CN110213185B - 一种基于原子范数最小化的三维信道参数估计方法 - Google Patents

一种基于原子范数最小化的三维信道参数估计方法 Download PDF

Info

Publication number
CN110213185B
CN110213185B CN201910483315.9A CN201910483315A CN110213185B CN 110213185 B CN110213185 B CN 110213185B CN 201910483315 A CN201910483315 A CN 201910483315A CN 110213185 B CN110213185 B CN 110213185B
Authority
CN
China
Prior art keywords
channel
matrix
anm
dimensional
frequency
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
Application number
CN201910483315.9A
Other languages
English (en)
Other versions
CN110213185A (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.)
Shanghai Zhuangyan Automation Technology Co ltd
Original Assignee
Beijing Jiaotong 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 Beijing Jiaotong University filed Critical Beijing Jiaotong University
Priority to CN201910483315.9A priority Critical patent/CN110213185B/zh
Publication of CN110213185A publication Critical patent/CN110213185A/zh
Application granted granted Critical
Publication of CN110213185B publication Critical patent/CN110213185B/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
    • 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
    • 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
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2602—Signal structure

Landscapes

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

Abstract

本发明属于通信技术领域，具体涉及一种基于原子范数最小化的三维信道参数估计方法。所述方法根据基站接收到的导频信号和已知的移动台发送的导频信号进行信道参数的估计；已知信道矩阵H[k]中包含了三个拥有范德蒙结构的参数：出发角度AoA、到达角度AoD、和时间延迟，根据毫米波信道的稀疏散射特性，将***的信道估计问题可以被看做三维的线性谱估计问题，利用ANM的方法以无网格的方式高精度恢复出以上参数；并且利用参数之间的独立性，将三维的ANM信道估计问题分解为二维ANM和一维ANM的估计问题，以降低三维ANM信道估计的算法复杂度。

Description

