基于MATLAB编程模拟光学实验

1、2022-2-5基于MATLAB编程模拟光学实验钱朝阳测试计量技术及仪器在物理学波动光学的部分，往往由于这部分知识较为抽象且要用到复杂的数学运算，有时很难掌握好这部分的知识。为了较好地解决这一问题，本文根据惠更斯-菲涅耳原理计算推导了夫琅禾费衍射场下单缝衍射的光强分布公式，建立了单缝衍射的Matlab实现模型，利用Matlab强大的计算功能，得出光强分布矩阵，然后用Matlab绘图函数得到衍射图样，同时绘制了相应的衍射光强分布图。并结合理论公式对衍射图样进行了分析。单缝衍射单缝衍射单缝衍射是光在传播过程中遇到障碍物，光波会绕过障碍物继续传播的一种现象。如果波长与缝、孔或障碍物尺寸相当或者更大时

2、，衍射现象最明显。把单色点光源放在透镜的焦点上，经过透镜后的单色平行光垂直照射衍射屏时，在屏后面不同距离上会观察到一些衍射现象，其中当屏远离到足够大的距离后，光斑中心出现一个较大的亮斑，外围是一些较弱的明暗相间的同心圆环，此后再往外移动，衍射花样出现稳定分布，中心处总是亮的，只是半径不断扩大而已，这种衍射称为夫琅禾费衍射。=2波程差 相位差光的衍射现象是光的波动性的一种表现，通过对光的各种衍射现象的研究可以深入具体地了解光的波动性。惠更斯原理认为：传播中的波面上任何一点都可以认为是一个新的次波源，由这些次波源发出的次波是球面波，这些次波的公共包络面就是下一个时刻的波面。菲涅尔将其进一步具体化，

3、认为可以将波面S分成许多面积元ds，每个ds在P点的振动可写为： 这便是计算衍射场的常用公式，即惠更斯-菲涅尔原理的表达式。 2(t)( )cosrdEkdsr2(t)()( )cosSrE Pkdsr单色平行光通过单缝,可看作N条光束的叠加来计算,单缝衍射的几何关系如图所示.把单缝看作一排等间隔的光源,N个光源分布,P处的光强为这N个光源照射结果的合成。22()rxyL设计成设计成MATLABMATLAB程序如下程序如下: :clearlamba=500e-9; %波长a=1e-3;L=1; %a是狭缝宽度，L屏与障碍物的距离xm=3*lamba*L/a; %屏幕上y方向的范围最大值n=51

4、; %屏幕上的点数xs=linspace(-xm,xm,n);n=51; %缝上的点数yp=linspace(-a/2,a/2,n);for i=1:n %对屏幕和光缝上各点循环，计算缝点到屏幕位置的距离r=sqrt(xs(i)-yp).2+L2);phi=2*pi.*(r-L)./lamba;%相位差sumcos=sum(cos(phi); %数组求和sumsin=sum(sin(phi);B(i,:)=(sumcos2+sumsin2)/n2; %屏幕上的归一化光强endN=256; %确定灰度的等级Br=(B/max(B)*N; %定标subplot(1,2,1)image(xm,xs,

5、Br);colormap(gray(N); %色调处理subplot(1,2,2)plot(B,xs,k);杨氏双缝干涉杨氏双缝干涉clearclearlam=500e-9lam=500e-9a=2e-3;D=1;a=2e-3;D=1;ym=5ym=5* *lamlam* *D/a;xs=ym;D/a;xs=ym;n=101;ys=linspace(-ym,ym,n);n=101;ys=linspace(-ym,ym,n);for i=1:nfor i=1:nr1=sqrt(ys(i)-a/2).2+D2);r1=sqrt(ys(i)-a/2).2+D2);r2=sqrt(ys(i)+a/2)

6、.2+D2); r2=sqrt(ys(i)+a/2).2+D2); phi=2phi=2* *pipi* *(r2-r1)./lam;(r2-r1)./lam;B(i,:)=sum(4B(i,:)=sum(4* *cos(phi/2).2);cos(phi/2).2);endendN=255;N=255;Br=(B/4.0)Br=(B/4.0)* *N Nsubplot(1,2,1)subplot(1,2,1)image(xs,ys,Br);image(xs,ys,Br);colormap(gray(N);colormap(gray(N);subplot(1,2,2)subplot(1,2,2

7、)plot(B,ys)plot(B,ys)1 光栅衍射光栅衍射clear clear lam=500e-9;N=2;lam=500e-9;N=2;a=2e-4;D=5;d=5a=2e-4;D=5;d=5* *a;a;ym=2ym=2* *lamlam* *D/a;xs=ym;D/a;xs=ym;n=1001;n=1001;ys=linspace(-ym,ym,n);ys=linspace(-ym,ym,n);for i=1:nfor i=1:nsinphi=ys(i)/D;sinphi=ys(i)/D;alpha=pialpha=pi* *a a* *sinphi/lam;sinphi/lam

8、;beta=pibeta=pi* *d d* *sinphi/lam;sinphi/lam;B(i,:)=(sin(alpha)./alpha).2.B(i,:)=(sin(alpha)./alpha).2.* *(sin(N(sin(N* *beta)./sin(beta).2;beta)./sin(beta).2;B1=B/max(B);B1=B/max(B);endendNC=255;NC=255;Br=(B/max(B)Br=(B/max(B)* *NC;NC;subplot(1,2,1)subplot(1,2,1)image(xs,ys,Br);image(xs,ys,Br);colormap(gray(NC);colormap(gray(NC);subplot(1,2,2)subplot(1,2,2)plot(B1,ys);plot(B1,ys);参考文献：参考文献：11：胡守信、