利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
vUqe.?5 @`5QG2 1、光疏射向光密
28 ja-1dB :MK:T JV clear
d3og?{i<}& 4{Q$!O> close all
lK "' n LL J`4Z<b53 n1=1,n2=1.45;
?n2C Wa (W&] theta=0:0.1:90;
FMoJ"6Q Fs?( UM a=theta*pi/180;
L^6"' # 6BUBk>A` rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
7F"3 <U@J e jP,29 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
&Q2N U$ Um/CR! tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
uPI v/&HA 1aBD^^Y ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
lQ+Ru8I >s;oOo+5 figure(1)
2j2 mW>Z :8](&B68gE subplot(1,2,1);
Qk?jG XB>^ &,^mM'
C plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
z`c%?_E K G%xb0%oi]% legend('r_p','r_s','|r_p|','|r_s|')
KX0<j AmP#'U5 xlabel('\theta_i')
-L+\y\F PthId aN@ ylabel('Amplitude')
oxUE79 ,Y`'myL8W title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
! x`;>0 );0<Odw%. axis([0 90 -1 1])
_{N0OX *!._Ais,\ grid on
Q&oC]u(="& RpBiE8F4 subplot(1,2,2);
lc\>DH\n6 \;nD)<)J plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
A ~qW. ~5lKL5w legend('t_p','t_s','|t_p|','|t_s|')
L!CX &
f$:7A0 xlabel('\theta_i')
kWz uz# jA20c(O ylabel('Amplitude')
L9E;Uii0 ZA *b9W title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
i_j9/k 4M;sD;3 axis([0 90 0 1])
ovhC4 2i \>*M Me grid on
}* f W!(* >i^y;5 Rp=abs(rp).^2;
:BewH?Ku t_[M & Rs=abs(rs).^2;
MgJ%26TZ 3q~":bpAp Rn=(Rp+Rs)/2;
( {;P#qCX fe&
t- Tp=1-Rp;
!Z<mrr;T@ @ A~B
, Ts=1-Rs;
`^ieT#(O Ns+)Y^(5 Tn=(Tp+Ts)/2;
[
4?cM\_u@ )@8'k]Glw. figure(2)
Ze$^ UR Hrd 5p+j subplot(1,2,1);
E5 ,%J `^X RrVX< plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
n9<roH Z^/z legend('R_p','R_s','R_n')
"4Vi=* 2V R`Ys;g/! xlabel('\theta_i')
:wv
:#EaH jQ .>2-;H9 ylabel('Amplitude')
OCbQB5k3 %fn'iKCB title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
PT uCN c t2_N axis([0 90 0 1])
"W?l R4 bl:a&<F grid on
7' :5
)o
" SB1 subplot(1,2,2);
h@PMCmf_ C,OB3 y plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
ohQz%?r r
48;_4d)D legend('T_p','T_s','T_n')
+H)!uLva B MYDf`0{$_a xlabel('\theta_i')
V:,3OLL* F'?5V0\he ylabel('Amplitude')
t K s4}vW ,z0~mN title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!~|"LA!jn $.r}g\43P axis([0 90 0 1])
oSu|Y n Gpxb_}P grid on
<M =W)2D7 o3TBRn, 
\0* yxSg,^ bY"eC i{K 
5JU(@}Db
