matlab生成曲线程序,Matlab绘制频散曲线程序代码(20210119130722).docx

Matlab绘制频散曲线

程序代码

-CAL-FENGHAI-(2020YEAR-YI

function disper

%绘制平板频散曲线

%tic

clc;clear;

cl=5790;%材料纵波波速(钢板)

cs=3200;%材料横波波速(钢板) dfd=*le3;

fdO=:dfd/le3:2O)*le3;%频厚积(MHz*mm) d_Q235二 6;

cps_mi n二2700;

cpa_min=100;

cp_max=10000;

mode=3;%绘制的模式数

precision=le-8;

cpa=zeros(length(fdO),mode);

cps=zeros(le ng th(fdO),mode);

for i=l:length(fdO)

fd=fdO(i);

[cpl2 n]=ss(cps_min/cp_max/fd/cl,cs,mode);

for j=l:n

cpl=cpl2(j,l);

cp2=cpl2(j,2);

cps(i,j)=serfe n(cpl,cp2,fctcl￡S'precisi on); end

[cpl2 n]=aa(cpa_min,cp_max/fd/cl/cs,mode); for j=l:n

cpl=cpl2(j,l);

cp2=cpl2(j,2);

cpa(ij)=aerfe n(cpbcp2,fd￡l‘cs,precisi on); end

<