基于分数低阶循环谱的图域通信信号调制识别方法

摘要

本发明公开了一种基于分数低阶循环谱的图域通信信号调制识别方法。利用接收信号的三维分数低阶循环谱，将被α稳定分布噪声干扰的调制信号转换到图域上，然后可以从图表示的稀疏邻接矩阵中提取有效特征参数行索引序列集合作为调制类型的特征，根据训练信号与接收信号的行索引序列集合汉明距离，来实现α稳定分布噪声干扰下，更稳定的更有效的通信信号调制类型的识别。

说明书

技术领域

本发明属于信号处理技术领域，更为具体地讲，涉及一种基于分数低阶循环谱的图域通信信号调制识别方法。

背景技术

自动调制分类(Automatic Modulation Classification，简称AMC)，也称通信信号调制识别可以在很少或没有先验知识的情况下识别接收信号的调制类型，是信号检测和解调之间必不可少的一个重要步骤，并广泛应用于许多军事和民用通信领域。

经典自动调制分类(AMC)方法，通常可分为两类：(i)基于似然的(LB)决策理论方法和(ii)基于特征的(FB)模式识别(PR)方法。然而，LB的方法不可避免地会有一些缺点，例如缺乏闭式解，难以忍受的高计算复杂度，概率模型不匹配。FB方法的性能不是最佳的，然而它们能非常有效地实现，因此，许多研究利用不同的特征和不同的分类算法以追求FB方法的鲁棒性能。

值得注意的是，LB方法和FB方法都是应用在高斯噪声信道的假设中，然而，各种各样的研究表明，在实际的无线通信中频道，通常是由明显的脉冲引起的多址干扰，低频大气噪声，电磁干扰等。这些物理噪声表现出尖锐的脉冲特性和具有重尾的概率密度分布。根据中心极限定理，这些在无线通信系统中主要误差来源的非高斯分布噪声可以被建模为α稳定分布噪声。在α稳定分布噪声出现的信道中，传统AMC方法性能会出现明显的恶化。

基于图域的自动调制分类(AMC_G)第一次将AMC技术引入图形域，并且已经实现了比现有PR和基于LB的决策理论算法更优的性能。但是该方法是对接收信号的二阶循环谱进行图域映射提取图域特征。然而在α稳定分布噪声中不存在二阶及更高阶的统计量，所以，现有的AMC_G方法在α稳定分布噪声中也失效，因此，新的更稳定的更有效的适用于α稳定分布噪声的AMC技术亟待被发现。

发明内容

本发明的目的在于克服现有技术的不足，提出一种基于分数低阶循环谱的图域通信信号调制识别方法，以适应α稳定分布噪声，实现更稳定的更有效的通信信号调制类型的识别。

为实现上述发明目的，本发明基于分数低阶循环谱的图域通信信号调制识别方法，其特征在于，包括以下步骤：

(1)、调制类型训练信号的特征提取

1.1)、基于分数低价循环谱的图域映射

对于无噪声的第k类调制类型的训练信号x^k(t)，k＝1,2,…,K，K为调制类型的类型数量；将其采样序列划分为L段，每一段进行一次图域映射：

采用FAM算法((Fast Fourier transform)Accumulation Method)：FFT累加算法，用于计算循环谱密度)计算出l段训练信号的FLOCS(Fractional Low-Order CyclicSpectrum，分数低阶循环谱)，得到图域集合：

其中，h＝1,2...H，表示第k类调制类型的训练信号的l段保留下来的循环频率ε_h所对应的时域平滑循环周期图，提取出H个循环频率ε_h所对应的时域平滑循环周期图的邻接矩阵，得到邻接矩阵集合：

其中，时域平滑循环周期图根据以下方式得到：

a1)、对计算出的FLOCS即分数低阶循环谱进行归一化和量化处理，得到最大值为1且离散的分数低阶循环谱其中，ε为循环频率，f为频率；

