一种大气湍流下无线光通信副载波调制星座图识别方法与流程

本发明属于通信技术领域，具体涉及一种大气湍流下无线光通信副载波调制星座图识别方法。

背景技术：

无线光通信普遍采用强度调制/直接检测系统，强度调制可分为脉冲调制和连续波调制两种。副载波调制是一种连续波调制，与ook调制相比，副载波调制不需要自适应阈值检测就可达到最佳差错性能，可以有效抑制大气湍流影响。和ppm调制相比较，也不需要增加系统带宽要求。大气信道是一种有记忆的时变信道，大气散射和大气湍流对副载波调制信号星座图相位影响较大，所导致的相位模糊、相位偏移等给星座图相位识别带来巨大困难。在信息对抗中，调制识别是从已调制信号的数据中提取信号的振幅、频率以及相位等特征信息，对非合作方通信来进行有效的信道、频谱监测、电子对抗以及干扰源识别等。采用信号谱线的调制识别一般是基于信号幅值分布进行，但是实际信号的调制信息是综合表征在幅值和相位分布上，仅利用幅值信息，其性能受限。调制信号星座图提供了信号相位的结构以及各种不同调制状态的关系，是数字调制方式设计与分析的经典工具。因此，采用星座图中的相位信息结合神经网络可有效地进行无线光副载波调制格式的自动识别，尤其可以有效抑制光强起伏对于高阶调制类型检测识别的影响，提高在较强光强振幅起伏方差下的高阶调制识别率。

技术实现要素：

本发明的目的在于提供一种大气湍流下无线光通信副载波调制星座图识别方法，解决了大气湍流光强起伏导致高阶调制识别率低的问题。

本发明所采用的技术方案是：一种大气湍流下无线光通信副载波调制星座图识别方法，具体按照以下步骤实施：

步骤1：建立大气湍流信道数学模型；

步骤2：根据步骤1建立的大气湍流信道数学模型，获得经过大气湍流信道后副载波调制信号的星座图，进行模糊c均值聚类并提取其硬趋势均值特征参数；

步骤3：设计改进三层bp神经网络分类器，将步骤2中提取的特征参数对改进bp网络进行训练，用训练好的网络最后进行调制星座图识别。

本发明的特点还在于，

步骤1建立大气湍流信道数学模型，具体按照以下步骤实施：设x(t)为发射端已调副载波信号，y(t)为接收机输出信号，则y(t)＝i(t)x(t)+n(t)，式中，i(t)是大气信道状态信息，表征为大气湍流强度，弱湍流情况则服从对数正态分布，强湍流情况服从负指数分布或k分布，n(t)为高斯白噪声，均值为0，方差为σ²。

步骤2中进行模糊c均值聚类并提取其硬趋势均值特征参数的具体实施步骤为：

步骤2.1：利用模糊c均值算法把n个待分类的星座图样本点xi(i＝1,2,...,n)划分为c个模糊组，求出每组的聚类中心，使得非相似性指标的价值函数达到最小；利用隶属度定义样本数据隶属于每一类的概率，隶属矩阵u取值为[0，1]之间的实数值，而且样本点隶属于每一类的取值和为1：

若选择欧几里德距离作为非相似性指标时，fcm的价值函数为：

式中的uij为隶属度矩阵，矩阵中元素属于[0，1]区间；ci为模糊组i的聚类中心；第i个聚类中心与第j个星座图数据点之间的欧几里德距离记为dij＝||ci-xj||；m∈[1,∞)是一个加权指数；对所有输入参量求导，可得到使fcm价值函数达到最优解的必要条件为：

通过反复修改隶属度矩阵和聚类中心，使jc(i+1)-jc(i)＜ε，ε取1e-5；

步骤2.2：获得迭代最终的聚类中心和隶属度矩阵后，采用硬趋势函数作为fcm算法的有效性函数，判断将接收信号点划分为c类是否合理，c取2、4、8、16；该函数求得的硬趋势均值特征参数作为三层改进bp神经网络的输入，从而区分不同的mpsk调制类型。

步骤2.2中fcm的硬趋势均值特征参数的计算步骤为：

步骤2.2.1：根据星座图样本点隶属度的最大值，把每个样本点划分到它所属的聚类集合中，那么样本点yi属于第i个聚类ci就定义为下式：

步骤2.2.2：计算每个样本点xj的隶属度矩阵中次大值与最大值之比并记为rj：

