利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
{H s""/sb l!e8=QlJ 1、光疏射向光密
`84pql, v<bq1QG clear
*0M#{HQ fVv#| close all
G3&ES3L 9PhdoREb n1=1,n2=1.45;
tVQfR*= T.="a2iS2 theta=0:0.1:90;
c/
%5IhX? ElAJR4'{*i a=theta*pi/180;
6'ye-}vD- ^zkTV_,cRp rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
w~M5)b ep<A d rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
{LTb-CB >|SB]'C| tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
HRi~TZ?\ xG0IA 7 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
{n%-^9b1{& FW&P`Iu figure(1)
9x.vz Pa2HFy2 subplot(1,2,1);
t+Mr1e 3P Twpq1 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
@8C^[fDL ,S;?3? a legend('r_p','r_s','|r_p|','|r_s|')
[YP8z~ ,R0@`t1 p xlabel('\theta_i')
*JGm b_ Sh#d& ylabel('Amplitude')
>JS\H6 n"Ec %n title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ba|x?kz K,tmh1 axis([0 90 -1 1])
%*OKhrM =O1CxsKt6 grid on
&5/`6-K DU$]e1 subplot(1,2,2);
&J^@TgqL^ '[JrP<~^o plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
,0[8/)$M cRz7.9-< legend('t_p','t_s','|t_p|','|t_s|')
{@g3AG% Tz58@VY V xlabel('\theta_i')
#0`2wuo
{ CU6rw+Vax ylabel('Amplitude')
t)g1ICt z_!P0` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(Z.K3 "m})~va axis([0 90 0 1])
TJ7on.; Svn|vH grid on
?yddr`?W 2j-|.l c Rp=abs(rp).^2;
aGNt?)8WPZ h+zJ"\ Rs=abs(rs).^2;
R|{AIa{} `y0ZFh1>X Rn=(Rp+Rs)/2;
/7|u2!#Ui 8gJ"7,}-' Tp=1-Rp;
JO5~Vj_" +Es3iE @
Ts=1-Rs;
NWwfNb> MR%M[SK1 Tn=(Tp+Ts)/2;
`kyr\+hp N4!YaQQ;} figure(2)
C2AP e{O mW subplot(1,2,1);
cg7NtY W5$jIQ}Bw plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
\% &QIe;:k $ePAsJ legend('R_p','R_s','R_n')
Mp?Ev. /-E>5 w U xlabel('\theta_i')
9+@_ZI- {i~qm4+o ylabel('Amplitude')
^w^cYM, k<"ZNQm$. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
x*i5g`jx =Z2U axis([0 90 0 1])
+%eMm.( Cv{rd##Y8 grid on
IyOujdKa 6;dB subplot(1,2,2);
y2\, L "FcA:7 + plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
"pdG%$ S#!PDg legend('T_p','T_s','T_n')
{\`#,[ 9g"a`a?c xlabel('\theta_i')
PQ@(p% PLg`\| ylabel('Amplitude')
Hh$D:ZO $&n!j'C: title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`iv,aQ ' T$GhE axis([0 90 0 1])
Da_g3z 7 -gt V# grid on
n0@e%=H)I nEQw6q~je
L2P~moVIi i4'?/UPc
s=~r. x