function disper
%绘制平板频散曲线
%tic
clc;clear;
df=0.01*1e3;
f0=(0.01:df/1e3:20)*1e6;%定义频率范围
cps_min=2700;
%cpa_min=100;
cp_max=10000;
mode=5;%绘制的模式数
precision=1e-8;
%cpa=zeros(length(fd0),mode);
cps=zeros(length(f0),mode);
for i=1:length(f0)
f=f0(i);
[cp12 n]=ss(cps_min,cp_max,f, mode);
for j=1:n
cp1=cp12(j,1);
cp2=cp12(j,2);
cps(i,j)=serfen(cp1,cp2,f, precision);
end
[cp12 n]=aa(cpa_min,cp_max,fd,cl,cs,mode);
for j=1:n
cp1=cp12(j,1);
cp2=cp12(j,2);
cpa(i,j)=aerfen(cp1,cp2,fd,cl,cs,precision);
end
end
h=zeros(mode,2);
%相速度
figure(1)
for j=1:2
if j==1
cp=cps;
color='b';
else
cp=cpa;
color='r';
end
for i=1:mode
cpp=cp(:,i);
ind=find(cpp==0);
if ~isempty(ind)
h(i,j)=plot((fd0(ind(end)+1:end))/d_Q235,cpp(ind(end)+1:end),color);
else
h(i,j)=plot(fd0/d_Q235,cpp,color);
评论3