[求助]Matlab循环语句模拟菲涅尔衍射算法求助 [复制链接]

看了不少帖子发现大多数都是用傅里叶变换法来模拟菲涅耳衍射的算法，我想问问有没有大神指导下如何通过循环语句来计算菲涅尔衍射的结果呢？ 10H)^p%3+  
=N 5z@;!  
思路是这样的，首先生成一个矩阵来模拟通过孔径光阑的光，以方孔为例，比如100*100的矩阵中心只有20*20的元素为1，其余为0。在光屏部分对应的分别计算模拟接收屏100*100矩阵每个点元素的衍射结果来描述衍射的结果。 p<:!)kt  
P (S>=,Y&  
从网上找到一个类似的例子是计算夫琅禾费衍射的，代码如下： qoAj] ")  
ES40?o*]x  
R=0.1; rb{P :MX  
lambda=1.064e-3; [|l?2j\  
k=2*pi/lambda; 2<}NB?f`N  
z=1.0e3; *YlV-C<}W"  
r=linspace(0,2*1.22*lambda/2/R*z,201); a@R]X5[O  
eta=linspace(0,2*pi,201); yeHDa+}  
[rho,theta]=meshgrid(r,eta); @}kv-*  
[x,y]=pol2cart(theta,rho); V <bd;m  
r0=linspace(0,R,201); Qqlup  
eta0=linspace(0,2*pi,201); RVeEkv[qp  
[rho0,theta0]=meshgrid(r0,eta0); `Qqk<o  
[x0,y0]=pol2cart(theta0 ,rho0); $HJwb-I  
deta=R/200*2*pi/200; : "1XPr  
E2=zeros(201,1); IUwm}9Q!  
for gk=1:201 8n>9;D5n  
   for m=1:200 `:axzCrCfR  
     for n=1:201 Zv#Ll@v  
E2(gk)=E2(gk)-j/lambda/z*exp(((x(1,gk)^2+y(1,gk)^2)/z/2+z)*j*k)*exp(j*k*(x(1,gk)*x0(m,n)+y(1,gk)*y0(m,n))/z)*deta*rho0(m,n); ~9j%Hm0ht  
     end GQ |Mr{.;  
   end K#"O a h  
end 5<w g8y  
Ie=conj(E2).*E2; plot(r,Ie,'k',-r,Ie,'k'); k?'B*L_Mzv  
:^(>YAyHj^  
不明白的地方是，计算一维分布如何带入严格的菲涅耳-基尔霍夫标量衍射公式呢？不用傅里叶变换法的话很难计算积分啊。。因为不是很懂坛子里大神发的fft法模拟衍射，所以想用简单的思路来解决问题。。还有怎么把一维的计算结果可视化的表现出来呢？刚学matlab好多都不太明白啊

速度怎么样啊

看不懂