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

..-----..-Px1,y1x,yz1x1y1,Matlab,,(x,y)(x1,y1)MATLABxy2003004..He-Ne=632.8nm12mm20mm50mm1r=12mm2r=20mm..3(r=50mm)15mm20mm30mm4564a=b=15mm..5(a=b=20mm)6(a=b=30mm)..MATLAB%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.^2+n.^2)/2/z);subplot(2,2,1);%22imagesc(I)%colormap([000;111])%axisimagetitle('')L=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);%H=fftshift(fft2(h));B=fftshift(fft2(I));%G=H.*B;%U=fftshift(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))r=12;a=1;b=1;I=zeros(N,N);[m,n]=meshgrid(linspace(-N/2,N/2-1,N));D=((m-a).^2+(n-b).^2).^(1/2);i=find(D=r);I(i)=1;%1