利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
-LEpT$v| VHTr;(]hk 1、光疏射向光密
'A9 U[| is}Y+^j. clear
Y(P <9 m: Q)}z$h55 close all
1IV
R4:a GYN Ly d) n1=1,n2=1.45;
I xE }v%& "3:TrM$|A
theta=0:0.1:90;
S#^-VZ~U4x SDICN0X* a=theta*pi/180;
P};GcV- %%f(R7n rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
b0R{cj=<[ \9s x_T rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Q1?0 ]5 wv_<be[?* tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Shb"Jc_i ,N`D{H"F ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
gPg2Ve0Qy Z;nUS,?om figure(1)
+a ((,wAN2 bdkxCt subplot(1,2,1);
7.tEi}O&_g 2x dN0S plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
#~`]eM5`J X3, + aL` legend('r_p','r_s','|r_p|','|r_s|')
7c.LyvM
]<?7Cp P xlabel('\theta_i')
WHvU|rJ .>{I S 4 ylabel('Amplitude')
F otHITw[ [u}2xsS x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'or8CGr^p j9/Ev]im|F axis([0 90 -1 1])
W05>\Rl q X>\*@ grid on
N(BCe\FV [*AWCV subplot(1,2,2);
g?d*cwtU ;(,1p i7| plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
#Do#e
{=+ * oU-V# legend('t_p','t_s','|t_p|','|t_s|')
p$*;> YKO k=):>} xlabel('\theta_i')
Y;6%pm $ #Bgq]6G2 ylabel('Amplitude')
hpO Uz% hCKx%&[^7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
hXx. ' 5%`[& axis([0 90 0 1])
W:z!fh- cPm- )/E)i grid on
njN]0l{p #-Ad0/ Rp=abs(rp).^2;
v9R"dc]0h DR w;.it2 Rs=abs(rs).^2;
37QX ML jwd {CN% Rn=(Rp+Rs)/2;
xRU ~h Q ~M4@hG! Tp=1-Rp;
bxA1fA; ,Y3wXmG Ts=1-Rs;
ie% _- Jf{
M[ z Tn=(Tp+Ts)/2;
d$(>=gzBQ XTOZ]H*^ figure(2)
oK3aW6 \<R .F
subplot(1,2,1);
3Ta<7tEM f8' $Mn, plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
HAr_z@#E xFY<
ns legend('R_p','R_s','R_n')
WgK |r~ AZTn!hrU xlabel('\theta_i')
:&oUI&(o 'o*:~n ylabel('Amplitude')
{k}EWV MlM2(/ok title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'F [ C 4 }_{QsPx9 axis([0 90 0 1])
sU P !'Av \O7Vo<B&D grid on
r\-25F<e5 * (4TasQu subplot(1,2,2);
k-LT'>CWl +.b@rU6H plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
%V CfcM}5I Dxwv\+7] legend('T_p','T_s','T_n')
qJ|n73yn P6E=*^^m( xlabel('\theta_i')
3oCw(Ff QF fK EMN ylabel('Amplitude')
*I(g~p k<M~co;L title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
;|Z;YK@20 s;0eD5b>x axis([0 90 0 1])
g}-Ch# ~',}]_'oR- grid on
"k${5wk#Fl !j3V'XU#Zn 
5w#*JK ; ^,2
Qs M 
|;-,(509
