利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
lof P$ -+[~eqRB 1、光疏射向光密
j937tn!Q q\xsXM clear
YRQ?: a{H m-]F]c=)w< close all
gO# %*
W b8**M'k n1=1,n2=1.45;
pgPm0 +N
{t| Q9& theta=0:0.1:90;
ce:wF#Qs .rQ cg.8/B a=theta*pi/180;
;g LOd5*0 j\iNag( rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
e!vWGnY XZrzG P( rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
w|f@sB>j JA% y{Wb tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
lib}dk C!Jy;Z=+u ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Wr`<bLq1vs +mN8uU~(kx figure(1)
9<.8mW^68 zEh&@{u? subplot(1,2,1);
@(IA:6 GN 5t|$Yt[ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
\+Y5b} -$I$z o legend('r_p','r_s','|r_p|','|r_s|')
! LCy:>i!d @ 6*eS+t\ xlabel('\theta_i')
G 2uM 6 bnN&E?{hF1 ylabel('Amplitude')
BwAmNW&i iD/+#UT Y title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
P!gY&>EU h6?o)Q>N axis([0 90 -1 1])
jOZ>^5} ` 5n^DP*X grid on
t ?05 yq!peFu subplot(1,2,2);
2OjU3z<J VVfTFi< plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
r @PVSH/ _o&NbD H legend('t_p','t_s','|t_p|','|t_s|')
@2` nBtk %vbov}R xlabel('\theta_i')
jI~$iDdOfs .g 94|P ylabel('Amplitude')
goND S5 } >8&fFq title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
eN`G2eE X9?0`6Li axis([0 90 0 1])
C3~O6<,Jh FGeKhA 8jT grid on
{REGoe=W% h-x~:$Z, Rp=abs(rp).^2;
[fY7| )~`UDaj_ Rs=abs(rs).^2;
Yp\Y]pym qRz /$|. Rn=(Rp+Rs)/2;
A\v53AT olKM0K Tp=1-Rp;
w-C%,1F,/ F I~=A/: Ts=1-Rs;
Ucz=\dO1 T7o7t5* Tn=(Tp+Ts)/2;
mQ9shdvt- P 6({wx figure(2)
7 0EH~ >CwI(vXn subplot(1,2,1);
;wTc_i tKds|0,j| plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
[.6bxK i 1Kq (7 legend('R_p','R_s','R_n')
GL,[32~C jv6>7@<G xlabel('\theta_i')
(Ffa{Tt! 8~T=p:z' ylabel('Amplitude')
qn'TIE. "<f?.l\+ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ZE1${QFkG C0 /G1\ axis([0 90 0 1])
J*m7
d 4^ SLp nVD:'1 grid on
s3'kzw X JIqg[Mao subplot(1,2,2);
S R_ -wD {,? Gj@$ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
O}j@+p%M Bwu?DK legend('T_p','T_s','T_n')
J Y,oXA6O 0k4XVd+Nv xlabel('\theta_i')
+}Mm5^6* ly6 dl ylabel('Amplitude')
B\qy:nr j nc1?c1s,f title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<fA}_BH%] E,cQ9} / axis([0 90 0 1])
J25/Iy*byG 8 qZbsZi4 grid on
;jO+<~YP! ]z | 2 
\ Co
Z+ N /<c ;"o 
LabI5+g
