[发明专利]一种基于状态优势统计的MQAM信号调制识别方法

权利要求书

1.一种基于状态优势统计的MQAM信号调制识别方法，其特征在于：

获得的MQAM信号的复基带数据为：

r(n)＝I(n)+j×Q(n)n＝1，2，...N

其中，n为当前的样点序号，N为总的采样点数；具体包括如下步骤：

步骤1、定义模板：若系统的MQAM调制类型共有M种，构成集合A，其对应的阶数分别为则确定模板k＝[k₁，k₂...k_M]，k_i的值为：

ki=k2mni=2mk2m+12n×k5ni=2m+1,ni≤5ni=5+2ni=1,2,...M,---(1)]]>

其中，k_2m为满足条件的最小正整数，k_2m+1为满足条件的最小正整数，m和n均为正整数；

步骤2、选取最低标准频次P_std和统计门限K_th；

步骤3、将r(n)映射到星座图，得到星座点坐标为：

R(n)＝I_R(n)+j×Q_R(n)n＝1，2，...N

然后判断I_R(n)与模板k中相匹配的值K_x，具体方法为：

a.令初始索引i_x＝1；

b.取统计I_R(n)n＝1，2，...N落入区间[k22(1-kth),k22(1+kth)],......,[kix2(1-kth),kix2(1+kth)]]]>内的样点数为NR={N1,N2,...Nix},]]>则各区间频次为p={p1,...,pix}={N1N,N2N,...NixN};]]>求全局频次P_g：

pg=1kixΣi=1ixpikix=2n1kix(Σi=12n+1pi+32Σi=2n+1+1ixpi)kix=2n+2n+1,---(2)]]>

其中，n为0或正整数；

c.计算tmp(i_x)：

为满足的最大整数；

d.若pg>pstdtmp(ix),]]>则判定kx=kix;]]>若pg≤pstdtmp(ix),]]>则令i_x＝i_x+1，重复

b至d的步骤；

步骤4、判断Q_R(n)与模板k中相匹配的值k_y，具体方法为：

a.令初始索引i_y＝1；

b.取统计Q_R(n)n＝1，2，...N落入区间[k22(1-kth),k22(1+kth)],......,[kiy2(1-kth),kiy2(1+kth)]]]>内的样点数为NR={N1,N2,...Niy},]]>则各区间频次为p={p1,...,piy}={N1N,N2N,...NiyN},]]>求全局频次P_g：

pg=1kiyΣi=1iypikiy=2n1kiy(Σi=12n+1pi+32Σi=2n+1+1iypi)kiy=2n+2n+1,---(4)]]>

其中，n为0或正整数；

c.计算tmp(i_y)：

为满足的最大整数；

d.若pg>pstdtmp(iy),]]>则判定kx=kiy,]]>若pg≤pstdtmp(iy),]]>则令i_y=i_y+1，重复b至d的步骤；

步骤5、利用步骤3和步骤4得到的k_x和k_y，求调制阶数M：

M=4×kx×kykx=ky=2n02tkx=ky=2n0+2n0+1Mpkx≠ky,---(6)]]>

其中，t为满足条件2t≤4×k_x×k_y的最大整数，n₀为0或正整数，M_p为一个与n_M有关的变量；求满足2nM<4×kx×ky<2nM+1]]>的n_M，若M1=2nM]]>和M2=2nM+1]]>均属于A，则M_p各以0.5的概率取M₁或M₂；否则，M_p＝min‖k_i-4×k_x×k_y‖²i=1，2，...M；

从而得到调制阶数M。

2.根据权利要求1所述的一种基于状态优势统计的MQAM信号调制识别方法，其特征在于：p_std和K_th是与信噪比有关的经验值，0＜p_std＜1，0＜k_th＜1。