物理光學(xué)matlab第二次小作業(yè)

1、電磁波界面反射1.要求 2.數(shù)學(xué)模型: 3仿真:3.1 光強(qiáng)曲線 3.2 偏振方向變化曲線: 3.3偏振光束相位變化曲線: 結(jié)果分析:光強(qiáng)曲線描述了光強(qiáng)隨入射角度變化的特點(diǎn),第二個(gè)圖象描繪了偏振方向隨角度變化的特點(diǎn),而最后的圖中第一個(gè)圖為p偏振光束相位變化,第二個(gè)圖為s偏振光束相位變化紅線表示的是入射角小于臨界角時(shí)偏振光束相位變化,藍(lán)線表示入射角大于臨界角時(shí)偏振光束相位變化?？偨Y(jié)和感悟:通過(guò)對(duì)電磁波界面反射的理解和應(yīng)用,我對(duì)物理光學(xué)有了更深的認(rèn)識(shí)和體會(huì),對(duì)matlab也有了進(jìn)一步的了解,當(dāng)真是受益匪淺。附件：matlab代碼 1. a1=0:0.1:70; a2=asin(1.54*sin(d

2、eg2rad(a1; rp=-tan(deg2rad(a1-a2./tan(deg2rad(a1+a2; rs=-sin(deg2rad(a1-a2./sin(deg2rad(a1+a2; I=2.5*(abs(rp.2+abs(rs.2; plot(a1,I; 2. 代碼： a1=0:0.1:70; a2=asin(1.54*sin(deg2rad(a1; rp=-tan(deg2rad(a1-a2./tan(deg2rad(a1+a2; rs=-sin(deg2rad(a1-a2./sin(deg2rad(a1+a2; pian=rad2deg(atan(rs./rp; L=length(

3、rp; for i=1:1:L if imag(rp(i=0 j=i; end end X=zeros(1,j; Y=zeros(1,j; for n=1:1:j; X(n=a1(n; Y(n=pian(n; end plot(X,Y; 3. 代碼： a1=0:0.1:70; a2=asin(1.54*sin(deg2rad(a1; rp=-tan(deg2rad(a1-a2./tan(deg2rad(a1+a2; rs=-sin(deg2rad(a1-a2./sin(deg2rad(a1+a2; pian=rad2deg(atan(rs./rp; L=length(rp; for pian=

4、1:1:L if imag(rp(i=0 j=i end end Phip1=zeros(1,j; Phip2=zeros(1,j; Phis1=zeros(1,j; Phip3=zeros(1,L-j; Phip4=zeros(1,L-j; Phis2=zeros(1,L-j; for h=1:1:L-j Phip3(h=a1(h+j; end for h=1:1:L-j Phip4(h=rad2deg(angle(rp(h+j; Phis2(h=rad2deg(angle(rs(h+j; end for h=1:1:j Phip1(h=a1(h; if rp(h0 Phip2(h=180; end if rs(h0 Phis1(h=180; end end subplot(2,1,1; plot(Phip1,Phip2,ro; hold on; plot(Phip3,Phip4; hold off; subp