24小时热门版块排行榜

返回列表

xiao又又

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 14
帖子: 4
在线: 2.2小时
虫号: 4394812
注册: 2016-02-04
专业: 河流海岸动力学与泥沙研究

[求助] 请教一下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;

回复此楼

» 猜你喜欢

1楼 2016-02-15 21:30:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 xiao又又的主题更新

返回列表