利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
#o yvsS8 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));
#+_OyZ* VQo7se1P 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 _BW$?:)9 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Gl1`Nx0 &+sO"j4<?r figure(1)
;=uHK'{ F ^E(AE 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&qX, 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 m9PcDhv subplot(1,2,2);
`I'=d4 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')
01w=;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 @D60 figure(2)
}e@j(*8 [J*)r8ys 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 <[} VK NCK grid on
.z{7
rH 8XY4 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