光课程设计光波在介质中界面上的反射及透射特性的仿真

1西安邮电大学光学报告学院：电子工程学生姓名：专业名称：光信息科学与技术班级：光信1103班2设计名称：光波在介质中界面上的反射及透射特性的仿真一、课程设计目的1．掌握反射系数及透射系数的概念；2．掌握反射光与透射光振幅和相位的变化规律；3．掌握布儒斯特角和全反射临界角的概念。二、任务与要求对n1=1、n2=1.52及n1=1.52、n2=1的两种情况下，分别计算反射光与透射光振幅和相位的变化，绘出变化曲线并总结规律。三、课程设计原理光在介质界面上的反射和折射特性与电矢量的振动方向密切相关。由于平面光波的横波特性，电矢量可在垂直传播方向的平面内的任意方向上振动，而它总可以分解成垂直于入射面振动的分量和平行于入射面振动的分量，一旦这两个分量的反射、折射特性确定，则任意方向上的振动的光的反射、折射特性也即确定。菲涅耳公式就是确定这两个振动分量反射、折射特性的定量关系式。（1）s分量和p分量3psmEEtEErimtmmimrmm,,,0000垂直入射面的振动分量--s分量平行入射面的振动分量--p分量定义：s分量、p分量的反射系数、透射系数分别为（2）反射系数和透射系数定义：s分量、p分量的反射系数、透射系数分别为psmEEtEErimtmmimrmm,,,0000（3）菲涅耳公式已知界面两侧的折射率21nn、和入射角1，就可由折射定律确定折射角2；进而可由菲涅耳公式求出反射系数和透射系数。绘出如下按光学玻璃（n=1.5）和空气界面计算，在21nn（光由光疏介质射向光密介质）和21nn（光由光密介质射向光疏介质）两种情况下，反射系数、透射系数随入射角1的变化曲线。4(a)光由光疏介质射向光密介质(b)光由光密介质射向光疏介反射光与入射光中s,p分量的相位关系：（1）n1＜n2时，光疏入射光密s分量的反射系数sr：反射光中的s分量与入射光中的s分量相位相反；反射光中的s分量相对入射光中的s分量存在一个相位突变(rs=)；p分量的反射系数pr：在1B范围内,pr＞0，反射光中的p分p量与入射光中的分量相位相同(rp=0)；在1B范围内,pr0，反射光中的p分量相对入射光中的p分量有π相位突变(rp=)；（2）n1n2时，光密入射光疏s分量的反射系数sr：入射角1在0到C(临界角，12/sinnnC）的范围内，s分量的反射系数sr＞0。反射光中的s分量与入射光中的s分量同相位，rs=0；入射角1C时，发生全反射，1212cossin2tannrs（21/nnn）；p分量的反射系数pr：在1＜B范围内，pr0，反射光中的p分量相对入射光中的p分量有π相位突5变（rp=)；在B1C范围内，pr0，反射光中的p分量与入射光中的p分量相位相同(rp=0)；入射角1C时，发生全反射，12122cossin2tannnrp；四、课程设计步骤（流程图）6（1）定义变量n1,n2,f1.（2）给变量赋值,其中n1=1,n2=1.52,还有一种情况其中n1=1.52,n2=1（3）设计for循环，使f1每循环一次加/1000，实现在f1每变化一次下，得出相应的反射系数，透射系数的值，从而得出程序的循环（4）根据程序仿真结果五、仿真结果分析0102030405060708090-1-0.500.51Fn1n2s/p分量与相位的关系010203040506070809001234Ffrs010203040506070809001234Ffrp0102030405060708090-101234n1n2s/p分量与相位的关系F010203040506070809001234Ffrs010203040506070809001234Ffrp结论：光在介质面上的反射、透射特性有三个因素决定：入射光的偏振态，入射角，界面两侧介质的折射率。(1)光波由光疏介质射向光密介质（n1n2）a.n1n2时，反射系数rs0,说明反射光中的s分量与入射光中的s分量相位相反。7(即frs=)b.而p分量的反射系数rp在f1fb范围内，rp0,说明反射光中的p分量与入射光中的p分量相位相同。（即frp=0）c.在f1fb范围内，rp0,说明反射光中的p分量与入射光中的p分量相位突变。（即frp=）(2)光波由光密介质射向光疏介质（n1n2）a.入射角f1在0到fc的范围内，s分量的反射系数rs0,说明反射光中s分量与入射光中的s分量同相位。（即frs=0）b.P分量的反射系数rp在f1fb范围内，rp0,说明反射光中的p分量相对入射光中的p分量有相位突变。（即frp=）c.在fbf1fc范围内，rp0,说明反射光中的p分量与入射光中的p分量相位相同。六、仿真小结光在介质界面上的反射、透射特性由三个因素决定：（1）入射光的偏振态；（2）入射角;（3）界面两侧介质的折射率。由rs、rp、ts、tp随入射角的变化曲线可知，在入射角从0度到90度的变化范围内，不论光波以什么角度入射至界面，也不论界面两侧折射率大小如何，s分量和p分量的透射系数t总是取正值，因此，折射光总是与入射光同相位。通过本次实验，掌握了反射系数及透射系数的概念，反射光与透射光振幅和相位的变化规律，布儒斯特角和全反射临界角的概念。8七、程序clearall;%n1=1;%n2=1.52;n1=1.52;n2=1;n=n2./n1;ifn1n2subplot(1,3,1)qa=0:pi/100:pi/2;qb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);rp=tan(qa-qb)./tan(qa+qb);ts=2.*cos(qa).*sin(qb)./sin(qa+qb);tp=2.*cos(qa).*sin(qb)./sin(qa+qb)./cos(qa-qb);plot(qa*180./pi,rs,'r',qa*180./pi,rp,'c',qa*180./pi,ts,'b',qa*180./pi,tp,'g')legend('rs','rp','ts','tp')%rssubplot(1,3,2)forqa=0:pi/1000:pi/2qb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);ifrs=0Frs=pi;elseFrs=0;endplot(qa*180./pi,Frs,'r')holdonendlegend('Frs')%rpsubplot(1,3,3)forqa=0:pi/1000:pi/2qb=asin(n1.*sin(qa)./n2);rp=tan(qa-qb)./tan(qa+qb);ifrp=0Frp=pi;elseFrp=0;9endplot(qa*180./pi,Frp,'b')holdonendlegend('Frp')elsesubplot(1,3,1)qc=asin(n2./n1);qa=0:0.0001:qc;qb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);rp=tan(qa-qb)./tan(qa+qb);ts=2.*cos(qa).*sin(qb)./sin(qa+qb);tp=2.*cos(qa).*sin(qb)./sin(qa+qb)./cos(qa-qb);plot(qa*180./pi,rs,'r',qa*180./pi,rp,'c',qa*180./pi,ts,'b',qa*180./pi,tp,'g')holdonqa=qc:0.0001:pi/2;tp=0;ts=0;rs=1;rp=1;plot(qa*180./pi,rs,'r',qa*180./pi,rp,'c',qa*180./pi,ts,'b',qa*180./pi,tp,'g')holdonlegend('rs','rp','ts','tp')%rsqc=asin(n2./n1);subplot(1,3,2)forqa=0:pi/1000:qcqb=asin(n1.*sin(qa)./n2);rs=-sin(qa-qb)./sin(qa+qb);ifrs=0Frs=pi;elseFrs=0;endplot(qa*180./pi,Frs,'r')holdonendqa=qc:pi/1000:pi/2;10Frs=2.*atan(sqrt(sin(qa).^2-(n.^2))./cos(qa));plot(qa*180./pi,Frs,'r')holdonlegend('Frs')%rpsubplot(1,3,3)forqa=0:pi/1000:qc;qb=asin(n1.*sin(qa)./n2);rp=tan(qa-qb)./tan(qa+qb);ifrp=0Frp=pi;elseFrp=0;endplot(qa*180./pi,Frp,'b')holdonendqa=qc:pi/1000:pi/2;Frp=2.*atan(sqrt(sin(qa).^2-(n.^2))./cos(qa)./n.^2);plot(qa*180./pi,Frp,'b')holdonlegend('Frp')end