利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
` #Qlr+X wc7gOrPpm 1、光疏射向光密
TVVL1wZ A&%vog]O clear
@~pIyy\_ /wplP+w2 close all
bt3v`q+V aewVq@ngq! n1=1,n2=1.45;
??& Q"6Oe ,B|~V 3)( theta=0:0.1:90;
jd<`W aO : wedfl a=theta*pi/180;
Le#>uWM Bw^*6P^l rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
$X1T!i[.X )}QtK+Rq rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
IC& >PwXb G9'Wo.$ t tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Jq 0aDf
f 13
`Or(>U ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*o<zo
` y;zp*(}f$h figure(1)
zu8 cMxu G'{=. subplot(1,2,1);
;Fw{p{7< VJW%y)_[ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
;a@%FWc _P6e%O8C# legend('r_p','r_s','|r_p|','|r_s|')
t[HfaW1W JC| j*x(k/ xlabel('\theta_i')
$ ?ayE o+{]&V->gN ylabel('Amplitude')
* E $& | Q0Wv8/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ph@hk0dgr/ 9FB k|g"U) axis([0 90 -1 1])
TmI~P+5w Mr/;$O{ grid on
\0gU)tVZ klkshlk d subplot(1,2,2);
|~)!8N.{ AQAZ+g(IK plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
'3B"@^] {O24:'K& legend('t_p','t_s','|t_p|','|t_s|')
`h%(ZG ~ v6uXik xlabel('\theta_i')
.|ZO2MCd ~kHWh8\b: ylabel('Amplitude')
D(bQFRBY6" Ife/:v title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{'O,G$Ldkr Y.>F fL axis([0 90 0 1])
S fl. &A( Cp!bsasj grid on
,3+ #?H 7T[Kjn^{Oj Rp=abs(rp).^2;
&E& _Z6# _]oNbcbt( Rs=abs(rs).^2;
_x +)Tv B!J~ t8 Rn=(Rp+Rs)/2;
X
?
eCK, iX]tL:,~i Tp=1-Rp;
$, I%g< <c; U 0! m Ts=1-Rs;
B EN
U
^T>P Tn=(Tp+Ts)/2;
Hl"^E*9x 86 $88`/2 figure(2)
5t=7- KE$ I!$zO subplot(1,2,1);
zE,1zBS< TzSEQ S{ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
& 9j*Y TUy
25E legend('R_p','R_s','R_n')
W!Qaa(o? pB(|Y]3A xlabel('\theta_i')
J 2H$ALl 8'<RPU} M ylabel('Amplitude')
7)-uYi]
dA y=HM]EH> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]a=n(`l? x->H~ / axis([0 90 0 1])
T<ka4 _&SST)Y| grid on
jNbU{Z%r [ 8 I*lsS subplot(1,2,2);
L9Z\|L5 xI>HY9i ) plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
eeVz Oq( i;l0)q legend('T_p','T_s','T_n')
%(}%#-X OCN:{ xlabel('\theta_i')
+T8h jOkC mb GL)NI ylabel('Amplitude')
r-e-2y7 y d=NafPM title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+~ n"@ / !n^7&Y[N; axis([0 90 0 1])
3UgPVCT uY;7&Lw
y1 grid on
\Dx5= Lh WupONrH1e 
-IE=?23Do? zq Cr'$ 
8+mu'RZ X
