[討論]黃轉換 Spectrum of energy
各位前輩安安~~ 小弟我最近剛開始接觸黃轉換
遇到一個百思不得其解的問題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)
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 61.70.78.141
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1435630773.A.2AC.html
→
06/30 13:14, , 1F
06/30 13:14, 1F
→
06/30 13:15, , 2F
06/30 13:15, 2F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
