利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
D# 9m\o_ P [G)sA_" 1、光疏射向光密
0I-9nuw,^; niMsQ clear
+ B,}Q r 60^`JVGWH close all
6fE7W>la e-})6)XgA n1=1,n2=1.45;
!,_u)4 K C*e/J theta=0:0.1:90;
PV.X z0@R 5zK4F raf a=theta*pi/180;
>mbHy<< XAD- 'i rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
D%[mW c@1I ih-#5M@ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
CCs%%U/= )J o: pkM tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<`8n^m* Y Vt% 0 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
(R,#a *CV nmee 'oEw figure(1)
\Gef \ r8t}TU>C subplot(1,2,1);
]6k\)# %2 hE:9{;Gf plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
9 &dtd WM Dl=6 legend('r_p','r_s','|r_p|','|r_s|')
>>4qJ%bL zF`0J xlabel('\theta_i')
&N$<e(K [Q~#82hBhY ylabel('Amplitude')
MVpGWTH@F X;+sUj8 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
9Z $"K- G IV~>I-rd axis([0 90 -1 1])
C$=%!w f q_: 4w$> grid on
3oj' y txN 4!{KWL `A subplot(1,2,2);
J'6PmPzY| tH@Erh|% plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
^cC,.Fdw l
K{hVqpt legend('t_p','t_s','|t_p|','|t_s|')
etDk35!h~, 1/B> XkCJ xlabel('\theta_i')
~Y[r`]X`"m >a<.mU|# ylabel('Amplitude')
AG
nxYV"p fC d&D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Uk wP 5:[0z5Hww axis([0 90 0 1])
3lL-)<0A( 5+0gR
&|j grid on
0@0w+&*"@ *->W^1eGM Rp=abs(rp).^2;
& GO}|W 'T;P;:!\ Rs=abs(rs).^2;
79rD7D&g IxN9&xa Rn=(Rp+Rs)/2;
kOrZv,qFG[ wYXQlxd y Tp=1-Rp;
un"Gozmt5 a#(?P .6 Ts=1-Rs;
B N5[,J w>&aEv/f Tn=(Tp+Ts)/2;
mkk6`,ov #4NaL figure(2)
`,*3[ oAVnK[EMq` subplot(1,2,1);
1p3z1_wrs y3Q sv plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
hp 5 0J e a2ayT legend('R_p','R_s','R_n')
u(.e8~s8 ,5p(T_V/ xlabel('\theta_i')
% g$o/A$ ,K s8*;#r ylabel('Amplitude')
uk:(pZ-uJ :K,i\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
cG zPI +F ,.8KN<A2]' axis([0 90 0 1])
dh iuI|?@ =U9*'EFr grid on
@ CL{D:d !X#OOqPr= subplot(1,2,2);
] IQ&>z}< # $07:UJ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
X=&ET)8-Y .p3,O6y2(F legend('T_p','T_s','T_n')
`:KY\ !sP {gi#= xlabel('\theta_i')
&-6Gc;f8 ;?i W%:_, ylabel('Amplitude')
20 h, ^ AM \'RHL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
BoWg0*5xb -zgI_u9=EB axis([0 90 0 1])
> uB# &Q ^A/k)x6 grid on
{$
JYw{a 5r|,CQ7o 
:Yh+>c}N ES7>H 
F^BS/Yag
