HHT10-12.docVIP

10总体经验模态分解（EEMD）、Fourier变换、HHT EEMD实际就是噪声分析法和EMD方法的结合，抑制模态混叠。Fourier变换是将任何信号分解为正弦信号的加权和,而每一个正弦信号对应着一个固定的频率(Fourier频率)和固定的幅值,因此,用Fourier 变换分析频率不随时间变化的平稳信号是十分有效的。但对于频率随时间变化的非平稳信号,Fourier 变换就无能为力了。HHT是历史上首次对Fourier变换的基本信号和频率定义作的创造性的改进。他们不再认为组成信号的基本信号是正弦信号,而是一种称为固有模态函数的信号,也就是满足以下两个条件的信号: (1) 整个信号中,零点数与极点数相等或至多相差1 ; (2) 信号上任意一点,由局部极大值点确定的包络线和由局部极小值点确定的包络线的均值均为零,即信号关于时间轴局部对称。无论Hilbert谱中的频率还是边际谱中的频率(即瞬时频率) ,其意义都与Fourier分析中的频率(即Fourier 频率) 完全不同,但在Fourier分析中,某一频率处能量的存在,代表一个正弦或余弦波在整个时间轴上的存在,而边际谱h中某一频率处能量的存在仅代表在整个时间轴上可能有这样一个频率的振动波在局部出现过,h越大,代表该频率出现的可能性越大。HHT时频灰度谱转黑白谱后面加上这个就可以了colormap(flipud(gray))小弟刚刚接触HHT，也不是学信号的，只是用HHT这个工具处理信号，在处理过程中遇到了这样的问题：对实测信号直接EMD，然后作HHT谱图如下： 2009-4-30 15:41 上传 下载附件 (45.25 KB) 然而对于实测信号的分析首先是要进行去噪处理的，我就试着去掉了两个高频IMF，然后作HHT谱图如下： 2009-4-30 15:41 上传 下载附件 (24.36 KB) 为什么在去噪之后的HHT谱中高频部分出现了很强烈的振幅，而去噪之前是没有的？请帮帮忙指点小弟该怎么做，谢谢大家！测得信号直接求边际谱： 2009-5-4 13:53 上传 下载附件 (16.53 KB) 去掉两个高频IMF后的边际谱： 2009-5-4 13:53 上传 下载附件 (20.44 KB) 为什么会这样呢？去掉高频IMF反倒在高频的地方出现了幅值！同样的程序，难道是我去高频的时候出错了？我去高频是直接相加的：IMF1=imf(3,:)+imf(4,:)+imf(5,:)+imf(6,:)+imf(7,:)+imf(8,:)+imf(9,:)+imf(10,:);? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?[A,fa,tt]=hhspectrum(IMF1);还是因为非正交性而不能直接相加？搞不懂？楼主你的程序“IMF1=imf(3,:)+imf(4,:)+imf(5,:)+imf(6,:)+imf(7,:)+imf(8,:)+imf(9,:)+imf(10,:);”不对。用这个语句后IMF1只有一行了。这样你再进行边际谱分析肯定不对。建议修改如下：for i=3:10IMF1(i-2,:)=imf(i,:);end谢谢了，这样是对的，太感谢了。