利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
.`p _vS9
Qi}LV"&L 1、光疏射向光密
e{O5y8, q8lK6p\:W clear
z~_\onC r]S9 z close all
IY:O? M +OmSR*fA0 n1=1,n2=1.45;
Rb0{t[IU -a[{cu{ theta=0:0.1:90;
O
o:jP6 r *l^'v9
a=theta*pi/180;
*LVM}| f Yp)U'8{h c rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L_r & 'B )-{~7@yqZ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
S#9SAX [ 6. jZy~ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
O+$70 LA+MX 0* ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
l.juys8s QhUr aZ figure(1)
H-W) Tq_?- ~}AP@t* subplot(1,2,1);
o~K 2K5I C] w< &o plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
h~@+M5r , c
v
9
6F legend('r_p','r_s','|r_p|','|r_s|')
OKVYpf *PEuaRDN xlabel('\theta_i')
Q(6(Scp{ w!RH*S ylabel('Amplitude')
\gkajY -? G= cx c_9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Anv8)J!9u bEKLameKv axis([0 90 -1 1])
mmL~`i/ <3 }l8Z grid on
G,-OH-M! }D+8 K subplot(1,2,2);
l,FoK76G Jf$wBPg plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
DcA'{21 g-FZel
legend('t_p','t_s','|t_p|','|t_s|')
>G2-kL_ P{eRDQ= xlabel('\theta_i')
]LTc )[5Zj dO,05?q| ylabel('Amplitude')
a`*WpP \+ !@ { [I:5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3L?a4,Q"k} a7\L-T+ axis([0 90 0 1])
C4tl4 df9 M[u3]dN grid on
Pfx71*u, QJWES%m` Rp=abs(rp).^2;
I;-5]/, ?w/nZQWi Rs=abs(rs).^2;
z|*6f FE 3 ?F@jEQk Rn=(Rp+Rs)/2;
*hs <Ez.cC gc3 U/
jM Tp=1-Rp;
:QV6 z*#zD -/c1qL dQ Ts=1-Rs;
/'6[*]IZP \)ZX4rs{8 Tn=(Tp+Ts)/2;
O^weUpe\ ' u};z:t figure(2)
#.+* G`m YwWTv subplot(1,2,1);
K)UOx#xe1 %@93^q[\2 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
j :Jdwf ?{,)XFck legend('R_p','R_s','R_n')
h_G|.7! KH>sCEt xlabel('\theta_i')
f^sb0nU ' 5 q L ylabel('Amplitude')
P8 ,jA<W %GCd?cFF title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
>ha Ixs`9 1vK(^u[ axis([0 90 0 1])
< VR&= YJ h;UdwmT grid on
x ETVt q ^% $W S, subplot(1,2,2);
2mU-LQ1 WN =tRe3o0( plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
O\F^@;] F6 `[ ` *@O(y legend('T_p','T_s','T_n')
.Xz"NyW I u~aTgHX% xlabel('\theta_i')
%802H%+ zHc 4e
ylabel('Amplitude')
b;`#Sea 4i\aW:_'i title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
MY*>) us\ r]//Q6|S axis([0 90 0 1])
C+ ]q 7U
)qC}( grid on
xKLcd+h CZ !na0 Y 
-kri3?Y, Wl RZ|. 
CE7pg&dJ)i
