电磁波界面反射特性理解+仿真分析

1、电磁波界面反射特性理解+仿真分析要求一束5W的线偏振光以* = 45度方位角振动，垂直入射到玻璃-空气表面，该光束波长0.6 m玻璃介质折射率1.54 0.6 m0-70度变化时，通过给定条件，分别完成如下要求：1建立反射光强 I （n,仇中）的数学模型；2画出该光束反射光的光强曲线（I -0 ）;3分析该反射光束的偏振方向或者偏振态变化，画出偏振方向变化曲线仲-日）和偏振光束相位变化曲线（*-6（注：光束从光密到光疏的界面，在入射角日从0-70度变化中，包括了临界角ec）.反射光的光强曲线（I - e ）数学模型:将这三种光波的电矢量振动方向都分解成两个分量，一个垂直于入射面，称为垂直分量s

2、;另一个平行于入射面，称为平行分量 p,这两个分量互相垂直。而任何偏振光都可以分解为互相垂直的两个分量，可以得出反射光强公式:I = Is +1 p, % = Ais/ Ais,rp = Ap/Aip2,2222I =A ,. I =& *Ais ”* Ap菲涅尔公式有：=-则一二21rp =则上=_)；因为偏振光以 =45度方位角振动，所以sin- 22tan 口 , 士As = Ap = JT/2,折射定律有n1sine =n2sin为，由此可以求出反射光的光强曲线(I-8仿真：反射光强曲残6 5 4 3 2 10-胆泥S吗10 2C3040入射角引60分析：入射角从0。增加，刚开始大部分

3、入射光发生折射，少数入射光发生反射，所以光强值很小，随着入 射角的增加，在接近临界角时大部分光发生反射，少部分光发生折射，此时反射光强快速增加，当入射角 大于临界角后发生全反射，反射光强与入射光强相等。.偏振方向变化曲线(*-0 )数学模型:反射光有-e-n匕，其中.为偏振方向，因为rs = As/As,p = Ap/Ap且As = Ap,可得tane=%/rp。可以得到偏振方向变化曲线的-8 )1Owo8。5040200建o曜-2040-60-SO,100051015202530354046入射角一分析：入射角大于临界角后为椭圆偏振光，没有相应的偏振方向。3.振光束相位变化曲线 仲-日）数学

4、模型：当入射角小于临界角时，rs，rP的正负值可以反应偏振光束的相位变化，例如rs为负值时，说明As和As的相位差相差180。当入射角大于临界角时，rs,rP变为复数，在全反射下 rs,rP有新的表达式：cos ? - j Esin2 口1 - n2n2 cosi - j - sin2Ui-n2rs =.2 .21rp =2.2 .2cos j . sin -1 -n n cos j . sin 11 - n由此可以求出偏振光束相位变化曲线-口p光相位变化曲线180iii r i16Q-140i2000 14020一20分析：入射角从0开始增加，rp开始为正，相位相差为 0。；当rp为负数时，

5、相位相差为 180。当入 射角大于临界角后，rp变为复数，其幅角为 p光相位变化。G光相位变化曲糕102030入射角厅rs为正数，此时相位相差 0。当入射角大分析：入射角从0开始增加，在入射角小于临界角的范围内, 于临界角时，rs为复数，其幅角为s光相位变化。总结通过此次matlab的仿真学习，并结合书本的知识我充分认识了电磁波界面反射中入射光和反射光的光强变化，以及偏振方向和偏振相位的变化关系。将偏振光分解为两个互相垂直的分量s光和p光，只要把s和p分量的振幅和相位关系弄清楚，任何偏振光的振幅和相位关系就可以完全确定。Matlab代码.光强曲线(I -8 )n=1.54;theta1=0:0

6、.1:70;I=5;theta2=asind(sind(theta1)*n);rs=-sind(theta1-theta2)./sind(theta1+theta2);rp=-tand(theta1-theta2)./tand(theta1+theta2);I1=0.5*I*(abs(rs).A2+abs(rp).A2);plot(theta1,I1,r);title( 反射光强曲线)；ylabel(反射光强I);xlabel(入射角 0 );.偏振方向变化曲线的-日)n=1.54;theta1=0:0.1:70;I=5;theta2=asind(sind(theta1)*n);rs=-sin

7、d(theta1-theta2)./sind(theta1+theta2);rp=-tand(theta1-theta2)./tand(theta1+theta2);phi=atand(rs./rp);n=length(rp);for i=1:1:nif imag(rp(i)=0%取入射角小于临界角的值m=i;endendx=zeros(1,m);phi1=zeros(1,m);for j=1:1:mx(j)=theta1(j);phi1(j)=phi(j);endplot（x,phi1）;title（偏振方向变化曲线）；ylabel（偏振角 4 ）;xlabel（入射角 0 ）;.偏振光束相

8、位变化曲线（1-日） theta1=0:0,1:70;theta2=asind(sind(theta1)*1.54);rs=-sind(theta1-theta2)./sind(theta1+theta2);rp=-tand(theta1-theta2)./tand(theta1+theta2);phi=atand(rs./rp);n=length(rp);for i=1:1:nif imag(rp(i)=0 m=i;endendx1=zeros(1,m);phi_p1=zeros(1,m);phi_s1=zeros(1,m);x2=zeros(1,n-m);phi_p2=zeros(1,n-

9、m); phi_s2=zeros(1,n-m);for j=1:1:m x1()=theta1();if rp(j)0 phi_p1(j)=180;endif rs(j)0 phi_s1(j)=180;endendfor h=1:1:n-m x2(h)=theta1(m+h);%取入射角小于临界角的值%右界角前相位变化%临界角后相位变化endfor k=1:1:n-mphi_p2(k)=rad2deg(angle(rp(m+k);phi_s2(k)=rad2deg(angle(rs(m+k);end plot(x1,phi_p1, r);hold on plot(x2,phi_p2);hold