利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
@{ *z1{ 1"? 3l`i 1、光疏射向光密
Q xZYy}2 ts% XjCN[ clear
9Jd{HI= C"}CD{<H]M close all
1Z%^U ? d/5i4g[q n1=1,n2=1.45;
z +,l"#Vv 12qX[39/ theta=0:0.1:90;
G x /sJ( n6 D9f~8" a=theta*pi/180;
eH,r% r, JAJo^}}{b rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
C^9G \s' 2f>G rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
H8.Aq\2S L A &W@ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
" #,]` ME; S.& =>
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
NwPGH= V s |B figure(1)
gS{hfDpk,h '-C%? *ku subplot(1,2,1);
Hl-!rP.?0 N fe plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
!qGE R. Sl,X*[ HGd legend('r_p','r_s','|r_p|','|r_s|')
N3/G6wn %]NaHf xlabel('\theta_i')
]^aOYtKX (=#[om( A ylabel('Amplitude')
:yg:sU ZqT8G title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'dJ#NT25 7G6XK axis([0 90 -1 1])
}/) vOUcEd Dxp8^VL grid on
+,oEcCi BpL7s
ej 7 subplot(1,2,2);
^n! j" Gg's gn
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
nE"##2X GM]" $ legend('t_p','t_s','|t_p|','|t_s|')
t<8vgdD O5rHN;\_ xlabel('\theta_i')
M$Sq3m`{! ]0[Gc
\h} ylabel('Amplitude')
~!V5Ug_2 +M )ep\j title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"gt-bo., V$ 8go#5 axis([0 90 0 1])
kHWW\?O ^Om0~)"q grid on
TRLz> m Q 8x)&4o@ Rp=abs(rp).^2;
^nO0/nqz] ;)n kY6- Rs=abs(rs).^2;
)Cas0~ RM ,Yp+&&p . Rn=(Rp+Rs)/2;
#'5C*RO iX&eQ{LB Tp=1-Rp;
n'?AZ4&z k<&zVV ' Ts=1-Rs;
Q>ZxJ!B<k w"PnN Tn=(Tp+Ts)/2;
g3 6:OK" kT@m*Etr { figure(2)
L'Yg$9 Vz )7aUDsu>4 subplot(1,2,1);
mRurGaR 9mmkFaBQ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
~M c'~:{O e7_.Xr~[ legend('R_p','R_s','R_n')
yjF ;%A/0 }O + a xlabel('\theta_i')
h&| S* 7Mj:bm &9 ylabel('Amplitude')
sQ}E4Iq1#S 2FV O@D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3>'TYXs- GXnrVI axis([0 90 0 1])
Ueb&<tS C>:,\=y% grid on
u ]z87#4 E]P7u" 1 subplot(1,2,2);
' 8"nXuL- &} rmDx plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
KSY HG Wt M1nnJp legend('T_p','T_s','T_n')
Be~ '@ ] :SbvsPm xlabel('\theta_i')
W9G1wU g7;OZ#\ ylabel('Amplitude')
4M )oA|1w _u9bZ' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
W,5_i7vr 4J{W8jX axis([0 90 0 1])
Ul Mc8 z iOB*K)U1 grid on
Y_FQB K U XE$;Z'Qhjm 
:I2spBx yTR5*{?j 
;3\'}2 ^|l
