本发明属于通信信号分析技术领域，尤其涉及一种mqam信号的调制阶数识别方法，可用于卫星空间通信系统、数字电视网络和数字微波系统。
背景技术：
qam信号的频谱利用率高，噪声容限能力强，被广泛的应用于卫星空间通信，数字电视网络高速数据传输，数字微波系统，vdsl等领域。目前，日本、美国和欧洲的数字有线电视标准均是采用qam调制方式。在通信技术的发展过程中，qam信号衍生出很多新的星座图类型，美国国防部在2011年底颁布的军标mil-std-188-110c附录d中提出一种圆形qam调制方式的星座映射，简称110c(吴天琳.mil-std-188-110c高阶qam信号解调关键技术研究[d].郑州：解放军信息工程大学，2014.)。farbodkayha在2016年提出一种由方形星座图映射而成的圆形星座图qci(farbodkayhan.qamtocircularisomorphicconstellations[c].20168thadvancedsatellitemultimediasystemsconferenceandthe14thsignalprocessingforspacecommunicationsworkshop(asms/spsc)，2016)。星座图综合反映了已调信号幅值和相位的分布信息，是调制阶数识别方面较为重要的依据之一。很多学者对qam信号的调制阶数识别作出研究，ningan等人利用混合矩构造统计量来区分不同调制阶数的qam信号，信噪比5db以上时识别率可以达到100％(ningan,bingbingli,minhuang.researchonmodulationclassificationofmqamsignalsusingjointmoments[c].2010secondinternationalconferenceoncommunicationsoftwareandnetworks，2010)；ningan等人利用混合矩构造统计量，利用fisher准则区分不同调制阶数的qam信号，信噪比5db以上时识别率可以达到98％(ningan,bingbingli,minhuang.modulationclassificationofhigherordermqamsignalsusingmixed-ordermomentsandfishercriterion[c].2010the2ndinternationalconferenceoncomputerandautomationengineering(iccae)，2010，vol3:150-153)；changyiyin等人利用信号的密度谱对信号进行识别，在信噪比10db以上时识别率可以达到100％(changyiyin,bingbingli,yanlingli等,modulationclassificationofmqamsignalsbasedondensityspectrumoftheconstellations[c].20102ndinternationalconferenceonfuturecomputerandcommunication，2010，vol3:57-61)。liyan-ling等人利用粒子群优化和减法聚类，得到不同信号的最佳聚类半径区间，从而对qam信号进行识别，信噪比5db以上时识别率可以达到100％(liyan-ling,libing-bing,yinchang-yi.modulationclassificationofmqamsignalsusingparticleswarmoptimizationandsubtractiveclustering[c].ieee10thinternationalconferenceonsignalprocessingproceedings，2010：1537-1540)；hishamabuella等人利用核密度估计对不同qam信号进行识别，在信噪比大于6db时识别率可以达到100％(hishamabuella,mehmetkemalozdemir.automaticmodulationclassificationbasedonkerneldensityestimation[j].canadianjournalofelectricalandcomputerengineering，2016，39(3):203-209)。这些研究方法在圆形qam信号调制阶数识别性能中降低。
综上所述，现有技术存在的问题是：目前大多数qam信号阶数识别方法仅仅局限于星座图形状为方形和十字形的qam信号，对于圆形qam信号调制阶数的识别将失效。
技术实现要素：
针对现有技术存在的问题，本发明提供了一种mqam信号的调制阶数识别方法、卫星空间通信系统。
本发明是这样实现的，一种mqam信号的调制阶数识别方法，所述mqam信号调制阶数识别方法包括以下步骤：
步骤一，对接收到的mqam信号做下变频处理恢复出基带信号，得到星座图并归一化，利用网格法估计星座图的概率密度函数；
步骤二，利用图像灰度化和改进的二值化对所估计的概率密度函数的俯视图进行图像处理，得到简化的星座图概率密度函数；
步骤三，利用多半径减法聚类和自校正的聚类有效性函数对简化的星座图的概率密度函数进行处理，识别出mqam信号的调制阶数。
进一步，星座图的概率密度函数定义为：经过功率归一化的第i个星座点ri，有ri＝si+vi，i＝1,2,...,l，其中si为发送端发送的星座点，坐标为(xk,yk)，vi为加性高斯白噪声对星座图的影响，星座点的偏移量，l为接收到的星座点个数，则ri的概率密度分布函数为：
各个星座点独立等概发射，当调制阶数为m时，每个星座点的传输个数为l/m，所有发射端发射的与si坐标相同的星座点，在接收端的概率密度分布函数为
化简得
进一步，所述网格法估计星座图的概率密度函数的方法为：根据接收到的星座点ri的坐标为(rii,riq)，构造方形网格，网格的边界由横纵坐标绝对值的最大值决定，ri,rq＝±max(rii,riq)；定义网格密度meshnum为单一方向上划分的网格数，在确定的区域内划分meshnum*meshnum的网格图，每个网格大小为
进一步，所述利用改进的二值化对所估计的概率密度函数的俯视图进行图像处理方法为：灰度空间，结合星座图的概率密度函数俯视图，对其灰度空间设定相应的阈值gth，将符合条件的像素点归一化灰度值为1，不符合条件的像素点归一化灰度值为0：
其中thresholdx,y为二值化后图片上(x,y)像素点的归一化灰度值；在得到的灰度图像上搜索灰度值最大的像素点(x,y)，归一化灰度值为gmax，令阈值gth＝pgmax。
进一步，自校正的聚类有效性函数：r＝{ri＝(rii,riq),i＝1,2,...,l}为n个点所组成的集合，c为计算出的聚类中心个数，mk为第k个聚类中心，k＝1,2,...,c，计算有效性函数s的步骤如下：
1)计算每个点ri的sik：
其中ri对应的聚类中心为mk，b(q,i)为第i个点ri和划分到第q个聚类中心mq的所有信号点的平均距离，其中q为变量且q≠k，a(i)为第i个点ri和在同一聚类中心下的其他信号点的平均距离；
2)计算mk中所有点的sik平均值sk：
其中nk为聚类中心mk中的样本点数：
3)计算所有聚类中心sk的均值s：
4)计算调制阶数mdec：
在第a个聚类半径下得到的聚类有效性函数为s(a),a＝1,2，3，实际计算出的s值越接近1表明接收信号划分为该类越合理，定义自校正的聚类有效性函数来对聚类得到的结果进行反馈和校正：
定义自校正的聚类有效性函数：
其中ma为第a次减法聚类得到的聚类中心个数，δsa＝|sa-smax|为第a次减法聚类的有效性函数sa和所有减法聚类得到的有效性函数最大值smax差的绝对值，ε为一个足够小的值，取ε＝0.01，n为累加的个数，roundms(x)表示在ms序列中取最接近x的值，其中ms＝[16,32,64]。
本发明的另一目的在于提供一种应用所述mqam信号的调制阶数识别方法的卫星空间通信系统。
本发明的另一目的在于提供一种应用所述mqam信号的调制阶数识别方法的数字电视网络。
本发明的另一目的在于提供一种应用所述mqam信号的调制阶数识别方法的数字微波系统。
本发明的优点及积极效果为：对mqam信号的星座图进行处理，识别其调制阶数；对于不同形状的16qam信号，当信噪比大于4db时，识别率达到100％；对于十字形32qam信号，当信噪比大于等于5db时，识别率达到100％；对于110c圆形32qam信号，当信噪比大于等于7db时，识别率达到100％；对于方形64qam信号，当信噪比大于等于10db时，识别率达到100％；对于110c圆形64qam信号，当信噪比大于等于10db时，识别率达到90％以上，由此可见，本发明的识别效果良好。
附图说明
图1是本发明实施例提供的mqam信号的调制阶数识别方法流程图。
图2是本发明实施例提供的16qam信号的调制阶数识别性能示意图。
图3是本发明实施例提供的32qam信号的调制阶数识别性能示意图。
图4是本发明实施例提供的64qam信号的调制阶数识别性能示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示，本发明实施例提供的mqam信号的调制阶数识别方法包括以下步骤：
s101：对接收到的mqam信号做下变频处理恢复出基带信号，得到星座图并归一化，并利用网格法估计星座图的概率密度函数；
s102：利用图像灰度化和改进的二值化对所估计的星座图概率密度函数的俯视图进行图像处理，得到简化的星座图概率密度函数；
s103：利用多半径减法聚类和自校正的聚类有效性函数对简化的星座图的概率密度函数进行处理，识别出不同星座图类型的mqam信号的调制阶数。
下面结合附图对本发明的应用原理作进一步的描述.
本发明实施例提供的mqam信号的调制阶数识别方法包括以下步骤：
(1)对接收到的mqam信号做下变频处理恢复出基带信号，得到星座图并归一化，利用网格法估计星座图的概率密度函数；
其中星座图的概率密度函数定义为：
针对经过功率归一化的第i个星座点ri，有ri＝si+vi，i＝1,2,...,l,其中si为发送端发送的星座点，其坐标为(xk,yk)，vi为加性高斯白噪声对星座图的影响，即星座点的偏移量，l为接收到的星座点个数，则ri的概率密度分布函数为：
假设各个星座点独立等概发射，当调制阶数为m时，每个星座点的传输个数为l/m，所有发射端发射的与si坐标相同的星座点，在接收端的概率密度分布函数为
化简得
通过网格法估计星座点的概率密度函数，包括：
根据接收到的星座点ri的坐标为(rii,riq)，构造方形网格，网格的边界由横纵坐标绝对值的最大值决定，即ri,rq＝±max(rii,riq)。定义网格密度meshnum为单一方向上划分的网格数，即在确定的区域内划分meshnum*meshnum的网格图，每个网格大小为
(2)对(1)中得到的星座图概率密度函数的俯视图利用图像灰度化和改进的二值化进行处理，得到简化的概率密度分布函数部分按以下进行：
首先对s1中得到的概率密度函数进行灰度化处理，灰度值加权公式为：
gray＝0.299*r+0.587*g+0.114*b；
其中，gray表示灰度值，r表示红色，g表示绿色，b表示蓝色。
再利用改进的二值化对灰度图像进行处理，得到简化的概率密度函数，处理方法为：
针对灰度空间，结合星座图的概率密度俯视图，对其灰度空间设定相应的阈值gth，将符合条件的像素点归一化灰度值设定为1，不符合条件的像素点归一化灰度值设定为0，即：
其中thresholdx,y为二值化后图片上(x,y)像素点的归一化灰度值。在得到的灰度图像上搜索灰度值最大的像素点(x,y)，其归一化灰度值为gmax，令阈值gth＝pgmax。
(3)对由(2)得到的简化的概率密度分布函数进行多半径减法聚类并利用自校正的聚类有效性函数来进行阶数识别，其中自校正的聚类有效性函数定义为：
r＝{ri＝(rii,riq),i＝1,2,...,l}为n个点所组成的集合，c为计算出的聚类中心个数，mk为第k个聚类中心，k＝1,2,...,c，则计算有效性函数s的步骤如下：
1)计算每个点ri的sik：
其中ri对应的聚类中心为mk，b(q,i)为第i个点ri和划分到第q个聚类中心mq的所有信号点的平均距离，其中q为变量且q≠k，a(i)为第i个点ri和在同一聚类中心下的其他信号点的平均距离。
2)计算mk中所有点的sik平均值sk：
其中nk为聚类中心mk中的样本点数。
3)计算所有聚类中心sk的均值s：
4)计算调制阶数mdec：
在第a个聚类半径下得到的聚类有效性函数为s(a),a＝1,2，3，实际计算出的s值越接近1表明接收信号划分为该类越合理，也可以理解为划分出的聚类中心越接近其真实调制阶数，因此，定义自校正的聚类有效性函数：
其中ma为第a次减法聚类得到的聚类中心个数，δsa＝|sa-smax|为第a次减法聚类的有效性函数sa和所有减法聚类得到的有效性函数最大值smax差的绝对值，ε为一个足够小的值，n为累加的个数，roundms(x)表示在ms序列中取最接近x的值，其中ms＝[16,32,64]。
下面结合仿真对本发明的应用效果作详细的描述。
为了评估本发明的性能，下面的仿真实验采用mqam信号，并进行2000次实验，mqam信号的参数设置如下：信号集φ＝{16qam的方形、110c、qci星座图；32qam的十字形、110c星座图；64qam的方形、110c、qci星座图}。码元速率为fb＝1mhz，采样频率fc＝8mhz。数据长度为50000点。其仿真结果如图2、图3和图4所示，对于不同形状的16qam信号，当信噪比大于4db时，识别率达到100％；对于十字形32qam信号，当信噪比大于等于5db时，识别率达到100％；对于110c圆形32qam信号，当信噪比大于等于7db时，识别率达到100％；对于方形64qam信号，当信噪比大于等于10db时，识别率达到100％；对于110c圆形64qam信号，当信噪比大于等于10db时，识别率达到90％以上，由此可见，本发明具有良好的识别效果。本发明可用于卫星空间通信系统、数字电视网络和数字微波系统。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。