一种RIS辅助的雷达通信一体化系统波形设计方法

本发明涉及一种ris辅助的雷达通信一体化系统波形设计方法，属于雷达通信一体化波束赋形技术领域。

背景技术：

雷达通信一体化(dual-functionalradarcommunications,dfrc)系统可以在作为基站服务用户的同时通过接收回波信号检测目标，其以集成雷达和通信功能的潜力在多个领域引起了广泛兴趣，例如车联网、无人机通信和传感以及多功能射频系统等等。为了同时支持信息传输和目标检测，许多工作致力于dfrc的波形设计。雷达通信联合波形设计的目标是实现预期雷达波形的同时满足信干噪比(signal-to-interference-plus-noiseratio,sinr)的要求。f.liu等研究了联合mimo雷达检测和多用户多输入单输出(mumiso)通信的波束成形设计。考虑到恒模约束，f.liu等提出了一种有效的分支定界方法和基于流形优化的方法来设计dfrc波形。进一步，f.liu等考虑了旁瓣抑制问题，开发了梯度下降方法来优化在每个天线功率预算下的通信和雷达指标的加权总和。上述研究表明可以通过适当设计dfrc波形来实现雷达和通信性能(例如雷达波束图和多用户干扰)之间的平衡。

最近，可重构智能表面(reconfigureintelligentsurface,ris)作为一项革命性的技术，因其在减少多用户干扰方面的潜力而备受关注。k.huang等设计了一种低复杂度的方法来优化ris的无源波束成形向量，以在多播和多用户下行链路传输中最大化用户的最小snr/sinr。j.wang等研究了在非精确信道状态信息下ris辅助的mu-miso系统的鲁棒波束成形和相移设计，并提出了一种基于交替优化的方法，以在中断概率约束下将发射功率降至最低。但是，在dfrc系统中使用ris缓解多用户干扰仍然是一个未解决的问题。

技术实现要素：

本发明的目的在于针对雷达通信一体化系统中存在的多用户干扰较大以及雷达和通信性能难以平衡、通信性能可能会受到严重损失的技术缺陷，提出了一种ris辅助的雷达通信一体化系统波形设计方法，通过基于流形优化的交替方法，联合设计dfrc波形和ris相位变换矩阵，并通过加权设计以发挥ris在缓解多用户干扰的作用，从而提高系统吞吐量；且能平衡雷达与通信功能的性能，使获得的波形更好地匹配所需的雷达波束图，并同时满足通信性能要求。

所述雷达通信一体化系统波形设计方法依托于一个ris辅助dfrc系统，该ris辅助dfrc系统包括一个n天线的均匀线性阵列组成的基站，k个单天线的用户以及一个有l个反射元件的可重构智能表面，在本系统中使用的通信帧的长度为m，即每一个通信帧包含m个符号；

具体包括以下步骤：

步骤1、初始化ris辅助的dfrc系统的组成架构及输入条件，具体包括如下子步骤：

步骤1.a：初始化ris辅助dfrc系统的组成架构；

其中，ris辅助dfrc系统的组成架构即通信系统模型，包括发射信号矩阵、ris辅助dfrc系统的输出、基站到用户间的基带等效信道、ris到用户间的基带等效信道、ris的相位变换矩阵以及用户处的噪声，通过公式(1)表示：

其中，x表示发射信号矩阵，y表示ris辅助dfrc系统的输出；gbu表示基站到用户间的基带等效信道；gru表示ris到用户间的基带等效信道；θ表示ris的相位变换矩阵，相位变换矩阵中对角线上的各元素满足模为1的约束；gbr表示基站到ris间的基带等效信道；w表示用户处的噪声；表示基站到用户间有ris辅助时的等效信道；

步骤1.b：初始化ris辅助的dfrc系统的输入条件；

其中，输入条件包括总功率限制、目标雷达波形的协方差矩阵及用户的信源符号矩阵；且输入条件包括总功率限制，记为p0；目标雷达波形的协方差矩阵，记为rd；用户的信源符号矩阵，记为s；

步骤2、初始化收敛精度阈值、黎曼梯度收敛阈值以及最大迭代次数；

其中，收敛精度阈值，记为η，黎曼梯度收敛阈值，记为∈，最大迭代次数，记为γ；

步骤3、初始化ris的相位变换矩阵θ为单位矩阵；

步骤4、将目标雷达波形的协方差矩阵rd进行cholesky分解、计算其分解得到的下三角复数矩阵的共轭转置与等效信道、信源符号矩阵乘后矩阵的奇异值分解，再计算波形样板矩阵u，具体包括如下子步骤：

步骤4.a：计算rd的cholesky分解rd＝ee^h，得到下三角复数矩阵e；

其中，[·]^t表示矩阵[·]的转置；

步骤4.b：计算得到的奇异值分解的左奇异矩阵和右奇异矩阵

