资源描述：
《光栅衍射的计算机仿真分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、1．光栅衍射的计算机仿真分析程序如下：clearL=0.000632;f=64;a=0.005;d=0.02;N=5;x=linspace(-3*pi,3*pi,1000);forj=1:1000u=(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);endNCLevels=255;Ir=NCLevels*I;figurecolormap(gray(NCLevels));subplot(2,1,2),image(x,I,Ir)subplo
2、t(2,1,1),plot(x,I(:)/max(I));N=5N=50N=11．光学干涉的计算机仿真实验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);axisoff;axissquare1．多光束干涉光的强度R=0.90;d=10;n=1.50;f=500;L=589.3e-6;
3、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),axissqure1.单缝衍射clc;clear;a=input('请输入单缝的宽度：a=');wavelength=input('请输入单色光的波
4、长：λ=');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.60.601])圆孔衍射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*
5、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);axissquare;axisoff双缝衍射f=64;b=0.005;d=0.02;L=632.8e-6;NCLevel=255;fori=1:4N=i^2+4;x=linspace(-3*pi,3*pi,1000);forj=1:1000u=(pi*b/L)*((x(j)/(sqrt(x(j)^2+f^2))));I(j)=(sin(u)/u
6、).^2*(sin(d*u*N/b)/sin(d*u/b)).^2;endIr=NCLevel*I;I1=Ir';colormap(gray(NCLevel));subplot(1,4,i),image(x,I,I1);axisoff;holdon;end三角形衍射图形A=imread('C:MATLAB6p5p1work5.bmp');B=fft2(A);C=abs(B).^2;D=fftshift(C);E=max(D(:));G=D/E;imagesc(sqrt(G))axisoff;axissquare;椭圆形衍射图形梯形衍射图形双圆孔衍射正六边形
7、的衍射方孔衍射矩形孔衍射双方孔衍射正方形分形垫片衍射图像正方形分形垫片衍射图像三角分形垫片衍射图像三角分形垫片衍射图像Koch雪花曲线Cantor集合垫片衍射图像布拉格晶格衍射图像
