利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
RR^I*kRH QT\ S>} 1、光疏射向光密
To DN^qE+ Q)H Vh[4 clear
xdp{y =,[ gwR ^Z{ close all
JT l
37j q^Oq:l$s n1=1,n2=1.45;
Vx-H W;, In=3#u
,M theta=0:0.1:90;
\w O)w@" )N`ia%p_] a=theta*pi/180;
R7Hn8;.. r 'pFHX rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
h Sr#/d w& rb&^ ei9B rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
]#))#-&1 uZ %b6+( tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L=4?vs "uqa~R{ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
B>#zrCD 8uS1HE \% figure(1)
# C4 LDr!d1A subplot(1,2,1);
qs 52)$ g|e^}voRM plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
B,|M
:#c? `>uV legend('r_p','r_s','|r_p|','|r_s|')
n4(w?,w } ;^ O^&< xlabel('\theta_i')
bo\|mvB~ "op1x to ylabel('Amplitude')
FhAuTZk X#1So .}c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
24 1YJ ,t61IU3" axis([0 90 -1 1])
R 5bt~U vy@;zrs grid on
iN+p>3w^l =EsKFt" subplot(1,2,2);
p5c'gziR X&
O
o1y plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Z]uN9c xgsD<3 legend('t_p','t_s','|t_p|','|t_s|')
~}z p}Pt B58H7NH ;G xlabel('\theta_i')
SECL(@0(^ RZm 5[n ylabel('Amplitude')
=@ gH$Q_1 p^ 9QYR title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:]=Y1*L\) ^X"G~#v=q axis([0 90 0 1])
0ydAdgD zu^?9k grid on
=n9adq
/$&~ 0pk Rp=abs(rp).^2;
T* -*U / 4x e:+sA.N Rs=abs(rs).^2;
</:f-J%U/ /=,^fCCN Rn=(Rp+Rs)/2;
9SC#N 5V @ g~kp Tp=1-Rp;
G/2@ Mn- P}DrUND Ts=1-Rs;
w v9s{I{P h7[VXE Tn=(Tp+Ts)/2;
1K0 9iB 1fViW^l_ figure(2)
JWlH(-U4| >`'#4!}G5j subplot(1,2,1);
iDp]l u pb_mW;JVu plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
~k|~Q\ tvf"w`H legend('R_p','R_s','R_n')
`: i|y 3vQ?vS |2 xlabel('\theta_i')
Crm](Z? 57v[b- SK ylabel('Amplitude')
p
Ic ;9 :NA cad title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
xf % ,UQ ,Xo9gn axis([0 90 0 1])
_Dwn@{[(8 Z9~~vf# grid on
*=~
9? <*2.B~ subplot(1,2,2);
4-ZiKM T/)$}#w0i plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
] bhzB w+2:eFi=/ legend('T_p','T_s','T_n')
wb"Jj &AUL]:<s xlabel('\theta_i')
84e)huAs F{b ET ylabel('Amplitude')
j
7a;g7. Y\dK- M{$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
F! c%& Z xO"5bj axis([0 90 0 1])
IDdhBdQ 1 p+2*c grid on
czdNqk.kh 8
6?D 
xfos>|0N O0WzDD 
3M+hjc.
