利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
G1/ 9E]7Etfw 1、光疏射向光密
S$i3/t }; ]f 3 clear
& BQ%df<y\ f}+8m .g2 close all
|BA<> WE p`3$NCJN n1=1,n2=1.45;
XnV$}T:?X #[NNb?`F theta=0:0.1:90;
hOIk6}r4X G>0 hi1 a=theta*pi/180;
Qb&gKQtt@ 3(>NS ?lX rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
JbEQ35r gq an]b_ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
O_y?5 3X bGa "r tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
#BX ^"J{~ ]LZ , >v ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
"3:TrM$|A
kK$*,]iCp figure(1)
pt-
1>Ui nN!R!tJPa subplot(1,2,1);
j-wz7B Af7&; 8pM plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
'.d]n(/lZd z+%74O"c legend('r_p','r_s','|r_p|','|r_s|')
U Zc%XZ`"V 2q* aq% xlabel('\theta_i')
9U~sRj=D Z|qUVD5Ic ylabel('Amplitude')
txXt<]N 4+15` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
f3HleA&& uQtwh08i axis([0 90 -1 1])
:pPn)j$ iZDZ/hohv grid on
3 eP7vy Z7Xic5PI{4 subplot(1,2,2);
)m[dfeqd + y]TNjLpo$ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
9V( esveq )90K^$93" legend('t_p','t_s','|t_p|','|t_s|')
Ug,23 %t<ba[9F xlabel('\theta_i')
A J"/T+g_ B[nkE+s ylabel('Amplitude')
%H'*7u2 <P4*7:jX title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2rJeON ![a~y`<K, axis([0 90 0 1])
=Frbhh57 o:"^@3 grid on
j: /cJt J{bNx8.& Rp=abs(rp).^2;
d65 t"U |=W=H6h* Rs=abs(rs).^2;
f:|O);n M (' VHL! Rn=(Rp+Rs)/2;
!wZ 9P Jh )K0>R Tp=1-Rp;
+1 j+%&). h1 ( MvEt Ts=1-Rs;
OlF5~VAbfb C ^ hCT Tn=(Tp+Ts)/2;
I~.d/!>Z ~ GW8|tw figure(2)
&\/b(|> du k:: |{F subplot(1,2,1);
uepL"%.@7| @Xb>GPVe#L plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
I_h{n{,sr lSk<euCYs legend('R_p','R_s','R_n')
@*rED6zH 7yK1Q_XY> xlabel('\theta_i')
SJdi*> c@1q8, ylabel('Amplitude')
#RA3 T[A /P3s.-sL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[M{EO) ?dcR!-3 axis([0 90 0 1])
9?_ybO~Oq 8 K/o / grid on
>7^+ag~& _noQk3N subplot(1,2,2);
w>W`8P_b@ ax)j$ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
}_{QsPx9 sU P !'Av legend('T_p','T_s','T_n')
'#b7Z?83C FTg4i\Wp xlabel('\theta_i')
* (4TasQu 6Q6l?!|W4 ylabel('Amplitude')
)5Bkm{v3 1xkU;no title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0y3< Ho,+$ H~P"uYKIZ axis([0 90 0 1])
A' \jaB -Jtx9P grid on
G2,r %|7ta h1 D#, 
M K<VjpP0( dTV:/QM 
g#ZuRL
