利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
p<*3mbgGO u;p{&\(] 1、光疏射向光密
Z;=G5O
uvQ kB { clear
b[< r+e8 P%
_cIR close all
)1wC].RFYm Nl ,M9 n1=1,n2=1.45;
i-w^pv' Q{T6t;eH theta=0:0.1:90;
j}3Avu% #g*U\y a=theta*pi/180;
]l7) F-v C']TO/2q rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Q;W[$yvW %^ f! = * rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
KWbnSL8 |1wfLJ4--l tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L|.q19b* $$ou qLu ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Rd@34"O i z!E1(z( figure(1)
l{ja2brX cih@: =Qy subplot(1,2,1);
s S(t
}$ /~[+' plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Dn~r~aR$g =NLsT.aa legend('r_p','r_s','|r_p|','|r_s|')
huFT_z_;; ~sj'GEhEg xlabel('\theta_i')
e f8_w6i / RU'~( ylabel('Amplitude')
yW 5 /Y02 07 [%RG title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)Be?axI 5Z"IM8? axis([0 90 -1 1])
I,;@\ VWDXEa9 grid on
X^;[X~g 3SI:su subplot(1,2,2);
"z Fv? ay Cq\1t plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
8w 1TX [ b !lQGoXQ'4 legend('t_p','t_s','|t_p|','|t_s|')
&c!d}pU} ZI. ;7G@| xlabel('\theta_i')
7{F(NJUO1 b-4g HW ylabel('Amplitude')
0 kC}qru' 5@l[!Jl0k title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:'1UX <&B |m{Q_zAB axis([0 90 0 1])
Z9P rw/8P }?[^q grid on
I#lvaoeN T}')QC&wQ Rp=abs(rp).^2;
s&(,_34 m BW E^ Rs=abs(rs).^2;
nADt8 [4qCW{x._ Rn=(Rp+Rs)/2;
>! u@> xW|^2k Tp=1-Rp;
~{69&T}9 }{FKs!(4 Ts=1-Rs;
CbPCj.MH ]MI> "hn Tn=(Tp+Ts)/2;
Twqkd8[ 'zYKG5A figure(2)
Bfh[C]yy dzA5l:5 subplot(1,2,1);
6 0C;J!D -anLp8G* plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
OPm ?kr "F_o%!l legend('R_p','R_s','R_n')
g~q+a- /JP]5M) xlabel('\theta_i')
x0t&hY >P! \Xy]z ylabel('Amplitude')
do3 BI4Q D+$ k title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ag Q5%t# }qf9ra axis([0 90 0 1])
$^&SEz 'Na \9b ( grid on
XD1 x*# /t "p^9!^ subplot(1,2,2);
M :m-i X n])-+[F plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
sOyWsXd+R' @( :ah legend('T_p','T_s','T_n')
o;#{N~4[$ e"jA#Y # xlabel('\theta_i')
Y~RZf /` HSt|Ua.c/h ylabel('Amplitude')
#
Sm M5% +{%@kX<V_ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[\9WqHs [V#"7O vl axis([0 90 0 1])
Oto pA) 4fDo }~ grid on
;bt@wgY cQ4TYr;? 
Ru%:
z>Y =hV-E
D 
Z=: oI Ae
