声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2740|回复: 7

[HHT] 怎么让hht变换的频率划分的更细啊

[复制链接]
发表于 2007-10-15 14:10 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我用riling的程序作hht变换,可是发现他的程序只是把频率归一化为400份,这跟小波包分析还有什么区别啊。
有没有直接给出原始频率信号的源程序啊?
回复
分享到:

使用道具 举报

发表于 2007-10-15 14:18 | 显示全部楼层

回复 #1 1021830113 的帖子

HHT从原理上和小波有本质的不同,分解出来的各阶和小波也不一样。楼主想得到什么样的源程序?
发表于 2007-10-15 14:23 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-7 16:08 编辑
原帖由 1021830113 于 2007-10-15 14:10 发表
我用riling的程序作hht变换,可是发现他的程序只是把频率归一化为400份

这句话从何说起?我之前没有留意过,能不能说得更详细一些:loveliness:
 楼主| 发表于 2007-10-15 14:25 | 显示全部楼层
我发现,hht变换最后把采样频率的一半进行400等分,而信号在每个时间点上对应的频率值都近视成划分好的频率值,而不是给出原始的信号频率值。我想得到一个程序,计算出来的结果中频率没有经过近似处理。我考虑这可能是riling为了便于处理数据而采取的手段,但是完全可以使用稀疏矩阵表征hht变换结果啊。
发表于 2007-10-15 14:36 | 显示全部楼层
 楼主| 发表于 2007-10-15 14:42 | 显示全部楼层
刚刚发现A和f应该是每个imf变换后的精确幅值和频率值。程序画图的时候才把频率近似了一下。
不知道这点发现对不对。
 楼主| 发表于 2007-10-15 14:48 | 显示全部楼层
%TOIMAGE  transforms a spectrum made of 1D functions in an 2D image
%
% [im,tt,ff] = TOIMAGE(A,f,t,splx,sply) transforms a spectrum made
% of 1D functions (e.g., output of "hhspectrum") in an 2D image
%
% inputs :   - A    : amplitudes of modes (1 mode per row of A)
%            - f    : instantaneous frequencies
%            - t    : time instants
%            - splx : number of columns of the output im (time resolution).
%                     If different from length(t), works only for uniform
%                     sampling.
%            - sply : number of rows of the output im (frequency resolution).
% outputs :  - im   : 2D image of the spectrum
%            - tt   : time instants in the image
%            - ff   : centers of the frequency bins
%
% Examples : [im,tt,ff] = toimage(A,f);[im,tt] = toimage(A,f,t);[im,tt,ff] = toimage(A,f,sply);
%            [im,tt,ff] = toimage(A,f,splx,sply);[im,tt,ff] = toimage(A,f,t,splx,sply);
%
%
% See also
%  emd, hhspectrum, disp_hhs
%
% G. Rilling, last modification 3.2007
% gabriel.rilling@ens-lyon.fr
function [im,tt,ff] = toimage(A,f,varargin)

DEFSPL = 400;
error(nargchk(2,5,nargin));

应该在TOIMAGE函数中修改红色部分,在这里可以选择画图时的频率划分程度。
发表于 2007-10-15 19:11 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-7 16:08 编辑
原帖由 1021830113 于 2007-10-15 14:48 发表
%TOIMAGE  transforms a spectrum made of 1D functions in an 2D image
%
%  = TOIMAGE(A,f,t,splx,sply) transforms a spectrum made
% of 1D functions (e.g., output of "hhspectrum") in an 2D image
...

这只能看他的文章有没有涉及这方面了,如果没有,可能是随便取的一个经验值。其实在其他函数也有这个情况,例如:

  1. T = 1:4:N;
  2. imagesc(flipud(rs1(1:128,:)))
复制代码

这些都是跟实际信号相关的值,因此,当你把它们弄懂以后,就可以自己修改了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 17:12 , Processed in 0.051707 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表