平板频散曲线.rar_lamb平板_lamb波频散曲线_导波 MATLAB_导波频散曲线_频散

%%%%%画频散曲线的主程序%%%%% clear all; close all; format long; c = 480; %扫描相速度的点数 f = 300; %扫描频率的点数 f0 = 2e4; %频率步长 ic = 1; jf = 1; tic for F = [1 : f]*f0 %扫描频率，步长2k w = 2*3.14*F; for Cp = [1 : c]*20 + 1000 %扫描相速度，步长20 E(ic,jf) = calsingledet(Cp,w); ic = ic + 1; end ic = 1; jf = jf + 1; end t1 = toc tic CF = zeros(c,f); %储存每个频率下行列式波谷对应的相速度值 for jf = 1:f w = 2*3.14*jf*f0; for ic = 1:c - 2 if(E(ic+1,jf) - E(ic,jf) <0 & E(ic+1,jf) - E(ic+2,jf)< 0) CF(ic+1,jf) = search_zero(ic*20+1000,(ic+2)*20+1000,w);%将找到的相速度点放入CF中 end end end t2 = toc for ic = c:-1:1 plot((1:f),CF(ic,:),'.'); hold on; end