| 查看: 1182 | 回复: 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; |
» 猜你喜欢
290求调剂
已经有9人回复
0856化工专硕求调剂
已经有6人回复
0856材料专业298分有科研经历 硕士研究生调剂自荐信
已经有6人回复
材料化工调剂
已经有7人回复
化工299分求调剂 一志愿985落榜
已经有3人回复
306分材料调剂
已经有4人回复
302材料工程求调剂
已经有5人回复
272求调剂
已经有5人回复
材料学调剂
已经有9人回复
321求调剂一志愿东北林业大学材料与化工英二数二
已经有6人回复













=instfreq(an(i,
回复此楼