分数阶傅立叶变换的最优阶数论文

《分数阶傅立叶变换的最优阶数论文》由会员分享，可在线阅读，更多相关《分数阶傅立叶变换的最优阶数论文（18页珍藏版）》请在金锄头文库上搜索。

1、现代数字信号处理学 号：140808040219学生所在学院：测试与光电工程学院学 生 姓 名 ：任 课 教 师 ：李 志 农教师所在学院：测试与光电工程学院2015年 1月分数阶傅立叶变换的最优阶数摘 要：传统傅立叶变换描述了信号时域或频域的特性，而不是描述信号时频特性，于是人们提出了一系列新的时频分析理论和方法来处理非平稳信号，分数阶傅立叶变换为其中的一种。本文主要介绍了它的定义、性质，还有它的离散算法，介绍了求最优阶数的方法，主要是峰值搜索算法。最后进行仿真验证，利用 MATLAB 对一个已知的 chirp 信号求解最优阶数。关键字：傅立叶变换；分数阶傅立叶变换；峰值搜索算法；MATLA

2、B；最优阶数1 引言传统的傅立叶变换在所有的信号处理工具中是应用最广泛、研究最成熟的数学工具，作为一种线性算子，传统傅立叶变换可视为在时频平面上，信号从时间轴逆时针旋转 到频率轴，而 FRFT 作为 FT 的广义形式可理解为对信2号旋转任意角度的线性算子，从而可以得到信号的任意阶次或者任意分数阶傅立叶域上的 FRFT 表示，并且在保留了传统的 FT 所有性质和优点的基础之上又增添了新的优势。2 FRFT 的定义及其性质1.1 FRFT 的定义如图 1.所示如果把信号的分数阶傅立叶变换看作是从时间-频率平面旋转的话，那么傅立叶变换就相当于在时频平面中逆时针旋转了 角度，从时间域2变换到频率域。令

3、 ， 是一个分数，那么就可以在时频平面内以任意2p角度的旋转定义线性算子 ，记作 ，我们就可以把傅立叶变换推广到2pRpF任意角度即分数阶傅立叶变换。图 1. 平面旋转 角度变成 平面(,)tw(,)uv分数阶傅立叶变换的定义为：(2-1)()()(,)pp pXuFxKutxd其中是核函数 ，,pKt(2-2)这里 ，p=1 或-1 时退化成为常规的傅立叶变换和逆变换。,02p分数阶傅立叶变换和经典傅立叶变换具有以下的关系： 1.分数阶傅立叶变换是线性算子 2.周期性(2-3)04()()Fxttx(2-4)12X1.2 分数傅立叶变换的性质（1） 线性分数阶傅立叶变换为线性变换，满足叠加原

4、理：若 和 分别是原函()PFft()Pgt数 和 的 阶分数傅立叶变换，则有()ftgp(2-5)1212()()()()PPPFcftgtcftt（2） 旋转相加行(2-6)PqF（3） 逆(2-7)1()P（4） 酋性(2-8)1()()PHF（5） 交换性(2-9)1221PP（6） 结合性(2-10)331212()()PPFF（7） 周期性(2-11)44()()()pkppk pfufxfxfu （8） 特征函数(2-12)e(/2)pl lFjl（9） 卷积、相乘、相关 函数 、 在阶次 p 分数阶傅立叶域的卷积记作 ()ftg(2-13)*()pppfgft函数 、 在阶次

5、p 分数阶傅立叶域的乘积记作()ftt(2-14)()pppFfft函数 、 在阶次 p 分数阶傅立叶域的相关定义为()ftgt(2-15)*(),(),pCORftgCPNVftg（10）时移特性(2-16)2sincosin()(c)jjuppFxtXue（11）频移特性(2-17)21(sincosin)()(si)jvupjt peXv（12）尺度特性(2-18)22cost(1)2t in() ()cosujp pjFxcteu其中 2arctn/（13）微分特性(2-19)()()cos()sinpppFxtXujXu（14）积分特性(2-20)22tantan()sec()uzl

