CN115834313A - 一种基于帧结构的gpu并行计算qpsk相干解调方法 - Google Patents

Abstract

本发明涉及GPU并行技术技术领域，具体涉及一种基于帧结构的GPU并行计算QPSK相干解调方法，主要包括位同步和载波同步以及符号判决，一帧数据采用并行处理方案，位同步误差和载波同步误差一帧更新一次，一帧采用相同的位同步插值位置和载波同步相位补偿。当信道条件良好，频偏波动很小时，帧长可以很大，并行度增大运行速度随之增大；由于gardner位同步算法可以独立于载波同步运行，故本方法先对IQ数据进行位同步，位同步后一个符号一个采样点，数据量大幅减小，随后给载波同步模块；载波同步采用鉴频和鉴相相结合的方案，使环路锁定更稳定，抗频偏能力更强。

Description

一种基于帧结构的GPU并行计算QPSK相干解调方法

技术领域

本发明涉及GPU并行技术技术领域，具体涉及一种基于帧结构的GPU并行计算QPSK相干解调方法。

背景技术

传统的QPSK解调大多采用专门的芯片或FPGA实现，随着近几年GPU芯片的发展，GPU应用范围越来越广，使得GPU应用到QPSK实时解调成为可能，凭借GPU强大的并行运算能力，在中低符号率的场合逐步得到应用。且GPU开发具有开发周期短，灵活修改的特点，近几年GPU用于解调越来越受到重视。

QPSK相干解调需要载波同步和位同步，载波同步原理框图如图1所示，AD数据经过正交下变频和匹配滤波器得到基带IQ数据，基带IQ数据经过载波同步模块得到误差电压，反馈到NCO，最终完成载波锁定，实现QPSK相干解调；位同步模块原理框图如图2所示，由于输入I和Q的采样数据采样率比较低，且最佳采样点未知，故需要位同步模块鉴别位同步误差，根据位同步误差调整I和Q的插值位置，得到最佳采样点。

论文“高性能并行QPSK软解调技术研究”中对GPU实现QPSK解调做了初步尝试，给出了解调框图，但未给出适用于GPU运算的并行算法。论文“基于CPU-GPU异构结构的流式MPSK软件无线电信号解调方法及系统”基于CPU-GPU异构结构的流式MPSK软件无线电信号解调，采用CPU进行载波同步和位同步和解调，用GPU完成帧同步和TPC译码以及解扰运算，由于载波同步和位同步存在实时反馈运算，故不适合并行计算，GPU并行运算能力得不到高效的发挥。

发明内容

针对上述技术的缺陷，本发明提供一种基于帧结构的GPU并行计算QPSK相干解调方法，主要包括位同步和载波同步以及符号判决，一帧数据采用并行处理方案，位同步误差和载波同步误差一帧更新一次，一帧采用相同的位同步插值位置和载波同步相位补偿，解决环路锁定更稳定，抗频偏能力更强。

为实现上述目的，本发明提供如下技术方案：

一种基于帧结构的GPU并行计算QPSK相干解调方法，包括增益控制模块、位同步模块和载波同步模块，AD数据经过DDC正交下变频后变为符号率8倍的基带IQ数据，DDC模块将数据分段为固定包长的一包数据传送到GPU解调模块，增益控制模块对一包数据做增益控制，增益控制后传输给位同步模块以及载波同步模块。

一种基于帧结构的GPU并行计算QPSK相干解调方法，包括如下步骤：

步骤1，位同步的数据拼接模块对本包数据和上次未残留的末尾数据进行拼接，拼接长度固定为L0×8+1，L0为一帧码元数，8为一个码元过采样8次，初始拼接IQ数据均为0+0i，i为虚数单位，拼接后的一包数据长度变为Len2＝Len+L0×8+1，记拼接数据向量为IQ_last,其长度为L0×8+1，记拼接后向量为IQ_all,长度为Len2；

步骤2，取数模块根据取数位置k以及参数L0进行取数，k的初始值为L0×8+2；如果位置k≥Len+1,则说明剩余的数据不足一帧长度，则保存本包末尾的一帧数据以便和第二包数据拼接，同时更新下包数据取数位置k＝L0×8-(Len2-k)+1；如果k<Len+1,则按照如下规则取数：

