利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
uD\?(LM KC/=TSSXd. 1、光疏射向光密
;XF:\<+ {_7 i8c<s= clear
(61EDKNd9 h]Wr [v close all
bSiYHRH.e _7H7
dV n1=1,n2=1.45;
$w)!3c4 `P@T$bC theta=0:0.1:90;
iIMd!Q.)@ 2;z b\d a=theta*pi/180;
"3Uv]F Mm xlp .l rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Gr7=:+0n|P [BWNRC1 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
O[I\A[* Yk)."r& ? tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ZI f KV]X@7`@ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
MLoYnR^ Y'1S`. figure(1)
3N dq> RP~|PtLw_ subplot(1,2,1);
hWM<
0= gnFr}L&j plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
`7 vHt` !ipR$ dM legend('r_p','r_s','|r_p|','|r_s|')
]^Qn 9`1O"R/ xlabel('\theta_i')
,I7E[LU -7*ET3NSI/ ylabel('Amplitude')
F<L
EQ7T
tasIDoo+!J title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
?\vh 9 [Na N>BZ? axis([0 90 -1 1])
P#R R9 >Q SW*"\X; grid on
OHeT,@(mh E8pB;\Z( subplot(1,2,2);
6A>dhU !D@ZYK; plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
S #MZV@nGF xCg52zkH# legend('t_p','t_s','|t_p|','|t_s|')
M@`;J jtSA $* AYcy7 xlabel('\theta_i')
H7d/X 8d O! ylabel('Amplitude')
gLE:g5v6 Jll-`b 1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uW 3 0ep' 7[> 6i axis([0 90 0 1])
t/cY=Wp U&Wt %U{ grid on
AfX}y+Ah =<{np Rp=abs(rp).^2;
?$*SjZt
{zn!vJX Rs=abs(rs).^2;
d{SG
Cr 9d )Qe~ 8u@? Rn=(Rp+Rs)/2;
}q $5ig {U1?Et# Tp=1-Rp;
z. kvX+7' )6q,>whI] Ts=1-Rs;
l8I /0 `_ fQ i4\m Tn=(Tp+Ts)/2;
V!lZ\) ]^lw*724'> figure(2)
)/pU.Z/ OW3sS+y subplot(1,2,1);
KVn []@# #73F}
tZ^ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
|G?htZF *v9 2 legend('R_p','R_s','R_n')
iPd[l {85Z 7J EbH?lEN xlabel('\theta_i')
-=~| ."O |N%fMPKa ylabel('Amplitude')
) L#i%)+ H@q?v+2 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6^eV"&+@ rEB @$C^ axis([0 90 0 1])
Xs~[ & lu<xv grid on
\Ta"} TF8 NYrQ$N" subplot(1,2,2);
IF44F3(V4 l|kSsP:GO plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
+kzo*zW$L ;1%-8f:lW legend('T_p','T_s','T_n')
qU#1i:(F* E`?3PA8 xlabel('\theta_i')
V:VO[e<e mj9 <%P ylabel('Amplitude')
aqj@Cjk4Z L*&p ! title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
czS+<
w )N7Y^CN~ axis([0 90 0 1])
uI1 q>[ ,N !o grid on
mt,OniU= Q ;[M}MFc/` 
umt.Um .m2 a1_GIM0 
;c4 gv,q@
