利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
r)U9u 0 r'xa' 6& 1、光疏射向光密
szmjp{g0 G=yQYsC$ clear
~)oC+H@{ 4fP>;9[F close all
sI*( MhU Ij#%Qu n1=1,n2=1.45;
V :d/;~ Kq-y1h]7H theta=0:0.1:90;
/ 91H! s Y`lC4 * g a=theta*pi/180;
Hb!Q}V+Kb8 $5il]D` rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
~POe 0!} :,GsbNKW rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
?;^_%XSQ* QD\S E tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
#-e3m/> 9;h 1;9sC| ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
-#)xe W.d ] +Gi~ figure(1)
G
V 0q? J
fsCkS subplot(1,2,1);
#n~/~*:i92 (">! vz plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
li/O&@ g` nwO;> Qr legend('r_p','r_s','|r_p|','|r_s|')
7sHtJr {&K#~[) xlabel('\theta_i')
`33h4G QKCk. 0Xe ylabel('Amplitude')
"b -KVZ
W-Hw%bwN/q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
xwK<f6H! y 0 B3*\ H}5 axis([0 90 -1 1])
I:mJWe oX?2fu- grid on
mE^6Zu "%>/rh2Iq subplot(1,2,2);
$lV0TCgba8 sVh!5fby& plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
1s!hl{n<~ m0 =CD legend('t_p','t_s','|t_p|','|t_s|')
LX*T<|c`' M%9PVePOe xlabel('\theta_i')
!^`ZHJ-3>; of{wZU\J+9 ylabel('Amplitude')
rBgLj,/`U/ Fnll&TF title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
q;B-np?U gDAA>U3|$ axis([0 90 0 1])
Gi;e Drgj~ _Vp9Y:mX2 grid on
tLV9b %i( x#Hq74H, Rp=abs(rp).^2;
T(3"bS., ! daXF&q Rs=abs(rs).^2;
7%)4cHZ^$? 6aMqU?- Rn=(Rp+Rs)/2;
;t* 45 }tj@* n_ Tp=1-Rp;
chfj|Ce]x <R]?8L0{h Ts=1-Rs;
dc%+ f 9>%f9 9n Tn=(Tp+Ts)/2;
Cww$ A %} \>9%=32u. figure(2)
ApS/,cV ^pZ(^ subplot(1,2,1);
> cSc
VN`2bp>5I plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Y. Gr(]tk $&lS7} legend('R_p','R_s','R_n')
atW;S99# v*]|1q%/ xlabel('\theta_i')
:K2N7 ?shA KvktC|~? ylabel('Amplitude')
mC J/gWDY ZJ+q<n_ 4} title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}bix+ /]
]km8M^P axis([0 90 0 1])
o t-!_w< 3go!P]) grid on
zy5@K) "C}nS=]8m subplot(1,2,2);
wf8vKl#Kfw s u![S T( plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
JUwP<C[ hDVD@b legend('T_p','T_s','T_n')
w*Kw#m 'U \?
MuORg xlabel('\theta_i')
{b]WLBy DB%=/ \U ylabel('Amplitude')
OL=IUg" (A R- 8 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ct"h.rD ] Ig5J_Z^]b axis([0 90 0 1])
D~2,0K 1V JE+3 grid on
183'1Z$KA ^B]M- XG 
gKS^-X{x
iM} cd$r{ 
n g<|lsZd
