Matlab钢板频散曲线程序代码(修改声速可绘制其他材料板材的频散曲线).docx

节选段落一：
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


节选段落二：
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);
end
hold on
end
if j==2
xlabel('f/(KHz)')
ylabel('C_{p}/(km·s^{-1})')
title('6mm钢板 相速度频散曲线


节选段落三：
*dcdf);
h(i,j)=plot(fd(1:end-1)/d_Q235,cg,color);
hold on
end
if j==2
xlabel('f/(KHz)')
ylabel('C_{g}/(km·s^{-1})')
title('6mm钢板 群速度频散曲线')
set(gca,'xtick',(0:0.6:20)*1e3/d_Q235,'xticklabel',(0:0.6:20)*1e3/d_Q235)
xlim([0, 1000]);%