利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
L< 3U)Gp u~kf z*hz 1、光疏射向光密
!YJfP@"e6r <HReh>)[ clear
46U?aHKW@| (N6 3k1M close all
*=+m;%]_ $"V gN ynq n1=1,n2=1.45;
_ ," -25a 'rWu}#Nb theta=0:0.1:90;
VU~
R Grot3a a=theta*pi/180;
k GaK(^w "'389*- rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
aI8k:FK" Z' cQ<
f rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
]#)1(ZE ARcPHV<(2 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
bwH l}3 ED9uKp<Wbv ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
n%N|?!rB hN(sz figure(1)
h( lkC[a& 6Xu^ cbD subplot(1,2,1);
j W a% vA &hciv\YT2W plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
g~z z[F 8U qx#k()E.U legend('r_p','r_s','|r_p|','|r_s|')
>FrF"u:kM &c;@u?:@S xlabel('\theta_i')
eVRFb#EU0e h>s|MZQ:* ylabel('Amplitude')
m(~5X0 }zA
kUt title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
# X~{p4Lr [A@K)A$f axis([0 90 -1 1])
hXxgKi% |~QHCg< grid on
UkO L 7M /^#8z(@B subplot(1,2,2);
.=y-T=} ;
E Nhy plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
*k&yD3br-V H
l'za legend('t_p','t_s','|t_p|','|t_s|')
N$Pi4 ifo^
M]v xlabel('\theta_i')
u!N Y@$Wc ~d+.w%Z ` ylabel('Amplitude')
yr p;G _ 1e Wl:S} title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
AsB ep SV-M8Im73z axis([0 90 0 1])
6 fP"I_c PS*=MyNa grid on
2(_+PQ6C= p&Os 5zw;| Rp=abs(rp).^2;
'Q R
@G B vXA9YQ3 Rs=abs(rs).^2;
N>Dr
z bnso+c A Rn=(Rp+Rs)/2;
FiN^}Kh *'b3Z3c,; Tp=1-Rp;
:>@6\ $d<vPpJ3 Ts=1-Rs;
80i-)a\n Y)X
'hk)5| Tn=(Tp+Ts)/2;
iX 3Y:
~Y-
!PZ figure(2)
\S] "nHX B@v\t pR subplot(1,2,1);
AF d 3_>h =qQH,{]c6 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
{km~,]N QkEIV<T&)l legend('R_p','R_s','R_n')
'N- nFc^ 5>E]C=maD xlabel('\theta_i')
8T:?C~ " zz4TJ (' ylabel('Amplitude')
z>\vYR$ \Lp|S:u title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
>8I?YT. ~EYsUC#B_ axis([0 90 0 1])
!B&OK&* 7Wd}H Z grid on
QD<GXPu?N *]L(,_:" subplot(1,2,2);
;WF3w N U>'$s plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
j. @CB` Ya%-/u legend('T_p','T_s','T_n')
[Pn(d [$z /7s^OkQ xlabel('\theta_i')
+#||
w9p ,Z2fVz~9 ylabel('Amplitude')
t`b !3U>I 5Op |="W. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:\]TA Qd- fyRSg B00$ axis([0 90 0 1])
-E*VF{IG1 q y]tuKZI grid on
K
*<+K<Tp 2j9+ f{ l 
3RwDI k?>% ,*y\b|<j 
z5V~m_RO
