求详细解释一下这个matlab源程序,并讲讲原理把单缝看作是np个分立的相干光源,屏幕上任意一点复振幅为np个光源照射结果的合成,对每个光源,光程差Δ=ypsinΦ,sinΦ=ys/D,光强I=I0(Σcosα)2+(Σsinα)2,
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/02 20:03:36
![求详细解释一下这个matlab源程序,并讲讲原理把单缝看作是np个分立的相干光源,屏幕上任意一点复振幅为np个光源照射结果的合成,对每个光源,光程差Δ=ypsinΦ,sinΦ=ys/D,光强I=I0(Σcosα)2+(Σsinα)2,](/uploads/image/z/6995984-32-4.jpg?t=%E6%B1%82%E8%AF%A6%E7%BB%86%E8%A7%A3%E9%87%8A%E4%B8%80%E4%B8%8B%E8%BF%99%E4%B8%AAmatlab%E6%BA%90%E7%A8%8B%E5%BA%8F%2C%E5%B9%B6%E8%AE%B2%E8%AE%B2%E5%8E%9F%E7%90%86%E6%8A%8A%E5%8D%95%E7%BC%9D%E7%9C%8B%E4%BD%9C%E6%98%AFnp%E4%B8%AA%E5%88%86%E7%AB%8B%E7%9A%84%E7%9B%B8%E5%B9%B2%E5%85%89%E6%BA%90%2C%E5%B1%8F%E5%B9%95%E4%B8%8A%E4%BB%BB%E6%84%8F%E4%B8%80%E7%82%B9%E5%A4%8D%E6%8C%AF%E5%B9%85%E4%B8%BAnp%E4%B8%AA%E5%85%89%E6%BA%90%E7%85%A7%E5%B0%84%E7%BB%93%E6%9E%9C%E7%9A%84%E5%90%88%E6%88%90%2C%E5%AF%B9%E6%AF%8F%E4%B8%AA%E5%85%89%E6%BA%90%2C%E5%85%89%E7%A8%8B%E5%B7%AE%CE%94%3Dypsin%CE%A6%2Csin%CE%A6%3Dys%2FD%2C%E5%85%89%E5%BC%BAI%3DI0%28%CE%A3cos%CE%B1%292%2B%28%CE%A3sin%CE%B1%292%2C)
求详细解释一下这个matlab源程序,并讲讲原理把单缝看作是np个分立的相干光源,屏幕上任意一点复振幅为np个光源照射结果的合成,对每个光源,光程差Δ=ypsinΦ,sinΦ=ys/D,光强I=I0(Σcosα)2+(Σsinα)2,
求详细解释一下这个matlab源程序,并讲讲原理
把单缝看作是np个分立的相干光源,屏幕上任意一点复振幅为np个光源照射结果的合成,对每个光源,光程差Δ=ypsinΦ,sinΦ=ys/D,光强I=I0(Σcosα)2+(Σsinα)2,其中α=2Δ/λ=πypys/λD
编写程序如下,得到图1
clear
lam=500e-9;
a=1e-3;D=1;
ym=3*lam*D/a;
ny=51;
ys=linspace(-ym,ym,ny);
np=51;
yp=linspace(0,a,np);
for i=1:ny
sinphi=ys(i)/D;
alpha=2*pi*yp*sinphi/lam;
sumcos=sum(cos(alpha));
sumsin=sum(sin(alpha));
B(i,:)=(sumcos^2+sumsin^2)/np^2;
end
N=255;
Br=(B/max(B))*N;
subplot(1,2,1)
image(ym,ys,Br);
colormap(gray(N));
subplot(1,2,2)
plot(B,ys);
求详细解释一下这个matlab源程序,并讲讲原理把单缝看作是np个分立的相干光源,屏幕上任意一点复振幅为np个光源照射结果的合成,对每个光源,光程差Δ=ypsinΦ,sinΦ=ys/D,光强I=I0(Σcosα)2+(Σsinα)2,
对于光学的东西不是太懂,能做到的只是帮助你更容易地理解程序每一句话,
最好将下面的程序复制粘贴到matlab的编辑窗口中再阅读
clear % 清除matlab工作内存
lam=500e-9; %定义一个常数λ:500乘以10的-9次方
a=1e-3;D=1; %定义常数a:10的-3次方;常数D=1
ym=3*lam*D/a; %定义ym
ny=51; %定义取点个数 ny=51
ys=linspace(-ym,ym,ny); %定义一个数组ys,从-ym到ym之间取ny个点
np=51; %定义取点个数 np=51
yp=linspace(0,a,np); %定义一个数组yp,从0到a之间取np个点
for i=1:ny %循环语句,一共循环51次
sinphi=ys(i)/D; %计算得到 sinΦ,
alpha=2*pi*yp*sinphi/lam; % 计算得到 α=2Δ/λ=πypys/λD(由于yp是一个向量,这里得到的α也是一个向量,包含51个数)
sumcos=sum(cos(alpha)); % 计算得到 Σcosα
sumsin=sum(sin(alpha)); % 计算得到 Σsinα
B(i,:)=(sumcos^2+sumsin^2)/np^2; % 计算得到光强 I=I0(Σcosα)2+(Σsinα)2,并存储为矩阵B的第i行
end
N=255;
Br=(B/max(B))*N; %计算得到光强分布矩阵
subplot(1,2,1) %将绘图区分成左右两部分,首先在左侧绘图区绘图
image(ym,ys,Br); %绘制光强分布图
colormap(gray(N)); %更改图像模式为灰度图像
subplot(1,2,2) %然后选择在右侧绘图区绘图
plot(B,ys); %绘制光强分布曲线