反射p波相位变化如下,横轴入射角。
所以在临界角之前,相位变化只有可能是0和pi,在临界角以后,全反射现象,折射角取复数,菲涅耳公式可以看成对应输入输出的传递函数。
或者说,rs,rp等本来就体现了输入输出之间的传递函数关系。
clear all;
n1=1.5;
n2=1;
theta1=0:0.01:pi/2;
theta2=asin((n1/n2)*sin(theta1));
rs=-sin(theta1-theta2)./sin(theta1+theta2);
ts=2*sin(theta2).*cos(theta1)./sin(theta1+theta2);
rp=tan(theta1-theta2)./tan(theta1+theta2);
delta=acos(real(rp)./(real(rp).*real(rp)+imag(rp).*imag(rp)));
% plot(theta1,sqrt(real(ts).*real(ts)+imag(ts).*imag(ts)))
% hold on;
% plot(theta1,ts)
% hold off;