(完整word版)圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学

1、工程光学综合练习-圆孔、矩孔的菲涅尔衍射模拟圆孔和矩孔的菲涅尔衍射模拟由惠更斯菲涅尔原理可知接收屏上的P点的复振幅可以表示为£(P)=c“£(Q)二处K(e)do(1)其中£(Q)为衍射屏上的复振幅分布，K(e)为倾斜因子。根据基尔霍夫对此公式的完善,有1C=iA1+cosOK(0)=-«1设衍射屏上点的坐标为(x1,yI)，接收屏上点的坐标为(x,y),衍射屏与接收屏间距离为Z1,当满足菲涅尔近似条件时，即3('xl)2+(7yi)2max«万(2)此时可得到菲涅尔衍射的计算公式户G,y)=exp(ikzjiAzidxrdyr把上式指

2、数项中的二次项展开，并改写成傅里叶变换的形式，可以写成E(x,y)=eX1Zi)=xp/(/+y2)jexp+yf)(4)上式为菲涅尔衍射的傅里叶变换表达式，它表明除了积分号前面的一个与xl、yl无关的振幅和相位因子外，菲涅尔衍射的复振幅分布是孔径平面的复振幅分布和一个二次相位因子乘积的傅里叶变换。相对于夫琅和费衍射而言，菲涅尔衍射的观察屏距衍射屏不太远。在菲涅尔衍射中，输入变量和输出变量分别为衍射孔径平面的光场分布和观察平面的光场以及光强分布，考虑到这三个量都是二维分布，而且Matlab主要应用于矩阵数值运算,所以本程序选择用二维矩阵来存储衍射孔径平面和观察平面的场分布，并分别以矩阵的列数和

3、行数来对应平面的直角坐标值(x,y)以及(xi,yi)o孔菲涅尔衍射用MATLAB分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同，考虑到运算量的要求，采样点数不能过多，所以每个屏的X和y方向各取200到300点进行运算。根据式(4),选取合适的衍射屏和接收屏尺寸和相距的距离，模拟结果如下：取典型的He-Ne激光器波长A=632.8nm,固定衍射屏和接收屏尺寸和相距的距离，分别取不同的圆孔半径，得到以下三组衍射图样，其圆孔半径分别为12mm,20mm,50mm同孔形状3103回2CD253vn衍封国近1DO2C030口图1(r=12mm)今射回在图2(r=20mm)300250200

4、15010050国孔形状1002DD3CO衍射图标图3(r=50mm)三、矩孔的菲涅尔衍射步骤与上述相同，仅需改变与衍射屏形状对应的矩阵。这里选择矩孔的长宽相等,分别为15mm,20mm,30mm,其衍射图样及强度分布如图4、5、6桁射屏刑状30025020015010050100200300衍射后的图样衍射屏形状30D25020015010D50100200300衍射后的图样-200.200图5(a=b=20mm)衍射所形状100200300衍射后的图样-3J0-200图6(a=b=30mm)四、MATLAB程序%圆孔图像画在2行2列的第一个位置%画衍射屏的形状%颜色以黑白区分%所仃长度单位

5、为毫米lamda=632.8e-6;k=2*pi/lamda;z=1000000;%先确定衍射屏N=300;%圆屏采样点数a=lS;b=15;m,n=meshgrid(linspace(-N/2/N/2-l/N);l=rect(m/(2*a).*rect(n/(2*b);q=exp(j*k*(m.A2+n.A2)/2/z);subplot(2,2,l);imagesc(l)colormap(000;111)axisimagetitle('衍射屏形状')若为圆孔，方框内替换为以下程序r=12;a=l;b=l;I二zeos(N,N);mrn=meshgrid(linspace(-N

6、/2/N/2-l/N);D=(m-a).A2+(n-b).A2).A(l/2);i=find(D<=r);l(i)=l;%孔半径范围内透射系数为1L=300;M=300;%取相同点数用于矩阵运算x,y=meshgrid(linspace(-L/2/L/2/M);h=exp(j*k*z)*exp(j*k*(x.八2+y.八2)/(2*z)/(j*lamda*z);%接4攵屏%网孔频谱%公式中为卷枳，空间域中相卷相当于频域中相乘%求逆变换，得到更振幅分布矩阵%归一化H=fftshift(fft2(h);B=fftshift(fft2(l);G=H.*B;U=fftshift(ifft2(G);Br=(U/max(U);subplot(2,2,2);imshow(abs(U);%画三维图形axisimage;colormap(hot)%f