利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
}MJy
+Z8& ~8'sBT 1、光疏射向光密
C)[,4wt, fa E t6 clear
(yH'{6g\ Q- cFtu-w close all
.?8;q A K0u |U` n1=1,n2=1.45;
*RI]?j%B u
1>2v theta=0:0.1:90;
0M_ DB= LdYB7T,
a=theta*pi/180;
R5N%e%[ +F R0(T rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
&|s0P x>**;#7) rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Y+kfBvxyf N H[k Ni' tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ih58 <Up5 Z+*9#!?J ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
$njUX SQ; u0JB\)(-/h figure(1)
: y'Ah# W>${zVu subplot(1,2,1);
\{Y 7FC~ I,E?h? 6Y plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Kr?TxhUHd !{ y@od@T legend('r_p','r_s','|r_p|','|r_s|')
2Z+Wu3# C'>|J9~Gz xlabel('\theta_i')
;;!yC GA$V0YQX ylabel('Amplitude')
OSRp0G20k\ Y4J 3-wK5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
M,j U}yD3 +Zb;Vn4 axis([0 90 -1 1])
$iN"9N%l 216 RiSr* grid on
)MtF23k)g 8EZ$g<} subplot(1,2,2);
S 8 kCp ; p\ Q5,eg plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
G~F b Q7HRzA^- legend('t_p','t_s','|t_p|','|t_s|')
Uf1!qP/H ? ^P[e1?SZG xlabel('\theta_i')
B al`y aHV;N#Lx3 ylabel('Amplitude')
?k+xSV +cJL7=V& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0CWvYC%e 1jx:; j axis([0 90 0 1])
h\$$JeSV] 4 R#c hQ grid on
A` N, iJ ($YvF4 Rp=abs(rp).^2;
OTvPU kp* Ab~3{Q]# Rs=abs(rs).^2;
_mj,u64 "SNn^p59k Rn=(Rp+Rs)/2;
'50OgF' ;o
6lf_ Tp=1-Rp;
WO7z N+9VYH"* Ts=1-Rs;
hXcyoZ8 L 8c0lx}Nn Tn=(Tp+Ts)/2;
e|g5=2(Pr& 7!e vm;A figure(2)
l^!
?@Kg,z ](Xb _xMf subplot(1,2,1);
j:Xq1f6a eln)BW# plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
w_ akn t T _jvxc'6 legend('R_p','R_s','R_n')
/{EP*,/* MOQ6 : xlabel('\theta_i')
n"h `5p5' ({ +!`}GY ylabel('Amplitude')
9# 23FK EZgq ?l~5O title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
GiJ *Wp -$t{>gO#Y axis([0 90 0 1])
C>]0YO
k2 k)i3
grid on
kq?Ms|h ^dI42 4 subplot(1,2,2);
?3/qz(bM ]N/=Dd+| plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
L HyB3V 3>3 Kwc~E legend('T_p','T_s','T_n')
zU,9T n#cC+>*>+ xlabel('\theta_i')
rEx nxQ<e ' Y?"{HZ ylabel('Amplitude')
*o]L|Vu ;RH;OE,A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
m1j*mtu M5i%jZ k axis([0 90 0 1])
.14~ J6 H(H<z,$}T grid on
k}S :RK 32(^Te]: 
@g|v;B|{ o m !!Sl 3 
{Je[ZQ$
