傅里叶变换相位解包裹程序

傅里叶变换相位解包裹程序，最原始的解包裹程序。

function varargout = fourier（ Uvarargin ）
%FOURIER Two-dimensional fourier transform
%  Syntax:
%  [Fdfxdfy] = fourier（Udxdy）
%  F = fourier（U）
%  fourier（U...）
%
%  U is the wavefront on the object plane
%  U must be a two-dimensional array
%  size of U must be even
%  F is the frequency distribution of U
%  dx is the sampling distance along the axis x in space domain
%  dy is the sampling distance along the axis y in space domain
%  if dx and dy are not inputted they are set to 1
%  dfx is the sampling distance along the axis fx in frequency domain
%  dfy is the sampling distance along the axis fy in frequency domain
%  
%  if there is no output frequency image will be displayed
%  else no image will be displayed
%
%  the origin of coordinates is at M/2+1N/2+1
%
%  the relation between sampling distance on the two planes is:
%  dfx=1/dx/N dfy=1/dy/M
%
error（nargchk（13nargin））
if nargout>3
    error（‘Too many output arguments‘）
end
[MN]=size（U）;
temp=zeros（MN）;
for m=1:M
    for n=1:N
        temp（mn）=（-1）^（m+n）;
    end
end
switch nargin
    case 1
        dx=1;
        dy=1;
    case 2
        dx=varargin{1};
        dy=dx;
    case 3
        dx=varargin{1};
        dy=varargin{2};
end
F=temp.*fft2（U.*temp）*（-1）^（（M+N）/2）*dx*dy;
switch nargout
    case 0
        f=abs（F）;
        OpticImage（f1/dx/N1/dy/M）;xlabel（‘fx‘）;ylabel（‘fy‘）;title（‘Amplitude of Frequency Spectrum‘）;
    case 1
        varargout{1}=F;
    case 2
        varargout{1}=F;
        varargout{2}=1/dx/N;
    case 3
        varargout{1}=F;
        varargout{2}=1/dx/N;
        varargout{3}=1/dy/M;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-07-04 11:05  傅里叶变换相位解包裹程序\
     文件        1611  2006-08-12 15:02  傅里叶变换相位解包裹程序\fourier.m
     文件        1655  2006-08-12 14:52  傅里叶变换相位解包裹程序\invfourier.m
     文件    26469438  2011-06-14 16:07  傅里叶变换相位解包裹程序\jiao.mat
     文件        2839  2006-11-17 15:01  傅里叶变换相位解包裹程序\OpticImage.m
     文件        1810  2006-08-12 14:59  傅里叶变换相位解包裹程序\paste.m
     文件         297  2011-05-24 14:57  傅里叶变换相位解包裹程序\Readme.txt
     文件        6513  2011-01-13 15:21  傅里叶变换相位解包裹程序\UBF.m
     文件        1465  2011-06-14 16:10  傅里叶变换相位解包裹程序\Unwrapping.asv
     文件        1466  2011-01-13 14:46  傅里叶变换相位解包裹程序\Unwrapping.m
     文件         615  2005-09-10 11:42  傅里叶变换相位解包裹程序\zero2one.m