2、光密射向光疏 , ~O>8VbF
xR/CP.dg
clear -AdDPWn
ERpAV-Zf
close all ,PTM'O@aU#
?<EzILM
n1=1.45,n2=1; w={q@.
g%
3' i6<
theta=0:0.1:90; =9GALoGL
%^IQ<
a=theta*pi/180; EfrQ~`\
Y 3BJ@sqz
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); qk2E>
Q[biy{(b8
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); Jr2yn{s=S
lFnls6dp
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ub-3/T
SIJ7Y{\.
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); [iub}e0
mgAjD.
figure(1) +'<PW+U$
#'_#t/u
subplot(1,2,1); yV(#z2|
}=[p>3Dd
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) s6,~JF^
*#T:
_
legend('r_p','r_s','|r_p|','|r_s|') dLiiJ6pl*
-tj#BEC[H(
xlabel('\theta_i') J\L'HIs
MJXnAIG?2
ylabel('Amplitude') >Ku4Il+36
!kovrvM6F
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) >G6kF!V
wk|+[Rl;L
axis([0 90 -1.5 1.5])
J=`
8
Lv%3 jj
grid on 6e-ME3!<l
P0l
fK}
subplot(1,2,2); NnLK!Q
pU7;!u:c4%
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) 72dRp!JU
4$xVm,n|
legend('t_p','t_s','|t_p|','|t_s|') ,a #>e
Q,3kaR@O
xlabel('\theta_i') Q,KNZxT,q
2{79,Js0
ylabel('Amplitude') yYP_TuNa
s2O()u-
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) A8'RM F1
f24W*#IX
axis([0 90 -0.5 3]) >+M[!;m}
MB)<@.A0
grid on yB UQ!4e
}Va((X w
Rp=abs(rp).^2; rN}pi@
d{S'6*`D
Rs=abs(rs).^2; }~
D
WB"
ZNBowZI
Rn=(Rp+Rs)/2; q} e#L6cM
7{m>W!
Tp=1-Rp; cq
I $9
|+
F ~zIu'
Ts=1-Rs; t2dsYU/
hEZvi
Tn=(Tp+Ts)/2; *tqD:hiF
rCPIz<
figure(2) :h(HKMSk1
<m-(B"FX
subplot(1,2,1); *Jsb~wta
w%])
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) n41#
>Sc yc-n
legend('R_p','R_s','R_n') ;Nn(
~+\=X`y
xlabel('\theta_i') b/S:&%E
s<YN*~
ylabel('Amplitude') |5~Oh`w
eK7A8\;e
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) x=DxD&I!J
=Y2 Rht
axis([0 90 0 1]) _Raf7 W
v2;E W p
grid on YiZk|K_
kY]"3a
subplot(1,2,2); H:DR?'yW
x# 0?$}f<
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) _[6sr7H!
EE"8s7ZF
legend('T_p','T_s','T_n') rj=as>6B
fu!T4{2
xlabel('\theta_i') rUKg<]&@
|+Wn5iT
ylabel('Amplitude') !<];N0nt#
%6Gg&Y$j!
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 2K:A4)jZ
irlFB#..
axis([0 90 0 1]) g
*,O
HT.,BF
grid on :?xH)J,imk
/r7xA}se^
yY8zTWji_
6:8s,a3&[k