Zemax光学设计从入门到精通	Zemax光学设计从基础到实践	Zemax光学设计实践基础(陈家璧)	ZEMAX光学系统设计实战
讯技光电:VirtualLab Fusion独家供应商	光学设计与光学元件	CODE V,LightTools,RSoft,LucidShape	计算光学带来的成像革命

发帖回复

2864阅读
2回复

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

上一主题下一主题

离线tzq7961

实习生

发帖: 3

光币: 5

光券: 0

只看楼主倒序阅读楼主发表于: 2017-03-09

看了不少帖子发现大多数都是用傅里叶变换法来模拟菲涅耳衍射的算法，我想问问有没有大神指导下如何通过循环语句来计算菲涅尔衍射的结果呢？

思路是这样的，首先生成一个矩阵来模拟通过孔径光阑的光，以方孔为例，比如100*100的矩阵中心只有20*20的元素为1，其余为0。在光屏部分对应的分别计算模拟接收屏100*100矩阵每个点元素的衍射结果来描述衍射的结果。

从网上找到一个类似的例子是计算夫琅禾费衍射的，代码如下：

R=0.1;
lambda=1.064e-3;
k=2*pi/lambda;
z=1.0e3;
r=linspace(0,2*1.22*lambda/2/R*z,201);
eta=linspace(0,2*pi,201);
[rho,theta]=meshgrid(r,eta);
[x,y]=pol2cart(theta,rho);
r0=linspace(0,R,201);
eta0=linspace(0,2*pi,201);
[rho0,theta0]=meshgrid(r0,eta0);
[x0,y0]=pol2cart(theta0 ,rho0);
deta=R/200*2*pi/200;
E2=zeros(201,1);
for gk=1:201
   for m=1:200
     for n=1:201
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);
     end
   end
end
Ie=conj(E2).*E2; plot(r,Ie,'k',-r,Ie,'k');

不明白的地方是，计算一维分布如何带入严格的菲涅耳-基尔霍夫标量衍射公式呢？不用傅里叶变换法的话很难计算积分啊。。因为不是很懂坛子里大神发的fft法模拟衍射，所以想用简单的思路来解决问题。。还有怎么把一维的计算结果可视化的表现出来呢？刚学matlab好多都不太明白啊