利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
9fr&Yb=_o@ DcFV^8O& 1、光疏射向光密
S%a}ip& v]d?6g clear
1*eWvYo1 rsq'60 close all
R, (+NT$ hH|XtQ.n^ n1=1,n2=1.45;
:#"OCXr k^vmRe<lk theta=0:0.1:90;
zL[U; c+\Gd}IJq a=theta*pi/180;
* Kp ^al 9Rt(G_' rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
}!J/ 9WKgU % 'L= rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
.bY1N5=sz 'E4(!H,k tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ZeDDH a7 '\* ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
sjVl/t`l ES:p^/ =* figure(1)
m SvTnd8 hF{mm(qyv subplot(1,2,1);
k+>p!1 m<VL19o>R plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
ROI$;B( u{y5'cJ{ legend('r_p','r_s','|r_p|','|r_s|')
BfO}4 T=)L5 Vuq< xlabel('\theta_i')
W0C$*oe!_i &`9bGO ylabel('Amplitude')
h11bK'TIv =K&q;;h title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
j(2tbWg9- hJtghG6v axis([0 90 -1 1])
jind!@}! v(WL 3[y; grid on
61 HqBa kv`3Y0R-" subplot(1,2,2);
I[YfF F^[Rwzv>c plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
zJym`NF (j\UoKLRt legend('t_p','t_s','|t_p|','|t_s|')
Xwn|. B\~3p4S xlabel('\theta_i')
r;s3(@[,@ i=+<7]Q ylabel('Amplitude')
Pp" )hFx I5AjEp title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
_a02# L?pvz} axis([0 90 0 1])
ckkM)|kK Lw78v@dY grid on
=I*ZOE3n N_~Wu Rp=abs(rp).^2;
MDXQj5s^ MZ#2WP)F Rs=abs(rs).^2;
UHm+5%ZC Y K 62#; Rn=(Rp+Rs)/2;
UmHb-uk ; DL/*t.)"et Tp=1-Rp;
B=qRZA!DQ? ' [M2Q"X Ts=1-Rs;
{`SGB;ho
a}8>(jtSt Tn=(Tp+Ts)/2;
w2_I/s6B nT6iS}h figure(2)
Kh_>V m/ +/ d8d subplot(1,2,1);
.*..pf|/ MhXm-<4
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
<UHf7:0V o]k]pNO legend('R_p','R_s','R_n')
~V)VGGOL$v K>E!W!-PJ xlabel('\theta_i')
%.} jF2GHyB ylabel('Amplitude')
i}12mjF 5 s2}nIe title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Y
.X-8 {.[,ee-)9 axis([0 90 0 1])
N1zrfn-VU n!,TBCNX grid on
@G&xq"Fg7 2G*#Czr" subplot(1,2,2);
M(gWd8?# G@6,O-Sj plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
:BiR6>1: ]uXsl0'`V legend('T_p','T_s','T_n')
<Gt{(is C F 0IP xlabel('\theta_i')
?DTP-#5Ba QRix_2+ ylabel('Amplitude')
k:yu2dQh G~j<I/)" title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
J:!Gf^/) Z4tq&^ :c= axis([0 90 0 1])
jtJ8r5j 1 }Bg<Fm grid on
[Cr~gd+q /1@m#ZxA: \r&@3a.> (9lx5 {N`<e>A]{