% 下面的x1 x2 y1 y2模拟的
x1 = 0:1.2:10
y1 = sin(x1)
x2 = 0:0.1:10
y2= sin(x2)+0.1
y1_linear = interp1(x1,y1,x2)
% y1_linear 就是y1数据根据x2步长插值派渗首出来的对应数据
clear
Lambda=input('输入光的波长(单位为nm):取500')
Lambda=Lambda*1e-9
d=input('输入两个缝的间距(单位为nm):(取2)')
d=d*0.001
z=input('输入缝到屏的距尺汪离(单位宴困码为m):(取1)')
ymax=5*Lambda/dxs=ymax
Ny=101ys=linspace(-ymax,ymax,Ny)
for i=1:Ny
L1=sqrt((ys(i)-d/2).^2+z.^2)
L2=sqrt((ys(i)+d/2).^2+z.^2)
phi=2*pi*(L2-L1)/Lambda
B(i,:)=4*cos(phi/2).^2
N1=11,dL=linspace(-0.1,0.1,N1)
Lambda1=Lambda*(1+dL')
phi1=2*pi*(L2-L1)./Lambda1
B(i,:)=sum(4*cos(phi1/2).^2)/N1
end
NCLevels=255
Br=(B/4.0)*NCLevels
subplot(1,4,3),image(xs,ys,Br)
set(gca,'color','w')
colormap(gray(NCLevels))
subplot(1,4,4),plot(B(:),ys)
pause,subplot(1,1,1)
你看看这个效果怎么样,因为我不知道这个的物理背景,我只是晌哪做了一些语法的检查
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)