一种基于原子范数最小化的三维信道参数估计方法
技术领域
本发明属于通信技术领域，具体涉及一种基于原子范数最小化的三维信道参数估计方法。
背景技术
毫米波通信是5G相关技术中的研究热点，毫米波通信工作在30-300GHz波段，拥有丰富的带宽资源和高速的数据传输速率。同时毫米波通信也带来了诸多挑战，首先高频段的毫米波通信在大气中的路径损耗极大，其次超高带宽的毫米波通信也造成了频率选择性衰落。为了解决这个问题，可以将毫米波与大规模MIMO(multiple inputmultiple output,MIMO)技术结合起来，使用大规模天线阵列的波束成形增益补偿路径损耗。而天线的尺寸与载波波长呈正比，所以毫米波的短波长利于发送端和接收端布置大规模天线阵列。为了克服频率选择性衰落，可以采用OFDM(orthogonal frequency selective fading,OFDM)技术将宽带信道分成多个并行的平坦子信道。所以毫米波大规模MIMO-OFDM***是5G无线通信的研究热点。
波束成形技术可以为毫米波大规模MIMO-OFDM***带来巨大的定向增益，然而波束成形技术需要根据信道状态信息来设计发射波束成形矢量和接收合并矢量。因此，如何准确地估计毫米波大规模MIMO-OFDM***中的信道状态信息成为值得研究的新挑战。
现有技术中，基于张量分解的方法用于毫米波大规模MIMO-OFDM***的信道估计，该方法需要在已知信道阶数的前提下才能进行信道估计，而信道阶数需要用其它算法获得，增加了计算量。并且，目前有二维ANM方法应用于窄带的毫米波MIMO-OFDM***，实际上毫米波带宽很宽，其中时延是一个不能忽略的问题。
综上，毫米波大规模MIMO-OFDM***的信道估计意义重大，值得研究。
发明内容
针对上述技术问题，本发明提供一种基于原子范数最小化的三维信道参数估计方法，运用可分解的三维ANM方法来解决三维的信道估计和信道参数估计问题，能够高精度地得到信道参数并降低计算复杂度。
其中，原子范数最小化(atomic norm minimization,ANM)是一种优化方法，其基本思想是将要估计的信号用已知原子集上的几个原子的简单线性组合表示，并且利用原子的结构提取估计信号中的角度信息。毫米波信道具有多路径稀疏特性，毫米波大规模MIMO-OFDM***的信道估计可被认为三维(出发角度，到达角度，时延)的线性频谱估计问题，该问题可以利用ANM方法恢复出稀疏信号——信道状态信息。
本发明是通过以下技术方案实现的：
一种基于原子范数最小化的三维信道参数估计方法，所述方法用于毫米波大规模MIMO-OFDM***中基站与单个移动台之间的信道参数估计；所述毫米波大规模MIMO-OFDM***使用子载波来发送数据流；
子载波对应的频率域的信道矩阵包括三个具有范德蒙结构的信道参数：出发角度AoA、到达角度AoD、和时间延迟；根据毫米波信道的稀疏散射特性，将所述毫米波大规模MIMO-OFDM***的信道参数估计看作三维的线性谱估计，利用原子范数最小化ANM的方法以无网格的方式高精度恢复出所述信道矩阵和信道矩阵的参数；其中，三维的ANM信道估计包括采用二维ANM估计出发角度、到达角度和信道矩阵；采用一维ANM的估计时间延迟和复增益，最后根据估计得到的出发角度、到达角度、时间延迟和复增益，获得完整的信道矩阵。
进一步地，所述毫米波大规模MIMO-OFDM***包括一个基站和多个移动台，所述基站和所述多个移动台的天线均采用均匀线性阵列；
所述基站配置有NBS根发射天线和
Figure BDA0002084546460000031
条射频链路，每个移动台配置有NMS根天线和
Figure BDA0002084546460000032
条RF链路；
所述基站和所述多个移动台采用混合模数预编码器。
进一步地，采用二维ANM估计出发角度和到达角度，具体为：
第k个子载波对应的频率域的信道矩阵为：
Figure BDA0002084546460000033
其中，L为基站和移动台之间的信道所拥有的散射路径的条数；βl,k为信道复增益，
Figure BDA0002084546460000034
αl是第l条路径上对应的复增益，l≤l≤L；j为虚数单位，τl为第l条路径上的时间延迟参数，fs是采样速率；
Figure BDA0002084546460000035
K是子载波的总数，
Figure BDA0002084546460000036
表示用于发送训练序列的子载波数目；φl∈[0,2π]和θl∈[0,2π]分别是第l条路径上的出发角度和到达角度，aMS(φl)和aBS(θl)分别对应移动台和基站的天线阵列响应矢量，(·)H表示对矩阵或向量进行共轭转置变化；
Figure BDA0002084546460000041
对表示aMS(φl)向量进行共轭转置后的向量；其中，因为τl和出发角度及到达角度是相互独立的，在出发角度及到达角度的估计中将τl视为常数；
建立矩阵形式原子集，将信道矩阵H[k]表示为原子集中几个原子线性组合的形式；给定矩阵形式的原子集
Figure BDA0002084546460000042
为：
Figure BDA0002084546460000043
其中，
Figure BDA0002084546460000044
表示矩阵形式的原子集中的元素；
Figure BDA0002084546460000045
Figure BDA0002084546460000046
基于矩阵形式原子集
Figure BDA0002084546460000047
的信道矩阵H[k]的原子范数定义为：
Figure BDA0002084546460000048
Inf{}表示H[k]的原子范数被最小化；
Figure BDA0002084546460000049
是原子集
Figure BDA00020845464600000410
中的元素，A()表示一种矩阵形式，括号内为矩阵的变量；
当矩阵形式原子集
Figure BDA00020845464600000411
中的原子接近真实的出发角度和到达角度时，信道矩阵H[k]的原子范数
Figure BDA00020845464600000412
将达到最小值；
假设信道矩阵H[k]中的φl与θl相互独立，且φl/θl对应的正弦值sin(φl)/sin(θl)满足充足的频率可分离条件，即：
Figure BDA0002084546460000051
Figure BDA0002084546460000052
式子中，p，q均表示序号；NBS为基站中配置的发射天线数量，NMS为每个移动台配置的天线数量；
满足充足的频率可分离条件的
Figure BDA0002084546460000053
的分解是唯一的和稀疏的，即
Figure BDA0002084546460000054
通过对信道矩阵H[k]的原子范数
Figure BDA0002084546460000055
进行求解，获得组成H[k]的原子，即确定φl/θl。
进一步地，在基站处信道矩阵H[k]是未知的，基站根据接收到的导频信号Y[k]和已知的移动台发送的导频信号s[k]恢复出信道矩阵H[k]和φl/θl，具体为：
根据给定的导频信号Y[k]，将信道矩阵H[k]的原子范数
Figure BDA0002084546460000056
重新写为去除噪声的ANM形式：
Figure BDA0002084546460000057
其中，
Figure BDA0002084546460000058
是噪声控制项，W是基站接收端BS的接收组合矩阵，F是移动台发送端MS的预编码矩阵，S是传输的导频序列；||·||F表示矩阵的范数，μ是与噪声功率相关的归一化参数，
Figure BDA0002084546460000059
其中N＝NBS+NMS；
为了避免无限规划问题，利用半定规划SDP求解，得到：
Figure BDA0002084546460000061
约束条件：
Figure BDA0002084546460000062
其中，uφ,k，uθ,k分别表示T(uφ,k)和T(uθ,k)托普利兹矩阵的第一行；求解得到T(uφ,k)和T(uθ,k)矩阵的第一行；tr(·)代表矩阵的迹，
Figure BDA0002084546460000063
T(uφ,k)是由它的第一行uφ,k确定的托普利兹矩阵；T(uθ,k)是由它的第一行uθ,k确定的托普利兹矩阵；
将T(uφ,k)和T(uθ,k)表示为范德蒙矩阵，即由天线阵列响应向量aMS(φl)和aBS(θl)构成的形式：
Figure BDA0002084546460000064
Figure BDA0002084546460000065
通过对T(uφ,k)和T(uθ,k)做一级范德蒙分解，获得φl/θl。
进一步地，采用一维ANM的估计时间延迟和复增益，具体为：
通过二维ANM信道估计获得φl与θl的估计值，进而获得移动台和基站的天线阵列响应
Figure BDA0002084546460000066
和
Figure BDA0002084546460000067
为了获得第l条路径上的时间延迟参数τl和复增益αl的估计值，将第k个子载波处的毫米波信道改写为OFDM信道形式：
Figure BDA0002084546460000068
其中，OFDM是正交频分复用技术；OFDM中包含了多个载波，h[k]是第k个载波处的信道；h[k]由公式
Figure BDA0002084546460000071
得到，[·]*表示矩阵或向量的伪逆；
给定OFDM信道向量
Figure BDA0002084546460000072
h中的元素与
Figure BDA0002084546460000073
中元素具有类似
Figure BDA0002084546460000074
Figure BDA0002084546460000075
映射关系；故，OFDM信道向量h表示为：
Figure BDA0002084546460000076
其中，时延响应向量aτ(fl)定义为：
Figure BDA0002084546460000077
其中，频率
Figure BDA0002084546460000078
假设τl满足
Figure BDA0002084546460000079
其中Lcp是循环前缀的长度；频率fl的取值范围是
Figure BDA00020845464600000710
的区间；
利用一维的ANM方法来估计信道状态信息时，将信道向量中包含的频率fl的取值范围应从
Figure BDA00020845464600000711
区间变成[0,1]的区间：
假设导频形式是梳状导频，
Figure BDA00020845464600000712
个导频子载波的位置是从等
Figure BDA00020845464600000713
间隔分布的位置集
Figure BDA00020845464600000714
中随机选择，用k′表示
Figure BDA00020845464600000715
个梳状导频子载波的位置索引，子载波位置索引与子载波序号满足
Figure BDA00020845464600000716
的关系，则第k′个子载波处的OFDM信道表示为：
Figure BDA00020845464600000717
此时，OFDM信道向量h中的频率
Figure BDA00020845464600000718
能够采用ANM方法估计时延；
基于原子集
Figure BDA0002084546460000081
的OFDM信道向量h的一维原子范数定义为：
Figure BDA0002084546460000082
其中，
Figure BDA0002084546460000083
假设频率fl满足充足的频率可分离条件，原子范数
Figure BDA0002084546460000084
表示为SDP规划：
Figure BDA0002084546460000085
Figure BDA0002084546460000086
其中，t表示一个常数；u表示托普利兹矩阵T(u)的第一行；hH表示OFDM信道向量h的共轭转置；
Figure BDA0002084546460000087
是由其第一行u构成的PSD托普利兹矩阵；T(u)能够表示为由范德蒙矩阵Aτ(f)组成的形式：
Figure BDA0002084546460000088
其中，Aτ(f)＝[aτ(f1),…,aτ(fL)]，D是半正定的对角矩阵；
采用MUSIC方法从T(u)中提取频率
Figure BDA0002084546460000089
当获得了L条传播路径的时间延迟参数
Figure BDA00020845464600000810
后，构成时延响应向量
Figure BDA00020845464600000811
复数路径增益对应的对角矩阵
Figure BDA00020845464600000812
的计算公式为：
Figure BDA00020845464600000813
对角矩阵
Figure BDA00020845464600000814
完整的信道矩阵
Figure BDA00020845464600000815
能够根据估计值出发角度、到达角度
Figure BDA00020845464600000816
时延
Figure BDA00020845464600000817
和路径增益
Figure BDA00020845464600000818
恢复出来。
本发明的有益技术效果：
(1)本发明提供的三维的基于ANM的信道状态信息估计方法，首次利用三维的ANM方法估计得到毫米波大规模MIMO-OFDM***的信道状态信息。ANM方法是无网格的信道估计方法，可以高精度地估计出多个信道的各项参数值，且不需要预知信道阶数。
(2)基于参数之间的独立性，本发明所述方法将三维的ANM信道估计分解为二维ANM和一维ANM的信道估计，在保持与三维ANM相当的精确度的同时明显降低了计算复杂度。
(3)本发明所述方法通过SDP求解ANM，获得由估计参数的阵列响应向量构成的托普利兹矩阵；并可以采用MUSIC算法来提取托普利兹矩阵中的信道参数，利用峰值搜索确定估计参数的取值。
附图说明
图1不同SNR下的不同参数的MSE；
图2不同SNR下的MSE性能比较；
图3本发明的应用场景模型。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细描述。应当理解，此处所描述的具体实施例仅仅用于解释本发明，并不用于限定本发明。
相反，本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步，为了使公众对本发明有更好的了解，在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
本发明实施例提供一种基于原子范数最小化的三维信道参数估计方法。该信道参数估计方法适用于毫米波大规模MIMO-OFDM***。毫米波大规模MIMO-OFDM***中包含一个基站和多个移动台，基站和移动台的天线采用均匀线性阵列(uniform linear array,ULA)，基站配置有NBS根发射天线和
Figure BDA0002084546460000101
条射频(radio frequency,RF)链路，移动台配置有NMS根天线和
Figure BDA0002084546460000102
条RF链路。
所述基站和所述多个移动台采用混合模数预编码器在本实施例中采用混合模数预编码器的原因是：模拟预编码器的性能和复杂度较低，数字预编码器的性能、复杂度和RF链的开销较高，在本实施例中基站和移动台均采用混合模数预编码器，能够提高预编码器的性能并降低复杂度和RF链开销。
其中，毫米波大规模MIMO-OFDM通信***使用K个子载波来发送NS条数据流。由于导频信号之间的正交性，本实施例所提供的方法仅估计基站与单个移动台之间的信道。
根据毫米波信道的有限散射特性，假设基站和移动台之间的信道拥有L条散射路径，在密集城市的非视线环境中，毫米波大规模MIMO-OFDM***对应的宽带频率选择性衰落信道矩阵为H(τ)：
Figure BDA0002084546460000103
其中，L为基站和移动台之间的信道所拥有的散射路径的条数；αl是第l条路径上对应的复增益，l≤l≤L；j为虚数单位；φl∈[0,2π]和θl∈[0,2π]是第l条路径上的出发角度(angle of departure,AoD)和到达角度(angle of arrival,AoA)，
Figure BDA0002084546460000111
和aBS(θl)分别对应移动台和基站的天线阵列响应矢量，另外，H(τ)表示时延域的信道矩阵，δ(·)表示狄拉克函数；τ是变量，τl是时延；τl为第l条路径上的时间延迟参数。由于基站和移动台都使用ULA天线，因此相应的天线阵列响应矢量具体表示为：
Figure BDA0002084546460000112
Figure BDA0002084546460000113
其中，λ和d分别表示载波波长和相邻天线之间的距离。
根据上面的时延域的信道模型，可以给出第k个子载波对应的频率域的信道矩阵：
Figure BDA0002084546460000114
其中，fs是采样速率，K是子载波的总数，
Figure BDA0002084546460000115
表示用于发送训练序列的子载波数目。
假设移动台在第k个子载波发送的连续M帧导频符号为
Figure BDA0002084546460000116
则基站接收到的导频信号
Figure BDA0002084546460000117
表示为：
Figure BDA0002084546460000118
其中，F[k]表示第k个子载波对应的移动台的预编码矩阵；W[k]表示第k个子载波对应的基站的接收组合矩阵；N[k]表示第k个子载波对应的表示高斯白噪声；为了方便信道估计算法的求解，假设不同子载波的预编码矩阵相同，导频符号也相同，所以省略了[k]。
Figure BDA0002084546460000121
分别是移动台的预编码矩阵和基站的接收组合矩阵。这两个预编码矩阵均由模拟RF预编码矩阵FRF/WRF和数字基带预编码矩阵FBB[k]/WBB[k]构成。为了方便信道估计算法的求解，假设不同子载波的预编码矩阵相同，导频符号也相同。
根据接收到的导频信号Y[k]和已知的发送导频信号s[k]进行信道
Figure BDA0002084546460000122
的估计。已知信道矩阵H[k]中包含了三个拥有范德蒙结构的参数(AoA，AoD和时间延迟)，根据毫米波信道的稀疏散射特性，该***的信道估计问题可以被看做三维的线性谱估计问题。这个三维的线性谱估计问题，可以利用ANM的方法以无网格的方式高精度恢复出以上参数。利用参数之间的独立性，将三维的ANM信道估计问题分解为二维ANM和一维ANM的估计问题，以降低三维ANM信道估计的算法复杂度。
本发明的信道估计方法分为两个步骤：二维ANM估计AoA\AoD和信道矩阵，一维ANM估计时间延迟和复增益。
二维ANM估计AoA和AoD：
采用二维ANM的无网格信道估计方法来获得AoDs/AoAs。因为时间延迟参数τl和AoDs/AoAsφl/θl是相互独立的，所以在AoDs/AoAs的估计中，将τl视为常数，设
Figure BDA0002084546460000123
因此(4)式另写为：
Figure BDA0002084546460000131
其中，L为基站和移动台之间的信道所拥有的散射路径的条数；βl,k为信道复增益，
Figure BDA0002084546460000132
αl是第l条路径上对应的复增益，l≤l≤L；j为虚数单位，τl为第l条路径上的时间延迟参数，fs是采样速率；
Figure BDA0002084546460000133
K是子载波的总数，
Figure BDA0002084546460000134
表示用于发送训练序列的子载波数目；φl∈[0,2π]和θl∈[0,2π]分别是第l条路径上的出发角度和到达角度，aMS(φl)和aBS(θl)分别对应移动台和基站的天线阵列响应矢量，(·)H表示对矩阵或向量进行共轭转置变化；
Figure BDA0002084546460000135
对表示aMS(φl)向量进行共轭转置后的向量；其中，因为τl和出发角度及到达角度是相互独立的，在出发角度及到达角度的估计中将τl视为常数；
ANM方法估计信道矩阵H[k]和AoDs/AoAs首先需要建立矩阵形式原子集，原子集是ANM方法的基础，信道矩阵H[k]可以表示为原子集中几个原子线性组合的形式。给定矩阵形式的原子集为：
Figure BDA0002084546460000136
其中，
Figure BDA0002084546460000137
表示矩阵形式的原子集中的元素；
Figure BDA0002084546460000138
Figure BDA0002084546460000139
H[k]具有可由原子集中的几个简单原子的线性组合表示的特性，所以选择矩阵形式原子集
Figure BDA00020845464600001310
中合适的原子，可以确定H[k]中的参数值(角度值)。基于原子集
Figure BDA0002084546460000141
的信道矩阵H[k]的原子范数定义为：
Figure BDA0002084546460000142
Inf{}表示H[k]的原子范数被最小化；
Figure BDA0002084546460000143
是原子集
Figure BDA0002084546460000144
中的元素，A()表示一种矩阵形式，括号内为矩阵的变量；
当矩阵形式原子集
Figure BDA0002084546460000145
中的原子接近真实的出发角度和到达角度时，信道矩阵H[k]的原子范数
Figure BDA0002084546460000146
将达到最小值；
假设信道矩阵H[k]中的φl与θl相互独立，且φl/θl对应的正弦值sin(φl)/sin(θl)满足充足的频率可分离条件，即：
Figure BDA0002084546460000147
Figure BDA0002084546460000148
式子中，p，q均表示序号；NBS为基站中配置的发射天线数量，NMS为每个移动台配置的天线数量；
满足充足的频率可分离条件的
Figure BDA0002084546460000149
的分解是唯一的和稀疏的，即
Figure BDA00020845464600001410
通过对信道矩阵H[k]的原子范数
Figure BDA00020845464600001411
进行求解，获得组成H[k]的原子，即确定φl/θl。
然而，实际上H[k]在基站处是未知的。基站要根据接收到的导频信号Y[k]和已知的发送导频信号s[k]恢复出H[k]和φl/θl。根据给定的Y[k]，将(8)重新写为去除噪声的ANM形式：
Figure BDA00020845464600001412
其中，
Figure BDA00020845464600001413
是噪声控制项，W是基站接收端BS的接收组合矩阵，F是移动台发送端MS的预编码矩阵，S是传输的导频序列；||·||F表示矩阵的范数，μ是与噪声功率相关的归一化参数，
Figure BDA0002084546460000151
其中N＝NBS+NMS；
为了避免无限规划问题，利用半定规划(semidefinite programming,SDP)求解(10)，
Figure BDA0002084546460000152
其中，tr(·)代表矩阵的迹，
Figure BDA0002084546460000153
Figure BDA0002084546460000154
分别是由它们的第一行
Figure BDA0002084546460000155
和
Figure BDA0002084546460000156
定义的一级的半正定(positivesemidefinite,PSD)托普利兹矩阵。PSD矩阵T(uφ,k)和T(uθ,k)也可表示为范德蒙矩阵(天线阵列响应向量aMS(φl)和aBS(θl))构成的形式，
Figure BDA0002084546460000157
Figure BDA0002084546460000158
综上，通过对T(uφ,k)和T(uθ,k)做一级范德蒙分解，可以获得φl/θl。
优选地，从T(uφ,k)和T(uθ,k)矩阵中提取角度的方法很多，例如矩阵树法和基于旋转不变技术的信号参数估计法等等。使用多信号分类(multiple single classification,MUSIC)法以峰值搜索的方式来获得φl/θl。然后，利用配对算法对得到的φl/θl进行配对。
采用一维ANM的估计时间延迟和复增益，具体为：
通过二维ANM信道估计获得φl与θl的估计值，进而获得移动台和基站的天线阵列响应
Figure BDA0002084546460000161
和
Figure BDA0002084546460000162
为了获得第l条路径上的时间延迟参数τl和复增益αl的估计值，将第k个子载波处的毫米波信道改写为OFDM信道形式：
Figure BDA0002084546460000163
其中，OFDM是正交频分复用技术；OFDM中包含了多个载波，h[k]是第k个载波处的信道；h[k]由公式
Figure BDA0002084546460000164
得到，[·]*表示矩阵或向量的伪逆；
给定OFDM信道向量
Figure BDA0002084546460000165
h中的元素与
Figure BDA0002084546460000166
中元素具有类似
Figure BDA0002084546460000167
Figure BDA0002084546460000168
映射关系；故，OFDM信道向量h表示为：
Figure BDA0002084546460000169
其中，时延响应向量aτ(f)定义为：
Figure BDA00020845464600001610
其中，频率
Figure BDA00020845464600001611
假设时延τl满足
Figure BDA00020845464600001612
其中Lcp是循环前缀的长度。因此，频率fl的取值范围是
Figure BDA00020845464600001613
的区间。
利用一维的ANM方法来估计信道状态信息，信道向量中包含的频率f的取值范围应从
Figure BDA0002084546460000171
区间变成[0,1]的区间。假设***的导频形式是梳状导频，
Figure BDA0002084546460000172
个导频子载波的位置是从等
Figure BDA0002084546460000173
间隔分布的位置集
Figure BDA0002084546460000174
中随机选择，用k′表示
Figure BDA0002084546460000175
个梳状导频子载波的位置索引，子载波位置索引与子载波序号满足
Figure BDA0002084546460000176
的关系。因此第k′个子载波处的OFDM信道表示为
Figure BDA0002084546460000177
现在OFDM信道向量h中的频率
Figure BDA0002084546460000178
因此可以采用ANM方法估计时延。基于原子集
Figure BDA0002084546460000179
的OFDM信道向量h的一维原子范数定义为：
Figure BDA00020845464600001710
其中
Figure BDA00020845464600001711
假设频率fl满足充足的频率可分离条件，随后原子范数
Figure BDA00020845464600001712
可以表示为SDP规划：
Figure BDA00020845464600001713
其中，t表示一个常数；u表示托普利兹矩阵T(u)的第一行；hH表示OFDM信道向量h的共轭转置；
Figure BDA00020845464600001714
是由其第一行u构成的PSD托普利兹矩阵；T(u)能够表示为由范德蒙矩阵Aτ(f)组成的形式，
Figure BDA0002084546460000181
其中，Aτ(f)＝[aτ(f1),…,aτ(fL)]，D是半正定的对角矩阵。同二维ANM中类似，采用MUSIC方法从T(u)中提取频率
Figure BDA0002084546460000182
当获得了L条传播路径的时延
Figure BDA0002084546460000183
后，可以构成时延响应向量
Figure BDA0002084546460000184
然后复数路径增益对应的对角矩阵
Figure BDA0002084546460000185
可以通过计算
Figure BDA0002084546460000186
得到，对角矩阵
Figure BDA0002084546460000187
最终，完整的信道矩阵
Figure BDA0002084546460000188
可以根据估计值AoAs,AoDs
Figure BDA0002084546460000189
时延
Figure BDA00020845464600001810
和路径增益
Figure BDA00020845464600001811
恢复出来。
对于毫米波MIMO-OFDM***信道估计，进行模拟仿真来评估本发明提出的可分解的三维的ANM方法的性能。将本发明所述方法同传统的矢量化三维的ANM和压缩感知(compressed sensing,CS)方法作对比。CS方法采用正交匹配追踪算法来解决估计问题，并在40×40×40(AoD-AoA-时间延迟)的网格上求解。均方误差(mean square error,MSE)作为评估所有方法性能的标准。仿真结果均是200次蒙特卡罗试验结果的平均值。如果没有特别说明，用D-ANM表示可分解的三维ANM方法，ANM表示传统的矢量化ANM法。仿真参数设置见表1。
表1***参数设置
Figure BDA00020845464600001812
Figure BDA0002084546460000191
首先比较信道参数(AoD,AoA和时延)估计的准确性。如图1(a)所示，本发明提供的方法可以精确的估计出角度值。D-ANM方法的MSE性能与ANM的MSE性能非常接近。而CS性能较差，且随着SNR的增加，其MSE性能无明显改变。图1(b)是关于时延估计性能的比较，ANM法中时延估计的性能和角度估计的性能相同，但D-ANM的方法对时延的估计性能略低于ANM，这是由于D-ANM中一维ANM的时延估计需要利用之前二维ANM所获得的信道矩阵，信道矩阵估计的精准度影响后续的时延估计。
接下来，比较不同信道估计方法的性能。从图2可以看出，所有方法的MSE随着SNR的增加而减小。其中，D-ANM方法的估计准确度与ANM相当，它们的MSE均高于CS。D-ANM方法的性能与ANM相当，且高于CS方法，D-ANM方法的计算复杂度远远小于ANM。
本发明所提供的方法使用ANM的方法来进行毫米波大规模MIMO-OFDM***的三维的信道估计，优势在于三点：其一，适用于宽带毫米波通信***；其二，获得的参数估计值精度极高；其三，计算复杂度较低。

