CN113395221A - 一种基于正交时频空联合信道估计与符号检测方法 - Google Patents

一种基于正交时频空联合信道估计与符号检测方法 Download PDF

Info

Publication number
CN113395221A
CN113395221A CN202110448874.3A CN202110448874A CN113395221A CN 113395221 A CN113395221 A CN 113395221A CN 202110448874 A CN202110448874 A CN 202110448874A CN 113395221 A CN113395221 A CN 113395221A
Authority
CN
China
Prior art keywords
symbol
matrix
channel
channel estimation
signal
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
CN202110448874.3A
Other languages
English (en)
Other versions
CN113395221B (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 CN202110448874.3A priority Critical patent/CN113395221B/zh
Publication of CN113395221A publication Critical patent/CN113395221A/zh
Application granted granted Critical
Publication of CN113395221B publication Critical patent/CN113395221B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0204Channel estimation of multiple channels
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/0224Channel estimation using sounding signals
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明公开了一种基于正交时频空联合信道估计与符号检测方法,包括以下步骤:步骤1,导频符号放置:根据MIMO‑OTFS系统输入与输出之间的关系建立数学模型,通过分析数学模型将信道估计问题转换为稀疏信号重构问题;步骤2,信道估计,在接收端,信号先经过预处理后得到时延‑多普勒域符号;步骤3,符号检测,在接收端,根据信道估计步骤所得到的信道参数进行符号检测;步骤4,迭代反馈,本发明的优越效果是:基于迭代联合的思想提出了多输入多输出正交时频空系统下的联合信道估计与符号检测方法,能够通过不断的反馈迭代来提高系统性能,对低于阈值概率的符号索引集合重新进行检测,形成信号检测与信道估计的联合迭代检测来提升系统性能。

Description

一种基于正交时频空联合信道估计与符号检测方法
技术领域
本发明属于高移动场景下无线通信技术领域,尤其涉及一种基于多输入多输出正交时频空(MIMO-OTFS)调制系统的联合信道估计与符号检测方法。
背景技术
目前,现有的5G通信技术采用的调制方式是正交频分复用(OFDM)方案,这种方案对频偏十分的敏感,而高速运动会产生较高的多普勒频移,会引起载波间干扰(Inter-Carrier Interference,ICI),从而导致子载波间的正交性被破坏,因此系统性能大幅衰减。Ronny Hadani等人近年来提出了正交时频空(OTFS)调制方案,技术在双选择性衰落信道中比OFDM具有显著的性能优势,尤其在对多普勒偏移具有较高的容忍性,原因在于OTFS调制方案是在时延-多普勒域上展开多路复用的,这样能够将时频域上的双选信道转换为时延-多普勒域上的恒不变信道,从而使得所有的时延-多普勒符号都经历相同的信道增益。由于OTFS调制优点,研究人员将其与多输入多输出(MIMO)系统相结合来实现高移动性场景下的高速率通信需求,然而,关于MIMO-OTFS系统的研究存在两大挑战,即信道估计与信号检测。对于OTFS系统的信道估计问题,许多现有方案将不适用于时延-多普勒域。同时,为保证高速移动的通信终端的性能,信号处理时间将大大减少,这使得盲信道估计方案也将不适用于高移动性场景。因此,OTFS调制方案的信道估计方案必须具备低复杂度,高准确度,低导频开销等特性。近年来的关于OTFS系统的信道估计方法分为线性和非线性这两类。其中基于嵌入式导频的信道估计方案就属于线性信道估计,先在时延-多普勒域中插入一个高功率的导频符号,并在其周围填充功率为零的符号来避免数据符号的干扰,然后在接收端通过阈值判定出时延-多普勒域信道抽头的非零值位置,最后利用这些位置所对应的接收信号与导频符号相除,从而得到信道的估计值。然而,这类信道估计方案准确度较低,而且导频功率至少高于符号功率20dB才有效,进而带来高峰均功率比这一缺陷。非线性信道估计方案则利用时延-多普勒域信道的稀疏特性,通过压缩感知,稀疏贝叶斯学习以及消息传递等理论来重构信道参数。但是非线性信道估计方案地缺点在于过高的计算复杂度。
对于OTFS系统的信号检测问题,最大的挑战在于系统性能与计算复杂度之间的权衡。根据输入信号与输出信号在时延多普勒域上的关系,每个接收信号都受到多个发送信号的影响,这使得最优的最大似然检测方法将不现实。目前信号检测方法也分为线性与非线性两类。其中线性方案就包括迫零(ZF)以及最小均方误差(MMSE)均衡,这类信号检测方案利用矩阵特性从而简化求解信道矩阵。线性信号均衡方案虽然降低了计算复杂度,但是也存在系统性能较低的缺点。非线性信号检测方案利用收发端间的关系,从而通过迭代的思想来逐步逼近最优解,其包括消息传递、变分贝叶斯以及反馈均衡等,这些方案极大地提高了系统性能,但是迭代使其计算复杂度也特别高。
发明内容
本发明的目的在于提供一种能够克服上述技术问题的基于正交时频空(MIMO-OTFS)联合信道估计与符号检测方法,本发明所述方法适用于正交时频空(MIMO-OTFS)系统的信道估计与符号检测。
本发明所述方法包括以下步骤:
步骤1,导频符号放置:根据MIMO-OTFS系统输入与输出之间的关系建立数学模型,通过分析数学模型将信道估计转换为稀疏信号重构,在时延-多普勒域放置少量的导频符号用于接收端的信道参数估计,其中收发天线数目Nt=Nr=Na,首先输入比特先经过编码,交织以及正交相移键控(QPSK)调制映射后得到时延-多普勒域符号,根据输入与输出之间的关系导频符号放置图案,设定系统的信道模型为以下公式(1):
Figure BDA0003038009300000021
公式(1)中,
Figure BDA0003038009300000022
分别表示第nt发送天线与第nr接收天线之间信号在第i径的时延以及多普勒,同理,
Figure BDA0003038009300000023
为信道复增益,则输入与输出之间的关系为以下公式(2)、(3):
Figure BDA0003038009300000024
其中:
Figure BDA0003038009300000025
Figure BDA0003038009300000026
公式(1)中,
Figure BDA0003038009300000027
表示第nr个接收天线的第k行,l列时延-多普勒域符号,M和N分别表示子载波和符号数目,
Figure BDA0003038009300000028
为第nr接收天线的系统噪声,其中,P表示信道多径数目,Q=5为分数多普勒近似参数,
Figure BDA0003038009300000029
分别表示第nt发送天线与第nr接收天线之间信号在第i径的时延以及多普勒抽头,同理,
Figure BDA00030380093000000210
为信道复增益,公式(2)表示相位偏移,其中,
Figure BDA00030380093000000211
是多普勒的分数抽头,取值范围为(-0.5,0.5],
Figure BDA00030380093000000212
代表分数多普勒抽头干扰项,根据公式(1)一个接收符号
Figure BDA00030380093000000213
由NaP(2Q+1)个发送信号
Figure BDA00030380093000000214
叠加以及高斯噪声所组成,其中kv和lτ分别是最大多普勒偏移和最大时延所对应的时延-多普勒域网格上多普勒轴以及时延轴位置;
步骤2,信道估计,在接收端,信号先经过预处理后得到时延-多普勒域符号,提取出导频符号并构建出求解稀疏信号恢复问题的观测矩阵,再通过压缩感知理论或者稀疏贝叶斯学习理论利用观测矩阵对信道进行估计得到信道参数,信号先经过预处理得到时延-多普勒域符号,根据时延-多普勒域信道的稀疏特性将信道估计转换为稀疏信号恢复,则接收端第p个导频位置处的公式(2)能够写为向量形式,如以下公式(4)所示:
Figure BDA0003038009300000031
其中:
Figure BDA0003038009300000032
分别为接收端导频信号和噪声,
Figure BDA0003038009300000033
为发送端导频信号,L=(2kv+2Q+1)(Nala+1)是待估计的时延多普勒域信道长,
Figure BDA0003038009300000034
为公式(5)中的
Figure BDA0003038009300000035
的向量形式,已知导频信息,通过α对导频符号进行预处理以提高信道估计准确度,
Figure BDA0003038009300000036
为待重构的时延多普勒域信道矩阵,h矩阵具有对称块稀疏特性,采用块正交匹配追踪(BOMP)算法来进行信道估计,具体包括以下步骤:
步骤2.1,首先初始化残差
Figure BDA0003038009300000037
初始化索引集合
Figure BDA0003038009300000038
迭代次数n=1;
步骤2.2,残差与字典矩阵Φp中的每一列进行内积运算,并记录绝对值最大的列向量的索引值λ,计算公式为
Figure BDA0003038009300000039
步骤2.3,由索引λ求解块稀疏所对应的连续索引位置ω=λ-Q,...,λ+Q;
步骤2.4,更新索引值集合Ωn=Ωn-1∪ω,即步骤2.3中求出的块稀疏位置集合与n-1次迭代索引集合Ωn-1的并集,生成新的索引集合Ωn
步骤2.5,获得新的估计值
Figure BDA00030380093000000310
其中Ωn为被选索引值选取的列向量,上标H表示共轭转置,上标-1表示求逆矩阵,
Figure BDA00030380093000000311
为最小均方估计表达式;
步骤2.6,更新残差
Figure BDA00030380093000000312
为第n次迭代的估计值,
Figure BDA00030380093000000313
为第nr接收天线在P导频处的接收信号,Φp为字典矩阵;
步骤2.7,判断是否达到迭代停止条件n>NaP,当达到停止条件则输出信道估计矩阵
Figure BDA00030380093000000314
否则返回步骤2.2;
步骤3,符号检测,在接收端,根据信道估计步骤所得到的信道参数进行符号检测,其中,符号检测采用低复杂度近似消息传递(AMP)算法,得到每个星座符号的概率并传递给下次迭代:
步骤3.1,根据步骤2所得的信道参数进行符号检测,由公式(2)推导输入与输出之间的关系矩阵表达式,表示为以下公式(5):
yMIMO=HMIMOxMIMO+wMIMO......(5),
其中:
Figure BDA0003038009300000041
分别为时延多普勒域的接收信号、发送信号以及系统噪声,
Figure BDA0003038009300000042
为时延-多普勒域的稀疏信道矩阵,根据公式(5)得到符号的最大后验概率准则为以下公式(6):
Figure BDA0003038009300000043
其中:A是Q个调制符号,由于HMIMO矩阵维度太大,公式(6)无法直接计算,则利用HMIMO每行以及每列仅有NaP(2Q+1)个非零元素这一特性来降低计算复杂度,利用信道参数提取出每行非零元素的索引矩阵为
Figure BDA0003038009300000044
以及其对应的元素矩阵Hrow,同理得到列元素对应的索引矩阵
Figure BDA0003038009300000045
和元素矩阵Hcol,对公式(6)进行逐符号的检测,得到第p个符号检测结果为以下公式(7):
Figure BDA0003038009300000046
在消息传递(AMP)信号检测算法中,干扰项的均值和方差作为观测节点到变量节点的传递信息,第p个变量节点传递给第q个观测节点的信息为概率质量函数(pmf)Pp,q{pp,q(aj)},第p个变量节点传递给S=NaP(2Q+1)个观测节点的信息近似相等,则简化为
Figure BDA0003038009300000047
矩阵HMIMO行非零元素干扰项的均值和方差的矩阵如以下公式(8)、(9)所示:
Figure BDA0003038009300000048
Figure BDA0003038009300000049
根据矩阵
Figure BDA00030380093000000410
推导出列维非零元素的均值和方差μcol
Figure BDA00030380093000000411
以及pmf的计算公式为以下公式(10):
Figure BDA0003038009300000051
步骤3.2,通过多次迭代得到最终的星座符号概率并根据星座符号概率判定接收信号的星座位置进行解映射操作,通过矩阵方式计算干扰项的均值和方差,在pmf更新过程中,当pp(aj)>0.9时,下次迭代将不再更新第p个符号的概率,即随着迭代次数的增加,需要重新计算的符号概率数目将不断减少以降低算法的复杂度;
步骤4,迭代反馈,优选部分数据符号与导频符号构成新的导频观测矩阵,利用新的导频观测矩阵重新进行信道估计,根据所估计出的新的信道参数对低于阈值概率的符号重新进行检测以形成信号检测与信道估计的联合迭代检测来提升系统性能:
步骤4.1,当没有达到迭代停止条件,如迭代次数,则将步骤3步所得到的数据符号经过解映射、解交织以及译码操作得到最终的信息比特并再将所述信息比特进行逆操作,即编码、交织以及映射得到时延-多普勒域符号,优选部分数据符号与已知的导频符号组成新的导频观测矩阵,即字典矩阵Φp,增加观测数目以提高信道重构的准确度;
步骤4.2,对低于阈值概率的符号重新进行检测以形成信号检测与信道估计的联合迭代检测来提升系统性能,在每次迭代中的所需要计算概率的符号数目都在不断减少以使得所提方案能够便于实现。
本发明所述方法的优越效果是:
1、本发明所述方法基于迭代联合的思想提出了多输入多输出正交时频空(MIMO-OTFS)系统下的联合信道估计与符号检测方法,能够通过不断的反馈迭代来提高系统性能;
2、本发明所述方法在发射端的时延-多普勒域上放置少量导频符号,在接收端,信号先经过预处理后得到时延-多普勒域符号,提取出导频符号后根据稀疏信号恢复理论对信道估计得到信道参数;
3、本发明所述方法根据所估计出的信道参数利用近似消息传递算法(AMP)对接收信号进行符号检测得到符号概率,优选部分数据符号与导频符号组成新的导频观测矩阵并重新进行信道估计;
4、本发明所述方法对低于阈值概率的符号索引集合重新进行检测,形成信号检测与信道估计的联合迭代检测来提升系统性能。
附图说明
图1是本发明所述方法的流程图;
图2是本发明所述方法的MIMO-OTFS系统收发端流程示意图;
图3是本发明所述方法的MIMO-OTFS系统收发端导频示意图。
具体实施方式
下面结合附图对本发明的实施方式进行详细描述。本发明所述方法包括以下步骤:
步骤1,导频符号放置:根据MIMO-OTFS系统输入与输出之间的关系建立数学模型,通过分析数学模型将信道估计转换为稀疏信号重构,在时延-多普勒域放置少量的导频符号用于接收端的信道参数估计,MIMO-OTFS系统收发端如图1所示,其中收发天线数目Nt=Nr=Na,首先输入比特先经过编码,交织以及正交相移键控(QPSK)调制映射后得到时延-多普勒域符号,根据输入与输出之间的关系导频符号放置图案如图2所示,设定系统的信道模型为以下公式(1):
Figure BDA0003038009300000061
公式(1)中
Figure BDA0003038009300000062
分别表示第nt发送天线与第nr接收天线之间信号在第i径的时延以及多普勒,同理,
Figure BDA0003038009300000063
为信道复增益,则输入与输出之间的关系为以下公式(2)、(3):
Figure BDA0003038009300000064
其中:
Figure BDA0003038009300000065
Figure BDA0003038009300000066
公式(1)中
Figure BDA0003038009300000067
表示第nr个接收天线的第k行,l列时延-多普勒域符号,M和N分别表示子载波和符号数目,
Figure BDA0003038009300000068
为第nr接收天线的系统噪声,其中P表示信道多径数目,Q=5为分数多普勒近似参数,
Figure BDA0003038009300000069
分别表示第nt发送天线与第nr接收天线之间信号在第i径的时延以及多普勒抽头,同理,
Figure BDA00030380093000000610
为信道复增益,公式(2)表示相位偏移,其中
Figure BDA00030380093000000611
是多普勒的分数抽头,取值范围为(-0.5,0.5],
Figure BDA00030380093000000612
代表分数多普勒抽头干扰项,根据公式(1)一个接收符号
Figure BDA00030380093000000613
由NaP(2Q+1)个发送信号
Figure BDA00030380093000000614
叠加以及高斯噪声所组成,导频符号放置位置如图2所示,其中kv和lτ分别是最大多普勒偏移和最大时延所对应的时延-多普勒域网格上多普勒轴以及时延轴位置;
步骤2,信道估计,在接收端,信号先经过预处理后得到时延-多普勒域符号,提取出导频符号并构建出求解稀疏信号恢复问题的观测矩阵,再通过压缩感知理论或者稀疏贝叶斯学习理论利用观测矩阵对信道进行估计得到信道参数,信号先经过预处理得到时延-多普勒域符号,根据时延-多普勒域信道的稀疏特性将信道估计问题转换为稀疏信号恢复问题,则接收端第p个导频位置处的公式(2)能够写为向量形式,如以下公式(4)所示:
Figure BDA0003038009300000071
其中:
Figure BDA0003038009300000072
分别为接收端导频信号和噪声,
Figure BDA0003038009300000073
为发送端导频信号,L=(2kv+2Q+1)(Nala+1)是待估计的时延多普勒域信道长,
Figure BDA0003038009300000074
为公式(5)中的
Figure BDA0003038009300000075
的向量形式,已知导频信息,通过α对导频符号进行预处理以提高信道估计准确度,
Figure BDA0003038009300000076
为待重构的时延多普勒域信道矩阵,h矩阵具有对称块稀疏特性,采用块正交匹配追踪(BOMP)算法来进行信道估计,具体包括以下步骤:
步骤2.1,首先初始化残差
Figure BDA0003038009300000077
初始化索引集合
Figure BDA0003038009300000078
迭代次数n=1;
步骤2.2,残差与字典矩阵Φp中的每一列进行内积运算,并记录绝对值最大的列向量的索引值λ,计算公式为
Figure BDA0003038009300000079
步骤2.3,由索引λ求解块稀疏所对应的连续索引位置ω=λ-Q,...,λ+Q;
步骤2.4,更新索引值集合Ωn=Ωn-1∪ω,即步骤2.3中求出的块稀疏位置集合与n-1次迭代索引集合Ωn-1的并集,生成新的索引集合Ωn
步骤2.5,获得新的估计值
Figure BDA00030380093000000710
其中Ωn为被选索引值选取的列向量,上标H表示共轭转置,上标-1表示求逆矩阵,
Figure BDA00030380093000000711
为最小均方估计表达式;
步骤2.6,更新残差
Figure BDA00030380093000000712
为第n次迭代的估计值,
Figure BDA00030380093000000713
为第nr接收天线在P导频处的接收信号,Φp为字典矩阵;
步骤2.7,判断是否达到迭代停止条件n>NaP,当达到停止条件则输出信道估计矩阵
Figure BDA00030380093000000714
否则返回步骤2.2;
步骤3,符号检测,在接收端,根据信道估计步骤所得到的信道参数进行符号检测,其中:符号检测采用低复杂度近似消息传递(AMP)算法,得到每个星座符号的概率并传递给下次迭代:
步骤3.1,根据步骤2所得的信道参数进行符号检测,由公式(2)推导输入与输出之间的关系矩阵表达式,表示为以下公式(5):
yMIMO=HMIMOxMIMO+wMIMO......(15),
其中
Figure BDA0003038009300000081
分别为时延多普勒域的接收信号、发送信号以及系统噪声,
Figure BDA0003038009300000082
为时延-多普勒域的稀疏信道矩阵,根据公式(5)得到符号的最大后验概率准则为以下公式(6):
Figure BDA0003038009300000083
其中:A是Q个调制符号,由于HMIMO矩阵维度太大,公式(6)无法直接计算,则利用HMIMO每行以及每列仅有NaP(2Q+1)个非零元素这一特性来降低计算复杂度,利用信道参数提取出每行非零元素的索引矩阵为
Figure BDA0003038009300000084
以及其对应的元素矩阵Hrow,同理得到列元素对应的索引矩阵
Figure BDA0003038009300000085
和元素矩阵Hcol,对公式(6)进行逐符号的检测,得到第p个符号检测结果为以下公式(7):
Figure BDA0003038009300000089
在消息传递(AMP)信号检测算法中,干扰项的均值和方差作为观测节点到变量节点的传递信息,第p个变量节点传递给第q个观测节点的信息为概率质量函数(pmf)Pp,q{pp,q(aj)},第p个变量节点传递给S=NaP(2Q+1)个观测节点的信息近似相等,则简化为
Figure BDA0003038009300000086
矩阵HMIMO行非零元素干扰项的均值和方差的矩阵如以下公式(8)、(9)所示:
Figure BDA00030380093000000810
Figure BDA00030380093000000811
根据矩阵μrow
Figure BDA0003038009300000087
推导出列维非零元素的均值和方差μcol
Figure BDA0003038009300000088
以及pmf的计算公式为以下公式(10):
Figure BDA0003038009300000091
步骤3.2,通过多次迭代得到最终的星座符号概率并根据星座符号概率判定接收信号的星座位置进行解映射操作,通过矩阵方式计算干扰项的均值和方差,在pmf更新过程中,当pp(aj)>0.9时,下次迭代将不再更新第p个符号的概率,即随着迭代次数的增加,需要重新计算的符号概率数目将不断减少以降低算法的复杂度;
步骤4,迭代反馈,优选部分数据符号与导频符号构成新的导频观测矩阵,利用新的导频观测矩阵重新进行信道估计,根据所估计出的新的信道参数对低于阈值概率的符号重新进行检测以形成信号检测与信道估计的联合迭代检测来提升系统性能:
步骤4.1,迭代流程图如图3所示,当没有达到迭代停止条件,如迭代次数,则将步骤3步所得到的数据符号经过解映射、解交织以及译码操作得到最终的信息比特并再将所述信息比特进行逆操作,即编码、交织以及映射得到时延-多普勒域符号,优选部分数据符号与已知的导频符号组成新的导频观测矩阵,即字典矩阵Φp,增加观测数目以提高信道重构的准确度;
步骤4.2,对低于阈值概率的符号重新进行检测以形成信号检测与信道估计的联合迭代检测来提升系统性能,在每次迭代中的所需要计算概率的符号数目都在不断减少以使得所提方案能够便于实现。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明公开的范围内,能够轻易想到的变化或替换,都应涵盖在本发明权利要求的保护范围内。

Claims (4)

1.一种基于正交时频空联合信道估计与符号检测方法,其特征在于,包括以下步骤:
步骤1,导频符号放置:根据MIMO-OTFS系统输入与输出之间的关系建立数学模型,通过分析数学模型将信道估计转换为稀疏信号重构,在时延-多普勒域放置少量的导频符号用于接收端的信道参数估计,其中收发天线数目Nt=Nr=Na,首先输入比特先经过编码,交织以及正交相移键控调制映射后得到时延-多普勒域符号,根据输入与输出之间的关系导频符号放置图案,设定系统的信道模型为以下公式(1):
Figure FDA0003038009290000011
公式(1)中
Figure FDA0003038009290000012
分别表示第nt发送天线与第nr接收天线之间信号在第i径的时延以及多普勒,同理,
Figure FDA0003038009290000013
为信道复增益,则输入与输出之间的关系为以下公式(2)、(3):
Figure FDA0003038009290000014
其中:
Figure FDA0003038009290000015
Figure FDA0003038009290000016
公式(1)中
Figure FDA0003038009290000017
表示第nr个接收天线的第k行,l列时延-多普勒域符号,M和N分别表示子载波和符号数目,
Figure FDA0003038009290000018
为第nr接收天线的系统噪声,其中P表示信道多径数目,Q=5为分数多普勒近似参数,
Figure FDA0003038009290000019
分别表示第nt发送天线与第nr接收天线之间信号在第i径的时延以及多普勒抽头,同理,
Figure FDA00030380092900000110
为信道复增益,公式(2)表示相位偏移,其中
Figure FDA00030380092900000111
是多普勒的分数抽头,取值范围为(-0.5,0.5],
Figure FDA00030380092900000112
代表分数多普勒抽头干扰项,根据公式(1)一个接收符号
Figure FDA00030380092900000113
由NaP(2Q+1)个发送信号
Figure FDA00030380092900000114
叠加以及高斯噪声所组成,其中kν和lτ分别是最大多普勒偏移和最大时延所对应的时延-多普勒域网格上多普勒轴以及时延轴位置;
步骤2,信道估计,在接收端,信号先经过预处理后得到时延-多普勒域符号,提取出导频符号并构建出求解稀疏信号恢复问题的观测矩阵,再通过压缩感知理论或者稀疏贝叶斯学习理论利用观测矩阵对信道进行估计得到信道参数,信号先经过预处理得到时延-多普勒域符号,根据时延-多普勒域信道的稀疏特性将信道估计转换为稀疏信号恢复,则接收端第p个导频位置处的公式(2)能够写为向量形式,如以下公式(4)所示:
Figure FDA0003038009290000021
其中
Figure FDA0003038009290000022
分别为接收端导频信号和噪声,
Figure FDA0003038009290000023
为发送端导频信号,L=(2kv+2Q+1)(Nala+1)是待估计的时延多普勒域信道长,
Figure FDA0003038009290000024
为公式(5)中的
Figure FDA0003038009290000025
的向量形式,已知导频信息,通过α对导频符号进行预处理以提高信道估计准确度,
Figure FDA0003038009290000026
为待重构的时延多普勒域信道矩阵,h矩阵具有对称块稀疏特性,采用块正交匹配追踪(BOMP)算法来进行信道估计;
步骤3,符号检测,在接收端,根据信道估计步骤所得到的信道参数进行符号检测,其中符号检测采用低复杂度近似消息传递算法,得到每个星座符号的概率并传递给下次迭代;
步骤4,迭代反馈,优选部分数据符号与导频符号构成新的导频观测矩阵,利用新的导频观测矩阵重新进行信道估计,根据所估计出的新的信道参数对低于阈值概率的符号重新进行检测以形成信号检测与信道估计的联合迭代检测来提升系统性能。
2.根据权利要求1所述的一种基于正交时频空联合信道估计与符号检测方法,其特征在于,所述步骤2包括以下步骤:
步骤2.1,首先初始化残差
Figure FDA0003038009290000027
初始化索引集合
Figure FDA0003038009290000028
迭代次数n=1;
步骤2.2,残差与字典矩阵Φp中的每一列进行内积运算,并记录绝对值最大的列向量的索引值λ,计算公式为
Figure FDA0003038009290000029
步骤2.3,由索引λ求解块稀疏所对应的连续索引位置ω=λ-Q,...,λ+Q;
步骤2.4,更新索引值集合Ωn=Ωn-1∪ω,即步骤2.3中求出的块稀疏位置集合与n-1次迭代索引集合Ωn-1的并集,生成新的索引集合Ωn
步骤2.5,获得新的估计值
Figure FDA00030380092900000210
其中Ωn为被选索引值选取的列向量,上标H表示共轭转置,上标-1表示求逆矩阵,
Figure FDA00030380092900000211
为最小均方估计表达式;
步骤2.6,更新残差
Figure FDA00030380092900000212
Figure FDA00030380092900000213
为第n次迭代的估计值,
Figure FDA00030380092900000214
为第nr接收天线在P导频处的接收信号,Φp为字典矩阵;
步骤2.7,判断是否达到迭代停止条件n>NaP,当达到停止条件则输出信道估计矩阵
Figure FDA00030380092900000215
否则返回步骤2.2。
3.根据权利要求1所述的一种基于正交时频空联合信道估计与符号检测方法,其特征在于,所述步骤3包括以下步骤:
步骤3.1,根据步骤2所得的信道参数进行符号检测,由公式(2)推导输入与输出之间的关系矩阵表达式,表示为以下公式(5):
yMIMO=HMIMOxMIMO+wMIMO……(5)
其中
Figure FDA0003038009290000031
分别为时延多普勒域的接收信号、发送信号以及系统噪声,
Figure FDA0003038009290000032
为时延-多普勒域的稀疏信道矩阵,根据公式(5)得到符号的最大后验概率准则为以下公式(6):
Figure FDA0003038009290000033
其中A是Q个调制符号,由于HMIMO矩阵维度太大,公式(6)无法直接计算,则利用HMIMO每行以及每列仅有NaP(2Q+1)个非零元素这一特性来降低计算复杂度,利用信道参数提取出每行非零元素的索引矩阵为
Figure FDA0003038009290000034
以及其对应的元素矩阵Hrow,同理得到列元素对应的索引矩阵
Figure FDA0003038009290000035
和元素矩阵Hcol,对公式(6)进行逐符号的检测,得到第p个符号检测结果为以下公式(7):
Figure FDA0003038009290000036
在消息传递(AMP)信号检测算法中,干扰项的均值和方差作为观测节点到变量节点的传递信息,第p个变量节点传递给第q个观测节点的信息为概率质量函数(pmf)Pp,q{pp,q(aj)},第p个变量节点传递给S=NaP(2Q+1)个观测节点的信息近似相等,则简化为
Figure FDA0003038009290000037
矩阵HMIMO行非零元素干扰项的均值和方差的矩阵如以下公式(8)、(9)所示:
Figure FDA0003038009290000038
Figure FDA0003038009290000039
根据矩阵μrow
Figure FDA00030380092900000310
推导出列维非零元素的均值和方差μcol
Figure FDA00030380092900000311
以及pmf的计算公式为以下公式(10):
Figure FDA0003038009290000041
步骤3.2,通过多次迭代得到最终的星座符号概率并根据星座符号概率判定接收信号的星座位置进行解映射操作,即通过矩阵方式计算干扰项的均值和方差,在pmf更新过程中,当pp(aj)>0.9时,下次迭代将不再更新第p个符号的概率,即随着迭代次数的增加,需要重新计算的符号概率数目将不断减少以降低算法的复杂度。
4.根据权利要求1所述的一种基于正交时频空联合信道估计与符号检测方法,其特征在于,所述步骤4包括以下步骤:
步骤4.1,当没有达到迭代停止条件,如迭代次数,则将步骤3步所得到的数据符号经过解映射、解交织以及译码操作得到最终的信息比特并再将所述信息比特进行逆操作,即编码、交织以及映射得到时延-多普勒域符号,优选部分数据符号与已知的导频符号组成新的导频观测矩阵,即字典矩阵Φp
步骤4.2,对低于阈值概率的符号重新进行检测以形成信号检测与信道估计的联合迭代检测来提升系统性能,在每次迭代中的所需要计算概率的符号数目都在不断减少以使得所提方案能够便于实现。
CN202110448874.3A 2021-04-25 2021-04-25 一种基于正交时频空联合信道估计与符号检测方法 Active CN113395221B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110448874.3A CN113395221B (zh) 2021-04-25 2021-04-25 一种基于正交时频空联合信道估计与符号检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110448874.3A CN113395221B (zh) 2021-04-25 2021-04-25 一种基于正交时频空联合信道估计与符号检测方法

Publications (2)

Publication Number Publication Date
CN113395221A true CN113395221A (zh) 2021-09-14
CN113395221B CN113395221B (zh) 2022-07-08

Family

ID=77617635

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110448874.3A Active CN113395221B (zh) 2021-04-25 2021-04-25 一种基于正交时频空联合信道估计与符号检测方法

Country Status (1)

Country Link
CN (1) CN113395221B (zh)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113507426A (zh) * 2021-06-16 2021-10-15 北京邮电大学 基于otfs调制的联合信道估计与信号检测方法及装置
CN113804981A (zh) * 2021-09-15 2021-12-17 电子科技大学 一种时频联合最优化多源多信道信号分离方法
CN113852580A (zh) * 2021-09-14 2021-12-28 电子科技大学 一种基于多级分离的mimo-otfs符号检测方法
CN113890796A (zh) * 2021-09-30 2022-01-04 成都工业学院 基于otfs系统调制解调的高速信道估计装置及方法
CN113890797A (zh) * 2021-09-30 2022-01-04 哈尔滨工业大学 一种基于短包通信传输过程的信道估计方法
CN114244670A (zh) * 2021-12-08 2022-03-25 北京理工大学 一种基于信道编码辅助的盲信道估计方法及系统
CN114337876A (zh) * 2021-12-31 2022-04-12 南京邮电大学 一种基于nsga2算法训练的amp检测算法和实施该算法的系统
CN114374447A (zh) * 2022-01-06 2022-04-19 中国人民解放军国防科技大学 一种信道探测的方法、装置及介质
CN114665932A (zh) * 2022-03-17 2022-06-24 东南大学 大规模mimo波束时延多普勒域统计信道信息获取方法
CN114690212A (zh) * 2022-03-18 2022-07-01 中国科学院国家空间科学中心 一种利用ddm数据质量监控进行干扰预警的方法及系统
CN114726688A (zh) * 2022-04-11 2022-07-08 桂林电子科技大学 一种信道估计方法、系统、设备及可读存储介质
CN114884777A (zh) * 2022-04-28 2022-08-09 中国科学院计算技术研究所 一种基于变换域的信道估计方法
CN115051894A (zh) * 2022-06-10 2022-09-13 西安电子科技大学 基于串行干扰消除的otfs系统符号检测方法
CN115086114A (zh) * 2022-06-10 2022-09-20 西安电子科技大学 基于分散式放置正交时频空otfs导频的信道估计方法
CN115174335A (zh) * 2022-07-19 2022-10-11 西安电子科技大学 基于交织置换机制的otfs安全传输方法
CN115296970A (zh) * 2022-07-07 2022-11-04 北京航空航天大学 基于逐元素外部信息的迭代正交时频空波形检测方法
CN115426223A (zh) * 2022-08-10 2022-12-02 华中科技大学 一种低轨卫星信道估计和符号检测方法及系统
CN115967595A (zh) * 2022-11-30 2023-04-14 西安电子科技大学 一种mimo-otfs系统的半盲信道估计方法
CN116055261A (zh) * 2023-01-17 2023-05-02 重庆邮电大学 一种基于模型驱动深度学习的otfs信道估计方法
CN116094876A (zh) * 2023-03-10 2023-05-09 南京邮电大学 一种基于非对称架构的正交时频空系统的信道估计方法
CN116346164A (zh) * 2023-03-13 2023-06-27 南京邮电大学 Otfs系统中连续并行干扰消除的最大似然检测方法
CN116633461A (zh) * 2022-05-31 2023-08-22 南方科技大学 一种信号检测的方法、装置和设备
CN114690212B (zh) * 2022-03-18 2024-05-14 中国科学院国家空间科学中心 一种利用ddm数据质量监控进行干扰预警的方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108156101A (zh) * 2017-12-18 2018-06-12 中国人民解放军空军工程大学 一种mimo-scfde系统联合迭代信道估计与迭代均衡方法
WO2020034394A1 (zh) * 2018-08-13 2020-02-20 南京邮电大学 基于压缩感知的大规模mimo信道反馈重构算法
CN111555780A (zh) * 2020-01-09 2020-08-18 北京邮电大学 基于正交时频空调制的多天线接收机设计
EP3761583A1 (en) * 2019-07-04 2021-01-06 Volkswagen AG Adaptive transmitter symbol arrangement for otfs channel estimation in the delay-doppler domain
CN112202479A (zh) * 2020-08-19 2021-01-08 北京邮电大学 多入多出-正交时频空系统的低复杂度信号检测方法
US20210036823A1 (en) * 2016-03-31 2021-02-04 Cohere Technologies, Inc. Channel acquisition using orthogonal time frequency space modulated pilot signals
US20210105155A1 (en) * 2018-02-08 2021-04-08 Cohere Technologies, Inc. Aspects of channel estimation for orthogonal time frequency space modulation for wireless communications

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210036823A1 (en) * 2016-03-31 2021-02-04 Cohere Technologies, Inc. Channel acquisition using orthogonal time frequency space modulated pilot signals
CN108156101A (zh) * 2017-12-18 2018-06-12 中国人民解放军空军工程大学 一种mimo-scfde系统联合迭代信道估计与迭代均衡方法
US20210105155A1 (en) * 2018-02-08 2021-04-08 Cohere Technologies, Inc. Aspects of channel estimation for orthogonal time frequency space modulation for wireless communications
WO2020034394A1 (zh) * 2018-08-13 2020-02-20 南京邮电大学 基于压缩感知的大规模mimo信道反馈重构算法
EP3761583A1 (en) * 2019-07-04 2021-01-06 Volkswagen AG Adaptive transmitter symbol arrangement for otfs channel estimation in the delay-doppler domain
CN111555780A (zh) * 2020-01-09 2020-08-18 北京邮电大学 基于正交时频空调制的多天线接收机设计
CN112202479A (zh) * 2020-08-19 2021-01-08 北京邮电大学 多入多出-正交时频空系统的低复杂度信号检测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
LEI ZHAO等: "Pilot Optimization for OFDM-based OTFS Systems over Doubly Selective Channels", 《GLOBECOM 2020 - 2020 IEEE GLOBAL COMMUNICATIONS CONFERENCE》 *
M. KOLLENGODE等: "MIMO-OTFS in High-Doppler Fading Channels: Signal Detection and Channel Estimation", 《2018 IEEE GLOBAL COMMUNICATIONS CONFERENCE (GLOBECOM)》 *
郭文彬等: "基于联合稀疏模型的OFDM压缩感知信道估计", 《北京邮电大学学报》 *
陈发堂等: "一种联合信道估计与多用户检测的新型算法", 《南京邮电大学学报(自然科学版)》 *

Cited By (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113507426A (zh) * 2021-06-16 2021-10-15 北京邮电大学 基于otfs调制的联合信道估计与信号检测方法及装置
CN113507426B (zh) * 2021-06-16 2022-12-06 北京邮电大学 基于otfs调制的联合信道估计与信号检测方法及装置
CN113852580A (zh) * 2021-09-14 2021-12-28 电子科技大学 一种基于多级分离的mimo-otfs符号检测方法
CN113804981A (zh) * 2021-09-15 2021-12-17 电子科技大学 一种时频联合最优化多源多信道信号分离方法
CN113890796A (zh) * 2021-09-30 2022-01-04 成都工业学院 基于otfs系统调制解调的高速信道估计装置及方法
CN113890797A (zh) * 2021-09-30 2022-01-04 哈尔滨工业大学 一种基于短包通信传输过程的信道估计方法
CN113890797B (zh) * 2021-09-30 2024-04-19 哈尔滨工业大学 一种基于短包通信传输过程的信道估计方法
CN114244670A (zh) * 2021-12-08 2022-03-25 北京理工大学 一种基于信道编码辅助的盲信道估计方法及系统
CN114244670B (zh) * 2021-12-08 2023-04-18 北京理工大学 一种基于信道编码辅助的盲信道估计方法及系统
CN114337876A (zh) * 2021-12-31 2022-04-12 南京邮电大学 一种基于nsga2算法训练的amp检测算法和实施该算法的系统
CN114337876B (zh) * 2021-12-31 2024-01-23 南京邮电大学 一种基于nsga2算法训练的amp检测方法和实施该方法的系统
CN114374447A (zh) * 2022-01-06 2022-04-19 中国人民解放军国防科技大学 一种信道探测的方法、装置及介质
CN114374447B (zh) * 2022-01-06 2024-02-02 中国人民解放军国防科技大学 一种信道探测的方法、装置及介质
CN114665932A (zh) * 2022-03-17 2022-06-24 东南大学 大规模mimo波束时延多普勒域统计信道信息获取方法
CN114665932B (zh) * 2022-03-17 2023-03-10 东南大学 大规模mimo波束时延多普勒域统计信道信息获取方法
CN114690212A (zh) * 2022-03-18 2022-07-01 中国科学院国家空间科学中心 一种利用ddm数据质量监控进行干扰预警的方法及系统
CN114690212B (zh) * 2022-03-18 2024-05-14 中国科学院国家空间科学中心 一种利用ddm数据质量监控进行干扰预警的方法及系统
CN114726688A (zh) * 2022-04-11 2022-07-08 桂林电子科技大学 一种信道估计方法、系统、设备及可读存储介质
CN114884777A (zh) * 2022-04-28 2022-08-09 中国科学院计算技术研究所 一种基于变换域的信道估计方法
CN116633461A (zh) * 2022-05-31 2023-08-22 南方科技大学 一种信号检测的方法、装置和设备
CN116633461B (zh) * 2022-05-31 2024-02-09 南方科技大学 一种信号检测的方法、装置和设备
CN115086114A (zh) * 2022-06-10 2022-09-20 西安电子科技大学 基于分散式放置正交时频空otfs导频的信道估计方法
CN115086114B (zh) * 2022-06-10 2023-08-15 西安电子科技大学 基于分散式放置正交时频空otfs导频的信道估计方法
CN115051894B (zh) * 2022-06-10 2024-02-06 西安电子科技大学 基于串行干扰消除的otfs系统符号检测方法
CN115051894A (zh) * 2022-06-10 2022-09-13 西安电子科技大学 基于串行干扰消除的otfs系统符号检测方法
CN115296970B (zh) * 2022-07-07 2024-01-02 北京航空航天大学 基于逐元素外部信息的迭代正交时频空波形检测方法
CN115296970A (zh) * 2022-07-07 2022-11-04 北京航空航天大学 基于逐元素外部信息的迭代正交时频空波形检测方法
CN115174335A (zh) * 2022-07-19 2022-10-11 西安电子科技大学 基于交织置换机制的otfs安全传输方法
CN115426223A (zh) * 2022-08-10 2022-12-02 华中科技大学 一种低轨卫星信道估计和符号检测方法及系统
CN115426223B (zh) * 2022-08-10 2024-04-23 华中科技大学 一种低轨卫星信道估计和符号检测方法及系统
CN115967595A (zh) * 2022-11-30 2023-04-14 西安电子科技大学 一种mimo-otfs系统的半盲信道估计方法
CN116055261A (zh) * 2023-01-17 2023-05-02 重庆邮电大学 一种基于模型驱动深度学习的otfs信道估计方法
CN116094876B (zh) * 2023-03-10 2023-08-29 南京邮电大学 一种基于非对称架构的正交时频空系统的信道估计方法
CN116094876A (zh) * 2023-03-10 2023-05-09 南京邮电大学 一种基于非对称架构的正交时频空系统的信道估计方法
CN116346164B (zh) * 2023-03-13 2023-11-24 南京邮电大学 Otfs系统中连续并行干扰消除的最大似然检测方法
CN116346164A (zh) * 2023-03-13 2023-06-27 南京邮电大学 Otfs系统中连续并行干扰消除的最大似然检测方法

Also Published As

Publication number Publication date
CN113395221B (zh) 2022-07-08

Similar Documents

Publication Publication Date Title
CN113395221B (zh) 一种基于正交时频空联合信道估计与符号检测方法
CN111404849B (zh) 一种基于深度学习的ofdm信道估计与信号检测方法
CN108833311B (zh) 联合时域聚类去噪与均衡判决的变换域二次估计方法
CN112600772B (zh) 一种基于数据驱动神经网络的ofdm信道估计与信号检测方法
CN111683024B (zh) 一种基于深度学习的时变ofdm系统信道估计方法
CN107682297A (zh) 一种移动水声通信方法
CN107332797B (zh) 一种电力线ofdm通信系统中的信道估计方法
CN113472706A (zh) 一种基于深度神经网络的mimo-ofdm系统信道估计方法
CN112968847B (zh) 一种基于深度学习和数据导频辅助的信道估计方法
CN113381951A (zh) 时变频选衰落信道下mftn联合信道估计与均衡方法
CN106453162A (zh) 针对多输入多输出正交频分复用系统的信道估计方法
CN114338305A (zh) 一种针对正交时频空调制系统的符号检测方法
CN101322365B (zh) 多载波系统中的噪声功率插值
Wang et al. Underwater acoustic sparse channel estimation based on DW-SACoSaMP reconstruction algorithm
CN101252560A (zh) 一种高性能的ofdm帧同步算法
CN102790746B (zh) Ofdm系统的信道估计方法
CN105812299A (zh) 基于联合块稀疏重构的无线传感网信道估计算法及系统
CN106911621B (zh) 一种基于v-ofdm的信道均衡和跟踪方法
CN111291511B (zh) 一种基于历史信息的软卡尔曼滤波迭代时变信道估计方法
CN112636855A (zh) 一种ofdm信号检测方法
CN115426224B (zh) 一种基于otfs信号的信道估计方法及系统
CN101662441B (zh) 信号估计方法、装置及检测系统
CN110113084A (zh) Mimo闭环传输系统的信道预测方法
CN113709075B (zh) 一种利用水声信道多径效应的水声通信接收机的实现方法
CN102231720B (zh) 样条函数Renyi熵与时间分集融合的小波盲均衡方法

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