2、光密射向光疏 fRiHs\+
Jps!,Mflc
clear FQ5# v{
\?j(U8mB>
close all e*tOXXY1
R&?p^!`%
n1=1.45,n2=1; HkrNt/]
BsJClKp/
theta=0:0.1:90; j'Fni4;
,-):&V:jF
a=theta*pi/180; aePk^?KbB
=#W6+=YN8
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); &:rf80`z.
NunT1ved
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); sS._N@f
$3So`8Bm[$
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); [8ih-k
>1YJETysO
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); qg_=5s
k]Yd4CC2
figure(1) [|vdr.
QgP
UP[
subplot(1,2,1); 2?&h{PA+
Na4\)({
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) +S$x}b'5q
TV} H
legend('r_p','r_s','|r_p|','|r_s|') r'&VH]m
T!8,R{V]4
xlabel('\theta_i') GE|V^_|i
DUaj]V{_^
ylabel('Amplitude') -0Ps.B
?Pa5skqR
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 2vynz,^ET
)gZ yW
axis([0 90 -1.5 1.5]) uKK+V6}!kj
`On%1%k8
grid on ]`LMyt0
_di[PU=Vh
subplot(1,2,2); \]zHM.E1
$. Ih-
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) V
V<Zl
iNcB6,++
legend('t_p','t_s','|t_p|','|t_s|') %5'6^bT
>mz<=n
xlabel('\theta_i') O9bIo]B
W 5-=,t
ylabel('Amplitude') |Gz(q4
,#nyEE
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) YH@^6Be9
H8X{!/,^
axis([0 90 -0.5 3]) 3:s!0ty"
i6aM}p<
grid on `2G 0B@
0`:0m/fsU
Rp=abs(rp).^2; T)MKhK9\Ab
+Y-Gp4"
Rs=abs(rs).^2; o!:
l4RZ!K*X_"
Rn=(Rp+Rs)/2; O|d"0P
W2'u]1bs
Tp=1-Rp; idEhxvAo
U<K)'l6#2n
Ts=1-Rs; J.$N<.
vkp_v1F%+
Tn=(Tp+Ts)/2; ",Mr+;;:[
;O+=
6>W
figure(2) N:_.z~>%
uWkW T.>$
subplot(1,2,1); 7*.nd
,?S1e#
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) 3VaL%+T$,
z#m ~}
legend('R_p','R_s','R_n') -e(e;e
G0)}?5L1J
xlabel('\theta_i') 3s;^p,9
Y
vhvFBx0
ylabel('Amplitude') DZ"'GQSg
#A\@)wJ
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) g%TOYZr!X
H}?"2jF
axis([0 90 0 1]) .~u[rc|<
DHQS7%)f`
grid on fN&@y$
JVydTvc
subplot(1,2,2); )V d^#p
a`I
\19p]
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) e>0gE`8A
- ({h @
legend('T_p','T_s','T_n') cDS\=Bf
m~04I~8vk
xlabel('\theta_i') xu\s2x$
R"W5R-
ylabel('Amplitude') Q<0X80w>
pmRm&VgE.
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) f YSH]!
F%P"T%|
axis([0 90 0 1]) Uo?4o*}
z^vfha
grid on ox*1F+Xri
w p\-LO~
<p/zm}?')
.d I".L