线调频小波分析.docVIP

线调频小波分析 一种自适应chirplet分解的快速算法的若干硏究 针对信号自适应chirplet分解未知参数多、实现起来比较困难的特点， 文献 ⑴提出了一种新的chirplet分解快速算法。该算法利用计算信号的二 次相位函数，得到其能量分布集中于信号的调频率曲线上的结论，此时通 过谱峰检测可同时获得chirplet调频率、时间中心和幅度的估计；然后通 过解线性调频技术获得其初始频率和宽度的估计，仿真结果验证了本文算 法的有效性。 如果线性调频信号具冇光滑的高斯包络，那么它就变成如下典型的调 副调频信号，称为chirplet函数: gn(t)?(??)2?0.25 n(t?tn)2exp{??jwn(t?tn)?j?(t?tn)2} (1) 22?n 其中?nmwn,?n分别表示chirplet函数的宽度，时间中心,初始频率 和调频率。此外，chirplet函数的Wigner-Ville分布(WVD)具冇如下形式: (t?tn)2 2WVDgn(t,w)?exp{??4?n[(w?wn)/2??n(t?tn)]2} (2) 2?2?nl 由上式可见，chirplet函数是唯一的WVD为肺腑能量的函数，因此它 在联合时频分析屮扮演着重耍的角色。 常用的联合吋频分析方法是将信号分解为一系列其函数的线性组合, 通过抑制函数的视频特性来了街特分析信号，其主要缺点是运算量大，尤 其是参数初值的估计，需要在一定范围内进行搜索，或者自适应迭代等, 不利于具体的应用；本文提岀一种新基于chirplet分解的时频分析快速算 法，该算法利用计算信号的二次相位函数，得到其能量分布集中于信号的 调频率曲线上的结论，通过谱峰检测可同时获得chirplet调频率、时间中 心和幅度的估计；然后通过解线性调频技 术化的起初始频率和宽度的估计。算法实现简单，即算量小，能够保 留信号更多的视频特性。 自适应chirplet分解法简介 自适应chirplet分解时以基函数与特分析信号最相似的原则來选择基 的，将特分析信号s(t)表示为一组线性调频小波基卯(t)的线性叠加： s(t)??Cngn(t) (3) n?0? 英中即(t)加式(1)所示，基函数按照下列准则逐一个自适应估计： Cn2?max?sn(t),gn(t)? (4) gn2 其中 sn(t)?sn?l(t)?Cn?lgn?l(t) (5) sO(t)?s(t) sn(t)是sn?l(t)向基函数gn?l(t)作正交投影后的剩余量，可以表示为 sn(t)?s(t)? M为基函数的个数。i?l,i?n?C￡(t) (6) iiM 采用(4) ~ (6)的方法，可完成对信号s(t)的分解，这个过程可描述 为：(a)在第一次分解过程中，按照式(4)估计与sO(t)最匹配的基函 数gO(t),利用式(5)得到剩余量sl(t)o 在第二次分解过程中，按照(4)、5)两式估计与sl(t)最匹配的 基函数gl(t)并得到剩余量s2(t)o 在以后每一步分解过程屮，重复(a), (b)两步，直到剩余能量 满足事先 给定的没摩一个条件，因而，设计最佳的基函数是分解的关键。 自适应chirplet分解快速算法 由式(4)可见，系数5的计算是一个多维非线性优化问题，通常没 有解析形式的解，因此使得己有算法的计算量非常大，不利于具体应用。 文献⑴提出一种新的快分解算法，过程如下； 考虑单分量信号 S(t)?C0(??)2?0.25 n(t?tn)2exp{??jw0(t?t0)?j?0(t?t0)2} (7) 22?n 首先定义其二次相位函数(Quadratic Phase, OP)： ? ■ OP(t,u)??s(t??)s(t??)ejw?d? (8) 02 将式(7)代入(8)可以得到 OP(t,u)?C(??)2 02?0.250(t?t0)2exp{??22?0 ?2jw0(t?t0)?2j?0(t?t0)2}*?exp{? 0??j?2(2?0?u)}d?2?02 (9) 利用积分公式 ? ? ? ? ???e?jmt2dt??me?j?4m?0 (10) 2?????2j?exp??j?(2??u)d??e (11) ??02???10??0?2(2?0?u)2?4?0 英中 11? (12) ??2?02?0(u?2?)4 由上式可以看出，当2?0?u?0时，有 OP(tzu)?C(??) ?202?0.2502(t?t0)2??2exp{?*?exp{?2}d??22?0?00(t?t0)12C0exp{?222?0 (13) 由上式可以得，而此相位函数OP(t,u)的峰值位于(02?0)处，而且峰值 大小为12C0,因此，参数{t0,2?0,C0}的估计方法可理解为：首先计算信号 的二次相位2 ?函数，然