第8讲衍 射.pptVIP

衍 射 信息光电子科技学院 李红燕 衍 射 1. 菲涅耳衍射 2. 夫琅和费衍射 1. 菲涅耳衍射近似和菲涅耳衍射计算公式 1. 菲涅耳衍射近似和菲涅耳衍射计算公式 2. 夫琅和费衍射近似和夫琅和费衍射计算公式 波长为632.8um的平行光照明边长为1mm的方形孔，记录 面的尺寸为5.12mm*5.12mm 令菲涅耳衍射近似 夫琅和费衍射近似 模拟矩孔的菲涅耳衍射 第一步：生成矩孔 clear clc t0=clock; %生成矩孔 r=512;c=512;L0=0.00512;%给出输入图像的像素数及图象尺寸，单位：米 a=zeros(r,c); a((r/2-5):(r/2+4),(c/2-5):(c/2+4))=1; figure;imshow(a,[ ]);title(矩孔图形);axis on; 第二步：计算物光传递到屏的衍射过程 lamda=6328*10^(-10); k=2*pi/lamda; x1=linspace(-L0/2,L0/2,r);y1=linspace(-L0/2,L0/2,c); [x1,y1]=meshgrid(x1,y1);%给出输入图像的坐标 z1=2; %图像到记录面的距离，单位：米 x=linspace(-L0/2,L0/2,r);y=linspace(-L0/2,L0/2,c); %记录面的尺寸 [x,y]=meshgrid(x,y);%给出记录面的坐标 F00=exp(j*k*z1)/(j*lamda*z1)*exp(j*k/2/z1*(x.^2+y.^2)); F=exp(j*k/2/z1*(x1.^2+y1.^2)); Fft=fftshift(fft2(a.*F)); %做FFT变换 Exy=F00.*Fft; %记录面上的复振幅分布 I=Exy.*conj(Exy); %记录面上的光强分布 I=I/max(max(I)); %强度归一化 figure;imshow(I,[0 0.2]);colormap(pink);title(矩孔的菲涅耳衍射);axis on; xlabel(像素数);ylabel(像素数); figure;plot(I(256,:));title(强度); t1=etime(clock,t0) 模拟多缝的菲涅耳衍射 第一步：生成多缝 clear clc t0=clock; %生成多缝 r=512;c=512;L0=0.00512;%给出输入图像的像素数及图象尺寸，单位：米 a=zeros(r,c); for k=1:512 if mod(k,6)==0 a(:,k)=1; else a(:,k)=0; end end figure;imshow(a,[ ]);title(多缝图形);axis on; 第二步：与矩孔的菲涅耳衍射相同 练 习 1. 编写一个单缝衍射的模拟程序 2. 编写一个圆孔衍射的模拟程序 * * 近似后的基尔霍夫衍射公式 当 大到使得上式第三项以后各项对相位kr的作用远小于 时，即： 第三项以后各项便可忽略，只取前两项来表示r , 即： (1) (2) (3) ▲ 所以菲涅耳衍射公式为 也可以写为 (4) (5) 当 很大，使得（3）式中第四项对相位的贡献远小于 ，即： 这时，第四项可略去 夫琅和费衍射计算公式为 (6) (7) (8) ▲ 计算出菲涅耳衍射距离为： 夫琅和费衍射距离为： *