在FAM算法中，FLOCS的频率分辨率为Δf＝f_s/N′，循环频率分辨率Δα＝1/Δt＝f_s/N，其中，f_s为采样频率，N′为复解调所用数据的点数，N为Δt时间内输入的数据点数，这样采用FAM算法计算出的FLOCS为(N′+1)×(2N+1)的矩阵；

a2)、由于FLOCS具有对称性，基于离散的分数低阶循环的四分之一象限建立相应的图域映射：

定义稳定的循环频率ε_p,p＝1,2...N，ε_p满足条件：

将稳定的循环频率ε_p,p＝1,2...N相应的频率值作为顶点，得到顶点集合：

将两个顶点之间的幅度差值作为边，得到边集合：

其中：

这样，在每一个稳定的循环频率ε_p下得到相应的图域映射，即时域平滑循环周期图为：

将分数低阶循环谱为0的循环频率删除，得到H个保留下来的循环频率ε_h所对应的时域平滑循环周期图：

1.2)、行索引序列的提取

对于每个邻接矩阵l＝1,2,…,L提取主对角线正上方的次对角线的非零条目(元素)，提取这些非零条目(元素)所对应的行索引序列行索引序列提取的原则如下：

b1)、检查次对角线的非零值，列出这些非零值所对应的行索引，并根据这些非零值的绝对值对这些行索引进行降序排列，然后，按降序依次提取行索引；

b2)、如果两个或多个非零条目具有相同的绝对值，则提取距离之前所提取的行索引距离最近的行索引，其他的丢弃；

b3)、如果两个或多个非零条目具有相同的绝对值，且最大，则选择最大的行索引，其他的丢弃；

这样得到循环频率ε_h所对应的得到L个行索引序列，选取在L个行索引序列中出现概率大于95％行索引构成一个稳定的行索引序列

对于第k类调制类型的训练信号，提取出H个循环频率ε_h稳定的行索引序列，构成稳定行索引序列集合：并作为第k类调制类型的特征；

(2)、通信信号调制类型的识别

对于接收信号，按照步骤(1)的方法获取其调制类型的特征，行索引序列集合其中，V是保留下来的循环频率个数；

计算行索引序列集合与第k类调制类型的特征的汉明距离，得到K个汉明距离k＝1,2,…,K，然后在其中找最小的汉明距离，其对应的调制类型即为接收通信信号的调制类型。

本发明的目的是这样实现的。

为应对α稳定分布噪声，本发明基于分数低阶循环谱的图域通信信号调制识别方法。利用接收信号的三维分数低阶循环谱，将被α稳定分布噪声干扰的调制信号转换到图域上，然后可以从图表示的稀疏邻接矩阵中提取有效特征参数行索引序列集合作为调制类型的特征，根据训练信号与接收信号的行索引序列集合汉明距离，来实现α稳定分布噪声干扰下，更稳定的更有效的通信信号调制类型的识别。

附图说明

图1是本发明应用的一种具体实施方式原理框图。

具体实施方式

下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。

为了方便描述，先对具体实施方式中出现的相关专业术语进行说明：

AMC(automatic modulation classification)：自动调制分类；

FB(feature-based)：基于统计特征

PR(pattern recognition)：模式识别

LB(Likelihood-based influence)：基于似然函数

