一种信道追踪方法和装置与流程
本发明涉及通信信道处理领域，尤其涉及一种信道追踪方法和装置。
背景技术：
随着全球无线通信技术和互联网技术的蓬勃发展，移动用户数量呈现爆炸式增长，各种移动新业务不断涌现。未来的5g移动通信系统不仅需要支持语音功能，还必须支持数字、图像、多媒体等数据传输，在传输效率和覆盖范围等方面较4g提高一个或多个量级满足10年内移动互联网流量增加1000倍的发展需求。毫米波频段(30–300ghz)具有丰富的频谱资源，毫米波与大规模多输入多输出(multipleinputmultipleoutput；mimo)相结合可以带来巨大的空间增益、能量效率和频谱效率，毫米波大规模mimo成为了下一代移动通信的关键技术。
毫米波大规模mimo系统的性能取决于信道状态信息的准确度，信道状态信息的获取是当前毫米波大规模mimo研究的重点问题。现有的信道追踪方法一般是针对天线阵列为均匀线型阵列的移动通信系统，但是对于毫米波大规模mimo技术对应的移动通信系统而言，其用户终端和和基站端的天线都采用的均匀矩形阵列，所以实际中的信道追踪方法如果应用在毫米波大规模mimo技术上会存在较大误差，不够精确。
技术实现要素：
本发明的实施例提供一种信道追踪方法和装置，用于对设置有均匀矩形阵列的通信系统进行信道追踪。
为达到上述目的，本发明的实施例采用如下技术方案：
第一方面，提供一种信道追踪方法，包括：获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角、以及预设增益信息；
根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道时变公式计算当前时间块的信道矩阵的期望值；
根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设角度时变公式和预设信道矩阵公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值；
根据当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值；
根据当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。
可选的，根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道时变公式计算信道矩阵的期望值包括：
根据预设发射端的水平角和俯仰角、预设接收端的水平角和俯仰角和预设增益信息，依据预设信道矩阵公式计算初始信道矩阵；
根据初始信道矩阵和预设信道时变公式计算当前时间块的信道矩阵的期望值。
可选的，根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设角度时变公式和预设信道矩阵公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值包括：
根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道矩阵公式计算初始信道矩阵；
根据初始信道矩阵建立信号传输方程；
根据初始发射水平角、初始发射俯仰角、初始接收水平角和初始接收俯仰角生成初始二维角度；
根据初始二维角度、预设角度时变公式和信号传输方程，依据卡尔曼滤波算法计算当前时间块的二维角度的估计值；
根据二维角度的估计值和预设增益信息，依据预设信道矩阵公式计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值。
可选的，根据当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值包括：
根据当前时间块的信道矩阵的期望值、当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值，依据预设目标函数计算当前时间块的目标函数值；
根据预设水平阵列迭代公式和预设垂直阵列迭代公式分别对当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值进行更新，以获取当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值；
根据当前时间块的信道矩阵的期望值、当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值，依据预设目标函数对当前时间块的目标函数值进行更新，以获取当前时间块的目标函数的更新值；
当当前时间块的目标函数的更新值和当前时间块的目标函数值之差小于预设阈值时，确定当前时间块的水平阵列的更新值为当前时间块的水平阵列的最优值、当前时间块的垂直阵列的更新值为当前时间块的垂直阵列的最优值。
优选的，卡尔曼滤波算法具体为扩展卡尔曼滤波算法。
可选的，当当前时间块的目标函数的更新值和当前时间块的目标函数值之差不小于预设阈值时，将当前时间块的水平阵列的估计值更新为当前时间块的水平阵列的更新值，将当前时间块的垂直阵列的估计值更新为当前时间块的垂直阵列的更新值。
示例性的，预设信道矩阵公式具体为：
其中，h为信道矩阵，u为水平阵列、v为垂直阵列；
α[τ+1]＝ρα[τ]+ξ[τ]；
at(v)为信道发射端垂直阵列响应矩阵，ar(u)为信道接收端水平阵列响应矩阵，ar(v)为信道接收端垂直阵列响应矩阵，α为路径增益，u为水平方向波程差，v为垂直方向的波程差，
示例性的，预设信道时变公式具体为：
其中，h1(τ)为时间块τ的信道矩阵的期望值，β为信道时变参数，h(τ-1)为时间块τ的前一时间块的信道矩阵的期望值，v(τ)为均值为零的高斯白噪声且其方差为σv2，v(τ)与h(τ-1)统计独立。
示例性的，预设角度时变公式具体为：
θ(τ)＝θ(τ-1)+u(τ)；
其中，θ(τ)为时间块τ的二维角度，θ(τ-1)为时间块τ的前一时间块的二维角度，u(τ)为服从均值为零，方差为qu的高斯分布；
示例性的，信号传输方程具体为：
其中，y(τ)为时间块τ的信道接收端接收的信号，x为时间块τ的信道发射端发射的信号，θ(τ)为时间块τ的二维角度，n(τ)为传输噪声且n(τ)～cν(0,qn)，w为时间块τ的信道接收端的组合器矩阵，f为时间块τ的信道的模拟波束赋形矩阵，
示例性的，预设目标函数具体为：
其中，l为目标函数值，λ为正则项参数，uik为水平阵列中第i行第k个元素，vkj为垂直阵列中第j列第k个元素，hij为信道矩阵的期望值中第i行第j列的元素。
示例性的，预设水平阵列迭代公式具体为：
所述预设垂直阵列迭代公式具体为：
其中，ui'k为uik的更新值，vk'j为vkj的更新值，β为迭代步长。
第二方面，提供一种信道追踪装置，包括：获取模块、期望计算模块、估计计算模块和算法模块；
获取模块，用于获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角、以及预设增益信息；
期望计算模块，用于根据获取模块获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道时变公式计算当前时间块的信道矩阵的期望值；
估计计算模块，用于根据获取模块获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设角度时变公式和预设信道矩阵公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值；
算法模块，用于根据估计计算模块计算的当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和期望计算模块计算的当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值；
算法模块还用于根据算法模块获取的当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。
可选的，期望计算模块具体包括：第一计算单元和第二计算单元；
第一计算单元，用于根据获取模块获取的预设发射端的水平角和俯仰角、预设接收端的水平角和俯仰角和预设增益信息，依据预设信道矩阵公式计算初始信道矩阵；
第二计算单元，用于根据第一计算单元计算的初始信道矩阵和预设信道时变公式计算当前时间块的信道矩阵的期望值。
可选的，估计计算模块具体包括：第三计算单元、方程建立单元、角度估计单元和第四计算单元；
第三计算单元，用于根据获取模块获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道矩阵公式计算初始信道矩阵；
方程建立单元，用于根据初始信道矩阵建立信号传输方程；
角度估计单元，用于根据获取模块获取的初始发射水平角、初始发射俯仰角、初始接收水平角和初始接收俯仰角生成初始二维角度；
角度估计单元还用于根据初始二维角度、预设角度时变公式和方程建立单元建立的信号传输方程，依据卡尔曼滤波算法计算当前时间块的二维角度的估计值；
第四计算单元，用于根据角度估计单元计算的二维角度的估计值和获取模块获取的预设增益信息，依据预设信道矩阵公式计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值。
可选的，算法模块具体包括：第五计算单元、更新单元和判断单元；
第五计算单元，用于根据期望计算模块计算的当前时间块的信道矩阵的期望值，以及估计计算模块计算的当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值，依据预设目标函数计算当前时间块的目标函数值；
更新单元，用于根据预设水平阵列迭代公式和预设垂直阵列迭代公式分别对估计计算模块计算的当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值进行更新，以获取当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值；
第五计算单元还用于根据期望计算模块计算的当前时间块的信道矩阵的期望值以及更新单元计算的当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值，依据预设目标函数对第五计算单元计算的当前时间块的目标函数值进行更新，以获取当前时间块的目标函数的更新值；
当判断单元确定第五计算单元获取的当前时间块的目标函数的更新值和第五计算单元计算的当前时间块的目标函数值之差小于预设阈值时，确定更新单元计算的当前时间块的水平阵列的更新值为当前时间块的水平阵列的最优值、更新单元计算的当前时间块的垂直阵列的更新值为当前时间块的垂直阵列的最优值。
可选的，当判断单元确定第五计算单元获取的当前时间块的目标函数的更新值和第五计算单元计算的当前时间块的目标函数值之差不小于预设阈值时，第五计算单元还用于将当前时间块的水平阵列的估计值更新为更新单元计算的当前时间块的水平阵列的更新值，同时将当前时间块的垂直阵列的估计值更新为更新单元计算的当前时间块的垂直阵列的更新值。
本发明实施例提供的信道追踪方法和装置，该方法包括：获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角，以及预设增益信息；根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道时变公式计算当前时间块的信道矩阵的期望值；根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设角度时变公式和预设信道矩阵公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值；根据当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值；根据当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。本发明实施例提供的信道追踪方法中，针对均匀矩形分布的天线阵列，将其对应信道的收发端之间的信道矩阵分为水平阵列和垂直阵列进行考虑，在分别计算出信道矩阵在当前时间块的期望值、水平阵列和垂直阵列在当前时间块的估计值后，利用矩阵分解算法，对水平阵列和垂直阵列在当前时间块的估计值不断迭代更新，以使水平阵列和垂直阵列合成的信道矩阵在当前时间块的估计值和期望值的差距缩至最小即目标函数值收敛，则此时根据水平阵列和垂直阵列各自在当前时间块的估计值更新得到的水平阵列和垂直阵列在当前时间块的更新值，则为最符合信道追踪要求的信道矩阵在的最优值分解后可得的水平阵列和垂直阵列在当前时间块的最优值，自此完成信道追踪，得到最优信道矩阵。因此本发明实施提供的技术方案，适用于采用均匀矩形排列的天线阵列的毫米波大规模mimo系统；进一步的，因为本发明实施例提供的技术方案采用了矩阵分解算法修正了信道追踪结果，所以相对于现有的信道追踪方法，最终得到的追踪结果更为精确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的均匀矩阵阵列排布的天线排列示意；
图2为本发明实施例提供的一种信道追踪方法流程示意图；
图3为本发明实施例提供的另一种信道追踪方法流程示意图；
图4为本发明实施例提供的一种信道追踪装置结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明的是，本发明实施例中，“示例性的”或者“例如”等词用于表示作例子、例证或说明。本发明实施例中被描述为“示例性的”或者“例如”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言，使用“示例性的”或者“例如”等词旨在以具体方式呈现相关概念。
还需要说明的是，本发明实施例中，“的(英文：of)”，“相应的(英文：corresponding，relevant)”和“对应的(英文：corresponding)”有时可以混用，应当指出的是，在不强调其区别时，其所要表达的含义是一致的。
为了便于清楚描述本发明实施例的技术方案，在本发明的实施例中，采用了“第一”、“第二”等字样对功能和作用基本相同的相同项或相似项进行区分，本领域技术人员可以理解“第一”、“第二”等字样并不是在对数量和执行次序进行限定。
现有的信道追踪方法大多适用于天线排列为均匀线型阵列的通信系统，但实际中越来越普及的毫米波大规模mimo系统的移动场景中，基站端和用户端都使用参照图1所示的均匀矩形阵列排布的天线阵列，现有信道追踪方法并不适用。
针对上述问题，参照图2所示，本发明实施例提供一种信道追踪方法，包括：
201、获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角、以及预设增益信息。
示例性的，参照图1所示，水平角如
202、根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道矩阵公式和预设信道时变公式计算当前时间块的信道矩阵的期望值。
示例性的，因为在研究过程中一般只考虑los(lineofsight视距)用以分析信道矩阵，所以信道矩阵与水平角和俯仰角都有关系，预设信道矩阵公式具体为：
其中，h为信道矩阵，v为水平阵列、u为垂直阵列；
α[τ+1]＝ρα[τ]+ξ[τ]；
示例性的，信道的时间相关性服从一阶ar(autoregressivemodel，自回归模型)模型，所以预设信道时变公式具体为：
其中，h1(τ)为时间块τ的信道矩阵的期望值，β为信道时变参数，h(τ-1)为时间块τ的前一时间块的信道矩阵的期望值，v(τ)为均值为零的高斯白噪声且其方差为σv2，v(τ)与h(τ-1)统计独立；上述时间块τ指时间段，因为实际中毫米波大规模mimo通信系统的信道矩阵是块衰落的，即其变化是稳定一段时间后变化，而其中稳定的时间一般称为时间块。
203、根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设信道矩阵公式和预设角度时变公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值。
示例性的，预设角度时变公式：
θ(τ)＝θ(τ-1)+u(τ)；
其中，θ(τ)为时间块τ的二维角度，θ(τ-1)为时间块τ的前一时间块的二维角度，u(τ)为服从均值为零，方差为qu的高斯分布；方差
需要说明的是，202步骤和203步骤可以同时进行也可以分先后进行，此处不对此做具体限制。
204、根据当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值。
具体的，矩阵分解算法主要是通过迭代更新信道矩阵的估计值(具体指水平阵列的估计值和垂直阵列的估计值)，使其不断靠近信道矩阵的期望值，达到修正信道信息的目的，从而提高信道追踪得到的信道矩阵的精确度。
205、根据当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。
上述实施例提供的信道追踪方法，针对均匀矩形分布的天线阵列，将其对应信道的收发端之间的信道矩阵分为水平阵列和垂直阵列进行考虑，在分别计算出信道矩阵在当前时间块的期望值、水平阵列和垂直阵列在当前时间块的估计值后，利用矩阵分解算法，对水平阵列和垂直阵列在当前时间块的估计值不断迭代更新，以使水平阵列和垂直阵列合成的信道矩阵在当前时间块的估计值和期望值的差距缩至最小即目标函数值收敛，则此时根据水平阵列和垂直阵列各自在当前时间块的估计值更新得到的水平阵列和垂直阵列在当前时间块的更新值，则为最符合信道追踪要求的信道矩阵在的最优值分解后可得的水平阵列和垂直阵列在当前时间块的最优值，自此完成信道追踪，得到最优信道矩阵。因此本发明实施提供的技术方案，适用于采用均匀矩形排列的天线阵列的毫米波大规模mimo系统；进一步的，因为本发明实施例提供的技术方案采用了矩阵分解算法修正了信道追踪结果，所以相对于现有的信道追踪方法，最终得到的追踪结果更为精确。
参照图3所示，本发明实施例还提供另一种信道追踪方法用于对上述实施例提供的信道追踪方法的补充说明，具体包括：
301、获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角、以及预设增益信息。
302、根据预设发射端的水平角和俯仰角、预设接收端的水平角和俯仰角和预设增益信息，依据预设信道矩阵公式计算初始信道矩阵。
303、根据初始信道矩阵和预设信道时变公式计算当前时间块的信道矩阵的期望值。
304、根据初始信道矩阵建立信号传输方程。
示例性的，信号传输方程具体为：
其中，y(τ)为时间块τ的信道接收端接收的信号，x为时间块τ的信道发射端发射的信号，θ(τ)为时间块τ的二维角度，n(τ)为传输噪声且n(τ)～cν(0,qn)，w为时间块τ的信道接收端的组合器矩阵，f为时间块τ的信道的模拟波束赋形矩阵，
305、根据初始发射水平角、初始发射俯仰角、初始接收水平角和初始接收俯仰角生成初始二维角度。
306、根据初始二维角度、预设角度时变公式和信号传输方程，依据卡尔曼滤波算法计算当前时间块的二维角度的估计值。
具体的，因为实际中二维角度与接收信号是非线性关系，因此此处使用扩展卡尔曼滤波方法能够提高得到的二维角度的估计值的准确性。
307、根据二维角度的估计值和预设增益信息，依据预设信道矩阵公式计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值。
308、根据当前时间块的信道矩阵的期望值、当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值，依据预设目标函数计算当前时间块的目标函数值。
示例性的，预设目标函数具体为：
其中，l为目标函数值，λ为正则项参数，uik为水平阵列中第i行第k个元素，vkj为垂直阵列中第j列第k个元素，hij为信道矩阵的期望值中第i行第j列的元素；
事实上，在矩阵分解算法中，目标函数值原则上只需要计算信道矩阵得估计值和信道矩阵的期望值的差值即可，然后逐步迭代缩小差值至最小，但是实际中这样会使得最后求得的信道矩阵存在过拟合的问题，所以在本发明实施例中通过对水平阵列和垂直阵列进行正则化后加入到目标函数表达式里用来解决过拟合问题，即上式中的
309、根据预设水平阵列迭代公式和预设垂直阵列迭代公式分别对当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值进行更新，以获取当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值。
310、根据当前时间块的信道矩阵的期望值、当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值，依据预设目标函数对当前时间块的目标函数值进行更新，以获取当前时间块的目标函数的更新值。
311、判断当前时间块的目标函数的更新值和当前时间块的目标函数值之差是否小于预设阈值。
当当前时间块的目标函数的更新值和当前时间块的目标函数值之差小于预设阈值时，执行312；当当前时间块的目标函数的更新值和当前时间块的目标函数值之差不小于预设阈值时，执行313。
312、确定当前时间块的水平阵列的更新值为当前时间块的水平阵列的最优值、当前时间块的垂直阵列的更新值为当前时间块的垂直阵列的最优值。
313、将当前时间块的水平阵列的估计值更新为当前时间块的水平阵列的更新值，将当前时间块的垂直阵列的估计值更新为当前时间块的垂直阵列的更新值。
313步骤执行308。
314、根据当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。
具体的，h3＝u′*v′，其中h3为信道矩阵的最优值，u′为水平阵列的最优值，v′为垂直阵列的最优值。
本发明实施例提供的信道追踪方法，该方法包括：获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角，以及预设增益信息；根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道时变公式计算当前时间块的信道矩阵的期望值；根据初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设角度时变公式和预设信道矩阵公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值；根据当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值；根据当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。本发明实施例提供的信道追踪方法中，针对均匀矩形分布的天线阵列，将其对应信道的收发端之间的信道矩阵分为水平阵列和垂直阵列进行考虑，在分别计算出信道矩阵在当前时间块的期望值、水平阵列和垂直阵列在当前时间块的估计值后，利用矩阵分解算法，对水平阵列和垂直阵列在当前时间块的估计值不断迭代更新，以使水平阵列和垂直阵列合成的信道矩阵在当前时间块的估计值和期望值的差距缩至最小即目标函数值收敛，则此时根据水平阵列和垂直阵列各自在当前时间块的估计值更新得到的水平阵列和垂直阵列在当前时间块的更新值，则为最符合信道追踪要求的信道矩阵在的最优值分解后可得的水平阵列和垂直阵列在当前时间块的最优值，自此完成信道追踪，得到最优信道矩阵。因此本发明实施提供的技术方案，适用于采用均匀矩形排列的天线阵列的毫米波大规模mimo系统；进一步的，因为本发明实施例提供的技术方案采用了矩阵分解算法修正了信道追踪结果，所以相对于现有的信道追踪方法，最终得到的追踪结果更为精确。
为了更好的实施上述实施例中提供的信道追踪方法，参照图4所示，本发明实施例还提供一种信道追踪装置01，包括：获取模块41、期望计算模块42、估计计算模块43和算法模块44；
获取模块41，用于获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角、以及预设增益信息；
期望计算模块42，用于根据获取模块41获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道时变公式计算当前时间块的信道矩阵的期望值；
估计计算模块43，用于根据获取模块41获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设角度时变公式和预设信道矩阵公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值；
算法模块44，用于根据估计计算模块43计算的当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和期望计算模块42计算的当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值；
算法模块44还用于根据算法模块44获取的当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。
可选的，期望计算模块42具体包括：第一计算单元421和第二计算单元422；
第一计算单元421，用于根据获取模块41获取的预设发射端的水平角和俯仰角、预设接收端的水平角和俯仰角和预设增益信息，依据预设信道矩阵公式计算初始信道矩阵；
第二计算单元422，用于根据第一计算单元421计算的初始信道矩阵和预设信道时变公式计算当前时间块的信道矩阵的期望值。
可选的，估计计算模块43具体包括：第三计算单元431、方程建立单元432、角度估计单元433和第四计算单元434；
第三计算单元431，用于根据获取模块41获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道矩阵公式计算初始信道矩阵；
方程建立单元432，用于根据初始信道矩阵建立信号传输方程；
角度估计单元433，用于根据获取模块41获取的初始发射水平角、初始发射俯仰角、初始接收水平角和初始接收俯仰角生成初始二维角度；
角度估计单元433还用于根据初始二维角度、预设角度时变公式和方程建立单元432建立的信号传输方程，依据卡尔曼滤波算法计算当前时间块的二维角度的估计值；
第四计算单元434，用于根据角度估计单元433计算的二维角度的估计值和获取模块41获取的预设增益信息，依据预设信道矩阵公式计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值。
可选的，算法模块44具体包括：第五计算单元441、更新单元442和判断单元443；
第五计算单元441，用于根据期望计算模块42计算的当前时间块的信道矩阵的期望值，以及估计计算模块43计算的当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值，依据预设目标函数计算当前时间块的目标函数值；
更新单元442，用于根据预设水平阵列迭代公式和预设垂直阵列迭代公式分别对估计计算模块43计算的当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值进行更新，以获取当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值；
第五计算单元441还用于根据期望计算模块42计算的当前时间块的信道矩阵的期望值以及更新单元442计算的当前时间块的水平阵列的更新值和当前时间块的垂直阵列的更新值，依据预设目标函数对第五计算单元441计算的当前时间块的目标函数值进行更新，以获取当前时间块的目标函数的更新值；
当判断单元确定第五计算单元441获取的当前时间块的目标函数的更新值和第五计算单元441计算的当前时间块的目标函数值之差小于预设阈值时，确定更新单元442计算的当前时间块的水平阵列的更新值为当前时间块的水平阵列的最优值、更新单元442计算的当前时间块的垂直阵列的更新值为当前时间块的垂直阵列的最优值。
可选的，当判断单元443确定第五计算单元441获取的当前时间块的目标函数的更新值和第五计算单元441计算的当前时间块的目标函数值之差不小于预设阈值时，第五计算单元441还用于将当前时间块的水平阵列的估计值更新为更新单元442计算的当前时间块的水平阵列的更新值，同时将当前时间块的垂直阵列的估计值更新为更新单元442计算的当前时间块的垂直阵列的更新值。
本发明实施例提供的信道追踪装置，因为该装置包括：获取模块，用于获取信道发射端发射信号的初始发射水平角和初始发射俯仰角、信道接收端接收信号的初始接收水平角和初始接收俯仰角、以及预设增益信息；期望计算模块，用于根据获取模块获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角和预设增益信息，依据预设信道时变公式计算当前时间块的信道矩阵的期望值；估计计算模块，用于根据获取模块获取的初始发射水平角、初始发射俯仰角、初始接收水平角、初始接收俯仰角、预设增益信息、预设角度时变公式和预设信道矩阵公式，计算当前时间块的水平阵列的估计值和当前时间块的垂直阵列的估计值；算法模块，用于根据估计计算模块计算的当前时间块的水平阵列的估计值、当前时间块的垂直阵列的估计值和期望计算模块计算的当前时间块的信道矩阵的期望值，依据矩阵分解算法获取当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值；算法模块还用于根据算法模块获取的当前时间块的水平阵列的最优值和当前时间块的垂直阵列的最优值，依据预设信道矩阵公式计算当前时间块的信道矩阵的最优值。所以本发明实施例提供的技术方案在对信道进行追踪时，可以针对均匀矩形分布的天线阵列，将其对应信道的收发端之间的信道矩阵分为水平阵列和垂直阵列进行考虑，在分别计算出信道矩阵在当前时间块的期望值、水平阵列和垂直阵列在当前时间块的估计值后，利用矩阵分解算法，对水平阵列和垂直阵列在当前时间块的估计值不断迭代更新，以使水平阵列和垂直阵列合成的信道矩阵在当前时间块的估计值和期望值的差距缩至最小即目标函数值收敛，则此时根据水平阵列和垂直阵列各自在当前时间块的估计值更新得到的水平阵列和垂直阵列在当前时间块的更新值，则为最符合信道追踪要求的信道矩阵在的最优值分解后可得的水平阵列和垂直阵列在当前时间块的最优值，自此完成信道追踪，得到最优信道矩阵。因此本发明实施提供的技术方案，适用于采用均匀矩形排列的天线阵列的毫米波大规模mimo系统；进一步的，因为本发明实施例提供的技术方案采用了矩阵分解算法修正了信道追踪结果，所以相对于现有的信道追踪方法，最终得到的追踪结果更为精确。
本申请实施例提供一种计算机设备，包括存储器、处理器。该存储器上存储有可在处理器上运行的计算机程序，该处理器执行上述计算机程序时实现前述的信道追踪方法。其中，上述存储介质包括：rom、ram、磁碟或者光盘等各种可以存储程序代码的介质。
本申请实施例提供一种计算机可读介质，其存储有计算机程序，该计算机程序被处理器执行时实现前述的信道追踪方法。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
- 还没有人留言评论。精彩留言会获得点赞！