利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
#o y vsS8 6e~+@S 1、光疏射向光密
06Uxd\E~ +;; fw |/ clear
M6DyOe< Y|B /( close all
@uH7GW}$g h)A+5^:^ n1=1,n2=1.45;
Jw
{:1 e,1Jxz4QH theta=0:0.1:90;
> O\-\L Pv0OoN*eJ{ a=theta*pi/180;
]]`+aF0 09x\i/nb rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
#+_Oy Z* VQo7 se1P rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
4r5,kOFWb mX |AptND tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
QAb[M\G _B W$?:)9 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Gl1`Nx0 &+sO"j4<?r figure(1)
;=u HK'{ F ^E(A E subplot(1,2,1);
9"V27"s pl"|NZz
7; plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
5~.\rcr% y?5*K legend('r_p','r_s','|r_p|','|r_s|')
v3]q2*`G# C2Y&q X, xlabel('\theta_i')
=20Q! wcu G_AAE#r` ylabel('Amplitude')
.s2d pr;L~$JW title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
gXH89n ~:+g+Mf~[ axis([0 90 -1 1])
c\;} ov+ ~*2PmD"+: grid on
$6'xRUx X m9PcD hv subplot(1,2,2);
`I'=d 4 z1V 0WDVm plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
/pyKTZ| w=^*)jZ8 legend('t_p','t_s','|t_p|','|t_s|')
UT^t7MY#O I1}{7-_t xlabel('\theta_i')
0 1w=;Q oY0*T9vv+ ylabel('Amplitude')
o[cKh7&+ z%;\q$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c6lEWC: aa".d[*1 axis([0 90 0 1])
z5{I3 Y!1 *#2`b%qh\M grid on
WVo%'DtF` r!x^P=f,MJ Rp=abs(rp).^2;
4&$G;?#W2 a7d782~ Rs=abs(rs).^2;
. upcUS8 {) .=G Rn=(Rp+Rs)/2;
J'7){C"G$ ucA6s:!={ Tp=1-Rp;
e=F'
O]
5 "0/OpT7h7 Ts=1-Rs;
h.7 1O"N uvmNQg
Tn=(Tp+Ts)/2;
m ' aakq @D 60 figure(2)
}e@j(*8 [J*)r8 ys subplot(1,2,1);
H$[--_dI{
-l,ib=ne plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Ot,sMRk' T+~~w'v0 legend('R_p','R_s','R_n')
#2WBYScW0 bMv9f
J xlabel('\theta_i')
mS );bs MWB uMF ylabel('Amplitude')
Q 5jP`<zWU @FbzKHdV/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\Xkx`C kv'n W axis([0 90 0 1])
6dinC <[} V K NC K grid on
.z{7
rH 8 XY4 subplot(1,2,2);
| ]I?^:I d-?~O~qD|! plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
Ne#nSx5, &O%Kj8)
legend('T_p','T_s','T_n')
60p1.;' /a yDyq. -Q xlabel('\theta_i')
Dv hK0L*Qr
: zZtZT! ylabel('Amplitude')
2>H\arEstR 0CDTj,eK title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
zY#U ]Is Sqb#U{E axis([0 90 0 1])
", |wG7N
K C&;'Pw9H grid on
-bj1y2)n @*_K#3 
Et- .[ =W4cWG?+ 
$ &^
,(z9
