function disper
%绘制平板频散曲线
%tic
clc;clear;
cl=5790;%材料纵波波速(钢板)
cs=3200;%材料横波波速(钢板)
dfd=0.01*1e3;
fd0=(0.01:dfd/1e3:20)*1e3;%频厚积(MHz*mm)
d_Q235=6;
cps_min=2700;
cpa_min=100;
cp_max=10000;
mode=3;%绘制的模式数
precision=1e-8;
cpa=zeros(length(fd0),mode);
cps=zeros(length(fd0),mode);
for i=1:length(fd0)
fd=fd0(i);
[cp12 n]=ss(cps_min,cp_max,fd,cl,cs,mode);
for j=1:n
cp1=cp12(j,1);
cp2=cp12(j,2);
cps(i,j)=serfen(cp1,cp2,fd,cl,cs,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);
评论0