利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
zH~P-MqC 6#Z] yk+p 1、光疏射向光密
]6M<c[H> + \Vw:~e clear
e$s&B!qJ !!Ww#x~k$[ close all
>s5}pkAv|e w=XIpWl n1=1,n2=1.45;
KNmU2-%l z6U'"T"a theta=0:0.1:90;
~T;:Tg* 8 ?82 p a=theta*pi/180;
[TA.|7& bIGcszWr rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
%j^QK >% cRs.@U\{R\ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
:lXY% [!6P ]AA|BeL?| tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
zd%f5L(' [ifw}( ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
CtMqE+j^ BlpyE[h
T figure(1)
ZY,$oFdsi 9~`#aQG T subplot(1,2,1);
8|@) #: U Y*`R plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
1av#u:jy~> f5F-h0HF`[ legend('r_p','r_s','|r_p|','|r_s|')
JL
G!;sov Tl yyJ{~ xlabel('\theta_i')
8T pYt)]S =)Hu(;Yv ylabel('Amplitude')
>bWpj8Kv ;@S'8 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
|kh{EUE
; }*bp4 <| axis([0 90 -1 1])
%ud-3u52M8 MUbKl X grid on
O\f`+Q`0 |a03S Zx subplot(1,2,2);
B Evt{q4 3U^Vz9LW plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
K 3\a~_0 2'=)ese legend('t_p','t_s','|t_p|','|t_s|')
Vj4 h#NN$ d;>G xlabel('\theta_i')
Jvc<j:{^w b4wT3 ylabel('Amplitude')
/1Gmga5 gpyio1V> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
B`)sc ~u d:iJU Vpr axis([0 90 0 1])
.f0qgmIyL Ws5N|g grid on
MJX4;nbl A-1K TD Rp=abs(rp).^2;
.aO6Y+Y ~x(|'` Rs=abs(rs).^2;
@+t|Aa^g !i8)si_ Rn=(Rp+Rs)/2;
rD21:1s o]tfvGvU* Tp=1-Rp;
syLdm 3d| ##''d||u Ts=1-Rs;
P9m LhKbZ oPp Tn=(Tp+Ts)/2;
nHeJ2 0 j]&Qai~}Y figure(2)
C.ynOo,W rb@[ Edj subplot(1,2,1);
68GH$ji JAKs [@: plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
OcMB)1uh\ | eCVq(R legend('R_p','R_s','R_n')
i 1w ]j zd 2_k 9 xlabel('\theta_i')
qJs_ahy( 9c4p9b! ylabel('Amplitude')
.?CaU +d LUq2 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
M
2|
k. n;.
M5}O axis([0 90 0 1])
cly p0`,7 E>
pr})^w grid on
v+x rn z `D= OEc subplot(1,2,2);
5"40{3 [4C_iaE plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
Hf H+U& &+0 2S n3A legend('T_p','T_s','T_n')
,F->*= 03)irq% l; xlabel('\theta_i')
K M )MUPr j<)$ [v6 ylabel('Amplitude')
#t Uhul/O :RIqA/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
" LDNkw' zqU$V~5;rG axis([0 90 0 1])
|FT.x9e- Zjn1,\(t~u grid on
%4rlB$x +%[,
m& 
Ff0V6j)ji 6(?@B^S>2 
% +
