大规模MIMO信道估计的信息几何方法及系统

大规模mimo信道估计的信息几何方法及系统
技术领域
1.本发明属于通信技术领域，涉及大规模mimo信道估计的信息几何方法及相关系统。


背景技术：

2.大规模多输入多输出(mimo，multiple-input multiple-ouput)技术是5g蜂窝系统及其持续演进系统的核心使能技术。在大规模mimo系统中，配备大量天线的基站(bs，base station)可以在相同的时间和频率资源上同时为数十个用户提供服务，从而潜在地提供了巨大的容量增益，并显著提高能源效率。在所有类型的天线阵列中，均匀平面阵列(upa，uniform plane array)具有紧凑的尺寸和三维覆盖能力，是实际应用的良好选择。正交频分复用(ofdm，orthogonal frequency division multiplexing)技术是一种多载波调制技术，可以减轻宽带无线通信中频率选择性衰落带来的影响。大规模mimo-ofdm在5g系统中发挥着重要作用，并在未来的6g系统中受到广泛关注。
3.在大规模mimo-ofdm系统中，信道估计起着至关重要的作用，这是因为系统性能高度依赖于系统所获取的信道信息的质量。在实际系统中，辅助导频的信道估计，即发射机周期性地发送导频信号，接收机根据接收到的导频信号获得信道状态信息(csi，channel state information)，是常用的信道估计方法。给定接收到的导频信号，信道估计的任务是获得信道参数的后验统计信息。当信道参数的先验分布为高斯分布时，其后验分布亦为高斯分布，其后验信息由后验均值和后验协方差矩阵给出。然而，由于大规模mimo-ofdm系统中信道的维度较大，后验信息的计算具有挑战性。传统信道估计算法，例如最小均方误差(mmse，minimum mean squared error)估计，由于计算过程中存在大维矩阵求逆，因此很难被应用在大规模mimo-ofdm系统中。随着基站侧天线单元个数及可支持的用户终端个数的进一步大幅增加，常规信道信息获取的计算复杂性和导频开销都将大幅提升，成为需要破解的瓶颈问题。
4.由后验概率密度函数定义的空间可以视作一具有黎曼结构的可微流形。因此，微分几何学中的定义和工具可以较好地应用于此，这正是信息几何的主题之一。因此，将信息几何应用于信道估计具备合理性。信息几何的主要思想是通过将概率密度函数(pdf，probability density function)的参数空间视为可微流形来研究特定pdf集合的内在几何结构。近年来，信息几何已经成功应用于多传感器估计融合、误报率检测和广义贝叶斯预测等方面。使用信息几何理论分析估计问题有以下几个优点：一是信息几何可以为现有算法的理论分析提供一个统一的框架；二是它从几何视角提供了对统计模型的直观理解，可以促进对现有问题的内在研究；此外，信息几何还可以从更一般和内在的角度改进算法。


技术实现要素：

