利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
a?Fz& BE Z .`+IN(>E 1、光疏射向光密
#pRbRT9 n. N0Nhd clear
rk=w~IZJ3 J qLPJUr close all
zz3{+1w] xM)P=y_!M+ n1=1,n2=1.45;
Wx/PD=Sf& lz0dt<8eP theta=0:0.1:90;
W{JR%Sq$ / tkV/ a=theta*pi/180;
]p( es,[ qtVgjT2#H rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
68~]_r.a 'GW~~UhdW rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
}c9RDpjh~ E\4ZUGy 0 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
@?tR-L<u 8}oe))b ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
: [?7,/w e#6H[t
figure(1)
ve/.q^JeJ + tkm,> s subplot(1,2,1);
#m7evb5eg* N["M "s(N plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
j5hQ;~Fa| R90chl legend('r_p','r_s','|r_p|','|r_s|')
JvT#Fxj k ]$)};8;7W xlabel('\theta_i')
)MN 6\v qoQ,3 &< ylabel('Amplitude')
ak}k e %MtaWZ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h /aG."U s *CBYzOm axis([0 90 -1 1])
q2Gm8>F1y. IH=% %AS grid on
9Z2aFW9 sN[<{;K4 subplot(1,2,2);
4 [r:DM|8 vKbGG plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
B& 3@b .Pe^u%J6F legend('t_p','t_s','|t_p|','|t_s|')
0}_1 ZU Kv5 !cll5 xlabel('\theta_i')
FGMYpapc~ l>Zp#+I- ylabel('Amplitude')
I*+*Wf }z-)!8vF title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
g{?{N ew Yk> axis([0 90 0 1])
>Co@K^' 5zJ#d}%}S" grid on
d$D3iv^hyx u"q!p5P%q Rp=abs(rp).^2;
L ~M6 ca" #=fd8} 9 Rs=abs(rs).^2;
y9GaxW* & ) Vy}oFT\ Rn=(Rp+Rs)/2;
6~Zq s\
]Rgi>w Tp=1-Rp;
9jUm 0B{? {M0pq3SL*t Ts=1-Rs;
wGpw+O JR /^Go$^ Tn=(Tp+Ts)/2;
^ @q $c {<^P YN>` figure(2)
5r\Rfma f,0oCBLP O subplot(1,2,1);
03y5$ kQ ' l'[ U plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
p~u11rH C Y8=prC legend('R_p','R_s','R_n')
1QDAfRx ld23 ^r xlabel('\theta_i')
?37Kc,o 8!dA1]2; ylabel('Amplitude')
2fJ2o[v dzc.s8T(0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}"9 jCxXL .S Sj=q4 ? axis([0 90 0 1])
!TuMrA * GfT`>M?QGK grid on
& AlX). !~aDmY 2 subplot(1,2,2);
k*xgF[ T
8 d\aU rsPn plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
?)#}Nj<R &mp@;wI6@ legend('T_p','T_s','T_n')
JS1''^G&. lo!_;`v=U xlabel('\theta_i')
owc#RW9 7 zpd Z. ylabel('Amplitude')
du4Q^-repC 'Sjcm@ILm title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Cy##+u,C pVm'XP axis([0 90 0 1])
t@EHhiBz p2~MJ
LK4 grid on
_[.3I1kG @RF !p 
JB= L\E} :X;' 37o#q 
'PrrP3lO_~
