利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
!j8FI Y'[ :(E@Gf 1、光疏射向光密
6+:iy'- e[{0)y>= clear
h/QXPdV KaLzg5is close all
Y U'k#\gi* F x] WCQo n1=1,n2=1.45;
E'f{i:O "~ rP'me2
B theta=0:0.1:90;
%HhBt5w [CTnXb a=theta*pi/180;
` _6C {<O Ucb F|vkI rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
d UD[e, ? CmWeY$Jb rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
,=uD ^n: ;_XFo&@ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
FBG4pb9=~ :ffY6L+ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
JB\UKZXw '"52uZ{ figure(1)
5f rX H[T?\Lq subplot(1,2,1);
`T1 V_)-#=J plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
)| cc X ufj ,T7g^ legend('r_p','r_s','|r_p|','|r_s|')
&Hrj3E (xycJ` N xlabel('\theta_i')
>:SHV W OU E (I3_ ylabel('Amplitude')
@ 8(q$ />C^WQI^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
pR<`H' C-xr" ]#] axis([0 90 -1 1])
1> ?M>vK Nd4f^Y grid on
t 6rRU~;} " > ypIR< subplot(1,2,2);
_=>He=v/ UhF-K#Z9 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
$B+ 8Of kb%;= t2 legend('t_p','t_s','|t_p|','|t_s|')
a9e>i U ';k5?^T xlabel('\theta_i')
]Ee?6]bN q~Hn -5H4Q ylabel('Amplitude')
%X]jaX 7 =}^9 wP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
TvoyZW\?w 1y &\5kB axis([0 90 0 1])
e)IzQ7Zex umBICC]CU grid on
[&[k^C5 Ep_HcX` Rp=abs(rp).^2;
5&g @3j] D>r&}6< Rs=abs(rs).^2;
/V8 #[9K `cn#B
BV Rn=(Rp+Rs)/2;
R+:yVi[F]U 2>9C-VL2 Tp=1-Rp;
;1W6 G=m T.F!+ Ts=1-Rs;
h8j .( $p8xEcQdU# Tn=(Tp+Ts)/2;
S>+|OCl"; PNhe figure(2)
1nM
#kJ" 68|E9^`l subplot(1,2,1);
lBE= (A`
G't$Qx,IC plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
]! &FKy n8 0?N}
legend('R_p','R_s','R_n')
M^Yh|%M ,,.QfUj/& xlabel('\theta_i')
7"D", 1h m 0C @G5 ylabel('Amplitude')
c ]-<vkpV =?* !"&h title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
N% B>M7-= Wb_J(!da axis([0 90 0 1])
(x|T+c"bAX ?<,l3pwqa grid on
h4}84}5d 0auYG>< = subplot(1,2,2);
1\m[$Gs: }tu C} plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
gg/-k;@ Rf :4w ?# legend('T_p','T_s','T_n')
T<>,lQs( a RC"MdcD:]y xlabel('\theta_i')
V &T~zh1
8W7J3{d ylabel('Amplitude')
Lb-OsKU 5;WH:XM title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
4 "ZP 'I; ^DLfY-F+j axis([0 90 0 1])
}>| s=uGW Y h7t"=o grid on
C): 1?@ H8jpxzXv 
**C R}
yV q
i;1L
Kc 
3}1u\(Mf
