利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
b-ZvEDCR R,t R{| 8 1、光疏射向光密
NK d8XQ=% Jln dypE clear
_w(SHWh2 37{mhU close all
3EAu#c@q" '*W/Bett n1=1,n2=1.45;
! k&< [t,7H theta=0:0.1:90;
Nmp1[/{J VTD'D+ t a=theta*pi/180;
i5|!M IY 2TY|)ltsF rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
xez~Yw2 {V8 v
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
voCQ_~*)9 eFbr1IV tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
:tNH Cx m ?LOd9 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
NfKi,^O O']-<E`1k figure(1)
2,'m]`;GNr K6_{AuL}4 subplot(1,2,1);
~- JkuRJ\ `AQv\@wp plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
t <x0?vfD *cJ GrLC legend('r_p','r_s','|r_p|','|r_s|')
>Xh 9{/o 3-srt^>w* xlabel('\theta_i')
;L++H5Kz6 ^#p S u ylabel('Amplitude')
sZ7{_}B !bS:!Il9= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
V # =o< 4)i/B9 9k axis([0 90 -1 1])
rl[&s\[ g&#.zJ[- grid on
K6{{\r ^. M*pe subplot(1,2,2);
vEOoG>'Zq >kd &>)9v plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
&Nt4dp`qj *h$Z:p-g legend('t_p','t_s','|t_p|','|t_s|')
-(ABQgSO] CFkM}`v0 xlabel('\theta_i')
sJ{NbN~`I
]8<;,} # ylabel('Amplitude')
1|Us"GQ (n =4 36/O`K title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Y FL9Q< 7lqj" o( axis([0 90 0 1])
t[Xx LG* "(> P= grid on
:"#EQq]ct lrPiaSO`I Rp=abs(rp).^2;
5\A[ra '!I^Lfz-Z Rs=abs(rs).^2;
!g-|@ W 4 jfkCU Rn=(Rp+Rs)/2;
H~^am }PTYNidlR Tp=1-Rp;
LkD$\i !U/iY%NE Ts=1-Rs;
6Xu8~%i )4vZIU # Tn=(Tp+Ts)/2;
y+hC !- y'(l]F1] figure(2)
h/j+ b.|
lU`]yL subplot(1,2,1);
!ZPaU11 mFC0f?nr plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
w O;\,zU T \5 5uQ legend('R_p','R_s','R_n')
KWq7M8mq 4W
&HUQ?^ xlabel('\theta_i')
N guJ[ N: [22`NP ylabel('Amplitude')
id="\12Bw ] ([:"j title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
k>I [U}h [#S}L(
axis([0 90 0 1])
[4KW64 %l g .ty#Z=: grid on
-
|n\
"E =\Vz subplot(1,2,2);
Bvj-LT=) r<,W{Va plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
2iYf)M C TO7 %TW{L legend('T_p','T_s','T_n')
B:Ec(USe 5RCZv\Wd& xlabel('\theta_i')
Hg\+:}k&9 :oh(M|;/2 ylabel('Amplitude')
6m"_=.k% =X6WK7^0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
t2d _XQOK {KYbsD axis([0 90 0 1])
Y>LgpO. zY6{ OP!# grid on
a "uO0LOb @| ! 9~F 
(,<&H;,8 (jv!q@@2C. 
iLiEh2%P
