| 查看: 1264 | 回复: 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; |
回复此楼» 猜你喜欢
科研孤儿太难了
已经有9人回复
-
时间戳又变了
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
40