光的干涉衍射计算机模拟仿真技术,实验报告之仿真(光的干涉与衍射)

实验报告之仿真(光的干涉与衍射)

(10页)

本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！

19.90 积分

大学物理创新性试验大学物理创新性试验实验项目：单缝实验项目：单缝﹑﹑双缝双缝﹑﹑多缝衍射现多缝衍射现 象仿真实验象仿真实验专业班级：材料成型及控制工程 0903 班 姓 名：曹惠敏 学 号：090201097目录目录 1 1 光的衍射光的衍射 2 2 衍射分衍射分类类 3 3 实验现实验现象象 4 4 仿仿真真模模拟拟 5 实验总结光的衍射光的衍射光在传播路径中，遇到不透明或透明的障碍物，绕过障碍物，产生偏离直线传播的现象称为光的衍射。光的衍射现象是光的波动性的重要表现之一.波动在传播过程中,只要其波面受到某种限制,如振幅或相位的突变等,就必然伴随着衍射的发生. 然而,只有当这种限制的空间几何线度与波长大小可以比拟时,其衍射现象才能显著地表现出来.所有光学系统,特别是成像光学系统,一般都将光波限制在一个特定的空间域内,这使得光波的传播过程实际上就是一种衍射过程.因此,研究各种形状的衍射屏在不同实验条件下的衍射特性,对于深刻理解衍射的实质,研究光波在不同光学系统中的传播规律分析复杂图像的空间频谱分布以及改进光学滤波器设计等具有非常重要的意义.随着计算机技术的飞速发展, 计算机仿真已深入各种领域。光的干涉与衍射既是光学的主要内容 , 也是人们研究与仿真的热点。由于光波波长较短,与此相应的复杂形状衍射屏的制作较困难,并且实验过程中对光学系统及环境条件的要求较高.因而在实际的实验操作和观察上存在诸多不便. 计算机仿真以其良好的可控性、无破坏、易观察及低成本等优点,为数字化模拟现代光学实验提供了一种极好的手段. 本次实验利用 MATLAB 软件实现对任意形状衍射屏的夫琅禾费衍射实验的计算机仿真。衍射分类衍射分类⒈⒈菲涅尔衍射菲涅尔衍射 菲涅尔衍射菲涅尔衍射：入射光与衍射光不都是平行光的衍射。惠更斯提出，媒质上波阵面上的各点，都可以看成是发射子波的波源，其后任意时刻这些子波的包迹，就是该时刻新的波阵面。菲涅尔充实了惠更斯原理，他提出波前上每个面元都可视为子波的波源，在空间某点 P 的振动是所有这些子波在该点产生的相干振动的叠加，称为惠更斯－菲涅尔原理。惠更斯－菲涅尔原理能定性地描述衍射现象中光的传播问题，成为我们解释光的各类衍射现象的理论依据。 2.夫琅禾费衍射夫琅禾费衍射 夫琅禾费衍射：夫琅禾费衍射：入射光与衍射光都是平行光的衍射。实验装置：实验装置：由基础光学可知,任意衍射屏的夫琅禾费衍射可借助两个透镜来实现.,如图 1 所示，位于透镜 L1 物方焦平面上的点源 S 所发出的单色球面光波经 L1 变换为一束平面光波,照射在衍射屏 AB 上. 按照平面波理论,衍射屏开口处的波前向各个方向发出次波方向彼此相同的衍射次波经透镜 L2 会聚到其像方焦平面的同一点 Pθ 上. 满足相长干涉条件时,该点为亮点;满足相消干涉时, 该点为暗点. 所有亮点和暗点的集合构成了该衍射屏的夫琅禾费衍射图样. 其次,从傅里叶光学角度, 任意衍射屏在单位振幅的单色平面波垂直照射下, 其夫琅禾费衍射光场复振幅即衍射屏透射系数的傅里叶变换, 而衍射图样实际上就是衍射屏的空间频谱强度分布。实验现象实验现象通过实验我们得到的衍射图像如下列图像所示： 图(1)——圆孔衍射 图(2)——单缝衍射 图(3)——双缝衍射 图(4)——光栅衍射图(1) 图(2)图(3) 图(4)仿真模拟仿真模拟1.Flash动画演示动画演示 2.计算机仿真计算机仿真 ⑴⑴单色光模拟单色光模拟1) 单缝. 设狭缝宽度为 a ,观察屏上点 Pθ 与透镜 L2 光心连线的方位角为 θ, 由几何成像理论, 此角正好也是相应平面波分量的方位角. 若取入射光波长为 λ, 透镜 L2 的焦距为f , 根据惠更斯- 菲涅耳原理[2 ] ,可得平面波垂直照射下的单缝夫琅禾费衍射图样的归一化强度为式中 α=πa sinθ/ λ, I ( P0) 为衍射图样中心点 P0 的强度. 若取Pθ 点到中心点 P0 的距离为 x′,则 sin θ= x′/ ( x′2 + f 2) 1/ 2 ,由此可以得到接收屏上任意位置 Pθ 到中心点 P0 的距离 x′与该的相对光强度之间的关系. 单缝单缝衍射衍射 MATLAB 程序程序设计设计： Clear lam=500e-9;a=1e-3;f=1;xm=3*lam*f/a; nx=51;xs=linspace(-xm,xm,nx);np=51;xp=linspace(0,a,np);for i=1:nxsinphi=xs(i)/falpha=2*pi*xp*sinphi/lam;sumcos=sum(cos(alpha));sumsin=sum(sin(alpha));B(i)=(sumcos^2+sumsin^2)/np^2;endN=255;Br=(B/max(B))*N;subplot(1,2,1)image(xm,xs,Br);colormap(bone(N));subplot(1,2,2) plot(B,xs);单缝衍射模拟图像 圆圆孔衍射：孔衍射：clear;f=700;a=0.04;wlr=700e-6;wlg=546.1e-6;wlb=435.8e-6;x=linspace(-30,30,800);[X,Y]=meshgrid(x);seta=atan(sqrt(X.^2+Y.^2)/f);aphr=2*pi*a*sin(seta)/wlr;aphg=2*pi*a*sin(seta)/wlg;aphb=2*pi*a*sin(seta)/wlb;Ir=(2*besselj(1,aphr)./aphr).^2;Ib=(2*besselj(1,aphb)./aphb).^2;Ig=(2*besselj(1,aphg)./aphg).^2;Iw=zeros(800,800,3);Iw=Ir;Iw=Ig;Iw=Ib;IO=255;Iw=IO*Iw;imshow(Iw);圆孔衍射模拟图像2) 多缝(朗琴光栅) . 设每个狭缝的宽度均为 a ,相邻两狭缝间不透明部分的宽度为 b , 则缝间距(光栅常量) 为 d = a + b. 同样取 θ 为衍射角, f 为透镜 L2 的焦距,衍射屏上透光的总狭缝数为 N ,则当平面光波垂直照射衍射屏时,沿 θ 方向的衍射光波在 Pθ 处的合振动的相对强度为式中 u =πa sin θ/λ, v =πd sin θ/λ= du / a , sin2 u/u2 称为单缝衍射因子, sin2 ( N v ) / sin2 v 称为缝间的干涉因子. 同样,若 Pθ 到 P0 的距离为 x′(θ) , 则由此可以得到相对衍射光强度与 x′(θ) 之间的关系.从以上两式可以看出,确定了仿真程序中的可调参数为入射光波长 λ、衍射屏的缝宽 a 、b ,透镜的焦距 f 后,只要求出接收屏上每一点的相对衍射光强度值,就可以绘出衍射图样. 杨杨氏双氏双缝缝干涉程序干涉程序设计设计：clear lam=600e-9 a= 2e-3; D= 1; ym= 5*lam* D /a; xs= ym; n= 101; ys= linspace( -ym, ym, n); for i= 1: n r1= sqrt( ( ys( i) - a /2).^2+ D^2); r2= sqrt( ( ys( i) + a /2).^2+ D^2); phi=2* pi* ( r2- r1)./lam;B(i) = sum( 4* cos( phi/2).^2); end N= 255; Br= ( B/4.0)* N subplot( 1, 2, 1) image( xs, ys, Br); colormap(bone(N) ); subplot( 1, 2, 2) plot(B, ys)杨氏双缝干涉模拟图像光光栅栅衍射：衍射：clear wavelength=0.0006328;f=100;b=0.005;d=0.02;n=5;angles=0; x=linspace(-4*pi,4*pi,1000); for j=1:1000u=(pi*b/wavelength)*(((x(j)/sqrt(x(j)^2+f^2)))+sin(angles));I(j)=((sin(u)/u).^2)*((sin(d*u*n/b)/sin(d*u/b)).^2); end hold on NCLevels=255; Ir=NCLevels*I; colormap(gray(NCLevels)); subplot(2,1,2),image(x,I,Ir) subplot(2,1,1),plot(x,I(:)/max(I));光栅衍射模拟图像⑵⑵白光模拟白光模拟实验以日光为仿真的白光光源。由于日光是一种连续谱的自然光, 因此, 要实现白光夫琅和费衍射实验的计算机仿真,须将日光的可见光谱均匀分割为若干份色光, 并计算出与各色光对应的色彩代码。本实验把日光的可见光谱均匀分割为 2 000 份色光。为便于利用 Matlab 软件进行仿真, 色光的色彩代码使用 RGB 模型表示。RGB 颜色模型立方体以红(R )、绿(G)、蓝(B)为坐标轴, 坐标值从 0~ 1 变化, 立方体范围内的每个彩色点都可以用 R、G、B 3 个参数表示, 分别表示红、绿、蓝基色的相对亮度。 1)白光夫琅和费矩孔衍射实验的计算机仿真夫琅和费矩孔衍射的衍射屏上任一点 P ( x, y )的光强为(1)式中, AU (axP)/(Kf), BU (bxP) /( Kf ), a、b 分别为衍射矩孔的长和宽(即矩孔沿 x、y 轴方向的长度)。依据式(1), 即可便捷地编写出白光夫琅和费矩孔衍射实验的仿真程序 Bgjkfz. m。 %白光矩孔衍射仿真程序Bgjkfz. m clear; clf %读取色光的波长及其RGB值 Myldrgb; lam= lambda* 1e- 9; RGB= rgb; %设置仿真参数及图像矩阵 a= 5. 5e- 4; b= 4e- 4; f= 1;%方孔的长与宽及透镜焦距 [ x, y] = meshgr id( - 0. 003: 6e- 6: 0. 003); Irgb= zeros( 1001, 1001, 3) ; Iw= Irgb; %计算相对光强及仿真结果RGB矩阵 for k= 1 B1 B2 000 a lpha= ( p*i a* x) /( *f lam( k) );be ta= ( p*i b* y) /( *f lam( k) ) ; Ia= sin( alpha). ^2. /( a lpha. ^2+ eps) ; Ib= s in( be ta). ^2. / ( be ta. ^2+ eps); I= Ia. * Ib;%相对光强I / I0 Iw(: , : , 1) = I* RGB( k, 1); %计算Iw矩阵R维 Iw(: , : , 2) = I* RGB( k, 2); %计算Iw矩阵G维 Iw(: , : , 3) = I* RGB( k, 3); %计算Iw矩阵B维 Irgb= Irgb+ Iw; Iw= [ ]; %计算仿真结果Irgb矩阵 end %显示仿真结果 B= 1 /max(max(max( Irgb) ) ); II= Irgb* B; Br= 28; imshow( II* Br) 2)白光夫琅和费单缝衍射实验的计算机仿真当 by 0 时, 式( 1)变为(2)即为夫琅和费单缝衍射的光强 分布。 由此说明, 单缝衍射是矩孔衍射的特例。因而, 只须适当调整矩孔的长和宽, 利用程序Bgjkfz. m, 即可仿真白光夫琅和费单缝衍射实验。实验总结实验总结本次实验一方面利用了 Flash 动画原理对衍射实验现象进行了动态描述，另一方面讨论了利用计算机及 MATLAB 语言仿真夫琅禾费衍射实验，对单色光﹑复色光分别进行了单缝﹑双缝﹑多缝衍射的仿真模拟，最终得到的计算机模拟图像与实验现象一致，并且仿真结果能够随参数的变化而相应变化。而在实际的干涉与衍射等光学实验中, 各种参数很难调节,利用计算机进行仿真, 则简单易行。因此,利用计算机仿真光学实验, 可以弥补实际实验的不足。 关 键 词： 实验 报告 仿真 干涉 衍射

 天天文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。