完整的EMD分解全过程，有Hilbert谱和边际谱

mning

运行提示：??? Undefined function or method 'emd' for input arguments of type 'double'.
emd.m哪里有

yangshunli

你的那个时间横轴是个采样点个数，不是时间啊，、如果用时间表示如何表示，单位用s.按你的程序应该是n/fs.但是如何加入到程序并显示出来，怎么弄啊？

a240612525

我是直接复制粘贴的，运行的时候为什么会出现下面的情况：
??? Error: File: D:\MATLAB\toolbox\hht2\toimage.m Line: 59 Column: 1
This statement is not inside any function.
(It follows the END that terminates the definition of the function "toimage".)
这个问题出在哪里啊

a240612525

上面的问题已经解决了，但是又有问题了：
??? Error using ==> toimage
Too many output arguments.
这是怎么回事？

a240612525

没问题了

z755924843

我这显示 没有EMD函数啊 楼主

雪缘

z755924843 发表于 2010-9-27 19:39

                               
登录/注册后可看大图

我这显示 没有EMD函数啊 楼主

需要安装emd工具箱，具体搜索一下本站

一夜孤舟100

18 19楼的问题我也遇到了 怎么修改啊 请指导一下

glwh

谢谢楼主，对我的帮助很大。先拷下了慢慢学

zhangjiuyun

我的也是在matlab直接运行出现
??? Undefined function or method 'instfreq' for input arguments of type 'double'.

Error in ==> hhspectrum at 79
  f(i,:)=instfreq(an(i,:)',tt,l)';

我把源程序拷过来，请大哥们帮帮忙。
%HHSPECTRUM  compute Hilbert-Huang spectrum
%
% [A,f,tt] = HHSPECTRUM(x,t,l,aff) computes the Hilbert-Huang spectrum
%
% inputs:
%   - x   : matrix with one signal per row
%   - t   : time instants
%   - l   : estimation parameter for instfreq (integer >=1 (1:default))
%   - aff : if 1, displays the computation evolution
%
% outputs:
%   - A   : instantaneous amplitudes
%   - f   : instantaneous frequencies
%   - tt  : truncated time instants
%
% calls:
%   - hilbert  : computes the analytic signal
%   - instfreq : computes the instantaneous frequency
%   - disprog : displays the computation evolution
%
%Examples:
%
%s = randn(1,512);
%imf = emd(s);
%[A,f,tt] = hhspectrum(imf(1:end-1,:));
%
%s = randn(10,512);
%[A,f,tt] = hhspectrum(s,1:512,2,1);
%
% rem: need the Time-Frequency Toolbox (http://tftb.nongnu.org)
%
% See also
%  emd, toimage, disp_hhs
%
% G. Rilling, last modification 3.2007
% gabriel.rilling@ens-lyon.fr

function [A,f,tt] = hhspectrum(x,t,l,aff)

error(nargchk(1,4,nargin));

if nargin < 2

  t=1:size(x,2);

end

if nargin < 3

  l=1;

end

if nargin < 4

  aff = 0;

end

if min(size(x)) == 1
        if size(x,2) == 1
                x = x';
                if nargin < 2
                        t = 1:size(x,2);
                end
        end
        Nmodes = 1;
else
        Nmodes = size(x,1);
end

lt=length(t);

tt=t((l+1):(lt-l));

for i=1:Nmodes

  an(i,:)=hilbert(x(i,:)')';
f(i,:)=instfreq(an(i,:)',tt,l)';----------就是这显示错误。
  A=abs(an(:,l+1:end-l));

  if aff
        disprog(i,Nmodes,max(Nmodes,100))
  end

end

hnlgzzk

楼主 你知道不知道 学习 HHT 这方面的matlab仿真书啊  给介绍介绍啊  
秋秋 ：494197178

youshen

非常好啊。。

wenyan1314

我的也是在matlab直接运行出现
??? Undefined function or method 'instfreq' for input arguments of type 'double'.

Error in ==> hhspectrum at 79
  f(i,:)=instfreq(an(i,:)',tt,l)';

这个问题到底是怎么解决的 请帮忙呀 急！

secondye

楼主这里是什么意思呢？？
%N=length(Cenf);%设置频率点数

另，不知道楼主有没有做过间断频率的设置呢？设置间断频率f，把频率高于f的信号分离出来

ncok123

谢谢楼主，看看