我这里是哪错了，各位大侠们进来给帮帮忙

bgpz2007

源程序是这样：
fs=1000;
N=1000;
n=0:N-1;
t=n/fs;
xn=randn(1,length(t));
x=sin(2*pi*120*t)+[1+0.2*sin(2*pi*7.5*t)].*cos(2*pi*30*t+0.5*sin(2*pi*15*t))+xn;
figure;
plot(t,x);
[thr,sorh,keepapp]=ddencmp('den','wv',x);
xd=wdencmp('gbl',x,'db3',2,thr,sorh,keepapp);
figure;
plot(t,xd);
imf=emd(xd);
plot_imf(xd,t,imf);%emdfenjie;
[A,f,t]=hhspectrum(imf(1:end-1, :));
[E,t,Cenf]=toimage(A,f,t,length(t));
figure;
set(gcf,'Color','w');
imagesc(t/N,[0,0.5*fs],E);
set(gca,'YDir','normal')
colormap('jet')
colorbar;
xlabel('时间 t/s');
ylabel('频率 f/Hz');axis([0 1 0 200])
title('Hilbert-Huang Spectrum') ;



运行到这：
plot_imf(xd,t,imf);%emdfenjie;
??? Undefined function or method 'plot_imf' for input arguments of type 'double'.
别人的没问题，怎么我一运行就错了呢？


还有一个
>> [A,f,t]=hhspectrum(imf(1:end-1, :));
??? Undefined function or method 'instfreq' for input arguments of type 'double'.
Error in ==> hhspectrum at 79
    f(i,:)=instfreq(an(i,:)',tt,l)';
这个是 只要一遇到  [A,f,t]=hhspectrum(imf(1:end-1, :));        就出现     f(i,:)=instfreq(an(i,:)',tt,l)';    怎么回事啊？

振动王子

bgpz2007

振动王子 发表于 2013-3-28 13:43

                               
登录/注册后可看大图

华电机械

你肯定是没有  [A,f,t]=hhspectrum(imf(1:end-1, :))这个子程序，，我也遇到了这样的问题

bgpz2007

华电机械 发表于 2013-5-7 16:36

                               
登录/注册后可看大图

你肯定是没有  [A,f,t]=hhspectrum(imf(1:end-1, :))这个子程序，，我也遇到了这样的问题

是的，解决了  ，谢谢啊！！！

华电机械

bgpz2007 发表于 2013-5-9 11:40

                               
登录/注册后可看大图

是的，解决了  ，谢谢啊！！！

我也是初学者，大家相互帮助吧。。

牛小贱

算是学习了！！