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

end

h=zeros(mode,2);

%相速度

figure(l)

for j=l:2

ifj==l

cp=cps;

color=,b,;

else

cp=cpa;

color二T;

end

for i=l:mode

cpp=cp(:,i);

in d=fi nd(cpp==0);

if ^isempty(ind)

h(i/j)=plot((fdO(ind(end)+l:end))/d_Q235/cpp(ind(end)+l:end),color);

else

h(i/j)=plot(fdO/d_Q235,cpp/color);

end

hold on

end

ifj==2

xlabel('f/(KHz)')

ylabel('C_{p}/(km-sA{-l})')

title('6mm钢板相速度频散曲线J

set(gca,,xtick,,(0::20)*le3/d_Q235,,xticklaber/(0::20)*le3/d_Q235)

xlim([0, 1000]);%

set(gca,'ylim,,[0 cp_max],'ytick,,(0:cp_max/le3)*le3,...

'yticklaber,0:cp_max/le3)

grid on

hSGroup = hggroup;%要在子对象构建之后构建，构建后立即使用，否则将失效 hAGroup = hggroup;

set(h(:/l)/parent,,hSGroup)

set(h(:,2)/parent,/hAGroup)

set(get(get(hSGroup,'A nn otationfLege ndlnformation'),..?

'Icon DisplayStyle'/on1);

set(get(get(hAGroup,'A nn otationfLege ndlnformatiorV),..?

'Icon DisplayStyle'/on1);

legend)'对称模式T反对称模式')

end

end

%群速度

figure(2)

for j=l:2

ifj==l

cp=cps;

color=,b,;

else

cp=cpa;

color=,r,;

end

for i=l:mode

cpp=cp(:,i);

in d=fi nd(cpp==0);

if ^isempty(ind)

fd=fdO(ind(end)+l:end)';

cpp=cpp(i nd(end)+l:end);

else

fd=fdO';

end

dcdf=diff(cpp)/dfd;

cg=cpp(l:end-l).A2./(cpp(l:end-l)-fd(l:end-l).*dcdf);

h(i,j)=plot(fd(l:end-l)/d_Q235,cg/color);

hold on

end

ifj==2

xlabel('f/(KHz)')

ylabel('C_{g}/(km-sA{-l})')

title('6mm钢板群速度频散曲线，) set(gca,,xtick,,(0::20)*le3/d_Q235,,xticklaber/(0::20)*le3/d_Q235) xlim([0,1000]);%

set(gca；ylim,/[0]*le3；ytick,/(0::*le3,