最近在搞hilbert滤波器,有少许收获,写下来共享。
希尔伯特(Hilbert)变换可以提供90°的相位变化而不改变频谱分量的幅度,即对信号进行希尔
设离散时间信号x(n)的希尔伯特变换是x^(n),希尔伯特变换器的单位抽样响应为h(n),由于其频率响应为:
单位抽样响应为:
离散时间信号的傅里叶变换为:
解析信号:
性质:
(1)希尔伯特变换器是幅频特性为1的全通滤波器,信号通过希尔伯特变换器后,其负频率成分作+90°相移,而正频率成分作-90°相移。信号通过希尔伯特变换器后
(2)对z(n)做傅里叶变换,由于x(n)的DFT为X(k),k=0,1,……,N-1,k=N/2,……,N-1对应负频率,则
Z(k)=
实现方法:
(1)希尔伯特变换器是幅频特性为1的全通滤波器.
希尔伯特变换中的h(n) 可以表示为:
则离散信号x(n)的离散希尔伯特变换可以表示为:
由此可见,离散希尔伯特变换器可以由FIR(有限冲击响应)滤波器来实现,可以用窗口法来设计FIR滤波器实现希尔伯特变换。利用矩形窗设计的55阶FIR滤波器幅频响应如图1所示
但要注意的是,利用FIR滤波器实现希尔伯特变换将会使输出信号延迟N/2(N为滤波器系数长度),而且输出信号的前N个数据和最后N个数据也是不对的,因为此时输入数据已经为0。
理想的Hilbert滤波器是非因果不可实现的,因为h(n)在n<0时不为0。
所以在求它的相频响应的时候,加了一个线性相位项在里面。假如你的滤波器系数h(n)的范围为-N<=n<=N,则相当于按照0<=n<=2N的因果可实现的滤波器得到的相频响应,所以在原本[0,pi]上为-90度。在[pi,2*pi]上为90度的相频响应的基础上叠加了一个exp(-jwN)的线性相位项。
(2)根据希尔伯特变换与解析信号的关系,可通过傅里叶变换实现。
应用领域:
1.希尔波特变换滤波器的一个主要作用就是构造解析信号。解析信号即没有负频率成分的复信号,将一个实信号加上它的希尔波特变换作为它的虚部即构成解析信号。通信中的单边带通信系统就是利用解析信号实现的,它可以节约一半的频带。
2.