利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
2lCgUe)N fm ^J- 1、光疏射向光密
9SBTeJ$RZ Q!!u=}GYK clear
n) j0h- = fy\W=c close all
1u8hnG %L|xmx!c n1=1,n2=1.45;
p'0X>>$ !Ln 'Mi_B theta=0:0.1:90;
Ou4 `#7FR (;;.[4,y a=theta*pi/180;
m1e Sn |)7 o-o - '0l rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
5\EnD, y *10qP?0H rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
?iZM.$![ ^?2zoS#iw tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
!rF1Remw N^M6*,F,J ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
#+N_wIP4 }U(bMo@; figure(1)
H#u N&^+H "L;@qCfhO subplot(1,2,1);
BIS ., yEos$/*u-N plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
jz~#K;3=, A i"MJ6) legend('r_p','r_s','|r_p|','|r_s|')
dW2 2v! .^ 2.h xlabel('\theta_i')
}d>Xh8:%) *kpP )\P ylabel('Amplitude')
+_{cq@c F!ph Tu title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"TP~TjXf q (?g+.]Dt, axis([0 90 -1 1])
;B;@MD, B 5
1N /XE k grid on
lU%oU&P/"S +'Y?K]zbt subplot(1,2,2);
P*B @it }]#z0'Aqsu plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Cn{v\Q~.4 jo0XF] legend('t_p','t_s','|t_p|','|t_s|')
6"+9$nFyW ^eyVEN xlabel('\theta_i')
]R>NmjAI h-|IZ}F7 ylabel('Amplitude')
ce#Iu#qT * >"NUHq title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
k6J&4?xZ Q"D5D
rj axis([0 90 0 1])
YK[2KTlo B=;kC#Emtf grid on
O$;#GpR }!{R;,5/n Rp=abs(rp).^2;
c ~~4eia) rM^2yr7H Rs=abs(rs).^2;
O F={k[ !h4A7KBYG Rn=(Rp+Rs)/2;
k \qFWFR #rF`Hk: Tp=1-Rp;
X.S<",a{qz &tvtL Ts=1-Rs;
9r+ 'DX?> >e&
L" Tn=(Tp+Ts)/2;
i NfAn& l*OR{!3H$ figure(2)
cD6 ^7QF #R:&Irh subplot(1,2,1);
.)$ MZyo 99`w'Nlk plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
} eW<P079 54Rp0o tv legend('R_p','R_s','R_n')
q^< H G] wggB^ }~ xlabel('\theta_i')
_tX=xAO9 $[Z~BfSQ ylabel('Amplitude')
MJqWc6{ n qPD(D{,f$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uV_% &P ;5L^)Nyd axis([0 90 0 1])
,b&h Lht KpYezdPF) grid on
~v;+-*t - e}(\ subplot(1,2,2);
$A9Pi" /*z P{UV3ZA% plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
$l"%o9ICG xSd&xwP legend('T_p','T_s','T_n')
z `t~N {pH# zs4Y xlabel('\theta_i')
YPI,u7- fn,hP_ ylabel('Amplitude')
C0/^ 6Lu"o ZSK_Lux> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(m,H 5 QEf@wv;T axis([0 90 0 1])
0&|0l>wy. YLQ0UeDN' grid on
C{pOGc@ cZ?$_;= 
#f) TAA rkR5>S( 2M 
tF} ^
