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