利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
rXVR X#Lh vn ^* 1、光疏射向光密
{#&j W <4%PT2R clear
PGA
`R m>+ e;5 close all
R?}<Cj I G9h B p n1=1,n2=1.45;
*[tLwl. TlJ'pG 4^ theta=0:0.1:90;
8s>OO& [[KIuW~ot a=theta*pi/180;
LR
y&/d J pKCux rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
zJG=9C ? CwsC)]{/o rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
K<f B]44Y iH)-8Q tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
w s7LDY&( jWh}cM= ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
d2 *uY., 0-8'. C1v figure(1)
rG{,8* Q&e*[l2M6 subplot(1,2,1);
nh>lDfJV< ykNPKzW: plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
77;|PKE / ;b ^" b{ legend('r_p','r_s','|r_p|','|r_s|')
@!%HEs!# # idNg&' xlabel('\theta_i')
n hG h5, pt~b=+bBm ylabel('Amplitude')
U Kf0cU K"VRHIhfg title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%Sw hNn `yrJ }f axis([0 90 -1 1])
U 46Z~B iTJE:[W"y grid on
YI),yj 9l}G{u9a subplot(1,2,2);
%Q|Hv jk=E [u7i)fn5? plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
{GS $7n I5w> *F legend('t_p','t_s','|t_p|','|t_s|')
L*1yK* U$j?2|v-x xlabel('\theta_i')
Uw->5 1D)=q^\I ylabel('Amplitude')
nmUMg QP!0I01 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%#@5(_' &#<>fT_ axis([0 90 0 1])
B^"1V{M U
)J/so ) grid on
NSQ#\:3:S [,ZHn$\ Rp=abs(rp).^2;
S"Drg m.
iKEHwm Rs=abs(rs).^2;
u51Lp | gP%8nh'C Rn=(Rp+Rs)/2;
X&cm)o%5Fe l&uBEYx Tp=1-Rp;
2q%vd =T sjwD x0(7= Ts=1-Rs;
GA2k g7 L'<.#(| Tn=(Tp+Ts)/2;
e2L4E8ST< d2O x:| <) figure(2)
b1-' q^M 3V@!}@y,F6 subplot(1,2,1);
i
E)Fo.H ;BYv&(#u1q plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Iz[wrtDI 1 `8G {-_ legend('R_p','R_s','R_n')
3J w}MFFV c_FnJ_+ +f xlabel('\theta_i')
v?`DP *&~wl(+O= ylabel('Amplitude')
' +E\-X `6 ?.ihV title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
jQ 9i<-zc i|QL6e*0 axis([0 90 0 1])
ZMmf!cKY:' ==Bxv:6 grid on
|RZI]H% P ze{5! subplot(1,2,2);
R]"3^k* 's 'H&sa plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
3Tz~DdB n_@cjO legend('T_p','T_s','T_n')
FP9FE `x XcM.<Dn3 xlabel('\theta_i')
|J ^$3RX ^6FU] ylabel('Amplitude')
)A: |8 m "wj-Qgz title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]%Z7wF</ %S]g8O[}nl axis([0 90 0 1])
GKa_6X_ 6'qu[ ~ }Q grid on
2*}qQ0J s[#ww
=T\ 
%mMPALN]{ Kld#C51X f 
nbdjk1E`~
