利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
}@^4,FKJ d`/tE?Gw 1、光疏射向光密
? mhs$g> tu5*Qp\ clear
d%0~c'D8a vC5n[0 close all
5A4&+rdU Y9`5G% n1=1,n2=1.45;
$/7pYl\n pm6>_K z theta=0:0.1:90;
:Pv*, qHE c-Pw]Ju a=theta*pi/180;
c?%(Dp E Dxk+P!!K rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
! XzF67 Z?O *'#yn rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
6A WKLFMV j6g[N4xr tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
oaK%Ww6~ c hEn |>~ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
W y%s1iu y7|P-3[ 4w figure(1)
SM^-Z|d? &< !U fa& subplot(1,2,1);
JXj8Br?Z@ G!XIc>F* plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
_fwb!T}$ 36n>jS& legend('r_p','r_s','|r_p|','|r_s|')
.&x}NYX4 )nd\7|5# xlabel('\theta_i')
X7g3 Rtjqx6-B; ylabel('Amplitude')
ZKdeB3D & V)6!,rb title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0 n/gd"M 9Mgq1Z axis([0 90 -1 1])
NxLXm, ?r2 #.W grid on
2z@\R@F 1lpwZ" subplot(1,2,2);
L.=w?%:H= )$Z=t-q plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
@EoZI~
E~kG2x{a legend('t_p','t_s','|t_p|','|t_s|')
^xZ
e2@ d;+[i xlabel('\theta_i')
Z,X'-7YkU W<< 9y ylabel('Amplitude')
5p
+ZD7jK nL?oTze*p title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
kHv[H]+v P%MfCpyj axis([0 90 0 1])
_e7-zg$/ P5gN #G grid on
:h1pBEiH Eaqca{%/^ Rp=abs(rp).^2;
[0UGuj Pe2w sR"_U Rs=abs(rs).^2;
Cb:gH}j RJeSi`19T) Rn=(Rp+Rs)/2;
-<#!DjV6( Ap$y%6 Tp=1-Rp;
tW;:- **>/}.%?K Ts=1-Rs;
\ sA*V%n mw^7oO# Tn=(Tp+Ts)/2;
{w
<+_++ vJ`.iRU | figure(2)
Xe+,wW3YF jn.C|9/mj subplot(1,2,1);
7 }_! #Z~C`n
u plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
_u}4j 9T `O6:t\d@ legend('R_p','R_s','R_n')
=?X$Yaw* ]Zf6Yw .Y xlabel('\theta_i')
4eH.9t <:|3rfm# ylabel('Amplitude')
O3o: qly! 8I,QD`
xu title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3CE[( N:" C+ a( axis([0 90 0 1])
o K 9' `)4a[thp grid on
s4H2/EC MP>n)!R[` subplot(1,2,2);
q*)+ K9LRk _KD5T4FZR plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
~svO*o Wa sm QVWs> legend('T_p','T_s','T_n')
JmpsQ,, #&,H"?" xlabel('\theta_i')
5f:DN\ ] <a le $[ ylabel('Amplitude')
r gcWRt M* pRv title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:1 )DqoAJ P6U%=xaC axis([0 90 0 1])
[XlB<P=|> TQ9D68
, grid on
<EO<x D=: #:ns64| 
!~ -6wN"k 8X&Ya = 
h F&}lPVtv
