我这里是哪错了,各位大侠们进来给帮帮忙
源程序是这样:fs=1000;
N=1000;
n=0:N-1;
t=n/fs;
xn=randn(1,length(t));
x=sin(2*pi*120*t)+.*cos(2*pi*30*t+0.5*sin(2*pi*15*t))+xn;
figure;
plot(t,x);
=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;
=hhspectrum(imf(1:end-1, :));
=toimage(A,f,t,length(t));
figure;
set(gcf,'Color','w');
imagesc(t/N,,E);
set(gca,'YDir','normal')
colormap('jet')
colorbar;
xlabel('时间 t/s');
ylabel('频率 f/Hz');axis()
title('Hilbert-Huang Spectrum') ;
运行到这:
plot_imf(xd,t,imf);%emdfenjie;
??? Undefined function or method 'plot_imf' for input arguments of type 'double'.
别人的没问题,怎么我一运行就错了呢?
还有一个
>> =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)';
这个是 只要一遇到=hhspectrum(imf(1:end-1, :)); 就出现 f(i,:)=instfreq(an(i,:)',tt,l)'; 怎么回事啊?
{:{10}:} 振动王子 发表于 2013-3-28 13:43 static/image/common/back.gif
{:{07}:} 你肯定是没有=hhspectrum(imf(1:end-1, :))这个子程序,,我也遇到了这样的问题 华电机械 发表于 2013-5-7 16:36 static/image/common/back.gif
你肯定是没有=hhspectrum(imf(1:end-1, :))这个子程序,,我也遇到了这样的问题
是的,解决了,谢谢啊!!! bgpz2007 发表于 2013-5-9 11:40 static/image/common/back.gif
是的,解决了,谢谢啊!!!
我也是初学者,大家相互帮助吧。。 算是学习了!!{:{39}:}
页:
[1]