我不太理解你的”取近似相当于高数上把图形离散为许多小矩形“的意思,是在Hilbert-Huang Spectrum上取吗?E是一个二维矩阵吧?函数[E,tt1]=toimage(A,fa,tt,length(tt))的意思是要根据瞬时幅度A,瞬时频率fa,时间tt转化成一个二维的图像吗?因为是对时间的积分,要求的是同一频率下H(w,t)对整个时间跨度的积分,所以要先对所有的瞬时频率排序,并将对应的瞬时幅度也排序,那E是不是满足这种条件了呀?然后用disp_hhs(E,tt1)就可以得到Hilbert-Huang Spectrum了,就是有颜色一张图片,看得眼花缭乱。。。
如果我不调用toimage函数,直接用以下的程序求Hilbert-Huang Spectrum可以吗?
imf=emd(y);%y是原始数据,列向量,9000个元素
[A,fa,tt]=hhspectrum(imf);%得到的A和fa是12*8998的矩阵,不知道为什么A和fa的长度减小了2?而tt是一个1*8998的行向量,范围是2~8999。。
figure;
for k=1:length(imf) %幅度-时间谱
plot(tt,A(k,:));
hold on;
end
figure;
for k=1:length(imf) %频率-时间谱
plot(tt,fa(k,:));
hold on;
end
我就是想把所有的瞬时幅度和瞬时频率叠加在一张图上,我不知道这样理解对不对,请不吝赐教!谢谢!
|