利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
.Mbz3;i0 !0+JbZ<%r| 1、光疏射向光密
9;{C IMg& CJ}%W# clear
?}oFg#m-<L 23PG q%R close all
dPlV>IM$z @ JMiO^ n1=1,n2=1.45;
.#gzP2 [q tmq OJ theta=0:0.1:90;
HJ"GnZp< C dn J&N{ a=theta*pi/180;
6x`t{g]f, )nkY_' BV rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
x 5Bk/e' Z@HEj_n rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^ 8WRqQdx oJ^P(] dw tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
q9"96({ \@ Wr
4,Y QM ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
l?e. 9o2- 7!1S)dup figure(1)
{: /}NpA$ ?,z}%p subplot(1,2,1);
cuX)8+ Nn6%9PX_) plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
M`_0C38
O - wzz legend('r_p','r_s','|r_p|','|r_s|')
*dQSw)R G c?a +T xlabel('\theta_i')
9up3[ F$ &C}*w2]0S ylabel('Amplitude')
4#D,?eA 7 00(\ZUj title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)0`C@um , 1`z"7\W axis([0 90 -1 1])
Yy8g( bU Rq -ZL{LR7 grid on
VbYdZCC /vt3>d%B; subplot(1,2,2);
z{q`G wW awRX1:T#;O plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Qs!5<)6
W?& %x(6M legend('t_p','t_s','|t_p|','|t_s|')
Eci\ a] ]L jf?tk xlabel('\theta_i')
UKGPtKE< F4QVAOM]U ylabel('Amplitude')
'/p4O2b, Wwo0%< 2y title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
u8^lB7! e/ [E_9 V%^ axis([0 90 0 1])
4+ n\k 42{~Lhxt grid on
btB% [] /r 5eWR1G Rp=abs(rp).^2;
BtZ yn7a 7<4qQ.deE Rs=abs(rs).^2;
Om&Dw |xG8 \V:^h [ad Rn=(Rp+Rs)/2;
H,J8 M{ XppOU Tp=1-Rp;
qs6 aB0ln *p d@.|^)m Ts=1-Rs;
]:;&1h3'7 xw%0>K[ Tn=(Tp+Ts)/2;
kfNWI#'9
2oW"'43X figure(2)
d9ihhqq3} fA-7VdR `R subplot(1,2,1);
reVgqYp{{- *Ex|9FCt$ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
CLSK'+ l Ac6=(B legend('R_p','R_s','R_n')
:Tc^y%b0
:&Nb w xlabel('\theta_i')
^Cmyx3O^
)jj0^f1!j ylabel('Amplitude')
w !CNRtM:~ GILfbN cd title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
4H g 9N} e!`i3KYn" axis([0 90 0 1])
C~[,z.FvO d-m7 } 2c grid on
@x1-!
~z# c,22*.V/ subplot(1,2,2);
E]6
6]+;0_ V :eD]zq5 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
- di o5a YqG7h,F legend('T_p','T_s','T_n')
5x de; d _
e WcI xlabel('\theta_i')
Si,6o!0k ,Q,^3*HX9} ylabel('Amplitude')
BY*Q_ Et >p/`;Kq@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8fb'yjIC 'S~5"6r axis([0 90 0 1])
#g=XUZ/" u>$t' grid on
JRFtsio* =xrv~ 
f)!Z~t & {$r[5%L\H 
] - .aL
