现代信号处理_08

1、1内内 容容2内内 容容3研究的必要性研究的必要性高阶统计量高阶统计量高阶谱高阶谱高阶累积量和多谱的性质高阶累积量和多谱的性质三阶相关和双谱及其性质三阶相关和双谱及其性质基于高阶谱的相位谱估计基于高阶谱的相位谱估计基于高阶谱的模型参数估计基于高阶谱的模型参数估计多谱的应用多谱的应用 参考:现代数字信号处理(184-199；204-205)4 研究高阶谱的必要性研究高阶谱的必要性v 关于关于模型参数估计问题模型参数估计问题所谓模型参数估计，就是根据有限长的数据序列所谓模型参数估计，就是根据有限长的数据序列(如模如模型输出端所观测到的信号型输出端所观测到的信号y(n)来估计图中随机信号模型来估计图

2、中随机信号模型的参数，的参数，)与前面所述不同之处在于：这里考虑了观测过程所引与前面所述不同之处在于：这里考虑了观测过程所引入的噪声入的噪声v(n). H ( z ) ( h ( n ) )v(n)y(n)x(n)u(n)5 研究高阶谱的必要性研究高阶谱的必要性v 基于二阶统计量的模型参数估计方法的缺陷基于二阶统计量的模型参数估计方法的缺陷 前述模型参数估计方法中，估计得到的模型参数仅与前述模型参数估计方法中，估计得到的模型参数仅与 信号的自相关函数或功率谱包络相匹配；其功率谱不信号的自相关函数或功率谱包络相匹配；其功率谱不 含信号的相位特性，亦称盲相。即含信号的相位特性，亦称盲相。即22)(

3、)(jueHS 这种模型只适合于高斯随机信号，因为高斯信号仅用这种模型只适合于高斯随机信号，因为高斯信号仅用 二阶统计量二阶统计量(均值和方差均值和方差)就能加以描述。就能加以描述。6 研究高阶谱的必要性研究高阶谱的必要性v二阶统计量方法的基本限制二阶统计量方法的基本限制 前面讨论的方法中，一般都假设：前面讨论的方法中，一般都假设： 信号模型中的系统信号模型中的系统H(z)是最小相位的。是最小相位的。 激励信号激励信号u(n)是均值为零，方差为是均值为零，方差为 的高斯白噪声。的高斯白噪声。 测量信号测量信号v(n)是均值为零，方差为是均值为零，方差为 的高斯白噪声；的高斯白噪声； 且且v(n

4、)与信号与信号x(n)统计无关，即统计无关，即v(n)不影响信号的谱形状不影响信号的谱形状 故有故有 2u2v)()()()()()()(22222mhmnynuEmRHSSuuyvuvxxyy7 研究高阶谱的必要性研究高阶谱的必要性v 二阶统计量方法存在的问题二阶统计量方法存在的问题 在许多实际应用在许多实际应用( (如地震勘探、水声信号处理、远程通如地震勘探、水声信号处理、远程通 信信) )中，往往不能满足上述假设；甚至系统是非线性的。中，往往不能满足上述假设；甚至系统是非线性的。 对于非高斯信号的模型参数，如仅仅考虑与自相关函数对于非高斯信号的模型参数，如仅仅考虑与自相关函数 匹配，就不

5、可能充分获取隐含在数据中的信息。匹配，就不可能充分获取隐含在数据中的信息。 若信号不仅是非高斯的，而且是非最小相位的，采用基若信号不仅是非高斯的，而且是非最小相位的，采用基 于自相关函数的估计方法所得到的模型参数，就不能反于自相关函数的估计方法所得到的模型参数，就不能反 映原信号的非最小相位特点。映原信号的非最小相位特点。 当测量噪声较大，尤其当测量噪声有色时，基于自相关当测量噪声较大，尤其当测量噪声有色时，基于自相关 函数的估计方法所得到的模型参数有较大的估计误差。函数的估计方法所得到的模型参数有较大的估计误差。8 研究高阶谱的必要性研究高阶谱的必要性v 解决问题的方法解决问题的方法 从观测

6、数据中提取相位信息从观测数据中提取相位信息 信号分析必须具有抗有色噪声干扰的能力信号分析必须具有抗有色噪声干扰的能力 因此，必须用高阶谱因此，必须用高阶谱( (高阶统计量高阶统计量) )来分析信号来分析信号9随机信号的高阶特征随机信号的高阶特征功率谱估计，功率谱估计，Wiener滤波器都是以信号的相关函数为工具滤波器都是以信号的相关函数为工具。 pijiqijieezAzBP121222211)()()(模型的多重性模型的多重性：考虑功率谱故由于,1)(111*jjee)(11)(121221221222PeePpiiqiipijiiqijii即不同即不同ARMA过程具有相同形状的功率谱。这一

7、特性过程具有相同形状的功率谱。这一特性 称为相关函数的多重性或模型的多重性。称为相关函数的多重性或模型的多重性。 相关函数的局限性相关函数的局限性10随机信号的高阶特征（续）随机信号的高阶特征（续） 两个具有零均值和相同方差的高斯白色噪声和两个具有零均值和相同方差的高斯白色噪声和 指数分布白色噪声显然是不同的随机过程，但它指数分布白色噪声显然是不同的随机过程，但它 们的功率谱相同。们的功率谱相同。 用这样两个白色噪声激励同一个用这样两个白色噪声激励同一个ARMAARMA模型，产生的模型，产生的 两个两个ARMAARMA过程显然是不同的随机过程，但它们的功过程显然是不同的随机过程，但它们的功 率

8、谱相同。率谱相同。 两个灰度图相同的图像有可能是不同的图像。两个灰度图相同的图像有可能是不同的图像。以上事实说明以上事实说明, 要准确地刻画随机信号要准确地刻画随机信号, 仅使用相关函仅使用相关函数数(二阶统计量二阶统计量)是不够的是不够的, 还必须使用更高阶的统计还必须使用更高阶的统计量。三阶和更高阶的统计量合称高阶统计量。量。三阶和更高阶的统计量合称高阶统计量。相关函数相关函数： 刻画信号的粗糙像高阶统计量高阶统计量：刻画信号的细节11v 特征函数与高阶矩特征函数与高阶矩 特征函数特征函数：随机变量：随机变量 x 的特征函数定义为的特征函数定义为)1 ()()(adxexfeEvjvxjv

9、x或或)1 ()()(bdxexfeEssxsx其中其中 f(x) 是随机变量是随机变量 x 的概率密度函数。的概率密度函数。 高阶矩高阶矩：对：对(1b)求求k 阶导数，得阶导数，得kksxkkdssdexEs)()(则随机变量则随机变量x 的的k 阶矩阶矩(即即k 阶原点矩阶原点矩)定义为定义为)2()()0(0skkkkkdssdxEm由于由于k 阶矩由阶矩由 生成，故特征函数生成，故特征函数 为随机变量为随机变量x的矩生成函数的矩生成函数(矩母函数矩母函数)，又成为第一特征函数，又成为第一特征函数。)(s)(s12v 累积量生成函数与高阶累积量累积量生成函数与高阶累积量(cumulan

10、t)累积量生成函数累积量生成函数)3()(ln)(avv或或)3()(ln)(bss称为累积量生成函数称为累积量生成函数(第二特征函数第二特征函数或或累积量母函数累积量母函数)。高阶累积量高阶累积量：随机变量：随机变量x x 的的k k 阶累积量阶累积量定义为定义为 ) 3()(0skkkdssdc即累积量生成函数的即累积量生成函数的k k 阶导数在原点的值。阶导数在原点的值。13v累积量生成函数与高阶累积量累积量生成函数与高阶累积量(cumulant)高阶矩与高阶累积量的关系高阶矩与高阶累积量的关系.)(61243)5()(23)(41412213122443131213321212211m

11、xEmmmmmmmcmxEmmmmcmxEmmcxEmc 关系关系:(注意注意：k 阶中心矩定义为阶中心矩定义为 ) 结论结论: 二、三阶累积量分别是二、三阶中心矩；均值为二、三阶累积量分别是二、三阶中心矩；均值为 零时零时, 就是二、三阶相关就是二、三阶相关(矩矩) 四阶以上的累积量不等于相应的中心矩四阶以上的累积量不等于相应的中心矩kxxE)(14v 累积量的物理意义累积量的物理意义 高斯随机变量的高阶矩与累积量高斯随机变量的高阶矩与累积量 高斯随机变量可用二阶矩完全描述。实际上高斯随机变量可用二阶矩完全描述。实际上,零均值高斯零均值高斯 随机变量的随机变量的k 阶矩阶矩(或零均值的或零均

12、值的k 阶中心矩阶中心矩)为为 高斯随机变量只有一阶和二阶累积量；其二阶以上的累高斯随机变量只有一阶和二阶累积量；其二阶以上的累 积量为零积量为零, 它不提供新的信息。即它不提供新的信息。即为奇数，为偶数kkkxEmkkk0,)1(,.,5 , 3 , 1 可见，其高阶矩仍然取决于二阶矩可见，其高阶矩仍然取决于二阶矩 。2 若若任一随机变量与高斯随机变量有相同的二阶矩任一随机变量与高斯随机变量有相同的二阶矩, 则累积则累积 量就是它们高阶矩的差。故有如下累积量的物理意义。量就是它们高阶矩的差。故有如下累积量的物理意义。) 3(0,221kccmck15v 累积量的物理意义累积量的物理意义 一一

13、阶累积量数学期望阶累积量数学期望:描述了概率分布的中心描述了概率分布的中心 二阶累积量方差：二阶累积量方差： 描述了概率分布的离散程度描述了概率分布的离散程度 三阶累积量三阶矩：三阶累积量三阶矩： 描述了概率分布的不对称程度描述了概率分布的不对称程度累积量累积量衡量衡量任意随机变量任意随机变量偏离偏离正态正态(高斯高斯)分布的分布的程度程度物理意义物理意义偏态与峰态偏态与峰态33csx 将三阶矩除以均方差的三次方将三阶矩除以均方差的三次方 ,得偏态系数或得偏态系数或偏态偏态：3 将四阶累积量除以均方差的四次方将四阶累积量除以均方差的四次方 ,得得峰态峰态:433344444422444mmmm

14、cx16高阶谱高阶谱功率谱的缺点功率谱的缺点：由功率谱只能恢复 ，不可能恢复自相关函数辨识系统，无法辨识非最小相位系统 “模型的多重性” “自相关函数等价性” “功率谱等价性”)()()()(*2fXfXfXfpx)( fX)( fX17高阶谱（续）高阶谱（续） 含义含义：高阶谱(Higher-order spectrum),又称多(polyspectrum), 是信号多个频率的能量谱。 111111)(1111),(),(krjkkxrkkxecS定义定义：高阶谱定义为：高阶谱定义为k阶累积量的阶累积量的k-1维维DFT，即，即 条件： “绝对可求和”11),(11krkkxrc通常将通常将

15、 的累积量谱称为的累积量谱称为高阶谱高阶谱或或多谱多谱。3k常用常用:常用的高阶谱是三阶谱常用的高阶谱是三阶谱(双谱双谱)和四阶谱和四阶谱(三谱三谱)。 18高阶谱（续）高阶谱（续）222111)(21321),()(rjxrxecB二阶谱二阶谱即为即为功率谱功率谱，它是单个频率的谱它是单个频率的谱。三阶谱三阶谱为为双谱双谱(bispectrum)，即两个频率的谱即两个频率的谱四阶谱四阶谱为为三谱三谱(trispectrum)，即三个频率的谱即三个频率的谱333221121)(3214321),()(rjxrrxecT19高阶谱（续）高阶谱（续）功率谱功率谱：2)()(Xpx双谱：)()()(

16、)()()(),(21*21212121XXXXXXBx三谱：)()()()()()()()(),(321*321321321321XXXXXXXXTx（1 1）双谱估计的直接方法双谱估计的直接方法：)()()(),()()(21*2121ffXfXfXffBfxnx20高阶谱（续）高阶谱（续）（2）双谱估计的间接方法双谱估计的间接方法：双谱),()(3nmcnxx2D-FT峰度峰度224)(3)(txEtxEk归一化峰度归一化峰度)()(2241txEtxEk 31k高斯信号亚高斯信号31k超高斯信号31k21高阶谱（续）高阶谱（续）归零化峰度3)()(2242txEtxEk高斯信号： 零峰

17、度亚高斯信号: 负峰度超高斯信号： 正峰度22v 主要性质主要性质 (8个性质) 最重要的性质如下最重要的性质如下: 和的累积量等于累积量之和，累积量因此得名。和的累积量等于累积量之和，累积量因此得名。 随机信号通过线性系统后的累积量等于该随机信号随机信号通过线性系统后的累积量等于该随机信号 的累积量与线性系统冲激响应累积量的卷积的累积量与线性系统冲激响应累积量的卷积信号的高阶累积量能够决定信号模型的冲激响应信号的高阶累积量能够决定信号模型的冲激响应h(n)， 即用信号模型的输出信号即用信号模型的输出信号(即观测到的信号即观测到的信号)y(n)的高的高 阶累积量就能决定阶累积量就能决定h(n)

18、。 23v 主要性质主要性质(续续) 确定性序列的多谱确定性序列的多谱: 确定性序列确定性序列h(1),h(k)的的k阶累量阶累量)7()().()(),.,(1111,knkhknhnhnhC其其 k 阶谱为阶谱为)8()()().()(),.,(11121121,kiikkhkHHHHS式中nnjenhH)()(24v用高阶累积量作为时间序列分析工具的原因用高阶累积量作为时间序列分析工具的原因 用高阶累量而不是高阶矩作为时间序列分析工具的原因：用高阶累量而不是高阶矩作为时间序列分析工具的原因： 理论上，使用高阶累积量可避免高斯有色噪声的影响，理论上，使用高阶累积量可避免高斯有色噪声的影响，

19、 高阶矩不能做到这一点。高阶矩不能做到这一点。高阶白噪声的高阶累积量是多维冲激函数高阶白噪声的高阶累积量是多维冲激函数, 其谱是多维其谱是多维 平坦的平坦的, 但高阶白噪声的高阶矩及其谱无此特性和优点；但高阶白噪声的高阶矩及其谱无此特性和优点；累积量问题的解具有唯一性累积量问题的解具有唯一性(因特征函数唯一地确定概因特征函数唯一地确定概 率密度函数率密度函数)，但矩问题不具有唯一性；，但矩问题不具有唯一性；两个统计独立的随机过程的累积量等于各随机过程累积两个统计独立的随机过程的累积量等于各随机过程累积 量之和，这一结论对高阶矩不成立。量之和，这一结论对高阶矩不成立。 25 三阶相关：三阶相关：

20、 设设x(n)为零均值的实平稳序列，其三阶相关函数为为零均值的实平稳序列，其三阶相关函数为)()()(),(2121mnxmnxnxEmmRx双谱双谱 R Rx x(m1,m2)的二维傅立叶变换就是双谱，其表达式为的二维傅立叶变换就是双谱，其表达式为2, 1)(2121,),(),(221112mmjmmxxemmRBv 性质性质 三阶相关函数的对称性三阶相关函数的对称性 双谱的对称性、周期性和共轭性双谱的对称性、周期性和共轭性v 定义定义26)()()(),(21*1121HHHBhv双谱中的相位信息双谱中的相位信息其中nnjenhH)()(这表明这表明双谱包含信号模型的相位信息 ；而功率谱

21、 不含相位信息 。 设设)(),(2121)()(),(),(21jjhheHHeBB则有则有)()()(),()()()(),(212121212121HHHBh且有)()(),(),(2121时当MnhnyBBhy)()(Sv确定性序列的双谱确定性序列的双谱 设设h(n)表示有限长确定性序列，其双谱可表示为表示有限长确定性序列，其双谱可表示为27v自相关函数丢失了信号的相位特性，而累积量可以得到自相关函数丢失了信号的相位特性，而累积量可以得到信号的相位谱。信号的相位谱。v实际应用中，基于三阶累积量的双谱和基于四阶累积量实际应用中，基于三阶累积量的双谱和基于四阶累积量的三谱已经够用。的三谱已

22、经够用。28v基本原理基本原理 与与AR功率谱估计功率谱估计(即单谱估计即单谱估计)相类似，相类似，AR过程的多谱过程的多谱 估计与已知的多谱相匹配的程度，也可用线性预测的多估计与已知的多谱相匹配的程度，也可用线性预测的多 谱来衡量，亦也可以用多谱的平坦度来衡量。说明如下：谱来衡量，亦也可以用多谱的平坦度来衡量。说明如下： 设用设用p 个值个值x(n)作线性预测，即作线性预测，即pkkknxanx1)()( 则预测误差则预测误差pkkknxanxnxne0)()( )()(其多谱为其多谱为)9(),.,()()().()(),.,(11,1112111,kxkkiikkekSAAAAS式中式中

23、nnjenhH)()(29v基本原理基本原理 ( (续续) ) 如果选择系数如果选择系数ak ，使得，使得)10()()().()(),.,(11111,11,kiiukkekAAAASukkekS,11,),.,( 式中式中 为一常量，则有为一常量，则有 uk,上式表明：上式表明：x(n)是由是由ukkneE,)(的非正态白噪声激励参数为的非正态白噪声激励参数为ak(k=1,p)的的AR过程产生的。过程产生的。结论结论：预测误差的多谱的平坦度可用作：预测误差的多谱的平坦度可用作AR过程多谱与实过程多谱与实际多谱接近程度的一种度量。际多谱接近程度的一种度量。30v 不稳定问题及其解决方法不稳定

24、问题及其解决方法 不稳定问题不稳定问题：用单谱(功率谱)和多谱估计AR模型参数时， 都存在稳定性问题。 解决办法解决办法 当用单谱估计AR模型时，只要把不稳定极点替换为其 倒数极点(反演技术)即可，这是因为 当用多谱估计AR模型时,不能作这种替换. 以双谱为例)()()()(1, 2111, 2zSzAzAzSxx)()()(),(12111211121, 3zzAzAzAzzSx)()()(),(2111211111211, 3zzAzAzAzzSx而故),(),(1211, 321, 3zzSzzSxx31v 多谱应用多谱应用: 用于信息学、海洋学、地球物理学、生物医学、用于信息学、海洋学

25、、地球物理学、生物医学、机械学和经济时间序列分析等学科领域机械学和经济时间序列分析等学科领域v对信号处理而言对信号处理而言，多谱可应用于自适应信号处理、阵列多谱可应用于自适应信号处理、阵列信号处理和多维信号处理信号处理和多维信号处理v信号处理中多谱的作用信号处理中多谱的作用 从正态信号中提取信息从正态信号中提取信息 检测和定性分析系统的非线性特征检测和定性分析系统的非线性特征 从有色正态噪声中提取信号（如水下信号、空间信号等）从有色正态噪声中提取信号（如水下信号、空间信号等） 提取非正态信号的相位信息提取非正态信号的相位信息32双谱在目标识别中的应用双谱在目标识别中的应用),(21)(212121222111),(),(),(jrjkyreBecB特性：特性：（1 1）保留了幅值