IQ2sym(2n-1)＝IQ_all(k+4(n-1))

IQ2sym(2n)＝IQ_all(k+4(n-1)+1),n＝1,2,3...2L0+1

其中，IQ2sym(2n-1)代表取数后的奇数位，IQ2sym(2n)表示取数的偶数位，IQ_all表示拼接后的IQ数据；

步骤3，线性插值模块根据插值间隔deta对IQ2sym线性插值，其中deta初始值为0.5,，插值结果记为IQ2sym_i公式可表示为：

IQ2sym_i(n)＝IQ2sym(2n-1)+deta(IQ2sym(2n)-IQ2sym(2n-1)),n＝1,2,3...2L0+1

对上述IQ2sym_i(n)从第一个数开始二倍抽取，得到符号序列IQsym，IQsym作为位同步输出，送入载波同步模块。IQsym公式表示为：

IQsym(n)＝IQ2sym_i(2n-1),n＝1,2,3,...,L0

步骤4，位同步误差模块：对上一步的插值结果利用gardner算法计算位同步误差e_t:

上式中Is表示插值结果IQ2sym_i的实部，Qs表示IQ2sym_i的虚部；

步骤5，二阶环路滤波器：对上式计算的误差进行环路滤波，更新积分部分cum_i：

cum_i＝cum_i+pid_i*e_t

积分加上比例：

deta＝deta+cum_i+pid_p*e_t

上式cum_i为二阶环路滤波的积分部分，初始值为0，deta为环路滤波器的输出，初始值为0.5，pid_i为滤波器的积分部分常数系数，pid_p为比例系数；更新整数位置k和插值间隔deta模块：本模块对上述deta做如下判断，如果deta>1,则deta＝deta-1，k＝k+L0*8+1；如果deta<0,则deta＝deta+1，k＝k+L0*8-1；如果0<deta<1,则deta不变，k＝k+L0*8；

步骤6，载波同步复数乘法模块：对位同步输出IQsym，复乘补偿载波，得到补偿后的IQ数据记为IQcom，公式表示为：

IQcom(n)＝IQsym(n)×IQc(n),n＝1,2,3,...,L0

IQc(n)为复数载波产生模块的输出，初始值均为0+0i。不难理解IQcom即为最终的解调输出；极性costas环鉴频和鉴相模块：本模块对IQcom向量计算当前频率和相位偏移方向，公式表示为：

相位符号输出：

记点积：dot(n)＝Ic(n)Ic(n+1)+Qc(n)Qc(n+1),n＝1,2,3,...,L0-1

记互积：cross(n)＝Ic(n)Qc(n+1)-Qc(n)Ic(n+1),n＝1,2,3,...,L0-1

则频率符号输出：

上式中Qc表示IQcom虚部，Ic表示IQcom实部，sign(x)表示符号函数，x≥0，则sign(x)＝1,否则sign(x)＝0；

步骤7，补偿相位计算模块：对补偿相位进行更新,

若pd＝0,fd＝0，

则ct_cum＝ct_cum+1.5*cos_i，ct_out＝ct_cum+1.5*cos_p

若pd＝0,fd＝1，

则ct_cum＝ct_cum-cos_i，ct_out＝ct_cum-cos_p

若pd＝1,fd＝0，

则ct_cum＝ct_cum+cos_i，ct_out＝ct_cum+cos_p

若pd＝1,fd＝1，

则ct_cum＝ct_cum-1.5*cos_i，ct_out＝ct_cum-1.5*cos_p

其中cos_i为积分因子，cos_p为比例因子，ct_out为补偿相位输出；

步骤8，复数载波产生模块：

其中，i是虚数单位，φ₀为初始相位，初始条件为φ₀＝0，更新φ₀：φ₀＝φ₀+L0·ct_out。

与现有技术相比，本发明提供的一种基于帧结构的GPU并行计算QPSK相干解调方法有益效果如下：

