利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
~mARgv 5K 2K'ZkI 1、光疏射向光密
(/N&_r4 x aY#?QjL clear
1 kKfFpN _1&Ar4: close all
xE
w\'tH 4|E^
#C n1=1,n2=1.45;
a!E22k?((z [)1vKaC theta=0:0.1:90;
Jp_#pV*}: >>,G3/ Zd* a=theta*pi/180;
GaG>0 x 4minzrKM\ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
8ZVQM7O *
l1*zaE rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
(X,i,qK/ j}eb
_K+I tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
E SIP+ *H/3xPh,* ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
twq~.:<o MCmb/.&wu figure(1)
f4b9o[,s2e v={{ $=/t subplot(1,2,1);
1wKXOy=v0 L|T?,^ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
R-S<7Q3E0= p8_^6wf g legend('r_p','r_s','|r_p|','|r_s|')
0H}tb}4 z 6C(?R xlabel('\theta_i')
=+Fb\HvX{ o+A1-&qhN ylabel('Amplitude')
kFWwz^x $TXxhd 6 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
31 a,i2Q4 fw jo? axis([0 90 -1 1])
)*;Tt @'y B/1j4/MS grid on
DFK@/. V M`H#Qo5/ subplot(1,2,2);
SQ_?4 s:: [~ s+,OO9) plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
jU j\<aW lD0-S0i legend('t_p','t_s','|t_p|','|t_s|')
)}`z<)3jP Ig `q[o xlabel('\theta_i')
dZ`Y>wH_ Ak9 {P` ylabel('Amplitude')
N | DI
k Qh1pX}X title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
W^tD6H; *#mmk1` axis([0 90 0 1])
AygdAg'\ 9:USxFM grid on
DQ Q]grU q !z"YpYB Rp=abs(rp).^2;
8(% F{&<; )A7^LLzG Rs=abs(rs).^2;
lPR=C0h}@ 3\XNOJH Rn=(Rp+Rs)/2;
w*VN = L/tpT?$fi Tp=1-Rp;
sm/a L^ 4 !E 5FU *s Ts=1-Rs;
AO R?2u =6[.||9 Tn=(Tp+Ts)/2;
L3, /7 rFg$7 figure(2)
f1 x&Fk T7 ,]^
1 subplot(1,2,1);
(u@:PiU/eP Ek)drt7c y plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
6!m#;8 4 Z6Fu~D2U y legend('R_p','R_s','R_n')
_ nMd \) ~d,M}kK xlabel('\theta_i')
PXMd=,} E )_n?>Ar ylabel('Amplitude')
d?*] /ZiR F/ s
n"2 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
k: P n.< X,k^p[Rcu axis([0 90 0 1])
]/=R ABi oe*1jR_J`[ grid on
8Sr' duY?LJ @g subplot(1,2,2);
4Hj)Av <O( z{rV |vQ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
QoZV 6 X0;u7g2Yz legend('T_p','T_s','T_n')
=NF0 E8O BDI|z/~& xlabel('\theta_i')
/+Z*)q+SbT k H.dtg_ ylabel('Amplitude')
{=Zy;Er Z =+Z96 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
1+FVM\<& 6g V*G axis([0 90 0 1])
Dkz/hg:q PK[mf\G\ grid on
su%(!XJQpg B0@
Tz39= 
U,=f}; pm6#azQ 
@_+aX.,
