高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于余弦符合度的自同步扰码盲识别

张立民 谭继远 钟兆根 吴昭军

张立民, 谭继远, 钟兆根, 吴昭军. 基于余弦符合度的自同步扰码盲识别[J]. 电子与信息学报, 2022, 44(4): 1412-1420. doi: 10.11999/JEIT210248
引用本文: 张立民, 谭继远, 钟兆根, 吴昭军. 基于余弦符合度的自同步扰码盲识别[J]. 电子与信息学报, 2022, 44(4): 1412-1420. doi: 10.11999/JEIT210248
Yang Guangzheng, Yang Xiangyu, Xu Lijuan. THEOREMS OF TSENG RUNS OF BINARY SEQUENCE[J]. Journal of Electronics & Information Technology, 1998, 20(4): 500-507.
Citation: ZHANG Limin, TAN Jiyuan, ZHONG Zhaogen, WU Zhaojun. Blind Recognition of Self-synchronous Scramblers Based on Cosine Conformity[J]. Journal of Electronics & Information Technology, 2022, 44(4): 1412-1420. doi: 10.11999/JEIT210248

基于余弦符合度的自同步扰码盲识别

doi: 10.11999/JEIT210248
基金项目: 国家自然科学基金(91538201);泰山学者工程专项(ts201511020);信息系统安全技术重点实验室基金(6142111190404)
详细信息
    作者简介:

    张立民:男,1966年生,教授,博士生导师,研究方向为卫星信号处理及应用

    谭继远:男,1994年生,博士生,研究方向为信道编码识别

    钟兆根:男,1984年生,博士,副教授,研究方向为扩频信号处理

    吴昭军:男,1992年生,博士,助理研究员,研究方向为信道编码识别

    通讯作者:

    谭继远 tjiyuan@outlook.com

  • 中图分类号: TN911.7

Blind Recognition of Self-synchronous Scramblers Based on Cosine Conformity