AMC_G(graph-based>

PDF(probability density function)：概率密度函数；

CF(characteristic function)：特征函数；

FLOCS(fractional low-order cyclic spectrum)：分数低阶循环谱；

FLOC(fractional low-order correlation)：分数低阶自相关函数；

FLOCC(fractional low-order cyclic correlation)：分数低阶循环自相关；

FAM(FFT(fast Fourier transform)accumulation method)：FFT累加算法，用于计算循环谱密度；

BPSK(binary phase-shift keying)：二进制相移键控；

QPSK(quadrature phase-shift keying)：正交相移键控；

OQPSK(offset quadrature phase-shift keying)：偏移四相相移键控；

2FSK(binary frequency-shift keying)：二进制频移键控；

4FSK(quadrature frequency-shift keying)：四进制频移键控；

MSK(minimum shift keying)：最小频移键控；

1、α稳定分布

α稳定分布又称为非高斯稳定分布、重尾分布，是一种广义的高斯分布，这种分布模型可以在实际的无线通信环境中，准确地模拟噪声的统计特性。

α稳定分布的模型是唯一满足稳定性和广义中心极限定理的模型，α稳定分布并不存在统一、封闭的概率密度函数(PDF)，但它存在统一的特征函数(CF)，其特征函数可以表示为：

ψ(u)＝exp{jau-γ|u|^α[1+jβsgn(u)ω(u,α)]}(9)；

其中，sgn(·)为符号函数。α(0＜α≤2)为特征指数，它决定该分布脉冲特性程度，α值越小，所对应分布的拖尾越厚，因此脉冲特性越显著；β(-1≤β≤1)为偏斜参数，用于确定分布的对称程度；γ(γ＞0)为分散系数，又称尺度参数，它是关于样本偏离其均值的分散程度的度量，类似于高斯分布中的方差；α(-∞＜a＜+∞)为位置参数，对应于稳定分布的均值或中值，u为特征函数的随机变量。

当α＝2时，α稳定分布退化为高斯分布；

当α＝1且β＝0时，α稳定分布为柯西分布；

当β＝0时，α稳定分布为关于均值α的对称分布，我们称这样的分布记为对称α稳定(SαS)分布。

2、分数低阶循环谱(FLOCS)分析

由于调制信号s(t)被服从α稳定分布的噪声n(t)所污染，因此接收信号x(t)可以被建模为：

x(t)＝s(t)+n(t)(11)；

其中，n(t)为服从SαS分布的噪声，由于α稳定分布的噪声具有显著的尖峰脉冲特性，不具有二阶或二阶以上统计量，传统的基于二阶或高阶循环统计量的AMC算法在α稳定分布的噪声会失效，对接收信号进行非线性变换得到的分数低阶循环谱(FLOCS)可以有效的抑制α稳定分布的噪声，因此，对于AMC技术，可以从接收信号的FLOCS中提取相应的信息进行调制信号的识别。

图1是本发明应用的一种具体实施方式原理框图。

在本实施例中，输入数据在发射机中经过调制器调制后，得到调制信号s(t)，然后在信道中混入的α稳定分布噪声n(t)，成为接收机的接收信号x(t)。

首先，在自动调制分类器中，对接收信号x(t)以采样频率F_s＝1/T_s进行均匀采样，采样后的离散信号x(n)的分数低阶自相关函数(FLOC)可以被表示为：

FLOC(n,m)＝E{[x(n+m)]^{b}[x^*(n)]^{b}}(12)；

x(n)^{b}＝|x(n)|^b-1x^*(n)(13)；

其中，式(12)是对散信号x(n)的b阶非线性变换，0＜b＜α/2；E(·)为期望，x^*(n)是x(n)的共轭。那么，信号的分数低阶循环自相关(FLOCC)为：

其中，<·>表示时间平均，值得注意的是，b阶非线性变换只改变了信号的幅度，没有改变周期信息，所以二阶循环相关下定义的循环频率同样适合分数低阶循环相关；若b＝1，则FLOCC退化为二阶循环自相关。FLOCS为FLOCC的傅里叶变换，可以被表示为：

实际上，可以利用时域平滑算法——FAM算法估计出，对于一个给定的频率f和循环频率ε，时域平滑循环周期图可以由下式表示：

其中g(n)是宽度为NT_s秒的统一权重函数，f₁和f₂是FAM算法中滤波器的中心频率，T_s是采样周期，其中，f₁＝f+α/2，f₂＝f-α/2，X_T(r,f₁)和X_T(r,f₂)是x(n)的复解调，可以由下式计算出。

其中a(r)是持续时间为T＝N′T_s秒的锥形数据窗，它的宽度即是FLOCS的频率分辨率Δf，如果a(r)是归一化的，FLOCS可以由时域平滑周期图实现无偏估计，如下式：

3、图域映射

采用FAM算法计算出的三维图的幅度为非负的，并对计算出的FLOCS进行归一化和量化处理，得到最大值为1且离散的分数低阶循环谱在FAM算法中，FLOCS的频率分辨率为Δf＝f_s/N′，循环频率分辨率Δα＝1/Δt＝f_s/N，其中，f_s为采样间隔，N′为复解调所用数据的点数，N为Δt时间内输入的数据点数。即采用FAM算法计算出的FLOCS矩阵(N′+1)×(2N+1)的矩阵。

由于FLOCS具有对称性，因此对离散谱的四分之一象限建立相应的图域映射。定义稳定的循环频率ε_p,p＝1,2...N，ε_p满足条件：

将稳定的循环频率相应的频率值作为顶点，设为：将两个顶点之间的幅度差值作为边，设为：q₁,q₂＝0,1,...,N′/2}，其中：

