HHT三维时频谱处理问题
在用surf函数显示HHT三维时频图时出现下列错误：??? Error using ==> surf at 78
Data dimensions must agree.
Error in ==> sanweishipinpu at 22
surf(FA,TT1,im)
请高手帮忙解决，程序如下：
clc;
clear all;
N=11000;
t=0.0001*(1:N);
fs=10000;
xx=load('sybx.csv');
x=xx(:,2)*1E3;
imf=emd(x);
=size(imf);
%=hhspectrum(imf);
if size(imf,1)>1
=hhspectrum(imf(1:end-1,:));
else
=hhspectrum(imf);
end
=toimage(A,fa,tt,length(tt),fs/2);
for i=1:m
faa=fa(i,:);
=meshgrid(faa,tt1);
figure(11);
surf(FA,TT1,im)
title('HHT显示三维时频谱')
end
怎么修改名字？？？ kukusuiyi 发表于 2015-4-11 10:00
怎么修改名字？？？
改名字？什么意思？ 这个是你输入的参数长度不一，看看是不是矩阵搞反了。
一是你搞错了，二是surf函数里面的要求，你看看它对输入要求。
解决你问题的办法是找一个网上的HHT三维谱程序，现在很少有人能修改程序的，大都用的现成的 你的问题解决了吗 qiuxianyangdedi 发表于 2015-4-11 10:57
改名字？什么意思？
s=wq(:,14);
imf=emd(s);
=size(imf);
=hhspectrum(imf(1:m-1,:));
=toimage(A,fa,tt,length(tt));
for i=1:m-1;
faa=fa(i,:);
=meshgrid(faa,tt1);%三维图显示HHT时频图
figure(11);
surf(TT1,FA,E);
hold on
title('HHT时频谱三维显示');
end
这段程序，不会导致你出现的问题，但是出来的时间和瞬时频率不是实际值，不知道你自己是怎么解决的 airazor 发表于 2015-6-13 19:50
这个是你输入的参数长度不一，看看是不是矩阵搞反了。
一是你搞错了，二是surf函数里面的要求，你看看它对 ...
你好，你研究过HHT吗 今天时频三维图没有画出来，我感觉不仅仅是画图的问题，还有一些问题没搞清楚。 武林低手 发表于 2015-6-15 16:58
你好，你研究过HHT吗
研究不敢当，把图画出来还是可以。 airazor 发表于 2015-6-17 16:47
研究不敢当，把图画出来还是可以。
加个扣扣，交流一下撒，扣扣：710482955 airazor 发表于 2015-6-17 16:47
研究不敢当，把图画出来还是可以。
我用的是西储大学的数据，采样频率42000hz,这是我用hht画出的三维时频图，这是正常轴承，故障样本图未给出。有几个个问题:1、瞬时频率出现12000hz，故障样本的也会出现这么大的，我感觉太大了，我怀疑是不是哪里出错了（震动信号用hht之前用小波变换滤波） 2、图中出现的蓝色平面是怎么回事 http://forum.vibunion.com/forum.php?mod=viewthread&tid=136440&page=1#pid805271 还是QQ吧，我已经发送给你了 http://forum.vibunion.com/forum.php?mod=viewthread&tid=136440&page=1#pid805271 我用的是LMD+HHT时频分析，但是也很有问题，用的相同数据长度，截取相同长度不同的数据部分，得到的瞬时频率随电流的变化也很问题
页:
[1]