利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
5S8>y7knQ SXP(C^?C 1、光疏射向光密
x|yJCs> < =~=IZ) clear
c9G%;U) (3K3)0fy close all
N,Z*d Z.&/,UU:4 n1=1,n2=1.45;
W+GC3W /7])] vZ_ theta=0:0.1:90;
E}KGZSj ilde<!? a=theta*pi/180;
m 94PFD@N <5D4h! rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
n807?FORB 1) ^\R(l rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
vQ* RrHG?c G^tazAEfo tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
P
JATRJ1. xxyc^\$ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
PDPK|FU m5iCvOP figure(1)
wgIm{;T[u {f\wIZ-K A subplot(1,2,1);
#2s}s<Sc; ;-8 .~Sm plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
JH{/0x#+ V]; i$ legend('r_p','r_s','|r_p|','|r_s|')
1Xo0(*O '5 Yzo^R; xlabel('\theta_i')
-N% V 5 TN &BLCP d ylabel('Amplitude')
>G8I X^*sG {ibu 0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h$kz3r;b," *ma
w`1 axis([0 90 -1 1])
]}PXN1( X5YOxMq grid on
:Rb\Ca
Nd RcA subplot(1,2,2);
i_Hm?Bi!F triU^uvh plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
e ,epKt L Et!J*{ s legend('t_p','t_s','|t_p|','|t_s|')
jQ;/ =9 cN0
*< xlabel('\theta_i')
:Bmn<2[Y; ttUK~%wSx ylabel('Amplitude')
\894 Jqh Sc]K-]1(H title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3s/1\m% ",3v%$ > axis([0 90 0 1])
?/OF=C# 5PO_qr= Hx grid on
= 4MTb_ <HoCt8>U Rp=abs(rp).^2;
!{r2`d09n) udq rHR5 Rs=abs(rs).^2;
KR#,6 z^U+ oG Rn=(Rp+Rs)/2;
w[
v {) U}5fjY Tp=1-Rp;
salC4z3
xcr2| Ts=1-Rs;
v1 8 <~ x6`mv8~9Db Tn=(Tp+Ts)/2;
a~7D 4G #+1|O;PB# figure(2)
u:f ]|Q `Y:]& w subplot(1,2,1);
i"}z9Ae~. 04-_ K plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Z?{\34lPj y017
B<Ou legend('R_p','R_s','R_n')
:oZ<[#p"* 3< &:av3 xlabel('\theta_i')
>8_y-74 UON=7}=$ & ylabel('Amplitude')
9>&zOITTaL K khuPBd2 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
OF7hp5 Cq!eAc axis([0 90 0 1])
ZU`9]7"87B S/X U4i:aV grid on
LG&BWs ! TI D gIK subplot(1,2,2);
Ab ,^y RqTO3K f plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
ML_VD*t9 m `-);y legend('T_p','T_s','T_n')
N1i pK9a "@&TC"YG0 xlabel('\theta_i')
ekhv.; N~ *)Qv;'U=rn ylabel('Amplitude')
RRt(%Wm* 23):OB>S` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^xa, r#N:V ZRxOXt&; axis([0 90 0 1])
aj`&ca8 :ZXd% grid on
=IC.FT} S[F06.(1 
.nD#:86M <IZt] P 
!kE5]<H\
