利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
>)8 <d3 m #m 6 W7_ 1、光疏射向光密
Nc\jA= ['DYP-1J clear
Ji e=/:& &2//\Qz close all
3dB{DuQ 1pJ?YV n1=1,n2=1.45;
U`kO<ztk ^wW{7Uq> theta=0:0.1:90;
=(P k7{ p.RSH$ ] a=theta*pi/180;
V\P
.uOI -5u. Ix3
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
IiZX IG4H :IRQouTf:, rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
W&p f%? UdK +,k~m/ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
t`G<}t @jT=SFf ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
hO.. j F ,h}HlU figure(1)
J 7]LMw7 3&5AbIZ subplot(1,2,1);
x>[f+Tc { PS|q? plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
I!Uj~jV P.'$L\ legend('r_p','r_s','|r_p|','|r_s|')
@i=_y+|d_ { 9|S,< 9 xlabel('\theta_i')
)O$T; U ^y2}C$1V ylabel('Amplitude')
drd5o Z
dEK bB title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
G^A }T 3 H':0 axis([0 90 -1 1])
#bCzWg z2god 1" grid on
-MuKeCg i VNHt ]Ewj subplot(1,2,2);
`(VVb@:o L] 3gHq plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
]6;oS-4gu? x_OZdI legend('t_p','t_s','|t_p|','|t_s|')
g#r,u5<*? 7-2,|(Xg xlabel('\theta_i')
r)Fd3)e G ; ylabel('Amplitude')
?jH u, C0-,<X title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h45RwQ5Z "= >8UR axis([0 90 0 1])
@h)X3X Jk,}3Cr/ grid on
O^<\] _l &C.m*^`^ Rp=abs(rp).^2;
aT}?-CUxx }`D-]/T8. Rs=abs(rs).^2;
w02t9vz J%lEyU Rn=(Rp+Rs)/2;
u#`'|ko \9 4<y|SI! Tp=1-Rp;
E9#.!re|^ =8 01nZJ Ts=1-Rs;
28 =L9q
^|lG9z%Foy Tn=(Tp+Ts)/2;
apd"p{ c%x. cbu> figure(2)
a 8.Xy])! {14sI*b16 subplot(1,2,1);
f <l .%B &~sk7iGi plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
y];@ M<<?e V(<(k,8=
legend('R_p','R_s','R_n')
@W\ H%VR CQ^ I;[=d xlabel('\theta_i')
~ Pi CA '^~3 8=FA ylabel('Amplitude')
Xr$hQbl5D *D;VZs0O title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
*[0)]|r g].v axis([0 90 0 1])
k=JT % 'QxPQ cU grid on
,7*-%05[\ b@yFqgJ_ subplot(1,2,2);
]_ &pIBp 5 @61=Au plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
IXt cHAgX R4Si{J*O legend('T_p','T_s','T_n')
P<s:dH" 62{[)jt{ xlabel('\theta_i')
MD=VR(P?eq l$eKV(CZ4 ylabel('Amplitude')
31n|ScXv &{(8EvuDd title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
u(P;) E"1 "U%jG`q axis([0 90 0 1])
ybgAy J{J< Cd51. Sk(l grid on
2 Ik@L, ljRR{HOl 
Al Q N6* v!M+ 
Mq2[^l!qu