步骤2.2.3：当聚类中心为c时，则第i个聚类中心的硬趋势即为所有隶属于ci的样本点，计算rj的均值其中ni为隶属于聚类中心ci的样本点数；

步骤2.2.4：硬趋势tci的平均为所有聚类中心下的均值tc：得到特征参数t2、t4、t8和t16。

步骤3中的分类器采用附加动量法和自适应学习速率相结合的三层改进bp神经网络，其中，附加动量法是在反向传播bp的基础上，给每一个权值的变化上加一项正比于前一次权值变化量的值，再由反向传播法产生新的权值变化；调整自适应学习速率的公式：

式中，k为迭代次数，sse(k+1)和sse(k)为第k+1次和第k次迭代得到的均方误差；若sse(k+1)＜sse(k)表明第k次迭代有效，乘以增量因子1.05来增大学习步长，反之，sse(k+1)＞1.04sse(k)，乘以因子0.7减小步长，加快网络学习速度，最终自适应改变学习速率；

之后设计网络输入层神经元个数取4，隐含层神经元个数取15，输出为2，共四种标识位，网络误差函数ee＝0.01，动量因子0.9，初始学习速率为0.02；利用训练好的改进bp网络再对经大气湍流信道后的调制星座图样本进行识别分类，网络输出设置阈值t1＝t2＝0.5，输出共有4种标识位，00为bpsk，01为qpsk，10为8psk，11为16psk。

本发明的有益效果是：本发明一种大气湍流下无线光通信副载波调制星座图识别方法，技术路线简单、经济而且切实可行，易于实现，提高了大气湍流导致的较大光强起伏下副载波高阶调制的正确识别率，能够很好的抑制大气湍流对星座图识别检测的影响。

附图说明

图1是本发明大气湍流下无线光副载波调制星座图识别系统框图；

图2是本发明所建立的大气湍流信道模型框图；

图3是本发明基于fcm硬趋势均的星座图值聚类特征提取流程图；

图4是本发明在σx²＝0.1时的16psk星座图及其fcm聚类中心图；

图5是本发明四种无线光副载波调制方式下的星座图所对应的聚类硬趋势均值；

图6是本发明改进bp网络训练、验证及测试的均方误差曲线。

具体实施方式

下面结合附图以及具体实施方式对本发明进行详细说明。

本发明提供了一种大气湍流下无线光通信副载波调制星座图识别方法，如图1所示，由光发射机发射的副载波mpsk调制光信号x(t)经过大气湍流信道后，对接收机输出y(t)进行星座图识别。首先采用模糊c均值聚类法对经过大气湍流信道的调制信号星座图进行分析，寻找聚类中心，并迭代出最终的隶属度矩阵uij，再根据隶属度矩阵最大值和次大值求取聚类硬趋势函数的均值，作为不同调制类型星座图的特征参数，设计改进三层bp神经网络分类器，将提取的特征参数对改进bp网络进行训练，用训练好的网络最后进行星座图识别。

步骤1：大气湍流信道模型的建立

大气湍流信道等效的数学模型框图如图2所示。对于无线激光通信系统，系统噪声主要考虑背景光噪声、apd噪声、元器件热噪声等高斯分布加性噪声以及大气湍流导致光强起伏的乘性噪声。图2中x(t)表示发射端已调副载波信号，y(t)表示接收机输出信号，有y(t)＝i(t)x(t)+n(t)，式中，i(t)是大气信道状态信息，表征为大气湍流强度，对于弱湍流情况则服从对数正态分布，强湍流情况服从负指数分布或k分布。i(t)相对于副载波信号频率，属于一种慢衰落，可看作低频的乘性噪声，n(t)为高斯白噪声，均值为0，方差为σ²。随着光强振幅起伏方差的增大，接收端副载波信号星座图样本点愈加弥散，相位信息越模糊，无法分辨。

步骤2：模糊c均值聚类及其硬趋势均值特征提取

步骤a：图3是基于模糊c均值聚类的提取硬趋势均值特征流程图。模糊c均值算法把n个待分类的星座图样本点xi(i＝1,2,...,n)划分为c个模糊组，求出每组的聚类中心，使得非相似性指标的价值函数达到最小。利用隶属度定义样本数据隶属于每一类的概率，隶属矩阵u取值为[0，1]之间的实数值，而且样本点隶属于每一类的取值和为1：

若选择欧几里德距离作为非相似性指标时，fcm的价值函数为：

