没看明白emd_online，请大虾指点。。
Rilling写的emd_online没看明白，在他的文章《ON EMPIRICAL MODE DECOMPOSITIONAND ITS ALGORITHMS》中提到ex_online.m可以帮助理解，但是也没找到这个，对于emd_online讲解瞒少，
程序读起来很吃力，用ecg.mat的数据输入： emd_online(ecg);
结果报错：??? Index exceeds matrix dimensions.
Error in ==> emd_online at 282
= extr(mps(k,max([(limpsl(k,i)-1),1]):finps(k,i),i));
可能是输入信号的问题，mps＝输入信号，而ecg是4170×1，不知道是不是输入参数有问题造成数组维数出问题，
另外请大虾讲解下emd_online这个程序的核心思想吧，感谢。。。。。。 刚找到了ex_online，准备来看～～～
看来eight说得对，要多利用搜索功能。。。:@$
如果仍然看不懂我再来问吧 在论坛上没找到ex_online，http://forum.vibunion.com/space/html/15/t-42415.html
有是有，但打开链接地时候：附件文件不存在或无法读入，请与管理员联系。
请管理员发一份过来吧，nos2006@163.com,谢谢～～
另外搜索了关键词emd_online，没有相关介绍哦
还得请大虾们指教，问题如上～ 原帖由 dudubear2009 于 2007-5-23 20:17 发表 http://www.chinavib.com/forum/images/common/back.gif
在论坛上没找到ex_online，http://www.chinavib.com/space/html/15/t-42415.html
有是有，但打开链接地时候：附件文件不存在或无法读入，请与管理员联系。
请管理员发一份过来吧，nos2006@163.com,谢谢～～
...
你说找不到 ex_online.m ？ 刚看了，每个版本都有，建议认真查找，实在不行的话，利用 windows 的搜索功能 online程序是建立在local的基础上的,你应该先看一下local.m,理解起来不难
新建EMD的qq群 主要是初学者的 希望大家参加 40652876
回复 #6 sharc 的帖子
你发一个帖子就可以了！这么多就不好了。:handshake 请以后注意。大家会看见的，放心好了 我最近也在看on-line emd程序，没看懂，急啊！！！哪位高手要是明白的话，请指点一下回复 7楼 zhangnan3509 的帖子
zhangnan3509 :你能否帮我看一下EX_ONLINE.m中的程序，我有个弱弱的问题，请你帮我解答一下程序如下：% elementary signal
N = 2000;% # of data samples
T = 1:4:N;
t = 1:N;
p = N/2;% period of the 2 sinusoidal FM's
% sinusoidal FM 1
fmin1 = 1/64;% min frequency
fmax1 = 1.5*1/8;% max frequency
x1 = fmsin(N,fmin1,fmax1,p,N/2,fmax1);
% sinusoidal FM 2
fmin2 = 1/32;% min frequency
fmax2 = 1.5*1/4;% max frequency
x2 = fmsin(N,fmin2,fmax2,p,N/2,fmax2);
% logon
f0 = 1.5*1/16;% center frequency
x3 = amgauss(N,N/2,N/8).*fmconst(N,f0);
a1 = 1;
a2 = 1;
a3 = 1;
x = real(a1*x1+a2*x2+a3*x3);
x = x/max(abs(x));
% periodized signal
y = ;
其中执行完x = x/max(abs(x));这个语句后，x是不是就是就是含有2000个数据的行向量了？然后y是含有8个列向量的矩阵？ 问题已经解决
回复 楼主 dudubear2009 的帖子
能问一下emd-online好干什么用马回复 11楼 haifeng1213 的帖子
我运行emd中的例子，怎么找不到FMSIN这个函数，忘指点 function =fmsin(N,fnormin,fnormax,period,t0,fnorm0,pm1);%FMSIN Signal with sinusoidal frequency modulation.
% =FMSIN(N,FNORMIN,FNORMAX,PERIOD,T0,FNORM0,PM1)
% generates a frequency modulation with a sinusoidal frequency.
% This sinusoidal modulation is designed such that the instantaneous
% frequency at time T0 is equal to FNORM0, and the ambiguity
% between increasing or decreasing frequency is solved by PM1.
%
% N : number of points.
% FNORMIN : smallest normalized frequency (default: 0.05)
% FNORMAX : highest normalized frequency (default: 0.45)
% PERIOD: period of the sinusoidal fm (default: N )
% T0 : time reference for the phase (default: N/2 )
% FNORM0: normalized frequency at time T0 (default: 0.25)
% PM1 : frequency direction at T0 (-1 or +1) (default: +1)
% Y : signal
% IFLAW : its instantaneous frequency law (optional).
%
% Example:
% z=fmsin(140,0.05,0.45,100,20,0.3,-1.0);plot(real(z));
%
% See also FMCONST, FMLIN, FMODANY, FMHYP, FMPAR, FMPOWER.
% F. Auger, July 1995.
% Copyright (c) 1996 by CNRS (France).
%
%This program is free software; you can redistribute it and/or modify
%it under the terms of the GNU General Public License as published by
%the Free Software Foundation; either version 2 of the License, or
%(at your option) any later version.
%
%This program is distributed in the hope that it will be useful,
%but WITHOUT ANY WARRANTY; without even the implied warranty of
%MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.See the
%GNU General Public License for more details.
%
%You should have received a copy of the GNU General Public License
%along with this program; if not, write to the Free Software
%Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA02110-1301USA
if (nargin == 0),
error ( 'At least 1 parameter required' ) ;
elseif (nargin == 1),
fnormin = 0.05; fnormax = 0.45; period = N;
t0= round(N/2); fnorm0 = 0.25; pm1=+1;
elseif (nargin == 2),
fnormax = 0.45; period = N;
t0 = round(N/2); fnorm0 = 0.5*(fnormin+fnormax); pm1=+1;
elseif (nargin == 3),
period = N; t0 = round(N/2);
fnorm0 = 0.5*(fnormin+fnormax); pm1=+1;
elseif (nargin == 4),
t0 = round(N/2); fnorm0 = 0.5*(fnormin+fnormax); pm1=+1;
elseif (nargin == 5),
fnorm0 = 0.5*(fnormin+fnormax); pm1=+1;
elseif (nargin == 6),
pm1= +1;
elseif (nargin==7),
if (abs(pm1)~=1),
error('pm1 must be equal to -1 or +1');
end;
end;
if (N <= 0),
error ('The signal length N must be strictly positive' );
elseif (abs(fnormin) > 0.5)|(abs(fnormax) > 0.5)|(abs(fnorm0) > 0.5),
error ('fnormin, fnormax and fnorm0 must be between -0.5 and 0.5') ;
elseif (fnormin > fnormax)
error ('fnormin must be lower than fnormax');
elseif (fnormin > fnorm0)|(fnorm0 > fnormax),
error ('fnorm0 must be between fnormin and fnormax') ;
else
fnormid=0.5*(fnormax+fnormin);
delta=0.5*(fnormax-fnormin);
phi=-pm1*acos((fnorm0-fnormid)/delta);
time=(1:N)-t0;
phase=2*pi*fnormid*time+delta*period*(sin(2*pi*time/period+phi)-sin(phi));
y = exp(j*phase).';
if (nargout==2)
iflaw=fnormid+delta*cos(2*pi*time'/period+phi);
end
end 楼上为广告帖请删除! 已禁止访问
一一报告还需等45秒! 没时间等:@Q
处理完记得也删除此帖
页:
[1]
2