作者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/cn.aspx?n=bbs/MATLAB/M.1435630773.A.2AC.html
1F：→ s4300026: hht 只是将时域图分割成数个时域图, 不能画的原因 06/30 13:14
2F：→ s4300026: 应该原本就不能画(code 有错), 跟hht 无关 06/30 13:15