www.pudn.com > pinsanquxian.rar > lamea0.m, change:2004-08-05,size:572b


%求解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