利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
44'=;/ KYnW7|* 1、光疏射向光密
#S@UTJa
3Dx Z#/! clear
n^* >a <I>%m, close all
j5PaS k&o= %T`4!:vy n1=1,n2=1.45;
>W>##vK /d{g lOk theta=0:0.1:90;
T r SN00 Zx}N Fcn a=theta*pi/180;
9=Y,["br$_ (:_%kmu rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
jHs<s` #h B.]qrS| rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
08zi/g2
3 {D;Xa`:O tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
g2+l@$W ."~7 \E> t ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
4y|xUO: T
[ T 6 figure(1)
DCEvr " ( T)b3N| ONB subplot(1,2,1);
"2)+) Db o>xxmyW| plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
mA&RN"+V Y&?|k'7 legend('r_p','r_s','|r_p|','|r_s|')
>/7KL2* T^ /Gj|N* xlabel('\theta_i')
^m6k@VM z5x ,fQw6O ylabel('Amplitude')
`U {o: p<Zs*
@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
- XE79 fQ
91-P)%? axis([0 90 -1 1])
vDu0 /Nj:!!
AN grid on
PphR4 sIM qT4I Y$h subplot(1,2,2);
8gVxiFjo J{nyo1A plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
s=H/b$v , aRJ!AZ legend('t_p','t_s','|t_p|','|t_s|')
l%sp[uqcg ,^MW)Gf< xlabel('\theta_i')
*=-__|t rK(x4]I
l" ylabel('Amplitude')
2+T 8Y,g s,UN'~e1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
&+E'1h10 |QAeQWP+1 axis([0 90 0 1])
4
|:Q1 "!%w9 grid on
3b g4# c W2r 6jm! Rp=abs(rp).^2;
GL;x:2XA ;
McIxvj Rs=abs(rs).^2;
<?q&PCAn^ R,d70w
(_ Rn=(Rp+Rs)/2;
<FUqD0sQ egx J 3. Tp=1-Rp;
$nGbT4sc U:mq7Rd8 Ts=1-Rs;
@\%)'WU 48 c
D3w Tn=(Tp+Ts)/2;
M5cOz|j/*R zCBtD_@ figure(2)
\p>]G[g an$ ]IN subplot(1,2,1);
WTV3p,;6a Vq .!(x plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
*! r\GGb |
Q1ub S legend('R_p','R_s','R_n')
Wvut) T E &iWtwkz xlabel('\theta_i')
(PsS E:r}+ RS ||KA])J ylabel('Amplitude')
uh C= Du O%B title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
A]tf>H#1 ~`G;=ITo axis([0 90 0 1])
YmO" EWb 6yu*a_ grid on
P x P?hk `L"p)5H subplot(1,2,2);
m]-v IUpb ;G4HMtL plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
gq/ePSa AjL?Qh4 legend('T_p','T_s','T_n')
aiR|.opIb (:fE _H2z xlabel('\theta_i')
Y6 ; 0khp 9U }MXY0 ylabel('Amplitude')
Cdt,//xrz T4H oSei title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
VJ6>3 j?f,~Y<k axis([0 90 0 1])
*&h XJJ[+ ^EuyvftZ grid on
/8$1[[[ 3I>S:|=K 
I*vj26qvg XZe p7d} 
w371.84