其中，[·]^h表示矩阵[·]的共轭转置；

步骤4.c：基于步骤4.a得到的下三角复数矩阵e、步骤4.b得到的左奇异矩阵右奇异矩阵一个通信帧中的符号数m以及n行m列的单位阵in×m，计算得到波形样板矩阵u，具体通过(1)实现：

其中，in×m为一个n行m列的单位阵；

步骤5、通过特征值分解和黄金分割搜索使加权的多用户干扰和雷达波形误差最小化，输出发射信号矩阵x；

其中，使加权的多用户干扰和雷达波形误差最小化通过公式(2)的优化函数表达：

其中，||·||f表示frobenius范数，即对||·||中矩阵的各元素求平方和再开方，ρ表示平衡雷达和通信性能的权重系数，u为波形样板矩阵；

步骤5，具体包括如下子步骤：

步骤5.a：将优化函数公式(2)变形为

其中，in表示n行n列的单位阵；

步骤5.b：计算a^ha的特征值分解，得到包含特征向量的正交矩阵v、包含特征值的对角矩阵λ以及包含特征向量的正交矩阵的共轭转置矩阵v^h；

步骤5.c：设定a^ha特征值λ的搜索范围，并使用黄金分割搜索寻找中特征值λ的最优解λopt；

其中，搜索范围为[-λmin,b]，λmin为a^ha的最小特征值，b≥0是搜索上界；[·]^-1表示矩阵[·]的逆；

步骤5.d：通过(3)计算发射信号矩阵x：

其中，表示矩阵的moore-penrose伪逆；

步骤6、固定步骤5中计算得到的发射信号矩阵x，将步骤4.a得到的下三角复数矩阵e的共轭转置与发射信号矩阵乘后矩阵进行奇异值分解，再计算得到新的波形样板矩阵，具体包括如下子步骤：

具体包括如下子步骤：

步骤6.a：将e^hx进行奇异值分解，得到e^hx的左奇异矩阵和右奇异矩阵

步骤6.b：基于下三角复数矩阵e、步骤6.a得到的左奇异矩阵右奇异矩阵一个通信帧中的符号数m以及n行m列的单位阵in×m，通过(4)计算波形样板矩阵u：

步骤7、固定发射信号矩阵x，利用黎曼最陡下降方法计算得到更新后的相位变换矩阵θ，具体包括如下子步骤：

步骤7.a：规定θ＝θ1l，其中，1l表示长度为l的全为1的列向量，初始化θ0为干扰最小化函数的一个可行解；

其中，d＝gbrx(gbux-s)^hgru，d＝[d1,1,·,dl,l]^t，b⊙c表示b与c的哈达玛积，且b⊙c为各个位置的元素等于两个输入矩阵相同位置元素的乘积的矩阵；[·]^*表示矩阵[·]的共轭；

步骤7.b：初始化k＝0；

步骤7.c：依据流形优化理论中的armijo准则计算步长αk；

步骤7.d：计算f(θk)的欧几里得梯度

步骤7.e：计算f(θk)的黎曼梯度

步骤7.f：更新其中(·)i表示对向量中的逐个元素进行操作；

步骤7.g：更新k＝k+1，判断黎曼梯度是否收敛，即f是否小于收敛精度阈值∈或达到最大迭代次数γ，若是，则更新后的θ为以θk中元素为对角线元素的对角阵，结束步骤7，否则跳至步骤7.c；

步骤8、判断优化函数(2)是否收敛，即优化后的优化函数值与更新前的优化函数值的差是否小于收敛精度阈值η或达到最大迭代次数γ，若是，则结束本方法，输出发射信号矩阵x、波形样板矩阵u以及相位变换矩阵θ，否则跳至步骤5。

有益效果

本发明提出的一种ris辅助的雷达通信一体化系统波形设计方法，与现有技术相比，具有以下有益效果：

1、所述方法，能够使雷达通信一体化系统中多用户干扰显著减少，从而获得更高的通信吞吐量；

2、所述方法，能够平衡雷达通信一体化系统中雷达和通信功能的性能，使获得的雷达波形更好地匹配所需雷达波束图的同时满足通信性能要求。

附图说明

图1为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1整体方法流程图；

图2为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中初始化波形样板矩阵的流程图；

图3为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中求解发射信号矩阵的流程图；

图4为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中求解新的波形样板矩阵的流程图；

图5为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中求解更新后的相位变换矩阵的流程图；

图6为本发明一种ris辅助的雷达通信一体化系统波形设计方法实施例1中实施方法前后随着发送信噪比变化时的通信速率仿真结果图；

图7为本发明一种ris辅助的雷达通信一体化系统波形设计方法实施例1中实施方法前后雷达波束仿真结果图。

具体实施方式

