利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
5lT lZRH1 8qT/1 b 1、光疏射向光密
JH 8^ZP:d' o!Y7y1$ clear
IMj{n.y4 h T<v8 close all
(u Sfr]89' [,VD^\ n1=1,n2=1.45;
d;`JDT @sP?@< C theta=0:0.1:90;
#<{sP 0v* 5( 3tPbm{ a=theta*pi/180;
$(BW |Pc ~MOIrF rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
HM`;%0T0( Kg 56.$ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
HJDM\j*5 ,a}+Jj{ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
8q_nOGd D_g+O"];P ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
~x2azY2DP d;K,2 figure(1)
cY5h6+ _ BV|LRB}G subplot(1,2,1);
GujmBb PA Jt M plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
o<Q~pd#Ip, )Zox ;}WK+ legend('r_p','r_s','|r_p|','|r_s|')
W 5-=,t |Gz(q4 xlabel('\theta_i')
,#nyEE YH@^6Be9 ylabel('Amplitude')
H8X{!/,^ 3:s!0ty" title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
i6a M}p< 04o (0 5K axis([0 90 -1 1])
qeypa ! X-=49) grid on
jx J5F3d V3
2 F subplot(1,2,2);
12.|E d*72 )
}(Po_ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
`ml _vm ~yKId legend('t_p','t_s','|t_p|','|t_s|')
b-,4< H8m ,XU<2jv] xlabel('\theta_i')
R10R,*6> iU+O(vi ylabel('Amplitude')
y2cYRHN[X} Dr)B0] KG title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0l;TZf=H ^? }-x axis([0 90 0 1])
0-6rIdDTM >4
VN1 ^ grid on
p C^d-Ii MR }=tO Rp=abs(rp).^2;
I;FHjnn( n&1q * Rs=abs(rs).^2;
}Y:V&4DW W^k95%zBM Rn=(Rp+Rs)/2;
k..AP<hH BlnR{Y Tp=1-Rp;
id+ ~ V #Pt_<?JtV Ts=1-Rs;
xa8;"Y~"bg ; Nk,bb K Tn=(Tp+Ts)/2;
Q`kV|
pjg $t0o*i{ figure(2)
X lLG/N
DaP,3>M subplot(1,2,1);
!y+uQ_IS@ 52ExRG S plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
F/V -@SF %tB7 &%ut legend('R_p','R_s','R_n')
y*=sboX T[Lz4;TRk5 xlabel('\theta_i')
'h R0JXy F_G .$a Cc ylabel('Amplitude')
(1 (~r"4I Hsn'" title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Wn2'uZ5If bb M^J axis([0 90 0 1])
sKCYGt$ ml@;ngmp. grid on
L I*=T Qo32oT[DM subplot(1,2,2);
#/_{(P >
a;iX.K plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
LAk Bf )pSA|Qt N legend('T_p','T_s','T_n')
bmu] zJ ]"} BqS0 xlabel('\theta_i')
C#-HWoSi Dj>eAO> ylabel('Amplitude')
)xq =V tk]D)+{u&c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
& .0A% Sl/]1[|mb axis([0 90 0 1])
,Qx]_gZ` }`kiULC'= grid on
BmKf%:l} ~m_{&,C A. 
nNq<x^@83 C$q};7b1N 
%N)B8A9kh