式中的uij为隶属度矩阵，矩阵中元素属于[0，1]区间；ci为模糊组i的聚类中心。第i个聚类中心与第j个星座图数据点之间的欧几里德距离记为dij＝||ci-xj||；m∈[1,∞)是一个加权指数。对所有输入参量求导，可得到使fcm价值函数达到最优解的必要条件为：

由上述两个必要条件可知，fcm算法是一个迭代过程，通过反复修改隶属度矩阵和聚类中心，使jc(i+1)-jc(i)＜ε，本发明中ε取1e-5。

步骤b：获得迭代最终的聚类中心和隶属度矩阵后，采用硬趋势函数作为fcm算法的有效性函数，判断将接收信号点划分为c(取2、4、8、16)类是否合理，该函数求得的硬趋势均值特征量作为三层改进bp神经网络的输入，从而区分不同的mpsk调制类型。fcm的硬趋势均值计算步骤：

(1)根据星座图样本点隶属度的最大值，把每个样本点划分到它所属的聚类集合中，那么样本点yi属于第i个聚类ci就定义为下式：

(2)计算每个样本点xj的隶属度矩阵中次大值与最大值之比并记为rj：

(3)当聚类中心为c时，则第i个聚类中心的硬趋势即为所有隶属于ci的样本点，计算rj的均值其中ni为隶属于聚类中心ci的样本点数。

(4)硬趋势tci的平均为所有聚类中心下的均值tc：

步骤3：设计改进的三层bp神经网络分类器

本发明识别分类器是采用附加动量法和自适应学习速率相结合的三层改进bp神经网络。此方法在反向传播bp的基础上，给每一个权值的变化上加一项正比于前一次权值变化量的值，再由反向传播法产生新的权值变化。其中调整自适应学习速率的公式：

式中，η(0)为迭代次数k＝0时的初始学习速率。sse(k+1)和sse(k)为第k+1次和第k次迭代得到的均方误差。若sse(k+1)＜sse(k)表明第k次迭代有效，乘以增量因子1.05来增大学习步长，反之，sse(k+1)＞1.04sse(k)，乘以因子0.7减小步长，加快网络学习速度，最终自适应改变学习速率。

因为星座图在经过模糊c均值聚类后计算硬趋势均值为t2、t4、t8和t16，一共四个特征参数，所以本发明设计网络输入层神经元个数取4，隐含层神经元个数取15，输出为2，共四种标识位。网络误差函数ee＝0.01，动量因子0.9，初始学习速率为0.02。利用训练好的改进bp网络再对经大气信道后的调制星座图样本进行识别分类，网络输出设置阈值t1＝t2＝0.5，输出共有4种标识位，00为bpsk，01为qpsk，10为8psk，11为16psk。

本发明的创新点在于：提出了一种新的大气湍流下调制星座图识别方法，并将其应用于副载波调制无线光通信系统中，针对无线光通信四种mpsk副载波调制信号进行了星座图聚类和神经网络识别，提高了在大气湍流下高阶调制识别率。

本发明对大气湍流下无线光mpsk调制信号星座图进行了识别实验。实验中经过大气湍流信道后四种调制(2psk、4psk、8psk和16psk)信号的星座图一共400个样本进行训练和测试，每种调制类型训练和测试样本数均各为50个，星座图的数据点长度为600。对200个训练样本信号星座图采用模糊c均值聚类算法首先进行聚类，其中在σx²＝0.1时的16psk星座图及其fcm聚类中心图如图4所示，然后在不同聚类中心数下迭代出最终隶属度矩阵，求出聚类硬趋势函数均值tc，其中c取2、4、8和16。相对应的tc如图5所示，将tc特征输入到所设计的改进bp网络进行分类识别，对4种共200个学习样本进行训练，经过5次和37次的权重修正后，网络的训练过程结束，训练好的网络再对200个测试样本进行识别实验。bp网络训练、验证和测试曲线见图6所示。其中两种光强起伏方差下的识别结果见表1和表2所示。

表1调制信号星座图识别结果(snr＝20db)

表2调制信号星座图识别结果(snr＝20db)

由上表可知，大气湍流较大光强起伏方差σx²＝0.1时，总正确识别率为100％，当时，总正确识别率也达到96％，且16psk调制识别率达到88％以上。实验结果说明了本发明可以有效抑制大气湍流光强起伏对副载波调制星座图识别的影响，能够提高较大起伏方差下的高阶调制星座图检测识别率。