信号瞬时频率估计的研究

《信号瞬时频率估计的研究》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、信号瞬时频率估计方法的研究：在信号处理中，信号本身有很多重要的属性，频率特性有：带宽、各频率分量的相对幅值、频率分量间的相对相位关系等；时域特性有信号时宽等。在很多时候，对信号的处理都涉及到需要对平稳或者非平稳信号的频率特性进行估计。平稳信号的频率特性是时不变的，而非平稳信号的频率特性往往是时变的，因此，瞬时频率的定义主要是针对非平稳信号而提出的。Ville给出了一种统一的瞬时频率的定义：其中，z(t)是实信号的解析信号。瞬时频率估计的方法可以分为时频分析和时域分析两类。就平稳信号而言，由于其功率谱

2、密度函数是不随时间变化的，因此可以直接用参数化或者非参数化谱估计的方法来得到其功率谱，将功率谱中峰值所对应的频率值作为组成该平稳信号的各频率分量的频率的估计值。但是，对于非平稳信号而言，由于其功率谱密度函数是时变的，因此如果要在频域估计其瞬时频率，最简单的方法就是先将其视为短时平稳的信号，每次都用足够短的时间内的数据来构建其功率谱密度函数，将估计得到的结果作为该短时间内的信号瞬时频率，这也就是时频分析中的短时傅立叶变换方法。当然，时频分析还有诸如小波变换等其他的性能更好的变换方法这里不再展开叙述。下

3、图是用短时傅立叶变换得到的一个非线性调频信号的时频分布图：时域处理方法则主要是根据信号瞬时频率的定义，先将实信号变换为复信号，再通过对复信号的相位进行求导（模拟）或者差分（数字）的方法来求得瞬时频率。时频分析处理的好处是对于有多个频率分量的信号可以根据功率谱密度函数的各个峰值点估计出对应分量的瞬时频率。而基于相位求导或者差分的时域处理方法却是无法对多频率分量的信号进行瞬时频率估计的。针对这一问题，HUANG.N.E提出了局域波分解方法，首先将复杂的信号分解成有限个基本模式分量，再对这些基本模式进行相

4、位求导或者差分以估计各分量的瞬时频率。通过局域波分解的方法可以很好的解决相位求导或差分方法的缺陷。时域处理的好处是计算量远小于时频分析处理。这里主要讨论时域的处理。而要进行时域处理，则通常要首先将物理上的实信号变换为复信号以便取其相位。现有的两种的方法分别是正交变换和hilbert变换。正交变换可以用模拟信号处理的方法实现，也可以用数字信号处理的方法实现。正交变换的实质就是用两个正交的本振进行同相和正交两路下变频混频，由这两路下变频的输出组成实信号包络的零中频复基带信号的实部和虚部。模拟信号处理由于

5、模拟电路的固有缺点（如本振频率不稳定、本振间不正交和IQ两路通道不均衡等）将会导致输出的两路正交信号会存在相位正交性误差和幅度一致性误差。而数字信号处理的方法则可以避免这些问题。因此，现在通常用数字正交处理的方法来实现正交变换。当然，数字正交处理由于滤波器处理速度的原因，其前置预滤波器和后置低通滤波器通常需要诸如多相滤波等特别设计以满足要求。而hilbert变换则是直接在载频上将实信号变换为复信号。hilbert变换的得到的复信号也叫解析信号。Hilbert变换可以用时域的方法实现，也可以用频域的方

6、法实现。在频域方法中：先对输入的一段数据进行FFT变换，再进行如下运算：之后再运行IFFT运算即可得到实信号的解析信号。这种在频域实现hilbert变换的方法，其计算量相比时域冲激响应的方法要大很多的。因此只适用于短数据的情况。时域的方法在理论上，只要将一个实信号通过一个冲激响应为：的hilbert滤波器就可以求出实信号的正交分量。但是在实际中，hilbert变换通常是在数字域实现的。因此，需要设计一个数字滤波器，使得其输出为数字信号的正交分量。这样的数字滤波器的冲激响应是：这是一个非因果IIR滤波

7、器。在实现的时候，考虑在实时要求不高的情况下，数据输出可以有一定的延迟，且的绝对值随n增大而减小，因此可以在足够小的时候将其截断。从而构成一个非因果FIR滤波器。但是，对于很短的一段非平稳信号数据而言，hilbert变换会存在的端点效应，尤其在用频域的方法实现hilbert变换时，这种端点效应会变得非常严重。即在hilbert变换得到的数据段的两端会存在I、Q不均衡现象。针对这一问题，有人提出通过对原数据序列进行数据延拓的方法，以期在有用的数据段内避免端点效应的出现。至于如何对数据进行有效地延拓，人

8、们也提出了很多的方法：HUANG.N.E在提出局域波分解法时也提出了一种有效的数据延拓方法，并在美国申请了专利；另外也有人用神经网络的方法来进行数据延拓。实际上，基带复信号只是解析信号的频移形式而已，两者反映的基本信息是相同的，但是在工程上多用的是正交变换的方法来建立复信号。在将实信号变换为复信号之后，就可以提取其瞬时相位进行瞬时频率估计了。其处理的流程如下：正交变换Q(n)Q(n)f(n-1)I(n)S(n)其中，在进行反正切运算求复信号相位的时候，由于反正切运算的