| 查看: 3368 | 回复: 0 | ||
[求助]
(matlab)角谱法程序解释。。完全不懂。。求解释。。越清晰越好
|
|
function varargout = ASP( U1,z,dxy1,lambda,varargin ) %ASP Diffraction calculation by Angular Spectrum Propagation method % Syntax: % [U2,dxy2] = ASP(U1,z,dxy1,'PropertyName','PropertyValue',...) % [U2,...] = ASP(U1,z,dxy1,...) % ASP(U1,z,dxy1,...) % % U1 is the wavefront of the object plane % U2 is the wavefront of the diffraction plane % U1 and U2 are all two-dimensional array % size of U1 and U2 are even % z is the distance between object plane and diffraction plane % dxy1 is the sampling distance of the object, dxy1=[dx1,dy1] % lambda is the wavelength of the laser % % if there is no output, image on diffraction plane will be displayed % else, no image will be displayed % % the origin of coordinates is at M/2+1,N/2+1 % %-------------------------------------------------------------------------- % PropertyName and PropertyValue: % % piz {off} | t % pad the input array with zeros % off - do not pad % t - paste U1 to an all zeros array whose size is t*size(U1) % % iiz {off} | t % interpolate the input array with zeros % off - do not interpolate % t - interpolate t-1 points between every two points % % pfz {off} | t % pad the frequency array with zeros % off - do not pad % t - paste the frequency array(F) to an all zeros array whose % size is t*size(F) % % FTmethod {fourier} | b,FSN,f0 % which method to be use to realize the fourier transform % fourier - traditional fourier transform using single fft % b,FSN,f0 - FAFS using 3 ffts, b,FSN,f0 are the parameters needed % % IFTmethod {invfourier} | b,SSN,xy0 % which method to be use to realize the inverse fourier transform % invfourier - traditional inverse fourier transform using single fft % b,SSN,xy0 - IFASS using 3 ffts, b,SSN,xy0 are the parameters needed %-------------------------------------------------------------------------- error(nargchk(4,14,nargin)) if nargout>2 error('Too many output arguments') end for n=1:length(varargin) if ~ischar(varargin{n}) error('Property names and values must be characters') end end % ----------------------construct the property array---------------------- parray=[0,0,0,0,0]; % initialize the property array, the property array % is used to store the property specified by the user freepv={'1','1','1','',''}; % initialize the free property value array, free property % values do not fall into those values in the property structure property=struct('name',{'piz','iiz','pfz','FTmethod','IFTmethod'}); % construct the property structure, property(1).value={'off',''}; % which stores all possible properties property(2).value={'off',''}; property(3).value={'off',''}; property(4).value={'fourier',''}; property(5).value={'invfourier',''}; for l=1:2:length(varargin) namefound=0; valuefound=0; m=0; n=0; while ~namefound && m if strcmp(varargin{l},property(m).name) namefound=1; end end while namefound && ~valuefound && n if strcmp(varargin{l+1},property(m).value{n}) valuefound=1; parray(m)=n; end end %---------------------------------------------- if namefound==1 && valuefound==0 parray(m)=3; freepv{m}=varargin{l+1}; valuefound=1; end %---------------------------------------------- if namefound==0 error('wrong property name') elseif valuefound==0 error('wrong property value') end end %-------------------------------------------------------------------------- % k=2*pi/lambda; if parray(1)==3 t=str2double(freepv{1}); U1=paste(zeros(t*size(U1)),U1); end if parray(2)==3 t=str2double(freepv{2}); temp=U1; U1=zeros(t*size(U1)); U1(1:t:end,1:t:end)=temp; clear temp; end %---------------------- implement fourier transform ---------------------- if parray(4)~=3 [F,dfxy(1),dfxy(2)]=fourier(U1,dxy1(1),dxy1(2)); else eval(['[F,dfxy] = FAFS(U1,dxy1,',freepv{4},');']); end clear U1 if parray(3)==3 t=str2double(freepv{3}); F=paste(zeros(t*size(F)),F); end %---------- construct the transfer function in frequency domain ---------- [M,N]=size(F); H=zeros(M,N); for m=1:M for n=1:N fx=(n-N/2-1)*dfxy(1); fy=-(m-M/2-1)*dfxy(2); H(m,n)=exp(1i*2*pi*z*sqrt((1/lambda)^2-fx^2-fy^2)); end end %------------------------------------------------------------------------- F=F.*H; clear H; %------------------ implement inverse fourier transform ------------------ if parray(5)~=3 && nargout~=0 [U2,dxy2(1),dxy2(2)]=invfourier(F,dfxy(1),dfxy(2)); elseif parray(5)==3 && nargout~=0 eval(['[U2,dxy2] = IFASS(F,dfxy,',freepv{5},');']); elseif parray(5)~=3 && nargout==0 invfourier(F,dfxy(1),dfxy(2)); else eval(['IFASS(F,dfxy,',freepv{5},');']); end %-------------------------------------------------------------------------- switch nargout case 1 varargout{1}=U2; case 2 varargout{1}=U2; varargout{2}=dxy2; end |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ASP.m
2012-07-12 20:21:26, 5.26 K
» 猜你喜欢
听说今天filecode变了
已经有30人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
基金中了
已经有4人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
找到一些相关的精华帖子，希望有用哦~
求助 谁帮我解释一下XPS结合能和XAS吸收边的关系
已经有7人回复
- 光纤拉曼增益谱
已经有7人回复
- 色谱出现一个未知大峰，请问是怎么回事？如图所示
已经有15人回复
- PCR出现不能解释的条带
已经有17人回复
- 关于这幅图，请给一个合理的解释，谢谢各位
已经有11人回复
- 求助：1，胶束颜色解释？2，为什么胶束经过微孔滤膜过滤后不蓝了呢？
已经有5人回复
- the co-occurring secologanic acid怎么解释？
已经有6人回复
- 看不懂一句英文，求解释，气质联用文章里的方法文章
已经有22人回复
- 刚刚做完了热失重谱图，但是看不懂，求解释
已经有15人回复
- pcr过程中遇到的诡异现象，求高手解释一下
已经有8人回复
- 最近出现的永动机视频，看不懂，求解释
已经有12人回复
- 【求助】求解释：荧光光谱法研究变性剂对蛋白质构象影响，峰强变而峰位不发生变化的可
已经有7人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50