作者dooksp (頹廢的人)
看板MATLAB
標題[討論]黃轉換 Spectrum of energy
時間Tue Jun 30 10:19:31 2015
各位前輩安安~~ 小弟我最近剛開始接觸黃轉換
遇到一個百思不得其解的問題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://webptt.com/m.aspx?n=bbs/MATLAB/M.1435630773.A.2AC.html
1F：→ s4300026: hht 只是將時域圖分割成數個時域圖, 不能畫的原因 06/30 13:14
2F：→ s4300026: 應該原本就不能畫(code 有錯), 跟hht 無關 06/30 13:15