利用matlab表示全反射下的s波与p波的位相跃变曲线

a_1=40:0.01:90;
n=1/1.5;
c_s=2*atand((-1)*(sind(a_1).^2.-n^2).^0.5./cosd(a_1))

c_p=2*atand((-1)*(sind(a_1).^2.-n^2).^0.5./(n^2*cosd(a_1)))
%axis ij
plot(a_1,c_p)
hold on
plot(a_1,c_s)
hold off

grid on
title('s波与p波全反射下的位相跃变')
xlabel('a_1角度变化')
ylabel('s波与p波的相位跃变')
text(53,-125,'p波')
text(65,-108,'s波')