5.发明目的：针对现有技术的不足，本发明公开了大规模mimo信道估计信息几何方法及相关系统，能够获得各用户终端信道的后验信息，在保证估计性能的同时，相比现有类
似技术手段可以进一步降低计算复杂度。
6.技术方案：为了达到上述目的，本发明提供如下技术方案：
7.大规模mimo信道估计方法，包括如下步骤：
8.基站侧/用户终端通过上行/下行探测获得信道的先验统计信息；
9.基站侧/用户终端通过接收到的上行/下行导频信号以及先验统计信息，利用信息几何方法获取信道的后验统计信息；所述信息几何方法将高斯分布的集合定义为原始流形，并且根据信道的后验分布构造目标流形和辅助流形，所述目标流形和辅助流形均为原始流形的子流形，迭代计算辅助流形中的分布在目标流形上的m-投影，并根据m-投影更新目标流形以及辅助流形中的分布，最后以目标流形上分布的均值和方差作为信道估计的后验均值和后验方差。
10.作为优选，所述目标流形是一类各个元素相互独立的高斯分布的集合，辅助流形是一类协方差矩阵为一对角阵与秩为1矩阵之和的逆矩阵的高斯分布的集合；m-投影通过最小化辅助流形中的分布与目标流形之间的kl散度得到。
11.作为优选，所述目标流形和辅助流形中的分布的均值及协方差矩阵通过各自辅助计算的参数矢量和参数实对角矩阵表示；其中目标流形中分布的协方差表示为先验方差的逆与参数实对角矩阵之差的逆矩阵，均值通过协方差矩阵与参数矢量的乘积表示；辅助流形中的分布的协方差矩阵表示为先验方差的逆与参数实对角矩阵之差再与秩为1矩阵之和的逆矩阵，其中秩为1矩阵由感知矩阵中相应行与噪声方差表示，均值通过协方差与参数矢量结合感知矩阵中相应行、接收导频信号矢量相应元素以及噪声方差组成矢量的乘积表示。
12.作为优选，将辅助流形的协方差矩阵使用sherman-morrison公式展开。
13.作为优选，利用信息几何方法获取信道估计后验均值和后验方差的步骤包括：
14.(1)建立大规模mimo信道的原始流形、目标流形以及辅助流形；
15.(2)初始化辅助流形以及目标流形上分布的参数；
16.(3)根据辅助流形上分布的参数、接收到的导频信号以及先验信道统计信息计算辅助流形中的分布在目标流形上的m-投影；
17.(4)根据m-投影更新辅助流形以及目标流形上分布的参数；重复步骤(3)-(4)直至预设迭代次数或目标流形上分布的参数收敛。
18.作为优选，基于空频波束基(space-frequency beam based)信道统计表征模型计算空频波束域信道的先验统计信息和后验统计信息，空频波束基信道统计表征模型中，空间频率域信道矩阵由空频波束域信道矩阵左乘采样空间舵矢量矩阵并右乘采样频率舵矢量矩阵的转置矩阵后得到，空频波束域信道各元素是统计独立的；对于基站侧，利用采样空间舵矢量矩阵和采样频率舵矢量矩阵将空频波束域信道的后验均值和后验方差转换为空间频率域信道的后验均值和后验方差；对于用户终端侧，将各自的空频波束域信道的后验统计信息并反馈给基站，基站侧利用采样空间舵矢量矩阵和采样频率舵矢量矩阵将所获得的空频波束域信道的后验均值和后验方差转换为空间频率域信道的后验均值和后验方差。
19.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的大规模mimo信道估计方法。
20.一种大规模mimo通信系统，包括基站和多个用户终端，所述基站用于：通过上行探
测获得各个用户终端的信道的先验统计信息；通过接收到的上行导频信号以及先验统计信息，利用信息几何方法获取各个用户终端的后验统计信息；所述信息几何方法将高斯分布的集合定义为原始流形，并且根据信道的后验分布构造目标流形和辅助流形，所述目标流形和辅助流形均为原始流形的子流形，迭代计算辅助流形中的分布在目标流形上的m-投影，并根据m-投影更新目标流形以及辅助流形中的分布，最后以目标流形上分布的均值和方差为信道估计后验均值和后验方差。
21.一种大规模mimo通信系统，包括基站和多个用户终端，所述用户终端用于：通过下行信道探测获得各自信道的先验统计信息；通过接收到的下行导频信号以及先验统计信息，利用信息几何方法以及信道预测方法获取各自信道的后验统计信息并反馈给基站；所述信息几何方法将高斯分布的集合定义为原始流形，并且根据信道的后验分布构造目标流形和辅助流形，所述目标流形和辅助流形均为原始流形的子流形，迭代计算辅助流形中的分布在目标流形上的m-投影，并根据m-投影更新目标流形以及辅助流形中的分布，最后以目标流形上分布的均值和方差为信道估计后验均值和后验方差。
22.一种大规模mimo通信系统，包括基站和多个用户终端，所述基站或用户终端包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的大规模mimo信道估计方法。
23.有益效果：与现有技术相比，本发明提出的大规模mimo信道估计的信息几何方法能够在保证信道估计准确度的前提下，以较低的计算复杂度和导频开销获得信道的后验均值和后验方差。所获得的后验均值和后验方差可以进一步应用于鲁棒预编码以及鲁棒检测中，提升系统性能，从而进一步提升系统的整体传输效率。
附图说明
24.图1为本发明实施例的信道估计方法流程图；
25.图2为本发明另一实施例的信道估计方法流程图；
26.图3为本发明实施例中信道估计的信息几何方法流程图；
27.图4为本发明实施例中信息几何方法与现有方法的信道估计性能比较示意图；
28.图5为本发明实施例中信噪比20db时信息几何方法与现有方法的收敛曲线比较示意图。
具体实施方式
29.以下将结合具体实施例对本发明提供的技术方案进行详细说明，应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
30.如图1所示，本发明实施例公开的一种大规模mimo信道估计方法，适用于基站侧，包括基站侧通过上行探测获得各个用户终端的信道的先验统计信息；基站侧通过接收到的上行导频信号以及先验统计信息，利用信息几何方法获取各个用户终端的信道的后验统计信息，包括后验均值和后验方差；其中信息几何方法将高斯分布的集合定义为原始流形，并且根据信道的后验分布构造目标流形和辅助流形，迭代计算辅助流形中的分布在目标流形上的m-投影，并根据m-投影更新目标流形以及辅助流形中的分布，最后以目标流形上分布的均值和方差为信道估计的后验均值和后验方差。
31.如图2所示，本发明另一实施例公开的一种大规模mimo信道估计方法，适用于用户终端侧，包括用户终端通过下行信道探测获得各自信道的先验统计信息；用户终端通过接收到的下行导频信号以及先验统计信息，利用信息几何方法获取各自信道的后验统计信息，后验信息包括后验均值和后验方差；其中信息几何方法与基站侧方法一致，只是将多用户退化为单用户。用户终端可以是手机、车载设备、智能装备等移动终端或固定终端。
32.图3示意了利用信息几何方法获取后验统计信息的具体步骤，包括：(1)建立大规模mimo信道的原始流形、辅助流形以及目标流形；(2)初始化辅助流形以及目标流形上分布的参数；(3)根据辅助流形上分布的参数、接收到的导频信号以及先验信道统计信息计算辅助流形中的分布在目标流形上的m-投影；(4)根据m-投影更新辅助流形以及目标流形上分布的参数；重复步骤(3)-(4)直至预设迭代次数或目标流形上分布的参数收敛，目标流形上分布的均值和方差即分别为信道的后验均值和后验方差。
33.本发明方法主要适用于基站侧配备大规模天线阵列以同时服务多个用户的大规模mimo系统。下面结合具体的通信系统实例对本发明涉及的信道估计信息几何方法的具体实现过程作详细说明，需要说明的是本发明方法不仅适用于下面示例所举的具体系统模型，也同样适用于其它配置的系统模型。
34.一、系统配置
35.考虑工作在时分复用模式(tdd，time division multiplexing)的大规模mimo-ofdm系统。基站侧配备upa天线阵列，其中天线数为nr＝n
r,v
×nr,h
，n
r,v
和n
r,h
分别为每列和每行的天线数，水平方向和垂直方向的天线间距分别记为δv和δh。基站同时为同一个小区内k个配备单天线的用户服务。在ofdm调制中，子载波个数为nc，系统采样间隔和循环前缀长度分别记记为ts和ng。记子载波集为其中用于上行训练额的子载波集记为其大小为在tdd模式下，由于信道互易性，上行训练获取的csi可以用于上行信号检测和下行预编码传输，因此实施例中考虑上行大规模mimo-ofdm信道估计。
36.二、空频波束基信道统计表征模型和信道估计问题陈述
37.下面对系统模型、空频波束基信道统计表征模型以及信道估计问题进行详细阐述。
38.1.系统模型
39.记为第k个用户在单个ofdm符号内的频域发送序列。在上行训练中，第n个子载波上的频域接收信号矢量可以表示为
[0040][0041]
其中为第k个用户在第n个子载波上的空域信道，为循环对称高斯噪声，噪声功率为进一步，定义第k个用户的空间频率域信道矩阵为
[0042][0043]
令其中xk＝[xk[n1]
…
xk[n2]]
t
，以及令
上标t表示矩阵或向量的转置。进一步，空间频率域接收信号模型可以写为
[0044][0045]
2.空频波束基信道统计表征模型
[0046]
定义方向余弦u＝sinθ,v＝cosθsinφ，其中θ,φ∈[-π/2,π/2]分别为极角和方位角。进一步，定义空间舵矢量v(u,v)为
[0047][0048][0049][0050]
其中表示kronecker积，λc为波长。定义频率舵矢量u(τ)为
[0051][0052]
其中τ为时延，δf＝1/ncts为子载波间隔。进一步，对方向余弦u,v∈[-1,1]和时延τ∈[0,ngts)采样量化为其中并且有nv,nh,n
τ
为采样倍数。则采样空间舵矢量矩阵和采样频率舵矢量矩阵可以分别表示为
[0053][0054][0055][0056][0057]
则所考虑空频波束基信道统计表征模型为：空间频率域信道矩阵由空频波束域信道左乘采样空间舵矢量矩阵并右乘采样频率舵矢量矩阵的转置矩阵后得到，其中空频波束域信道各元素是统计独立的。具体表达式为
[0058]gk
＝vh
kft
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12)
[0059]
其中定义为空频波束域信道矩阵。假设采样倍数nv≥n
r,v
,nh≥n
r,h
,n
τ
≥nf，并定义精细化因子当精细化因子为整数时，可以验证采样舵矢量矩阵有dft结构。
[0060]
本领域技术人员可以理解的是，上述模型中的具体矢量表示仅以均匀平面阵列为例，对于采用如均匀线阵、均匀圆阵等不同天线阵列的系统，上述空频波束基信道统计表征
模型依然适用，仅需将v改为对应的空间采样舵矢量矩阵即可。
[0061]
3.问题陈述
[0062]
结合空间频率域接收信号模型(3)和空频波束基信道统计表征模型(12)，可以得到
[0063]
y＝vham z
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0064]
其中进一步，对等式两边矢量化并去除空频波束域中的零元可以得到
[0065]
y＝ah z
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0066]
其中为根据ha非零元位置从抽取相关列得到的矩阵,n＝nrn
p
,m为ha中的非零元个数；y,z分别为对y,z矢量化后的矢量；h为ha矢量化并抽取非零元得到的矢量。信道估计问题为利用接收导频信号y，获取各用户终端的空频波束域信道h的后验统计信息。空频波束域信道的后验统计信道信息为空频波束域信道的后验均值和后验方差；信道后验均值和后验方差包括：基站侧在给定接收到的上行导频信号条件下的空频波束域信道的条件均值和条件方差。假设其中d为正定实对角矩阵，为各用户的先验方差，i为单位阵，为噪声方差，表示均值为μ、协方差矩阵为σ的循环对称高斯分布。于是，空频波束域信道后验分布p(h|y)也属于高斯分布，其均值和协方差分别为
[0067][0068][0069]
其中上标h为矩阵或向量的共轭转置。可以验证后验均值(15)和最小均方误差(mmse)估计结果是等价的。空频波束域信道的后验信息，包括后验均值(15)和后验协方差(16)的计算复杂度为以下将基于信息几何方法，设计低复杂度的空频波束域信道估计方法。
[0070]
本领域技术人员可以理解的是，接收信号模型(14)中的感知矩阵a在其它信道统计表征模型中可以有不同的形式，任意满足接收信号模型(14)以及相同高斯先验的统计推断问题，均可以用信息几何方法求解，而不仅仅局限于波束域信道估计问题。接收信号模型(14)中去除零元是为了减少计算量，是优选方案而非必要。
[0071]
三、空频波束域信道估计的信息几何方法
[0072]
1.原始流形、目标流形和辅助流形的建立
[0073]
首先将高斯分布的集合定义为原始流形，具体为
[0074]mor
:p(h)＝pg(h；μ,σ)＝exp{-(h-μ)σ-1
(h-μ)h}
ꢀꢀꢀꢀꢀꢀꢀ
(17)
[0075]
其中μ,σ分别为原始流形上分布pg(h；μ,σ)的均值和协方差。进一步，在高斯假设的基础上，后验分布p(h|y)可以表示为
[0076][0077]
其中hi为h第i个元素，yn为y第n个元素，为a的第n行；dh＝f(0,-d-1
)，th＝f(h,i
⊙
(hh
t
))，其中
⊙
为hadamard积，函数f(a,a)＝[a
t
,vec(a)
t
]
t
；符号表示同维度矢量的操作符号，c,ψq为归一化因子；cn(h)为
[0078][0079]
进一步，定义目标流形为一类各个元素相互独立的高斯分布的集合，具体为
[0080][0081]
其中分别为辅以计算的参数矢量和参数实对角矩阵，为m维实对角矩阵的集合。上述目标流形为原始流形的子流形，目标流形上分布的均值μ0及方差σ0和θ0,θ0的关系为
[0082][0083]
σ0＝(d-1-θ0)-1
ꢀꢀꢀꢀꢀꢀꢀ
(22)
[0084]
进一步，定义n个辅助流形，辅助流形是一类协方差矩阵为一对角阵与秩为1矩阵之和的逆矩阵的高斯分布的集合，其中第n个辅助流形包含cn(h)，具体为
[0085][0086]
其中上述辅助流形为原始流形的子流形，辅助流形上分布的均值μn及协方差σn和θn,θn的关系为
[0087][0088][0089]
由于辅助流形上分布的协方差矩阵σn为一对角阵与秩为1的矩阵之和的逆矩阵，可以使用sherman-morrison公式展开，因此σn可以写为
[0090]
[0091]
2.信道估计方法
[0092]
信道估计信息几何方法分为四个步骤，分别为建立大规模mimo空频波束域信道的原始流形、辅助流形以及目标流形，初始化辅助流形以及目标流形上分布的参数，根据辅助流形上分布的参数、接收到的导频信号以及空频波束域先验信道信息计算辅助流形中的分布在目标流形上的m-投影，以及根据m-投影更新辅助流形以及目标流形上分布的参数。其中建立大规模mimo空频波束域信道的原始流形、辅助流形以及目标流形这一步骤见上一节，下面详细介绍后续步骤。
[0093]
(1)初始化辅助流形以及目标流形上分布的参数
[0094]
设置迭代次数初始值t＝0，初始化目标流形参数以及辅助流形参数
[0095]
(2)计算辅助流形中的分布在目标流形上的m-投影
[0096]
根据辅助流形上分布的参数、接收到的导频信号以及空频波束域先验信道信息计算辅助流形中的分布在目标流形上的m-投影。将辅助流形上的分布向目标流形m-投影，得到的分布记为所述m-投影为在目标流形上寻找一个分布使和辅助流形上分布之间的kl散度最小，即
[0097][0098]
其中kl散度定义为
[0099][0100]
其中表示关于概率分布p(x)的期望。给定辅助流形参数kl散度可以表示为其中c
p
为常数。可以表示为
[0101][0102]
进一步分别对θ
on
求偏导
[0103]
[0104][0105]
其中上横线表示取共轭。令偏导等于0可以得到
[0106][0107][0108]
由于原始流形为e-平坦子流形，辅助流形上的分布属于原始流形，目标流形为原始流形的子流形，因此辅助流形上分布向目标流形m-投影是唯一的，于是通过上述一阶充分条件获得的参数即为m-投影点的参数。结合式(32)(33)和式(24)(26)，可以进一步整理得到m-投影结果为
[0109][0110][0111]
(3)更新辅助流形和目标流形上分布的参数
[0112]
信息几何方法的目标在于分别在辅助流形和目标流形上寻找近似后验分布的概率分布。所以目标流形上分布的是对的近似，辅助流形上分布的是对的近似。在目标流形上的投影结果为
[0113][0114]
其中并根据辅助流形和目标流形的定义可以得到是对cn(h)的近似。进一步可以得到辅助流形和目标流形上分布的参数更新方式
[0115][0116][0117]
为了使算法更加稳定，对上述迭代式进行松弛，辅助流形和目标流形上分布的参数的更新可以进一步表示为
[0118][0119]
[0120][0121]
其中0≤α≤1。
[0122]
进一步更新迭代次数t＝t 1，迭代式(34)(35)(39)(40)(41)至预设迭代次数或目标流形上分布的参数收敛。在上述迭代过程中，乘法运算仅涉及n维对角矩阵和n维向量的乘法、n维对角矩阵和n维对角矩阵的乘法以及标量和n维向量的乘法，上述运算的乘法数均为m，所以每次迭代的计算复杂度为远低于mmse估计的复杂度，从而支持同时估计大量用户的信道，有效地降低了导频开销。目标流形上分布的均值和方差即分别为空频波束域信道的后验均值和后验方差，具体为
[0123][0124][0125]
进一步，通过各用户空频波束域信道的非零元位置获取每个用户的空频波束域信道的后验均值和后验方差进一步通过采样空间舵矢量矩阵和采样频率舵矢量矩阵将各用户空频波束域信道的后验均值和后验方差转换为空间频率域信道的后验均值和后验方差具体为
[0126][0127][0128]
本领域技术人员可以理解的是，当信息几何方法应用于下行信道估计时，其计算过程与上行估计基本一致。此时，接收信号模型(3)中退化为仅包含单个用户，即下标k仅取1。接着，信道标准模型、信号处理过程与信息几何方法的应用与上行信道估计完全一致，得到信道的后验统计信息后，结合信道预测等方法，各用户终端将所获取的信道后验统计信息发送至基站。基站侧利用采样空间舵矢量矩阵和采样频率舵矢量矩阵将所获得的空频波束域信道的后验均值和后验方差转换为空间频率域信道的后验均值和后验方差。
[0129]
四、实施效果
[0130]
为了使本技术领域的人员更好地理解本发明方案，下面给出两种具体系统配置下的本实施例中信息几何方法与现有方法的信道估计性能结果比较。
[0131]
首先，给出本实施例中信息几何方法与现有方法的估计性能结果比较。对比的方法有mmse估计，文献“generalized approximate message passing for estimation with random linear mixing,in ieee isit,st.petersburg,russia,july 31-august 5,2011,pp.2168
–
2172.”中提出的gamp方法以及文献“unifying message passing algorithms under the framework of constrained bethe free energy minimization,ieee trans.wireless commun.,vol.20,no.7,pp.4144
–
4158,jul.2021.”中的算法2(下文简称vep)。考虑一配置为nr＝128，k＝48和n
t
＝1的大规模mimo系统，其中基站天线配置为n
r,v
＝8，n
r,h
＝16。图4给出了在所考虑大规模mimo系统上行链路下，本实施例中信息几何方法与gamp、vep方法在不同信噪比下的信道估计性能比较。信息几何方法、gamp以及vep的最大迭代次数均设置为100次。从图4中，可以看出在所有信噪比下，信息几何方法(iga)均可以获
得与mmse估计几乎相同的信道估计性能。当信道估计性能为-29db时，信息几何方法相较于gamp以及vep的信噪比增益大约为5db。这表明与gamp和vep方法相比，本实施例中的信息几何方法能够更加获得更加准确的信道估计性能。
[0132]
接着，给出本实施例中信息几何方法与现有方法的收敛曲线比较示意图。保持所考虑大规模各项参数不变，信噪比设置为20db，以mmse估计作为性能基线供参考，作出了信息几何方法、gamp以及vep方法的收敛曲线。从图5中，可以发现信息几何方法(iga)需要约200次迭代收敛，并且收敛后其估计性能与mmse估计几乎完全一致。而gamp以及vep方法需要超过550次迭代收敛。这表明相较于现有算法，信息几何方法的收敛速度更快。
[0133]
基于相同的发明构思，本发明实施例公开的一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，该计算机程序被加载至处理器时实现上述的适用于基站或用户终端的大规模mimo信道估计方法。
[0134]
在具体实现中，该设备包括处理器，通信总线，存储器以及通信接口。处理器可以是一个通用中央处理器(cpu)，微处理器，特定应用集成电路(asic)，或一个或多个用于控制本发明方案程序执行的集成电路。通信总线可包括一通路，在上述组件之间传送信息。通信接口，使用任何收发器一类的装置，用于与其他设备或通信网络通信。存储器可以是只读存储器(rom)或可存储静态信息和指令的其他类型的静态存储设备，随机存取存储器(ram)或者可存储信息和指令的其他类型的动态存储设备，也可以是电可擦可编程只读存储器(eeprom)、只读光盘(cd-rom)或其他光盘存储、盘存储介质或者其他磁存储设备、或者能够用于携带或存储具有指令或数据结构形式的期望的程序代码并能够由计算机存取的任何其他介质，但不限于此。存储器可以是独立存在，通过总线与处理器相连接。存储器也可以和处理器集成在一起。
[0135]
其中，存储器用于存储执行本发明方案的应用程序代码，并由处理器来控制执行。处理器用于执行存储器中存储的应用程序代码，从而实现上述实施例提供的信道估计方法。处理器可以包括一个或多个cpu，也可以包括多个处理器，这些处理器中的每一个可以是一个单核处理器，也可以是一个多核处理器。这里的处理器可以指一个或多个设备、电路、和/或用于处理数据(例如计算机程序指令)的处理核。
[0136]
基于相同的发明构思，本发明实施例公开的一种大规模mimo通信系统，包括基站和多个用户终端，其中基站通过上行探测获得各个用户终端的信道的先验统计信息，并通过接收到的上行导频信号以及先验统计信息，利用信息几何方法获取各个用户终端的后验统计信息。具体的信道估计信息几何方法参见前述实施例，此处不再赘述。
[0137]
基于相同的发明构思，本发明实施例公开的一种大规模mimo通信系统，包括基站和多个用户终端，其中用户终端通过下行信道探测获得各自信道的先验统计信息，并通过接收到的下行导频信号以及先验统计信息，利用信息几何方法以及信道预测方法获取各自信道的后验统计信息并反馈给基站。具体的信道估计信息几何方法参见前述实施例，此处不再赘述。
[0138]
基于相同的发明构思，本发明实施例公开的一种大规模mimo通信系统，包括基站和多个用户终端，其中基站或用户终端包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，该计算机程序被加载至处理器时实现前述的大规模mimo信道估计方法。在本技术所提供的实施例中，应该理解到，所揭露的方法，在没有超过本技术的精神
和范围内，可以通过其他的方式实现。当前的实施例只是一种示范性的例子，不应该作为限制，所给出的具体内容不应该限制本技术的目的。例如，一些特征可以忽略，或不执行。
[0139]
本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段，还包括由以上技术特征任意组合所组成的技术方案。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。