[发明专利]一种基于双谱和稀疏矩阵的数字调制信号分类方法

权利要求书

1.一种基于双谱和稀疏矩阵的数字调制信号分类方法，其特征在于，包括以下步骤：

（1）、稀疏矩阵字典的建立

1.1）、分别计算确知的四种数字调制信号2FSK、4FSK、QPSK、16QAM的训练序列x_2FSK(t),x_4FSK(t),x_QPSK(t),x_16QAM(t)的双谱，得到其双谱表示，对每个序列进行以下处理：

将长度为N的训练序列分割成M段，计算每个分段的双谱

Y(i)(λ)=1KΣk=0K-1xi(k)exp(-j2πkλ/K),]]>λ=λ₁+k₁、λ₂+k₂或λ₁+λ₂+k₁+k₂；

其中L₁是预设整数，Δ₀是双谱域的频率采样间隔，K=N/M，xⁱ(·)表示数字调制信号的第i个部分；

训练序列的双谱采用均值定义，表示为

这样四种数字调制信号2FSK、4FSK、QPSK、16QAM的训练序列x_2FSK(t),x_4FSK(t),x_QPSK(t),x_16QAM(t)的双谱为：以及

1.2）、对四种数字调制信号的双谱分别进行稀疏表示，得到各自的稀疏矩阵ψ：

B‾(λ1,λ2)=ψs0,]]>其中：

ψ=[v₁,v₂,...v_n,]∈Rⁿ

Rⁿ是n维实数域，vj是第j个数据，j=1,2…,K；

s0=.[0,···,0,sd1,sd2,···,sdm,0,···,0]nT]]>

由于维度为n相对较高，计算复杂度高，需要进行降维处理：

用矩阵Ф同时对两边进行相乘，得到：

y=.ΦB‾(λ1,λ2)=Φψs0∈Rf]]>

R^f表示f(f<n)维的实数域；

这样得到四种数字调制信号2FSK、4FSK、QPSK、16QAM的稀疏矩阵为：

y_2FSK,y_4FSK,y_QPSK,y_16QAM，并构成稀疏矩阵字典；

（2）、数字调制信号调制方式的分类

2.1）对输入待分类的数字调制信号进行步骤（1）同样的变换，得到变换结果

2.2）、计算输入待分类数字调制信号变换后结果与字典中每个稀疏矩阵的残差：

rd(y~x)=||y~x-yd||=||y~x-ΦB‾(λ1,λ2)||2,d=1,2,3,4]]>

d=1的时候y_d=y_2FSK,d=2的时候y_d=y_4FSK,

d=3的时候y_d=y_QPSK,d=4的时候y_d=y_16QAM,

为了提高分类的稳定性，采用残差均值,得到其残差均值为：

E[rd]=1lΣj=1,···,lrd]]>

其中，l表示对同一待分类数字调制信号的测试次数；

2.3）、鉴别待分类数字调试信号的调制方式：对哪种调制方式的稀疏矩阵的残差均值E[r_d]最小，就确定为该调制类别。