1、本发明提供一种基于帧结构的GPU并行计算QPSK相干解调方法，主要包括位同步和载波同步以及符号判决，一帧数据采用并行处理方案，位同步误差和载波同步误差一帧更新一次，一帧采用相同的位同步插值位置和载波同步相位补偿，解决环路锁定更稳定，抗频偏能力更强。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为背景技术中的载波同步原理框图；

图2为背景技术中的位同步模块原理框图；

图3为本发明提供一种基于帧结构的GPU并行计算QPSK相干解调方法的框图；

图4为本发明提供一种基于帧结构的GPU并行计算QPSK相干解调方法的位同步模块和载波同步模块框图；

图5为本发明提供一种基于帧结构的GPU并行计算QPSK相干解调方法的数据分帧与拼接框图。

具体实施方式

下面将通过具体实施方式对本发明的技术方案进行清楚、完整地描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图3-5所示，本实施例提供一种基于帧结构的GPU并行计算QPSK相干解调方法，包括增益控制模块、位同步模块和载波同步模块，AD数据经过DDC正交下变频后变为符号率8倍的基带IQ数据，DDC模块将数据分段为固定包长的一包数据传送到GPU解调模块，增益控制模块对一包数据做增益控制，增益控制后传输给位同步模块以及载波同步模块。

本申请的并行计算主要体现在以下三个方面：1、增益控制模块：对整包IQ数据做幅度均值归一化。2、位同步模块：对一帧数据为单位进行线性插值，插值后对插值数据采用gardner算法进行位同步误差鉴别，插值算法和gardner算法架构固定，采用GPU并行计算实现。3、载波同步模块：载波同步算法主要包括三个子模块，复数向量相乘；极性costas环鉴频和鉴相；补偿相位计算；其中复数向量相乘和costas环鉴频和鉴相的计算以一帧数据为单位采用并行计算。

如图3所示，该方法包括增益控制模块、位同步模块和载波同步模块。AD数据经过DDC正交下变频后变为符号率8倍的基带IQ数据。DDC模块将数据分段为固定包长的一包数据传送到GPU解调模块。增益控制模块对一包数据做增益控制，增益控制后传输给位同步模块以及载波同步模块。位同步模块和载波同步模块详细框图如图4所示，详细步骤如下：

步骤1，位同步的数据拼接模块对本包数据和上次未残留的末尾数据进行拼接，拼接长度固定为L0×8+1，L0为一帧码元数，8为一个码元过采样8次，初始拼接IQ数据均为0+0i，i为虚数单位，拼接后的一包数据长度变为Len2＝Len+L0×8+1，记拼接数据向量为IQ_last,其长度为L0×8+1，记拼接后向量为IQ_all,长度为Len2；

步骤2，取数模块根据取数位置k以及参数L0进行取数，其示意图如下图5所示，k的初始值为L0×8+2；如果位置k≥Len+1,则说明剩余的数据不足一帧长度，则保存本包末尾的一帧数据以便和第二包数据拼接，同时更新下包数据取数位置k＝L0×8-(Len2-k)+1；如果k<Len+1,则按照如下规则取数：

IQ2sym(2n-1)＝IQ_all(k+4(n-1))

IQ2sym(2n)＝IQ_all(k+4(n-1)+1),n＝1,2,3...2L0+1

IQ2sym(2n-1)代表取数后的奇数位，IQ2sym(2n)表示取数的偶数位，IQ_all表示拼接后的IQ数据；

步骤3，线性插值模块根据插值间隔deta对IQ2sym线性插值，其中deta初始值为0.5。插值结果记为IQ2sym_i公式可表示为：

IQ2sym_i(n)＝IQ2sym(2n-1)+deta(IQ2sym(2n)-IQ2sym(2n-1)),n＝1,2,3...2L0+1

对上述IQ2sym_i(n)从第一个数开始二倍抽取，得到符号序列IQsym，IQsym作为位同步输出，送入载波同步模块。IQsym公式表示为：

IQsym(n)＝IQ2sym_i(2n-1),n＝1,2,3,...,L0

