利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
6EY W:o hha!uD~ ( 1、光疏射向光密
ZBxV& .9/ mzH3Q56 4 clear
$?M$^ - (e K[s!3.u close all
C<"b99\2` # _9Jam%M n1=1,n2=1.45;
8;+dlWp X6 SqOb\(a theta=0:0.1:90;
71wtO f)9{D[InM^ a=theta*pi/180;
rpn&.#KS >Dk Rl rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
i0py5Q b 8~7C4 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
<HMmsw {Ip)%uR tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
34s>hm=0. &O;' ?/4
S ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
$geDB~ 2> _a&Mk figure(1)
M?G4k] FF'Ul 4y subplot(1,2,1);
2 i:tPe& biy[h3b plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
1Uf8ef1, o)SA^5 legend('r_p','r_s','|r_p|','|r_s|')
-5.~POO Ps=<@,dks xlabel('\theta_i')
#1VejeTi y>iot e~ ylabel('Amplitude')
z>9gt l>{+X ) title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
-gR
}^D MB3 0.V/\ axis([0 90 -1 1])
,`. `}' V(6GM+ grid on
)uxXG `,h kO'_g1f<[ subplot(1,2,2);
; +%| !~ <rAWu\d; plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
YdiXj |k+ [&H?--I legend('t_p','t_s','|t_p|','|t_s|')
QoTjKck. \r^*4P,, xlabel('\theta_i')
6S6E
1~ t}A n: ylabel('Amplitude')
yMCd5%=M\ RQt\_x 7P title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
avwhGys# r5(-c]E7 axis([0 90 0 1])
<h<4R Rj uU#7SX(uu grid on
9<Kc 9Z zm`^=cV Rp=abs(rp).^2;
8j%hxA V$ *oP&'$P Rs=abs(rs).^2;
Tt*n.HA /m+q!yi & Rn=(Rp+Rs)/2;
o])2_e5 &] e uL:C Tp=1-Rp;
tW7 *(D F G 5e{ Ts=1-Rs;
B}A7Usm L?~-<k Tn=(Tp+Ts)/2;
SYZS@o c N^,-~U figure(2)
hp6%zUR .f~x*@ subplot(1,2,1);
;.wWw" ) ee5QZ, plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
~.'NG?
%7P i][af legend('R_p','R_s','R_n')
S/l?wwD tE>hj:p xlabel('\theta_i')
5Zdx n> Tl ?]K ylabel('Amplitude')
SpG^kI # ;q6FdS title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
d,0Yi
u.p rBUdHd9 axis([0 90 0 1])
j )[
w X =5yI >A0 grid on
USrg,A bWAVBF subplot(1,2,2);
{V)Z!D 5lMm8<v plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
I!;&#LT+b L>qLl_. legend('T_p','T_s','T_n')
xo46L\ 5$ra4+k0 xlabel('\theta_i')
(8*& 42W @'yD(ZMAz ylabel('Amplitude')
m+J3t @$ k 8Swra?j title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^KsiTVY Jc:gNQCsP axis([0 90 0 1])
= os%22* qq+fUfB2: grid on
5$|wW}SA Yy"05V. 
TCd1JF0 ?r,lgaw 
Z ?{;|Z5
