CN116094876B - 一种基于非对称架构的正交时频空系统的信道估计方法 - Google Patents
一种基于非对称架构的正交时频空系统的信道估计方法 Download PDFInfo
- Publication number
- CN116094876B CN116094876B CN202310228849.3A CN202310228849A CN116094876B CN 116094876 B CN116094876 B CN 116094876B CN 202310228849 A CN202310228849 A CN 202310228849A CN 116094876 B CN116094876 B CN 116094876B
- Authority
- CN
- China
- Prior art keywords
- frequency
- time
- delay
- channel
- angle
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 10
- 238000004891 communication Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000005516 engineering process Methods 0.000 claims description 8
- 238000004378 air conditioning Methods 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000011478 gradient descent method Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 2
- 101150057833 THEG gene Proteins 0.000 claims 1
- 230000003247 decreasing effect Effects 0.000 claims 1
- 230000001934 delay Effects 0.000 claims 1
- 230000033001 locomotion Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0802—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using antenna selection
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Radio Transmission System (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于非对称架构的正交时频空系统的信道估计方法包括,基于互质阵列设计接收天线阵列元素,弥补非对称架构中阵列孔径损失,以实现更高精度的角度参数估计;基于互质阵列设计针对性的角度估计算法;在完成角度估计的基础上,实现对不同路径的时延、多普勒频移和信道增益估计,以恢复上行信道,重建下行信道;本发明可实现对基于非对称架构的正交时频空系统信道状态信息估计的同时,提升了参数估计的精度,大大降低了算法复杂度,提高了非对称架构系统的性能。
Description
技术领域
本发明涉及无线通信技术领域,尤其涉及一种基于非对称架构的正交时频空系统的信道估计方法。
背景技术
在例如自动驾驶汽车、无人机、低地球轨道卫星、高速列车等场景中,由于收发器之间的相对运动,无线通信会受到严重的多普勒频移效应的影响。在高移动性场景中,无线信道通常在时间频率域中是双色散的。具体来说,多径效应会导致时间色散,而频率色散是由多普勒频移引起的。然而,作为5G中的关键技术之一的正交频分复用调制技术,容易受频移的影响,导致正交频分复用技术子载波之间的完美正交性被破坏,可能无法支持在高载波频率下的高移动性场景中的稳定高速通信。因此,有研究人员提出了一种名为正交时频空的二维调制技术,以作为支持6G无线系统要求的潜在解决方案之一。
相较于经典的正交频分复用技术,正交时频空具有优异的抗多普勒频移鲁棒性。正交时频空将每个比特信息符号调制到一组跨越信号传输带宽与时间周期的二维正交基函数上,称正交时频空调制信号所处的域为时延多普勒域。与在时间频率域中相比,双色散信道在时延多普勒域中相对更加稳定,表现出稀疏性,仅仅需少量参数即可表征信道状态信息。
与此同时,考虑到全数字大规模MIMO技术在提高系统频谱效率与系统信道容量方面的巨大优势,通过将正交时频空技术与全数字大规模MIMO技术结合起来,构建正交时频空系统为用户提供在高载波频率下的高移动性场景中的高质量高速率通信。然而,关于全数字正交时频空系统存在着算法复杂度过高、系统硬件成本开销过大等问题,极大的阻碍了将大规模MIMO应用到实际通信的脚步。通过设计全数字非对称收发架构,即仅允许上行链路使用部分天线进行接收,在满足用户通信需求的情况下,极大的降低系统成本与复杂度,为此,也需要解决上下行链路天线阵列孔径不一致所带来的路径参数估计精度差异的问题。但是目前的研究大多基于对称架构的正交时频空系统所设计,并没有考虑到非对称架构的正交时频空系统的特殊性,现有的信道估计方法将不再适用。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于上述现有存在的问题,提出了本发明。因此,本发明提供了一种基于非对称架构的正交时频空系统的信道估计方法,用来解决实际问题中,非对称架构中的阵列孔径不一致所带来的路径参数估计精度低、系统成本高和算法复杂的问题。
为解决上述技术问题,本发明提供如下技术方案:
本发明提供了一种基于非对称架构的正交时频空系统的信道估计方法,包括:
基于大规模均匀线性阵列,依据互质阵列的排布选择接收射频链的位置,设计上行链路的线性接收天线阵列 ,构建非对称系统的上行信道模型 /> ;
基于正交时频空调制方案,获得非对称架构中上行链路的输出信号表达式 ,通过所述输出信号表达式推导获得上行链路正交时频空的输入输出关系,并根据此输入输出关系,从中提取角度信息,重构虚拟的均匀线性接收天线阵列 /> ,得到到达角的估计值 />;
于获得的到达角估计值参数,利用信号在时延多普勒维度的稀疏性,设计时延多普勒二维搜索网格 ,筛选与角度相对应的最匹配的路径时延 />和多普勒频移/> ,并采用最小二乘法计算复信道增益/> ;
将所得到的上行信道参数:到达角、时延、多普勒频移与复信道增益,基于信道互易性,重构恢复出下行信道 。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:基于非对称架构的正交时频空系统,包括:
所述系统由一个配备了 根天线的基站与/>个高速移动的单天线用户组成,系统采用正交时频空调制方案,工作在时分双工模式;
发送端用户首先是将放置在时延多普勒域服从复高斯分布的导频序列,通过逆辛有限傅里叶变换将其变换到时间频率域,然后每个二维正交时频空符号经过正交频分复用技术调制后获得一维的时域信号,最后用户在不同时隙将信号发送给基站;
首先所述基站依据根天线的接收信号,恢复出 /> 根天线的虚拟阵列的接收信号,再从中提取角度信息,然后执行与发送端相对应的逆操作来解调信号,最后估计出到达角、时延、多普勒频移和复移信道增益参数,并重构出下行信道以辅助下行通信。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:基于大规模均匀线性阵列,依据互质阵列的排布选择接收射频链的位置,设计上行链路的线性接收天线阵列 ,包括:
基站配备了一个均匀线性阵列,共根天线;
设计上行链路,使得只有 个接收射频链链接到线性阵列,得到上行导向矢量表示为:
;
其中, 表示第/> 个用户的第/> 条路径的到达角;接收天线的索引是通过互质阵列确定的,表示为/>且 /> ,
, /> 与 /> 都为质数,/>为虚数。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:构建非对称系统的上行信道模型 ,包括:
上行信道模型公式表示为:
;
其中, 分别表示复信道增益、时延与多普勒频移, />为路径数,为时延抽头索引,/> 为采样频率,/>表示狄拉克函数,n为时隙。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:基于正交时频空调制方案,获得非对称架构中上行链路的输出信号表达式 ,包括:
输出信号是由发送信号经历时频双选信道后获得的,因此上行链路的输出信号在时域的公式表示为:
;
其中, 表示为一个正交时频空帧沿着时延与多普勒索引方向的长度; 与/> 分别表示路径的物理频移与时延在时延多普勒域的映射;
对于实际应用场景中的用户其起始时间已知,/> 表示有效点的个数,导频序列/> 服从复高斯分布, /> 为噪声。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:通过所述输出信号表达式推导获得上行链路的正交时频空的输入输出关系,根据此输入输出关系,从中提取角度信息,重构虚拟的均匀线性接收天线阵列,包括:
计算接收信号 的协方差矩阵公式表示为:
;
其中,表示导频信号的方差,/>表示噪音的方差,/>为单位向量;
从所述的协方差矩阵中提取角度信息,进行选择、重排得到虚拟的上行均匀线性阵列 的观测信号,公式表示为:
;
其中,观测信号 是维度为 /> 的列向量,因为路径数P<<A,此时的角度估计问题等效于一个稀疏信号恢复问题;
设计字典矩阵为:
;
其中, 表示可能的角度;
通过从字典矩阵 中选择约束问题/> 中最小一列作为估计的角度,/>为所有路径的噪声功率,/>为惩罚因子。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:得到到达角的估计值 ,包括:
得到到达角的估计值,采用近端梯度下降法;
令 , /> ,计算函数/> 的下降梯度,公式表示为:
;
其中, 表示第/> 次迭代, />表示步长;然后将 /> 代入梯度算子中,迭代更新软阈值函数得到/> ,直到循环收敛,于是得到角度的估计结果/> 。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:基于获得的到达角估计值参数,利用信号在时延多普勒维度的稀疏性,包括:
在获得角度估计值 之后,将接收信号沿数据方向堆叠,获得维度为的列向量表示为:
;
其中,多普勒分量与导频的时延分量分别为:
,和
。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:设计时延多普勒二维搜索网格 ,筛选与角度相对应的最匹配的路径时延 与多普勒频移/> ,并采用最小二乘法计算复信道增益/> ,包括:
创建不考虑分数时延的二维时延多普勒网格 ,令时延方向网格精度为1,多普勒方向网格精度为 /> ,表示为:
;
初始化残差 ,从/>中依次选择网格点进行匹配筛选,计算公式为:
;
其中, ,选择最大值对应的网格点作为时延与多普勒的估计值/> ,并使用最小二乘法计算/> 更新复信道增益,同时更新每一步的残差结果 /> ,重复上述过程/> 次后,获得信道剩余参数的所有估计值 /> 。
作为本发明所述的基于非对称架构的正交时频空系统的信道估计方法的一种优选方案,其中:基于信道互易性,重构恢复出下行信道 ,包括:
设计下行链路,使得每根天线都链接到一个发送射频链,得到下行导向矢量表示为:
;
其中,表示第/> 个用户的第 />条路径的离开角;
下行信道模型公式表示为:
;
基于信道的互易性,上行信道的信道参数可以直接用于重建下行信道,以辅助下行通信。
与现有技术相比,发明有益效果为:本发明所述方法能够通过合理选择接收天线阵列分布,消除上下行天线阵列孔径不一致所带来的估计精度下降的影响,实现准确的角度估计,基于所估计的角度,从接收信号中提取有关路径时延、多普勒频移参数的精确结果;并将角度、时延与多普勒频移的三维搜索空间转换为一维角度搜索与二维时间频率搜索网格,显著降低了算法的复杂度,同时提升了在低信噪比时的参数估计性能,更有效准确的重建下行信道;同时,基于非对称架构的正交时频空系统,适用于用户快速移动和静止的复杂场景,具有较低的算法时延与系统硬件开销,更加适用于实际系统。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为本发明一个实施例所述的基于非对称架构的正交时频空系统的信道估计方法的整体流程图;
图2为本发明一个实施例所述的基于非对称架构的正交时频空系统的信道估计方法的正交时频空系统调制解调示意图;
图3为本发明一个实施例所述的基于非对称架构的正交时频空系统的信道估计方法的在不同信噪比下下信道重建的均方根误差的仿真示意图;
图4为本发明一个实施例所述的基于非对称架构的正交时频空系统的信道估计方法的在不同网格精度下算法运行时间的仿真示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
本发明结合示意图进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的剖面图会不依一般比例作局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。此外,在实际制作中应包含长度、宽度及深度的三维空间尺寸。
同时在本发明的描述中,需要说明的是,术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一、第二或第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
本发明中除非另有明确的规定和限定,术语“安装、相连、连接”应做广义理解,例如:可以是固定连接、可拆卸连接或一体式连接;同样可以是机械连接、电连接或直接连接,也可以通过中间媒介间接相连,也可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
参照图1,为本发明第一个实施例,该实施例提供了一种基于非对称架构的正交时频空系统的信道估计方法,包括:
S1、基于大规模均匀线性阵列,依据互质阵列的排布选择接收射频链的位置,设计上行链路的线性接收天线阵列 ,构建非对称系统的上行信道模型/> ;
基站配备了一个均匀线性阵列,共根天线;
设计上行链路,使得只有 个接收射频链链接到线性阵列,得到上行导向矢量表示为:
;
其中, 表示第/> 个用户的第/> 条路径的到达角;接收天线的索引是通过互质阵列确定的,表示为/> ,
,/> ,/>与 /> 都为质数,/>为虚数;
进一步的,上行信道模型公式表示为:
;
其中, 分别表示复信道增益、时延与多普勒频移, />为路径数,为时延抽头索引, />为采样频率,/>表示狄拉克函数,n为时隙;
S2、基于正交时频空调制方案,获得非对称架构中上行链路的输出信号表达式 ,通过输出信号表达式推导获得上行链路正交时频空的输入输出关系,并根据此输入输出关系,从中提取角度信息,重构虚拟的均匀线性接收天线阵列/> ,得到到达角的估计值 />;
输出信号是由发送信号经历时频双选信道后获得的,因此上行链路的输出信号在时域的公式表示为:
;
其中, 表示为一个正交时频空帧沿着时延与多普勒索引方向的长度; 与/> 分别表示路径的物理频移与时延在时延多普勒域的映射;
对于实际应用场景中的用户其起始时间已知, /> 表示有效点的个数,
导频序列 服从复高斯分布,/>为噪声;
进一步的,计算接收信号 的协方差矩阵公式表示为:
;
其中,表示导频信号的方差,/>表示噪音的方差,/>为单位向量;
从协方差矩阵中提取角度信息,进行选择、重排得到虚拟的上行均匀线性阵列 的观测信号,公式表示为:
;
其中,观测信号 是维度为/> 的列向量,因为路径数P<<A,此时的角度估计问题等效于一个稀疏信号恢复问题;
进一步的,设计字典矩阵为:
;
其中, 表示可能的角度;
通过从字典矩阵 中选择约束问题/> 中最小一列作为估计的角度,/>为所有路径的噪声功率,/>为惩罚因子;
进一步的,得到到达角的估计值,采用近端梯度下降法;
令 , /> ,计算函数/> 的下降梯度,公式表示为:
;
其中,表示第/> 次迭代, />表示步长;然后将 /> 代入梯度算子中,迭代更新软阈值函数得到/> ,直到循环收敛,于是得到角度的估计结果 />;
S3、在所估计获得的角度参数的基础上,并利用信号在时延多普勒维度的稀疏性,设计时延多普勒二维搜索网格 ,从中筛选与角度相对应的最匹配的路径时延/> 与多普勒频移/> ,最后利用最小二乘法计算复信道增益/> ;
进一步的,在获得角度估计值 之后,将接收信号沿数据方向堆叠,获得维度为 /> 的列向量表示为:
;
其中,多普勒分量与导频的时延分量分别为:
,和
;
进一步的,创建不考虑分数时延的二维时延多普勒网格 ,令时延方向网格精度为1,多普勒方向网格精度为/> ,表示为:
;
初始化残差 ,从/> 中依次选择网格点进行匹配筛选,计算公式为:
;
其中, ,选择最大值对应的网格点作为时延与多普勒的估计值 ,并使用最小二乘法计算 /> 更新复信道增益,同时更新每一步的残差结果 /> ,重复上述过程/> 次后,获得信道剩余参数的所有估计值 /> ;
S4、将所得到的上行信道参数:到达角、时延、多普勒频移与复信道增益,基于信道互易性,重构恢复出下行信道;
进一步的,设计下行链路,使得每根天线都链接到一个发送射频链,得到下行导向矢量表示为:
;
其中,表示第 />个用户的第 />条路径的离开角;
下行信道模型公式表示为:
;
基于信道的互易性,上行信道的信道参数 可以直接用于重建下行信道,以辅助下行通信。
实施例2
参照图2,为本发明第二个实施例,该实施例提供了一种基于非对称架构的正交时频空系统的信道估计方法,包括:
正交时频空系统上行链路的调制解调过程如图2所示;
发送端用户首先是将放置在时延多普勒域服从复高斯分布的导频序列,通过逆辛有限傅里叶变换将其变换到时间频率域,然后每个二维正交时频空符号经过正交频分复用技术调制后获得一维的时域信号,最后用户在不同时隙将信号发送给基站;
所述基站依据现实场景交互的需要,设置23根天线的接收信号,通过正交频分复用技术调制后恢复出144根天线的虚拟阵列的接收信号,再从中提取角度信息,然后执行与发送端相对应的逆操作来解调信号,最后估计出到达角、时延、多普勒频移与复信道增益,并重构出下行信道以辅助下行通信。
实施例3
参照图3,为本发明第三个实施例,该实施例提供了一种基于非对称架构的正交时频空系统的信道估计方法,包括:
OTFS的帧大小被设置为,上行链路的接收天线数量设置为/> ,下行链路的发送天线数量设置为/> ,系统载频设置为 /> GHz,子载波间隔设置为 /> kHz;设定用户的最大运动速度为540 km/h,则相应的多普勒频谱范围为[-2kHz,2 kHz],因此最大多普勒索引/> ,最大时延索引 /> ,角度设为在范围/> 内均匀分布;
由图3可见,随着信噪比的增加,各类方法的信道估计误差都在逐步降低,从整个信噪比区间可以看出,本发明方法的误差均小于其他两种方法,特别是在低信噪比下,本发明的方法与其他两种算法的性能差距更大,这表明采用互质阵列能够有效消除上下行阵列孔径差异,改善信道参数估计性能,体现了本发明方法的有效性。
实施例4
参照图4,为本发明第四个实施例,该实施例提供了一种基于非对称架构的正交时频空系统的信道估计方法,包括:
与实施例3进行同样的设置,将OTFS的帧大小被设置为,上行链路的接收天线数量设置为/> ,下行链路的发送天线数量设置为/> ,系统载频设置为 GHz,子载波间隔设置为/> kHz;设定用户的最大运动速度为540 km/h,则相应的多普勒频谱范围为[-2 kHz,2 kHz],得到最大多普勒索引为 /> ,最大时延索引/> ,角度设为在范围/> 内均匀分布;
由图4可见,随着网格精度的降低,各方法的时间开销都在减小,但在相同的运行时间内,本发明的方法可以采用更小的网格精度,即在相同的时间、硬件开销条件下,本发明的方法能够实现更精确的估计,更好的性能,这体现了本发明方法的低复杂性,能够有效降低算法时延,提升效率。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (2)
1.一种基于非对称架构的正交时频空系统的信道估计方法,其特征在于,包括:
基于大规模均匀线性阵列,依据互质阵列的排布选择接收射频链的位置,设计上行链路的线性接收天线阵列aUL,包括:
基站配备了一个均匀线性阵列,共At根天线;
设计上行链路,使得只有Ar个接收射频链链接到线性阵列,得到上行导向矢量表示为:
其中,θq,p表示第q个用户的第p条路径的到达角或离开角;θq,p在上行导向矢量中表示为到达角,在下行导向矢量中表示离开角;表示接收天线的索引,接收天线的索引是通过互质阵列确定的,表示为{L1i|0≤i≤(L-1)/2}∪{L2j|0≤j≤(L-1)/2}且At=L1L2+1,Ar=L1+L2-1,L1与L2都为质数,/>为虚数;
构建非对称系统的上行信道模型hUL,包括:
上行信道模型公式表示为:
其中,hq,p,τq,p,νq,p分别表示复信道增益、时延与多普勒频移,P为路径数,l为时延抽头索引,Ts为采样频率,prc表示狄拉克函数,n为时隙;
基于正交时频空调制方案,获得非对称架构中上行链路的输出信号表达式yg,包括:
输出信号是由发送信号经历时频双选信道后获得的,因此上行链路的输出信号在时域的公式表示为:
其中,M,N表示为一个正交时频空帧沿着时延与多普勒索引方向的长度;kp=νpNT与lp=τpMΔf分别表示路径的物理频移与时延在时延多普勒域的映射;hp,θp,νp和τp分别表示第p条路径的信道增益、离开角、多普勒频移和传播时延,Δf表示子载波间隔,T表示一个符号的持续时间;对于实际应用场景中的用户其起始时间已知,/>表示有效点的个数,c=[c0,c1,...,cG-1]T表示服从复高斯分布的导频序列,wg为噪声;w表示噪声向量;表示有效点的总个数,[]T表示向量的转置;
通过所述输出信号表达式推导获得上行链路正交时频空的输入输出关系,并根据此输入输出关系,从中提取角度信息,重构虚拟的均匀线性接收天线阵列包括:
计算接收信号yg的协方差矩阵公式表示为:
其中,σc表示导频信号的方差,σw表示噪音的方差;IG为一个全零列向量;
从所述的协方差矩阵中提取角度信息,进行选择、重排得到虚拟的上行均匀线性阵列的观测信号,公式表示为:
其中,观测信号是维度为At×1的列向量,因为路径数P<<A,此时的角度估计问题等效于一个稀疏信号恢复问题;
设计字典矩阵为:D=[D1,...,Dη,D0]
其中,表示可能的角度;
通过从字典矩阵D中选择约束问题中最小一列作为估计的角度,r为所有路径的噪声功率,λc为惩罚因子;
得到到达角的估计值包括:
得到到达角的估计值,采用近端梯度下降法;
令g(r)=λc‖r‖1,计算函数f(r)的下降梯度,公式表示为:
其中,(i)表示第i次迭代,μc表示算法每次迭代的步长;然后将z(i)代入梯度算子中,迭代更新软阈值函数得到r(i+1)=argminr{g(r)+1/2μc‖z(i)-r‖2},直到循环收敛,于是得到角度的估计结果
基于获得的到达角估计值参数,利用信号在时延多普勒维度的稀疏性,包括:
在获得角度估计值之后,将接收信号沿数据方向堆叠,获得维度为/>的列向量表示为:
其中,多普勒分量与导频的时延分量分别为:和 表示导频序列等同于cG;
设计时延多普勒二维搜索网格Γ,筛选与角度相对应的匹配的路径时延l和多普勒频移并采用最小二乘法计算复信道增益/>包括:
创建不考虑分数时延的二维时延多普勒网格Γ,令时延方向网格精度为1,多普勒方向网格精度为Δδ,表示为:
Γ={(d1,d2),d1=0,1,...,lmax;d2=-kmax,-kmax+Δδ,...,kmax}
初始化残差yr=y,从Γ中依次选择网格点进行匹配筛选,计算公式为:
其中,选择最大值对应的网格点作为时延与多普勒的估计值/>并使用最小二乘法计算/>更新复信道增益,同时更新每一步的残差结果/>重复上述过程P次后,获得信道剩余参数的所有估计值/>
将所得到的上行信道参数:到达角、时延、多普勒频移与复信道增益,基于信道互易性,重构恢复出下行信道hDL,包括:
设计下行链路,使得每根天线都链接到一个发送射频链,得到下行导向矢量表示为:
下行信道模型公式表示为:
基于信道的互易性,上行信道的信道参数可以直接用于重建下行信道,以辅助下行通信。
2.如权利要求1所述的基于非对称架构的正交时频空系统的信道估计方法,其特征在于,基于非对称架构的正交时频空系统,包括:
所述系统由一个配备了At根天线的基站与Q个高速移动的单天线用户组成,系统采用正交时频空调制方案,工作在时分双工模式;
发送端用户首先是将放置在时延多普勒域服从复高斯分布的导频序列,通过逆辛有限傅里叶变换将其变换到时间频率域,然后每个二维正交时频空符号经过正交频分复用技术调制后获得一维的时域信号,最后用户在不同时隙将信号发送给基站;
首先所述基站依据Ar根天线的接收信号,恢复出At根天线的虚拟阵列的接收信号,再从中提取角度信息,然后执行与发送端相对应的逆操作来解调信号,最后估计出到达角、时延、多普勒频移和复移信道增益参数,并重构出下行信道以辅助下行通信;
其中,Q表示几个高速移动的单天线用户,At表示几根天线的基站,Ar表示根天线的接收信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310228849.3A CN116094876B (zh) | 2023-03-10 | 2023-03-10 | 一种基于非对称架构的正交时频空系统的信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310228849.3A CN116094876B (zh) | 2023-03-10 | 2023-03-10 | 一种基于非对称架构的正交时频空系统的信道估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116094876A CN116094876A (zh) | 2023-05-09 |
CN116094876B true CN116094876B (zh) | 2023-08-29 |
Family
ID=86210326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310228849.3A Active CN116094876B (zh) | 2023-03-10 | 2023-03-10 | 一种基于非对称架构的正交时频空系统的信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116094876B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016161009A1 (en) * | 2015-04-01 | 2016-10-06 | Verasonics, Inc. | Method and system for coded excitation imaging by impulse response estimation and retrospective acquisition |
WO2018032016A1 (en) * | 2016-08-12 | 2018-02-15 | Cohere Technologies | Localized equalization for channels with intercarrier interference |
WO2018049595A1 (zh) * | 2016-09-14 | 2018-03-22 | 深圳大学 | 一种基于交替方向乘子法的稳健稀疏恢复stap方法及其系统 |
CN112087247A (zh) * | 2020-08-04 | 2020-12-15 | 西安电子科技大学 | 一种基于大规模mimo-otfs的径分多址接入方法 |
CN112887233A (zh) * | 2021-01-21 | 2021-06-01 | 中国科学技术大学 | 一种基于2维簇结构的稀疏贝叶斯学习信道估计方法 |
CN112882005A (zh) * | 2021-01-14 | 2021-06-01 | 东南大学 | 一种基于贝叶斯学习的otfs雷达目标参数估计方法 |
CN113395221A (zh) * | 2021-04-25 | 2021-09-14 | 北京邮电大学 | 一种基于正交时频空联合信道估计与符号检测方法 |
CN113507426A (zh) * | 2021-06-16 | 2021-10-15 | 北京邮电大学 | 基于otfs调制的联合信道估计与信号检测方法及装置 |
CN113765581A (zh) * | 2021-09-27 | 2021-12-07 | 北京理工大学 | 基于压缩感知与波束对齐的ris快时变信道估计方法 |
CN114401171A (zh) * | 2022-01-07 | 2022-04-26 | 南京邮电大学 | 一种基于互质阵列的非对称大规模mimo信道估计方法 |
-
2023
- 2023-03-10 CN CN202310228849.3A patent/CN116094876B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016161009A1 (en) * | 2015-04-01 | 2016-10-06 | Verasonics, Inc. | Method and system for coded excitation imaging by impulse response estimation and retrospective acquisition |
WO2018032016A1 (en) * | 2016-08-12 | 2018-02-15 | Cohere Technologies | Localized equalization for channels with intercarrier interference |
WO2018049595A1 (zh) * | 2016-09-14 | 2018-03-22 | 深圳大学 | 一种基于交替方向乘子法的稳健稀疏恢复stap方法及其系统 |
CN112087247A (zh) * | 2020-08-04 | 2020-12-15 | 西安电子科技大学 | 一种基于大规模mimo-otfs的径分多址接入方法 |
CN112882005A (zh) * | 2021-01-14 | 2021-06-01 | 东南大学 | 一种基于贝叶斯学习的otfs雷达目标参数估计方法 |
CN112887233A (zh) * | 2021-01-21 | 2021-06-01 | 中国科学技术大学 | 一种基于2维簇结构的稀疏贝叶斯学习信道估计方法 |
CN113395221A (zh) * | 2021-04-25 | 2021-09-14 | 北京邮电大学 | 一种基于正交时频空联合信道估计与符号检测方法 |
CN113507426A (zh) * | 2021-06-16 | 2021-10-15 | 北京邮电大学 | 基于otfs调制的联合信道估计与信号检测方法及装置 |
CN113765581A (zh) * | 2021-09-27 | 2021-12-07 | 北京理工大学 | 基于压缩感知与波束对齐的ris快时变信道估计方法 |
CN114401171A (zh) * | 2022-01-07 | 2022-04-26 | 南京邮电大学 | 一种基于互质阵列的非对称大规模mimo信道估计方法 |
Non-Patent Citations (1)
Title |
---|
"R1-2101522 Joint channel estimation for PUSCH".3GPP tsg_ran\wg1_rl1.2021,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN116094876A (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111698182B (zh) | 一种基于压缩感知的时频分块稀疏的信道估计方法 | |
Liu et al. | Uplink-aided high mobility downlink channel estimation over massive MIMO-OTFS system | |
CN112737985A (zh) | 基于深度学习的大规模mimo信道联合估计和反馈方法 | |
CN113395221A (zh) | 一种基于正交时频空联合信道估计与符号检测方法 | |
CN103780521B (zh) | 一种稀疏度自适应的ofdm系统信道估计方法 | |
US20160315791A1 (en) | Methods and devices for channel estimation and ofdm receiver | |
CN111865843B (zh) | 大规模mimo-ofdm系统混合消息传递信道估计方法 | |
CN109005133B (zh) | 双稀疏多径信道模型及基于此模型的信道估计方法 | |
CN106972875B (zh) | 一种mimo系统下多维联合估计动态稀疏信道的方法 | |
CN103428154B (zh) | 基于Vector OFDM的双选择性信道的变换域复用方法 | |
CN104869086A (zh) | 基于二维压缩感知的mimo-ofdm通信系统下行信道估计方法、装置 | |
CN104539562A (zh) | 基于多输入多输出正交频分复用的宽带短波信道估计方法 | |
CN111355677A (zh) | 一种基于滤波器组的多载波水下高速通信系统 | |
Uwaechia et al. | Compressed channel estimation for massive MIMO-OFDM systems over doubly selective channels | |
CN105847192B (zh) | 一种动态稀疏信道的联合估计方法 | |
CN112769462B (zh) | 一种基于联合参数学习的毫米波mimo宽带信道估计方法 | |
CN112887233B (zh) | 一种基于2维簇结构的稀疏贝叶斯学习信道估计方法 | |
CN116094876B (zh) | 一种基于非对称架构的正交时频空系统的信道估计方法 | |
CN103346984A (zh) | 基于bsl0的局部聚类稀疏信道估计方法 | |
CN117749241A (zh) | 一种基于改进gomp的低轨卫星系统的信道估计方法 | |
CN115412416B (zh) | 一种面向高速移动场景的低复杂度otfs信号检测方法 | |
CN114448535B (zh) | 基于fdd部分互易性的信道预测方法、网络侧设备及终端 | |
CN114268346B (zh) | 非高斯噪声下电力线载波通信压缩感知信道估计方法 | |
Zheng et al. | Joint position, orientation and channel estimation in hybrid mmWave MIMO systems | |
CN111970216A (zh) | 一种基于vamp算法的mimo-ofdm系统稀疏信道估计的设计方法 |
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 |