Claims (3)

1.一种基于原子范数最小化的三维信道参数估计方法，其特征在于，所述方法用于毫米波大规模MIMO-OFDM***中基站与单个移动台之间的信道参数估计；所述毫米波大规模MIMO-OFDM***使用子载波来发送数据流；
子载波对应的频率域的信道矩阵包括三个具有范德蒙结构的信道参数：出发角度AoA、到达角度AoD、和时间延迟；根据毫米波信道的稀疏散射特性，将所述毫米波大规模MIMO-OFDM***的信道参数估计看作三维的线性谱估计，利用原子范数最小化ANM的方法以无网格的方式高精度恢复出所述信道矩阵和信道矩阵的参数；其中，三维的ANM信道估计包括采用二维ANM估计出发角度、到达角度和信道矩阵；采用一维ANM的估计时间延迟和复增益，最后根据估计得到的出发角度、到达角度、时间延迟和复增益，获得完整的信道矩阵；
其中，采用二维ANM估计出发角度和到达角度，具体为：
第k个子载波对应的频率域的信道矩阵为：
Figure FDA0002379143690000011
其中，L为基站和移动台之间的信道所拥有的散射路径的条数；βl，k为信道复增益，
Figure FDA0002379143690000012
αl是第l条路径上对应的复增益，l≤l≤L；j为虚数单位，τl为第l条路径上的时间延迟参数，fs是采样速率；
Figure FDA0002379143690000013
K是子载波的总数，
Figure FDA0002379143690000014
表示用于发送训练序列的子载波数目；φl∈[0，2π]和θl∈[0，2π]分别是第l条路径上的出发角度和到达角度，aMS(φl)和aBS(θl)分别对应移动台和基站的天线阵列响应矢量，(·)H表示对矩阵或向量进行共轭转置变化；
Figure FDA0002379143690000021
对表示aMS(φl)向量进行共轭转置后的向量；其中，因为τl和出发角度及到达角度是相互独立的，在出发角度及到达角度的估计中将τl视为常数；
建立矩阵形式原子集，将信道矩阵H[k]表示为原子集中几个原子线性组合的形式；给定矩阵形式的原子集
Figure FDA0002379143690000022
为：
Figure FDA0002379143690000023
其中，
Figure FDA0002379143690000024
表示矩阵形式的原子集中的元素；
Figure FDA0002379143690000025
Figure FDA0002379143690000026
基于矩阵形式原子集
Figure FDA0002379143690000027
的信道矩阵H[k]的原子范数定义为：
Figure FDA0002379143690000028
Inf{}表示H[k]的原子范数被最小化；
Figure FDA0002379143690000029
是原子集
Figure FDA00023791436900000210
中的元素，A()表示一种矩阵形式，括号内为矩阵的变量；
当矩阵形式原子集
Figure FDA00023791436900000211
中的原子接近真实的出发角度和到达角度时，信道矩阵H[k]的原子范数
Figure FDA00023791436900000212
将达到最小值；
假设信道矩阵H[k]中的φl与θl相互独立，且φl/θl对应的正弦值sin(φl)/sin(θl)满足充足的频率可分离条件，即：
Figure FDA00023791436900000213
Figure FDA00023791436900000214
式子中，p，q均表示序号；NBS为基站中配置的发射天线数量，NMS为每个移动台配置的天线数量；
满足充足的频率可分离条件的
Figure FDA0002379143690000031
的分解是唯一的和稀疏的，即
Figure FDA0002379143690000032
通过对信道矩阵H[k]的原子范数
Figure FDA0002379143690000033
进行求解，获得组成H[k]的原子，即确定φl/θl；
采用一维ANM的估计时间延迟和复增益，具体为：
通过二维ANM信道估计获得φl与θl的估计值，进而获得移动台和基站的天线阵列响应
Figure FDA0002379143690000034
和
Figure FDA0002379143690000035
为了获得第l条路径上的时间延迟参数τl和复增益αl的估计值，将第k个子载波处的毫米波信道改写为OFDM信道形式：
Figure FDA0002379143690000036
其中，OFDM是正交频分复用技术；OFDM中包含了多个载波，h[k]是第k个载波处的信道；h[k]由公式
Figure FDA0002379143690000037
得到，[·]*表示矩阵或向量的伪逆；
给定OFDM信道向量
Figure FDA0002379143690000038
h中的元素与
Figure FDA0002379143690000039
中元素具有类似
Figure FDA00023791436900000310
Figure FDA00023791436900000311
映射关系；故，OFDM信道向量h表示为：
Figure FDA00023791436900000312
其中，时延响应向量aτ(fl)定义为：
Figure FDA0002379143690000041
其中，频率
Figure FDA0002379143690000042
假设τl满足
Figure FDA0002379143690000043
其中Lcp是循环前缀的长度；频率fl的取值范围是
Figure FDA0002379143690000044
的区间；
利用一维的ANM方法来估计信道状态信息时，将信道向量中包含的频率fl的取值范围应从
Figure FDA0002379143690000045
区间变成[0，1]的区间：
假设导频形式是梳状导频，
Figure FDA0002379143690000046
个导频子载波的位置是从等
Figure FDA0002379143690000047
间隔分布的位置集
Figure FDA0002379143690000048
中随机选择，用k′表示
Figure FDA0002379143690000049
个梳状导频子载波的位置索引，子载波位置索引与子载波序号满足
Figure FDA00023791436900000410
的关系，则第k′个子载波处的OFDM信道表示为：
Figure FDA00023791436900000411
此时，OFDM信道向量h中的频率
Figure FDA00023791436900000412
能够采用ANM方法估计时延；
基于原子集
Figure FDA00023791436900000413
的OFDM信道向量h的一维原子范数定义为：
Figure FDA00023791436900000414
其中，
Figure FDA00023791436900000415
假设频率fl满足充足的频率可分离条件，原子范数
Figure FDA00023791436900000416
表示为SDP规划：
Figure FDA0002379143690000051
Figure FDA0002379143690000052
其中，t表示一个常数；u表示托普利兹矩阵T(u)的第一行；hH表示OFDM信道向量h的共轭转置；
Figure FDA0002379143690000053
是由其第一行u构成的PSD托普利兹矩阵；T(u)能够表示为由范德蒙矩阵Aτ(f)组成的形式：
Figure FDA0002379143690000054
其中，Aτ(f)＝[aτ(f1)，…，aτ(fL)]，D是半正定的对角矩阵；
采用MUSIC方法从T(u)中提取频率
Figure FDA0002379143690000055
当获得了L条传播路径的时间延迟参数
Figure FDA0002379143690000056
后，构成时延响应向量
Figure FDA0002379143690000057
复数路径增益对应的对角矩阵
Figure FDA0002379143690000058
的计算公式为：
Figure FDA0002379143690000059
对角矩阵
Figure FDA00023791436900000510
完整的信道矩阵
Figure FDA00023791436900000511
能够根据估计值出发角度、到达角度
Figure FDA00023791436900000512
时延
Figure FDA00023791436900000513
和路径增益
Figure FDA00023791436900000514
恢复出来。
2.根据权利要求1所述一种基于原子范数最小化的三维信道参数估计方法，其特征在于，所述毫米波大规模MIMO-OFDM***包括一个基站和多个移动台，所述基站和所述多个移动台的天线均采用均匀线性阵列；
所述基站配置有NBS根发射天线和
Figure FDA00023791436900000515
条射频链路，每个移动台配置有NMS根天线和
Figure FDA00023791436900000516
条RF链路；
所述基站和所述多个移动台采用混合模数预编码器。
3.根据权利要求1所述一种基于原子范数最小化的三维信道参数估计方法，其特征在于，在基站处信道矩阵H[k]是未知的，基站根据接收到的导频信号Y[k]和已知的移动台发送的导频信号s[k]恢复出信道矩阵H[k]和φl/θl，具体为：
根据给定的导频信号Y[k]，将信道矩阵H[k]的原子范数
Figure FDA0002379143690000061
重新写为去除噪声的ANM形式：
Figure FDA0002379143690000062
其中，
Figure FDA0002379143690000063
是噪声控制项，W是基站接收端BS的接收组合矩阵，F是移动台发送端MS的预编码矩阵，S是传输的导频序列；||·||F表示矩阵的范数，μ是与噪声功率相关的归一化参数，
Figure FDA0002379143690000064
其中N＝NBS+NMS；
为了避免无限规划问题，利用半定规划SDP求解，得到：
Figure FDA0002379143690000065
约束条件：
Figure FDA0002379143690000066
其中，uφ，k，uθ，k分别表示T(uφ，k)和T(uθ，k)托普利兹矩阵的第一行；求解得到T(uφ，k)和T(uθ，k)矩阵的第一行；tr(·)代表矩阵的迹，
Figure FDA0002379143690000067
T(uφ，k)是由它的第一行uφ，k确定的托普利兹矩阵；T(uθ，k)是由它的第一行uθ，k确定的托普利兹矩阵；
将T(uφ，k)和T(uθ，k)表示为范德蒙矩阵，即由天线阵列响应向量aMS(φl)和aBS(θl)构成的形式：
Figure FDA0002379143690000071
Figure FDA0002379143690000072
通过对T(uφ，k)和T(uθ，k)做一级范德蒙分解，获得φl/θl。
CN201910483315.9A 2019-06-04 2019-06-04 一种基于原子范数最小化的三维信道参数估计方法 Active CN110213185B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910483315.9A CN110213185B (zh) 2019-06-04 2019-06-04 一种基于原子范数最小化的三维信道参数估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910483315.9A CN110213185B (zh) 2019-06-04 2019-06-04 一种基于原子范数最小化的三维信道参数估计方法

