利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
T/3LJGnY )XDBK* ! 1、光疏射向光密
e P @#I^_ jw:z2:0~ clear
`Eh>E, 4u:SE close all
:hZM $4 r'OqG^6JFN n1=1,n2=1.45;
mH$tG
$ CT[9=wV)m% theta=0:0.1:90;
F0<)8{s ?XHQdN3e a=theta*pi/180;
V'&;r'#O Ab@ G^SLX rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
4/b.;$ \_`qon$9 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
61S;M8tNv e 'K~WNT tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
5skN'*oG /TpTR-\I0 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
<gF]9 %2E A9.TRKb=8 figure(1)
1p}H,\o SboHo({5VA subplot(1,2,1);
qdix@ @ ;bE/(nz M plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
q 6> eb . $&^yp legend('r_p','r_s','|r_p|','|r_s|')
:0r ,.) RsZj xlabel('\theta_i')
FH=2, "A K W .4 9 ylabel('Amplitude')
Oh! {E5!) J~2SGXH)^? title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
5%I3eL %s =Wcv b?;* axis([0 90 -1 1])
7 9t E w5%Yi { grid on
D~C'1C&W ab6I*DbF subplot(1,2,2);
[J,.?'V zS%XmS\ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
OD8
fn [ ~9UsHfH legend('t_p','t_s','|t_p|','|t_s|')
h!w::c V UMGiJO\yH xlabel('\theta_i')
s)E \ <w9~T TS ylabel('Amplitude')
d VbFMQ&
^ }7O|Y7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(uC8M,I\ !eF(WbU0 axis([0 90 0 1])
@"7S$@cO p [Po*c.b grid on
yD3 }USw ~ XOmxz0 Rp=abs(rp).^2;
l bHgxZ !l%: Rs=abs(rs).^2;
8<L{\$3HP| jo e)b Rn=(Rp+Rs)/2;
b > D fmW{c mr| Tp=1-Rp;
Jy(G
A yx]9rD1cz Ts=1-Rs;
YlrN ^rO ZwUBeyxS=c Tn=(Tp+Ts)/2;
jYp!?%! i7#4&r figure(2)
Bv<g Vt L8`v subplot(1,2,1);
0ID9=:J =~;~hZj plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
0/GBs~P ng%[yY legend('R_p','R_s','R_n')
-EiTP:A ]rv\sD`[ xlabel('\theta_i')
e0`z~z]6& cB uuq ylabel('Amplitude')
^-"Iw y b? );
D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\bARp z?a A6]:BuP;c axis([0 90 0 1])
&ksuk9M gzs \C{4D grid on
"6jt$-? 3x;UAi+& subplot(1,2,2);
KfiSQ!{ &>\;4E.O5 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
;\pVc)\4" l 2Sar1~1 legend('T_p','T_s','T_n')
'-v:"%s| (h0@;@@7hW xlabel('\theta_i')
R/~ !km ^2k j O/ ylabel('Amplitude')
~[Tcl =ANr|d title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
z x-[@G <U3X4)r axis([0 90 0 1])
=K kqk
'FX M7D grid on
w
% Hj' V}s/knd 
/WnCAdDgZ Eg)24C R 4 
&Sd5]r@+
