利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
fA" VLQE L7OFZ|gUz 1、光疏射向光密
:V
ZXI#([ y\@INA^ clear
#2*6e sP H%G|8,4 close all
Dg'BlrwbR 0$:jZ/._ n1=1,n2=1.45;
wNQ hg t~kh?u].j theta=0:0.1:90;
W+`T: Mgh 7 D W_G a=theta*pi/180;
?my2dd,| C|-QU rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`g^b Q x Dt glPo_( rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
[I2vg<my X6G 2$| tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
wHE1Jqpo Qa@]
sWcM ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
R>y/Y<5= QUK v :; figure(1)
<}( 'w / 1 ;"DIsz@d subplot(1,2,1);
gH:+$FA Ux+UcBKm- plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Sgq" 3(+%, {N'<_%cu legend('r_p','r_s','|r_p|','|r_s|')
2f0qfF H C(7,3 xlabel('\theta_i')
O%*:fd,o- JN> h: ylabel('Amplitude')
6S&OE k )JXy > q# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!~ fy".|x 0@/C 5 v axis([0 90 -1 1])
(g3@3.Kk) ,?(U4pzX grid on
g66x;2Q ^k5# {?I subplot(1,2,2);
9^1li2z k{ bT c^ huP plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
>B$J y7U?nP ')+ legend('t_p','t_s','|t_p|','|t_s|')
Kdr} 7#c b A02)?L xlabel('\theta_i')
;&`6b:ug ^7iP!-w/ ylabel('Amplitude')
NEW0dF&) vE C#W43l title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
OCv,EZ p@[n(?duC. axis([0 90 0 1])
Z}yd` 7 k\Z@B!VAq grid on
{Y1&GO; gPcOm
b Rp=abs(rp).^2;
Q3rLCg,; yII+#?D Rs=abs(rs).^2;
D8\9n HUD` Y{g[LG`U Rn=(Rp+Rs)/2;
}S 6h1X NuR3]Ja\0 Tp=1-Rp;
Z=9gok\ EqF>=5* Ts=1-Rs;
K8{ef 9=. 7[-6i9 Tn=(Tp+Ts)/2;
:Ny^-4-N 9\!=i figure(2)
VU;98 VfkQc$/ subplot(1,2,1);
vY);7 C=x70Y / plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
,*Wp$ / 5y _ < legend('R_p','R_s','R_n')
iiv`ji 9v?rNJs xlabel('\theta_i')
[E)&dl_k ? $.x%G+ ylabel('Amplitude')
qflOi8 ]e(\<R 6Gf title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
iM M s3 N14Q4v-*x axis([0 90 0 1])
GC7W7B Xa_:B\ic grid on
?G 'sb}. mNKcaM?h subplot(1,2,2);
+zZ]Txb( S~U5xM^s plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
O:Wd
,3_ 2Ws'3Jz legend('T_p','T_s','T_n')
rm4t lw _@(E]E xlabel('\theta_i')
iz3H oj \d-9Ndp
nf ylabel('Amplitude')
J~)JsAXAI =Y*zF>#lP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
TecWv@. ce7 $#
# f axis([0 90 0 1])
>OKc\m2%Q 4@=[r Zb9 grid on
T' a& ?;GXFKy 
8=u88?Bh {v* X}`.h 
"qoJIwl#q