下面结合附图和具体实施例对本发明所述的一种ris辅助的雷达通信一体化系统波形设计方法进行详细说明。

实施例1

本实施例详细阐述了本发明所述的一种ris辅助的雷达通信一体化系统波形设计方法具体实施时的步骤。

本实例考虑一个ris辅助的雷达通信一体化系统，信道矩阵的每项都服从标准复数高斯分布，基站的天线数设置为n＝20，总功率限制设置为p0＝20dbm，通过基于流形优化的交替方法，实现联合设计dfrc波形和ris的相位变换矩阵，通过加权设计以发挥ris在缓解多用户干扰的作用，从而提高系统吞吐量；并且平衡了雷达与通信功能的性能，使获得的波形更好地匹配所需的雷达波束图的同时满足通信性能要求；

图1为本发明所述的一种ris辅助的雷达通信一体化系统波形设计方法及实施例1整体方法流程图；

从图1可以看出本发明所述的一种ris辅助的雷达通信一体化系统波形设计方法的详细流程，具体到本实施例，方法操作流程如下：

步骤1、初始化ris辅助的dfrc系统的组成架构及输入条件；

步骤1.a：初始化ris辅助dfrc系统的组成架构，本实例考虑一个ris辅助dfrc系统，其组成架构如下所述，基站天线数量n为20，单天线用户个数k为4个，可重构智能表面反射元件的个数l为16个，设定每一帧包含m＝30个符号；

通信系统模型表示为

其中，x表示发射信号矩阵，y表示ris辅助dfrc系统的输出；

gbu＝[gbu,1,…,gbu,k]^t为一个k×n(k＝4，n＝20，即4×20)的复数矩阵，表示基站到用户间的基带等效信道；gru＝[gru,1,…,gru,k]^t为一个k×l(k＝4，l＝16，即4×16)的复数矩阵，表示ris到用户间的基带等效信道；gbr＝[gbr,1,…,gbr,l]^t为一个l×n(l＝16，n＝20，即16×20)的复数矩阵，表示基站到ris间的基带等效信道；w＝[w1,…,wk]^t为一个k×m(k＝4，m＝30，即4×30)的复数矩阵，表示用户处的噪声，其中wk服从均值为0的高斯分布，表示基站到用户间有ris辅助时的等效信道；

步骤1.b：初始化ris辅助的dfrc系统的输入条件；

初始化总功耗预算p0为20dbm，初始化目标雷达波形的协方差矩阵rd为n×n(n＝20，即20×20)的矩阵，初始化用户的信源符号矩阵s为k×m(k＝4，m＝30，即4×30)的矩阵，其中各元素为随机生成的qpsk符号；

步骤2、初始化收敛精度阈值η为0.01，黎曼梯度收敛阈值∈为0.01，最大迭代次数γ为100；

步骤3、初始化ris的相位变换矩阵θ为单位矩阵；

步骤4、将目标雷达波形的协方差矩阵rd进行cholesky分解、计算其分解得到的下三角复数矩阵的共轭转置与等效信道、信源符号矩阵乘后矩阵的奇异值分解，再计算波形样板矩阵u，具体包括如下子步骤：

图2为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中初始化波形样板矩阵的流程图；

从图2可以看出步骤4中通过cholesky分解和奇异值分解初始化波形样板矩阵u的方法，具体到本实施例，方法操作流程如下，具体包括如下子步骤：

步骤4.a：计算rd的cholesky分解rd＝ee^h，其中e是一个n×n的下三角复数矩阵；

步骤4.b：计算的奇异值分解，得到的左奇异矩阵和右奇异矩阵

步骤4.c：计算得到波形样板矩阵的初始值

步骤5、将相位变换矩阵θ和波形样板矩阵u固定，通过特征值分解和黄金分割搜索使加权的多用户干扰和雷达波形误差最小化，即将上式称为优化函数，其中ρ表示平衡雷达和通信性能的权重系数，从而计算得到发射信号矩阵x；

具体到本实施例，ρ分别取0.2，0.01，0；

图3为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中求解发射信号矩阵的流程图；

从图3可以看出步骤5中通过特征值分解和黄金分割搜索使加权的多用户干扰和雷达波形误差最小化，从而得到发射信号矩阵的方法，具体到本实施例，方法操作流程如下，具体包括如下子步骤：

步骤5.a：计算

步骤5.b：计算a^ha的特征值分解a^ha＝vλv^h，得到包含特征向量的正交矩阵v、包含特征值的对角矩阵λ以及包含特征向量的正交矩阵的共轭转置矩阵v^h；

步骤5.c：设定搜索范围是[-λmin,b]，其中λmin为a^ha的最小特征值，b≥0是搜索上界，使用黄金分割搜索寻找中特征值λ的最优解λopt；

其中，[·]^-1表示矩阵[·]的逆；

具体到本实施例，搜索上界b设为100；