Publications (2)

Publication Number Publication Date
CN110213185A CN110213185A (zh) 2019-09-06
CN110213185B true CN110213185B (zh) 2020-06-12

Family

ID=67790727

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910483315.9A Active CN110213185B (zh) 2019-06-04 2019-06-04 一种基于原子范数最小化的三维信道参数估计方法

Country Status (1)

Country Link
CN (1) CN110213185B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112737649B (zh) * 2020-12-25 2022-04-01 杭州电子科技大学 一种基于角度网格优化与范数约束的毫米波信道估计方法
CN112887250A (zh) * 2021-01-12 2021-06-01 电子科技大学 一种基于张量分解的毫米波宽带mimo***的信道估计方法
CN113376605B (zh) * 2021-05-20 2022-07-26 南京航空航天大学 一种基于原子范数最小化的mmwave-MIMO雷达超分辨DOA估计方法
CN114286307B (zh) * 2022-01-07 2024-03-26 深圳泓越信息科技有限公司 一种基于矩阵束的信道状态信息参数估计方法
CN115022129B (zh) * 2022-03-08 2024-03-22 东南大学 基于anm的多用户上行传输ris辅助***的信道估计方法
CN114928518B (zh) * 2022-05-18 2024-01-30 东南大学 毫米波mimo-ofdm***中基于3d-music算法的信道估计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106031068A (zh) * 2014-04-02 2016-10-12 华为技术有限公司 一种基于波束成形的通信方法及装置
CN107817465A (zh) * 2017-10-12 2018-03-20 中国人民解放军陆军工程大学 超高斯噪声背景下的基于无网格压缩感知的doa估计方法
CN108768585A (zh) * 2018-04-27 2018-11-06 南京邮电大学 基于深度学习的上行免信令noma***多用户检测方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10263763B2 (en) * 2015-02-26 2019-04-16 New York University Systems, methods, and computer-accessible media for measuring or modeling a wideband, millimeter-wave channel and methods and systems for calibrating same
CN107181705B (zh) * 2017-05-24 2020-04-14 深圳先进技术研究院 一种半盲信道估计方法及***
CN107171985B (zh) * 2017-06-07 2020-01-10 深圳先进技术研究院 一种信道估计方法及***

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106031068A (zh) * 2014-04-02 2016-10-12 华为技术有限公司 一种基于波束成形的通信方法及装置
CN107817465A (zh) * 2017-10-12 2018-03-20 中国人民解放军陆军工程大学 超高斯噪声背景下的基于无网格压缩感知的doa估计方法
CN108768585A (zh) * 2018-04-27 2018-11-06 南京邮电大学 基于深度学习的上行免信令noma***多用户检测方法

