matlab仿真光束的传输特性

sinth=sqrt(r./(r+f^2));%角度 fai=2*pi*a*sinth./lam;úi=k*a*sinth hh=(2*BESSELJ(1,fai)).^2./fai.^2;%贝塞尔函数 b(:,i)=hh.*5000; end

subplot(1,2,1) image(xs,ys,b) colormap(gray(n)) subplot(1,2,2) b(:,m/2)

plot(ys,b(:,m/2))

（3）夫朗和费单缝衍射 clear all lam=500e-9; a=3;f=1; xm=3*lam*f/a; nx=50;

x=linspace(-xm,xm,nx); ny=50;

y=linspace(0,a,ny); for i=1:ny

sinphi=x/f;%角

af=(pi*a*sin(sinphi))/lam; I(i,:)=5*(sin(af)./af).^2; end

N=255;%确定灰度等级 Br=(I/max(I(1,:)))*N; subplot(1,2,1) image(x,y,Br);

colormap(gray(N));%颜色 subplot(1,2,2) plot(x,I(1,:));

（4）夫朗和费多缝衍射 clear all;

lamda=500e-9; %波长

N=2; %缝数，可以随意更改变换 a=2e-4;f=5;d=5*a;

ym=2*lamda*f/a;%选择坐标范围 xs=ym; n=1001;

ys=linspace(-ym,ym,n); for i=1:n sinphi=ys(i)/f;

alpha=pi*a*sinphi/lamda; fai=2*pi*d*sinphi/lamda;

I1=(sin(alpha)./alpha).^2;%单缝衍射因子

B(i,:)=I1*(sin(N*fai/2)./sin(fai/2)).^2;%多缝衍射光强的计算公式

B1=B/max(B);%归一化光强 end

NC=256; %确定灰度的等级 Br=(B/max(B))*NC; subplot(1,2,1) image(xs,ys,Br);

colormap(gray(NC)); %色调处理 subplot(1,2,2) plot(B1,ys,'k');

五、画出仿真图形 ?

?