利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
`$ZX]6G }tF/ca:XPQ 1、光疏射向光密
C~ t?< L)a8W
clear
bTHKMaGWC h8Q+fHDYv close all
pzbR.L}'D $HF. 02{| n1=1,n2=1.45;
8;ke,x >N@tInE theta=0:0.1:90;
+(x^5~QX .X\p;~H
5 a=theta*pi/180;
%%`Nq&' n1 =B rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
g34<0%6jd 8*!<,k="9 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
8i!AJF9IQ} :6S!1roi tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
iZ<^p1i Yz=(zj ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
%'a%ynFs 8W#/=Xh? figure(1)
CL.JalR`b &PaqqU. subplot(1,2,1);
ns[v.YDL ']Gqa$(YC plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
XAwo~E bXF>{%(}E legend('r_p','r_s','|r_p|','|r_s|')
/x,gdZPX lT#&\JQ
xlabel('\theta_i')
C93BK)$} {e\Pd!D?| ylabel('Amplitude')
gKeqf-UWKJ 8]skAh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,(dg]7 v".q578
0B axis([0 90 -1 1])
,yus44w[ zI$^yk-vn grid on
i]J*lM7' [Oen{c9A subplot(1,2,2);
)b`Xc+{> 5Ds/^fA plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
.)
uUpY%K^ "X0"=1R~ legend('t_p','t_s','|t_p|','|t_s|')
lffp\v{w 7[h_"@_A7 xlabel('\theta_i')
T6m#sVq L9Sd4L_e ylabel('Amplitude')
[ -"o5!0< d~F`q7F'?] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
tvXoF;Yq nwswy]e8/ axis([0 90 0 1])
TP o%zZo O4RNt,?l grid on
:f%FM&b (XA=d
4 Rp=abs(rp).^2;
yTzP{I e8g"QDc Rs=abs(rs).^2;
uMVM- (g% xFxl9oM." Rn=(Rp+Rs)/2;
CT,caa E}AOtY5a Tp=1-Rp;
=sy>_ @Gp=9\L Ts=1-Rs;
t\~lGG-p b?9c\-} Tn=(Tp+Ts)/2;
Ot}fGiio Atl`J.;G figure(2)
IN=pki|. pm$2*!1F( subplot(1,2,1);
n@n608 ,K9\;{C plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Vij P; ?~~sOf AP legend('R_p','R_s','R_n')
>2h|$6iWP %x@
D i`; xlabel('\theta_i')
NbOeF7cq+ rt">xVl ylabel('Amplitude')
RrrW0<Ed t`NZ_w / title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
_svEPHU f!6oW( r-L axis([0 90 0 1])
+Go(yS 5v"r>q[
X grid on
)006\W|t9 pRI<L' subplot(1,2,2);
}L
Q9db1 RtVy^~=G plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
~3byAL O@JgVdgf legend('T_p','T_s','T_n')
P^m+SAAB $YZsaw xlabel('\theta_i')
@X`~r8& AA][}lU:5 ylabel('Amplitude')
[MSLVTR 9~+A<X]Hd title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
*hVbjI$ mF!/8qk axis([0 90 0 1])
D|-^}I4 d:/8P985 grid on
V]b1cDx{ 5.gM]si
<viIpz2jh% ,w~0U
;>^oe:@