Funds: The National Natural Science Foundation of China (91538201), Taishan Scholar Special Foundation (ts201511020), The Chinese National Key Laboratory of Science and Technology on Information System Security (6142111190404)
  • 摘要: 为克服现有非合作自同步扰码识别算法在低信噪比下识别率低、适应性差的缺点,该文提出一种基于余弦符合度的自同步扰码盲识别方法。首先基于信源不平衡性和自同步扰码解扰原理建立自同步扰码的含错校验方程,然后将接收到的软判决序列转化为信息码元的后验概率序列,遍历可能的生成多项式,在遍历的过程中引入余弦符合度作为统计量,通过分析符合度的统计特性求解出最优的判别门限,根据统计量和判别门限的关系完成自同步扰码生成多项式的识别。仿真结果表明:该算法能有效识别出生成多项式,且在低信噪比下识别率优于现有算法,具有较好的低信噪比适应能力,在信源不平衡度ε为0.1,截获扰码序列长度为800 bit和ε为0.05,截获扰码序列长度为3000 bit时,能够有效完成生成多项式的识别,与目前算法相比,该文算法识别性能优于现有的硬判决算法,且比硬判决算法性能提升1~2 dB。
  • 在通信领域中,为了改善信号的连0或连1特性,提高信息传输的随机性,会对发送信号进行加扰处理,加扰后的信息序列0和1趋于平衡,通过加扰不仅便于提取定时信息和功率谱的平坦化,而且可以增强信息传输的安全性[1,2]。然而,对于非合作领域,侦察方需要从截获的扰码数据中得到相关信息,此时便需要对扰码数据进行解扰得到原始信源序列,因此扰码数据的参数盲识别技术对后续的信号分析处理有着重要的意义和实用价值[3]

    扰码可分为自同步扰码和同步扰码两种。两者的构造基础都为线性反馈移位寄存器(Linear Feedback Shift Register, LFSR)。同步扰码的盲识别包括生成多项式和扰码初态的识别,近年来得到了广泛的研究[4-7]。自同步扰码的盲识别是根据截获的扰码数据进行生成多项式的确定。自同步扰码由于其保密性强,且需要将扰码序列作为LFSR的输入,因此其生成多项式的识别更加复杂。目前针对自同步扰码盲识别的研究主要集中为两大类:基于求解方程或基于统计量。文献[8]通过构造解扰序列的随机变量,当正确解扰和错误解扰时,随机变量服从不同的分布,从而设立判别门限,识别出生成多项式。但是判别门限的设立需要知道信源不平衡度,算法实用性不强。文献[9]通过对扰码数据的重码进行统计,根据重码统计的极性分布与阶数关系判定多项式阶数,然后根据组合不平衡性判定多项式抽头位置,从而完成生成多项式的识别。该方法阶数估计时所需数据量会随着阶数呈指数增加,计算复杂量较大。文献[10]在已知信源不平衡度的情况下,提出了一种基于游程统计的自同步扰码多项式阶数估计方法,通过信源不平衡度和游程统计极值的关系完成多项式阶数的估计。文献[11]根据解扰后的比特状态不平衡性,以统计概率分布和均匀分布间的修正平方欧几里得距离为准则,完成了自同步扰码的识别。文献[12,13]分别提出了基于卷积码和RS码的自同步扰码盲识别。以上方法都是采用硬判决加扰序列进行盲识别的,但是通常截获的扰码序列都为调制信息和信道噪声的软信息序列,因此以上算法普遍存在低信噪比下识别率低的问题。软判决信息由于含有可靠度信息,近几年在编码识别领域的应用取得了较好效果[14-16]。文献[17]采用了硬判决和软判决相结合的方法,采用软判决度量过程中,算法采用了简单的近似处理,其算法性能在低信噪比条件下仍然具有较大的损失。同时未设置相应的判别门限。从现有算法来看,低信噪比下自同步扰码的识别性能还有待进一步提高。

    针对现有算法在低信噪比下识别率较低的缺点,本文提出一种低信噪比下自同步扰码的盲识别算法。该算法首先根据扰码的信源不平衡性和解扰原理,构造满足扰码约束关系的校验方程。然后将接收到的软判决序列转化为信息码元的后验概率序列,遍历可能的生成多项式。在校验过程中,引入余弦符合度,该符合度能较好地反映校验方程成立的概率。利用错误解扰时的随机序列和正确解扰时的不平衡序列在平均余弦符合度下的统计特性差异,设立相应的判决门限,完成对生成多项式的识别。与以往方法相比,提出的算法在低信噪比下的识别性能得到了较大的提升。

    自同步扰码的加扰和解扰过程的基础单元为LFSR,通过将信源序列xk和反馈逻辑的输出sk进行模二求和,得到扰码序列yk。加扰过程和LFSR反馈结构如图1所示。

    图 1  自同步扰码加扰器结构

    加扰过程可以表示为

    yk=xkLi=1ciyki
    (1)

    其中,表示模2累加。ciGF(2),i=1,2,,L为LFSR的反馈系数,GF(2)表示2元域,ci的取值为0或1,c1=cL=1。扰码的生成多项式可以表示为

    f(x)=1+c1x+c2x2++cLxL
    (2)

    扰码的解扰过程如图2所示。

    图 2  自同步扰码解扰

    解扰过程表达式为

    xk=ykLi=1ciyki=Li=0ciyki
    (3)

    实际的通信系统信息传输中,截获的扰码序列是加扰序列经过调制和信道传输后的信息序列,在调制过程中会进行相应的映射,信道传输过程中会不可避免地混入噪声。因此本文研究的自同步扰码的模型如图3所示。

    图 3  自同步扰码模型

    模型表达式可以表示为

    rk=zk+nk
    (4)

    其中,zkyk调制后的调制信息序列,nk为信道噪声。有偏信源中P(xk=0)=1/2+εε通常取值为0.05或0.1。本文假设接收到的扰码序列采用BPSK调制,信号幅度为A,噪声为均值为0,方差为σ2的高斯白噪声。信噪比定义为

    SNR=10lg(A22σ2)
    (5)

    对于接收方来说,截获的信息序列为rk,需要获取的信息序列为xk,此时,需要已知生成多项式系数ci以完成发送信源序列的恢复。自同步扰码的识别问题即为根据接收到的信息序列rk识别出生成多项式ci的过程。在实际通信中,自同步扰码使用的LFSR阶数绝大多数为3~100阶,并且其生成多项式多为稀疏多项式,自同步扰码由于误码扩散率和生成多项式的项数成正比,因此一般都是用2项式或3项式[11]。本文的算法也建立在此基础上。

    由式(3)可知,正确解扰的情况下,当输入信息为0时,有表达式(6)成立

    ykLi=1ciyki=Li=0ciyki=0
    (6)

    输入信息为1时,有表达式(7)成立

    ykLi=1ciyki=Li=0ciyki=1
    (7)

    综合式(6)和式(7),解扰过程可以转化为如式(8)所示的方程组

    [000y1y1y20y1y1y2yLyL+1y2y3yL+2yNsLyNsL+1yNs][cLcL1c11]=[x1x2xNs]
    (8)

    其中,Ns为截获的扰码序列长度,L为生成多项式阶数。

    此时根据式(8)构造校验方程为

    [000y1y1y20y1y1y2yLyL+1y2y3yL+2yNsLyNsL+1yNs][cLcL1c11]=[000]
    (9)

    当遍历的多项式fk(x)为正确的生成多项式时,解扰序列就是加扰前的序列,因此解扰序列x1,x2,,xNs0,1不平衡序列,此时校验方程组中方程成立的个数和方程不成立的个数不均衡;当遍历的多项式fk(x)不是正确的生成多项式时,假设自同步扰码多项式为fm(x)=xr+xm+1,m=1,2,,r1,,,mk,则有式(10)成立

    xj=yjyjkyjr=(ˆxjyjmyjr)yjkyjr=ˆxjyjmyjk
    (10)

    其中,ˆxj为加扰前的序列。

    由加扰序列的白化性可知[9]P((yjmyjk)=1)=1/2,进而推导得到

    P(xj=1)=P((ˆxjyjmyjk)=1)1/2(1p)+(11/2)p=1/2
    (11)

    因此,解扰序列x1,x2,,xNs0,1趋于平衡的随机序列,此时校验方程组中方程成立的个数和不成立的个数近似相等。

    传统的自同步扰码的识别方法主要采用硬判决扰码序列或者软硬判决相结合的方式。这些方法忽略了软判决中的可靠度信息,在复杂的信道环境中性能较差。为了尽可能利用来自信道的调制和噪声等信息,本文采用直接截获的软判决序列,完成生成多项式的识别。为了方便说明,将校验方程的第k行单独列出讨论,即

    L+1l=1yk,lcLl=0
    (12)

    其中,1kNs

    为了衡量方程成立的可能性大小,建立余弦符合度为

    F(k)=Ll=0cos(πP(yk,l|rk,l)cLl)
    (13)

    其中,yk,0,yk,1,,yk,l为硬判决序列,对应的获取到的软判决序列为rk,0,rk,1,,rk,l,记P(yk,l|rk,l)为截获扰码软判决序列为rk,l的条件下,硬判决序列yk,l取1的概率。

    P(yk,l|rk,l)cLl<0.5时,yk,lcLl0的概率大,此时余弦符合度为正值;当P(yk,l|rk,l)cLl>0.5时,yk,lcLl取1的概率大,此时余弦符合度为负值。因此当校验方程成立时,式(12)等于0,式(12)中模2加等于1的个数为偶数个,此时对应的F(k)为正值;当校验方程不成立时,式(12)等于1,式(12)中模2加等于1的个数为奇数个,此时对应的F(k)为负值。式(13)中后验概率P(yk,l|rk,l)的推导如下。

    假设截获的加扰软判决为rk,l,对应的硬判决扰码为yk,lyk,l的条件概率密度可以表示为

    P(rk,l|yk,l=0)=12πσe(rk,l+A)2/(2σ2)
    (14)
    P(rk,l|yk,l=1)=12πσe(rk,lA)2/(2σ2)
    (15)

    根据贝叶斯公式可知

    P(yk,l|rk,l)=P(rk,l|yk,l=1)P(yk,l=1)P(rk,l)
    (16)

    其中,P(rk,l)可以表示为

    P(rk,l)=P(rk,l|yk,l=1)P(yk,l=1)+P(rk,l|yk,l=0)P(yk,l=0)
    (17)

    由于没有先验信息,此时令P(yk,l=0)=P(yk,l=1)=0.5,由式(14)—式(17)可以得到

    P(yk,l|rk,l)=e2Ark,l/σ2e2Ark,l/σ2+1
    (18)

    将式(18)代入式(13)可以得到余弦符合度的表达式为

    F(k)=Ll=0cos(πe2Ark,lσ2e2Ark,lσ2+1cLl)
    (19)

    由式(19)可知,当校验方程式(12)成立时,对应的解扰信息为0,0与SNR成正比,SNR越大,F(k)越趋近于1;当校验方程不成立时,对应的解扰信息为0,F(k)与SNR成反比,SNR越大,F(k)越趋近于–1。此时考虑式(9)中所有的校验方程,对余弦符合度取统计平均,得到平均余弦符合度为

    ¯F=1NsNsk=1F(k)
    (20)

    当遍历的多项式为生成多项式时,解扰序列的0,1比例不均衡,此时F(k),k=0,1,,Ns中取正值和取负值的比例不均衡,解扰序列中0多于1时,¯F为正值;解扰序列中1多于0时,¯F为负值,且信噪比越大,¯F越趋近于1或–1。当遍历的多项式不是扰码的生成多项式时,解扰序列的0,1比例均衡,为随机序列。此时F(k),k=0,1,,Ns中取正值和取负值的比例均衡,¯F值在0附近徘徊。以此为依据可以判定遍历的多项式是否为扰码的生成多项式。

    对于平均余弦符合度的统计特性,需要引入相应的判别门限以完成生成多项式的识别。

    G=[cLcL1c11],设其中等于1的抽头个数为m,对应的抽头位置集合为G1=[cLcnc11],此时参与校验方程的加扰硬判决序列为:[yk,1  yk,Ln+1    yk,L  yk,L+1]

    当校验方程成立时,参与校验方程的加扰硬判决序列[yk,1yk,Ln+1yk,Lyk,L+1]中码元1模2加等于0,此时对应的码元1的个数为偶数个,此时所有可能的组成情况可表示为

    Vk,0=m/2j=0C2jm
    (21)

    其中, 表示向下取整,Cnm=n(n1)(nm+1)m(m1)1

    根据均值和方差的定义,对每种可能情况的均值和方差取统计平均,结果为

    μk,0=m/2j=0C2jmVk,0(cos(πe2Ax/σ2e2Ax/σ2+1)p(r|y=1)dx)2j(cos(πe2Ax/σ2e2Ax/σ2+1)p(r|y=0)dx)m2j
    (22)
    σ2k,0=m/2j=0C2jmVk,0((cos(πe2Ax/σ2e2Ax/σ2+1))2p(r|y=1)dx)2j((cos(πe2Ax/σ2e2Ax/σ2+1))2p(r|y=0)dx)m2jμ2k,0
    (23)

    当校验方程不成立时,参与校验方程的加扰硬判决序列[yk,1  yk,Ln+1    yk,L  yk,L+1]中码元1模2加等于1,此时对应的码元1的个数为奇数个,所有可能的组成情况可表示为

    Vk,1=m/2j=0C2j1m
    (24)

    根据均值和方差的定义,对每种可能情况的均值和方差取统计平均,结果为

    μk,1=m/2j=0C2j1mVk,1(cos(πe2Ax/σ2e2Ax/σ2+1)p(r|y=1)dx)2j1(cos(πe2Ax/σ2e2Ax/σ2+1)p(r|y=0)dx)m2j+1
    (25)
    σ2k,1=m/2j=0C2j1mVk,1((cos(πe2Ax/σ2e2Ax/σ2+1))2p(r|y=1)dx)2j1((cos(πe2Ax/σ2e2Ax/σ2+1))2p(r|y=0)dx)m2j+1μ2k,1
    (26)

    当遍历的多项式是扰码生成多项式时,校验方程成立的个数l0和不成立的个数l1符合信源不平衡度,即

    l0l1=0.5+ε0.5ε
    (27)

    当遍历的多项式不为生成多项式时,参与校验方程的加扰硬判决序列[yk,1  yk,Ln+1    yk,L  yk,L+1]中码元1的个数奇偶随机,此时所有组成情况可表示为

    V1=mj=0Cjm
    (28)

    根据均值和方差的定义,对每种可能情况的均值和方差取统计平均,结果为

    μ1=mj=0CjmV(cos(πe2Ax/σ2e2Ax/σ2+1)p(r|y=1)dx)j(p(r|y=0)dx)mj
    (29)
    σ21=mj=0CjmV1((cos(πe2Ax/σ2e2Ax/σ2+1))2p(r|y=1)dx)j((cos(πe2Ax/σ2e2Ax/σ2+1))2p(r|y=0)dx)mjμ21
    (30)

    式(21)—式(30)中的积分表达式不存在解析解,为了完成快速运算和达到较高的计算精度,本文采用数值积分方式进行求解。

    首先给出以下两种假设检验:

    H0:遍历的多项式不是自同步扰码生成多项式;

    H1:遍历的多项式是自同步扰码生成多项式。

    根据大数定律可知,当截获的扰码序列长度Ns较大,遍历的多项式不是自同步扰码生成多项式时,¯F服从均值为μ2,方差为σ22=σ21/Ns的高斯分布,在非生成多项式的情况下,均值μ的上下浮动会降低生成多项式的检测性能,因此,适当提升H0假设中的均值,仍能保证检测性能[18],但是过度提升均值会导致算法所需数据量的增加,本算法令μ2=1.5σ2/Ns,此时大部分(92%)的测试值在此范围内。

    H0:¯FN(μ2,σ22)
    (31)

    同理,当遍历的多项式是自同步扰码生成多项式时,¯F服从均值为μ3,方差为σ23=σ2/Ns的高斯分布,即

    H1:¯FN(μ3,σ23)
    (32)

    其中,μ3=(0.5+ε)μk,0+(0.5ε)μk,1, σ23=(0.5+ε)2σ2k,0+(0.5ε)2σ2k,1

    设判别门限为Λ,虚警概率Pf和漏警概率Pn分别为

    Pf=Λ12πσ2e(yμ2)2/σ22dx=2[1Φ(Λ|μ2|σ2)]
    (33)
    Pn=Λ12πσ3e(yμ3)2/σ23dx=Φ(Λ|μ3|σ3)
    (34)

    在实际运用中,首先需要确定虚警概率Pf,然后根据式(33)得到判别门限为

    Λ=σ2(Φ1(1Pf/2)+1.5)
    (35)

    ¯F>Λ时,假设H1成立,判定遍历的多项式为生成多项式;反之,假设H0成立,判定遍历的多项式不为生成多项式。

    基于余弦符合度的自同步扰码盲识别算法步骤如下:

    步骤1 将截获的扰码软判决序列根据式(16)转化为后验概率序列;

    步骤2 构造并储存2抽头和3抽头的3~100阶多项式;

    步骤3 确定虚警概率Pf

    步骤4 遍历步骤2中的多项式,同时根据式(9)构造校验矩阵;

    步骤5 根据式(19)和式(20)计算平均余弦符合度¯F

    步骤6 根据式(35),计算判别门限Λ。若¯FΛ,则识别出生成多项式;否则跳转到步骤4,直到¯FΛ

    本节主要进行了如下的仿真验证:仿真1:给定虚警概率时,在固定信源下,得出平均余弦符合度和判别门限分布图;仿真2:截获扰码序列长度对识别性能的影响;仿真3:生成多项式阶数对识别性能的影响;仿真4:信源不平衡度对识别性能的影响;仿真5:与其他算法的对比。在不做特殊说明的情况下,以下仿真验证默认给定虚警概率Pf=103,生成多项式为ITUV.34中使用的生成多项式:f(x)=1+x18+x23,蒙特卡罗仿真次数为1000次。

    SNR=6dB,在信源不平衡度ε=0.1,截获扰码数据Ns=1000ε=0.05,Ns=3000两种情况下,遍历所有多项式后,统计量和遍历多项式以及判别门限分布图分别如图4(a)图4(b)所示。

    图 4  统计量和判别门限关系图

    图4可以看出,遍历的多项式正确时,统计量有一个明显的突出谱线,且遍历正确时,平均余弦符合度谱线明显高于判别门限;遍历错误时,平均余弦符合度谱线低于判别门限。因此本文所提算法能有效识别生成多项式。

    仿真设定的截获的扰码序列信源不平衡度分别为0.1和0.05。在ε=0.1时设定截获的扰码序列长度分别为400, 600, 800, 1000;在ε=0.05时设定截获的扰码序列长度分别为1500, 2000, 2500, 3000;信噪比范围为–1~10 dB,间隔1取值。统计在不同的信噪比情况下,截获扰码序列长度与生成多项式正确识别率的关系,结果如图5(a)图5(b)所示。

    图 5  截获扰码序列长度对算法识别率的影响

    图5可以看出,通过增加截获的扰码序列长度可以明显提高算法的识别率。当信源不平衡度较低时,可以通过增加扰码序列长度来克服由于ε下降造成的算法识别性能下降。从仿真结果可以看出,算法在低信噪比下识别性能较好。在ε=0.1,截获扰码数据为800,以及ε=0.05,截获扰码数据为3000,信噪比为2dB的情况下,生成多项式的识别率能达到90%以上。

    在不同阶数生成多项式情况下验证算法的性能。分别采用8阶生成多项式:f(x)=1+x2+x8,23阶生成多项式:f(x)=1+x18+x23和31阶生成多项式:f(x)=1+x11+x31。信源不平衡度ε=0.1,生成的扰码序列长度为600时,信噪比范围为–1~10 dB,间隔1取值。统计在不同信噪比情况下,生成多项式阶数与算法正确识别率之间的关系,结果如图6所示。

    图 6  生成多项式阶数对算法识别率的影响

    图6中可以看出,相同截获序列长度和ε下,同一信噪比时,生成多项式的阶数越低,算法的识别率越高。且三者的识别率均随信噪比的提升而升高,阶数越低,算法识别率越早达到100%。在信噪比为4 dB时,算法的正确识别率均能达到90%以上,具有较好的多项式阶数容错性能。

    为了验证信源不平衡度对算法性能的影响。采用23阶生成多项式:f(x)=1+x18+x23,序列长度2500,信源不平衡度取值0.04, 0.06, 0.08。信噪比范围为–1~10 dB,间隔1取值。统计在不同信噪比下信源不平衡度与算法正确识别率的关系,结果如图7所示。

    图 7  信源不平衡度对算法识别率的影响

    图7中可以看出,在相同信噪比下,随着ε的增大,算法的正确识别率得到了提高,信源越不平衡,算法的识别率越高。随着信噪比的增大,ε越大,算法识别率越快达到100%。

    与本文算法进行对比的是基于比特不均衡算法[11],Walsh软判决算法[17]和Cluzeau算法[8]。其中Cluzeau算法需要知道信源不平衡度,仿真该算法时ε已知。假设首先对这4种算法的容错性能进行对比,设定生成多项式:f(x)=1+x18+x23,截获序列长度分别为800,ε=0.1和3000, ε=0.05。信噪比范围为–1~10 dB,间隔1取值。通过蒙特卡罗仿真实验,统计在不同信噪比下4种算法的正确识别率,结果如图8(a)图8(b)所示。

    图 8  4种算法的容错性能对比

    图8中可以看出,本文所提算法在低信噪比约有1~2 dB的信噪比增益。主要原因是本文较好利用软判决序列构造符合度,极大利用了软序列中的信息可靠度。在低信噪比条件下性能提升较为明显。

    其次,对这4种算法的复杂度进行了比较。假设截获的扰码序列长度为N,扰码多项式的项数为d,扰码多项式阶数为L,遍历的多项式个数为S。本文算法的计算复杂度主要集中在余弦符合度的计算上,第1次遍历多项式的过程中,需要进行d次余弦运算,d次乘法运算和N次加法运算。为方便分析,将1次余弦运算等效为3次乘法运算,因此遍历1次多项式需要进行4d次乘法运算和N次加法运算。考虑到最差的情况,遍历到最后一个多项式才识别出生成多项式,则最大计算量为4Sd次乘法运算和SN次加法运算。比特不均衡算法比特状态间隔为L,一次遍历需要进行NL次比特操作,故需要(NL)S次比特运算。Walsh软判决算法计算复杂度主要集中在Hadamard变换部分,共需要2L+1(L+1)S次乘法运算和SN次加法运算;Cluzeau算法1次遍历时遍历部分的计算量为1/ε2次乘法和NL次模二和运算(等效为加法运算),总计算量为S/ε2乘法运算和S(NL)次加法运算。ε越小,文献[8]的算法复杂度越大。文献[18]的算法复杂度随着扰码阶数L的增加呈指数增加,算法复杂度远大于本文算法。虽然本文算法复杂度较文献[12]略有提升,但是本文算法复杂度在可承受范围内,同时识别性能和算法容错性得到了较大的改善。

    本文首先利用自同步扰码解扰原理和信源的不平衡性,构建出了符合自同步扰码约束关系的校验方程,然后引入平均余弦符合度,将校验方程成立的可能性以概率的形式表示,当遍历正确的多项式和遍历错误的多项式时,平均余弦符合度服从不同的分布;最后通过分析平均余弦符合度的统计特性设立最优的判别门限,完成生成多项式的识别。本文算法在低信噪比下的适应性能较好,具有较好的识别性能和算法容错性,同时计算复杂度低于传统的Walsh和Cluzeau等算法,在非合作通信领域具有较强的工程实用性。

  • 图  1  自同步扰码加扰器结构

    图  2  自同步扰码解扰

    图  3  自同步扰码模型

    图  4  统计量和判别门限关系图

    图  5  截获扰码序列长度对算法识别率的影响

    图  6  生成多项式阶数对算法识别率的影响

    图  7  信源不平衡度对算法识别率的影响

    图  8  4种算法的容错性能对比

  • [1] 马钰, 张立民. 基于实时检测的扰码重建算法[J]. 电子与信息学报, 2016, 38(007): 1794–1799.

    MA Yu, ZHANG Limin. Scrambling code reconstruction algorithm based on real-time detection[J]. Journal of Electronics and Information Technology, 2016, 38(007): 1794–1799.
    [2] KIM D, SONG J, and YOON D. On the estimation of synchronous scramblers in direct sequence spread spectrum systems[J]. IEEE Access, 2020, 8: 166450–166459. doi: 10.1109/ACCESS.2020.3023425
    [3] 邱钊洋, 李天昀, 查雄. 非周期长码直扩信号同步及伪码序列盲估计[J]. 电子与信息学报, 2021, 43(8): 2171–2180. doi: 10.11999/JEIT201096

    QIU Zhaoyang, LI Tianyun, and CHA Xiong. Aperiodic Long Code DSSS Signal Synchronization and Blind Estimation of Pseudo Code Sequence[J]. Journal of Electronics and Information Technology, 2021, 43(8): 2171–2180. doi: 10.11999/JEIT201096
    [4] 钟兆根, 孙雪丽, 马钰. 线性扰码重建算法性能分析[J]. 系统工程与电子技术, 2019, 41(2): 205–210.

    ZHONG Zhaogen, SUN Xueli, and MA Yu. Performance analysis of linear scrambling code reconstruction algorithm[J]. Systems Engineering and Electronics, 2019, 41(2): 205–210.
    [5] 解辉, 韩壮志, 丁爽. 基于传播算子算法的扰码序列估计方法[J]. 系统工程与电子技术, 2017, 39(10): 2327–2332.

    XIE Hui, HAN Zhuangzhi, and DING Shuang. Scrambling sequence estimation method based on propagation operator algorithm[J]. Systems Engineering and Electronics, 2017, 39(10): 2327–2332.
    [6] HAN Shunan, ZHANG Min, and LI Xinhao. Reconstruction of Feedback Polynomial of Synchronous Scrambler Based on Triple Correlation Characteristics of M-sequences[J]. IEICE Transactions on Communications, 2018, ElO1.B(7): 1723–1732.
    [7] 韩树楠, 张旻, 李歆昊. 基于m序列三阶相关性的同步扰码反馈多项式重构[J]. 电子学报, 2019, 47(3): 552–559.

    HAN Shunan, ZHANG Min, and LI Xinhao. Synchronous scrambling feedback polynomial reconstruction based on the third-order correlation of m-sequences[J]. Acta Electronica Sinica, 2019, 47(3): 552–559.
    [8] CLUZEAU M. Reconstruction of a linear scrambler[J]. IEEE Transactions on Computers, 2007, 56(9): 1283–1291. doi: 10.1109/TC.2007.1055
    [9] 吕喜在, 苏绍璟, 黄芝平. 一种新的自同步扰码多项式盲恢复方法[J]. 兵工学报, 2011, 32(6): 680–685.

    LV Xizai, SU Shaojin, and HUANG Zhiping. A new blind recovery method of self-synchronizing scrambling polynomial[J]. Acta Armamentarius, 2011, 32(6): 680–685.
    [10] 黄芝平, 周靖, 苏绍璟, 等. 基于游程统计的自同步扰码多项式阶数估计[J]. 电子科技大学学报, 2013(4): 541–545.

    HUANG Zhiping, ZHOU Jing, SU Shaojing, et al. Self-synchronizing scrambling polynomial order estimation based on run-length statistics[J]. Journal of University of Electronic Science and Technology of China, 2013(4): 541–545.
    [11] 廖红舒, 袁叶, 甘露. 自同步扰码的盲识别方法[J]. 通信学报, 2013, 34(1): 136–143.

    LIAO Hongshu, YUAN Ye, and GAN Lu. Blind recognition method of self-synchronizing scrambling code[J]. Journal on Communications, 2013, 34(1): 136–143.
    [12] HAN S, and ZHANG M. A method for blind identification of a scrambler based on matrix analysis[J]. Communications Letters, IEEE, 2018, 22(11): 2198–2201. doi: 10.1109/LCOMM.2018.2868681
    [13] 尹瑾, 王建新. RS码的自同步扰码盲识别方法[J]. 计算机工程与应用, 2017, 53(22): 77–81.

    YIN Jin, and WANG Jianxin. Self-synchronizing scrambling code blind recognition method of RS code[J]. Computer Engineering and Applications, 2017, 53(22): 77–81.
    [14] WU Zhaojun, ZHANG Limin, and ZHONG Zhaogen. A maximum cosinoidal cost function method for parameter estimation of RSC Turbo codes[J]. IEEE Communications Letters, 2019, 23(3): 390–393. doi: 10.1109/LCOMM.2018.2890224
    [15] 吴昭军, 张立民, 钟兆根, 等. 一种软判决下的RS码识别算法[J]. 电子与信息学报, 2020, 42(9): 2150–2157. doi: 10.11999/JEIT190592

    WU Zhaojun, ZHANG Limin, ZHONG Zhaogen, et al. A RS code recognition algorithm under soft decision[J]. Journal of Electronics and Information, 2020, 42(9): 2150–2157. doi: 10.11999/JEIT190592
    [16] 欧世峰, 赵艳磊, 宋鹏,等. 基于概率耦合的双直接判决先验信噪比估计算法[J]. 电子学报, 2020, 48(8): 151–160.

    OU Shifeng, ZHAO Yanlei, SONG Peng, et al. Dual direct decision prior signal-to-noise ratio estimation algorithm based on probabilistic coupling[J]. Acta Electronica Sinica,, 2020, 48(8): 151–160.
    [17] 尹瑾, 王建新. 基于软判决求解含错方程的自同步扰码盲识别[J]. 探测与控制学报, 2017(2): 46–50, 55.

    YIN Jin and WANG Jianxin. Self-synchronizing scrambling code blind recognition based on soft decision to solve error equations[J]. Journal of Detection and Control, 2017(2): 46–50, 55.
    [18] MA Yu, ZHANG Limin, and WANG Haotong. Reconstructing synchronous scrambler with robust detection capability in the presence of noise[J]. IEEE Transactions on Information Forensics and Security, 2015, 10(2): 397–408. doi: 10.1109/TIFS.2014.2378143
  • 加载中
图(8)
计量
  • 文章访问数:  727
  • HTML全文浏览量:  446
  • PDF下载量:  60
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-26
  • 修回日期:  2021-06-29
  • 网络出版日期:  2021-07-08
  • 刊出日期:  2022-04-18

目录

/

返回文章
返回