利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
u5T \_0 4 S03W
1、光疏射向光密
m,)Re8W- B.{0,b W?
clear
"FG6R' hQHV ]xW close all
T[ky7\ "$
u"Py n1=1,n2=1.45;
ap_(/W c;(}Ih(# theta=0:0.1:90;
8kO|t!?:U $,'r}
% a=theta*pi/180;
c5O1h8 =\oNu&Q^ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
?uh7m 2l0D Ur5X~a\y rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
c;wt9J.f EY> %#0 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
%;n y E/5w
H/ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
(lt/ t #(XP=PUj figure(1)
|urohua t%N #Yh! subplot(1,2,1);
+a;: 7[%& sA}=o.\j: plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
&w LI:x5 k(' 2 K2P legend('r_p','r_s','|r_p|','|r_s|')
C]eb=rw$ X*d,z~k%*d xlabel('\theta_i')
]sbu9O ^"f ydoCoD
w ylabel('Amplitude')
. 5_w^4`b ]Fs Plxk6 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Eepy%-\ 1gCp/m2r7 axis([0 90 -1 1])
gIRFqEz@o h9SS
o0]F grid on
MUVp8! *@ okr'=iDg subplot(1,2,2);
W}JJaZR*X >/evL
/ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
ow] n)Te s )_sLt8? legend('t_p','t_s','|t_p|','|t_s|')
cl1 4FrpYu %o*afd xlabel('\theta_i')
MhJq~G p
]Alu~ Dw ylabel('Amplitude')
0'aZ*ozk X~n Kuo title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#WfJz}P,! `M p]iD { axis([0 90 0 1])
vmW 4a3 $6ITa }o grid on
qdO^)uJJ BKV vu}V(o Rp=abs(rp).^2;
WYI? M ZLo3
0* Rs=abs(rs).^2;
_ mFb+8C .6 Rn=(Rp+Rs)/2;
D~8f6Ko"m N b(se*Y# Tp=1-Rp;
:'.-*Ew M.H4 ud Ts=1-Rs;
ilpg() v)rN] b] Tn=(Tp+Ts)/2;
N!Kd VDdT| pOw4H67 figure(2)
+ ZGOv,l z?FZu,h} subplot(1,2,1);
Awe\KJ^` CbK7="48 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
]Jv Z:'g} @{XN}tWDOp legend('R_p','R_s','R_n')
x"_f$,:! gY;N>Yq,C xlabel('\theta_i')
C, jPr )6) tHhY1[A8m ylabel('Amplitude')
0.&gm@A~c$ )pJ}o&J title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,CwhpW\Y p 4(- axis([0 90 0 1])
x"U/M ?l Gf]oRNP,N grid on
"ve?7&G7U >-+M Wu= subplot(1,2,2);
s7T=/SC54 28Q`O$=v plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
5F&i/8Ib s7df<dBC legend('T_p','T_s','T_n')
;#f_e; ^W #[6]S xlabel('\theta_i')
2ZLK`^S mX,#|qLf ylabel('Amplitude')
D {>, 2hC ^k u~m5v title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
=oiY'}%(i j >0S3P, axis([0 90 0 1])
yf_<o xp> <7{ grid on
-c|O !Lc- BAqu@F\): 
TSE(K t B(b[Dbb 
4+l7v?:Pr
