各位前辈安安~~ 小弟我最近刚开始接触黄转换
遇到一个百思不得其解的问题QQ
我使用EEMD分析出了11个IMF
接下来想要使用nnspe画出时频图
但始终跑出以下Error:
Error in fa (line 181)
[f, a] = FAimpHilbert(data,dt);
Error in nnspe (line 200)
[freq,amp] = fa(data,dt,ifmethod,normmethod,nfilter);
以下是我的code恳请各位指点一下clear all
samplerate=500;
N=1024;
sinefreq1=30;
sinefreq2=100;
noverlap=round(nfft*0.5);
ts=[1:N/4]/samplerate; % time axis for each segment
data=[zeros(1,N/4) sin(2*pi*sinefreq1*ts) sin(2*pi*sinefreq2*ts)
zeros(1,N/4)];
taxis=[1:N]/samplerate; % time axis for whole data length
figure,
plot(taxis,data),xlim([taxis(1) taxis(end)])
xlabel('Time (s)')
noiselevel=0;
Nensemble=1;
EMDIMF=eemd(data,noiselevel,Nensemble);
figure(2)
anoiselevel=num2str(noiselevel);
aNensemble=num2str(Nensemble);
strips(EMDIMF);
title(['EMD result ,noise level=',anoiselevel,' Number of
ensemble=',aNensemble]);
[nte,tae,fae]=nnspe(EEMDIMF)