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
<