Cdma反向数据信道调制类型识别方法
【技术领域】
[000。 本发明涉及信号处理技术领域，具体针对CDMA20001X EV-D0 Rev. A系统（W下简 称CDMA2000系统）反向业务信道中数据信道的调制类型的快速识别方法。
【背景技术】
[0002] 在对通信信号进行调制模式识别时很难找到一个通用的调制分类特征和方法， 针对每种分类问题都必须依据所需分类的调制类型来寻找特定的方法和特征。正如 Lamontagne所说；"调制识别是一个非常直觉的领域，分类特征的选取依赖于作者的知识 背景和想象力"。因此，在通信对抗领域，调制识别一直是被当作一个具有挑战性的问题来 进行研究。
[0003] 通过对信号的解调和参数提取来构造信号的幅度直方图、频率直方图、差分相位 直方图、幅度方差和频率方差等分类特征，然后通过模式识别的分类方法，将选取分类特征 与理想样本的特征参数相比较，按最近原则进行信号自动分类。CDMA2000反向业务信道中 的数据子信道根据用户对传输速率的不同要求采用不同的调制和扩频的形式，因此直接恢 复数据子信道的数据需要遍历各种可能性，增加了计算的复杂性。对CDMA2000反向业务信 道中数据信道调制类型进行识别时，传统的方法需要首先解调出数据信道数据，然后才能 对数据信道的调制类型进行识别。在不知道数据信道使用哪种扩频方式和物理结构的前提 下，需要对各种情况进行遍历才能恢复出数据信道的数据，送样就大量增加了识别过程的 运算量。
【发明内容】
[0004] 本发明针对现有技术存在的上述技术问题，针对传统方法对CDMA2000反向业务 信道中数据信道调制类型进行识别时，需要对各种情况进行遍历才能恢复数据信道的数 据，大量增加了识别过程的运算量的问题，提出一种直接对反向业务信道CDM信号直接进 行处理的识别方法，计算量远小于传统方法，可W实现数据信道调制类型的快速识别。
[0005] 本发明解决上述问题的技术方案是，根据数据信道调制类型与其所使用的扩频码 周期相对应的特点，W CDMA2000系统反向业务信道CDMA信号为基础，使用确定长度对信号 进行分段，求取每段信号的相关矩阵，然后对每段信号的相关矩阵累加取平均获得累加矩 阵，再对累加矩阵使用奇异值分解。根据产生的奇异值的数目判断数据信道使用的扩频码 周期，进而识别出B4、Q2和Q4调制类型；最后结合奇异值分解得到的数据信道扩频码和系 统特点，恢复出数据信道数据，使用高阶累积量识别算法识别出Q4Q2和E4E2类型。
[0006] 一种CDMA2000调制类型识别方法，反向业务信道对CDMA信号求取自相关矩阵进 行奇异值分解，根据奇异值的数目判断数据信道所使用的调制类型和所使用的物理结构； 根据奇异值分解结果和各调制类型结构恢复数据信道数据；对数据信道数据采用高阶累积 量识别算法进一步区分是Q4Q2或E4E2调制类型。
[0007] 所述进行奇异值分解进一步包括；将CDMA信号x(t)按照反向业务信道码片速率 = 1. 2288 Mchip/s)的整数倍采样形成一个时间序列，并w Τ。= 32E/fs为周期进行 连续不重叠分段（其中32是反向业务信道中各信道所使用的Walsh-Hadamard(WH)扩频码 中最长的码序列长度，E为正整数），每段序列构成一个数据向量Xi，取出Μ个数据向量（M 为分段数，其值越大估计出的自相关矩阵越准确），使用送些分段数据向量组成数据矩阵A =[Xi，X2,…，Xm]，求出相关矩阵的估计值
对相关矩阵的估计值 进行奇异值分解。所述判断数据信道所使用的调制类型和所使用的物理结构的方法具体包 括：当数据信道奇异值的数目为0个、8个或16个时，数据信道对应采用B4、Q4、Q2调制类 型；当奇异值数目为24个时，数据信道采用Q4Q2或E4E2调制类型。
[0008] 图1中各信道数据扩频W后的速率均为1. 2288Mchip/s，因此I路各信道在扩频之 前共有3组速率（使用相同扩频码长度的速率相同归为一组），每组速率下又存在1到2个 用户（扩频码长度相同而码序列不同）。由此可见，反向业务信道属于多速率多用户Direct Sequence-Code Division Multiple Access(直接序列码分多址，DS-CDMA)系统模型，由于 I路和Q路模型相同，W下仅W其中一路为例对本发明的实施进行分析说明。
[0009] 可W表达为：
其中也!^为第r组速率的 第k个用户每段数据所包含的信息码数量，R为速率组数，Kf为第r组速率所包含的用户数 目，
为第i段数据的第j个信息码，
巧 应的扩频码序列，P,k为第r组速率的第k个用户采样后的一周期WH码序列，P' 为包 含了 Pf,k的后段的序列，P' f,k,2为包含了 Pf,k的前段的序列。构造四阶累积量特征参数f =C4。I / I C421区分QPSK信号和8PSK信号；当f〉= 0. 6时为QPSK信号；当f <0. 6时为8PSK 信号。
[0010] 本发明提出的一种直接对反向业务信道CDMA信号直接进行处理的识别方法，计 算量远小于传统方法，可W实现数据信道调制类型的快速识别。
【附图说明】
[0011] 图1反向业务信道的物理结构框图；
[0012] 图2B4调制类型结构框图；
[0013] 图3Q4调制类型结构框图；
[0014] 图4Q4Q2调制类型结构框图；
[0015] 图5本发明算法流程框图；
[0016] 图6单个用户数据分段模型；
[0017] 图7各种调制类型的数据x(t)奇异值
[0018] (a)B4调制类型
[001引 （b)Q4调制类型
[0020] (c)Q2调制类型
[0021] (d) Q4Q2/E4E2 调制类型；
[0022] 图8各种调制类型的数据y (t)奇异值
[0023] (a)B4调制类型
[0024] (b)Q4调制类型 [00巧](c)Q2调制类型
[0026] (d) Q4Q2/E4E2 调制类型；
[0027] 图9Q4Q2和E4E2类型对应的特征参数f值。
【具体实施方式】
[0028] CDMA2000反向业务信道的物理结构如图1所示，各子信道经不同的扩频码扩频后 叠加在一起形成I路信道x(t)和Q路信道y(t)。其中，导频信道、辅助导频信道、反向速 率指示信道（RRI信道）、确认/数据源信道（ACK/DSC信道）和数据速率控制信道值RC信 道）的扩频码已知，而数据信道使用的调制类型和扩频码未知。数据信道存在W下5种调 制类型：
[0029] B4调制类型，如图2数据子信道数据经过Binary Phase Shift Ifeying(二相相移 键控，BPSK)调制器后由WH码W/扩频，然后传输到Q路；没有数据传输到I路。Q4调制类 型，如图3数据子信道数据经过如a化ature化ase Shift Ifeying(四相相移键控，QPSK)调 制器后由WH码W/扩频，然后一路传输到I路，另外一路传输到Q路。Q2调制类型与Q4类 型相似，只是Q2类型采用WH码Wi 2进行扩频调制。Q4Q2调制类型，如图4信号W 6比特为 一组进入调制器，前2比特使用QPSK调制后使用WH码W24扩频调制，而后4比特使用QPSK 调制后使用WH码Wi2扩频调制。E4E2调制类型与Q4Q2调制类型相似，信号W 9比特为一 组进入调制器，前3比特使用8如a化ature化ase Shift K巧ing(8相相移键控，8PSK)调制 后使用WH码W24扩频调制，而后6比特使用8PSK调制后使用WH码Wi 2扩频调制。
[0030] 如图5所示为本发明CDMA2000反向业务信道调制类型识别方法处理流程示意图。 将信号x(t)采样形成一个时间序列，并进行连续不重叠分段，使用送些分段数据组成数据 矩阵A，（A = [XI，而，…，Xm])，根据数据矩阵求出相关矩阵的估计值房?Λ/)。对相关矩阵的估 计值Μ)进行奇异值分解，根据产生的奇异值的数目L判断数据信道所使用的调制类型 和所使用的物理结构。数据信道分别采用B4、Q4或Q2调制类型时，分别对应0个、8个或16 个奇异值（图7化）中点1到8为0、8个奇异值，图7 (C)中点1到16为16个奇异值）；而 辅助导频信道、ACK/DSC信道和RRI信道奇异值分解后对应的奇异值数目都是1个（图7 (a) 中点1，图7(b)中点9,图7(c)中点17)、1个（图7(a)中点2,图7(b)中点10,图7(c)中 点18)和2个（图7 (a)中点3到4,图7化）中点11到12,图7 (C)中点19到20)。所此 L的取值分别为4,12, 20。数据信道采用Q4Q2和E4E2调制类型时，由于它们使用相同的扩 频码和类似的物理结构，因此都是对应24个奇异值（此时L为28)，需要使用其它方法识 别。
[0031] 当确定调制类型为Q4Q2或者E4E2时，此时数据信道使用的扩频码和对应的物理 结构就成了已知的，再结合其它信道的扩频码，可W恢复出数据信道的数据，最后使用高阶 累积量识别数据信道的数据是QPSK信号还是8PSK信号，QPSK对应Q4Q2调制类型，8PSK对 应E4E2调制类型。
[0032] 由标准可知，图1中各信道数据扩频W后的速率均为1. 2288Mchip/s，因此I路各 信道在扩频之前共有3组速率（使用相同扩频码长度的速率相同归为一组），每组速率下又 存在1到2个用户（扩频码长度相同而码序列不同）。由此可见，反向业务信道属于多速率 多用户DS-CDM系统模型，由于I路和Q路模型相同，W下仅W其中一路为例对本发明的