lames.m

基于瑞利波及lamb波的频散曲线程序,可用于根据实验绘制Rali波及lamb波的频散曲线,检验其在不同材料中的频散特性

format long
syms p;
kh=6.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;
ezplot(y1,[0,1.5]);
grid on
hold on
ezplot(y2,[0.,1.5]);
zoom on
pause;
[xx,yy]=ginput(3);  %从图上获取曲线零点所对应的p值点
if xx(1)>=1
   c1=sqrt(xx(1)^2+1)  %求出所对应的c/cs值
else
   c1=sqrt(-xx(1)^2+1)
end
if xx(2)>=1
   c2=sqrt(xx(2)^2+1)  %求出所对应的c/cs值
else
   c2=sqrt(-xx(2)^2+1)
end
if xx(1)>=1
   c3=sqrt(xx(3)^2+1)  %求出所对应的c/cs值
else
   c3=sqrt(-xx(3)^2+1)
end
[x(1),fval,exitflag]=fzero(inline(y1-y2),xx(1),[]);
[x(2),fval,exitflag]=fzero(inline(y1-y2),xx(2),[]);
[x(3),fval,exitflag]=fzero(inline(y1-y2),xx(3),[]);
if x(1)>=1
   c11=sqrt(x(1)^2+1)  %求出所对应的c/cs值
else
   c11=sqrt(-x(1)^2+1)
end
if x(2)>=1
   c22=sqrt(x(2)^2+1) %求出所对应的c/cs值
else
   c22=sqrt(-x(2)^2+1)
end
if x(3)>=1
   c33=sqrt(x(3)^2+1)  %求出所对应的c/cs值
else
   c33=sqrt(-x(3)^2+1)
end
%检验结果
if c11>=1
   p1=sqrt(c11^2-1);
else
    p1=sqrt(-c11^2+1);
end  
y11=-4*tan(kh/2*((p1^2+1)*m^2-1)^(1/2))*p1*((p1^2+1)*m^2-1)^(1/2);
y21=(p1^2-1)^2*tan(kh/2*p1);
y01=y11-y21
if c22>=1
   p2=sqrt(c22^2-1);
else
    p2=sqrt(-c22^2+1);
end 
y12=-4*tan(kh/2*((p2^2+1)*m^2-1)^(1/2))*p2*((p2^2+1)*m^2-1)^(1/2);
y22=(p2^2-1)^2*tan(kh/2*p2);
y02=y12-y22
if c33>=1
   p3=sqrt(c33^2-1);
else
    p3=sqrt(-c33^2+1);
end 
y13=-4*tan(kh/2*((p3^2+1)*m^2-1)^(1/2))*p3*((p3^2+1)*m^2-1)^(1/2);
y23=(p3^2-1)^2*tan(kh/2*p3);
y03=y13-y23