lame1.m

来自「基于瑞利波及lamb波的频散曲线程序,可用于根据实验绘制Rali波及lamb波的」· M 代码 · 共 7 行

function y=lam1(p)
syms p;
kh=1.0;
cd=6100;cs=3300;m=cs/cd;
y1=-4*tan(kh/2*((p^2+1)*m^2-1)^(1/2))*p*((p^2+1)*m^2-1)^(1/2);
y2=(p^2-1)^2*tan(kh/2*p);
y=y1-y2;