各位前輩安安~~ 小弟我最近剛開始接觸黃轉換
遇到一個百思不得其解的問題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)