利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
vErbX3RY2 {ILQ
CvP* 1、光疏射向光密
Q'vIeG "o c"jhbH!u4 clear
l% 3 Q=c s` dkEaS close all
$txWVjR?\ *8(t y%5F0 n1=1,n2=1.45;
X]f#w \p_8YC theta=0:0.1:90;
`^@g2c+d
A*?/F:E a=theta*pi/180;
&vGEz*F NC#kI3 { rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
vFkyfX( m78MWz] Yo rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
knj,[7uh S
-mz xj tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
aF_ZV bS KfN`ZZ< ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
R&d _ WB4w s`7
_J9 figure(1)
tue%L]hc -t706(#k subplot(1,2,1);
g#bfY=C A+Pm "| plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
EKI+Dq, ps4W wk( legend('r_p','r_s','|r_p|','|r_s|')
]((
>i%%~ p{pzOMi6 xlabel('\theta_i')
azATKH+j I'wk / ylabel('Amplitude')
4F P~+ e#Zf>hlAz title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,1.([%z+r 3C5D~9v axis([0 90 -1 1])
w1/T>o 2ucsTh@ grid on
(Os
OPTp I@3c QxI subplot(1,2,2);
6d?2{_} , bm]dz;ljh plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
dzC&7
9$ y{}, {~FA" legend('t_p','t_s','|t_p|','|t_s|')
?tM]. \ SwmPP- n xlabel('\theta_i')
X ;[zfEB #p:jKAc3 ylabel('Amplitude')
7ump:| 1*Fvx-U' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8=_| qy}l/ kl<B*:RqH axis([0 90 0 1])
b "3T(#2<* =;?PVAdu%# grid on
5 @r Zm4U =D:R'0YH Rp=abs(rp).^2;
uL{~(?U $ R)d 7b,_Yd Rs=abs(rs).^2;
QcVtv7+*v 7Mb t*[n Rn=(Rp+Rs)/2;
XIW: Nk!S :FgRe,D Tp=1-Rp;
b.u8w2( vrEaNT$J- Ts=1-Rs;
+53zI |I K*i1! "w Tn=(Tp+Ts)/2;
r H_:7#.E 8$xKg3-3M figure(2)
hx;kEJ g"AfI subplot(1,2,1);
>Ti2E+}[M 9^h% }> plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
qc4 "0Ap' hb\Y )HSp/ legend('R_p','R_s','R_n')
X!+ a;wr =id $ xlabel('\theta_i')
'|R@k_nx D{d$L9. ylabel('Amplitude')
~oR&0e t ') cgx9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7CN[Z9Y^} 1X}Tp\e axis([0 90 0 1])
wLp
t2 b8S 9@a;1Wr/f grid on
,^Srd 20 9<r}s subplot(1,2,2);
<R8Z[H:bV PKs%-Uk plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
2cL )sP} A0k>Nb\c3 legend('T_p','T_s','T_n')
ZvKMRW 4gN Rln- xlabel('\theta_i')
~0{Kga )GKgK;=~ ylabel('Amplitude')
y$N o o)Z I*R$*/) title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Q g .: w PGhZ`nl axis([0 90 0 1])
e[dRHl 1@dx (_ grid on
~ J{{n_G{ TRok4uc 
ABDUp: bbkI}d%(Ng 
Ax*~[$$~%
