CN112087247A - 一种基于大规模mimo-otfs的径分多址接入方法 - Google Patents

一种基于大规模mimo-otfs的径分多址接入方法 Download PDF

Info

Publication number
CN112087247A
CN112087247A CN202010771007.9A CN202010771007A CN112087247A CN 112087247 A CN112087247 A CN 112087247A CN 202010771007 A CN202010771007 A CN 202010771007A CN 112087247 A CN112087247 A CN 112087247A
Authority
CN
China
Prior art keywords
delay
channel
user
otfs
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.)
Granted
Application number
CN202010771007.9A
Other languages
English (en)
Other versions
CN112087247B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202010771007.9A priority Critical patent/CN112087247B/zh
Publication of CN112087247A publication Critical patent/CN112087247A/zh
Application granted granted Critical
Publication of CN112087247B publication Critical patent/CN112087247B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/0413—MIMO systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/382—Monitoring; Testing of propagation channels for resource allocation, admission control or handover
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • 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/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • 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/0891—Space-time diversity
    • H04B7/0897—Space-time diversity using beamforming per multi-path, e.g. to cope with different directions of arrival [DOA] at different multi-paths
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Power Engineering (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明属于无线通信技术领域，公开了一种基于大规模MIMO‑OTFS的径分多址接入方法，使用OTFS调制，得到大规模MIMO‑OTFS角度‑延迟‑多普勒域接收信号模型和角度‑延迟‑多普勒域上行主、次3D信道模型；利用大规模MIMO‑OTFS信道在角度‑延迟‑多普勒域的弥散特性，设计用户角度域分组和路径调度算法，在用户端合理分配角度域资源；在调度好的接收格点位置接收不同用户的导频和有效数据，使用最小二乘估计器估计其信道，并用最大比合并的方法进行数据恢复；使用上行传输中的用户分组和路径调度算法指导下行接入过程。本发明能降低多用户间的干扰，实现用户大规模接入，有效提升频谱利用率。

Description

一种基于大规模MIMO-OTFS的径分多址接入方法
技术领域
本发明属于无线通信技术领域，尤其涉及一种基于大规模MIMO-OTFS的径分多址接入方法。
背景技术
目前：大规模MIMO是下一代通信网络的一项有前途的技术。通过在基站上使用大量天线，可以实现在能量和容量效率方面的显着增益。利用准确的信道状态信息，可以分析系统的可达速率、量化网络干扰、优化资源分配并得出系统的能量效率。进而，可以优化大规模MIMO的信号处理单元，例如波束形成、预编码、用户调度和数据检测。但是，以当用户以较高的速度运动时，信道会随时间变化。在有限的信道相干时间内获取时变信道将更具挑战性。
为了突破以上瓶颈，研究者们提出了两类解决方法：基于自回归(AR)模型的方案和基于基扩展模型(BEM)的方案。对于时间选择性大规模MIMO信道，基于AR的方案中包括以下内容：(1)根据时变信道在相邻时间块内的关联性，建立了时变信道的AR模型；(2)开发基于期望最大化(EM)算法的稀疏贝叶斯学习(SBL)框架以学习时间和空间信道参数；(3)使用卡尔曼滤波进行上行时变信道跟踪；(4)利用用户信道的空间信息，对下行信道参数进行恢复并降维，跟踪下行时变信道。在双选择性大规模MIMO信道衰落情况下，基于BEM的方案中包括以下内容：(1)根据空间基扩展模型(SBEM)，提出了广义的空间基扩展模型(G-SBEM)以使用数个非正交基表示信道；(2)利用复指数BEM(CE-BEM)和提出的G-SBEM，建立稀疏信道估计模型；(3)提出两个基于压缩感知理论的算法：准块同时正交匹配跟踪(QBSO)算法和自适应QBSO算法，以恢复信道估计模型中的准块稀疏系数矩阵。
此外，一种被称为正交时频空间(OTFS)调制新型的二维调制技术应运而生。在传统的OFDM方案中，有效数据是在时频域上进行映射和处理的。因此，通常会采用插入导频的方式估计时频域内的时变信道。但是，OTFS提供了新的二维信号空间，即延迟-多普勒域。相应地，可在延迟-多普勒域上实现数据映射和解映射。有趣的是，时变信道可以用延迟-多普勒域上的几个常参数来等效描述。有研究通过考虑时变物理散射模型，推导了延迟-多普勒域上的显式输入输出模型。重要的是，证明了延迟-多普勒域上的等效信道是由每个散射路径的参数确定的，例如延迟和多普勒频率。另外一种在OTFS系统的上行多址方法被提出，其中为用户分配了延迟-多普勒域资源块，该资源块在域中以相等的间隔间隔开。针对时变大规模MIMO的多用户接入问题，对于信道估计，基于AR模型的框架通常假定大规模MIMO信道处于块衰落状态，并且块内的时变特性将丢失。尽管基于BEM的方案可以表征给定块内的大规模MIMO信道动态特征，但BEM的参数对用户移动性条件敏感。此外，两个模型都存在模型近似误差，特别是在高速条件下，误差将变得不可忽略。出现上述现象的原因是，这两种方案并非来自时变物理散射模型，并且无法捕获内在因素，从而有效地量化和描述时变信道。对于信号恢复，OTFS提供了时变信道的新观点。然而,OTFS的重要特征是，在延迟-多普勒域上，将以不同的双循环移位结构接收许多信号块。因此，低复杂度的数据检测方法不易设计。此外，在大规模MIMO-OTFS系统中接入大量用户，由于不同用户的信道具有不同的弥散特性，其用户间干扰管理问题较为复杂。若对每个用户等间隔的分配延迟-多普勒域中的资源块，又将造成频谱资源的严重浪费。
通过上述分析，现有技术存在的问题及缺陷为：
(1)针对时变大规模MIMO的多用户接入问题，对于信道估计，基于AR模型的框架通常假定大规模MIMO信道处于块衰落状态，并且块内的时变特性将丢失。尽管基于BEM的方案可以表征给定块内的大规模MIMO信道动态特征，但BEM的参数对用户移动性条件敏感。此外，两个模型都存在模型近似误差，特别是在高速条件下，误差将变得不可忽略。
(2)在大规模MIMO-OTFS系统中接入大量用户，由于不同用户的信道具有不同的弥散特性，其用户间干扰管理问题较为复杂；若对每个用户等间隔的分配延迟-多普勒域中的资源块，又将造成频谱资源的严重浪费。
解决以上问题及缺陷的难度为：基于大规模MIMO-OTFS的径分多址接入方法的难点在于：构建大规模MIMO-OTFS信道模型，并基于该模型，利用信道在角度-延迟-多普勒域的弥散特性，设计用户角度域分组和路径调度算法给接入的用户分配角度-延迟-多普勒域资源块，进行大规模MIMO-OTFS 3D信道跟踪和数据恢复，实现大量用户无用户间干扰的多址接入，即径分多址接入。
解决以上问题及缺陷的意义为：设计基于大规模MIMO-OTFS的新型径分多址接入方法，能面向高移动性场景，大规模地接入用户，减小用户间干扰，显著提高资源利用率和服务用户数，为下一代移动通信系统中超密集场景的高移动用户接入提供一种有效的解决思路。
发明内容
针对现有技术存在的问题，本发明提供了一种基于大规模MIMO-OTFS的径分多址接入方法。
本发明是这样实现的，一种基于大规模MIMO-OTFS的径分多址接入方法，所述基于大规模MIMO-OTFS的径分多址接入方法包括：
基于大规模均匀线性阵列天线，构建大规模MIMO时变上行信道模型；
基于OTFS调制，得到大规模MIMO-OTFS角度-延迟-多普勒域接收信号模型和角度-延迟-多普勒域上行主、次3D信道模型；
以信道模型为基础，利用大规模MIMO-OTFS信道在角度-延迟-多普勒域的弥散特性，设计用户角度域分组和路径调度算法，在用户端合理分配角度域资源，确保不同用户的观测区域在接收端的角度-延迟-多普勒域3D立方体区域上正交；
在调度好的接收格点位置接收不同用户的导频和有效数据，使用最小二乘LS估计器估计其信道，并用最大比合并MRC的方法进行数据恢复；
使用上行传输中的用户分组和路径调度算法指导下行接入过程，设计一种低复杂度的下行波束形成策略。
进一步，所述基于大规模MIMO-OTFS的径分多址接入方法的上行接入流程；首先，用户将数据映射在调度算法给本用户分配的延迟-多普勒域的格点上，经过OTFS调制转换到时频域发送出去；在基站端经过OTFS解调器，将接收信号变换到天线-延迟-多普勒域，再沿天线做DFT，将信号变换到角度-延迟-多普勒域；基站在期望接收的角度-延迟-多普勒域格点分别解映射并接收每个用户的信号，使用导频估计每个用户的角度-延迟-多普勒域信道，然后使用MRC恢复每个用户的数据；下行过程与上行类似，但在发送数据之前增加了波束形成向量的设计，最后在用户端估计等效信道并恢复数据。
进一步，所述基于大规模MIMO-OTFS的径分多址接入方法的大规模MIMO-OTFS信道模型的构建包括：
步骤一：考虑高移动场景下单小区毫米波大规模MIMO系统，基站服务K个随机分布的单天线用户。基站配备有Nr根天线的均匀线性阵列，且Nr＞＞K。空间中存在着散射体，用户信道由多条传播路径组成。无线信号能沿着视线路径到达基站，也可以由多个散射体反射，这意味着基站与用户之间的信道链路会受到频率选择性衰落的影响。由于用户的高移动性，信道变化很快，经历时间选择性衰落。假设特定用户和基站之间的信道有P条散射径，每条散射径对应着一个到达方向(DOA)、一个多普勒频移、一个时延。
定义θk，p为第k个用户第p条路径的DOA，相应的空间导向矢量可表示为：
Figure BDA0002616601930000041
其中λ为载波波长，天线间距d设置为半波长。因此，从几何信道模型，用户k在时隙n的时变信道表示为
Figure BDA0002616601930000042
其中l表示延迟域索引，hk，p、τk，p和νk，p分别表示第k个用户第p条路径信道增益、延迟和多普勒频率。δ(·)表示狄拉克函数，Ts是系统采样周期。因为Ts的分辨率足以在典型宽带系统沿着延迟域捕获详细的信道信息，所以假设τk，p＝nτ，pTs，nτ，p为整数。
步骤二：对第k个用户，将长度为LDND的数据序列重新排列成延迟-多普勒域的二维OTFS数据块
Figure BDA0002616601930000051
其中LD＝512和ND＝128分别表示子载波数和OFDM符号数；
在预处理块，应用ISFFT得到时频域数据块
Figure BDA0002616601930000052
其中
Figure BDA0002616601930000053
和
Figure BDA0002616601930000054
为归一化DFT矩阵；对
Figure BDA0002616601930000055
的每一列做LD点逆DFT得到传输信号块
Figure BDA0002616601930000056
其中sk，j表示一个OFDM符号；
通过给每个OFDM符号加循环前缀CP，本发明可以得到时域的一维传输信号
Figure BDA0002616601930000057
其中Lcp＝32为CP长度，sk将在NDT时间内占据LDΔf的带宽，其中Δf和T＝(Lcp+LD)Ts分别为子载波间隔和OFDM符号周期；
在基站侧，在后处理块连续地应用与用户侧对称的操作，如重新排列、去CP、LD点DFT和LD×ND维SFFT操作，在第nr个天线处，可以得到延迟-多普勒域二维接收数据块
Figure BDA0002616601930000058
Figure BDA0002616601930000059
的第(i，j+ND/2)元素表示为：
Figure BDA00026166019300000510
其中
Figure BDA00026166019300000511
为i-i′除以LD的余数，i＝0，1，...，LD-1，j＝-ND/2，...，0，...，ND/2-1，wi，j为均值为0，方差为
Figure BDA00026166019300000512
的复高斯噪声；
步骤三：利用步骤一中建立的大规模MIMO信道模型，将
Figure BDA00026166019300000513
的第(i，j+ND/2)元素重新表示为：
Figure BDA00026166019300000514
其中
Figure BDA0002616601930000061
且：
Figure BDA0002616601930000062
Figure BDA0002616601930000063
分别表示天线-延迟-多普勒域的主信道和次信道；此外l是基站接收信号格点的延迟域坐标，
Figure BDA0002616601930000064
步骤四：进一步考虑由大规模天线阵列带来的信道稀疏性，对信道按天线维度做DFT得到角度-延迟-多普勒域主、次信道模型：
Figure BDA0002616601930000065
Figure BDA0002616601930000066
其中
Figure BDA0002616601930000067
Figure BDA0002616601930000068
定义
Figure BDA0002616601930000069
通过对其做DFT可得角度-延迟-多普勒域接收信号模型：
Figure BDA00026166019300000610
进一步，所述基于大规模MIMO-OTFS的径分多址接入方法的角度-延迟-多普勒域的上下行径分多址接入方法包括：
步骤一：通过已有的算法在上行链路估计信道参数
Figure BDA00026166019300000611
步骤二：由大规模MIMO-OTFS角度-延迟-多普勒域信道模型可得，3D信道
Figure BDA0002616601930000071
和
Figure BDA0002616601930000072
在
Figure BDA0002616601930000073
内有主导值，其中：
Figure BDA0002616601930000074
(ik，p，jk，p，qk，p)与第k个用户的第p条散射径相关，认为是该路径的角度-延迟-多普勒域标签，第p条散射径几乎包含了格点(ik，p，jk，p，qk，p)处所有的能量，对于格点(ik，p，jk，p，qk，p)，有以下近似：
Figure BDA0002616601930000075
Figure BDA0002616601930000076
相应地，可导出以下关系式：
Figure BDA0002616601930000077
步骤三：对第k个用户，其角度标签集定义为
Figure BDA0002616601930000078
角度标签不重叠的用户被分配到同一组
Figure BDA0002616601930000079
满足：
Figure BDA00026166019300000710
其中
Figure BDA00026166019300000711
Dθ是角度域的保护间隔，对
Figure BDA00026166019300000712
分配给他们相同的延迟-多普勒域资源，即
Figure BDA00026166019300000713
和不同的接收端角度-延迟-多普勒域格点，即
Figure BDA00026166019300000714
步骤四：对不同用户组
Figure BDA00026166019300000715
给他们分配延迟-多普勒域可区分的资源满足以下限制：
Figure BDA00026166019300000716
其中
Figure BDA00026166019300000717
表示
Figure BDA00026166019300000718
保护间隔Dτ和Dν用于对抗3D信道
Figure BDA00026166019300000719
和
Figure BDA00026166019300000720
在延迟和多普勒域的弥散效应，定义延迟和多普勒域的最大弥散长度：
Figure BDA00026166019300000721
分别将Dτ和Dν设置为
Figure BDA00026166019300000722
和
Figure BDA00026166019300000723
调度后，不同的用户可以将其各自的数据映射到调度的延迟-多普勒域格点，同时在同一OTFS块中将数据发送到基站，并在基站处占用不同的3D资源；每个用户的有效数据块的延迟和多普勒域宽度分别为Wd和WD；基站并行地对不同用户的数据进行解映射和解码，而不会受到用户之间的干扰；
步骤五：第k个用户的发送数据块Xk中的导频的位置表示为
Figure BDA0002616601930000081
且
Figure BDA0002616601930000082
其在基站的观测信号
Figure BDA0002616601930000083
可写为：
Figure BDA0002616601930000084
其中
Figure BDA0002616601930000085
通过LS估计器，从上式恢复
Figure BDA0002616601930000086
得：
Figure BDA0002616601930000087
相应地，次信道：
Figure BDA0002616601930000088
为了加强数据恢复的性能，根据主信道
Figure BDA0002616601930000089
的近似表示和已知的
Figure BDA00026166019300000810
估计
Figure BDA00026166019300000811
得：
Figure BDA00026166019300000812
通过已获得的当前OTFS块内信道参数集
Figure BDA00026166019300000813
的准确信息，通过大规模MIMO-OTFS主、次信道模型重构所有格点的3D信道
Figure BDA00026166019300000814
和
Figure BDA00026166019300000815
步骤六：考虑用户发送有效数据块
Figure BDA00026166019300000816
内的一个数据符号，其元素索引为
Figure BDA00026166019300000817
由大规模MIMO-OTFS信道模型，导出基站处关于
Figure BDA00026166019300000818
的接收信息在集合
Figure BDA00026166019300000819
内，将数据符号
Figure BDA00026166019300000820
的P个接收格点列成一个P×1向量
Figure BDA00026166019300000821
得：
Figure BDA0002616601930000091
其中
Figure BDA0002616601930000092
是P×1噪声向量，式右边第三项为干扰的和；使用MRC的策略进行数据恢复，第k个用户第p个路径，即
Figure BDA0002616601930000093
的信干噪比表示为：
Figure BDA0002616601930000094
根据MRC原理，第p个路径的最优加权因子为：
Figure BDA0002616601930000095
定义权向量βk＝[βk，1，...，βk，P]T，并将P个接收信号用其本身的加权因子结合起来，MRC接收信号表示为：
Figure BDA0002616601930000096
可以导出合并信号的最优信噪比为
Figure BDA0002616601930000097
获得了最优信噪比，使用经典的算法恢复传输的数据，根据LS的特性，数据恢复的均方误差MSE为
Figure BDA0002616601930000098
步骤七：对于TDD系统，由于上下行信道在角度-延迟-多普勒域的互易性，下行信道参数
Figure BDA0002616601930000099
与上行的相同；上行的路径调度策略可用于下行的多用户服务，过与步骤二中相同的角度分组，第k个用户的接收信号主要来自于P层，即
Figure BDA0002616601930000101
相应地，接收信号
Figure BDA0002616601930000102
可写为：
Figure BDA0002616601930000103
其中Nk为噪声矩阵，其元素服从0均值、方差
Figure BDA0002616601930000104
的高斯分布。
不失一般性，将第k个用户的下行有效观测区域
Figure BDA0002616601930000105
设置为上行接入时分配的发送区域；在基站的3D传输资源空间内，分配给第k个用户的有效发送区域为
Figure BDA0002616601930000106
对给定的
Figure BDA0002616601930000107
其不同于
Figure BDA0002616601930000108
根据步骤三中的路径调度结果，确保：
Figure BDA0002616601930000109
在径分之后，第k个用户的格点
Figure BDA00026166019300001010
将观测到来自
Figure BDA00026166019300001011
的信号成分，其将分别经历自己的3D信道
Figure BDA00026166019300001012
在操作前，给出以下公式：
Figure BDA00026166019300001013
其中P×1向量bk，l，n表示对
Figure BDA00026166019300001014
中观测点(l，n)的波束形成操作，uk，l，n表示第k个用户期望的有效数据，得用户接收信号：
Figure BDA00026166019300001015
其中bk，l，n可用MMSE波束形成框架设计为：
Figure BDA00026166019300001016
参数
Figure BDA00026166019300001017
通过上下行参数的互易性获得，而
Figure BDA00026166019300001018
可通过最近的上行OTFS块获得，在波束形成后，
Figure BDA00026166019300001019
内所有格点的等效信道
Figure BDA00026166019300001020
相同，仅需
Figure BDA00026166019300001021
内发送一个格点的导频即可实施等效信道
Figure BDA00026166019300001022
的估计，估计方法与步骤四相同。
本发明的另一目的在于提供一种实施所述基于大规模MIMO-OTFS的径分多址接入方法的径分多址接入系统，所述基于大规模MIMO-OTFS的径分多址接入方法包括：
第一模型构建模块，用于基于大规模均匀线性阵列天线，构建大规模MIMO时变上行信道模型；
第二模型构建模块，用于基于OTFS调制，得到大规模MIMO-OTFS角度-延迟-多普勒域接收信号模型和角度-延迟-多普勒域上行主、次3D信道模型；
用户角度域分组和路径调度算法设计模块，用于以信道模型为基础，利用大规模MIMO-OTFS信道在角度-延迟-多普勒域的弥散特性，设计用户角度域分组和路径调度算法，在用户端合理分配角度域资源，确保不同用户的观测区域在接收端的角度-延迟-多普勒域3D立方体区域上正交；
数据恢复模块，用于在调度好的接收格点位置接收不同用户的导频和有效数据，使用最小二乘LS估计器估计其信道，并用最大比合并MRC的方法进行数据恢复；
下行波束形成策略模块，用于使用上行传输中的用户分组和路径调度算法指导下行接入过程，设计一种低复杂度的下行波束形成策略。
本发明的另一目的在于提供一种无线通信系统，所述无线通信系统搭载所述的径分多址接入系统。
本发明的另一目的在于提供一种移动通信系统，所述移动通信系统搭载所述的径分多址接入系统。
结合上述的所有技术方案，本发明所具备的优点及积极效果为：本发明基于大规模MIMO-OTFS的径分多址接入方法，大规模MIMO提供了角度维度的信号空间，而OTFS调制提供了延迟-多普勒域的二维信号空间。在高移动性场景下，利用大规模MIMO-OTFS多径信道的弥散特性，提出新型的径分多址接入技术。通过用户分组和路径调度，合理分配角度-延迟-多普勒域资源，实现用户的大规模接入。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对本申请实施例中所需要使用的附图做简单的介绍，显而易见地，下面所描述的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的基于大规模MIMO-OTFS的径分多址接入方法流程图。
图2是本发明实施例提供的基于大规模MIMO-OTFS的径分多址接入方法系统结构示意图；
图2中：1、第一模型构建模块；2、第二模型构建模块；3、用户角度域分组和路径调度算法设计模块；4、数据恢复模块；5、下行波束形成策略模块。
图3是本发明实施例提供的基于大规模MIMO-OTFS的径分多址接入方法上行部分的流程图。
图4是本发明实施例提供的基于大规模MIMO-OTFS的径分多址接入方法中上行大规模MIMO-OTFS信道角度-延迟多普勒域弥散特性示意图。
图5是本发明实施例提供的上行信道重构对比用户不同速度情况下的MSE仿真图。
图6是本发明实施例提供的上行数据恢复对比用户不同速度情况下、不同数据块的MSE仿真图。
图7是本发明实施例提供的下行数据恢复对比不同信道条件下的MSE仿真图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
针对现有技术存在的问题，本发明提供了一种基于大规模MIMO-OTFS的径分多址接入方法，下面结合附图对本发明作详细的描述。
如图1所示，本发明提供的径分多址接入方法包括以下步骤：
S101：基于大规模均匀线性阵列天线，构建大规模MIMO时变上行信道模型；
S102：基于OTFS调制，得到大规模MIMO-OTFS角度-延迟-多普勒域接收信号模型和角度-延迟-多普勒域上行主、次3D信道模型；
S103：以该信道模型为基础，利用大规模MIMO-OTFS信道在角度-延迟-多普勒域的弥散特性，设计用户角度域分组和路径调度算法，在用户端合理分配角度域资源，确保不同用户的观测区域在接收端的角度-延迟-多普勒域3D立方体区域上正交；
S104：在调度好的接收格点位置接收不同用户的导频和有效数据，使用最小二乘(LS)估计器估计其信道，并用最大比合并(MRC)的方法进行数据恢复；
S105：使用上行传输中的用户分组和路径调度算法指导下行接入过程，设计一种低复杂度的下行波束形成策略。
本发明提供的基于大规模MIMO-OTFS的径分多址接入方法业内的普通技术人员还可以采用其他的步骤实施，图1的本发明提供的基于大规模MIMO-OTFS的径分多址接入方法仅仅是一个具体实施例而已。
如图2所示，本发明提供的基于大规模MIMO-OTFS的径分多址接入系统包括：
第一模型构建模块1，用于基于大规模均匀线性阵列天线，构建大规模MIMO时变上行信道模型；
第二模型构建模块2，用于基于OTFS调制，得到大规模MIMO-OTFS角度-延迟-多普勒域接收信号模型和角度-延迟-多普勒域上行主、次3D信道模型；
用户角度域分组和路径调度算法设计模块3，用于以信道模型为基础，利用大规模MIMO-OTFS信道在角度-延迟-多普勒域的弥散特性，设计用户角度域分组和路径调度算法，在用户端合理分配角度域资源，确保不同用户的观测区域在接收端的角度-延迟-多普勒域3D立方体区域上正交；
数据恢复模块4，用于在调度好的接收格点位置接收不同用户的导频和有效数据，使用最小二乘(LS)估计器估计其信道，并用最大比合并(MRC)的方法进行数据恢复；
下行波束形成策略模块5，用于使用上行传输中的用户分组和路径调度算法指导下行接入过程，设计一种低复杂度的下行波束形成策略。
下面结合附图对本发明的技术方案作进一步的描述。
本发明具体涉及基于大规模MIMO-OTFS的径分多址接入方法，可应用于大规模MIMO无线网络，有效满足高移动性场景下网络的密集多用户接入需求，提升系统整体性能。对于上行链路通过用户角度域分组和路径调度分配角度-延迟-多普勒域资源，使用LS估计器实现当前OTFS块内角度-延迟-多普勒域3D信道的估计，并使用MRC的策略恢复数据。对于下行，使用与上行相同的路径调度策略，使用MMSE的方法设计波束形成向量，进而可实现下行等效信道估计和数据恢复。
下面结合附图对本发明的应用原理作详细的描述。
如图3所示，本发明实例提供的基于大规模MIMO-OTFS的径分多址接入方法包括以下步骤：
在大规模MIMO系统中使用OTFS调制，根据大规模均匀线性天线阵列和OTFS调制的特性，构建角度-延迟-多普勒域信道模型，得到上行信道，通过现有算法估计上行信道参数。
基站端依据大规模MIMO-OTFS信道在角度域的稀疏性和在延迟-多普勒域的弥散特性，使用路径调度算法，对用户进行空间分组，并给用户分配角度-延迟-多普勒域资源使其在基站端的接收区域正交。
每个用户在同一OTFS块内各自分配的延迟-多普勒域资源块内发送导频与数据，使用导频估计当前OTFS块内的角度-延迟-多普勒域信道，然后用其进行基于MRC的数据恢复。
使用上行路径调度算法指导下行，设计波束形成向量，进而可实现下行等效信道估计和数据恢复。
如图4所示，本发明实施例提供的基于大规模MIMO-OTFS的径分多址接入方法的上行接入流程。首先，用户将数据映射在调度算法给本用户分配的延迟-多普勒域的格点上，经过OTFS调制转换到时频域发送出去。在基站端经过OTFS解调器，将接收信号变换到天线-延迟-多普勒域，再沿天线做DFT，将信号变换到角度-延迟-多普勒域。基站在期望接收的角度-延迟-多普勒域格点分别解映射并接收每个用户的信号，使用导频估计每个用户的角度-延迟-多普勒域信道，然后使用MRC恢复每个用户的数据。下行过程与上行类似，但在发送数据之前增加了波束形成向量的设计，最后在用户端估计等效信道并恢复数据。
本发明实施例提供的基于大规模MIMO-OTFS的径分多址接入方法的应用场景为：
步骤一：考虑高移动场景下单小区毫米波大规模MIMO系统，基站服务K＝32个随机分布的单天线用户。基站配备有Nr＝128根天线的均匀线性阵列。空间中存在着散射体，用户信道由多条传播路径组成。无线信号能沿着视线路径到达基站，也可以由多个散射体反射，这意味着基站与用户之间的信道链路会受到频率选择性衰落的影响。用户移动速度范围为vs∈[120，306]km/h，最大多普勒频率为2kHz。由于用户的高移动性，信道变化很快，经历时间选择性衰落。假设特定用户和基站之间的信道有P条散射径，每条散射径对应着一个到达方向(DOA)、一个多普勒频移、一个时延。
定义θk，p为第k个用户第p条路径的DOA，相应的空间导向矢量可表示为：
Figure BDA0002616601930000151
其中λ为载波波长，天线间距d设置为半波长。因此，从几何信道模型，用户k在时隙n的时变信道表示为：
Figure BDA0002616601930000161
其中l表示延迟域索引，hk，p、τk，p和νk，p分别表示第k个用户第p条路径信道增益、延迟和多普勒频率。δ(·)表示狄拉克函数，
Figure BDA0002616601930000162
是系统采样周期。假设τk，p＝nτ，pTs，nτ，p为整数。τk，p分布在{0，Ts，2Ts，...，15Ts}中随机选择。θk，p和vk，p分别均匀分布在[-90°，90°]和[-2kHz，2kHz]内。
在本发明的优选实施例中，大规模MIMO-OTFS信道模型的构建包括：
步骤一：对第k个用户，将长度为LDND的数据序列重新排列成延迟-多普勒域的二维OTFS数据块
Figure BDA0002616601930000163
其中LD＝512和ND＝128分别表示子载波数和OFDM符号数。
在预处理块，应用ISFFT得到时频域数据块
Figure BDA0002616601930000164
其中
Figure BDA0002616601930000165
和
Figure BDA0002616601930000166
为归一化DFT矩阵。然后对
Figure BDA0002616601930000167
的每一列做LD点逆DFT得到传输信号块
Figure BDA0002616601930000168
其中sk，j表示一个OFDM符号。
通过给每个OFDM符号加循环前缀CP，本发明可以得到时域的一维传输信号
Figure BDA0002616601930000169
其中Lcp＝32为CP长度。sk将在NDT时间内占据LDΔf的带宽，其中Δf和T＝(Lcp+LD)Ts分别为子载波间隔和OFDM符号周期。
在基站侧，在后处理块连续地应用与用户侧对称的操作，如重新排列、去CP、LD点DFT和LD×ND维SFFT操作。相关地，在第nr个天线处，可以得到延迟-多普勒域二维接收数据块
Figure BDA00026166019300001610
的第(i，j+ND/2)元素可表示为：
Figure BDA00026166019300001611
其中
Figure BDA00026166019300001612
为i-i′除以LD的余数，i＝0，1，...，LD-1，j＝-ND/2，...，0，...，ND/2-1。此外，wi，j为均值为0，方差为
Figure BDA00026166019300001613
的复高斯噪声。
步骤二：利用步骤一中建立的大规模MIMO信道模型，将
Figure BDA00026166019300001614
的第(i，j+ND/2)元素重新表示为：
Figure BDA0002616601930000171
其中
Figure BDA0002616601930000172
且：
Figure BDA0002616601930000173
Figure BDA0002616601930000174
分别表示天线-延迟-多普勒域的主信道和次信道。此外l是基站接收信号格点的延迟域坐标，
Figure BDA0002616601930000175
步骤三：进一步考虑由大规模天线阵列带来的信道稀疏性，本发明可对信道按天线维度做DFT得到角度-延迟-多普勒域主、次信道模型：
Figure BDA0002616601930000176
Figure BDA0002616601930000177
其中
Figure BDA0002616601930000178
Figure BDA0002616601930000179
定义
Figure BDA00026166019300001710
通过对其做DFT可得角度-延迟-多普勒域接收信号模型：
Figure BDA00026166019300001711
在本发明的优选实施例中，角度-延迟-多普勒域的上下行径分多址接入方法包括：
步骤一：通过已有的算法在上行链路估计信道参数
Figure BDA0002616601930000181
步骤二：由大规模MIMO-OTFS角度-延迟-多普勒域信道模型可得，3D信道
Figure BDA0002616601930000182
和
Figure BDA0002616601930000183
在
Figure BDA0002616601930000184
内有主导值，其中：
Figure BDA0002616601930000185
(ik，p，jk，p，qk，p)与第k个用户的第p条散射径相关，可认为是该路径的角度-延迟-多普勒域标签。进一步地，第p条散射径几乎包含了格点(ik，p，jk，p，qk，p)处所有的能量。因此，对于格点(ik，p，jk，p，qk，p)，有以下近似：
Figure BDA0002616601930000186
Figure BDA0002616601930000187
相应地，可导出以下关系式：
Figure BDA0002616601930000188
步骤三：对第k个用户，其角度标签集定义为
Figure BDA0002616601930000189
角度标签不重叠的用户被分配到同一组
Figure BDA00026166019300001810
满足：
Figure BDA00026166019300001811
其中
Figure BDA00026166019300001812
Dθ是角度域的保护间隔。对
Figure BDA00026166019300001813
分配给他们相同的延迟-多普勒域资源，即
Figure BDA00026166019300001814
和不同的接收端角度-延迟-多普勒域格点，即
Figure BDA00026166019300001815
步骤四：对不同用户组
Figure BDA00026166019300001816
给他们分配延迟-多普勒域可区分的资源满足以下限制：
Figure BDA00026166019300001817
其中
Figure BDA00026166019300001818
表示
Figure BDA00026166019300001819
保护间隔Dτ和Dν用于对抗3D信道
Figure BDA00026166019300001820
和
Figure BDA00026166019300001821
在延迟和多普勒域的弥散效应。定义延迟和多普勒域的最大弥散长度：
Figure BDA00026166019300001822
典型地，分别将Dτ和Dν设置为
Figure BDA0002616601930000191
和
Figure BDA0002616601930000192
调度后，不同的用户可以将其各自的数据映射到调度的延迟-多普勒域格点，同时在同一OTFS块中将数据发送到基站，并在基站处占用不同的3D资源。每个用户的有效数据块的延迟和多普勒域宽度分别为Wd＝48和WD＝50。然后，基站可以并行地对不同用户的数据进行解映射和解码，而不会受到用户之间的干扰。
步骤五：第k个用户的发送数据块Xk中的导频的位置表示为
Figure BDA0002616601930000193
且
Figure BDA0002616601930000194
其在基站的观测信号
Figure BDA0002616601930000195
可写为：
Figure BDA0002616601930000196
其中
Figure BDA0002616601930000197
通过LS估计器，可从上式恢复
Figure BDA0002616601930000198
得：
Figure BDA0002616601930000199
相应地，次信道：
Figure BDA00026166019300001910
为了加强数据恢复的性能，根据主信道
Figure BDA00026166019300001911
的近似表示和已知的
Figure BDA00026166019300001912
可估计
Figure BDA00026166019300001913
得：
Figure BDA00026166019300001914
通过已获得的当前OTFS块内信道参数集
Figure BDA00026166019300001915
的准确信息，可通过大规模MIMO-OTFS主、次信道模型重构所有格点的3D信道
Figure BDA00026166019300001916
和
Figure BDA00026166019300001917
步骤六：考虑用户发送有效数据块
Figure BDA00026166019300001918
内的一个数据符号，其元素索引为
Figure BDA00026166019300001919
由大规模MIMO-OTFS信道模型，可导出基站处关于
Figure BDA00026166019300001920
的接收信息在集合
Figure BDA00026166019300001921
内。然后，将数据符号
Figure BDA00026166019300001922
的P个接收格点列成一个P×1向量
Figure BDA00026166019300001923
得：
Figure BDA0002616601930000201
其中
Figure BDA0002616601930000202
是P×1噪声向量，式右边第三项为干扰的和。使用MRC的策略进行数据恢复。第k个用户第p个路径，即
Figure BDA0002616601930000203
的信干噪比可表示为：
Figure BDA0002616601930000204
根据MRC原理，第p个路径的最优加权因子为：
Figure BDA0002616601930000205
定义权向量βk＝[βk，1，...，βk，P]T，并将P个接收信号用其本身的加权因子结合起来，MRC接收信号可表示为：
Figure BDA0002616601930000206
可以导出合并信号的最优信噪比为
Figure BDA0002616601930000207
获得了最优信噪比，可以使用经典的算法恢复传输的数据，如LS估计器。根据LS的特性，数据恢复的均方误差(meansquareerror，MSE)为
Figure BDA0002616601930000208
步骤七：对于TDD系统，由于上下行信道在角度-延迟-多普勒域的互易性，下行信道参数
Figure BDA0002616601930000209
与上行的相同。进一步地，上行的路径调度策略可用于下行的多用户服务。通过与步骤二中相同的角度分组，第k个用户的接收信号主要来自于P层，即
Figure BDA00026166019300002010
相应地，接收信号
Figure BDA00026166019300002011
可写为：
Figure BDA0002616601930000211
其中Nk为噪声矩阵，其元素服从0均值、方差
Figure BDA0002616601930000212
的高斯分布。
不失一般性，将第k个用户的下行有效观测区域
Figure BDA0002616601930000213
设置为上行接入时分配的发送区域。在基站的3D传输资源空间内，分配给第k个用户的有效发送区域为
Figure BDA0002616601930000214
对给定的
Figure BDA0002616601930000215
其不同于
Figure BDA0002616601930000216
根据步骤三中的路径调度结果，可确保:
Figure BDA0002616601930000217
在径分之后，第k个用户的格点
Figure BDA0002616601930000218
将观测到来自
Figure BDA0002616601930000219
的信号成分，其将分别经历自己的3D信道
Figure BDA00026166019300002110
在操作前，给出以下公式：
Figure BDA00026166019300002111
其中P×1向量bk，l，n表示对
Figure BDA00026166019300002112
中观测点(l，n)的波束形成操作，uk，l，n表示第k个用户期望的有效数据。进而可得用户接收信号：
Figure BDA00026166019300002113
其中bk，l，n可用MMSE波束形成框架设计为：
Figure BDA00026166019300002114
注意到，参数
Figure BDA00026166019300002115
通过上下行参数的互易性获得，而
Figure BDA00026166019300002116
可通过最近的上行OTFS块获得。进一步地，在波束形成后，
Figure BDA00026166019300002117
内所有格点的等效信道
Figure BDA00026166019300002118
相同。因此，仅需
Figure BDA00026166019300002119
内发送一个格点的导频即可实施等效信道
Figure BDA00026166019300002120
的估计，估计方法与步骤四相同。
图5(a)和图5(b)展示了较低速度和较高速度的信道重构性能。带三角形的曲线是在不考虑次信道的情况下估计的三维主信道，带圆曲线和带叉曲线代表估计的主信道和次信道。可以看出，不考虑次信道的主信道的MSE比其他信道的MSE要高得多。此外，在高移动性场景下，其性能变差，而其他两条曲线与低移动性场景下的曲线更接近。这种特性可以解释为次信道的影响随着多普勒频移频率的增大而增大。在高移动性情况下，多普勒频移变大，因此，三角形曲线变差。对于其他两条曲线，此外，由于估计方程的近似形式，次信道的MSE总是高于主信道的MSE。
图6描述了两种移动情况下的数据恢复性能，其中考虑了有效数据块的稀疏性。正如预期的那样，MSE随信噪比的增加而减小，并趋于收敛；后一种情况是由于数据块内的干扰。当信噪比很低时，与噪声相比，干扰很小。但在高信噪比情况下，干扰远大于噪声，严重影响数据恢复性能。数据块稀疏调度的MSE与数据块密集分布的MSE之间的差距值得注意。这是由于主弥散接收格点的距离越远，由能量泄露引起的干扰就越小。
图7显示了不同信道质量和用户速度下下行数据恢复性能与信噪比的对比。可以看到，完美信道下的MSE比不匹配信道下的MSE要低得多，而不匹配信道下的MSE是由于3D信道估计不准确造成的。此外，在高信噪比下，MSE趋于收敛，原因与图6相同。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种基于大规模MIMO-OTFS的径分多址接入方法，其特征在于，所述基于大规模MIMO-OTFS的径分多址接入方法包括：
基于大规模均匀线性阵列天线，构建大规模MIMO时变上行信道模型；
基于OTFS调制，得到大规模MIMO-OTFS角度-延迟-多普勒域接收信号模型和角度-延迟-多普勒域上行主、次3D信道模型；
以信道模型为基础，利用大规模MIMO-OTFS信道在角度-延迟-多普勒域的弥散特性，设计用户角度域分组和路径调度算法，在用户端合理分配角度域资源，确保不同用户的观测区域在接收端的角度-延迟-多普勒域3D立方体区域上正交；
在调度好的接收格点位置接收不同用户的导频和有效数据，使用最小二乘LS估计器估计其信道，并用最大比合并MRC的方法进行数据恢复；
使用上行传输中的用户分组和路径调度算法指导下行接入过程，设计一种低复杂度的下行波束形成策略。
2.如权利要求1所述的基于大规模MIMO-OTFS的径分多址接入方法，其特征在于，所述基于大规模MIMO-OTFS的径分多址接入方法的上行接入流程；首先，用户将数据映射在调度算法给本用户分配的延迟-多普勒域的格点上，经过OTFS调制转换到时频域发送出去；在基站端经过OTFS解调器，将接收信号变换到天线-延迟-多普勒域，再沿天线做DFT，将信号变换到角度-延迟-多普勒域；基站在期望接收的角度-延迟-多普勒域格点分别解映射并接收每个用户的信号，使用导频估计每个用户的角度-延迟-多普勒域信道，然后使用MRC恢复每个用户的数据；下行过程与上行类似，但在发送数据之前增加了波束形成向量的设计，最后在用户端估计等效信道并恢复数据。
3.如权利要求1所述的基于大规模MIMO-OTFS的径分多址接入方法，其特征在于，所述基于大规模MIMO-OTFS的径分多址接入方法的大规模MIMO-OTFS信道模型的构建包括：
步骤一：考虑高移动场景下单小区毫米波大规模MIMO系统，基站服务K个随机分布的单天线用户；基站配备有Nr根天线的均匀线性阵列，且Nr＞＞K；空间中存在着散射体，用户信道由多条传播路径组成；无线信号能沿着视线路径到达基站，也可以由多个散射体反射，这意味着基站与用户之间的信道链路会受到频率选择性衰落的影响；由于用户的高移动性，信道变化很快，经历时间选择性衰落；假设特定用户和基站之间的信道有P条散射径，每条散射径对应着一个到达方向(DOA)、一个多普勒频移、一个时延；
定义θk，p为第k个用户第p条路径的DOA，相应的空间导向矢量可表示为：
Figure FDA0002616601920000021
其中λ为载波波长，天线间距d设置为半波长；因此，从几何信道模型，用户k在时隙n的时变信道表示为
Figure FDA0002616601920000022
其中l表示延迟域索引，hk，p、τk，p和
Figure FDA0002616601920000023
分别表示第k个用户第p条路径信道增益、延迟和多普勒频率；δ(·)表示狄拉克函数，Ts是系统采样周期；因为Ts的分辨率足以在典型宽带系统沿着延迟域捕获详细的信道信息，所以假设τk，p＝nτ，pTs，nτ，p为整数；
步骤二：对第k个用户，将长度为LDND的数据序列重新排列成延迟-多普勒域的二维OTFS数据块
Figure FDA0002616601920000024
其中LD＝512和ND＝128分别表示子载波数和OFDM符号数；
在预处理块，应用ISFFT得到时频域数据块
Figure FDA0002616601920000025
其中
Figure FDA0002616601920000026
和
Figure FDA0002616601920000027
为归一化DFT矩阵；对
Figure FDA0002616601920000028
的每一列做LD点逆DFT得到传输信号块
Figure FDA0002616601920000029
其中sk，j表示一个OFDM符号；
通过给每个OFDM符号加循环前缀CP，本发明可以得到时域的一维传输信号
Figure FDA00026166019200000210
其中Lcp＝32为CP长度，sk将在NDT时间内占据LDΔf的带宽，其中Δf和T＝(Lcp+LD)Ts分别为子载波间隔和OFDM符号周期；
在基站侧，在后处理块连续地应用与用户侧对称的操作，如重新排列、去CP、LD点DFT和LD×ND维SFFT操作，在第nr个天线处，可以得到延迟-多普勒域二维接收数据块
Figure FDA0002616601920000031
Figure FDA00026166019200000314
的第(i，j+ND/2)元素表示为：
Figure FDA0002616601920000032
其中
Figure FDA0002616601920000033
为i-i′除以LD的余数，i＝0，1，...，LD-1，j＝-ND/2，...，0，...，ND/2-1，wi，j为均值为0，方差为
Figure FDA00026166019200000313
的复高斯噪声；
步骤三：利用步骤一中建立的大规模MIMO信道模型，将
Figure FDA00026166019200000315
的第(i，j+ND/2)元素重新表示为：
Figure FDA0002616601920000035
其中
Figure FDA0002616601920000036
且：
Figure FDA0002616601920000037
Figure FDA0002616601920000038
分别表示天线-延迟-多普勒域的主信道和次信道；此外l是基站接收信号格点的延迟域坐标，
Figure FDA0002616601920000039
步骤四：进一步考虑由大规模天线阵列带来的信道稀疏性，对信道按天线维度做DFT得到角度-延迟-多普勒域主、次信道模型：
Figure FDA00026166019200000310
Figure FDA00026166019200000311
其中
Figure FDA00026166019200000312
Figure FDA0002616601920000041
定义
Figure FDA0002616601920000042
通过对其做DFT可得角度-延迟-多普勒域接收信号模型：
Figure FDA0002616601920000043
4.如权利要求1所述的基于大规模MIMO-OTFS的径分多址接入方法，其特征在于，所述基于大规模MIMO-OTFS的径分多址接入方法的角度-延迟-多普勒域的上下行径分多址接入方法包括：
步骤一：通过已有的算法在上行链路估计信道参数
Figure FDA0002616601920000044
步骤二：由大规模MIMO-OTFS角度-延迟-多普勒域信道模型可得，3D信道
Figure FDA0002616601920000045
和
Figure FDA0002616601920000046
在
Figure FDA0002616601920000047
内有主导值，其中：
Figure FDA0002616601920000048
(ik，p，jk，p，qk，p)与第k个用户的第p条散射径相关，认为是该路径的角度-延迟-多普勒域标签，第p条散射径几乎包含了格点(ik，p，jk，p，qk，p)处所有的能量，对于格点(ik，p，jk，p，qk，p)，有以下近似：
Figure FDA0002616601920000049
Figure FDA00026166019200000410
相应地，可导出以下关系式：
Figure FDA00026166019200000411
步骤三：对第k个用户，其角度标签集定义为
Figure FDA00026166019200000412
角度标签不重叠的用户被分配到同一组
Figure FDA00026166019200000413
满足：
Figure FDA00026166019200000414
其中
Figure FDA0002616601920000051
Dθ是角度域的保护间隔，对k1，
Figure FDA0002616601920000052
分配给他们相同的延迟-多普勒域资源，即
Figure FDA0002616601920000053
和不同的接收端角度-延迟-多普勒域格点，即
Figure FDA0002616601920000054
k1≠k2；
步骤四：对不同用户组
Figure FDA0002616601920000055
给他们分配延迟-多普勒域可区分的资源满足以下限制：
Figure FDA0002616601920000056
其中
Figure FDA0002616601920000057
表示
Figure FDA0002616601920000058
保护间隔Dτ和
Figure FDA0002616601920000059
用于对抗3D信道
Figure FDA00026166019200000510
和
Figure FDA00026166019200000511
在延迟和多普勒域的弥散效应，定义延迟和多普勒域的最大弥散长度：
Figure FDA00026166019200000512
分别将Dτ和
Figure FDA00026166019200000513
设置为
Figure FDA00026166019200000514
和
Figure FDA00026166019200000515
调度后，不同的用户可以将其各自的数据映射到调度的延迟-多普勒域格点，同时在同一OTFS块中将数据发送到基站，并在基站处占用不同的3D资源；每个用户的有效数据块的延迟和多普勒域宽度分别为Wd和WD；基站并行地对不同用户的数据进行解映射和解码，而不会受到用户之间的干扰；
步骤五：第k个用户的发送数据块Xk中的导频的位置表示为
Figure FDA00026166019200000516
且
Figure FDA00026166019200000517
其在基站的观测信号
Figure FDA00026166019200000518
可写为：
Figure FDA00026166019200000519
其中
Figure FDA00026166019200000520
通过LS估计器，从上式恢复
Figure FDA00026166019200000521
得：
Figure FDA00026166019200000522
相应地，次信道：
Figure FDA0002616601920000061
为了加强数据恢复的性能，根据主信道
Figure FDA0002616601920000062
的近似表示和已知的
Figure FDA0002616601920000063
估计
Figure FDA0002616601920000064
得：
Figure FDA0002616601920000065
通过已获得的当前OTFS块内信道参数集
Figure FDA0002616601920000066
的准确信息，通过大规模MIMO-OTFS主、次信道模型重构所有格点的3D信道
Figure FDA0002616601920000067
和
Figure FDA0002616601920000068
步骤六：考虑用户发送有效数据块
Figure FDA0002616601920000069
内的一个数据符号，其元素索引为
Figure FDA00026166019200000610
由大规模MIMO-OTFS信道模型，导出基站处关于
Figure FDA00026166019200000611
的接收信息在集合
Figure FDA00026166019200000612
内，将数据符号
Figure FDA00026166019200000613
的P个接收格点列成一个P×1向量
Figure FDA00026166019200000614
得：
Figure FDA00026166019200000615
其中
Figure FDA00026166019200000616
是P×1噪声向量，式右边第三项为干扰的和；使用MRC的策略进行数据恢复，第k个用户第p个路径，即
Figure FDA00026166019200000617
的信干噪比表示为：
Figure FDA00026166019200000618
根据MRC原理，第p个路径的最优加权因子为：
Figure FDA0002616601920000071
定义权向量βk＝[βk，1，...，βk，P]T，并将P个接收信号用其本身的加权因子结合起来，MRC接收信号表示为：
Figure FDA0002616601920000072
可以导出合并信号的最优信噪比为
Figure FDA00026166019200000718
获得了最优信噪比，使用经典的算法恢复传输的数据，根据LS的特性，数据恢复的均方误差MSE为
Figure FDA0002616601920000074
步骤七：对于TDD系统，由于上下行信道在角度-延迟-多普勒域的互易性，下行信道参数
Figure FDA0002616601920000075
与上行的相同；上行的路径调度策略可用于下行的多用户服务，过与步骤二中相同的角度分组，第k个用户的接收信号主要来自于P层，即
Figure FDA0002616601920000076
相应地，接收信号
Figure FDA0002616601920000077
可写为：
Figure FDA0002616601920000078
其中Nk为噪声矩阵，其元素服从0均值、方差
Figure FDA0002616601920000079
的高斯分布；
不失一般性，将第k个用户的下行有效观测区域
Figure FDA00026166019200000710
设置为上行接入时分配的发送区域；在基站的3D传输资源空间内，分配给第k个用户的有效发送区域为
Figure FDA00026166019200000711
对给定的
Figure FDA00026166019200000712
其不同于
Figure FDA00026166019200000713
根据步骤三中的路径调度结果，确保：
Figure FDA00026166019200000714
在径分之后，第k个用户的格点
Figure FDA00026166019200000715
将观测到来自
Figure FDA00026166019200000716
的信号成分，其将分别经历自己的3D信道
Figure FDA00026166019200000717
在操作前，给出以下公式：
Figure FDA0002616601920000081
其中P×1向量bk，l，n表示对
Figure FDA0002616601920000082
中观测点(l，n)的波束形成操作，uk，l，n表示第k个用户期望的有效数据，得用户接收信号：
Figure FDA0002616601920000083
其中bk，l，n可用MMSE波束形成框架设计为：
Figure FDA0002616601920000084
参数
Figure FDA0002616601920000085
通过上下行参数的互易性获得，而
Figure FDA0002616601920000086
可通过最近的上行OTFS块获得，在波束形成后，
Figure FDA0002616601920000087
内所有格点的等效信道
Figure FDA0002616601920000088
相同，仅需
Figure FDA0002616601920000089
内发送一个格点的导频即可实施等效信道
Figure FDA00026166019200000810
的估计，估计方法与步骤四相同。
5.一种实施权利要求1～4任意一项所述基于大规模MIMO-OTFS的径分多址接入方法的径分多址接入系统，其特征在于，所述径分多址接入系统包括：
第一模型构建模块，用于基于大规模均匀线性阵列天线，构建大规模MIMO时变上行信道模型；
第二模型构建模块，用于基于OTFS调制，得到大规模MIMO-OTFS角度-延迟-多普勒域接收信号模型和角度-延迟-多普勒域上行主、次3D信道模型；
用户角度域分组和路径调度算法设计模块，用于以信道模型为基础，利用大规模MIMO-OTFS信道在角度-延迟-多普勒域的弥散特性，设计用户角度域分组和路径调度算法，在用户端合理分配角度域资源，确保不同用户的观测区域在接收端的角度-延迟-多普勒域3D立方体区域上正交；
数据恢复模块，用于在调度好的接收格点位置接收不同用户的导频和有效数据，使用最小二乘LS估计器估计其信道，并用最大比合并MRC的方法进行数据恢复；
下行波束形成策略模块，用于使用上行传输中的用户分组和路径调度算法指导下行接入过程，设计一种低复杂度的下行波束形成策略。
6.一种无线通信系统，其特征在于，所述无线通信系统搭载权利要求5所述的径分多址接入系统。
CN202010771007.9A 2020-08-04 2020-08-04 一种基于大规模mimo-otfs的径分多址接入方法 Active CN112087247B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010771007.9A CN112087247B (zh) 2020-08-04 2020-08-04 一种基于大规模mimo-otfs的径分多址接入方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010771007.9A CN112087247B (zh) 2020-08-04 2020-08-04 一种基于大规模mimo-otfs的径分多址接入方法

Publications (2)

Publication Number Publication Date
CN112087247A true CN112087247A (zh) 2020-12-15
CN112087247B CN112087247B (zh) 2021-10-08

Family

ID=73735510

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010771007.9A Active CN112087247B (zh) 2020-08-04 2020-08-04 一种基于大规模mimo-otfs的径分多址接入方法

Country Status (1)

Country Link
CN (1) CN112087247B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112882005A (zh) * 2021-01-14 2021-06-01 东南大学 一种基于贝叶斯学习的otfs雷达目标参数估计方法
CN112929316A (zh) * 2021-01-25 2021-06-08 南京邮电大学 基于otfs调制的交错式时频多址方式调制解调方法及装置
CN113098818A (zh) * 2021-04-01 2021-07-09 北京交通大学 一种正交扩频数据的交织与映射方法
CN113114393A (zh) * 2021-03-26 2021-07-13 华南理工大学 针对基于索引调制的环境反向散射通信在多接入信道下的低复杂度检测方法、读写器和系统
CN113286355A (zh) * 2021-05-20 2021-08-20 南京邮电大学 基于otfs-noma跨域传输系统中的功率分配方法
CN113507426A (zh) * 2021-06-16 2021-10-15 北京邮电大学 基于otfs调制的联合信道估计与信号检测方法及装置
CN113556306A (zh) * 2021-07-19 2021-10-26 上海交通大学 离散傅里叶变换扩展正交时频空调制方法
CN113890796A (zh) * 2021-09-30 2022-01-04 成都工业学院 基于otfs系统调制解调的高速信道估计装置及方法
WO2022171084A1 (zh) * 2021-02-10 2022-08-18 维沃移动通信有限公司 接入方法、装置、通信设备及可读存储介质
CN115002793A (zh) * 2022-06-17 2022-09-02 东南大学 均衡用户角度间距的毫米波大规模mimo系统用户分组方法
CN115051894A (zh) * 2022-06-10 2022-09-13 西安电子科技大学 基于串行干扰消除的otfs系统符号检测方法
WO2023274126A1 (zh) * 2021-06-30 2023-01-05 维沃移动通信有限公司 传输处理方法、装置、通信设备及可读存储介质
CN116094876A (zh) * 2023-03-10 2023-05-09 南京邮电大学 一种基于非对称架构的正交时频空系统的信道估计方法
CN116094875A (zh) * 2023-02-09 2023-05-09 重庆邮电大学 一种超大规模mimo系统中基于上行辅助的otfs下行信道估计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1691562A (zh) * 2004-04-19 2005-11-02 北京三星通信技术研究有限公司 基于可变扩频和跳频的正交时频域的码分多址接入方法
WO2018121460A1 (zh) * 2016-12-26 2018-07-05 华为技术有限公司 非正交传输数据的方法和设备
CN108964725A (zh) * 2018-07-20 2018-12-07 西安电子科技大学 时变大规模mimo网络中信道参数的稀疏估计方法
CN110024343A (zh) * 2016-11-30 2019-07-16 瑞典爱立信有限公司 用于多用户mimo系统中参考信号的自适应子载波间隔
EP3610582A1 (en) * 2017-04-11 2020-02-19 Cohere Technologies, Inc. Digital communication using dispersed orthogonal time frequency space modulated signals

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1691562A (zh) * 2004-04-19 2005-11-02 北京三星通信技术研究有限公司 基于可变扩频和跳频的正交时频域的码分多址接入方法
CN110024343A (zh) * 2016-11-30 2019-07-16 瑞典爱立信有限公司 用于多用户mimo系统中参考信号的自适应子载波间隔
WO2018121460A1 (zh) * 2016-12-26 2018-07-05 华为技术有限公司 非正交传输数据的方法和设备
EP3610582A1 (en) * 2017-04-11 2020-02-19 Cohere Technologies, Inc. Digital communication using dispersed orthogonal time frequency space modulated signals
CN108964725A (zh) * 2018-07-20 2018-12-07 西安电子科技大学 时变大规模mimo网络中信道参数的稀疏估计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YUSHAN LIU: "Uplink-Aided High Mobility Downlink Channel Estimation Over Massive MIMO-OTFS System", 《IEEE》 *
李建东等: "面向5G无线网络的智能干扰管理技术", 《电信科学》 *
杨维等: "一种改进基于门限的稀疏码多址接入低复杂度多用户检测算法", 《电子与信息学报》 *

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112882005A (zh) * 2021-01-14 2021-06-01 东南大学 一种基于贝叶斯学习的otfs雷达目标参数估计方法
CN112882005B (zh) * 2021-01-14 2024-01-26 东南大学 一种基于贝叶斯学习的otfs雷达目标参数估计方法
CN112929316A (zh) * 2021-01-25 2021-06-08 南京邮电大学 基于otfs调制的交错式时频多址方式调制解调方法及装置
WO2022171084A1 (zh) * 2021-02-10 2022-08-18 维沃移动通信有限公司 接入方法、装置、通信设备及可读存储介质
CN113114393B (zh) * 2021-03-26 2022-04-22 华南理工大学 针对基于索引调制的环境反向散射通信在多接入信道下的低复杂度检测方法、读写器和系统
CN113114393A (zh) * 2021-03-26 2021-07-13 华南理工大学 针对基于索引调制的环境反向散射通信在多接入信道下的低复杂度检测方法、读写器和系统
CN113098818A (zh) * 2021-04-01 2021-07-09 北京交通大学 一种正交扩频数据的交织与映射方法
CN113098818B (zh) * 2021-04-01 2022-04-22 北京交通大学 一种正交扩频数据的交织与映射方法
CN113286355A (zh) * 2021-05-20 2021-08-20 南京邮电大学 基于otfs-noma跨域传输系统中的功率分配方法
CN113286355B (zh) * 2021-05-20 2022-07-29 南京邮电大学 基于otfs-noma跨域传输系统中的功率分配方法
CN113507426A (zh) * 2021-06-16 2021-10-15 北京邮电大学 基于otfs调制的联合信道估计与信号检测方法及装置
WO2023274126A1 (zh) * 2021-06-30 2023-01-05 维沃移动通信有限公司 传输处理方法、装置、通信设备及可读存储介质
CN113556306B (zh) * 2021-07-19 2022-04-08 上海交通大学 离散傅里叶变换扩展正交时频空调制方法及系统
CN113556306A (zh) * 2021-07-19 2021-10-26 上海交通大学 离散傅里叶变换扩展正交时频空调制方法
CN113890796B (zh) * 2021-09-30 2022-08-12 成都工业学院 基于otfs系统调制解调的高速信道估计装置及方法
CN113890796A (zh) * 2021-09-30 2022-01-04 成都工业学院 基于otfs系统调制解调的高速信道估计装置及方法
CN115051894A (zh) * 2022-06-10 2022-09-13 西安电子科技大学 基于串行干扰消除的otfs系统符号检测方法
CN115051894B (zh) * 2022-06-10 2024-02-06 西安电子科技大学 基于串行干扰消除的otfs系统符号检测方法
CN115002793A (zh) * 2022-06-17 2022-09-02 东南大学 均衡用户角度间距的毫米波大规模mimo系统用户分组方法
CN115002793B (zh) * 2022-06-17 2024-01-02 东南大学 均衡用户角度间距的毫米波大规模mimo系统用户分组方法
CN116094875A (zh) * 2023-02-09 2023-05-09 重庆邮电大学 一种超大规模mimo系统中基于上行辅助的otfs下行信道估计方法
CN116094876A (zh) * 2023-03-10 2023-05-09 南京邮电大学 一种基于非对称架构的正交时频空系统的信道估计方法
CN116094876B (zh) * 2023-03-10 2023-08-29 南京邮电大学 一种基于非对称架构的正交时频空系统的信道估计方法

Also Published As

Publication number Publication date
CN112087247B (zh) 2021-10-08

Similar Documents

Publication Publication Date Title
CN112087247B (zh) 一种基于大规模mimo-otfs的径分多址接入方法
CN112511201B (zh) 天波大规模mimo通信方法及模型和系统
CN109691050B (zh) 一种ofdm系统中频偏估计方法、装置及系统
CN103166685A (zh) Lte中一种基于联合功率分配的干扰对齐方法
WO2017121095A1 (zh) 一种实现测量参考符号传输的方法及装置
US11477053B2 (en) Method and apparatus for estimating channel in wireless communication system
Bartolome et al. MMSE techniques for space diversity receivers in OFDM-based wireless LANs
CN102571659A (zh) 一种干扰噪声估计和干扰抑制方法及相应系统
CN110881010B (zh) 统计csi辅助的多用户noma下行传输方法
CN102891815A (zh) 一种时分双工多基站协作系统中的低复杂度信道估计方法
Ge et al. OTFS signaling for SCMA with coordinated multi-point vehicle communications
CN115065432A (zh) 天波大规模mimo三重波束基信道建模及信道信息获取
CN102215073B (zh) 多点—多用户上行低复杂度mimo检测方法
CN113645165B (zh) 5g下行链路的分组插值-加权合并信道估计方法及系统
Li et al. Multiple access for massive MIMO-OTFS networks over angle-delay-doppler domain
Yin et al. Design and performance analysis of AFDM with multiple antennas in doubly selective channels
CN103220025A (zh) 一种应用于vmimo系统的多用户配对算法
Ji et al. Multipath extraction based UL/DL channel estimation for FDD massive MIMO-OFDM systems
CN109194375B (zh) Fd-mimo多小区下行干扰协调方法
Darsena et al. Reception strategies for sky-ground uplink non-orthogonal multiple access
CN113938360A (zh) 一种基于指纹定位的分布式mimo系统协方差矩阵估计方法
Shao et al. Angle-delay-doppler domain NOMA over massive MIMO-OTFS networks
CN101958874B (zh) 基于角分复用的双正交频分多址蜂窝系统
KR101234245B1 (ko) 다중 안테나를 사용하는 직교 주파수 분할 시스템에서의 블라인드 채널 추정을 위한 선형 프리코딩과 통신 용량을 최대화시키는 프리코딩과의 호환 방법
CN115996158A (zh) 一种宽带毫米波大规模mimo波束空间信道估计方法

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