利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
d&5c_6oW .<HC[ls 1、光疏射向光密
6<<'bi </|)"OD9 clear
})KJ60B M5F(<,n; close all
|7 ]? >- hEWx. n1=1,n2=1.45;
LSSW.Oz2L F. }l(KuJ theta=0:0.1:90;
cxY$LY!zX >dyhox2*" a=theta*pi/180;
6$;L]<$W> uC- A43utv rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
[gpOu TW ob)Q,;8R rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Q H4k!^ ujiZM tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
'/ihL ^^@L 7Pwg+| ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*}>Bkq9h J6eJIKK figure(1)
kkT3 wP O+p]3u subplot(1,2,1);
Mm "Wk ,.[.SU#V plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Cy<T Vk8 TH<fbd legend('r_p','r_s','|r_p|','|r_s|')
`b#/[3 .F4oo = xlabel('\theta_i')
z<n" {% *e%Dg{_ ylabel('Amplitude')
3T" #T&eL nrF%wH/5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'BqrJfv xy$vYD AFw axis([0 90 -1 1])
@55bE\E?@ MJ1W*'9</W grid on
~ l(G6/R |^Y*~d<H subplot(1,2,2);
Hr
/W6 C {-o7w0d_ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
y> @v>S mgo'MW\ legend('t_p','t_s','|t_p|','|t_s|')
)!N2'Ld y= -{Q xlabel('\theta_i')
tceIA8d6
W" W@WG9X0 ylabel('Amplitude')
BHF{-z \H,V 9!B title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
9@( O\ xr 's =Q.s axis([0 90 0 1])
C ]yvK} 2~#ZO?jE6 grid on
nsZDZ/jx ^:qpa 5^" Rp=abs(rp).^2;
T {hyt L v Rs=abs(rs).^2;
J CM)N8~i \#uqD\DE Rn=(Rp+Rs)/2;
-$AjD?; "B3iX@C Tp=1-Rp;
;X+G6F' hCU)W1q# Ts=1-Rs;
gcX5Q^`a= } 9MW! Ss Tn=(Tp+Ts)/2;
{[l'S 7'-)/Pk figure(2)
{fAh@:{@ z2rQ$O -# subplot(1,2,1);
;6DR .2}?> ~Y1"k]J plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
tfi2y]{A wlm3~B\64 legend('R_p','R_s','R_n')
j)6@q@P/ Q.j-C}a xlabel('\theta_i')
M3hy5 j(b sL!;hKK ylabel('Amplitude')
&@mvw=d ^JYF1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
>g5T;NgH9 0-8ELX[# axis([0 90 0 1])
$=\oJ-(!@S @/_XS4 grid on
x(C]O, /.pa
??u subplot(1,2,2);
[,\'V0 T+ &x{+gZ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
6LSPPMM v&t`5-e-A legend('T_p','T_s','T_n')
U92B+up- [i,5>YIk xlabel('\theta_i')
sW>P- D6 C -x ylabel('Amplitude')
9Q
SUCN_ n r ' YWW title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
y>)mSl@1y
y)N.LS axis([0 90 0 1])
b&hF')_UOz Y-}hNZn"{ grid on
TE*> a5C| ]1/W8z% 
9m<X-B&P lB5 [#z 
Ev GU j$