步骤4，位同步误差模块：对上一步的插值结果利用gardner算法计算位同步误差e_t:

上式中Is表示插值结果IQ2sym_i的实部，Qs表示IQ2sym_i的虚部。

步骤5，二阶环路滤波器：对上式计算的误差进行环路滤波，更新积分部分cum_i：

cum_i＝cum_i+pid_i*e_t

积分加上比例：

deta＝deta+cum_i+pid_p*e_t

上式cum_i为二阶环路滤波的积分部分，初始值为0，deta为环路滤波器的输出，初始值为0.5，pid_i为滤波器的积分部分常数系数，pid_p为比例系数，更新整数位置k和插值间隔deta模块：本模块对上述deta做如下判断，如果deta>1,则deta＝deta-1，k＝k+L0*8+1；如果deta<0,则deta＝deta+1，k＝k+L0*8-1；如果0<deta<1,则deta不变，k＝k+L0*8；

步骤6，载波同步复数乘法模块：对位同步输出IQsym，复乘补偿载波，得到补偿后的IQ数据记为IQcom，公式表示为：

IQcom(n)＝IQsym(n)×IQc(n),n＝1,2,3,...,L0

IQc(n)为复数载波产生模块的输出，初始值均为0+0i。不难理解IQcom即为最终的解调输出；极性costas环鉴频和鉴相模块：本模块对IQcom向量计算当前频率和相位偏移方向，公式表示为：

相位符号输出：

记点积：dot(n)＝Ic(n)Ic(n+1)+Qc(n)Qc(n+1),n＝1,2,3,...,L0-1

记互积：cross(n)＝Ic(n)Qc(n+1)-Qc(n)Ic(n+1),n＝1,2,3,...,L0-1

则频率符号输出：

上式中Qc表示IQcom虚部，Ic表示IQcom实部。sign(x)表示符号函数，x≥0，则sign(x)＝1,否则sign(x)＝0。

步骤7，补偿相位计算模块：对补偿相位进行更新，

若pd＝0,fd＝0，

则ct_cum＝ct_cum+1.5*cos_i，ct_out＝ct_cum+1.5*cos_p

若pd＝0,fd＝1，

则ct_cum＝ct_cum-cos_i，ct_out＝ct_cum-cos_p

若pd＝1,fd＝0，

则ct_cum＝ct_cum+cos_i，ct_out＝ct_cum+cos_p

若pd＝1,fd＝1，

则ct_cum＝ct_cum-1.5*cos_i，ct_out＝ct_cum-1.5*cos_p

其中cos_i为积分因子，cos_p为比例因子，ct_out为补偿相位输出。

步骤8，复数载波产生模块：

其中i是虚数单位。φ₀为初始相位，初始条件为φ₀＝0，更新φ₀：φ₀＝φ₀+L0·ct_out。

本实施例提供的该QPSK相干解调方法，详细步骤如下：

1.对本包数据和上次未残留的末尾数据进行拼接，拼接长度固定为L0×8+1，L0为一帧码元数，8为一个码元过采样8次。初始拼接IQ数据均为0。拼接后的一包数据长度变为Len2＝Len+L0×8+1，记拼接数据向量为IQ_last,其长度为L0×8+1，记拼接后向量为IQ_all,长度为Len2。拼接示意图如图5所示。

2.取数模块根据取数位置k以及参数L0进行取数，其示意图如图5所示，k的初始值为L0×8+2；如果k≥Len+1,则保存本包末尾的一帧数据以便和第二包数据拼接，同时更新下包数据取数位置k＝L0×8-(Len2-k)+1；如果k<Len+1,则按照如下规则取数：

IQ2sym(2n-1)＝IQ_all(k+4(n-1))

IQ2sym(2n)＝IQ_all(k+4(n-1)+1),n＝1,2,3...2L0+1

IQ2sym(2n-1)代表取数后的奇数位，IQ2sym(2n)表示取数的偶数位，IQ_all表示拼接后的IQ数据。

3.根据插值间隔deta对IQ2sym线性插值，其中deta初始值为0.5。插值结果记为IQ2sym_i表示为：

