[发明专利]基于纠错纠删RS‑Turbo级联码的跳频抗干扰方法

摘要

本发明公开一种基于纠错纠删RS‑Turbo级联码的跳频抗干扰方法，主要解决现有跳频通信系统在载波频偏和部分频带干扰下无法正常通信的问题。其实现步骤是(1)对一帧数据进行RS编码；(2)将RS码各个码元分配到相应的跳上；(3)分别对每跳数据进行Turbo编码；(4)设置跳频组帧格式，并进行QPSK调制；(5)对调制信号加入载波频偏、相偏和高斯白噪声，得到复基带信号；(6)对复基带信号进行载波同步；(7)对同步信号进行Turbo译码，给出RS码的删除位；(8)根据RS码的删除位进行纠错纠删译码。本发明具有误比特性能好，抗干扰性能强的优点，可用于移动通信、卫星通信、深空通信及遥测。

主权项

一种基于纠错纠删RS‑Turbo级联码的跳频抗干扰方法，包括如下步骤：(1)对一帧数据分组进行RS编码：(1a)将长度为N比特的一帧数据分为B组，每组含N/B比特数据；(1b)分别对每组数据进行RS(n,k0,d)编码，得到B个RS码字，其中n为码字长度，k0为信息码元长度，每个码元包含的比特数为m，d为最小码距，d＝2t+1，t为纠错个数，为向下取整运算；(2)将每个RS码字的各个码元分配到跳频系统相应的跳上：(2a)将B个RS码字并行排放成B列，每列含n个码元，每行含mB个编码比特，再将B列RS码字划分为B0块，每块n行，每行含mB/B0个编码比特；(2b)将第1块第1行B/B0个码元即mB/B0个编码比特分配到跳频系统的第1跳，依次类推，第1块第n行B/B0个码元分配到跳频系统的第n跳，同理，将第B0块第1行B/B0个码元分配到跳频系统的第[n(B0‑1)+1]跳，依次类推，第B0块第n行B/B0个码元分配到跳频系统的第nB0跳；(3)分别对每跳数据进行Turbo编码：(3a)分别在每跳mB/B0个编码比特末尾添加M个冲洗比特，得到(mB/B0+M)比特数据；(3b)分别对每跳(mB/B0+M)比特数据进行码率为1/3的Turbo编码，最后得到长为Da的Turbo编码比特，其中Da＝3(mB/B0+M)；(4)设置跳频组帧格式：(4a)将长度为P比特的导频序列作为跳频系统的导频跳；(4b)设第1,2,…,nB0跳的数据格式相同，其每一跳均将长为Le比特的引导码置于长为Da比特的Turbo编码数据之前，组成长为Da+Le比特的数据，其中Da+Le＝P；(4c)将导频跳置于第1,2,…,n跳之前组成跳频系统的第1帧，依次类推，将导频跳置于第[n(B0‑1)+1],[n(B0‑1)+2],…,nB0跳之前组成跳频系统的第B0帧，一帧数据帧长Nb＝(Le+Da)n；(5)对每帧数据进行正交相移键控QPSK调制，得到发送的复基带信号sk；(6)对复基带信号sk加入载波频偏、载波相偏和相位抖动，通过高斯白噪声信道，在接收端将接收到的信号经过滤波、下变频、数字采样后得到复基带信号rk：rk=sk·ej(2πkΔfT+θ+Δθ)+nk,k=0,1,...,Nb2-1,]]>其中，Δf、θ和Δθ分别为载波频偏、载波相偏和相位抖动，T为符号周期，ΔfT为归一化频偏，|ΔfT|<1，nk是均值为零方差为σ2的复高斯随机变量，N0为噪声的单边功率谱密度，Nb为一帧数据帧长；(7)载波同步：(7a)从复基带信号rk中提取导频序列dk，利用导频序列由旋转平均周期图RPA算法对载波偏差进行粗估计，得到频偏粗估计值fRPA，利用该粗估计值fRPA对复基带信号rk进行频偏校正，得到频偏校正信号rk′；(7b)从频偏校正信号rk′提取引导码，利用引导码由最大似然算法得到相位初始估计值θ0；(7c)从频偏校正信号rk′中提取数据序列xk，由科斯塔斯环法迭代计算相偏校正序列yk：yk=xke-jθk,k=0,1,...,Da2-1,]]>其中，θk为第k个数据码元的相位校正量，当k＝0时，θ0为相位初始估计值，Da为引导码的长度；(8)将相位校正序列yk送入Turbo译码器，采用基于对数的最大后验概率LOG‑MAP算法进行迭代译码，输出所有RS码元每个比特的软值；(9)纠错纠删RS译码：(9a)设输入至纠错纠删RS译码器的码字C＝(c1,c2,…,ch,…,cn)，其中，ch为第h个码元，h＝1,2,…,n，设码元ch的软值度量为βh＝(βh1,βh2,…,βhj,…,βhm)，其中，βhj为码元ch第j个比特的软值，j＝1,2,…,m，找出软值度量βh中绝对值的最小值γh，并将γh作为码元ch的可靠性度量值，则码字C的可靠性度量γ＝(γ1,γ2,…,γh,…,γn)；(9b)对码字C的可靠性度量γ从小到大排序，排列结果为(γp1,γp2,…,γph,…,γpn)，ph∈1,2,…,n，其中，ph为排列在第h位的可靠性度量值对应码元的下标，相应的码元排列结果为(cp1,cp2,…,cph,…,cpn)；(9c)设置码字C的删除个数为e，即cp1,cp2,…,cpe位为删除位；(9d)对码字C进行纠错纠删译码，如果错误数在码字C的纠删能力范围内，则将译码码字作为输出，译码结束，否则，执行步骤(9e)；(9e)设置计数i＝1、码字存储矩阵U为空、搜索区间长度为s，在e+1,e+2,…,e+s中不重复的选出q个数，共有种组合，用vi1,vi2,…,viq，i∈1,2,…,w表示其中的一种组合，构造w行、q列的删除矩阵V，将vi1,vi2,…,viq，i∈1,2,…,w存储在删除矩阵V的第i行；(9f)在删除位cp1,cp2,…,cpe以外再选取删除矩阵V中第i行的q位作为删除位，即将位作为删除位，重新设置码字C的删除个数e′＝e+q；(9g)对码字C进行纠错纠删译码，如果错误数在码字C的纠删能力范围内，则译码正确，并将译码码字保存在码字存储矩阵U中，其中，h＝1,2,…,n，第h个码元为码元第j个比特的硬判决值，其中，(9h)将计数i加1后重新赋给i；(9i)如果赋值后的计数i小于等于组合种类w，则返回到步骤(9f)，否则，执行步骤(9j)；(9j)如果码字存储矩阵U为空则将接收的码字作为译码输出，译码结束，否则，分别将码字存储矩阵U中每个码字与码字C的软值按位进行相关运算，译码码字Di与码字C的软值按位进行相关运算计算公式为：Cori=Σh=1nΣj=1mβhjdhji;]]>(9k)选取最大相关值所对应的码字作为译码输出，译码结束。