利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
/DQaGq/Ld Sz- J y:j 1、光疏射向光密
Oy U Bp8'pj;~ clear
(q+)'H%iK 6f')6X'x close all
| 7yAX+ $YN6<5R) n1=1,n2=1.45;
r\qj! OQW#a[=WQ theta=0:0.1:90;
?7<JQh)"e S;$-''o?9 a=theta*pi/180;
s l]_M t2,A@2DU 2 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
QFYWA1<pDh ~0beuK&p rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
d
,4]VE &boOtl^
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
N"nd*? xx[9~z=d ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
ZovW0 Q)m O8 B\{T1 figure(1)
$ s9Vrw0Z dZi "$ g subplot(1,2,1);
Y.g59X!Ub2 :<}=e@/~| plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
*:ZDd g !rQ4#4 legend('r_p','r_s','|r_p|','|r_s|')
xuq v 6b. NR` C(^} xlabel('\theta_i')
{NHdyc$ SQX:7YF~ ylabel('Amplitude')
rg^'S1x| `DV.+>O-1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
SHe49!RA'{
TWA-.>c axis([0 90 -1 1])
V5UF3'3;} _f$^%?^ grid on
_d5QbTe i\ ,-o O subplot(1,2,2);
N@t|7~ Wk)OkIFR plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
R}O_[ DGS $Ukz&T legend('t_p','t_s','|t_p|','|t_s|')
vN`klDJgW[ 8fl`r~bqZ xlabel('\theta_i')
n*2UnKaJ #ZB~ x6i6 ylabel('Amplitude')
kqFP)!37 wB.&}p9p title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
02 c ':a=7 3kybLO G axis([0 90 0 1])
W"k"I vTW} <J) ]mh dm grid on
As'=tIro o]4*|ARPs Rp=abs(rp).^2;
| ?,A]|j sB7#
~p A Rs=abs(rs).^2;
.+$ Q <L 8WXQ Oo8 Rn=(Rp+Rs)/2;
Sc;BCl {=| .(K)?r-g5 Tp=1-Rp;
OCe!.` KVc lhT<F Ts=1-Rs;
fp`;U_-&0 k>;`FFQU> Tn=(Tp+Ts)/2;
F1*>y ZOh`(})h y figure(2)
r>>%2Z-P Mk"^?%PxT subplot(1,2,1);
`dq,>HdW %)1y AdG
8 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
~%<X0s| 8\+uec]k legend('R_p','R_s','R_n')
-t!~%_WCv m| n xlabel('\theta_i')
< ^#,_o,! ~vm%6CABM ylabel('Amplitude')
]cHgleHQ t}r ' k/[ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
f6hnTbJ marQ NZ axis([0 90 0 1])
V?6a 8lJ -r `. #c4 grid on
gb[5&> (# b RFLcM subplot(1,2,2);
3l rT3a3vV 'j#*6xD plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
~Y^+M* n] . _uza legend('T_p','T_s','T_n')
n(1l}TJy }T(D7| ^R xlabel('\theta_i')
&* M!lxDN T<n ylabel('Amplitude')
X _q\S g ,0 M_ Bk" title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
n@<YI 03$mYS_? axis([0 90 0 1])
`V}q-Z dy ?+))}J5N\ grid on
Z F!h<h&, 0"jY.*_EW 
Wf+cDpK ~[t[y~Hup 
3#LlDC_WC
