| 查看: 1174 | 回复: 0 | ||
[求助]
请教一下emd程序问题出在哪里，得到的许多都是负频
|
|
我的时间序列n=432，想运用emd分解多个IMF，并求出振幅与频率，但出现很多负频，新手一枚，请各位指点迷津，小女不胜感激。 z= Book1S2; t=(1:1:432); figure plot(t, z) title('原始信号') % EMD imf=emd(z); emd_visu(z,t,imf) [A,f,tt]=hhspectrum(imf); [im,tt]=toimage(A,f，tt,length(tt)); disp_hhs(im); xlabel('时间（月）') ylabel('含沙量'); figure disp_hhs(A,tt) 未定义与 'double' 类型的输入参数相对应的函数 'instfreq'。 出错 hhspectrum (line 79) f(i,  =instfreq(an(i,',tt,l)';z=imf1; >> hx=hilbert(z); xr=real(hx);xi=imag(hx); %计算瞬时振幅 sz=sqrt(xr.^2+xi.^2); %计算瞬时相位 sx=atan(xi./xr);sx=angle(hx); %计算瞬时频率 dt=diff(t); dx=diff(sx); sp=dx./dt; |
回复此楼» 猜你喜欢
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有6人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有6人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
遇见不省心的家人很难过
已经有22人回复