利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
6SN$El 0|G Qvhz$W[P> 1、光疏射向光密
%] [6TZ} fW8whN clear
XI58Cy*! OIdoe0JR:O close all
8I,/ysT: 6V6,m4e n1=1,n2=1.45;
\R (Yf!> dEG1[QG theta=0:0.1:90;
$qy ST |^$?9Dn9.L a=theta*pi/180;
K1[(% <Gp kCZxv "Ts rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
*-.,QpgTX w>uo-88 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
vK,.P:n !=rJ~s
F/{ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
(=/ }i' Rq RyZ*n ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
e{7 "7wn= Q9NKQu Su figure(1)
Xwt}WSdF`k ZIikDi h1 subplot(1,2,1);
cSWn4-B@l TxXX}6 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
)w'GnUqWz h ; S? legend('r_p','r_s','|r_p|','|r_s|')
a5@XD_ b +W[NgUrGJ xlabel('\theta_i')
]d -U eL{6;.C ylabel('Amplitude')
y?.l9
T@x_}a:g title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
NG?- dkD J!@`tR- axis([0 90 -1 1])
,ou&WI yC "E}38 grid on
(YKkJ zaMKwv}BR subplot(1,2,2);
Yi+$g c},wW@SF2W plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
G+zI h}9
uhO-0H legend('t_p','t_s','|t_p|','|t_s|')
RI#o9d"x} ISALR{Aq xlabel('\theta_i')
S^*ME*DDz [ %:%C]4 ylabel('Amplitude')
DZ5QC aA G*\U'w4w|* title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
fe$O Pl~ gO,2:, axis([0 90 0 1])
8lfKlXR78 Z z@wbhMV grid on
B96"|v$ p{S #>JTr Rp=abs(rp).^2;
P2>Y0"bY .:V4> Rs=abs(rs).^2;
V/W{d[86G 4VrL@c
@ Rn=(Rp+Rs)/2;
3?:?dy(3z E{W(5.kb;i Tp=1-Rp;
+!Lz]@9K 3}25=% ;[ Ts=1-Rs;
>P[BwL] F=l. 2t*9 Tn=(Tp+Ts)/2;
TDR2){I kQQhZ8Ch figure(2)
w6FV SU]sY nMU[S + subplot(1,2,1);
h(MS>= L qdz qq plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
A
^U`c'$ C3GI?| b legend('R_p','R_s','R_n')
l_z@.</8P@ TSHH=`cx xlabel('\theta_i')
Jl|^ (X( c.Jj ylabel('Amplitude')
!:|[ ?M.` (3fU2{sm title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:Yy8Ie# 1H]E:Bq axis([0 90 0 1])
5KvqZ1L Xb MAcgS grid on
2#g 4R d 0CFMy6 subplot(1,2,2);
n , .t~ j3yz"-53e plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
'W]oQLD^R v'X=|$75 legend('T_p','T_s','T_n')
KM jnY2 ;|H (_J=6k xlabel('\theta_i')
^":Dk5gl >
g=u Y{Rf ylabel('Amplitude')
!-Br? 9 &p;2/H title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
bhg
OLh# 9HJA:k*k| axis([0 90 0 1])
[V _?`M sksop4gu5 grid on
_E < B@@j- 
;?'=*+'> q26%Z)'nf 
j~:N8(=
