利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
81n%2G _<~05Eh 1、光疏射向光密
<AIsNqr Mfz5:' clear
"s*{0'jo q{@Wn]!k close all
vbFAS:Y:+ 8t^"1ND n1=1,n2=1.45;
f>'7~69 "2h#inS theta=0:0.1:90;
2KG j !w ZD<,h`
lZ a=theta*pi/180;
K4rr.f6 9w(j2i
q rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
a=J?[qrx 2Fh_ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
QZfPd\Q5 :vr,@1c tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
2ReulL8j kj8zWG4KH ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
\uYUX~}i" /}VQzF figure(1)
i" )_M|
J8%|Gd0#4 subplot(1,2,1);
nqH^%/7)A@ EW!$D plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Brl6r8LGi /X:lt^?%I legend('r_p','r_s','|r_p|','|r_s|')
zPmVECS Zu$f[U)X xlabel('\theta_i')
Dux`BKl ]`NbNr]K ylabel('Amplitude')
f4_\F/ J!H5{7.efN title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3UaP7p+d Ao\Vh\rQkq axis([0 90 -1 1])
bXW)n<y q!,do2T grid on
~;+vF-]R Y1~SGg7(@ subplot(1,2,2);
7mjj% d#*5U9\z plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
zm:=d>D.. e!8_3BE legend('t_p','t_s','|t_p|','|t_s|')
(HF,p,h_ h<M7[p= xlabel('\theta_i')
Lr^xp,_ n t2:c@) ylabel('Amplitude')
W
wE)XE A)#sh)
}Q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
w(U/(C7R ,dp?'_q{ axis([0 90 0 1])
K8Y/XEK 2%8Y-o? grid on
IX(yajc[~M I5A^/=bf& Rp=abs(rp).^2;
]~zJ7I pd1m/: Rs=abs(rs).^2;
)eEvyU
{|gJC>f@ Rn=(Rp+Rs)/2;
U{_s1 >2:S v1T Tp=1-Rp;
^zJ.W `|w#K28t" Ts=1-Rs;
"{k
)nr+7U 8_m9CQ6 i Tn=(Tp+Ts)/2;
t/ 1NTa PPPwDsJ figure(2)
wUeOD.;#F 9/M!S[N9 subplot(1,2,1);
>>cd3)b DRoxw24 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
<`3(i\-X qlJOb}$ I legend('R_p','R_s','R_n')
[J`G`s! Zsogx}i- xlabel('\theta_i')
B|=maz:_ 5r<(Z0 ylabel('Amplitude')
eW)I}z+{ S7/v,E title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
UK>=y_FYO P`
F'Nf2U axis([0 90 0 1])
)T5h\ZO`; @>IjfrjV grid on
"T
u[n\8 ->&VbR) subplot(1,2,2);
1$VI\} /A.i5=k plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
`_A?a_[* m tQ{6u
legend('T_p','T_s','T_n')
ls\E%d t)Q@sKT6 xlabel('\theta_i')
!#I/be] U_;J.{n ylabel('Amplitude')
]~Vu-@
/} 'F?Znd2L title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'?gIcWM )x x/di axis([0 90 0 1])
[0ffOTy TDE1z>h+" grid on
>Mz|e(6 |K;Txe_
(OcNC/9 RhNaYO
B;]5,`#!