|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我导入一组数据,然后就会出来这样的错误。好像意思是索引超过了超出了模型范围。我用的是论坛里有位前辈的程序。
tspan=30; %时长
fs=50; %采样频率
t=(1/fs:1/fs:tspan);
N=length(t); %采样点数
x=data(1:N,2048);
figure(10);
plot(t,x);
[imf,ort,nbits]=emd(x); %对信号进行EMD分解
%emd_visu(x,1:length(x),imf)
%%% 求边际谱
[A,fa,tt]=hhspectrum(imf(1:end-1,:)); % A--瞬时幅值,fa--瞬时频率,tt--
[E,tt1,ff]=toimage(A,fa,tt,length(tt));%% ff--中心频率
错误所在
disp_hhs(E);colorbar;
% colormap(flipud(gray)) %显示灰度图
for k=1:size(E,1)
bjp(k)=sum(E(k,:))*1/fs*1/tspan;
end
f=(0:N-3)/N*(fs/2); % fs/(2*N)是时频图上y方向的分辨率
figure(60)
plot(f,bjp); %%% 作边际谱图
xlabel('频率 / Hz');
ylabel('幅值');
我导入的数据
请问是不是程序参数出了问题??
|
|