利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
K=::)/{ P ?~e3 &ux 1、光疏射向光密
o{:xp r=( / ~ %KVe clear
Kv&g5&N, } T2xXbU close all
8 IQtz2 (2ot5x}`j n1=1,n2=1.45;
%r\n% $@_ =wWpP-J& theta=0:0.1:90;
:M3oUE{ D/y bFk a=theta*pi/180;
rqPo)AL 9pn>-1NJ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
f|7\DeY9U Uu
s. rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^2D1`,| N c_grPk2O4 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
> et-{(G I=)hWC/ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
(IqZ@->nw B(g _Gm< figure(1)
S{t +>/ B X ms;[ subplot(1,2,1);
Kb#4ILA !LMN[3M_ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
j l.p'$Fbn kl!wVLE legend('r_p','r_s','|r_p|','|r_s|')
{6;9b-a] Ks^6.) xlabel('\theta_i')
tsOrt3 z}QwP~Z ylabel('Amplitude')
b
\KL;H/ GG064zPq7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8
;d$54
b Ix@&$!'k axis([0 90 -1 1])
=uk0@hy9b z<sg0K8z63 grid on
H`bSYjgM! RoA?p;]< subplot(1,2,2);
U @ ?LP (s<Dd2&.H plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
^H{R+} @bM2{Rh: legend('t_p','t_s','|t_p|','|t_s|')
+y7;81ND }VS3L_
;}/ xlabel('\theta_i')
s^0/"j |7 El_wdbbT ylabel('Amplitude')
n~"$^Vr Ee)[\Qjn title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
B[=( #W fH`P[^N axis([0 90 0 1])
Wt)Drv{@ { 0N[&3Ee8 grid on
(Fq5IGs K 8n4oz#z Rp=abs(rp).^2;
T{V/+ RM 3~,d+P Rs=abs(rs).^2;
q"O.Cbk LTNj| u Rn=(Rp+Rs)/2;
XBd>tdEP i|w8.}0 Tp=1-Rp;
xq-17HKs [Vd[- Ts=1-Rs;
NaVQ9ku7VW EXP%Mk/ Tn=(Tp+Ts)/2;
R1nJUOE4w^ /V3=K Y`_J figure(2)
bCref$|
RZM"~ 0 subplot(1,2,1);
@L.82p{ h A +y plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
tz#gClo xe7O/',pa= legend('R_p','R_s','R_n')
]| N3eu q@b| F- xlabel('\theta_i')
+!QJTn"3 u@ jX+\ ylabel('Amplitude')
lb'GXd % Hzrtl et title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3qY K_M^[ 0ye!R
axis([0 90 0 1])
J]m{ b09F da1]mb=4 5 grid on
k >t )g-,2 ? uYu`Ojzr subplot(1,2,2);
SyAvKd`g 763+uFx^ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
e]jH+IR:> [((;+B legend('T_p','T_s','T_n')
8$(Dz]v|[& {] "]uT# xlabel('\theta_i')
; 7N
Z<k D")_;NLE1 ylabel('Amplitude')
.{;Y'Zc14S ^Rx9w!pAN title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
m*$| GW9 ?<G]&EK~~] axis([0 90 0 1])
J8 Yd1.Qj P"<U6zM\sP grid on
5 H#W[^s " `"qP 
L`th7d" p7:{^ 
rP!#RzL
