[发明专利]基于聚类实现任意星座图的MQAM信号调制方式识别的方法

权利要求书

1.基于聚类实现任意星座图的MQAM信号调制方式识别的方法，其特征在于包括如下步骤：

步骤一：信号预处理：主要是通过通用解调算法，从中频信号中解调出符号同步的复基带信号；

步骤二：自适应减法聚类重构星座图，这是第一步聚类算法，主要目的是给出聚类中心的初始个数和初始位置，

符号同步的复基带信号r_k为

rk=Rkejθk+nk,]]>k＝1，2，…，N.                (1)

其中，其中代表信号的调制星座点，n_k代表了高斯噪声，N表示采样序列长度；信噪比为SNR，并且这个值已通过盲的信噪比估计得到；由SNR，算出平均噪声功率P_n。进而算出每个采样点的“密度”指标D_k；

Dk=Σj=1Nexp(-||rk-rj||2Ka*Pn).---(2)]]>

其中，K_a是一个调整系数。根据式(2)的结果，选出D_k的最大值D_c，1对应的采样点r_c，1为第一个聚类中心；然后，对上面的“密度”指标修正；

Dk=Dk-Dc,1exp(-||rk-rc,1||2Kb*Pn).---(3)]]>

其中，K_b是一个调整系数。用同样的方法选出第二个聚类中心r_c，2；

从第三个聚类中心r_c，3开始，每次选出新聚类中心r_c，l+1后，判断聚类中心是否已全部选出；

min{‖r_c，k-r_c，l+1‖²}＞K_c*P_n，k＝1，2，…，l.        (4)

其中，K_b是一个调整系数。如果式(4)成立，则继续进行减法聚类，否则结束；

步骤三、使用模糊C均值聚类算法或C均值聚类算法重构星座图，这是第二步聚类算法，主要目的是基于第一步聚类结果完成星座图的重构；

{r_k，k＝1，2，…，N}是N个基带复信号样本组成的样本集合，为减法聚类估计出的聚类中心个数，m_i，i＝1，2，…为每个聚类的中心，μ_i(r_k)是第k个样本对于第i个聚类的隶属度函数；

根据已知的初始聚类中心，计算样本对所有初始聚类点的隶属度大小，隶属度函数定义如下

μi(rk)=(1/||rk-mi||2)1/(b-1)Σl=1C^(1/||rk-ml||2)1/(b-1),]]>k＝1，2，…，N，i＝1，2，…，C^.---(5)]]>

隶属度计算完成后，计算在此隶属度下的聚类中心值

mi=Σk=1N[μi(rk)]brkΣk=1N[μi(rk)]b,]]>i＝1，2，…，C^.---(6)]]>

在有了聚类中心和隶属度后，计算用隶属度函数定义的聚类代价函数

Jf=Σi=1C^Σk=1N[μi(rk)]b||rk-mi||2.---(7)]]>

其中，b＞1是一个可以控制聚类结果的模糊程度的常数。通过上面的公式(5)(6)(7)反复迭代，使得聚类中心逐渐收敛到聚类代价函数最小的点，此时隶属度值趋于稳定，迭代的结束条件是

J_f(k+1)≈J_f(k).                    (8)

通过以上两部聚类算法可完成星座图重构；

步骤四、使用基于重构星座图的广义似然比测试(GLRT)完成调制方式识别：

首先根据重构出的星座点个数，判断调制方式阶数。如果待识别调制方式集合中有两种调制方式的阶数是一样的，那么继续使用GLRT分类；

通过重构星座点m_i和调制星座点间的欧式距离最近，给出调制星座点的极大似然估计上标c代表了不同的调制类型；然后使用GLRT分类同阶调制类型；表达式如下

其中，表示估计出的调制阶数，表示使l_GLR，m取得最小值的调制方式c。通过式(10)完成了最终的调制方式识别。