利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
k'BLos1W . #;;pu7W 1、光疏射向光密
iM Xl}3 a]4 65FY clear
M#qZ0JT4 uj-q@IKe close all
n8M/Y}mH p5~;8Q7 n1=1,n2=1.45;
8P 3EQ Y - N3nk\ )V\E theta=0:0.1:90;
/y,~? 9zkR)C a=theta*pi/180;
K`@GN T& 15r <n rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
1!/
U#d" |400N
+MK rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Z<|ca T]Q( G2k r~FG tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
9b !+kJD A]{8 = ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
i)0*J?l= %xR;8IO figure(1)
>:s.` jV< fwI Zr~l subplot(1,2,1);
w<awCp zk4yh%Cd_ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Fa9]!bW ^AD/N|X^ legend('r_p','r_s','|r_p|','|r_s|')
q5 [%B K ;1cX|N= xlabel('\theta_i')
"$#x+|PyC #4''C s ylabel('Amplitude')
79a{Zwdd9j KTk% N p title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
e@ L'H)w, T.HI
$(d axis([0 90 -1 1])
mB#`{|1[ [6
!/ grid on
5RTAM o"v>
BhpC subplot(1,2,2);
D|Z, eench ;=6 ++Oq plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
y6; '?.Y1 '&3Sl?E legend('t_p','t_s','|t_p|','|t_s|')
jo<[|ZD `ivr$b# xlabel('\theta_i')
Uz H)fB [(`T*c.#.X ylabel('Amplitude')
.d9VV& i[^?24~ c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
DSy, #yA [ 8SW0wsk axis([0 90 0 1])
*Yvfp{B X<(h)&E grid on
4%p5X8|\ih TRAs5I% Rp=abs(rp).^2;
d+ih]? S53 [Ja Rs=abs(rs).^2;
uHQJ& f<WnPoV Rn=(Rp+Rs)/2;
Z[AJat@H Ajq;\- : Tp=1-Rp;
Y. i<7pBt ^=D77 jS Ts=1-Rs;
eJ%~6c`@! Z5-"a?{Y Tn=(Tp+Ts)/2;
S5v>WI^0h cWp
n/.a figure(2)
w_,. D Z ~|yH subplot(1,2,1);
ED"5y }.2pR*W plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
ERGDo=j =t&B8+6 legend('R_p','R_s','R_n')
$|6Le;
K HC4ad0Gs+{ xlabel('\theta_i')
cGsxfwD .(CzsupY_q ylabel('Amplitude')
I LF"m; ) Ah title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
?_W "=WpC ;csAhkf:S axis([0 90 0 1])
5&2=;?EO 5:CC\!&QBV grid on
Ej 'a
G A~nq4@uj subplot(1,2,2);
V[+ Pb] |m k $W$h plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
lUCdnp;w' N.xmHv Pk legend('T_p','T_s','T_n')
kc|`VB8L }
%S1OQC xlabel('\theta_i')
fBD 5K3 f>ED ylabel('Amplitude')
gA2\c5F< 8xg:ItJaA0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Mk|h ><Q" "H%TOk7l axis([0 90 0 1])
Q` AJR$L -Q 6W`*8 grid on
$CL =M xE;fM\7pu 
a2Q9tt>Q ] wa?~;1^& 
@ %ChPjN
