光栅衍射的计算机仿真分析

1、1 光栅衍射的计算机仿真分析程序如下：clearL=0.;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+f2); 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)subplot(2,1,1),plot(x,I(:)/max(I); N=

2、5 N=50 N=12 光学干涉的计算机仿真实验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 square3 多光束干涉光的强度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

3、(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 squre1. 单缝衍射clc;clear;a=input( 请输入单缝的宽度：a= );wavelength=input( 请输入单色光的波长：= );t=- 0.6:0.001:0.6; %sina1=(1000*a/wavelength)*t;y10=sinc(a1)y1=y10.*

4、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

5、,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=i2+4; x=linspace(-3*pi,3*pi,1000); for j=1:1000 u=(pi*b/L)*(x(j)/(sqrt(x(j)2+f2); 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:MATLAB6p5p1work5.bmp);B=fft2(A);C=abs(B).2;D=fftshift(C);E=max(D(:);G=D/E;imagesc(sqrt(G)axis off;axis square;椭圆形衍射图形梯形衍射图形双圆孔衍射正六边形的衍