MATLAB 实现物理光学模拟.pdfVIP

班级：光电0705班 姓名：罗钦 学号：U200713786 物理光学MATLAB作业 光谱扩展条纹分布特性 1.要求 杨氏干涉仪：a=20mm,D=1000mm,I =I ,d=1.5mm,λ=1µm 1 2 0 点光源位于光轴上，屏上任意点：当光谱范围扩展 ∆λ=0.4µm 时，观 察屏上的条纹图样。 2 技术 分析：对于谱线宽度为 ∆λ 的非单色光，干涉条纹消失的位置是 波长为 λ+∆λ/2 的成分的 k 级明纹与波长为 λ-∆λ/2 的成分的 k+1 级 明纹重合的位置。由于两成分在此位置上有同一光程差，根据光程差 与 明 纹 级 次 的 关 系 可 知 ， 条 纹 消 失 时 ， 应 满 足 ： （λ+∆λ/2 ） k =（ λ-∆λ/2 ）（ k+1 ）。 由 于 ∆λλ, 于 是 有 k λ/∆λ 。由该表达式可知，∆λ 越大，光的单色性越差，能观察到 的级次就越少。 3.操作 使用 MATLAB 模拟，对于单色光干涉时，其相邻的明纹与暗纹之 间的间距是Dλ/d,为了在屏上显示除了零级外共10 条明纹，程序中y 方向的计算范围为±5.5Dλ/d，用linspace 求出此范围内的空间点数 值。针对单色光和非单色光分别模拟仿真。由题意知光谱的宽度为中 心波长的 40%，近似的取100 根谱线，分别求每根谱线的总光强，然 后再对计算出的每根谱线产生的光强取平均值，这个平均值就是非单 色光的总光强，即可模拟出非单色光的干涉问题。 我并没有用彩色图标记其干涉图样，而是和教科书一样用灰度等 级表现出来，设灰度等级为 255 级，并将最大光强设为最大灰度 255， 即最大光强为白色。 我所做绘制的干涉图是二维的，并非老师给出的三维图样，但我 给出了干涉光强图。 4.结果截图 单色光干涉图及干涉光强 分析：这与实验现象完全相符，出现等间距的干涉条纹。 波长展宽为 ∆λ=0.4µm 分析：由于中间比较白（即灰度值大），故其对应的干涉光强就越大， 这在右面的图中可以看出来，这进一步验证了该程序的正确性。同时 从干涉图中看出，当 ∆λ较大时，光的单色性就较弱，所能见到的干 涉条纹就越少。 5．代码清单 clear all lamd=1e-6; %设定入射波长 d=1.5e-3; %缝间距 D=1; %屏缝间距 ym=5.5*lamd*D/d;%设定屏幕范围 xs=ym; n=201; ys=linspace(-ym,ym,n); for i=1:n r1=sqrt((ys(i)-d/2).^2+D^2); r2=sqrt((ys(i)+d/2).^2+D^2); phi=2*pi*(r2-r1)/lamd; I(i,:)=4*cos(phi/2).^2; %干涉光强 N1=100; %近似取100 根谱线 d1=linspace(-0.4,0.4,N1);%非单色光谱线宽度 lamd1=lamd*(1+d1); phi1=2*pi*(r2-r1)./lamd1; I1(i,:)=sum(4*cos(phi1/2).^2)/N1;%非单色光干涉强度 end N=255; %设定灰度等级 Ir=(I/4.0)*N; %最大光强为最大灰度（白色） I1r=(I1/4.0)*N; figure subplot(1,2,1) image(xs,ys,