The computation of two-dimensional Fourier transform for point spread function of diffraction limited pupil can be transformed to the computation of one-dimensional Fourier and Hankel transform by use of the relationship between exponential function and Bessel function and property of an expansion of a optical symmetrical pupil function into azimuthal Fourier series. The detailed numerical procedure based on one-dimensional digital algorithm is completed on MATLAB. In comparison with the classic two-dimensional Fourier transform algorithm, the one-dimensional method of Fourier and Hankel transform can obtain intensity samples along arbitrary radial direction at the image plane. And compared with Fortran, Basic and C language etc., the MATLAB program has advantages of succinct language, great reductions of demand of mathematic basis and computer language, and it is a efficient tool for scientific and engineering research.