利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
sUfYEVjr {x+"Ru~7, 1、光疏射向光密
L+c7.l.yT !F |mCEU clear
;}tEU'& #9{9T"ed close all
vSt7&ec J<
E"ZoY n1=1,n2=1.45;
[Ipg",Su;f Tk^J#};N theta=0:0.1:90;
P~>E {E oRY/] a=theta*pi/180;
UogkQ& B :N8 26_q rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
C/JeD-JG H9x,C/r, rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
N 34.Bt Y=%SK8]Q; tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
D*>EWlZ aX oD{zA ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Q Eh_2 \WDL?(G< figure(1)
y U-^w^4 )wmG&"qsP subplot(1,2,1);
m#D+Yh/y{n d ,Fj|}S plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
' THcO*< uN:|4/;{& legend('r_p','r_s','|r_p|','|r_s|')
Wz}8O]#/. .H M1c xlabel('\theta_i')
uB(16|W>S OX"^a $ ylabel('Amplitude')
hnQDm$k J3]W2m2Zw title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6I$laHx? 9@Iz:!oqb axis([0 90 -1 1])
>q'xW=Y
j\ _f~(g1sE grid on
$`2rtF +<G |Ru- subplot(1,2,2);
;g3z?Uz) 19Cs
3B \4 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
@R5jUPUVV Bf72 .gx{0 legend('t_p','t_s','|t_p|','|t_s|')
pJ` M5pF B! `Dj, _ xlabel('\theta_i')
hi{#HX a yGNZw7^( ylabel('Amplitude')
K3jPTAw=# Ub0hISA title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/Hox]r]'e y:U'3G- axis([0 90 0 1])
(,5oqU9s@ r/X4Hy0!lT grid on
w/ 8`]q uHBEpqC % Rp=abs(rp).^2;
K[wOK DCJmk6 p%0 Rs=abs(rs).^2;
z (N3oBW E8TJ*ZU Rn=(Rp+Rs)/2;
+`EF0sux `EV"
/&` Tp=1-Rp;
yI&{8DCCw |- WoR u Ts=1-Rs;
]L'FYOfrpx dQoZh E Tn=(Tp+Ts)/2;
\9U4V>p 9;Z 2.P"w figure(2)
}PZz(Ms @%4MFc0`! subplot(1,2,1);
M*DF tp< \JJ > y plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
b-/zt Z@u =j)y. x( legend('R_p','R_s','R_n')
T6X%.tR>` j3Ng] @N xlabel('\theta_i')
] gN]Cw\L wzw`9^B ylabel('Amplitude')
64Ot`=A" Hp o/CY/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]dX HjOpA omxBd#;F$ axis([0 90 0 1])
A),nkw0X 2<d l23 grid on
br!:g]Vh r ]XXN2[jO subplot(1,2,2);
T5m dC ; Nw . plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
mhp&;
Q9 } 3t bqFiH legend('T_p','T_s','T_n')
?/mk FDN ryz
[A:^G xlabel('\theta_i')
O"otzla DVu_KT[H d ylabel('Amplitude')
\z}/=Qgc m oQ><>/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^y .e
Fz 9m2_zfO[ w axis([0 90 0 1])
>".,=u ' )ca^%(25!z grid on
E"j u<q/Q 3Lxk7D>0c 
@PXb^x#k ap;tggi(H 
zj!&12w%3
