马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 baishp 于 2011-9-17 23:16 编辑
%Data是工程实测的一列随机信号
imf=emd(Data);
运行后得到一个15*58320的矩阵imf。
%Data的IMF分量图
for k=1:15
subplot(5,3,k)
plot(imf(k,:))
end
求imf各分量振幅与瞬时频率
[A_imf,f_imf]=hhspectrum(imf);
%A_imf各分量振幅向量图
for k=1:15
subplot(5,3,k)
plot(A_imf(k,:))
end
%f_imf各分量瞬时频率向量图
for k=1:15
subplot(5,3,k)
plot(f_imf(k,:))
end
可以看出,在低频imf分量中,瞬时频率有多处出现频率的最大值0.5HZ。这是为什么?这个肯定与实际不符啊。
单独将第14个分量的瞬时频率作图:
plot(f_imf(14,:))
频率轴用对数坐标再作图:
plot(log10(f_imf(14,:)))
这下看的更清楚了吧?
请教各位能人,这是为什么?怎么解决? |