lamea0.m

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

%求解Lame波反对称部分的基本模态Ao
format long
syms c;   %定义速度c为符号函数
kh=1.5    %k为波数，h为总的板厚
cd=5790;cs=3200; %cd为纵波波速，cs为横波波速
y1=4*tanh(kh/2*(1-(c/cs)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2);
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cd)^2)^(1/2));
ezplot(y1,[0,3500]); %作图求解
hold on
grid on
ezplot(y2,[0,3500]);
[xx,yy]=ginput(1);
[x(1),fval,exitflag]=fzero(inline(y1-y2),xx(1),[]);
c=x(1)
y1=4*tanh(kh/2*(1-(c/cs)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2);
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cd)^2)^(1/2));
e1=y1-y2