利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
Tt%} 4{"
OT;cfkf7 1、光疏射向光密
K\aAM;)- O k*Z clear
VXwPdMy*L <ZVZ$ZW~D close all
#52NsVaT@ xHe^"LL n1=1,n2=1.45;
KJ dz v!l= GQ[pG{ _+ theta=0:0.1:90;
K#wK1 Sv @701S(0 '7 a=theta*pi/180;
=U
c$D* UGCox-W" rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
8kS~ENe?o {@45?L( ' rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
e C;!YG Z RG&6FRoq tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
{>h97}P }PZ= `w*O ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
'W(xgOP1 !UcOl0"6 figure(1)
4w;~4#ZPp T
.hb#oO subplot(1,2,1);
$kl$D"*0 nP ;;MX:B plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
-X8eabb LipxAE?O legend('r_p','r_s','|r_p|','|r_s|')
k}U
JVH21k )88nMH- xlabel('\theta_i')
Um\0i;7 ~4 ;s }3e#$L ylabel('Amplitude')
$ rB6< i. M2E$b| title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
lQL:3U0DjU R8 jovr axis([0 90 -1 1])
($ S Lb6 dyO E6Ex grid on
GS4_jvD- Pr/]0 <s subplot(1,2,2);
F@EJtwLd5y Uouq>N plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Nxi)Q$ I7zn>^0} legend('t_p','t_s','|t_p|','|t_s|')
s i?HkJv5 Onw24 & xlabel('\theta_i')
]Uh 1l.O [O 1|75 ylabel('Amplitude')
KoF
iQ? W+hV9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
lkwh'@s. V Dpxk$a axis([0 90 0 1])
%E.S[cf%8& U0IE1_R grid on
N,|r1u 9X# ~K` 1 Rp=abs(rp).^2;
O,i rpQ qd8pF!u|# Rs=abs(rs).^2;
no|Gq>Xp Q7(eq0na Rn=(Rp+Rs)/2;
v|&s4x?D ,C:o`fQ\ Tp=1-Rp;
]EN&S Wh O] ZC+]}/ Ts=1-Rs;
0H+c4IW g_M ^E-3 Tn=(Tp+Ts)/2;
s#P:6]Ar -l[jEJS} figure(2)
iaRR5D- M5O'=\+,F subplot(1,2,1);
zo44^=~% 6#|qg*OS plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
rJV?) =Z vi|Zit legend('R_p','R_s','R_n')
Ay%:@j(E j )";:v xlabel('\theta_i')
PpFsp ( )x yoU2AMH2D^ ylabel('Amplitude')
choL %g} ]3+`` vL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!g2 a|g 2GUupnQkD axis([0 90 0 1])
Ux_<d? p j+Zt .KXjT grid on
9wMEvX70 tW(+xu36 subplot(1,2,2);
[#Fg\2bq_y pl \ b- plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
jsdBd2Gdc \ JU{xQMB legend('T_p','T_s','T_n')
>?V<$>12 Z"D W 2k xlabel('\theta_i')
0!_*S ) (3O1?n [n ylabel('Amplitude')
(Y rR8 f3t. T=S title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~S;! T b0YNac.l axis([0 90 0 1])
/RqhykgZ =GTD"*vwr grid on
u-39r^`5 LzE/g)> 
.[,6JU% +V1}@6k
: 
n^Vxi;F
