利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
1RCXc> }/ RP,A!pa@ 1、光疏射向光密
@>M8 Pe UCe ,2v% clear
s.VUd R" D LNa 6 close all
zt-'SY +fC#2%VnU n1=1,n2=1.45;
IRlN++I! 1P(%9 theta=0:0.1:90;
wCV>F- UedvA9$&; a=theta*pi/180;
I/^q+l.=`{ , D exJ1 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
A
Z4|&i T F9Ifw><XM rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Oyi;bb<# Sg/: n,68 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
}l,T~Pjb <P+G7!KZ& ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Z[ a O_6L ;[;)P tFz\ figure(1)
,Adus M di8W2cwz subplot(1,2,1);
@PT`CK} V<7R_}^_7 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
o} WB(WsG q,<n ,0)K legend('r_p','r_s','|r_p|','|r_s|')
IW5*9)N ? 66I |0_ xlabel('\theta_i')
i0,%}{` Rf)'HT ylabel('Amplitude')
*Gg1h@& KU1+<OCh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
zkjPLeX eNI kiJ$uS axis([0 90 -1 1])
GCcwEl!K^ ?R|fS*e2EB grid on
X)`(nj |HaU3E*R subplot(1,2,2);
4 ?9soc *4(/t$)pEl plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
^/_\etV r!{w93rPX legend('t_p','t_s','|t_p|','|t_s|')
9F 2w.(m PWf{aHsr xlabel('\theta_i')
:N^@a- hKk\Y{wv' ylabel('Amplitude')
L
LYHr .yctE:n title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}4b B7,j t 3(%UB axis([0 90 0 1])
N oX_? @D. R0uM grid on
v YRt2({}Z Z]mM Rp=abs(rp).^2;
pRQ fx^ On JVJ1Ay/be Rs=abs(rs).^2;
|@ o]X?^ 6MLN>)t Rn=(Rp+Rs)/2;
>>oASo v$gMLu= Tp=1-Rp;
Bq$e|t)' HI"!n$p Ts=1-Rs;
"TBQNWZ ~pd1 ) Tn=(Tp+Ts)/2;
]wkSAi5z* 9B!im\]O figure(2)
>wg9YZ~8 ^D W# subplot(1,2,1);
<|KK v5[ '(3Nopl plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Q|j@# @O 1 G 1#Bb5q: legend('R_p','R_s','R_n')
.o i}S G Uy$U8b-ov xlabel('\theta_i')
Dyouk+08x , 6EZb[;g^ ylabel('Amplitude')
"Cyo<| VgFF+Eg title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
X/`#5<x `_ J^g&y~ axis([0 90 0 1])
y~]I Vl" Y^c, mK^ grid on
%# Wg^l
' c-s `>m subplot(1,2,2);
Kc JP^ :Fi%Cef| plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
ecY ^C3+S 'K;4102\ legend('T_p','T_s','T_n')
.WL\:{G8; Q
!RVD*( xlabel('\theta_i')
( l3UNP Uf7F8JZmM ylabel('Amplitude')
4,w{rm j e\d5SKY title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Xv A0nEi DSDl[;3O{s axis([0 90 0 1])
UALg!M# fncwe ' ;? grid on
d}wa[WRv
[/+dHW| 
9?g]qy,1) MUCJ/GF* 
7+D'W 7Yx
