资源描述：
《圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、工程光学综合练习圆孔、矩孔的菲涅尔衍射模拟圆孔和矩孔的菲涅尔衍射模拟一、原理由惠更斯-菲涅尔原理可知接收屏上的P点的复振幅可以表示为exp(??????软?？=?????Q)K(9)d(T(1)????其中？?Q)为衍射屏上的复振幅分布，K(8为倾斜因子。根据基尔霍夫对此公式的完善,有1C二 ????1+????????K(0)=〜1设衍射屏上点的坐标为(X1,y1)，接收屏上点的坐标为(x,y),衍射屏与接收屏间距离为当满足菲涅尔近似条件时，即??C8??[(??-?如+(??-??)2]?????????(2)此时可得到菲涅尔衍射的计算公式-exp(??????8-????.°??x
2、，y)=???????，??)exp{2??[(?2??)2+(????)2]}???????(3)I八1・・-OO2.?软x,y)exp(??????FV?exp[把上式指数项中的二次项展开，并改写成傅里叶变换的形式，可以写成缪???+??)]f储??,？?)exp[考冷+??)]}2??2??上式为菲涅尔衍射的傅里叶变换表达式，它表明除了积分号前面的一个与x1、y1无关的振幅和相位因子外，菲涅尔衍射的复振幅分布是孔径平面的复振幅分布和一个二次相位因子乘积的傅里叶变换。相对于夫琅和费衍射而言，菲涅尔衍射的观察屏距衍射屏不太远。在菲涅尔衍射中，输入变量和输出变量分别为衍射孔径平面的光场分
3、布和观察平面的光场以及光强分布，考虑到这三个量都是二维分布，而且Matlab主要应用于矩阵数值运算，所以本程序选择用二维矩阵来存储衍射孔径平面和观察平面的场分布，并分别以矩阵的列数和行数来对应平面的直角坐标值(x,y)以及(xi,y1)。二、圆孔菲涅尔衍射用MATLAB分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同，考虑到运算量的要求，采样点数不能过多，所以每个屏的x和y方向各取200到300点进行运算。根据式(4),选取合适的衍射屏和接收屏尺寸和相距的距离，模拟结果如下：取典型的He-Ne激光器波长入=632.8nm,固定衍射屏和接收屏尺寸和相距的距离，分别取不同的圆孔半径
4、，得到以下三组衍射图样，其圆孔半径分别为12mm,20mm,50mm区九档状宿射国祥图2(r=20mm)图2(r=20mm)图1(r=12mm)国孔刑状衍射固岸1G0200300图2(r=20mm)图3(r=50mm)三、矩孔的菲涅尔衍射步骤与上述相同，仅需改变与衍射屏形状对应的矩阵。这里选择矩孔的长宽相等,分别为15mm,20mm,30mm,其衍射图样及强度分布如图4、5、6衍射屏解状衍射后的图杠图4(a=b=15mm)衍射屏彤状箱射后的国以图5(a=b=20mm)和附居系状JUO25TJ2in15010050100200300-20C-200200图6(a=b=30mm)四、MATL
5、AB®序%所有长度单位为毫米lamda=632.8e-6;k=2*pi/lamda;z=1000000;%先确定衍射屏N=300;%圆屏采样点数a=15;b=15;[m,n]=meshgrid(linspace(-N/2,N/2-1,N));I=rect(m/(2*a)).*rect(n/(2*b));q=exp(j*k*(m.A2+n.A2)/2/z);subplot(2,2,1);imagesc(I)colormap([000;111])axisimagetitle('衍射屏形状')L=300;M=300;[x,y]=meshgrid(linspace(-L/2,L/2,M));若为
6、圆孔，方框内替换为以下程序r=12;a=1;b=1;I=zeros(N,N);[m,n]=meshgrid(linspace(-N/2,N/2-1,N));D=((m-a).A2+(n-b).A2).A(1/2);i=find(D<=r);I(i)=1;%孔半径范围内透射系数为1%圆孔图彳t画在2行2列的第一个位置%画衍射屏的形状%颜色以黑白区分%取相同点数用于矩阵运算h=exp(j*k*z)*exp((j*k*(x.A2+y.A2))/(2*z))/(j*lamda*z);%接收屏H=fftshift(fft2(h));B=fftshift(fft2(I));G=H.*B;U=ffts
7、hift(ifft2(G));Br=(U/max(U));subplot(2,2,2);imshow(abs(U));axisimage;colormap(hot)%figure,imshow(C);title('衍射后白内图样');subplot(2,2,3);mesh(x,y,abs(U));subplot(2,2,4);plot(abs(Br))%圆孔频谱%公式中为卷积，空间域中相卷相当于频域中相乘%求逆变换，得到复振幅分布矩
