Matlab数字衍射光学实验(二)分析.doc

clear; close all; lamda=0.6328;%入射光波长,单位 um z=200.0e3;%传播距离,单位 um h=2e3;w=2e3;%模拟光波的大小,单位 um dx=10.0;dy=10.0;%抽样间隔,单位 um M=h/dx;N=w/dy;%抽样点数 a0=1;%初始振幅大小为1 x0=0;y0=0;%模拟光源位于坐标原点 k=2*pi/lamda;%波数 %---------球面波表示-------% [x,y]=meshgrid(-M/2*dx:(M-1)/2*dx,-N/2*dy:(N-1)/2*dy);%M个点,两点之间的距离为dx SphFunct=a0/z.*exp(i*k*z).*exp(i*pi/(lamda*z).*((x-x0).^2+(y-y0).^2));%球面波 I=abs(SphFunct).^2;%光强度 I1=I./max(max(I));%归一化强度,二维所以求两次max() Ph=angle(SphFunct);%相位,函数angle()表示取相位 figure(1); imagesc(I1);%做出球面波的强度分布图，没有强度变化，均匀光斑,这里做了彩图 figure(2); P=(Ph+pi)/(2*pi);%归一化相位 imshow(P);%画出球面波的相位分布图 figure(3); imagesc(P);