利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
[%:NR p { 3|W< 1、光疏射向光密
XzsK^E0R K2Zy6lGOZ clear
|{Q,,<C *^ BE1- close all
Jsl,r+'H RR|Eqm 3) n1=1,n2=1.45;
Wd# 6Y}: [XK"$C]jHJ theta=0:0.1:90;
5Tq 3L[T5; kE h# 0 a=theta*pi/180;
_i#Z'4?2E ,zaveQ~l rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L X'US-B.! fCR;Fk2B rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
MiRB*eA #S)*MT4ke tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
DGnswN%n1 oGZ9@Y)(T ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
gu!A:Q HS7!O figure(1)
o hCPNm H Vy^^$ subplot(1,2,1);
I( e>ff cae }dHG2 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
[A47O R [(mq8Nb legend('r_p','r_s','|r_p|','|r_s|')
|%} ?*|- Z[VKB3Pb8 xlabel('\theta_i')
zoU. \]#C #0c`"2t&M ylabel('Amplitude')
JQ qD Ud "O`;zC title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Hw
I s7 ~A)$= " axis([0 90 -1 1])
$O#h4L_ s0u$DM 2 grid on
}`(k X] ][ `}bUf epMJ subplot(1,2,2);
t j0vB]c g7pFO cV plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
F-6*
BUqJ V-31x ) legend('t_p','t_s','|t_p|','|t_s|')
)Jc>l;G(M ? g&]*zc^\ xlabel('\theta_i')
b#M<b.R) h$!qb'| ylabel('Amplitude')
jL# ak V =% p"oj]: title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
5^k#fl2 >Et?7@
axis([0 90 0 1])
{:Q2 Itsy VXa]L4jJ9 grid on
<i5^izg A28w/ =e7 Rp=abs(rp).^2;
gDv]n^ & R8E<;^?j Rs=abs(rs).^2;
v#6.VUAw D7 @10;F}[ Rn=(Rp+Rs)/2;
YANEdH`d cA*%K[9 Tp=1-Rp;
p4[W@JV R8 KL4g-d Ts=1-Rs;
!\m.&lk'^ ru&RL
HFV Tn=(Tp+Ts)/2;
1l i`+~L
F ; ?rW`e2 figure(2)
m H?hzxa+ z R 4hu o subplot(1,2,1);
I4*N o TL "]3`' plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
l*eA
?Qz Z j JD@,j legend('R_p','R_s','R_n')
L=$P %tMfOW xlabel('\theta_i')
IH0Uq_ \8pbPo=x ylabel('Amplitude')
eZv0"FK
X 4eKJ\Q=nX5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
YcI]_[ ldxUq,p axis([0 90 0 1])
Io X(Pa gp}S 1 grid on
CR/LV]G ,kP{3.#Q subplot(1,2,2);
T&c[m!}X|t
ql8:s>1T plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
C{Fo^-3 ~$@I <=L legend('T_p','T_s','T_n')
y`N1I 8_ uDxd xlabel('\theta_i')
`8Om*{xg QhZ%<zN ylabel('Amplitude')
<Z5prunov z/TRqD title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
BP7_o63/G }X$>84s>[P axis([0 90 0 1])
`S VR_ B)(A#&nr b grid on
2@H~ nw 0 s)C.e# xl 
us]ah~U6A q\~7 z1 
?]})Xf.A
