利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
wf\"&xwh? Xxg |01 1、光疏射向光密
sm/a L^ 4 !E 5FU *s clear
:W *yfhLt u /F!8# close all
F?Lt-a + avR tYL n1=1,n2=1.45;
f1 x&Fk .`OU\LA theta=0:0.1:90;
(u@:PiU/eP @Z~YFnEJi a=theta*pi/180;
>&P M' k 2LtDS?)@ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
c4tw )O-X 1!v{#w{u7 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
ka9@7IFM R5uG.Oj-2 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
6nW)2LV /4an@5.\C ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
%GhI0F # ~XTC:6ts figure(1)
$gUlM+sK S0^ a)#D & subplot(1,2,1);
t eY@) F ,UY1.tR( plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
i/9iM\2 c;VqEpsbl legend('r_p','r_s','|r_p|','|r_s|')
mJUM#ry )zr*Ecz xlabel('\theta_i')
{rKC4 : hC\
l
\y ylabel('Amplitude')
+8Lbz^# Q HR,p/p title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
B)JMughq_ Ek{Q NlQ]4 axis([0 90 -1 1])
$tmdE )"& Ln h =y2 grid on
<YaT r9%w 9J3fiA_ subplot(1,2,2);
> yC=@Uq+ Bh3F4k2bg7 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
pm6#azQ Tv=lr6t8 legend('t_p','t_s','|t_p|','|t_s|')
\Bo%2O%4 8o~
NJ 6 xlabel('\theta_i')
[YOH'i&X O4R\] B#Xu ylabel('Amplitude')
lfgJQzi
G FzInIif title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
r<vMp'u =1IK"BA2? axis([0 90 0 1])
_SBbd9 W@d&X+7e grid on
2aYBcPFQh# g:Hj1!' Rp=abs(rp).^2;
^&>B,;Wu ? o "
Vkc: Rs=abs(rs).^2;
h&0zR# t [^qT?se{ Rn=(Rp+Rs)/2;
P 45Irir | :eTo<
Tp=1-Rp;
e5' I W__ 42dv3bE " Ts=1-Rs;
. &^p@A~ }+@9[Q
L Tn=(Tp+Ts)/2;
[HILK `@@ 0kE[=#'.' figure(2)
j?K$w` J2z /XHS subplot(1,2,1);
<*( Z}p i~';1
.g plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
| v!N1+v0 oZBD.s legend('R_p','R_s','R_n')
!xsfhLZK Pq /5Dy xlabel('\theta_i')
hYJzF.DW<$ -#Z df | ylabel('Amplitude')
)3h\QE!z (D>_O$o title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2k .VTGak @ Xo*TJB axis([0 90 0 1])
@](vFb JEto_&8, C grid on
}A 3/( Q#eMwM# ~ subplot(1,2,2);
@c|=onx5 ~o'1PAW7 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
<%)vl P#@ H*W>v[> legend('T_p','T_s','T_n')
dNe!X0[ ~c)&9' xlabel('\theta_i')
TQ"XjbhU;X Q "&Mr + ylabel('Amplitude')
3TwjC:Yhv2 5Tl5T& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
KF rsXf "^]gI Qc axis([0 90 0 1])
[q9B" @X Hx.|5n,5 grid on
!l [;,l 2 `=jKt 
Bu:%trlgV si0}b~t 
w#.3na
