利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
}NX9"}/ Zq w xi1 1、光疏射向光密
::bK{yZm {@j 0?s clear
f`Km ctI UykOQ-2-n close all
*a }NRf}W 15PFnk6E| n1=1,n2=1.45;
R7axm<PR= w
D|p' N theta=0:0.1:90;
LObS
7U FCO5SX#-g a=theta*pi/180;
=apcMW(zn bpq2TgFj rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
fa#xEWaFr YbS$D rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
mDJF5 I &n:3n tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L0ZgxG3:g M#v#3:&5 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
&=bWXNU. 8NY $Iw figure(1)
,P G d 9OIX5$,S; subplot(1,2,1);
<PBrW#:' |Ld/{&Q r plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Bvbv ~7g ( 2^-Z17Z} legend('r_p','r_s','|r_p|','|r_s|')
v0pyyUqS 8<cD+Jt j xlabel('\theta_i')
^T}6o Ud t- Rp_2t ylabel('Amplitude')
U}9B
wr^ c ;M7 [y& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Rg/*)SKj 4]\ f} axis([0 90 -1 1])
1UR ;} uFd.2 ,XNP grid on
&r0b~RwUv V w j^h subplot(1,2,2);
TF)OBN~/ CuE>=y- "I plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
P%Fkd3e+ .SdHFWx legend('t_p','t_s','|t_p|','|t_s|')
;d"F'd PA5_ xlabel('\theta_i')
p h[
^ve XwlA W7lU= ylabel('Amplitude')
FB<#N+L \ .r~!d| title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
A aLj.HR &<</[h/B/F axis([0 90 0 1])
?y( D_Nt L zh2gU@" grid on
`N 0Mm7 ZVo%ssVt Rp=abs(rp).^2;
9Vp|a&Ana $Iuf(J-5[ Rs=abs(rs).^2;
1( V>8}zn (<?6X9F:N Rn=(Rp+Rs)/2;
np|3 os oi P8~ Tp=1-Rp;
m&)5QX hX#s3)87 Ts=1-Rs;
G([8Q8B4 + [zp v3Uw Tn=(Tp+Ts)/2;
F,0 @z/8a 9W$d'IA figure(2)
2mS3gk MS{pur D subplot(1,2,1);
!6 _lD 0 h>v;1Q O9D plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
b-e3i;T!}~ rK%<2 i legend('R_p','R_s','R_n')
b#p)bcz!I /6 0 `"xH xlabel('\theta_i')
ne nYP0 Uq /#\7/rL ylabel('Amplitude')
mQt';|X @ nOxCni~ T title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
j&(2ze:=*$ l,Un7]* axis([0 90 0 1])
l,}{Y4\G [Af&K22M(X grid on
ZN!OM)@:! T^h;T{H2 subplot(1,2,2);
LjUBV_J BC V<( @c plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
Q+|8|V} w zAEq)9Y"l' legend('T_p','T_s','T_n')
hNB;29r~ kYBTmz} z xlabel('\theta_i')
4 fBgmL @AUx%:}0Y: ylabel('Amplitude')
Z; shFMu q+XL,E title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
rDkAeX0 8 [i#x|`g axis([0 90 0 1])
U#G[#sd> K :kS A^w8 grid on
9_'xq.uP jI@0jxF 
B]dHMLzl ]xQPSs_ 
quc?]rb
