ex1.m

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

function y=ex1(c) 
h=0.00276;cd=6100;cs=3300;j=1
for w=1000000:1000000:1000000;
    for c=1000;
    b=w/c;
     kdt=sqrt((w/cd)^2-b^2);
     kst=sqrt((w/cs)^2-b^2);
     y=tan(kst*h/2)/tan(kdt*h/2)+4*b^2*kst*kdt/(kst^2-b^2)^2;
     if y<=0.001
        z(j)=c; 
        j=j+1;
     end
     end
end