利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
.`p_vS9
Qi}LV"&L 1、光疏射向光密
e{O5y8, q8lK6p\:W clear
z~_\onC r]S9z close all
IY:O? M +OmSR*fA0 n1=1,n2=1.45;
Rb0{t[IU -a[{cu{ theta=0:0.1:90;
O
o:jP6r *l^'v9
a=theta*pi/180;
*LVM}| f Yp)U'8{h c rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L_r &'B )-{~7@yqZ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
S#9SAX [ 6.jZy~ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
O+$70 LA+MX0* ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
l.juys8s QhUraZ figure(1)
H-W)Tq_?- ~}AP@t* subplot(1,2,1);
o~K 2K5I C] w< &o plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
h~@+M5r, c
v
9
6F legend('r_p','r_s','|r_p|','|r_s|')
OKVYpf *PEuaRDN xlabel('\theta_i')
Q(6(Scp{ w!RH*S ylabel('Amplitude')
\gkajY-? G= cxc_9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Anv8)J!9u bEKLameKv axis([0 90 -1 1])
mmL~`i/ <3}l8Z grid on
G,-OH-M! }D+8K subplot(1,2,2);
l,FoK76G Jf$wBPg plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
DcA'{21 g-FZel
legend('t_p','t_s','|t_p|','|t_s|')
>G2-kL_ P{eRDQ= xlabel('\theta_i')
]LTc)[5Zj dO,05?q| ylabel('Amplitude')
a`*WpP \+ !@{[I:5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3L?a4,Q"k} a7\L-T+ axis([0 90 0 1])
C4tl4df9 M[u3]dN grid on
Pfx71*u, QJWES%m` Rp=abs(rp).^2;
I;-5]/, ?w/nZQWi Rs=abs(rs).^2;
z|*6fFE 3 ?F@jEQk Rn=(Rp+Rs)/2;
*hs<Ez.cC gc3 U/
jM Tp=1-Rp;
:QV6z*#zD -/c1qLdQ Ts=1-Rs;
/'6[*]IZP \)ZX4rs{8 Tn=(Tp+Ts)/2;
O^weUpe\ ' u};z:t figure(2)
#.+*G`m YwWTv subplot(1,2,1);
K)UOx#xe1 %@93^q[\2 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
j :Jdwf ?{,)XFck legend('R_p','R_s','R_n')
h_G|.7! KH>sCEt xlabel('\theta_i')
f^sb0nU ' 5 qL ylabel('Amplitude')
P8 ,jA<W %GCd?cFF title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
>ha Ixs`9 1vK(^u[ axis([0 90 0 1])
<VR&=YJ h;UdwmT grid on
x ETVtq ^%$W S, subplot(1,2,2);
2mU-LQ1WN =tRe3o0( plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
O\F^@;]F6 `[ ` *@O(y legend('T_p','T_s','T_n')
.Xz"NyW I u~aTgHX% xlabel('\theta_i')
%802H%+ zHc 4e
ylabel('Amplitude')
b;`#Sea 4i\aW:_'i title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
MY*>)us\ r]//Q6|S axis([0 90 0 1])
C+]q 7U
)qC}( grid on
xKLcd+hCZ !na0 Y
-kri3?Y, WlRZ|.
CE7pg&dJ)i