至此，可以在每一个稳定的循环频率下得到相应的图域映射p＝0,1,...,N，显然每个循环频率下的图具有循环性，因此可以提取相应图的邻接矩阵作为不同信号的判别特征。

4、提取特征

设调制类型集合为其中，表示第k类调制类型，k＝1,2,...,K。在本实施例中，可以对6类调制类型信号进行识别，即BPSK，2FSK，4FSK，QPSK，OQPSK，MSK，对于无噪声的第k类调制类型训练信号，可以计算其FLOCS，根据第3部分的方法构建图域集。

将训练信号的采样序列划分为L段，可以建立L次图域映射，对于每一次图域映射，可以得到H个图，对于第l次图域映射，图域的集合可以表示为其中h＝1,2...H，表示第k种调制类型的训练信号保留下来的循环频率ε_h所对应的图，其对应的图提取出的邻接矩阵集合表示为

因为FLOCS在图域中代表一个加权的有向环，任意邻接矩阵是如下性质的稀疏矩阵

其中，为邻接矩阵的第(u,v)个条目，对于每个邻接矩阵提取邻接矩阵主对角线正上方的次对角线的非零条目，提取这些非零条目所对应的行索引序列行索引序列提取的原则如下：

b1)、检查次对角线的非零值，列出这些非零值所对应的行索引，并根据这些非零值的绝对值对这些行索引进行降序排列，然后，按降序依次提取行索引；

b2)、如果两个或多个非零条目具有相同的绝对值，则提取距离之前所提取的行索引距离最近的行索引，其他的丢弃；

b3)、如果两个或多个非零条目具有相同的绝对值，且最大，则选择最大的行索引，其他的丢弃；

这样得到循环频率ε_h所对应的得到L个行索引序列，选取在L个行索引序列中出现概率大于95％行索引构成一个稳定的行索引序列

对于第k类调制类型的训练信号，提取出H循环频率ε_h稳定的行索引序列，构成稳定行索引序列集合：并作为第k类调制类型的特征。

注意，这些行索引序列不必具有相同个数的元素，因为每个序列的长度由相对应的邻接矩阵的非零元素决定。

5、通信信号调制类型的识别

对于接收信号，按照第3、4部分的方法获取其调制类型的特征，行索引序列集合其中，V是保留下来的循环频率个数；

计算行索引序列集合与第k类调制类型的特征的汉明距离，得到K个汉明距离k＝1,2,…,K，然后在其中找最小的汉明距离，其对应的调制类型即为接收通信信号的调制类型。

在本实施例中，如图1所示，接收信号x(t)进行预处理后送入分类器中按照前述第5部分的方法进行通信信号调制识别，并把调制类型送入解调器中，按照对应的调制类型对预处理后的接收信号进行解调，得到输出数据。

如图1所示，本发明通过计算分数低阶循环谱FLOCS将被α稳定分布噪声干扰的调制信号转换到图域上，然后通过图域映射及特征提取，得到调制类型训练信号的特征，然后根据特征进行图域分类，实现α稳定分布噪声干扰下，更稳定的更有效的通信信号调制类型的识别。

尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。