6、 j jpxded3 分数阶傅立叶变换的离散的离散算法分数阶傅立叶变换的出现引起了各个领域研究人员的重视，在工程上也有十分广阔的应用前景。在数字信号处理的应用中，必须采用离散形式的分数阶傅立叶变换（DFRFT） ，这使得离散分数阶傅立叶变换及其快速算法的研究显得十分重要。目前，DFRFT 的离散化算法主要有四种：1. 利用 来计算离散 FRFT 的核矩阵，从而利用 FFT 来320()iiFW计算 DFRFT其中 W 是离散傅立叶变换核矩阵。这种方法实际上缺乏理论基础，而且其离散 FRFT 矩阵不满足连续 FRFT 的旋转相加性，因此不能用相同的方法计算逆 FRFT。实际计算所产生的误差比较大

7、，与连续 FRFT 没有相似的输出结果。其计算复杂度与传统傅立叶变换相同，为 。2(log)ON2. 分解方法根据 FRFT 的表达式，将 FRFT 分解为信号的卷积形式，从而利用 FFT 来计算 FRFT。这种方法思想比较直观，计算出的记过与连续 FRFT 的输出比较接近。但它要经过一次 2 倍内插和 2 倍抽取，而且还要进行坐标的无量纲化，实现起来较为烦琐。其计算复杂度为 。2(log)ON3. 利用矩阵的特征值和特征向量来计算 DFRFT这种方法保持了连续 FRFT 的特征值-特征函数的关系，克服了第一种方法中特征值和特征向量不匹配的缺点。采用了两种正交映射的办法 DFT 的 Hermi

8、te 特征向量，由此开发出两种快速方法，即 OPA 方法和 GSA 方法，这两种方法都有和连续 FRFT 相近的输出结果，可逆性好。其计算复杂度均为。2()ON4. 直接对 FRFT 进行离散化来计算 DFRFT。 这种方法采用直接将连续 FRFT 离散化的方法来获得离散 FRFT 的核矩阵，避开了烦琐的特征值和特征向量匹配问题以及矩阵的正交归一化运算，与连续 FRFT 有相似的输出结果，该算法的计算复杂度为 。2()ON在目前的研究中，采用的最多的是分解方法和矩阵特征值和特征向量两种方法，下面的内容重点介绍这两种方法。3.1 分解方法所谓分解方法是指根据 FRFT 的表达式，将 FRFT 分

9、解为信号的卷积形式，从而利用 FFT 来计算 FRFT。这种算法由 H.M.Ozaktas 等人提出，其计算速度几乎与 FFT 相当，被公认为目前计算速度最快的一种 FRFT 数字计算方法，非常适合于对信号进行实时的 FRFT 计算。但这种快速算法的运算机理决定了在进行 FRFT 数值计算之前必须对原始信号进行量纲归一化处理。3.1.1 量纲归一化原理如果一个信号在时间轴或频率轴的一个子集取非零值，并且取非零值的条件限定在有限区间内，则称该信号在时间轴或频率轴上是紧凑的。从理论上讲，任何信号和它的傅立叶变换不可能是同时紧凑的。然而我们实际中需要处理的信号往往是时限的和带限的。信号的时宽带宽积可

10、以用来确定信号的采样频率和采样点数，用于唯一地从离散化的信号中恢复原始信号。假定原始连续信号在时间轴和频率轴上都是紧凑的，其时域表示限定在区间，而其频域表示限定在区间 ， 和 分别表示/2,t/2,ftf信号的时宽和带宽。信号的时宽带宽积为 ，根据不确定性定理， 的NtN值应当大于 1。由于时域和频域具有不同的量纲，为了 FRFT 计算处理方便，需要将时域和频域分别转换成无量纲的域。具体方法是引入一个具有时间量纲的尺度因子 ，并定义新的尺度化坐标S,/xtS/vf新的坐标系 实现了无量纲化。信号在新坐标系中被限定在区间()x和 内 。 为 使 两 个 区 间 的 长 度 相 /2,/tt/2,

