Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
高级搜索

留言板

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

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

Twofish算法中密钥相关S-盒的差分性质分析及其改进

周旋 李超

张祖凡, 杨作为, 王国仲. 智能反射表面辅助太赫兹信道估计的低复杂度算法[J]. 电子与信息学报, 2023, 45(10): 3640-3647. doi: 10.11999/JEIT221476
引用本文: 周旋, 李超. Twofish算法中密钥相关S-盒的差分性质分析及其改进[J]. 电子与信息学报, 2004, 26(6): 912-916.
ZHANG Zufan, YANG Zuowei, WANG Guozhong. A Low Complexity Algorithm for Intelligent Reflective Surface-assisted Tera Hertz Channel Estimation[J]. Journal of Electronics & Information Technology, 2023, 45(10): 3640-3647. doi: 10.11999/JEIT221476
Citation: Zhou Xuan, Li Chao. Differential Analysis and Modification of the Key-Dependent S-Boxes of Twofish[J]. Journal of Electronics & Information Technology, 2004, 26(6): 912-916.

Twofish算法中密钥相关S-盒的差分性质分析及其改进

Differential Analysis and Modification of the Key-Dependent S-Boxes of Twofish

  • 摘要: 该文从理论上证明了Twofish算法中,密钥越长,密钥相关S-盒的差分概率就越小,提出了一种新的与密钥作用的方式来产生密钥相关S-盒的方法,理论与测试结果表明新的S-盒的异或差分概率和模加差分概率比原算法的差分概率要小.
  • 作为下一代移动通信候选关键技术,太赫兹(Tera Hertz, THz)通信和智能反射表面(Intelligent Reflecting Surface, IRS)成为专家学者研究的热点。简单来说,THz通信相比于毫米波(millimeter Wave, mmWave),THz通信可以提供更丰富的带宽(0.1~10 THz)、更高的数据速率(Tbit/s级)和更低的延迟(微秒级)[1,2]。但是,由于THz信号波长较短,导致信号衰减严重、衍射能力差[3],这些缺点严重限制了THz频段的应用,IRS能够有效地弥补这些缺点。IRS是一个由大量反射元件组成的平面阵列,每个反射元件可以控制入射信号的幅度和相位,通过调整和反射发送端的信号,从而达到重构信道环境的作用[4,5]。将IRS与THz通信相结合,可以增加THz信号的传输距离,减少信号阻塞以及提高通信的可靠性。因此,许多科研人员针对IRS在通信的应用进行了大量的研究工作。

    IRS相关基础理论研究方面,文献[6]讨论了IRS硬件实现的工作原理,并且讨论了在传统无线通信系统中引入IRS时的基本问题以及和传统的中继的比较。在通过部署IRS提高通信系统性能方面,文献[7]讨论了在传统多小区通信场景中,通过在多小区的边界部署IRS,协助下行小区边缘用户通信,减轻小区间干扰。在IRS辅助多用户通信部分信道未知的场景下,文献[8]提出一种混合波束赋形算法,提高系统中所有用户的和速率。在物理层安全通信方面,为了减少窃听者对合法接收者的信息接收影响和窃听者的有用信息接收量,文献[9]通过IRS实现在无线供电通信网络中窃听端的波束抵消。

    面对IRS辅助THz通信的场景中,复杂的传播环境给信道估计带来了挑战,如何充分发挥THz巨大的通信潜力,以及如何通过IRS解决THz通信中的信道不稳定的缺点成为讨论的重点。文献[10]提出了一种基于压缩感知的低复杂度的信道估计算法,通过消除迭代过程中感知矩阵的冗余列,减少了计算复杂度。同样地,文献[11]利用信道传输过程中的稀疏性,解决了级联信道的估计问题,很好地解决了信道估计过程中复杂度过高和训练开销过大的缺点。文献[12]提出一种基于发送端和IRS协同的波束训练方案来进行信道估计,根据得到的信道信息进一步进行混合波束赋形,减少了射频链路的消耗,并且提出的方案能够在完美信道状态信息的条件下达到接近全数字波束赋形的性能。虽然上述研究成果可以解决无线通信系统中的信道估计问题,但是,随着THz通信中天线数目的增加、IRS的大量部署以及IRS元件数的增加,接收信号的维度也会大大增加,从而导致解决信道估计问题的成本增加。

    因此,本文针对上述缺点,提出一种低复杂度的信道估计算法。在分析THz信道特点的基础上对IRS进行分组设计,并且将IRS辅助THz通信的信道模型表示为统一数学表达式,有效避免在估计IRS级联信道中由IRS元件数增加带来开销过高的问题。通过仿真和复杂度分析表明本文提出的算法满足归一化均方误差(Normalized Mean Square Error, NMSE)指标的前提下进一步降低信道估计开销。

    考虑一个THz通信系统,假设直视路径受阻,通过引入IRS缓解通信过程中的多径衰落和阻塞。发送端发射的信号经过环境中散射体的散射或者经过IRS反射后到达接收端。IRS辅助THz通信系统模型如图1所示,基站配备NtBS根天线,用户端配备NrUE根天线,IRS具有N个反射阵元。其中,基站和用户分别配备NtRFNuRF个射频链路。则用户接收到的信号可以表示为

    图 1  系统模型
    y=WHHx+n (1)

    其中,W=WRF×WB表示用户端的预编码矩阵,WRFCNrUE×NrRFWBCNrRF×r分别表示用户端数字和模拟预编码矩阵;HCNrUE×NtBS表示信道矩阵;xCNtBS×1表示经过发送端混合波束赋形后的信号,其中x=FRFFBs, FRFCNtBS×NtRFFBCNtRF×r分别表示基站端数字和模拟预编码矩阵,sCr×1为基站发送的数据流;nCr×1表示加性高斯白噪声。

    为了不失一般性,用散射路径有限的几何信道模型表示THz信道,根据文献[13]的描述,THz信号经过多次反射后的能量衰减超过20 dB,因此,本文只考虑1次反射情况。根据几何信道模型,每个散射路径决定1个单一的传播路径[14],可以用式(2)来描述

    Hs=L1l1=1gl1(f,d)aUE(φl1)aHBS(θl1) (2)

    其中,L1表示路径数量;gl1(f,d)表示路径损耗。在THz频段,本文主要考虑信号的自由空间损耗和分子吸收损耗,满足

    gl1(f,d)=c4πfde12τ(f)d (3)

    其中,c表示光速,τ(f)表示分子吸收损耗,详细的参数可以参考文献[15]。aBS(θ), aUE(φ)分别表示基站和用户的阵列响应向量。假设基站和用户采用均匀线性阵列,则阵列响应可以表示为

    aBS(θ)=1NtBS[1,ejkdtsin(θ),,ejkdt(NtBS1)sin(θ)]T (4)
    aUE(φ)=1NrUE[1,ejkdrsin(φ),,ejkdr(NrUE1)sin(φ)]T (5)

    其中,k=2π/λλ表示波长,dtdr表示收发端天线的间距,θφ表示收发端角度。

    在如图1所示的系统模型中,信号经过IRS反射的级联信道可以描述为

    HI=HtΦHr (6)

    其中,HtHr分别表示基站-智能反射表面(Base Station-Intelligent Reflecting Surface, BS-IRS)和智能反射表面-用户终端(Intelligent Reflecting Surface-User Equipment, IRS-UE)的信道矩阵;Φ表示IRS的相移矩阵,Φ=diag([β1ejθ1,β2ejθ2,,βnejθn]),幅度βn{0,1},相位θn(0,2π],简单起见,假设βn=1。本文将图1中完整的IRS划分为M个子阵列,每个子阵列具有Mi=N/M个反射元件,简单起见,假设每个子阵列为线性阵列,则级联信道可以进一步表示为

    HI=Mm=1gBIm(f,d)aUE(φ)aHIRS,m(ψ)ΦmgImU(f,d)aIRS,m(ψ)aHBS(θ) (7)

    gm(f,d)=gBI(f,d)gIU(f,d), am=aHIRS,m(ψ)ΦmaIRS,m(ψ)。与之前相关工作文献中的信道不同的是,本文将IRS视作受控散射体,根据信号在空间中的传播路径,重新建立一个信道模型。则式(7)可以进一步表示为

    HI=Mm=1gmamaUE(φm)aHBS(θm) (8)

    可以清楚地看到式(8)和式(2)完成了形式上的统一。同时,也可以很清晰地看出IRS在环境中的作用为增加BS到UE的路径数量,并且通过控制IRS元件的振幅和相位,能够有效地调整每一条路径的增益,并且增益是已知可控的,从而实现控制信道环境的作用。完整的信道可以表述为

    H=HI+Hs=Kk=1αkaUE(φk)aHBS(θk) (9)

    其中,K=L1+M表示整个信道环境中的总路径数。从式(9)可以很清楚地看出IRS在无线通信中的作用为增加信号传播路径,即使在传播环境较为复杂的情况下也能提供较为稳定的传播环境。将信道模型表示成统一的数学表达式可以使得在信道估计的过程中,不需要对不同链路的信道进行单独估计。

    CP分解最早由Hitchcock提出[16],其主要思想是将张量分解为一系列低秩张量外积之和。假设XRI×J×K是一个3阶张量,根据CP分解可以表示成

    XNn=1anbncn (10)

    其中,表示向量的外积;aRI×n, bRJ×n, cRK×n。将所有的向量a组合成因子矩阵A=[a1,a2,,an],同理可得到因子矩阵BC。通过因子矩阵可以将张量表示成如式(11)的矩阵形式[17]

    X(1)A(CB)TX(2)B(CA)TX(3)C(BA)T} (11)

    其中,表示矩阵的Khatri–Rao积;X(n)表示张量的模式-n(n=1,2,3)展开形式。综上所述,张量的CP分解有如式(12)的表达形式

    X[[A,B,C]]Nn=1anbncn (12)

    对于下行估计信道,本文考虑一个信道训练过程。简单来说,将信道相干时间Ts分成T个块,每个块有L个时隙,从而使得Ts=TL,在训练过程中,假设发送信号在每个块重复发送,并且IRS在一个块内的相移保持不变。根据信道训练过程将用户在第t个块l个时隙内接收到的信号表示为

    y[t,l]=Kk=1gkαt,l,kWHaUE(φk)aHBS(θk)Fx[t,l]+n[t,l] (13)

    将式(13)表示为第t块总的接收信号

    Y[t] = Kk=1akaUE(φk)aHBS(θk)+N[t] (14)

    其中,Y[t]=[y[t,1],y[t,2],,y[t,l]],  aUE(φk)WHaUE(φk), N[t]=[n[t,1],n[t,2],,n[t,l]], akgkαt,k,  aBS(θk)(Fx[l])HaBS(θk)。式(14)满足文献[17]中CP分解的要求,则进一步将式(11)写成CP分解形式,即秩一张量叠加和的形式

    Y=Kk=1aUE(ϕk)aBS(θk)ak(t)+N=[[A,B,C]]+N (15)

    其中,YCr×L×T,并且

    A=[aUE(φ1),aUE(φ2),,aUE(φk)]Cr×KB=[aBS(θ1),aBS(θ2),,aBS(θk)]CL×KC=[a1,a2,,ak]CT×K} (16)

    由于THz信号具有稀疏性,即路径数L1很小,并且IRS子阵的数量M是人为可控的,在操作时,可以将总路径数控制在一定的范围内,接收信号张量Y具有低秩结构,这样的低秩结构可以保证张量的分解具有唯一性。

    假设THz信道的路径数量和IRS反射元件个数是已知的先验信息,CP分解的目标为

    min (17)

    其中,\hat {\boldsymbol{A}} = \left[ {{{\hat a}_1},{{\hat a}_2},\cdots,{{\hat a}_k}} \right], \hat {\boldsymbol{B}} = \left[ {{{\hat b}_1},{{\hat b}_2},\cdots,{{\hat b}_k}} \right], \hat {\boldsymbol{C}} = \left[ {{{\hat c}_1},{{\hat c}_2},\cdots,{{\hat c}_k}} \right]。上述优化问题可以通过交替最小二乘算法解决。在该算法中,通过固定两个变量矩阵求解另外一个变量矩阵,重复进行该过程,直到满足收敛准则。在CP分解唯一性的条件下,将张量 {\boldsymbol{\mathcal{Y}}} 展开成模式-n( n{\text{ = 1,2,3}} )的形式

    \left. \begin{gathered} {{\boldsymbol{Y}}_{(1)}} \approx {\boldsymbol{A}}{({\boldsymbol{C}} \odot {\boldsymbol{B}})^{\text{T}}} + {{\boldsymbol{N}}_{(1)}} \\ {{\boldsymbol{Y}}_{(2)}} \approx {\boldsymbol{B}}{({\boldsymbol{C}} \odot {\boldsymbol{A}})^{\text{T}}} + {{\boldsymbol{N}}_{(2)}} \\ {{\boldsymbol{Y}}_{(3)}} \approx {\boldsymbol{C}}{({\boldsymbol{B}} \odot {\boldsymbol{A}})^{\text{T}}} + {{\boldsymbol{N}}_{(3)}} \\ \end{gathered} \right\} (18)

    根据上述展开公式,采用交替最小二乘法求解得到

    {\hat {\boldsymbol{A}}^{(t + 1)}} = {\rm{arg}}\mathop {\min}\limits_A \left\| {{\boldsymbol{Y}}_{(1)}^{\text{T}} - \left( {{{\hat {\boldsymbol{C}}}^{(t)}} \odot {{\hat {\boldsymbol{B}}}^{(t)}}} \right){{\hat {\boldsymbol{A}}}^{\text{T}}}} \right\| (19)
    {\hat {\boldsymbol{B}}^{(t + 1)}} = {\rm{arg}}\mathop {\min}\limits_B \left\| {{\boldsymbol{Y}}_{(2)}^{\text{T}} - \left( {{{\hat {\boldsymbol{C}}}^{(t)}} \odot {{\hat {\boldsymbol{A}}}^{(t + 1)}}} \right){{\hat {\boldsymbol{B}}}^{\text{T}}}} \right\| (20)
    {\hat {\boldsymbol{C}}^{(t + 1)}} = {\rm{arg}}\mathop {\min}\limits_C \left\| {{\boldsymbol{Y}}_{(3)}^{\text{T}} - \left( {{{\hat {\boldsymbol{B}}}^{(t + 1)}} \odot {{\hat {\boldsymbol{A}}}^{(t + 1)}}} \right){{\hat {\boldsymbol{C}}}^{\text{T}}}} \right\| (21)

    然后对张量进行重构,得到{\hat {{\boldsymbol{\mathcal{Y}}}}^{(t + 1)}} = \left[\kern-0.15em\left[ {{\hat {\boldsymbol{A}}}^{(t + 1)}}, {{\hat {\boldsymbol{B}}}^{(t + 1)}},{{\hat {\boldsymbol{C}}}^{(t + 1)}} \right]\kern-0.15em\right]。当条件\left\| {{{\hat {{\boldsymbol{\mathcal{Y}}}}}^{(t + 1)}} - {{\hat {{\boldsymbol{\mathcal{Y}}}}}^{(t)}}} \right\|_{\text{F}}^2\Bigr/ \left\| {{{\hat {{\boldsymbol{\mathcal{Y}}}}}^{(t)}}} \right\|_{\text{F}}^2 \le \varepsilon满足时,认为收敛。接下来根据估计的因子矩阵对信道参数进行估计。

    通过上节对因子矩阵 {\boldsymbol{A}} , {\boldsymbol{B}} , {\boldsymbol{C}} 进行估计后,本节将会对具体的信道参数 {\varphi _k} , {\theta _k} , {g_k} 进行估计。在张量分解满足唯一性的条件下,分解得到的因子矩阵相较于原因子矩阵具有一定的误差,他们之间的关系为

    {{\hat {\boldsymbol{A}}}} = {\boldsymbol{A}}{{\boldsymbol{\varLambda}} _1}{\boldsymbol{\varPi}} + {{\boldsymbol{E}}_1} (22)
    {{\hat {\boldsymbol{B}}}} = {\boldsymbol{B}}{{\boldsymbol{\varLambda}} _2}{\boldsymbol{\varPi}} + {{\boldsymbol{E}}_2} (23)
    {{\hat {\boldsymbol{C}}}} = {\boldsymbol{C}}{{\boldsymbol{\varLambda}} _3}{\boldsymbol{\varPi}} + {{\boldsymbol{E}}_3} (24)

    其中, \{ {{\boldsymbol{\varLambda}} _1},{{\boldsymbol{\varLambda}} _2},{{\boldsymbol{\varLambda}} _3}\} 为未知的非奇异对角矩阵,满足 {{\boldsymbol{\varLambda}} _1}{{\boldsymbol{\varLambda}} _2}{{\boldsymbol{\varLambda}} _3} = {\boldsymbol{I}} {\boldsymbol{\pi}} 为未知的置换矩阵,在估计的过程中, {\boldsymbol{\pi}} 对所有信道参数的作用是一致的,可以不用对其进行估计; \{ {{\boldsymbol{E}}_1},{{\boldsymbol{E}}_2},{{\boldsymbol{E}}_3}\} 对应3个因子矩阵的估计误差。估计因子矩阵{{\hat {\boldsymbol{A}}}} {{\hat {\boldsymbol{B}}}} 的每一列包含收发端天线的角度信息 {\varphi _k} {\theta _k} ,可以通过一个基于相关性的估计器进行估计

    {\hat \varphi _l} = \mathop {\arg \max }\limits_{{\varphi _l}} \frac{{\left| {{\boldsymbol{\hat a}}_l^{\text{H}}{{{\boldsymbol{\tilde a}}}_{{\text{UE}}}}\left( {{\varphi _j}} \right)} \right|}}{{{{\left\| {{{{\boldsymbol{\hat a}}}_l}} \right\|}_2}{{\left\| {{{{\boldsymbol{\tilde a}}}_{{\text{UE}}}}\left( {{\varphi _j}} \right)} \right\|}_2}}},1 \le j \le J (25)
    {\hat \theta _l} = \mathop {\arg \max }\limits_{{\theta _l}} \frac{{\left| {{\boldsymbol{\hat b}}_l^{\text{H}}{{{\boldsymbol{\tilde a}}}_{{\text{BS}}}}\left( {{\theta _j}} \right)} \right|}}{{{{\left\| {{{{\boldsymbol{\hat b}}}_l}} \right\|}_2}{{\left\| {{{{\boldsymbol{\tilde a}}}_{{\text{BS}}}}\left( {{\theta _j}} \right)} \right\|}_2}}},1 \le j \le J (26)

    其中, {{\boldsymbol{\hat a}}_l} {{\boldsymbol{\hat b}}_l} 分别表示因子矩阵 {{\hat {\boldsymbol{A}}}} {{\hat {\boldsymbol{B}}}} 的第 l 列; J 表示搜索维度。

    根据求解的角度信息 {\varphi _k} {\theta _k} 重构出原始的阵列响应向量矩阵 {{\boldsymbol{\tilde a}}_{{\text{UE}}}}\left( \varphi \right) {{\boldsymbol{\tilde a}}_{{\text{BS}}}}\left( \theta \right) ,进一步重构出 {{\tilde {\boldsymbol{A}}}} = [{{{\tilde {\boldsymbol{a}}}}_1},{{{\tilde {\boldsymbol{a}}}}_2},\cdots,{{{\tilde {\boldsymbol{a}}}}_k}] {{\tilde {\boldsymbol{B}}}} = [{{{\tilde {\boldsymbol{b}}}}_1},{{{\tilde {\boldsymbol{b}}}}_2},\cdots,{{{\tilde{\boldsymbol{ b}}}}_k}] ,然后根据式(22)和式(23)得到 {{\boldsymbol{\varLambda}} _1} = {{{\tilde {\boldsymbol{A}}}}^\dagger }{{\hat {\boldsymbol{A}}}} {{\boldsymbol{\varLambda}} _2} = {{{\tilde {\boldsymbol{B}}}}^\dagger }{{\hat {\boldsymbol{B}}}} ,再利用 {{\boldsymbol{\varLambda}} _1}{{\boldsymbol{\varLambda}} _2}{{\boldsymbol{\varLambda}} _3} = {\boldsymbol{I}} 求出 {{\boldsymbol{\varLambda}} _3} ,接着将 {{\boldsymbol{\varLambda}} _3} 代入式(24)中求得 {{\tilde {\boldsymbol{C}}}} = [{{{\tilde {\boldsymbol{c}}}}_1},{{{\tilde {\boldsymbol{c}}}}_2},\cdots,{{{\tilde {\boldsymbol{c}}}}_k}] ,最后根据已知的智能反射表面矩阵 {\alpha _{t,k}} 、式(14)中的定义和式(16),使用式(27)可以求出信道路径增益

    [{\tilde g_1},{\tilde g_2},\cdots,{\tilde g_k}] = {\text{diag}}\left( {{{[{{\boldsymbol{a}}_1},{{\boldsymbol{a}}_2},\cdots,{{\boldsymbol{a}}_k}]}^\dagger }{{\tilde {\boldsymbol{C}}}}} \right) (27)

    至此,信道参数 {\varphi _k} , {\theta _k} , {g_k} 估计完成,代入式(9)中便可完成对信道的估计。

    本节讨论张量具有唯一分解需要满足的条件。根据文献[18],需要满足Kruskal’s条件,则平行因子分解具有唯一性,即满足

    {\text{rank}}({\boldsymbol{A}}) + {\text{rank}}({\boldsymbol{B}}) + {\text{rank}}({\boldsymbol{C}}) \ge 2K + 2 (28)

    首先讨论矩阵 {\boldsymbol{A}} 的秩。因为

    {\boldsymbol{A}} = {{\boldsymbol{W}}^{\text{T}}}[{\boldsymbol{a}}({\varphi _1}),{\boldsymbol{a}}({\varphi _2}),\cdots,{\boldsymbol{a}}({\varphi _k})] \triangleq {{\boldsymbol{W}}^{\text{T}}}{{\boldsymbol{A}}_{\text{R}}} \in {\mathbb{C}^{r \times K}} (29)

    其中, {{\boldsymbol{A}}_{\text{R}}} \in {\mathbb{C}^{N_{{\text{UE}}}^{\text{r}} \times K}} 是线性阵列响应向量的组合,波束成形矩阵 {\boldsymbol{W}} 是从单位圆中随机选择,根据文献[19]的证明,可以得到

    {\text{rank}}({\boldsymbol{A}}) = \min (r,K) (30)

    同理

    {\text{rank}}({\boldsymbol{B}}) = \min (L,K) (31)

    对于矩阵 {\boldsymbol{C}} 来说,{\boldsymbol{C}} = [{{\boldsymbol{a}}_1},{{\boldsymbol{a}}_2},\cdots,{{\boldsymbol{a}}_k}]{\text{diag}}({g_1}, {g_2},\cdots,{g_k}) \in {\mathbb{C}^{T \times K}}。根据文献[20],本文将IRS相移矩阵设计成离散傅里叶矩阵,则

    {\text{rank}}({\boldsymbol{C}}) = {\text{min}}(T,K) (32)

    通常来说,发送信号数目大于路径数,所以有 {\text{rank}}({\boldsymbol{A}}) = K 。为了满足式(28),则只需让{\text{rank}}({\boldsymbol{B}}) + {\text{rank}}({\boldsymbol{C}}) \ge K + 2,可以使得张量分解满足唯一性。

    实验仿真的参数设置如下:本文考虑收发端天线皆为线性阵列,天线数 {N_{\text{r}}} = {N_{\text{t}}} = 64 ,天线间距为信号波长的1/2。信道模型采用基于几何的THz信道模型,收发端角度随机分布在 \text{0}~\text{2}\pi ,THz信号频率 f = 0.14{\text{ THz}} ,分子吸收损耗[15] \tau (0.14T) = 1.83 \times {10^{ - 5}}/m , L = T = 16,IRS元件数设置为N{\text{=32}},将IRS分成多个子阵列并且对相移进行单独设计。收发两端的波束形成矩阵 {\boldsymbol{F}} 和组合矩阵 {\boldsymbol{W}} 从单位圆中随机选取,射频链路数 N_{{\text{RF}}}^{\text{r}} = N_{{\text{RF}}}^{\text{t}} = 16 。信噪比(SNR)定义为信号分量与噪声分量的比值。在本文仿真实验中,信噪比范围设置为(–10~20 dB)。本文采用均方误差(Mean Square Error, MSE)指标来统计信道参数的估计误差,具体的计算公式为

    \left. \begin{gathered} {\text{MSE}}({\boldsymbol{\varphi}} ) = \left\| {{\boldsymbol{\varphi}} - \hat {\boldsymbol{\varphi}} } \right\|_2^2 \\ {\text{MSE}}({\boldsymbol{\theta}} ) = \left\| {{\boldsymbol{\theta}} - \hat {\boldsymbol{\theta}} } \right\|_2^2 \\ {\text{MSE}}({\boldsymbol{g}}) = \left\| {{\boldsymbol{g}} - \hat {\boldsymbol{g}}} \right\|_2^2 \\ \end{gathered} \right\} (33)

    其中,{\boldsymbol{\varphi}} \triangleq \left[ {{\varphi _1},{\varphi _2},\cdots,{\varphi _k}} \right], {\boldsymbol{\theta }} \triangleq \left[ {{\theta _1},{\theta _2},\cdots,{\theta _k}} \right], {\boldsymbol{g}} \triangleq [{g_1},{g_2},\cdots,{g_k}]。并且采用NMSE指标统计信道的估计误差。

    基于正则平行因子分解的信道估计结果将会展示在本节中。首先讨论路径数 K 的大小对信道估计的影响。由于THz信号传播主要由直视路径控制,并且直视路径增益要高出散射路径15 dB左右,因此假设整个系统中不存在散射路径即 {L_1} = 0 ,当 N{\text{ = 36}} ,具体的仿真结果如图2所示。随着路径数的增加,信道估计的准确度有所下降,这是因为 K 的增大会导致信道矩阵的秩增加,从而导致估计的准确度降低。此外,随着信噪比的增加,由路径数增加导致估计性能降低的问题可以得到缓解。根据图2的结果,再结合实际信道情况,选取 K = 4 进行后续的仿真实验。

    图 2  不同路径数和信噪比下的信道估计的NMSE

    在收发端天线数均为64,IRS反射元件数32,并且将IRS分成4个子块, T = L = 16 ,在不同的信噪比下信道参数估计的均方误差如图3所示。从中可以看出,当信噪比增加时,估计误差能有效的降低。特别地,当信噪比范围为(–5~0 dB),信道参数估计会有一个急速下降的趋势,这表明本文提出的信道估计算法能够在信噪比稍差的场景中使用。

    图 3  信道参数估计的MSE

    图4比较了不同反射元件数下信噪比和导频开销的NMSE比较情况。可以看出随着信噪比的增大,信道的NMSE估计误差下降,当信噪比固定时,随着导频开销的增大,信道估计的准确度也会增加,并且当IRS反射元件数增加时,估计的性能会进一步降低,详细的结果可以参考图5。在信道条件较为恶劣的情况下,利用大量的导频训练,会改善估计的准确度,这为现实情况中当信号传输环境不足以满足信道估计条件时,可以采用增加导频开销的方法缓解信道估计准确度不足的问题。

    图 4  不同反射元件数下信噪比和导频开销的NMSE比较
    图 5  不同信噪比下反射元件和导频开销的NMSE比较

    图5比较了当 {\text{SNR}} = 10{\text{ dB}} {\text{SNR}} = 15{\text{ dB}} 时,IRS元件数变化时不同导频开销下信道估计的均方误差,信道估计的误差随着IRS元件数增多而降低,这是因为增加IRS元件数会增加传输信息的准确性。特别地,当导频开销一定时,随着IRS元件数的增加,信道估计的NMSE会趋于一个稳定值。通过观察图4图5可知,在信道估计的过程中,需要选取合适的导频开销和IRS反射元件数目。

    图6对比了文献[20]提出的最小二乘Khatri-Rao分解(Least Squares Khatri-Rao Factorization, LSKRF)、文献[21]提出的基于迭代的信道估计算法(Iterative Channel Estimation, ICE)和本文提出的基于正则平行因子分解的信道估计算法。具体的参数设置为 {N_{\text{r}}} = 16 , {N_{\text{t}}} = 32 , N = 16 , T = L = 32 , K = 4 。从图6可以看出本文提出的信道估计算法优于其他两种算法,具体来说,本文提出的算法在信噪比为20 dB时相比于文献[21]有4.28 dB左右的提升,相比于文献[20]有7.12 dB左右的提升。之所以出现这种差距,是因为在没有直视路径的条件下,IRS能够提供相对稳定的传输环境。

    图 6  3种不同算法的NMSE对比

    本文提出的信道估计算法的复杂度主要为求解式(19)—式(21)中因子矩阵的最小二乘问题。对于因子矩阵 {\boldsymbol{A}} 来说

    \begin{split} {\boldsymbol{A}} =& { {\boldsymbol{Y}}_{(1)}}{\left( {{{\left( {{\boldsymbol{C}} \odot {\boldsymbol{B}}} \right)}^{\text{T}}}} \right)^\dagger }{\text{ = }}{{\boldsymbol{Y}}_{{\text{(1)}}}}\left( {{\boldsymbol{C}} \odot {\boldsymbol{B}}} \right)\\ & \cdot {\left( {{{\boldsymbol{C}}^{\text{T}}}{\boldsymbol{C}} \times {{\boldsymbol{B}}^{\text{T}}}{\boldsymbol{B}}} \right)^\dagger } \end{split} (34)

    因为 LT > K ,所以 \left( {{\boldsymbol{C}} \odot {\boldsymbol{B}}} \right) \in {\mathbb{C}^{LT \times K}} 是列满秩矩阵,且 {{\boldsymbol{Y}}_{(1)}} \in {\mathbb{C}^{r \times LT}} ,则求解因子矩阵 {\boldsymbol{A}} 所需要的复杂度为 \mathcal{O}\left( {rLTK + rT{K^2} + {K^3}} \right) ,同理,求解因子矩阵 {\boldsymbol{B}} {\boldsymbol{C}} 也需要同样的复杂度,因此,本文所提出的算法复杂度为\mathcal{O}\left( 3\left( rLTK + rT{K^2} + {K^3} \right) \right)。文献[21]提出的ICE算法估计整个信道需要的复杂度为 \mathcal{O}\left( {4{N^3}\left( {{N_t} + {N_r}} \right) + 2{N^2}} \right) 。可以看出,文献[21]的复杂度量级最高为立方阶,且随着IRS元件数增加,复杂度是难以接受的。本文提出的算法复杂度量级最高为线性阶,对比发现,本文提出的算法具有很低的复杂度。

    本文针对IRS辅助THz无线通信场景下的信道估计问题进行研究,将IRS辅助通信信道写成统一的表达形式,利用不同块不同时隙将接收信号构建成3维矩阵,最后用基于正则平行因子分解的算法对信道进行准确的估计。仿真结果表明,所提出的信道分解估计算法不仅能有效地对信道参数进行估计,通过实验对比,本文提出的算法在估计精度和复杂度上要优于其他算法。

  • 加载中
计量
  • 文章访问数:  2517
  • HTML全文浏览量:  143
  • PDF下载量:  686
  • 被引次数: 0
出版历程
  • 收稿日期:  2003-02-27
  • 修回日期:  2003-08-08
  • 刊出日期:  2004-06-19

目录

/

返回文章
返回