ex01.m

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

clear
h=0.00276;cd=6100;cs=3300;
for w=1000:1000:2000000;
    i=1;j=1;
for c=100:100:9000
b=w/c;
kdt=sqrt(abs((w/cd)^2-b^2));
kst=sqrt(abs((w/cs)^2-b^2));
y1(i)=(kst^2-b^2)^2*tan(kst*h/2);
y2(i)=4*b^2*kst*kdt*tan(kdt*h/2);
 if abs(y1(i)+y2(i))<=0.000001
   y(w/1000,j)=c;
   j=j+1;
 end
   i=i+1;
end
end