Also Published As

Publication number Publication date
CN110213185A (zh) 2019-09-06

Similar Documents

Publication Publication Date Title
CN110213185B (zh) 一种基于原子范数最小化的三维信道参数估计方法
CN108933745B (zh) 一种基于超分辨率角度和时延估计的宽带信道估计方法
CN111106859B (zh) 毫米波/太赫兹网络大规模mimo无线传输方法
Lin et al. A new view of multi-user hybrid massive MIMO: Non-orthogonal angle division multiple access
CN106559367B (zh) 基于低秩张量分解的mimo-ofdm***毫米波信道估计方法
CN103152293B (zh) Ofdm***中信道估计的方法和设备
CN108886826A (zh) 用于无线多天线和频分双工***的混合波束成形方法
CN107483091B (zh) 一种fdd大规模mimo-ofdm***下的信道信息反馈算法
CN107465636B (zh) 一种毫米波大规模阵列空频双宽带***的信道估计方法
CN109412983B (zh) 一种基于dft域的无网格化大规模mimo信道估计算法
CN110099016A (zh) 一种基于深度学习网络的毫米波稀疏阵面信道估计方法
CN105915473B (zh) 一种基于压缩感知技术的ofdm***参数化信道估计及均衡方法
CN104052691A (zh) 基于压缩感知的mimo-ofdm***信道估计方法
WO2022104993A1 (zh) 天波大规模mimo通信方法、模型及***
CN108599825A (zh) 一种基于mimo-ofdm毫米波结构的混合编码方法
CN108881076A (zh) 一种基于压缩感知的mimo-fbmc/oqam***信道估计方法
CN108881074A (zh) 一种低精度混合架构下宽带毫米波信道估计方法
CN110034916B (zh) 基于终端反馈的天线相位同步与信道互易性校准的方法
Shafin et al. DoA estimation and RMSE characterization for 3D massive-MIMO/FD-MIMO OFDM system
CN113411107A (zh) 一种基于波束束的毫米波大规模mimo***宽带信道估计方法
WO2023185373A1 (zh) 天波大规模mimo三重波束基信道建模及信道信息获取
CN107276726B (zh) 一种Massive MIMO FBMC波束空时编码下行传输方法
CN101155164B (zh) 一种dft扩频的广义多载波***的sinr估计方法
CN112769462B (zh) 一种基于联合参数学习的毫米波mimo宽带信道估计方法
CN114143896A (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240115

Address after: 200439, 9th Floor, Building 1, No. 2816 Yixian Road, Baoshan District, Shanghai

Patentee after: Shanghai Zhuangyan Automation Technology Co.,Ltd.

Address before: 100044 Beijing city Haidian District Xizhimen Shangyuan Village No. 3

Patentee before: Beijing Jiaotong University