11、/fSf等 ， 选 择 ，则两个区间长度都等于无量纲量 ，即两t xtf个区间归一化为 。归一化以后信号的 Wigner-Ville 分布限定在以/,x原点为中心，直径 的圆内，如图 3-2 所示。最后根据采样定理对归一化后的信号进行采样，采样间隔为 ，采样点数为 .1/x2Nx图 3.1 归一化后的时频支撑区域3.1.2 分解算法可以把（2-1 ）式改写为(3-1)221cotcotcot()exp()expexpcs2PjuXujtjjtudt 式可以具体分解为以下几个步骤：（1） 用 chirp 信号调制信号 ：()xt(3-2)2eptan(/)gj（2） 调制信号与另一个 chirp

12、 信号卷积：(3-3)2()x()cs()uAjutgtd（3） 用 chirp 信号调制卷积后的信号：(3-4)2()tan(/)PXgj要实现 FRFT 的数值计算必须对以上每个分解步骤都进行离散化处理，具体的实现过程如下：（1）对 与 chirp 信号的乘积进行采样 ()xt假定分数阶次 ，chirp 信号的调频率 ， 经 chirp 1,ptan(/2)1()xt信号调制后所得的信号时宽带宽积可以是原信号的时宽带宽积的两倍，所以要求 的采样间隔为 ，如果原来的采样间隔是 ，可通过插值的方()xt/2x/法获得样本值，然后再 chirp 信号的离散采样相乘，以得到 的采样值。()gt（2

13、）实现 与 chirp 信号的卷积()gt由于 是带限信号，所以 chirp 信号也可取其带限形式。所以有：(3-5)2()exp()cs()()uAjutgtdAhutgd 其中 2(,xhHvdju而 则是 chirp 信号傅立叶变换()v(3-6)231()exp()4cscsjjvv于是，式的离散形式为(3-7)22NmngAhgxx这一离散卷积可利用 FFT 来计算。（3）计算分数阶傅立叶变换 的采样值 ()pXu(/)pXmx由于在第一步操作时对信号作了 2 倍内插操作，所以在最后的结果需要对再进行 2 倍抽取，以得到离散采样 。(/)pXmx(/)p总之上述方法从连续信号 的 N

14、 个离散采样 开始，最后得到由()xtx的 N 个离散样本 值得注意的是上述方法只适用于 的()pu/pXm1p情况，对位于该区间外的情况，可以利用分数阶傅立叶变换的周期性将阶次变到 后再进行计算。13.1.3 特征值和特征向量方法分解方法虽然计算复杂度较低，但不严格遵守分数阶傅立叶变换的旋转特性，因此不能从变换后的信号通过其逆变换精确地恢复原始信号。为了克服这种方法的缺点和不足，Pei Soo-Chang 等人提出一种新的离散化方法，该方法具有与连续情况相似的变换性质和结果，并可以通过其逆变换恢复原始信号。 这种方法就是特征值和特征向量方法，它从连续傅立叶变换的特征函数为Hermite 函数

15、出发，通过对 Hermite 函数的离散化近似和正交投影，得到一组与 Hermite 函数形状相似的 DFT 矩阵的正交化离散 Hermite 特征向量，然后，按照连续分数阶傅立叶变换的核函数谱分解表达式，构造离散分数阶傅立叶变换矩阵。3.1.3.1 DFT 矩阵特征值和特征向量傅立叶变换的特征函数 为 Hermite-Gaussian 函数，其表达式为：()nu(3-8)2()unAHe其中 ， 为 阶 Hermite 多项式。1/42!nA0,; )nDFT 矩阵的特征值及重复度如下表：DFT 矩阵 F 的特征值是 ，共有 1、-j、-1、j 四个值，每个特征值(/2)jke对应的特征向量全体组成一个特征子空间，记为 ，每个特征值的重0123,E复度决定了子空间的秩。矩阵 S 可用于计算 DFT 矩阵 F 的特征向量，S 的表达式为：(3-9)210cos00210cos(1)SN 可以证明矩阵 S 和 F 满足乘法交换性，即 SF = FS。因此矩阵 S 的特征向量也是矩阵 F 的特征向量，但它们对应不同的特征值。3.1.3.2 DFT