利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
AA&5wDMV> BEfP#h=hr 1、光疏射向光密
xE0'eC5n^ @xq jAcfg clear
`A\|qH5`W t XbMP close all
7uI~Xo ?N :!cNkJa n1=1,n2=1.45;
!_a @autj xDsB%~ theta=0:0.1:90;
4ayZ.`aK /'g /yBY a=theta*pi/180;
n/zTS3< ML)5nJD rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
1( nK | oiKY2.yW rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Bh%Yu*.f I<&(Dg|XQ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
,cF
$_7M >3,t`Z : ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
6u6, 9VG, 2Nau]y]= figure(1)
"QV?C $Fr>'H+i subplot(1,2,1);
5Mb5t;4b Vs:x3)m5j plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
7)FI_uW 1>"Yw|F-|3 legend('r_p','r_s','|r_p|','|r_s|')
q_G O;-b{ LOPw 0@ xlabel('\theta_i')
YU (|i}b `#Yv(a2TY ylabel('Amplitude')
k2@|fe {~=Z%Cj2 Q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c_R AtM<n 0Xe?{!@a axis([0 90 -1 1])
]?/7iM UbDRE[^P grid on
9Z_!}eY2mc j 56Dt_ subplot(1,2,2);
@qaK5 ymr# OP$<S plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
/z(;1$Ld6{ ^j)0&}fB legend('t_p','t_s','|t_p|','|t_s|')
aEJds}eE6) cCWk^ lF], xlabel('\theta_i')
>oOZ Duj K[G =J ylabel('Amplitude')
U 5f<4I \5ZDP3I title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Scfk] DT [$+N"4 axis([0 90 0 1])
|."thTO Obl,Qa:5 grid on
'_`O&rbT +bC=yR Rp=abs(rp).^2;
_go1gf7 {B d 0 Rs=abs(rs).^2;
PRpW*#"EI m~x O;_m Rn=(Rp+Rs)/2;
]u(EEsG/ y G{; kJ P Tp=1-Rp;
^|+;~3<J )K"7= TvY Ts=1-Rs;
sfuA
{c'v gb,X"ODq Tn=(Tp+Ts)/2;
omEnIfQSO zK@DQ5 figure(2)
m@ 2;9 d0"Xlle ld subplot(1,2,1);
rERHfr`OU QaAWO plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
#383W)n
Bl6> y/ legend('R_p','R_s','R_n')
0**.:K <i nTd[-3o xlabel('\theta_i')
<dA D-2O+ nYF;. k ylabel('Amplitude')
O*%@(w 6 }vGW lNd#g title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
G&?,L:^t fSL' +l3 axis([0 90 0 1])
$^XCI%DH ra[*E4P9L* grid on
}wkZ\q[ LaolA qU subplot(1,2,2);
<Jwx| `r$c 53|<u plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
1P17]j2C Tz"Xm/Gy legend('T_p','T_s','T_n')
/O$7A7Tl ^o@N.+`&< xlabel('\theta_i')
[$bK%W{f |[ lmW% ylabel('Amplitude')
wm<` 0} ztRe\(9bL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
=8!FY"c* ]3/_?n-"` axis([0 90 0 1])
g
AZe&"K SJXA grid on
j$u=7Z&E m+ww 
"dN4EA&QJ lz-t+LD@ST 
-{tB&V~+v