IQ2sym_i(n)＝IQ2sym(2n-1)+deta(IQ2sym(2n)-IQ2sym(2n-1)),n＝1,2,3...2L0+1，

4.对上述IQ2sym_i(n)从第一个数开始二倍抽取，得到符号序列IQsym：

IQsym(n)＝IQ2sym_i(2n-1),n＝1,2,3,...,L0

5.上一步的插值结果利用gardner算法计算位同步误差，同步误差e_t:

6.上式中Is表示插值结果IQ2sym_i的实部，Qs表示IQ2sym_i的虚部，对上式计算的误差进行环路滤波，更新积分部分cum_i：

cum_i＝cum_i+pid_i*e_t

deta＝deta+cum_i+pid_p*e_t

上式cum_i为二阶环路滤波的积分部分，初始值为0，deta为环路滤波器的输出，初始值为0.5，pid_i为滤波器的积分部分常数系数，pid_p为比例系数。

7.如果deta>1,则deta＝deta-1，k＝k+L0*8+1；如果deta<0,则deta＝deta+1，k＝k+L0*8-1；如果0<deta<1,则deta不变，k＝k+L0*8；

8.对步骤4的输出IQsym，做载波补偿载波，得到补偿后的IQ数据记为

IQcom：IQcom(n)＝IQsym(n)×IQc(n),n＝1,2,3,...,L0

IQc(n)为复数载波向量，初始值均为0+0i。不难理解IQcom即为最终的解调输出。

9.对IQcom向量计算当前频率和相位偏移方向，公式表示为：

相位输出

记点积：dot(n)＝Ic(n)Ic(n+1)+Qc(n)Qc(n+1),n＝1,2,3,...,L0-1

记互积：cross(n)＝Ic(n)Qc(n+1)-Qc(n)Ic(n+1),n＝1,2,3,...,L0-1

频率输出

上式中Qc表示IQcom虚部，Ic表示IQcom实部。sign(x)表示符号函数，x≥0，则sign(x)＝1,否则sign(x)＝0。

10.对补偿相位进行更新，

若pd＝0,fd＝0，

则ct_cum＝ct_cum+1.5*cos_i，ct_out＝ct_cum+1.5*cos_p

若pd＝0,fd＝1，

则ct_cum＝ct_cum-cos_i，ct_out＝ct_cum-cos_p

若pd＝1,fd＝0，

则ct_cum＝ct_cum+cos_i，ct_out＝ct_cum+cos_p

若pd＝1,fd＝1，

则ct_cum＝ct_cum-1.5*cos_i，ct_out＝ct_cum-1.5*cos_p

其中cos_i为积分因子，cos_p为比例因子，ct_out为补偿相位输出。

11.复数载波产生：

其中i是虚数单位。φ₀为初始相位，初始条件为φ₀＝0，对φ₀更新：

φ₀＝φ₀+L0·ct_out。

以上所述的实施例仅仅是对本发明的优选实施方式进行描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案做出的各种变形和改进，均应落入本发明装置权利要求书确定的保护范围内。

Claims (2)

1.一种基于帧结构的GPU并行计算QPSK相干解调方法，其特征在于，包括增益控制模块、位同步模块和载波同步模块，AD数据经过DDC正交下变频后变为符号率8倍的基带IQ数据，DDC模块将数据分段为固定包长的一包数据传送到GPU解调模块，增益控制模块对一包数据做增益控制，增益控制后传输给位同步模块以及载波同步模块。

2.根据权利要求1所述一种基于帧结构的GPU并行计算QPSK相干解调方法，其特征在于，包括如下步骤：

步骤1，位同步的数据拼接模块对本包数据和上次未残留的末尾数据进行拼接，拼接长度固定为L0×8+1，L0为一帧码元数，8为一个码元过采样8次，初始拼接IQ数据均为0+0i，i为虚数单位，拼接后的一包数据长度变为Len2＝Len+L0×8+1，记拼接数据向量为IQ_last,其长度为L0×8+1，记拼接后向量为IQ_all,长度为Len2；

