利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
l8n[8AT1 S r#fyr 1、光疏射向光密
^Ih dq89 t B #V 4 clear
V44sNi hcqmjqJ close all
` a1R "A Dm`U|<o n1=1,n2=1.45;
_$jJpy cE8 _keR~ theta=0:0.1:90;
(k HQK Qmq ~>-;(YU"t a=theta*pi/180;
1J% qbh V/[,1W[B rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
fm2M i~}0
0{ Ll4 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
=rtA{g$)+ 0 dX= tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
(R
2P<
Zr x [uXD ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^DCv-R+ p L" ^366M! figure(1)
(9q {J(44 Y GO 7lar subplot(1,2,1);
,ce^"yG -v+^x`HR plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
pxnUe1= x r-`i legend('r_p','r_s','|r_p|','|r_s|')
%_=R&m'n` r0uXMr=Z96 xlabel('\theta_i')
(c|qX-%rC -$|X\#R ylabel('Amplitude')
s=@Ce V@4W ~acK$.# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^3 s&90 M[N.H9 axis([0 90 -1 1])
:Ty*i KcF+!;: grid on
bvRGTOxO .`#R%4Xl subplot(1,2,2);
Gv3Fg[MA@c 6
&Aa b 56 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
k_ ^/ ^ nZ2p$ legend('t_p','t_s','|t_p|','|t_s|')
f0bV]<_9 |VEAzY|[# xlabel('\theta_i')
_Vf|F R2ZQBwB ylabel('Amplitude')
2 3 w{h d n L20}"$E title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
__%E!*m"<_ JJ3JULL2 axis([0 90 0 1])
tBUQf*B >^HTghgRD grid on
y%l#lz=6 eyjUNHeh# Rp=abs(rp).^2;
[t) i\ }V ryLNMh Rs=abs(rs).^2;
&V{,D))6[ ZTC1t_ Rn=(Rp+Rs)/2;
RteTz _ z{ B;je|M!d Tp=1-Rp;
jc`',o'[+ {S(T1ua Ts=1-Rs;
<s3( DA@hf Tn=(Tp+Ts)/2;
jn Y3G ^{bEq\5& figure(2)
^}\R]})w" K 8c #/o subplot(1,2,1);
^i1 :PlW] bj{f[nZ d plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
$ zi\ /Yw +Z]%@"S? legend('R_p','R_s','R_n')
kovzB] a/#,Y<kJ xlabel('\theta_i')
5)<jPyC t +q`h3 ylabel('Amplitude')
l);8y5 N)Kr4GC title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
zD-8#H35X" Ih"Ol(W axis([0 90 0 1])
[
Ulo; #P @1Lc`;Wd grid on
2oASz| ]~ UkD*Ct subplot(1,2,2);
S`KCVQ>V 9%6`ZS~3 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
.CpO+z h h}E Dnx legend('T_p','T_s','T_n')
-y<uAI g "r+<=JU>OV xlabel('\theta_i')
h_!"CF <n DL_ \luh ylabel('Amplitude')
eO G%6C%a CU_06A|} title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
&X(-C9'j oFp&j@`k8j axis([0 90 0 1])
it vdzPO KZNyp% q grid on
*[n^6 ) i[# Tn52D 
!o!04_ v%[mt` I 
;v\n[
