声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1029|回复: 2

[HHT] 关于导入数据的问题。。请教前辈一下。。

[复制链接]
发表于 2010-10-12 20:45 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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('幅值');

我导入的数据

我导入的数据

请问是不是程序参数出了问题??  

pigwenti.png
回复
分享到:

使用道具 举报

 楼主| 发表于 2010-10-13 09:45 | 显示全部楼层
哪位前辈帮忙解决一下啊。。谢谢啊。。
发表于 2011-4-16 16:38 | 显示全部楼层
回复 2 # youshen 的帖子

是不是数据有问题 我用的是N*1的
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-16 11:43 , Processed in 0.098687 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表