利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
]Q"T8drL lsax.uG5x 1、光疏射向光密
uU^DYgs <@Q27oEuA clear
HTL6;87w+] &qbEF3p^@ close all
it}h8:^< Wep^He\: n1=1,n2=1.45;
72;'8 f\
P0% theta=0:0.1:90;
=F@
+~)_ :|bL2T@>[ a=theta*pi/180;
uZld9u TnQ>v{Rx rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
H@(O{ 9Yl; 2^B_iyF; rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
X,49(-~\ x'L=p01 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
naR< &@6xu{o ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
jG1(Oe;# 8c) eaDu figure(1)
RsDSsux qxI$F subplot(1,2,1);
"w?0f[" %V=%ARP| plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
X(tx8~z =K:[26 legend('r_p','r_s','|r_p|','|r_s|')
q|ce7HnK Sn[xI9}O xlabel('\theta_i')
DQ8/]Z{H d}O\:\}y ylabel('Amplitude')
b|_e):V| '#c#.O title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Q>$B.z XHA|v^ axis([0 90 -1 1])
qiet<F fd4gB6> grid on
J!om"h L"jA#ULg subplot(1,2,2);
-Mr_Ao`E Mst%]@TG plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
P7w
RX F{ $]LS!@ Rm legend('t_p','t_s','|t_p|','|t_s|')
R([zlw~B5 bkdXBCBx? xlabel('\theta_i')
"" UyfC[ rfonM~3?' ylabel('Amplitude')
)M<+?R$]; \~8W0q.4M title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
W_@ b. 1 u43Mo\"<&% axis([0 90 0 1])
:WRD<D_4 &S(>L[)9 grid on
V%4P.y waMV6w)< Rp=abs(rp).^2;
]?]M5rP _=0Ja
S>M. Rs=abs(rs).^2;
!BVCuuM>w >8/Otg+h Rn=(Rp+Rs)/2;
-G>J /}1|'?P Tp=1-Rp;
-o~zb-E c!HGiqp Ts=1-Rs;
2D%2k Do(G;D`h+_ Tn=(Tp+Ts)/2;
7KiraKb| n#}@|"J figure(2)
9gA@D%0 HKbV@NW subplot(1,2,1);
vk|f"I 6{;6~?U plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
E"*E[> <N<0 ?GQ legend('R_p','R_s','R_n')
n=j)M XcNL\fl1 xlabel('\theta_i')
D 6'd&U{_ :JSxsA6k ylabel('Amplitude')
7[=G;2< ZNH-0mk title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^; /~$ "yH?df24 axis([0 90 0 1])
{3KY:%6qj D?y-Y
grid on
^nZ=B>Yn2 mD>
J,E subplot(1,2,2);
'?$@hqQn ~H#c-B plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
J8\l'}?& U{dK8~ legend('T_p','T_s','T_n')
xppnBnu$7 Up%XBA xlabel('\theta_i')
Z?S?O#FED bCP2_h3* ylabel('Amplitude')
,>#\aO1n o,j_eheAM title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
;m]
n l_vg 6v{&, q axis([0 90 0 1])
hfJ&o7Dt PJ:!O?KVq grid on
b,:^\HKC r1q'+i
Ruq;:5u %_n%-Qn
zdem}kBIe