HHT三维时频谱处理问题

qiuxianyangdedi

在用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);
[m,n]=size(imf);
%[A,fa,tt]=hhspectrum(imf);
if size(imf,1)>1
   [A,fa,tt]=hhspectrum(imf(1:end-1,:));
else
    [A,fa,tt]=hhspectrum(imf);
end
[im,tt1]=toimage(A,fa,tt,length(tt),fs/2);

for i=1:m
    faa=fa(i,:);
[FA,TT1]=meshgrid(faa,tt1);
figure(11);
surf(FA,TT1,im)
title('HHT显示三维时频谱')
end

kukusuiyi

怎么修改名字？？？

qiuxianyangdedi

kukusuiyi 发表于 2015-4-11 10:00
怎么修改名字？？？

改名字？什么意思？

airazor

这个是你输入的参数长度不一，看看是不是矩阵搞反了。
一是你搞错了，二是surf函数里面的要求，你看看它对输入要求。

解决你问题的办法是找一个网上的HHT三维谱程序，现在很少有人能修改程序的，大都用的现成的

武林低手

你的问题解决了吗

武林低手

qiuxianyangdedi 发表于 2015-4-11 10:57
改名字？什么意思？

s=wq(:,14);
imf=emd(s);
[m,n]=size(imf);
[A,fa,tt]=hhspectrum(imf(1:m-1,:));
[E,tt1]=toimage(A,fa,tt,length(tt));
for i=1:m-1;
faa=fa(i,:);
[FA,TT1]=meshgrid(faa,tt1);%三维图显示HHT时频图
figure(11);
surf(TT1,FA,E);
hold on
title('HHT时频谱三维显示');
end
这段程序，不会导致你出现的问题，但是出来的时间和瞬时频率不是实际值，不知道你自己是怎么解决的

武林低手

airazor 发表于 2015-6-13 19:50
这个是你输入的参数长度不一，看看是不是矩阵搞反了。
一是你搞错了，二是surf函数里面的要求，你看看它对 ...

你好，你研究过HHT吗

cufflink

今天时频三维图没有画出来，我感觉不仅仅是画图的问题，还有一些问题没搞清楚。

airazor

武林低手 发表于 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、图中出现的蓝色平面是怎么回事

airazor

http://forum.vibunion.com/forum. ... mp;page=1#pid805271

airazor

还是QQ吧，我已经发送给你了

airazor

http://forum.vibunion.com/forum. ... mp;page=1#pid805271

luogangtju

我用的是LMD+HHT时频分析，但是也很有问题，用的相同数据长度，截取相同长度不同的数据部分，得到的瞬时频率随电流的变化也很问题