基于MATLAB的光学实验仿真研究和交互式界面设计

《基于MATLAB的光学实验仿真研究和交互式界面设计》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、西安工业大学北方信息工程学院毕业设计（论文）中期报告题目：基于MATLAB的光学实验仿真研究与交互式界面设计系别专业班级代写论文姓名18192322624学号导师设计(论文)进展状况1.编写杨氏双缝干涉仿真程序clear%基本的初始化设置lam=580e-9%单色光波长a=2e-3;%缝宽D=l;%缝到屏的距离ym=5*lam*D/a;xs=ym;n=101;ys=linspace(-ym,ym,n);fori=l:nrl=sqrt((ys(i)-a/2).A2+DA2);r2=sqrt((ys(i)+a/2

2、).A2+DA2);phi=2*pi*(r2-r1)./lam;B(i,:)=sum(4*cos(phi/2).A2);%函数End%仿真图样的显示N=255;Br=(B/4.0)*Nsubplot(121)%显示比例image(xs,ys,Br);%图样colormap(gray(N));%对应图的曲线subplot(l22)%显示比例plot(B,ys)运行的得到的图像如图1：图1杨氏双缝干涉仿真模拟图样1.编写等厚干涉仿真程序Clearformatlong%参数设置N=301;R=5;%曲率半径lam=

3、5893e-10;%单色光波长fork=0:l:60xm=0.00300+k*2e-5;ym=0.00300+k*2e-5;x=linspace(-xm,xm,N);y=linspace(-xm,xm,N);fori=l:Nforj=l:Nr(i,j)=sqrt(x(i)*x(i)+y(j)*y(j));B(i,j)=4*cos(pi*(r(i,j)A2/R+lam/2)/lam).A2;endendi=(N+l)/2;forj=l:NBB(i,j)=4*cos(pi*(r(i,j)A2/R+lam/2)/l

4、am).A2;End%图样显示nclevels=255;Br=0.25*B*nclevels;subplot(2,l,l);image(x,y,Br)；axis(卜xm,xm,-ym,ym]);colormap(gray(255));subplot(2,l,2);axis(卜xm，xm,-4,4]);plot(y,BB);N=N+2;pauseend运行得到的图像如图2：-3X10图2牛顿坏干涉仿真模拟图样1.编写等倾干涉仿真程序Clear%参数设置lambda=500*10*9);%单色光的波长f=0.2;

5、%透镜焦距(1=2.5*10A(-4);%薄膜厚度theta=0.15;%角度rMax=f*tan(theta/2);N=501;fori=l:Nx(i)=(i-l)*2*rMax/(N・1)-rMax;forj=l:Ny(j)=(j-1)*2*rMax/(N-1)-rMax;r(i,j)=sqrt(x(i)A2+y(j)A2);delta(i,j)=2*d/sqrt(1+r(i,j)A2/fA2);Phi(i,j)二2*pi*delta(i,j)/lambda;B(ij)=4*cos(Phi(ij)/2)A

6、2;endend%图样显示NCLevels=255;Br=(B/4.0)*NCLevels;figure(1);image(x,y,Br)；colormap(gray(NCLevels));axissquare;figure(2);mesh(x,y,Br)；colormap(gray(NCLevels));axissquare;运行得到的图像如图3：-0015・001-0005000050010015-0015-0014)005000050010015图3边克尔逊干涉仿真模拟图样1.编写多光束干涉的仿真程序c

7、learlam二[600610570550460440410]*le-5;%七色光的波长RGB=[1,0,0;1,0.5,0;1,1,0;0,1,0;0,1,1;0,0,1;0.67,0,1];%七色光的三原色光模式(1=0.05;%两镀银平面间隔距离pc=zeros(500,500,3)；%设置光屏%可调参数「0.54;%反射系数n=1.5;%折射率D二1.7;%薄膜厚度[x,y]=meshgrid(linspace(-0.005,0.005,500));z二x+li*y;theta=atan(abs(z)

8、/d);%入射角设置fork=l:7phi=4*pi*n*D*cos(theta)/lam(k);%相位差It=17(1+4*rA2*sin(phi/2).A2/(l-rA2)A2);%光强p(:,:,l)=It*RGB(k,l);p(:,:,2)=It*RGB(k,2);p(:,:,3)二It*RGB(k,3);%将包含颜色信息的光强用矩阵的形式保存起来pc=pc+p;p=[];end;%显示Br