基于高阶累积量短波通信电台细微特征提取与分析

张广吉，袁仕继，刘志华，黄文晶

(解放军63888 部队，河南 济源 454650)

1 问题描述

人们对信号进行处理，其最终目的就是要对单通道或多通道的有限长观测数据进行必要的处理，从而检测出隐含于这些观测数据中的各种重要信息并加以利用。高阶统计量在信号处理与系统分析中扮演着一个极为重要的角色。根据最近的资料显示，在通信、生物医学工程、语音处理、地震信号分析、图象处理、雷达、声纳等领域都进行了有关高阶统计量处理的研究，具体应用于时延估计、系统辨识、自适应滤波及阵处理等方面。高阶统计量能提供比功率谱更多的有用信息，能够有效地检测信号幅度以外的其他信息，这就具有明显的优点。高阶统计量作为非高斯信号处理的主要分析工具，不仅提供高阶相关信息，而且能够衡量随机序列偏离正态的程度，并对加性高斯噪声(白色或有色)不敏感。这些特点都非常有利于特征参数的提取。

2 基于高阶统计量的信号特征分析

所谓高阶统计量，通常理解为高阶矩、高阶累积量以及它们的谱。本文主要讨论高阶矩和高阶累积量，对它们的谱暂不作分析。首先将利用特征函数来定义高阶矩和高阶累积量，然后引出它们的有关性质，最后推导出高阶矩和高阶累积量之间的变换关系。

2.1 通信信号高阶统计量的算法实现

通过接收机，在不同的时间里对几个不同电台所发射的不同调制样式的短波通信信号进行了不同采样速率的接收。对接收到的信号进行去均值处理，去其直流成分，然后再对信号的高阶统计量(二阶、三阶、四阶)进行计算，并画出曲线图进行比较。

如果，x={xi}，i =1，…，N 表示随机序列或信号的采样，则序列的矩为:，其中E{}表示期望值的运算符。其中p(x)表示xn的概率密度函数。当采集得到的信号时间长度比较长时，认为是各态历经和平稳的，同时进行了去均值处理。

这样则采样信号的矩可由下式估计:

同样它的n 阶累计量:

这样，前4 个累计量可以用下式计算:

2.2 信号的高阶统计量实验计算

对同一电台、同一种调制样式信号进行分析，由于比较得到AM 信号的高阶统计量参数较稳定，分别将两个时间段内对同一电台AM 信号的采样进行了分析画图比较发现，测量的时间不同对同一电台AM 信号的高阶统计量参数影响不大，它们的各阶统计量的波动范围都几乎在同一个很小的区间内，这反应了该电台的同一调制样式的信号的相对稳定性。现以在不同时间内对AM 666 kHz 电台的二阶矩和二阶累计量为例，如图1 所示。

对4 个不同电台的AM 信号进行了采样，它们是800 kHz、864 kHz、936 kHz、1 592.5 kHz。发现它们各自的高阶统计量参数仍是较平稳，不同的是各自的中心值和波动范围不同，以此可以区分各自电台。现将它们的各阶统计量画出，如图2 所示。

对所有电台的高阶统计量参数计算、画曲线图后发现一个共同特点:对一定采样速率得到的采样信号，当对它们合并数目越多时(即新的采样信号的时间长度越长)，信号的高阶统计量参数曲线越平稳，即参数越集中;而对它们截的越短时(即新的采样信号的时间长度越短)，信号的高阶统计量参数曲线越杂乱，抖动越大，即参数越分散。现以电台800 kHz 的三阶统计量为例，如图3 所示。

图1 AM 666 kHz 电台的二阶矩和二阶累计量

图2 不同信号的高阶统计量

以电台864 kHz 的三阶统计量为例，如图4 所示。通过实验结果的比较，得出了上述结论。在对同一电台的同一种调制样式(AM)信号在不同时间内的比较中得到，高阶统计量的值相对稳定，区别也只是随机噪声带来的。这就为利用高阶统计量去识别这一信号的电台提供了实验依据。

同样，在对不同采样时间长度的比较中得到:采样时间越长，高阶统计量越稳定。对于这种现象，可以这样理解:因为在计算信号高阶统计量的时候，所采用的算法的前提条件是信号为各态历经和平稳的，而实际试验中采样的信号，由于噪声的影响和采样时间较短等因素，分析的信号是不太满足上述条件。但是如果采样信号时间越长，它就是越接近各态历经和平稳的，也就越满足前提条件。当合并信号，实际上是在加长采样信号的时间长度，而截短信号，是在减少采样信号时间。而同一电台发射的同一调制样式的信号基本是稳定的，这样当越接近条件时，对实验的效果越好，高阶统计量参数曲线也就越显得平稳。

图3 800 kHz 的三阶统计量

图4 864 kHz 的三阶统计量

3 结束语

由于诸多因素的影响，对通信信号高阶统计量特征的研究还没有深入下去。这篇论文仅是对该方面的一个浅显的探索，只是其中的一小部分。对于下一步工作，可以从以下几个方面开展:

1)解决算法。由于计算信号高阶统计量的计算量太大，总是寻找简便的算法，这样就会出现不少的前提条件，如何使采样得到的信号更满足前提条件，这是需要考虑的;同时，对由于前提条件的不同而产生的高阶统计量的误差，还要进行分析校正。

2)信号样本。由于本文只着重对AM 信号进行了分析，而没对其他的信号作进一步研究，这样所得到的参数还不充足。可以在对各种调制样式的信号采样时，不断改变采样信号的时间长度，看看在什么样的一个时间范围内，信号的高阶统计量可以进一步的稳定，这对信号识别很有借鉴价值。

3)高阶矩谱和高阶累积量谱。在这次分析中，只是针对于信号的时域分析，在对信号归一化(去均值)后利用算法进行的计算。而作为工程上所研究的高阶统计量还包括高阶矩谱和高阶累积量谱，即站在频谱的角度上对高阶统计量进行分析。高阶矩谱和高阶累积量谱会含有对提取信号特征更为有用的信息和规律。

［1］樊昌信. 通信原理［M］. 5 版. 北京: 国防工业出版社，1999.

［2］Richard A.通信电子战系统导论［M］.北京:电子工业出版社，2003.

［3］王炳武.MATLAB5.3 教程［M］.北京:中国水利水电出版社，2000.

［4］张贤达. 现代信号处理［M］. 北京: 清华大学出版社，1995.

［5］姚文冰.基于高阶累积量的抗噪语音识别［D］.武汉:华中科技大学，2001.

［6］梁清泉.基于高阶统计量的波达方向估计技术的研究［D］.西安:西北工业大学，2002.

［7］宋文武，丁志尧，张凯，等. 短波通信系统阻塞干扰控制方法研究［J］.装备环境工程，2010（2）:68-70