光栅衍射的计算仿真分析.docVIP

光栅衍射的计算机仿真分析 程序如下： clear L=0.000632; f=64; a=0.005; d=0.02; N=5; x=linspace(-3*pi,3*pi,1000); for j=1:1000 u=(pi*a/L)*((x(j)/(sqrt(x(j)^2+f^2)))); I(j)=((sin(u)/u).^2)*((sin(d*u*N/a)/sin(d*u/a)).^2); end NCLevels=255; Ir=NCLevels*I; figure colormap(gray(NCLevels)); subplot(2,1,2),image(x,I,Ir) subplot(2,1,1),plot(x,I(:)/max(I)); N=5 N=50 N=1 光学干涉的计算机仿真实验 3．牛顿环 L=589.3e-6;R=3000;NCLevels=256; x=-4:8/1000:4;y=x; [X,Y]=meshgrid(x,y); I=cos((pi*(X.^2+Y.^2)/R+L/2)/L).^2; colormap(gray(NCLevels)); Ir=256-I*NCLevels; image(x,y,Ir); axis off; axis square 多光束干涉光的强度 R=0.90;d=10;n=1.50;f=500; L=589.3e-6;NCLevels=255; x=linspace(-5,5,1000);y=x; [X,Y]=meshgrid(x,y); i=asin(n*sin(atan(sqrt(X.^2+Y.^2)/f))); Zeta=4*pi*d.*cos(i)/L; I=1./(1+4*R.*sin(Zeta/2).^2/(1-R)^2); colormap(gray(NCLevels)); Ir=I*NCLevels; image(x,y,Ir),axis squre 单缝衍射 clc;clear; a=input( 请输入单缝的宽度：a= ); wavelength=input( 请输入单色光的波长：λ= ); t=- 0.6:0.001:0.6; %sinθ a1=(1000*a/wavelength)*t; y10=sinc(a1) y1=y10.*y10; subplot(2,1,1) ,image(2,t,y1*255) colormap(gray(3)) subplot(2,1,2), plot(t,y1) axis([- 0.6 0.6 0 1]) 圆孔衍射 f=600;a=0.03;L=600e-6; x=linspace(-20,20,100);y=x; [X,Y]=meshgrid(x,y); Zeta=atan(sqrt(X.^2+Y.^2)/f); alpha=2*pi*a*sin(Zeta)/L; I=(2*bessel(1,alpha)./alpha).^2; NCLevel=255;Ir=NCLevel*I; colormap(jet); plot3(X,Y,I) image(x,y,Ir*5);axis square;axis off 双缝衍射 f=64;b=0.005;d=0.02;L=632.8e-6;NCLevel=255; for i=1:4 N=i^2+4; x=linspace(-3*pi,3*pi,1000); for j=1:1000 u=(pi*b/L)*((x(j)/(sqrt(x(j)^2+f^2)))); I(j)=(sin(u)/u).^2*(sin(d*u*N/b)/sin(d*u/b)).^2; end Ir=NCLevel*I;I1=Ir; colormap(gray(NCLevel)); subplot(1,4,i),image(x,I,I1); axis off; hold on; end 三角形衍射图形 A=imread(C:\MATLAB6p5p1\work\5.bmp); B=fft2(A);C=abs(B).^2; D=fftshift(C); E=max(D(:)); G=D/E; imagesc(sqrt(G)) axis off; axis square; 椭圆形衍射图形 梯形衍射图