CN115396263A - 一种otfs通信感知一体化信号目标参数估计方法 - Google Patents

一种otfs通信感知一体化信号目标参数估计方法 Download PDF

Info

Publication number
CN115396263A
CN115396263A CN202210903235.6A CN202210903235A CN115396263A CN 115396263 A CN115396263 A CN 115396263A CN 202210903235 A CN202210903235 A CN 202210903235A CN 115396263 A CN115396263 A CN 115396263A
Authority
CN
China
Prior art keywords
domain
doppler
target
representing
index
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210903235.6A
Other languages
English (en)
Other versions
CN115396263B (zh
Inventor
刘喜庆
熊礼亮
闫实
彭木根
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN202210903235.6A priority Critical patent/CN115396263B/zh
Publication of CN115396263A publication Critical patent/CN115396263A/zh
Application granted granted Critical
Publication of CN115396263B publication Critical patent/CN115396263B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/024—Channel estimation channel estimation algorithms
    • H04L25/0256—Channel estimation using minimum mean square error criteria
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/0212—Channel estimation of impulse response
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/024—Channel estimation channel estimation algorithms
    • H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

本发明公开了一种OTFS通信感知一体化信号目标参数估计方法。包括以下步骤：步骤1，生成由发射天线发送的OTFS通信感知一体化信号：步骤2，在一体化信号感知到目标后，由目标反射回波信号至通信感知一体化接收机：步骤3，进行下变频处理与OTFS解调，获得DD域接收调制符号矩阵：步骤4，完成雷达信道与目标参数估计：步骤5，对OTFS通信感知一体化信号目标参数估计的结果进行校准：本发明优点是当目标的距离和速度分别为距离分辨率和速度分辨率的整数倍时能够准确估计目标参数信息且提出的DD域能量质心谱校正法能够将分数倍情况下产生的参数估计误差大大降低，通过校准获得更准确的参数估计结果并且极大地简化了计算过程。

Description

一种OTFS通信感知一体化信号目标参数估计方法
技术领域
本发明属于无线通信和雷达信号处理技术领域，具体涉及一种基于MMSE准则与DD域能量质心谱校正的OTFS通信感知一体化信号目标参数估计方法。
背景技术
现有的无线通信系统提供通信服务以解决用户的话音及数据业务需求。未来第六代移动通信系统(6G)除了实现传统通信网络的功能，还将为新兴业务如扩展现实、智能车联网提供高精度传感服务，因而无线通信与雷达感知技术的融合受到了广泛的关注。长期以来无线通信和雷达感知并行发展、交叉有限，然而随着两者技术的演进，通信系统与感知系统逐渐具备联合设计的可能与潜力。在系统组成结构上，通信与感知在信号处理算法、设备甚至系统架构上有许多共同点；在工作频段兼容方面，通信系统频段的抬高将导致两者的频段出现越来越多的交叠，两者共享频谱将具有潜在可能；此外将通信系统与感知系统集成于一体能够减小设备尺寸与成本。因此为了解决新兴业务的通感双功能需求并提高无线电频谱效率及降低硬件成本，需要联合设计通信系统与感知系统从而共享相同的频段和硬件设备，最终实现通信感知一体化。通信与感知的深度融合依赖于通信感知一体化信号处理，对通感一体化信号波形进行优化与融合设计，使得一种波形能同时满足通信传输和雷达感知的需求，这是实现目标感知和数据通信两种功能集成的关键挑战。
目前，一种新型二维调制技术即正交时频空(Orthogonal Time FrequencySpace，OTFS)，其在时延-多普勒(Delay-Doppler，DD)域调制符号信息能同时获得时间域和频率域全分集增益，因而对多普勒频移并不敏感，即使在高移动性场景下仍能实现良好的通信性能。雷达感知的目标参数是距离和速度的元组形式，分别对应于时延和多普勒信息，这与OTFS在DD域进行信号调制和信道表示的思路相契合，目前，现有技术的解决思路是基于发送调制符号与接收调制符号估计雷达感知信道完成对目标距离和速度的参数估计。但目前现有技术在感知性能方面还有待提升，如峰值旁瓣比和雷达图像信噪比，此外当目标的实际距离和速度为分数倍的感知分辨率时，如何降低感知参数估计误差并获得更加精确的估计结果也成为一个急待解决的技术问题。
发明内容
针对现有技术存在的问题，本发明提供了一种基于MMSE准则与DD域能量质心谱校正的OTFS通信感知一体化信号目标参数估计方法，本发明所述方法包括以下步骤：
步骤1，通信感知一体化发射机对通信信息序列进行星座映射、OTFS调制与射频调制，生成由发射天线发送的OTFS通信感知一体化信号：
步骤1.1，通信信源生成由二进制比特构成的数据信息序列，将信息序列输入至星座映射器并采用正交振幅调制(Quadrature Amplitude Modulation，QAM)进行数字调制获得发送调制符号向量；
Figure BDA0003771656370000011
表示长度为MN的发送调制符号向量，M表示OTFS帧的子载波数量，N表示OTFS帧的符号数量；
步骤1.2，将待传输的发送调制符号向量排列成N行M列的符号矩阵并摆放于二维时延-多普勒(Delay-Doppler，DD)域平面栅格Γ获得DD域发送调制符号矩阵
Figure BDA0003771656370000021
DD域平面栅格定义为Γ如以下公式(1)所示：
Figure BDA0003771656370000022
以上公式(1)中，k和l分别表示多普勒域的索引和时延域的索引，T和Δf分别表示OTFS符号周期和子载波间隔，并满足TΔf＝1以保持正交性，(k)N表示对k取多普勒频移运算，如以下公式(2)所示：
Figure BDA0003771656370000023
步骤1.3，对DD域发送调制符号进行0TFS调制，获得时域上的基带发射信号：
步骤1.3.1，通过辛有限傅里叶逆变换(Inverse Symplectic Finite FourierTransform，ISFFT)将DD域发送调制符号矩阵映射为时间-频率(Time-Frequency，TF)域发送调制符号矩阵
Figure BDA0003771656370000024
占据TF域平面栅格Λ上的N个符号和M个子载波对应的时频格点，映射过程如以下公式(3)所示：
Figure BDA0003771656370000025
其中，ISFFT(·)表示应用辛有限傅里叶逆变换，{XTF[n，m]∈XTF，n＝0，…，N-1，m＝0，…，M-1}表示TF域平面栅格Λ上的第n个时间和第m个频率格点处的发送调制符号，DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的发送调制符号表示为{XDD[k，l]∈XDD，k＝0，…，N-1，l＝0，…，M-1}；
TF域平面栅格定义为Λ，如以下公式(4)所示：
Λ＝{(nT，mΔf)，n＝0，1，...，N-1，m＝0，1，...，M-1}......(4)，
其中，n和m分别为时间域的索引和频率域的索引；
步骤1.3.2，对TF域发送调制符号矩阵进行海森堡变换(Heisenberg Transform)，通过矩形脉冲成形滤波器获得时域基带发射信号，变换公式如以下公式(5)所示：
Figure BDA0003771656370000026
其中，s(t)表示时域基带发射信号，gtx(t)为发送端持续时间为T的矩形脉冲成形滤波器，当0≤t≤T时有gtx(t)＝1，否则gtx(t)＝0；
步骤1.4，对时域上的基带发射信号进行射频调制，产生的由发射天线发送的OTFS通信感知一体化信号如以下公式(6)所示：
Figure BDA0003771656370000027
其中，sRF(t)表示由发射天线发送的OTFS通信感知一体化信号，fc表示射频载波频率；
步骤2，通信感知一体化发射机发送OTFS通信感知一体信号进行目标感知，在一体化信号感知到目标后，由目标反射回波信号至通信感知一体化接收机：
步骤2.1，根据所感知目标的特征建立DD域雷达信道脉冲响应的模型；
DD域雷达信道脉冲响应的模型如以下公式(7)所示：
Figure BDA0003771656370000031
其中，h(τ，v)表示DD域雷达信道，τ和v分别表示时延和多普勒频移，P表示感知目标的数量，hp表示第p个目标的复信道增益，p＝1，…，P，τp和vp分别表示通信感知一体化发射机与第p个目标的往返时延和多普勒频移，τp满足τp＝2Rp/c，vp满足vp＝2fcVp/c，Rp和Vp分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度，c代表光速，δ(·)为狄利克雷函数；
步骤2.2，将发送的OTFS通信感知一体化信号经感知目标反射后，产生回波并返回至通信感知一体化接收机的这一过程，等效于发送的0TFS通信感知一体化信号经过步骤2.1所述DD域雷达信道的射频等效信道到达通信感知一体化接收机处，通信感知一体化接收机处收到回波信号的等效过程如以下公式(8)所示：
rRF(t)＝∫∫hRF(τ，v)sRF(t-τ)ej2πv(t-τ)dτdv......(8)，
其中，rRF(t)表示通信感知一体化接收机处收到的回波信号，hRF(τ，v)表示步骤2.1所述的DD域雷达信道的射频等效信道；其中，DD域雷达信道的射频等效信道如以下公式(9)所示：
Figure BDA0003771656370000032
步骤3，通信感知一体化接收机在接收到由目标反射的回波信号后，进行下变频处理与OTFS解调，获得DD域接收调制符号矩阵：
步骤3.1，通信感知一体化接收机对回波信号进行下变频处理，获取时域基带接收信号；
时域基带接收信号如以下公式(10)所示：
Figure BDA0003771656370000033
其中，r(t)表示通信感知一体化接收机处的时域基带接收信号；
步骤3.2，对时域基带接收信号进行OTFS解调，获得DD域接收调制符号矩阵
Figure BDA0003771656370000034
步骤3.2.1，对时域基带接收信号应用魏格纳变换(Wigner Transform)，得到TF域接收调制符号矩阵
Figure BDA0003771656370000035
变换公式如以下公式(11)所示：
Figure BDA0003771656370000036
其中，{YTF[n，m]∈YTF，n＝0，…，N-1，m＝0，…，M-1}表示TF域平面栅格Λ上的第n个时间和第m个频率格点处的接收调制符号，grx(t)为接收端的矩形脉冲成形滤波器，(·)*表示求复共轭函数；
步骤3.2.2，采用辛有限傅里叶变换(Symplectic Finite Fourier Transform，SFFT)将TF域接收调制符号矩阵反映射，得到DD域接收调制符号矩阵
Figure BDA0003771656370000041
映射过程如以下公式(12)所示：
Figure BDA0003771656370000042
其中，SFFT(·)表示应用辛有限傅里叶变换，{YDD[k，l]∈YDD，k＝0，…，N-1，l＝0，…，M-1}表示DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的接收调制符号；
步骤4，通信感知一体化接收机根据DD域接收调制符号与DD域发送调制符号进行雷达感知信号处理，完成雷达信道与目标参数估计：
步骤4.1，将步骤2.1中所述的DD域雷达信道使用离散形式表达，如以下公式(13)所示：
Figure BDA0003771656370000043
其中，h[k，l]表示DD域雷达信道离散响应，即多普勒域索引k和时延域索引l处的复信道增益值，多普勒域索引k和时延域索引l分别对应于多普勒频移v＝(k)N/(NT)和时延τ＝l/(MΔf)，1/(NT)和1/(MΔf)分别表示多普勒分辨率和时延分辨率；
其中，第p个目标的往返时延τp与时延分辨率的关系为τp＝lp/(MΔf)，多普勒频移νp与多普勒分辨率的关系为vp＝kp/(NT)且关系表示第p个目标位于DD域平面栅格Γ上第lp个时延域索引和第kp个多普勒域索引处，当k＝kp且l＝lp时则h[k，l]＝hp，当在第l个时延域索引和第k个多普勒域索引处未感知到目标时则h[k，l]＝0；
步骤4.2，根据步骤4.1所述的离散形式DD域雷达信道，获得DD域发送调制符号与DD域接收调制符号之间的输入输出关系，DD域发送调制符号与DD域接收调制符号之间的输入输出关系如以下公式(14)所示：
Figure BDA0003771656370000044
其中，<·>N和<·>M分别表示模N和模M运算，ω[k，l]表示DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的加性高斯白噪声，服从均值为零且方差(即平均功率)为σ2的复高斯分布，α[k，l]表示由时域样本的循环移位引起的多普勒域线性相移，α[k，l]如以下公式(15)所示：
Figure BDA0003771656370000045
步骤4.3，根据DD域接收调制符号和DD域发送调制符号进行基于最小均方误差准则(Minimum Mean Square Error，MMSE)的DD域雷达信道与目标参数估计；
步骤4.3.1，将步骤4.2所述的输入输出关系使用矩阵形式表达式如以下公式(16)所示：
Figure BDA0003771656370000051
其中，
Figure BDA0003771656370000052
表示DD域接收调制符号向量，yDD的第k+Nl个元素为YDD[k，l]，且0≤k+Nl≤MN-1，
Figure BDA0003771656370000053
表示DD域雷达信道向量，hDD的第k+Nl个元素为h[k，l]，
Figure BDA0003771656370000054
表示DD域噪声向量，wDD的第k+Nl个元素为w[k，l]，符号矩阵
Figure BDA0003771656370000055
的第i行第j列元素如以下公式(17)所示：
Figure BDA0003771656370000056
其中，i和j分别满足0≤i＝k+Nl≤MN-1，O≤j＝k′+Nl′≤MN-1；
步骤4.3.2，根据步骤4.3.1.所述的矩阵形式输入输出关系，基于MMSE准则进行DD域雷达信道估计，DD域雷达信道向量估计值如以下公式(18)所示：
Figure BDA0003771656370000057
其中，
Figure BDA0003771656370000058
表示基于MMSE准则估计的DD域雷达信道向量，(·)-1表示矩阵求逆，(·)H表示复共轭转置，Ps和σ2分别表示发送的OTFS通信感知一体化信号的平均功率和噪声的平均功率，IMN表示维度为MN×MN的单位矩阵；
步骤4.3.3，根据步骤4.3.2中基于MMSE准则估计的DD域雷达信道向量获取感知目标在DD域平面栅格Γ上的时延域索引和多普勒索引的估计结果；
其中，
Figure BDA0003771656370000059
表示估计的DD域雷达信道离散响应，即多普勒域索引k和时延域索引l处的复信道增益估计值，
Figure BDA00037716563700000510
等于DD域雷达信道向量估计值
Figure BDA00037716563700000511
的第k+Nl个元素，估计的DD域雷达信道离散响应其峰值对应的多普勒域索引和时延域索引分别为感知目标在DD域平面栅格Γ上的多普勒域索引和时延域索引的估计结果，分别表示为
Figure BDA00037716563700000512
和
Figure BDA00037716563700000513
步骤4.3.4，由感知目标在DD域平面栅格Γ上的多普勒域索引和时延域索引的估计结果进行目标参数估计，估计结果如以下公式(19)、公式(20)所示：
Figure BDA00037716563700000514
Figure BDA00037716563700000515
以上公式(19)、公式(20)中，
Figure BDA00037716563700000516
和
Figure BDA00037716563700000517
分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度的参数估计结果，
Figure BDA0003771656370000061
和
Figure BDA0003771656370000062
分别表示通信感知一体化发射机与第p个目标的往返时延和多普勒频移的参数估计结果，
Figure BDA0003771656370000063
满足
Figure BDA0003771656370000064
Figure BDA0003771656370000065
满足
Figure BDA0003771656370000066
Figure BDA0003771656370000067
和
Figure BDA0003771656370000068
分别表示第p个目标在DD域平面栅格Γ上的时延域索引和多普勒域索引的估计结果；
步骤5，在DD域平面栅格Γ上采用DD域能量质心谱校正法对步骤4中OTFS通信感知一体化信号目标参数估计的结果进行校准：
步骤5.1，对步骤4.3.3中估计的DD域雷达信道离散响应进行样值比选，分别实施行比选和列比选，根据比选结果求向量元素中的最大值,即DD域雷达信道离散响应的峰值；
行比选表示按行比较选取每一行的最大值，行比选的结果为向量
Figure BDA0003771656370000069
满足
Figure BDA00037716563700000610
Figure BDA00037716563700000611
表示DD域雷达信道离散响应
Figure BDA00037716563700000612
中第k行的最大值，k为多普勒域索引，列比选表示按列比较选取每一列的最大值，列比选的结果为向量
Figure BDA00037716563700000613
满足
Figure BDA00037716563700000614
Figure BDA00037716563700000615
表示DD域雷达信道离散响应
Figure BDA00037716563700000616
中第l列的最大值，l为时延域索引，峰值
Figure BDA00037716563700000617
对应于时延域索引
Figure BDA00037716563700000618
和多普勒域索引
Figure BDA00037716563700000619
步骤5.2，以DD域雷达信道离散响应的峰值为中心，向上取X个样值、向下取个X样值，同时向左取Y个样值、向右取Y个样值，并取样值的时延域索引与多普勒域索引；
向上取X个样值的结果为
Figure BDA00037716563700000620
下取X个样值的结果为
Figure BDA00037716563700000621
左取Y个样值的结果为
Figure BDA00037716563700000622
右取Y个样值的结果为
Figure BDA00037716563700000623
步骤5.3，根据DD域样值的总能量与样值在时延域维度、多普勒域维度的加权能量，计算时延域维度的质心，并行计算多普勒域维度的质心；
步骤5.3.1，计算步骤5.2中DD域样值的总能量，并行计算样值在时延域维度的加权能量及多普勒域维度的加权能量，计算公式如以下公式(21)至公式(23)所示：
Figure BDA00037716563700000624
Figure BDA00037716563700000625
以上公式(21)至公式(23)中，E表示DD域样值的总能量，El表示时延域维度的加权能量，Ek表示多普勒域维度的加权能量，
Figure BDA00037716563700000626
表示多普勒域索引k和时延域索引l处的能量，
集合
Figure BDA00037716563700000627
表示由步骤5.2中峰值与旁瓣值对应的时延域索引、多普勒域索引构成的集合，即
Figure BDA00037716563700000628
步骤5.3.2，计算时延域维度的质心，并行计算多普勒域维度的质心，计算公式如以下公式(24)、公式(25)所示：
Figure BDA0003771656370000071
Figure BDA0003771656370000072
其中，
Figure BDA0003771656370000073
表示时延域维度的质心，
Figure BDA0003771656370000074
表示多普勒域维度的质心；
步骤5.4，根据步骤5.3得到的时延域维度的质心与多普勒域维度的质心，对目标参数估计结果进行校准；
校准结果如以下公式(26)、公式(27)所示(26)、(27)：
Figure BDA0003771656370000075
Figure BDA0003771656370000076
其中，
Figure BDA0003771656370000077
和
Figure BDA0003771656370000078
分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度在经过DD域能量质心谱校正法校准后得到的参数估计结果，
Figure BDA0003771656370000079
和
Figure BDA00037716563700000710
分别表示经过DD域能量质心谱校正法校准后得到的时延域维度的质心和多普勒域维度的质心。
本发明具有的显著的有益效果在于：
1.本发明所述方法使用基于OTFS调制的通信感知一体化信号完成了目标感知与参数估计功能且实现了对目标距离和速度的参数估计，本发明所述方法能够与现有的OTFS通信系统相兼容并在支持高速移动通信的同时也能解决目标感知的需求，从而实现了通信感知一体化双功能。
2.本发明所述方法提出了针对OTFS通信感知一体化信号的基于MMSE准则的感知处理算法，当目标的距离和速度分别为距离分辨率和速度分辨率的整数倍时能够准确估计目标参数信息且提出的DD域能量质心谱校正法能够将分数倍情况下产生的参数估计误差大大降低，通过校准获得更准确的参数估计结果并且极大地简化了计算过程，即大大降低了计算的复杂度。
附图说明
图1是本发明所述方法基于MMSE准则与DD域能量质心谱校正的OTFS通信感知一体化信号目标参数估计方法的整体框图；
图2是本发明所述方法步骤5中DD域能量质心谱校正法的原理示意图；
图3是本发明所述方法实施例中归一化距离像性能曲线的仿真示意图；
图4是本发明所述方法实施例中归一化速度像性能曲线的仿真示意图；
图5是本发明所述方法的流程图。
具体实施方式
下面结合说明书附图，进一步详细描述本发明的具体实施方式。如图1、图2所示，
本发明所述方法包括以下步骤：
步骤1，通信感知一体化发射机对通信信息序列进行星座映射、OTFS调制与射频调制，生成由发射天线发送的OTFS通信感知一体化信号：
步骤1.1，通信信源生成由二进制比特构成的数据信息序列，将信息序列输入至星座映射器并采用正交振幅调制(Quadrature Amplirude Modulation，QAM)进行数字调制获得发送调制符号向量；
Figure BDA0003771656370000081
表示长度为MN的发送调制符号向量，M表示OTFS帧的子载波数量，N表示OTFS帧的符号数量；
步骤1.2，将待传输的发送调制符号向量排列成N行M列的符号矩阵并摆放于二维时延一多普勒(Delay-Doppler，DD)域平面栅格Γ获得DD域发送调制符号矩阵
Figure BDA0003771656370000082
DD域平面栅格定义为Γ如以下公式(1)所示：
Figure BDA0003771656370000083
以上公式(1)中，k和l分别表示多普勒域的索引和时延域的索引，T和Δf分别表示OTFS符号周期和子载波间隔，并满足TΔf＝1以保持正交性，(k)N表示对k取多普勒频移运算，如以下公式(2)所示：
Figure BDA0003771656370000084
步骤1.3，对DD域发送调制符号进行OTFS调制，获得时域上的基带发射信号：
步骤1.3.1，通过辛有限傅里叶逆变换(Inverse Symplectic Finite FourierTransform，ISFFT)将DD域发送调制符号矩阵映射为时间-频率(Time-Frequency，TF)域发送调制符号矩阵
Figure BDA0003771656370000085
占据TF域平面栅格Λ上的N个符号和M个子载波对应的时频格点，映射过程如以下公式(3)所示：
Figure BDA0003771656370000086
其中，ISFFT(·)表示应用辛有限傅里叶逆变换，{XTF[n，m]∈XTF，n＝0，…，N-1，m＝0，…，M-1}表示TF域平面栅格Λ上的第n个时间和第m个频率格点处的发送调制符号，DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的发送调制符号表示为{XDD[k，l]∈XDD，k＝0，…，N-1，l＝0，…，M-1}；
TF域平面栅格定义为Λ，如以下公式(4)所示：
Λ＝{(nT，mΔf)，n＝0，1，...，N-1，m＝0，1，...，M-1}......(4)，
其中，n和m分别为时间域的索引和频率域的索引；
步骤1.3.2，对TF域发送调制符号矩阵进行海森堡变换(Heisenberg Transform)，通过矩形脉冲成形滤波器获得时域基带发射信号，变换公式如以下公式(5)所示：
Figure BDA0003771656370000091
其中，s(t)表示时域基带发射信号，gtx(t)为发送端持续时间为T的矩形脉冲成形滤波器，当0≤t≤T时有gtx(t)＝1，否则gtx(t)＝0；
步骤1.4，对时域上的基带发射信号进行射频调制，产生的由发射天线发送的OTFS通信感知一体化信号如以下公式(6)所示：
Figure BDA0003771656370000092
其中，sRF(t)表示由发射天线发送的0TFS通信感知一体化信号，fc表示射频载波频率；
步骤2，通信感知一体化发射机发送0TFS通信感知一体信号进行目标感知，在一体化信号感知到目标后，由目标反射回波信号至通信感知一体化接收机：
步骤2.1，根据所感知目标的特征建立DD域雷达信道脉冲响应的模型；
DD域雷达信道脉冲响应的模型如以下公式(7)所示：
Figure BDA0003771656370000093
其中，h(τ，v)表示DD域雷达信道，τ和v分别表示时延和多普勒频移，P表示感知目标的数量，hp表示第p个目标的复信道增益，p＝1，…，P，τp和vp分别表示通信感知一体化发射机与第p个目标的往返时延和多普勒频移，τp满足τp＝2Rp/c，vp满足vp＝2fcVp/c，Rp和Vp分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度，c代表光速，δ(·)为狄利克雷函数；
步骤2.2，将发送的OTFS通信感知一体化信号经感知目标反射后，产生回波并返回至通信感知一体化接收机的这一过程，等效于发送的0TFS通信感知一体化信号经过步骤2.1所述DD域雷达信道的射频等效信道到达通信感知一体化接收机处，通信感知一体化接收机处收到回波信号的等效过程如以下公式(8)所示：
rRF(t)＝∫∫hRF(τ，v)sRF(t-τ)ej2πv(t-τ)dτdv......(8)，
其中，rRF(t)表示通信感知一体化接收机处收到的回波信号，hRF(τ，v)表示步骤2.1所述的DD域雷达信道的射频等效信道；其中，DD域雷达信道的射频等效信道如以下公式(9)所示：
Figure BDA0003771656370000094
步骤3，通信感知一体化接收机在接收到由目标反射的回波信号后，进行下变频处理与OTFS解调，获得DD域接收调制符号矩阵：
步骤3.1，通信感知一体化接收机对回波信号进行下变频处理，获取时域基带接收信号；时域基带接收信号如以下公式(10)所示：
Figure BDA0003771656370000095
其中，r(t)表示通信感知一体化接收机处的时域基带接收信号；
步骤3.2，对时域基带接收信号进行OTFS解调，获得DD域接收调制符号矩阵
Figure BDA0003771656370000101
步骤3.2.1，对时域基带接收信号应用魏格纳变换(Wigner Transform)，得到TF域接收调制符号矩阵
Figure BDA0003771656370000102
变换公式如以下公式(11)所示：
Figure BDA0003771656370000103
其中，{YTF[n，m]∈YTF，n＝0，…，N-1，m＝0，…，M-1}表示TF域平面栅格Λ上的第n个时间和第m个频率格点处的接收调制符号，grx(t)为接收端的矩形脉冲成形滤波器，(·)*表示求复共轭函数；
步骤3.2.2，采用辛有限傅里叶变换(Symplectic Finite Fourier Transform，SFFT)将TF域接收调制符号矩阵反映射，得到DD域接收调制符号矩阵
Figure BDA0003771656370000104
映射过程如以下公式(12)所示：
Figure BDA0003771656370000105
其中，SFFT(·)表示应用辛有限傅里叶变换，{YDD[k，l]∈YDD，k＝0，…，N-1，l＝0，…，M-1}表示DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的接收调制符号；
步骤4，通信感知一体化接收机根据DD域接收调制符号与DD域发送调制符号进行雷达感知信号处理，完成雷达信道与目标参数估计：
步骤4.1，将步骤2.1中所述的DD域雷达信道使用离散形式表达，如以下公式(13)所示：
Figure BDA0003771656370000106
其中，h[k，l]表示DD域雷达信道离散响应，即多普勒域索引k和时延域索引l处的复信道增益值，多普勒域索引k和时延域索引l分别对应于多普勒频移v＝(k)N/(NT)和时延τ＝l/(MΔf)，1/(NT)和1/(MΔf)分别表示多普勒分辨率和时延分辨率；
其中，第p个目标的往返时延τp与时延分辨率的关系为τp＝lp/(MΔf)，多普勒频移νp与多普勒分辨率的关系为vp＝kp/(NT)且关系表示第p个目标位于DD域平面栅格Γ上第lp个时延域索引和第kp个多普勒域索引处，当k＝kp且l＝lp时则h[k，l]＝hp，当在第l个时延域索引和第k个多普勒域索引处未感知到目标时则h[k，l]＝0；
步骤4.2，根据步骤4.1所述的离散形式DD域雷达信道，获得DD域发送调制符号与DD域接收调制符号之间的输入输出关系，DD域发送调制符号与DD域接收调制符号之间的输入输出关系如以下公式(14)所示：
Figure BDA0003771656370000111
其中，<·>N和<·>M分别表示模N和模M运算，ω[k，l]表示DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的加性高斯白噪声，服从均值为零且方差(即平均功率)为σ2的复高斯分布，α[k，l]表示由时域样本的循环移位引起的多普勒域线性相移，α[k，l]如以下公式(15)所示：
Figure BDA0003771656370000112
步骤4.3，根据DD域接收调制符号和DD域发送调制符号进行基于最小均方误差准则(Minimum Mean Square Error，MMSE)的DD域雷达信道与目标参数估计；
步骤4.3.1，将步骤4.2所述的输入输出关系使用矩阵形式表达式如以下公式(16)所示：
Figure BDA0003771656370000113
其中，
Figure BDA0003771656370000114
表示DD域接收调制符号向量，yDD的第k+Nl个元素为YDD[k，l]，且0≤k+Nl≤MN-1，
Figure BDA0003771656370000115
表示DD域雷达信道向量，hDD的第k+Nl个元素为h[k，l]，
Figure BDA0003771656370000116
表示DD域噪声向量，wDD的第k+Nl个元素为w[k，l]，符号矩阵
Figure BDA0003771656370000117
的第i行第j列元素如以下公式(17)所示：
Figure BDA0003771656370000118
其中，i和j分别满足0≤i＝k+Nl≤MN-1，0≤j＝k′+Nl′≤MN-1；
步骤4.3.2，根据步骤4.3.1.所述的矩阵形式输入输出关系，基于MMSE准则进行DD域雷达信道估计，DD域雷达信道向量估计值如以下公式(18)所示：
Figure BDA0003771656370000119
其中，
Figure BDA00037716563700001110
表示基于MMSE准则估计的DD域雷达信道向量，(·)-1表示矩阵求逆，(·)H表示复共轭转置，Ps和σ2分别表示发送的OTFS通信感知一体化信号的平均功率和噪声的平均功率，IMN表示维度为MN×MN的单位矩阵；
步骤4.3.3，根据步骤4.3.2中基于MMSE准则估计的DD域雷达信道向量获取感知目标在DD域平面栅格Γ上的时延域索引和多普勒索引的估计结果；
其中，
Figure BDA00037716563700001111
表示估计的DD域雷达信道离散响应，即多普勒域索引k和时延域索引l处的复信道增益估计值，
Figure BDA00037716563700001112
等于DD域雷达信道向量估计值
Figure BDA00037716563700001113
的第k+Nl个元素，估计的DD域雷达信道离散响应其峰值对应的多普勒域索引和时延域索引分别为感知目标在DD域平面栅格Γ上的多普勒域索引和时延域索引的估计结果，分别表示为
Figure BDA0003771656370000121
和
Figure BDA0003771656370000122
步骤4.3.4，由感知目标在DD域平面栅格Γ上的多普勒域索引和时延域索引的估计结果进行目标参数估计，估计结果如以下公式(19)、公式(20)所示：
Figure BDA0003771656370000123
Figure BDA0003771656370000124
以上公式(19)、公式(20)中，
Figure BDA0003771656370000125
和
Figure BDA0003771656370000126
分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度的参数估计结果，
Figure BDA0003771656370000127
和
Figure BDA0003771656370000128
分别表示通信感知一体化发射机与第p个目标的往返时延和多普勒频移的参数估计结果，
Figure BDA0003771656370000129
满足
Figure BDA00037716563700001210
Figure BDA00037716563700001211
满足
Figure BDA00037716563700001212
Figure BDA00037716563700001213
和
Figure BDA00037716563700001214
分别表示第p个目标在DD域平面栅格Γ上的时延域索引和多普勒域索引的估计结果；
步骤5，在DD域平面栅格Γ上采用DD域能量质心谱校正法对步骤4中OTFS通信感知一体化信号目标参数估计的结果进行校准：
步骤5.1，对步骤4.3.3中估计的DD域雷达信道离散响应进行样值比选，分别实施行比选和列比选，根据比选结果求向量元素中的最大值,即DD域雷达信道离散响应的峰值；
行比选表示按行比较选取每一行的最大值，行比选的结果为向量
Figure BDA00037716563700001215
满足
Figure BDA00037716563700001216
Figure BDA00037716563700001217
表示DD域雷达信道离散响应
Figure BDA00037716563700001218
中第k行的最大值，k为多普勒域索引，列比选表示按列比较选取每一列的最大值，列比选的结果为向量
Figure BDA00037716563700001219
满足
Figure BDA00037716563700001220
Figure BDA00037716563700001221
表示DD域雷达信道离散响应
Figure BDA00037716563700001222
中第l列的最大值，l为时延域索引，峰值
Figure BDA00037716563700001223
对应于时延域索引
Figure BDA00037716563700001224
和多普勒域索引
Figure BDA00037716563700001225
步骤5.2，以DD域雷达信道离散响应的峰值为中心，向上取X个样值、向下取个X样值，同时向左取Y个样值、向右取Y个样值，并取样值的时延域索引与多普勒域索引；
向上取X个样值的结果为
Figure BDA00037716563700001226
下取X个样值的结果为
Figure BDA00037716563700001227
左取Y个样值的结果为
Figure BDA00037716563700001228
右取Y个样值的结果为
Figure BDA00037716563700001229
步骤5.3，根据DD域样值的总能量与样值在时延域维度、多普勒域维度的加权能量，计算时延域维度的质心，并行计算多普勒域维度的质心；
步骤5.3.1，计算步骤5.2中DD域样值的总能量，并行计算样值在时延域维度的加权能量及多普勒域维度的加权能量，计算公式如以下公式(21)至公式(23)所示：
Figure BDA0003771656370000131
Figure BDA0003771656370000132
以上公式(21)至公式(23)中，E表示DD域样值的总能量，El表示时延域维度的加权能量，Ek表示多普勒域维度的加权能量，
Figure BDA0003771656370000133
表示多普勒域索引k和时延域索引l处的能量，
集合
Figure BDA0003771656370000134
表示由步骤5.2中峰值与旁瓣值对应的时延域索引、多普勒域索引构成的集合，即
Figure BDA0003771656370000135
步骤5.3.2，计算时延域维度的质心，并行计算多普勒域维度的质心，计算公式如以下公式(24)、公式(25)所示：
Figure BDA0003771656370000136
Figure BDA0003771656370000137
其中，
Figure BDA0003771656370000138
表示时延域维度的质心，
Figure BDA0003771656370000139
表示多普勒域维度的质心；
步骤5.4，根据步骤5.3得到的时延域维度的质心与多普勒域维度的质心，对目标参数估计结果进行校准；
校准结果如以下公式(26)、公式(27)所示(26)、(27)：
Figure BDA00037716563700001310
Figure BDA00037716563700001311
其中，
Figure BDA00037716563700001312
和
Figure BDA00037716563700001313
分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度在经过DD域能量质心谱校正法校准后得到的参数估计结果，
Figure BDA00037716563700001314
和
Figure BDA00037716563700001315
分别表示经过DD域能量质心谱校正法校准后得到的时延域维度的质心和多普勒域维度的质心。
为了验证本实施例中本发明所述方法的技术效果，采用计算机仿真的方式，从归一化距离像和归一化速度像的角度来评估本实施例所述方法的性能。在仿真中，雷达信道模型中的目标个数为1，即单目标，感知目标的实际距离为455m，实际速度为86m/s；射频载波频率为28GHz；带宽为5MHz；OTFS子载波数量和符号数量分别为128和64；子载波间隔为39kHz；数字调制方式采用4进制的正交幅度调制(4QAM)；距离分辨率为30m，速度分辨率为3.27m/s。图3中归一化感知目标距离像的仿真结果显示，当OTFS通信感知一体化信号在采用MMSE准则进行目标参数估计时，未使用DD域能量质心谱校正方法校准前目标距离参数估计的结果为450m，而目标实际距离为455m，在DD域能量质心谱校正方法的作用下距离参数估计的结果被校准为454.09m。因此当目标实际距离为距离分辨率的分数倍时，即目标实际距离455m是距离分辨率30m的分数倍，本发明所述方法会将其距离感知为最近的整数时延域索引对应的距离即450m，由此带来距离感知误差；在此基础上，DD域能量质心谱校正法能够对距离参数估计的结果进行进一步地校准，在距离分辨率有限且固定时能够有效降低距离参数估计误差，提升感知精度。从图4能够看出，仿真得到的归一化感知目标速度像显示在未校准前目标速度参数估计的结果为85.01m/s，其对应整数多普勒域索引，而目标实际速度为86m/s，是速度分辨率的分数倍，DD域能量质心谱校正法将速度参数估计的结果校准为85.76m/s。因此，本发明所述方法提出的OTFS通信感知一体化信号目标参数估计及DD域能量质心谱校正法方案能够有效地解决目标参数估计的技术问题且DD域能量质心谱校正法能够大大降低目标感知时的参数估计误差。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明公开的范围内，能够轻易想到的变化或替换，都应涵盖在本发明权利要求的保护范围内。

Claims (10)

1.一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，包括以下步骤：
步骤1，通信感知一体化发射机对通信信息序列进行星座映射、OTFS调制与射频调制，生成由发射天线发送的OTFS通信感知一体化信号；
步骤2，通信感知一体化发射机发送OTFS通信感知一体信号进行目标感知，在一体化信号感知到目标后，由目标反射回波信号至通信感知一体化接收机；
步骤3，通信感知一体化接收机在接收到由目标反射的回波信号后，进行下变频处理与OTFS解调，获得DD域接收调制符号矩阵；
步骤4，通信感知一体化接收机根据DD域接收调制符号与DD域发送调制符号进行雷达感知信号处理，完成雷达信道与目标参数估计；
步骤5，在DD域平面栅格Γ上采用DD域能量质心谱校正法对步骤4中OTFS通信感知一体化信号目标参数估计的结果进行校准。
2.根据权利要求1所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤1包括以下步骤：
步骤1.1，通信信源生成由二进制比特构成的数据信息序列，将信息序列输入至星座映射器并采用正交振幅调制(Quadrature Amplitude Modulation，QAM)进行数字调制获得发送调制符号向量；
Figure FDA0003771656360000011
表示长度为MN的发送调制符号向量，M表示OTFS帧的子载波数量，N表示OTFS帧的符号数量；
步骤1.2，将待传输的发送调制符号向量排列成N行M列的符号矩阵并摆放于二维时延-多普勒(Delay-Doppler，DD)域平面栅格Γ获得DD域发送调制符号矩阵
Figure FDA0003771656360000012
DD域平面栅格定义为Γ如以下公式(1)所示：
Figure FDA0003771656360000013
以上公式(1)中，k和l分别表示多普勒域的索引和时延域的索引，T和Δf分别表示OTFS符号周期和子载波间隔，并满足TΔf＝1以保持正交性，(k)N表示对k取多普勒频移运算，如以下公式(2)所示：
Figure FDA0003771656360000014
步骤1.3，对DD域发送调制符号进行OTFS调制，获得时域上的基带发射信号：
步骤1.4，对时域上的基带发射信号进行射频调制，产生的由发射天线发送的OTFS通信感知一体化信号如以下公式(6)所示：
Figure FDA0003771656360000015
其中，sRF(t)表示由发射天线发送的OTFS通信感知一体化信号，fc表示射频载波频率。
3.根据权利要求1所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤2包括以下步骤：
步骤2.1，根据所感知目标的特征建立DD域雷达信道脉冲响应的模型；
DD域雷达信道脉冲响应的模型如以下公式(7)所示：
Figure FDA0003771656360000021
其中，h(τ，v)表示DD域雷达信道，τ和v分别表示时延和多普勒频移，P表示感知目标的数量，hp表示第p个目标的复信道增益，p＝1，…，P，τp和vp分别表示通信感知一体化发射机与第p个目标的往返时延和多普勒频移，τp满足τp＝2Rp/c，vp满足vp＝2fcVp/c，Rp和Vp分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度，c代表光速，δ(·)为狄利克雷函数；
步骤2.2，将发送的OTFS通信感知一体化信号经感知目标反射后，产生回波并返回至通信感知一体化接收机的这一过程，等效于发送的OTFS通信感知一体化信号经过步骤2.1所述DD域雷达信道的射频等效信道到达通信感知一体化接收机处，通信感知一体化接收机处收到回波信号的等效过程如以下公式(8)所示：
rRF(t)＝∫∫hRF(τ，v)sRF(t-τ)ej2πv(t-τ)dτdv......(8)，
其中，rRF(t)表示通信感知一体化接收机处收到的回波信号，hRF(τ，ν)表示步骤2.1所述的DD域雷达信道的射频等效信道；其中，DD域雷达信道的射频等效信道如以下公式(9)所示：
Figure FDA0003771656360000022
4.根据权利要求1所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤3包括以下步骤：
步骤3.1，通信感知一体化接收机对回波信号进行下变频处理，获取时域基带接收信号；时域基带接收信号如以下公式(10)所示：
Figure FDA0003771656360000023
其中，r(t)表示通信感知一体化接收机处的时域基带接收信号；
步骤3.2，对时域基带接收信号进行OTFS解调，获得DD域接收调制符号矩阵
Figure FDA0003771656360000024
5.根据权利要求1所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤4包括以下步骤：
步骤4.1，将步骤2.1中所述的DD域雷达信道使用离散形式表达；DD域雷达信道使用离散形式的表达式如以下公式(13)所示：
Figure FDA0003771656360000025
其中，h[k，l]表示DD域雷达信道离散响应，即多普勒域索引k和时延域索引l处的复信道增益值，多普勒域索引k和时延域索引l分别对应于多普勒频移v＝(k)N/(NT)和时延τ＝l/(MΔf)，1/(NT)和1/(MΔf)分别表示多普勒分辨率和时延分辨率；
其中，第p个目标的往返时延τp与时延分辨率的关系为τp＝lp/(MΔf)，多普勒频移vp与多普勒分辨率的关系为vp＝kp/(NT)且关系表示第p个目标位于DD域平面栅格Γ上第lp个时延域索引和第kp个多普勒域索引处，当k＝kp且l＝lp时则h[k，l]＝hp，当在第l个时延域索引和第k个多普勒域索引处未感知到目标时则h[k，l]＝0；
步骤4.2，根据步骤4.1所述的离散形式DD域雷达信道，获得DD域发送调制符号与DD域接收调制符号之间的输入输出关系，DD域发送调制符号与DD域接收调制符号之间的输入输出关系如以下公式(14)所示：
Figure FDA0003771656360000031
其中，<·>N和<·>M分别表示模N和模M运算，w[k，l]表示DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的加性高斯白噪声，服从均值为零且方差(即平均功率)为σ2的复高斯分布，α[k，l]表示由时域样本的循环移位引起的多普勒域线性相移，α[k，l]如以下公式(15)所示：
Figure FDA0003771656360000032
步骤4.3，根据DD域接收调制符号和DD域发送调制符号进行基于最小均方误差准则(Minimum Mean Square Error，MMSE)的DD域雷达信道与目标参数估计。
6.根据权利要求1所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤5包括以下步骤：
步骤5.1，对步骤4.3.3中估计的DD域雷达信道离散响应进行样值比选，分别实施行比选和列比选，根据比选结果求向量元素中的最大值，即DD域雷达信道离散响应的峰值；
行比选表示按行比较选取每一行的最大值，行比选的结果为向量
Figure FDA0003771656360000033
满足
Figure FDA0003771656360000034
Figure FDA0003771656360000035
表示DD域雷达信道离散响应
Figure FDA0003771656360000036
中第k行的最大值，k为多普勒域索引，列比选表示按列比较选取每一列的最大值，列比选的结果为向量
Figure FDA0003771656360000037
满足
Figure FDA0003771656360000038
Figure FDA0003771656360000039
表示DD域雷达信道离散响应
Figure FDA00037716563600000310
中第l列的最大值，l为时延域索引，峰值
Figure FDA00037716563600000311
对应于时延域索引
Figure FDA00037716563600000312
和多普勒域索引
Figure FDA00037716563600000313
步骤5.2，以DD域雷达信道离散响应的峰值为中心，向上取X个样值、向下取个X样值，同时向左取Y个样值、向右取Y个样值，并取样值的时延域索引与多普勒域索引；
向上取X个样值的结果为
Figure FDA0003771656360000041
x＝1，…，X，下取X个样值的结果为
Figure FDA0003771656360000042
左取Y个样值的结果为
Figure FDA0003771656360000043
y＝1，…，Y，右取Y个样值的结果为
Figure FDA0003771656360000044
步骤5.3，根据DD域样值的总能量与样值在时延域维度、多普勒域维度的加权能量，计算时延域维度的质心，并行计算多普勒域维度的质心；
步骤5.4，根据步骤5.3得到的时延域维度的质心与多普勒域维度的质心，对目标参数估计结果进行校准；
校准结果如以下公式(26)、公式(27)所示(26)、(27)：
Figure FDA0003771656360000045
Figure FDA0003771656360000046
其中，
Figure FDA0003771656360000047
和
Figure FDA0003771656360000048
分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度在经过DD域能量质心谱校正法校准后得到的参数估计结果，
Figure FDA0003771656360000049
和
Figure FDA00037716563600000410
分别表示经过DD域能量质心谱校正法校准后得到的时延域维度的质心和多普勒域维度的质心。
7.根据权利要求2所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤1.3包括以下步骤：
步骤1.3.1，通过辛有限傅里叶逆变换(Inverse Symplectic Finite FourierTransform，ISFFT)将DD域发送调制符号矩阵映射为时间-频率(Time-Frequency，TF)域发送调制符号矩阵
Figure FDA00037716563600000411
占据TF域平面栅格Λ上的N个符号和M个子载波对应的时频格点，映射过程如以下公式(3)所示：
Figure FDA00037716563600000412
其中，ISFFT(·)表示应用辛有限傅里叶逆变换，{XTF[n，m]∈XTF，n＝0，…，N-1，m＝0，…，M-1}表示TF域平面栅格Λ上的第n个时间和第m个频率格点处的发送调制符号，DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的发送调制符号表示为{XDD[k，l]∈XDD，k＝0，…，N-1，l＝0，…，M-1}；
TF域平面栅格定义为Λ，如以下公式(4)所示：
Λ＝{(nT，mΔf)，n＝0，1，...，N-1，m＝0，1，...，M-1}......(4)，
其中，n和m分别为时间域的索引和频率域的索引；
步骤1.3.2，对TF域发送调制符号矩阵进行海森堡变换(Heisenberg Transform)，通过矩形脉冲成形滤波器获得时域基带发射信号，变换公式如以下公式(5)所示：
Figure FDA0003771656360000051
其中，s(t)表示时域基带发射信号，gtx(t)为发送端持续时间为T的矩形脉冲成形滤波器，当0≤t≤T时有gtx(t)＝1，否则gtx(t)＝0。
8.根据权利要求4所述的一种0TFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤3.2包括以下步骤：
步骤3.2.1，对时域基带接收信号应用魏格纳变换(Wigner Transform)，得到TF域接收调制符号矩阵
Figure FDA0003771656360000052
变换公式如以下公式(11)所示：
Figure FDA0003771656360000053
其中，{YTF[n，m]∈YTF，n＝0，…，N-1，m＝0，…，M-1}表示TF域平面栅格Λ上的第n个时间和第m个频率格点处的接收调制符号，grx(t)为接收端的矩形脉冲成形滤波器，(·)*表示求复共轭函数；
步骤3.2.2，采用辛有限傅里叶变换将TF域接收调制符号矩阵反映射，得到DD域接收调制符号矩阵
Figure FDA0003771656360000054
映射过程如以下公式(12)所示：
Figure FDA0003771656360000055
其中，SFFT(·)表示应用辛有限傅里叶变换，{YDD[k，l]∈YDD，k＝0，…，N-1，l＝0，…，M-1}表示DD域平面栅格Γ上的第k个多普勒和第l个时延格点处的接收调制符号。
9.根据权利要求5所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤4.3包括以下步骤：
步骤4.3.1，将步骤4.2所述的输入输出关系使用矩阵形式表达式如以下公式(16)所示：
Figure FDA0003771656360000056
其中，
Figure FDA0003771656360000057
表示DD域接收调制符号向量，yDD的第k+Nl个元素为YDD[k，l]，且0≤k+Nl≤MN-1，
Figure FDA0003771656360000058
表示DD域雷达信道向量，hDD的第k+Nl个元素为h[k，l]，
Figure FDA0003771656360000059
表示DD域噪声向量，wDD的第k+Nl个元素为w[k，l]，符号矩阵
Figure FDA00037716563600000510
的第i行第j列元素如以下公式(17)所示：
Figure FDA0003771656360000061
其中，i和j分别满足0≤i＝k+Nl≤MN-1，0≤j＝k′+Nl′≤MN-1；
步骤4.3.2，根据步骤4.3.1.所述的矩阵形式输入输出关系，基于MMSE准则进行DD域雷达信道估计，DD域雷达信道向量估计值如以下公式(18)所示：
Figure FDA0003771656360000062
其中，
Figure FDA0003771656360000063
表示基于MMSE准则估计的DD域雷达信道向量，(·)-1表示矩阵求逆，(·)H表示复共轭转置，Ps和σ2分别表示发送的OTFS通信感知一体化信号的平均功率和噪声的平均功率，IMN表示维度为MN×MN的单位矩阵；
步骤4.3.3，根据步骤4.3.2中基于MMSE准则估计的DD域雷达信道向量获取感知目标在DD域平面栅格Γ上的时延域索引和多普勒索引的估计结果；
其中，
Figure FDA0003771656360000064
表示估计的DD域雷达信道离散响应，即多普勒域索引k和时延域索引l处的复信道增益估计值，
Figure FDA0003771656360000065
等于DD域雷达信道向量估计值
Figure FDA0003771656360000066
的第k+Nl个元素，估计的DD域雷达信道离散响应其峰值对应的多普勒域索引和时延域索引分别为感知目标在DD域平面栅格Γ上的多普勒域索引和时延域索引的估计结果，分别表示为
Figure FDA0003771656360000067
和
Figure FDA0003771656360000068
步骤4.3.4，由感知目标在DD域平面栅格Γ上的多普勒域索引和时延域索引的估计结果进行目标参数估计，估计结果如以下公式(19)、公式(20)所示：
Figure FDA0003771656360000069
Figure FDA00037716563600000610
以上公式(19)、公式(20)中，
Figure FDA00037716563600000611
和
Figure FDA00037716563600000612
分别表示通信感知一体化发射机与第p个目标的相对距离和相对运动速度的参数估计结果，
Figure FDA00037716563600000613
和
Figure FDA00037716563600000614
分别表示通信感知一体化发射机与第p个目标的往返时延和多普勒频移的参数估计结果，
Figure FDA00037716563600000615
满足
Figure FDA00037716563600000616
Figure FDA00037716563600000617
满足
Figure FDA00037716563600000618
Figure FDA00037716563600000619
和
Figure FDA00037716563600000620
分别表示第p个目标在DD域平面栅格Γ上的时延域索引和多普勒域索引的估计结果。
10.根据权利要求6所述的一种OTFS通信感知一体化信号目标参数估计方法，其特征在于，所述步骤5.3包括以下步骤：
步骤5.3.1，计算步骤5.2中DD域样值的总能量，并行计算样值在时延域维度的加权能量及多普勒域维度的加权能量，计算公式如以下公式(21)至公式(23)所示：
Figure FDA0003771656360000071
Figure FDA0003771656360000072
以上公式(21)至公式(23)中，E表示DD域样值的总能量，El表示时延域维度的加权能量，Ek表示多普勒域维度的加权能量，
Figure FDA0003771656360000073
表示多普勒域索引k和时延域索引l处的能量，
集合
Figure FDA0003771656360000074
表示由步骤5.2中峰值与旁瓣值对应的时延域索引、多普勒域索引构成的集合，即
Figure FDA0003771656360000075
步骤5.3.2，计算时延域维度的质心，并行计算多普勒域维度的质心，计算公式如以下公式(24)、公式(25)所示：
Figure FDA0003771656360000076
Figure FDA0003771656360000077
其中，
Figure FDA0003771656360000078
表示时延域维度的质心，
Figure FDA0003771656360000079
表示多普勒域维度的质心。
CN202210903235.6A 2022-07-29 2022-07-29 一种otfs通信感知一体化信号目标参数估计方法 Active CN115396263B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210903235.6A CN115396263B (zh) 2022-07-29 2022-07-29 一种otfs通信感知一体化信号目标参数估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210903235.6A CN115396263B (zh) 2022-07-29 2022-07-29 一种otfs通信感知一体化信号目标参数估计方法

Publications (2)

Publication Number Publication Date
CN115396263A true CN115396263A (zh) 2022-11-25
CN115396263B CN115396263B (zh) 2023-06-30

Family

ID=84116405

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210903235.6A Active CN115396263B (zh) 2022-07-29 2022-07-29 一种otfs通信感知一体化信号目标参数估计方法

Country Status (1)

Country Link
CN (1) CN115396263B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011054279A1 (zh) * 2009-11-05 2011-05-12 Gao Xiqi 基于多相分解的多天线信道估计方法
US20150117395A1 (en) * 2010-05-28 2015-04-30 Cohere Technologies, Inc. Otfs methods of data channel characterization and uses thereof
CN112751798A (zh) * 2021-03-05 2021-05-04 上海交通大学 用于太赫兹通信感知一体化的nmw-ofdm系统
EP3826253A1 (en) * 2019-11-19 2021-05-26 Volkswagen AG Otfs embedded pilot estimation extension
CN112882005A (zh) * 2021-01-14 2021-06-01 东南大学 一种基于贝叶斯学习的otfs雷达目标参数估计方法
CN113507426A (zh) * 2021-06-16 2021-10-15 北京邮电大学 基于otfs调制的联合信道估计与信号检测方法及装置
WO2021254636A1 (en) * 2020-06-19 2021-12-23 NEC Laboratories Europe GmbH User localization in a wireless cellular network
CN114124238A (zh) * 2021-10-20 2022-03-01 北京邮电大学 基于时分体制的otfs通信雷达一体化波形设计方法
CN114785644A (zh) * 2022-05-19 2022-07-22 江苏科技大学 一种移动水声otfs通信稀疏信道估计方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011054279A1 (zh) * 2009-11-05 2011-05-12 Gao Xiqi 基于多相分解的多天线信道估计方法
US20150117395A1 (en) * 2010-05-28 2015-04-30 Cohere Technologies, Inc. Otfs methods of data channel characterization and uses thereof
EP3826253A1 (en) * 2019-11-19 2021-05-26 Volkswagen AG Otfs embedded pilot estimation extension
WO2021254636A1 (en) * 2020-06-19 2021-12-23 NEC Laboratories Europe GmbH User localization in a wireless cellular network
CN112882005A (zh) * 2021-01-14 2021-06-01 东南大学 一种基于贝叶斯学习的otfs雷达目标参数估计方法
CN112751798A (zh) * 2021-03-05 2021-05-04 上海交通大学 用于太赫兹通信感知一体化的nmw-ofdm系统
CN113507426A (zh) * 2021-06-16 2021-10-15 北京邮电大学 基于otfs调制的联合信道估计与信号检测方法及装置
CN114124238A (zh) * 2021-10-20 2022-03-01 北京邮电大学 基于时分体制的otfs通信雷达一体化波形设计方法
CN114785644A (zh) * 2022-05-19 2022-07-22 江苏科技大学 一种移动水声otfs通信稀疏信道估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
C. YOU: "Channel estimation and passive beamforming for intelligent reflecting surface: Discrete phase shift and progressive refinement", 《IEEE J. SELECT. AREAS COMMUN.,》 *
SHUANGYANG LI等: "A Novel ISAC Transmission Framework Based on Spatially-Spread Orthogonal Time Frequency Space Modulation", 《IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS》 *
闫实等: "通信-感知-计算融合:6G愿景与关键技术", 《北京邮电大学学报》 *
黄瑞;杜小勇;胡卫东;: "OFDM雷达多目标运动参数的近似最大似然估计", 雷达学报, no. 04 *

Also Published As

Publication number Publication date
CN115396263B (zh) 2023-06-30

Similar Documents

Publication Publication Date Title
Schweizer et al. Stepped-carrier OFDM-radar processing scheme to retrieve high-resolution range-velocity profile at low sampling rate
Wei et al. Integrated sensing and communication signals toward 5G-A and 6G: A survey
Liu et al. Joint range and angle estimation for an integrated system combining MIMO radar with OFDM communication
Xu et al. A bandwidth efficient dual-function radar communication system based on a MIMO radar using OFDM waveforms
JP6586129B2 (ja) Ｍｉｍｏ方式システムの試験装置および試験方法
CN113765581A (zh) 基于压缩感知与波束对齐的ris快时变信道估计方法
CN113612716B (zh) 基于恒定包络正交频分复用的多用户通信雷达一体化系统
CN104486266A (zh) 一种基于mimo-ofdm系统的信道估计方法及装置
CN113542162A (zh) 基于块稀疏贝叶斯算法的上下链路通信感知一体化方法
Han et al. THz ISAC: A physical-layer perspective of terahertz integrated sensing and communication
CN112968722A (zh) 基于fd-mimo的雷达通信一体化信号设计方法
Zhang et al. Radar sensing via otfs signaling: A delay doppler signal processing perspective
Wu et al. Accurate channel estimation for frequency-hopping dual-function radar communications
CN111641575B (zh) 一种正交时频二维空间调制信号接收方法及接收器
CN115396263B (zh) 一种otfs通信感知一体化信号目标参数估计方法
Wu et al. A non-uniform multi-wideband OFDM system for terahertz joint communication and sensing
Delamou et al. An Efficient OFDM-Based Monostatic Radar Design for Multitarget Detection
CN1351456A (zh) 获得一个传输增益函数的方法
CN116047503A (zh) 基于时延-多普勒联合调制的分布式雷达通信一体化方法
CN116405348A (zh) 通信感知一体化系统阵列天线与5g新空口估角方法
Xu et al. A dual-function radar communication system with OFDM waveforms and subcarrier sharing
Xiao et al. A novel joint angle-range-velocity estimation method for MIMO-OFDM ISAC systems
Li et al. A novel design of multi-user sequence set for joint vehicular radar-communication based on Oppermann family
CN114884787A (zh) 适用快时变信道的波形可控多载波通信方法、装置及系统
Jacome et al. Multi-antenna dual-blind deconvolution for joint radar-communications via soman minimization

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant