计算光波由光密入射到光疏介质时相位变化
o;5 ns =d M'n}@U
clear
,\U c/w R /c:78@ close all
1@S(v L3a e=u ?-8 n1=1.45,n2=1;
!/RL.`!> :.bBV] 6q theta=0:0.1:90;
sf5 F$ ;\48Q; a=theta*pi/180;
=jIP 29+ G~ nQR
qv rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
*P0sl( & thK4@C|X4 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
'74*- y d >o,l/# z tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
=Hf`yH\# ' |Ia-RbX ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
8qF OO3c\V 5|_El/G arp=angle(rp);
$WICyI{$ Mv\]uAT` ars=angle(rs);
&X&msEM vMv?
fE" atp=angle(tp);
X 3q2XU %s* F~E ats=angle(ts);
(F=q/lK$ 1}E R+;If figure(1)
q:ah%x[ mGP&NOR0^y subplot(1,2,1);
~b;u1;ne WinwPn+9 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
L)yc_ d5 7Q>bJ Ek7 legend('arg(r_p)','arg(r_s)')
26B+q XEt !  * xlabel('\theta_i')
]gjB%R[.m 8'| _O ylabel('\phi')
'9XSz? b; jr;I title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
y*H rv ywY[g{4+ axis([0 90 -3.5 3.5])
ni~1)"U. &MJ cLM] grid on
!cNw 8"SIU 4#Cm5xAt6 subplot(1,2,2);
Cc&SHG*R }!*CyO* plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
CX3yIe~u , .]1N :
legend('arg(t_p)','arg(t_s)')
/ ei(Q'pc[ IrMxdF~c xlabel('\theta_i')
_;'<}a [oD u3Qn ylabel('\phi')
OKV/=]GS MHPh! title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
R4g;-Ci-> vR*p1Kq: axis([0 90 -3.5 3.5])
6Ij'z9nJw E'+?7ZGWj grid on
'LMMo4o3 2 F#DJN# 
