利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
p%CcD]o Bj\ oo+L/ 1、光疏射向光密
h/#s\>)T ':T6m=yv clear
+*$@ K'VL {`[u XH?3d close all
z%L\EP;o} >2<
Jb!f& n1=1,n2=1.45;
j}$Q`7-wB1 .^uYr^( |[ theta=0:0.1:90;
QdP)-Fx n )>nfnh a=theta*pi/180;
&?y VLft ) Ft>X9$ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
}9(:W </} ^2!l/(? rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Vu1X@@z A(qy>x-BI tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
0D48L5kH#' %%k`+nK~ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
~,P." uVV;"LVK~ figure(1)
rfc;
>t}0o$\?E subplot(1,2,1);
n/ - d56 m=%WA5c? plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
u6u1> C<\O;-nHH legend('r_p','r_s','|r_p|','|r_s|')
% wS5 m#n EU-=\Y xlabel('\theta_i')
qUF}rl S=r *ZA.O ylabel('Amplitude')
rz3!0P!"K :g=z}7!s title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
N6_<[` S, TK;g axis([0 90 -1 1])
R} aHo0r X3;|h93.a grid on
7tr;adjs <xQHb^: subplot(1,2,2);
EBPm7{&0| ~Fd<d[b ? plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Z &R{jQ, 5n'C6q " legend('t_p','t_s','|t_p|','|t_s|')
7'9~Kx&+ <UcbBcW, xlabel('\theta_i')
bC@9
*/i L0O },O ylabel('Amplitude')
Q>}*l|Ci @ \(*pa title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
bl; C =n NbtNu $%t axis([0 90 0 1])
h&}XG\ioNA GY %48}7 grid on
N\OeWjA F ~L. )<{? Rp=abs(rp).^2;
^q 0O x&X A12 #v, Rs=abs(rs).^2;
3G7Qo r!vSYgee Rn=(Rp+Rs)/2;
kK[m=rTx1$ ,U#$Qb 12 Tp=1-Rp;
h)qapC5z, E%vG# Ts=1-Rs;
.?YLD+\A oX9rpTi Tn=(Tp+Ts)/2;
L-lDvc?5c P]4C/UDS-~ figure(2)
{b^JH2,
v<z%\`y subplot(1,2,1);
{-( B x xh(VQdg plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
_f8<t=R *hp3w legend('R_p','R_s','R_n')
N| dwuBW lxpi xlabel('\theta_i')
eVZ/3o 9oL/oL-J/ ylabel('Amplitude')
5};Nv{km^2 Ih]'OaE title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Jm|eZDp 8Ilg[Drj* axis([0 90 0 1])
}-:s9Lt ZCZY gf@ grid on
85<zl|ZD IG1+_-H: subplot(1,2,2);
E5P?(5Nv |7V:~MTkk& plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
$ 4\,a^ _-^Lr
/`G! legend('T_p','T_s','T_n')
TM8WaH |.x |BJ xlabel('\theta_i')
z
(,%<oX ei>8{v&g
ylabel('Amplitude')
xG05OqKpE gu[ 3L title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
& >I4-D[ }kaU0 P axis([0 90 0 1])
\Ze"Hv C< GS._V& grid on
e 'I13)
opK =Z 
y5L%_
{n aO9a G*9T 
t'bzhPQO)f
