|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
源程序是这样:
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)'; 怎么回事啊?
|
|