步骤5.d：计算得到最优发射信号矩阵其中表示矩阵的moore-penrose伪逆；

步骤6、固定步骤5中计算得到的发射信号矩阵x，将步骤4.a得到的下三角复数矩阵e的共轭转置与发射信号矩阵乘后矩阵进行奇异值分解，再计算得到新的波形样板矩阵；

图4为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中求解新的波形样板矩阵的流程图；

从图4可以看出步骤6中计算得到新的波形样板矩阵u的方法，具体到本实施例，方法操作流程如下，具体包括如下子步骤：

步骤6.a：计算e^hx的奇异值分解，得到e^hx的左奇异矩阵和右奇异矩阵

步骤6.b：计算得到新的波形样板矩阵

步骤7、将发射信号矩阵x固定，利用黎曼最陡下降方法计算得到更新后的相位变换矩阵θ；

图5为本发明一种ris辅助的雷达通信一体化系统波形设计方法及实施例1中求解更新后的相位变换矩阵的流程图；

从图5可以看出步骤7中通过黎曼最陡下降方法计算得到更新后的相位变换矩阵θ的方法，具体到本实施例，方法操作流程如下，具体包括如下子步骤：

步骤7.a：规定θ＝θ1l，其中，1l表示长度为l的全为1的列向量，初始化θ0为干扰最小化函数的一个可行解，其中|θl|＝1,d＝gbrx(gbux-s)^hgru，d＝[d1,1,·,dl,l]^t，b⊙c表示b与c的哈达玛积，且b⊙c为各个位置的元素等于两个输入矩阵相同位置元素的乘积的矩阵；[·]^*表示矩阵[·]的共轭；

步骤7.b：初始化k＝0；

步骤7.c：依据流形优化理论中的armijo准则计算步长αk；

步骤7.d：计算欧几里得梯度

步骤7.e：计算黎曼梯度

步骤7.f：更新其中(·)i表示对向量中的逐个元素进行操作；

步骤7.g：更新k＝k+1，判断黎曼梯度是否收敛，即f是否小于收敛精度阈值∈或达到最大迭代次数γ，若是，则更新后的θ为以θk中元素为对角线元素的对角阵，结束步骤7，否则跳至步骤7.c；

步骤8、判断优化函数是否收敛，即优化后的优化函数值与更新前的优化函数值的差是否小于收敛精度阈值η或达到最大迭代次数γ，若是，则结束本方法，输出发射信号矩阵x、波形样板矩阵u以及相位变换矩阵θ，若否则跳至步骤5；

至此，经过步骤1到步骤8，完成了一种可重构智能表面辅助的雷达通信一体化系统波形设计。

图6为本发明一种ris辅助的雷达通信一体化系统波形设计方法实施例1中实施方法前后随着发送信噪比变化时的通信速率仿真结果图；

图6横坐标为发送信噪比，依次取值为0db，2db，4db，6db，8db，10db，纵坐标为通信速率，仿真实验对六种情况进行了对比分析：1)有ris辅助，平衡雷达和通信性能的权重系数ρ为0.2，即本方法；2)无ris辅助，平衡雷达和通信性能的权重系数ρ为0.2；3)有ris辅助，平衡雷达和通信性能的权重系数ρ为0，严格满足雷达波束要求，即本方法；4)无ris辅助，平衡雷达和通信性能的权重系数ρ为0，严格满足雷达波束要求；5)有ris辅助，平衡雷达和通信性能的权重系数ρ为0.01，即本方法；6)无ris辅助，平衡雷达和通信性能的权重系数ρ为0.01；

由图6可以看出在ris的辅助下，通信速率可以显著地提高，有效提高了通信吞吐量；同时，随着平衡雷达和通信性能的权重系数ρ的增大，通信速率也随之增加，但代价是雷达波形匹配度下降；

图7为本发明一种ris辅助的雷达通信一体化系统波形设计方法实施例1中实施方法前后雷达波束仿真结果图；

仿真实验对四种情况进行了分析：1)目标波束，即平衡雷达和通信性能的权重系数ρ为0，严格满足雷达波形要求；2)无ris辅助，平衡雷达和通信性能的权重系数ρ为0.01；3)有ris辅助，平衡雷达和通信性能的权重系数ρ为0.01，即本方法；4)有ris辅助，平衡雷达和通信性能的权重系数ρ为0.2，即本方法；

由图7可以看出在ris的辅助下，雷达波形匹配度明显比没有ris辅助时更高，使获得的雷达波形更好地匹配所需雷达波束图；

由图6和图7可以看出在ris的辅助下，可以实现雷达和通信功能更好的平衡。

以上所述为本发明的较佳实施例而已，本发明不应该局限于该实施例和附图所公开的内容。凡是不脱离本发明所公开的精神下完成的等效或修改，都落入本发明保护的范围。