利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
u?>8`]r 9A|9:OdG1 1、光疏射向光密
n;:C{5 =+[`9 clear
~at:\h4: 0bSnD|#I close all
v_pFI8Cz) I=
cayR n1=1,n2=1.45;
t8.3 8.7lc2aX theta=0:0.1:90;
~0}d=d5g e)87
&
7 a=theta*pi/180;
gkca{BJ 41Ab, rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<7gv<N6BQf b?,=|H rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
R+=wSG ] `}FZ;q3DP tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
[81q 0@ nV8iYBBym ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
!HP=Rgh I%B\Wy/j^ figure(1)
x`2du/
C Q|U
[|U subplot(1,2,1);
]*Kv[%r07c PR,8c plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
lvLz){ wLvM<p7OX legend('r_p','r_s','|r_p|','|r_s|')
4[Wwm ][YC.J xlabel('\theta_i')
k[@/N+;")` eF' l_* ylabel('Amplitude')
JBLh4c3 +c__U
Qx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
hf7[<I,jov x,fL656t axis([0 90 -1 1])
A8:eA ~)6EH`- grid on
IA`Lp3Z |=V~CQ] subplot(1,2,2);
+/" \.wYv j[dgY1yE: plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
n8`WU3& Ry? f; s legend('t_p','t_s','|t_p|','|t_s|')
J6<O|ng:: &)_
z! xlabel('\theta_i')
#]Jg> . lNf.x#u ylabel('Amplitude')
P'*Fd3B#A= }XqC'z title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
vI,T1%llu Tk`|{Ph0 axis([0 90 0 1])
2(Aw s=556 grid on
W! J@30 #q xo1uV(c Rp=abs(rp).^2;
U*go}dt"5 m1xR uj] Rs=abs(rs).^2;
;kT~&.,y *.Z~f"SZy* Rn=(Rp+Rs)/2;
91nw1c! a>Zp?*9 Tp=1-Rp;
J"TF@7{p J93xxj Ts=1-Rs;
qVjMflVoay o/oLL w Tn=(Tp+Ts)/2;
~qK/w0=j kv;P2:"| figure(2)
[ugr<[6 G^eXJusOv subplot(1,2,1);
(%"9LYv k|U2Mp plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
xV=Tmu6l {-%8RSK=< legend('R_p','R_s','R_n')
v-/vj/4> !(Y,2{ xlabel('\theta_i')
;k,@^f8 v *`M3jb ylabel('Amplitude')
IIAmx[ b ydAiH*> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
E]=>@EX qwO@>wQ}~ axis([0 90 0 1])
g!9|1z u(\O@5a grid on
X;
6=WqJj NB+/S ;` subplot(1,2,2);
3xiDt?&H ZDov2W plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
tBX71d
T d`~~Ww1 legend('T_p','T_s','T_n')
2U(qyC o$rF-? xlabel('\theta_i')
HxcL3Bh$~} &~2IFp ylabel('Amplitude')
PC|ul{[*} (OA4H1DL^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
oYt 34@{? =itQ@``r axis([0 90 0 1])
t[@>u'YKt 5pK
_-:? grid on
n9n)eI)R A7|L|+ ? AWKJ@&pA9m 3dSC`K )ax>*