利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
oKyl2jg+, 03N|@Tu 1、光疏射向光密
1. A@5* Q m\*;Fx clear
<&w(%<; qm mQH S close all
oZ]^zzoEcg t =LIkwD n1=1,n2=1.45;
{7d(B1[1 VT ikLuH theta=0:0.1:90;
up _Qv#`Q Z G(. Q:1 a=theta*pi/180;
e^!>W %.7Z #su R[K*S rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
X|&v]mJ Y @(i zC&h rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
\7V[G6'{ r0^ *|+
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
~ eNKu d.e_\]o<@ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
eewlK] %*bGW' Cw figure(1)
F TgqE@ ZuLW%z. subplot(1,2,1);
op8[8 pt% g9~QNA plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
pNFVa<D q"LE6?hs legend('r_p','r_s','|r_p|','|r_s|')
2I[(UMI$7 evSr?ys xlabel('\theta_i')
1S{AGgls5 "J(T?|t ylabel('Amplitude')
O 'rz ]0o_-
NI title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
;$.^ m'n<.1;1{j axis([0 90 -1 1])
yDapl( zp'Vn 7 grid on
MA1y@ 99GK6}~TGm subplot(1,2,2);
@| 2sF #1)#W6 h\ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Hb&C; lk :pJK Z2B, legend('t_p','t_s','|t_p|','|t_s|')
@@W-]SR T`$!/BlZ xlabel('\theta_i')
aN5"[& 0<[g 7BbR ylabel('Amplitude')
FZBdQhYF JZup } {a title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
vqhu%ZyP <Z
j>} axis([0 90 0 1])
u>\u}c {Gvv^.H 7 grid on
i'YM9* yN X ..<U}e Rp=abs(rp).^2;
RW`+F|UbE nP R*mbW Rs=abs(rs).^2;
'H8b+ GR,gCtG+L Rn=(Rp+Rs)/2;
c#zx" ,K 'wA4yJ < Tp=1-Rp;
]KM 3G bLsN?_jy Ts=1-Rs;
fR&x5Ika0 o?|
]ciY Tn=(Tp+Ts)/2;
;q33t% j e_+SBN1`P& figure(2)
jZgCDA8Mr! *@ eZt*_ subplot(1,2,1);
VQ7A"&hh dZ kj|Ua~ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
"Z]z9( B8`!A legend('R_p','R_s','R_n')
)|,Zp`2/ z P8rW5/ xlabel('\theta_i')
LI9
Uc\ #i 5@G* ylabel('Amplitude')
h4 E[\<? % 7A?gY81 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:Aa5,{v _ O`f[9^fN axis([0 90 0 1])
tjJi| BIxjY!!" grid on
MZvxcr{x PUEEfq!% subplot(1,2,2);
CwaW> (`v xM:9XhH1 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
?wVq5^ e J'C% legend('T_p','T_s','T_n')
| Zj=E$ =`g@6S xlabel('\theta_i')
Zvkb= =gA n;~ ylabel('Amplitude')
4W<8 u( ,In}be$: title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ey&gZ$|& `Y,<[ Lnr axis([0 90 0 1])
TG 6E^3a P *6wt+twH grid on
)@_ 5}8 -r%k)4 _ 
mJ '5!G &%;K_asV; 
{-l:F2i
