利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
Ma =6kX] ^A "lkV7 1、光疏射向光密
bqQO E 4; >6l ;/ J clear
3ES[ N.V# KjwY'aYwr: close all
&QOWW} i
B!h Ebz n1=1,n2=1.45;
H (NT| k#E z theta=0:0.1:90;
4$zFR}f $]H= a=theta*pi/180;
`f 6)Q`n <f.>jjwFE rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
2^ ,H_PS Y(
$Ji12 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
|j~EV~A J Y7kb1UG tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
{NY~JFM Rg?{?qK\K ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
OSa}8rlr' .qIy7_^ figure(1)
~C"k$;(n c.8((h/
subplot(1,2,1);
:(l $^
M Y1fy2\<' plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
b$goF
}b'g j
FPU
zB" legend('r_p','r_s','|r_p|','|r_s|')
oGJ*Rn)Z rW FcIh 5 xlabel('\theta_i')
5w,Z 7I8 #6N+5Yx_[ ylabel('Amplitude')
{C/L5cZ]J xMNNXPz( title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
.L^pMU+!^ YXX36 axis([0 90 -1 1])
YA"Ti9-EV >d{dZD} grid on
ws>WA{]gq b.R!2]T]i^ subplot(1,2,2);
g= FDm* wc7mJxJxA plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
QY~<~<d+G ?E+: ]j_ legend('t_p','t_s','|t_p|','|t_s|')
. # 6n Me gE--h xlabel('\theta_i')
WxVn&c\ .:{h{@a ylabel('Amplitude')
j\`EUC t['k%c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ew
%{ i(d EjR_-8@ FK axis([0 90 0 1])
PPoI>J 9;0V
/y grid on
7Wwp )D c= A(o Rp=abs(rp).^2;
.KLm39j( >~I
xyQp Rs=abs(rs).^2;
Pz]bZPHn 3h9Sz8 Rn=(Rp+Rs)/2;
eyeN rk*2o q&XCX$N Tp=1-Rp;
oA4D\rn8" -}{%Q?rYj Ts=1-Rs;
Zm%VG(l qwq/ Xcv Tn=(Tp+Ts)/2;
(Wd_G-da r] A" Og_U figure(2)
lLu ID uY^v"cw/F subplot(1,2,1);
xS6(K #ZG3|#Q=L plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
"?ap gx 6 9=t#5J#O legend('R_p','R_s','R_n')
<^lJr82 ([:]T$0 # xlabel('\theta_i')
qbS'|--wH v5(q) h ylabel('Amplitude')
;i<$7MR.e g%`i=s&N% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ecr88 6 bTZ>@~$ axis([0 90 0 1])
^"3\i A : )^4 ko grid on
LnP3z5d( wgCvD subplot(1,2,2);
e8$l0gzaD TT'Ofvdc plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
N>+ P WE$ Lltc 4Mzw legend('T_p','T_s','T_n')
&^ V~cJ t?gJNOV xlabel('\theta_i')
YiD-F7hf.* _p\6 29` ylabel('Amplitude')
z0#-)AeS Lp_$?MCD. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ls&+XlrX8 'eD J@4Xm axis([0 90 0 1])
UQ/qBbn rkkU"l$v grid on
94\t1fE &~RR&MdZ2 
>D 97c|?c (=1)y'. 
]1m"V;vZ