步骤2，取数模块根据取数位置k以及参数L0进行取数，k的初始值为L0×8+2；如果位置k≥Len+1,则说明剩余的数据不足一帧长度，则保存本包末尾的一帧数据以便和第二包数据拼接，同时更新下包数据取数位置k＝L0×8-(Len2-k)+1；如果k<Len+1,则按照如下规则取数：

IQ2sym(2n-1)＝IQ_all(k+4(n-1))

IQ2sym(2n)＝IQ_all(k+4(n-1)+1),n＝1,2,3...2L0+1

其中，IQ2sym(2n-1)代表取数后的奇数位，IQ2sym(2n)表示取数的偶数位，IQ_all表示拼接后的IQ数据；

步骤3，线性插值模块根据插值间隔deta对IQ2sym线性插值，其中deta初始值为0.5,，插值结果记为IQ2sym_i公式可表示为：

IQ2sym_i(n)＝IQ2sym(2n-1)+deta(IQ2sym(2n)-IQ2sym(2n-1)),n＝1,2,3...2L0+1

对上述IQ2sym_i(n)从第一个数开始二倍抽取，得到符号序列IQsym，IQsym作为位同步输出，送入载波同步模块。IQsym公式表示为：

IQsym(n)＝IQ2sym_i(2n-1),n＝1,2,3,...,L0

步骤4，位同步误差模块：对上一步的插值结果利用gardner算法计算位同步误差e_t:

上式中Is表示插值结果IQ2sym_i的实部，Qs表示IQ2sym_i的虚部；

步骤5，二阶环路滤波器：对上式计算的误差进行环路滤波，更新积分部分cum_i：

cum_i＝cum_i+pid_i*e_t

积分加上比例：

deta＝deta+cum_i+pid_p*e_t

上式cum_i为二阶环路滤波的积分部分，初始值为0，deta为环路滤波器的输出，初始值为0.5，pid_i为滤波器的积分部分常数系数，pid_p为比例系数；更新整数位置k和插值间隔deta模块：本模块对上述deta做如下判断，如果deta>1,则deta＝deta-1，k＝k+L0*8+1；如果deta<0,则deta＝deta+1，k＝k+L0*8-1；如果0<deta<1,则deta不变，k＝k+L0*8；

步骤6，载波同步复数乘法模块：对位同步输出IQsym，复乘补偿载波，得到补偿后的IQ数据记为IQcom，公式表示为：

IQcom(n)＝IQsym(n)×IQc(n),n＝1,2,3,...,L0

IQc(n)为复数载波产生模块的输出，初始值均为0+0i。不难理解IQcom即为最终的解调输出；极性costas环鉴频和鉴相模块：本模块对IQcom向量计算当前频率和相位偏移方向，公式表示为：

相位符号输出：

记点积：dot(n)＝Ic(n)Ic(n+1)+Qc(n)Qc(n+1),n＝1,2,3,...,L0-1

记互积：cross(n)＝Ic(n)Qc(n+1)-Qc(n)Ic(n+1),n＝1,2,3,...,L0-1

则频率符号输出：

上式中Qc表示IQcom虚部，Ic表示IQcom实部，sign(x)表示符号函数，x≥0，则sign(x)＝1,否则sign(x)＝0；

步骤7，补偿相位计算模块：对补偿相位进行更新,

若pd＝0,fd＝0，

则ct_cum＝ct_cum+1.5*cos_i，ct_out＝ct_cum+1.5*cos_p

若pd＝0,fd＝1，

则ct_cum＝ct_cum-cos_i，ct_out＝ct_cum-cos_p

若pd＝1,fd＝0，

则ct_cum＝ct_cum+cos_i，ct_out＝ct_cum+cos_p

若pd＝1,fd＝1，

则ct_cum＝ct_cum-1.5*cos_i，ct_out＝ct_cum-1.5*cos_p

其中cos_i为积分因子，cos_p为比例因子，ct_out为补偿相位输出；

步骤8，复数载波产生模块：

其中，i是虚数单位，φ₀为初始相位，初始条件为φ₀＝0，更新φ₀：φ₀＝φ₀+L0·ct_out。

Images