一种基于容积卡尔曼滤波的毫米波波束跟踪方法技术

本发明专利技术公开了一种基于容积卡尔曼滤波的毫米波波束跟踪方法，步骤一：建立通信信道模型；步骤二：给定滤波器初始条件，完成初始化；步骤三：计算容积点集；步骤四：进行状态预测；步骤五：更新状态；步骤六：阈值判断，检测跟踪的有效性。本发明专利技术在模拟波束成形系统下，利用容积变换解决均值和协方差的非线性传递问题，具有完全的数值稳定性，提高了波束跟踪精度。本发明专利技术提出的基于容积卡尔曼滤波的波束跟踪方法，在每次迭代过程中都会通过一定方式生成采样点集，与无迹卡尔曼滤波相比，点集更小，计算复杂度更低，效率更高。

【技术实现步骤摘要】
一种基于容积卡尔曼滤波的毫米波波束跟踪方法
本专利技术属于毫米波通信
，具体涉及一种基于容积卡尔曼滤波的毫米波波束跟踪方法。
技术介绍
5G作为新一代通信标准，与之前的4G相比，数据速率需要提高10倍，平均小区吞吐量需要提高25倍，系统容量需要提高1000倍。然而数据吞吐量的增长需要提供更大的带宽。5G之前使用的频段集中在6GHz以下，目前，该频段的频谱资源已十分紧张，因此，需要寻找新的频谱资源用于5G。毫米波凭借其拥有的巨大的免许可连续带宽(30～300GHz)以及高速数据传输能力引起人们的广泛关注。然而，毫米波的绕射能力较差，路径损耗较大，尤其在雨雪等恶劣天气下会产生极大衰减。鉴于此，人们提出采用大规模阵列天线下的波束成形技术来弥补毫米波传播中的路径损耗。由于该技术依靠特定方向上的高增益窄波束进行通信，发射机和终端用户间的相对运动会导致波束失配，影响通信质量。为了保持良好的通信性能，需要不断调整收发端的波束角度，实现波束实时对准。因此，有效的波束跟踪方法对于将毫米波应用在移动环境中具有重要意义。目前已有一些毫米波波束跟踪方法，如通过卡尔曼滤波算法进行跟踪。扩展卡尔曼滤波算法通过扫描所有可能的波束组合，建立扩展卡尔曼滤波算法中的测量矩阵，完成波束跟踪。但该算法需要进行波束全扫描，测量时间较长，很难在快速变化的场景中进行波束实时跟踪。一些基于扩展卡尔曼滤波算法的改进算法相对于普通的扩展卡尔曼滤波算法，只需要训练一个波束，极大降低了开销；同时在每个步骤中只使用一个测量方程，更加适合快速移动场景下的波束跟踪。然而两者采用的扩展卡尔曼滤波是通过对非线性系统的一阶线性化来解决非线性状态跟踪问题的，这会在线性化的真实后验均值和协方差中引入较大误差，导致滤波器的次优性能，有时还会出现滤波器发散。无迹卡尔曼滤波算法通过放弃以往对非线性函数进行线性化的方法，采用无迹变换处理均值和协方差的非线性传递克服了该问题。由于不需要线性化，不用忽略高阶项，因此具有较高的计算精度。但是无迹卡尔曼滤波的数值稳定性较差，随着状态维数增加，滤波精度会大幅下降。且无迹卡尔曼滤波算法的计算复杂度仍较高。因此，需要一种跟踪精度高，数值稳定性好，同时计算复杂度较低的算法来解决上述算法的问题。
技术实现思路
本专利技术所要解决的技术问题是针对上述现有技术的不足，提供一种基于容积卡尔曼滤波的毫米波波束跟踪方法，针对毫米波窄波束场景，以高精度和实时性为目标，实现快速准确的毫米波波束跟踪。为实现上述技术目的，本专利技术采取的技术方案为：一种基于容积卡尔曼滤波的毫米波波束跟踪方法，包括：步骤一：建立通信信道模型；步骤二：给定滤波器初始条件，完成初始化；步骤三：计算容积点集；步骤四：进行状态预测；步骤五：更新状态；步骤六：阈值判断，检测跟踪的有效性。为优化上述技术方案，采取的具体措施还包括：上述的步骤一建立的通信信道模型，发射端和接收端均配置均匀线性阵列，其信道模型表示为：其中，L表示信道路径的个数，αl[k]是第k个时刻第l条路径的信道复增益系数,服从一阶复高斯马尔科夫模型，表示为αl[k]＝ραl[k-1]+ζ[k-1]；其中ρ表示信道相关系数；和分别是第k个时刻第l条路径对应的接收端和发射端的阵列响应向量，其中和分别表示第l条路径的信号到达角和信号发射角，[·]H表示矩阵的共轭转置变换。上述的发射端和接收端的接收信号表示为：y[k]＝wHhHf+n(1)其中，y[k]表示第k个时刻的接收信号，h表示信道矩阵，f表示发射端模拟波束成形系统的波束成形向量，w表示接收端模拟波束成形系统的波束合并向量，n为均值为0、方差为R的加性高斯白噪声。上述的对于有Nr个阵元的接收端的均匀线性阵列，阵列响应向量表示为：对于有Nt个阵元的发射端的均匀线性阵列，阵列响应向量表示为：其中，表示阵元间距，λ为信号波长，[·]T表示矩阵的转置变换；波束成形向量f和波束合并向量w分别表示为上述的步骤二初始化中，首先定义容积卡尔曼滤波中的状态向量为：x[k]＝[αR[k],αI[k],θr[k],θt[k]]T(8)其中，θr[k],θt[k]分别表示信号到达角和信号发射角，α[k]＝αR[k]+jαI[k]，即αR[k],αI[k]别表示路径增益的实部和虚部；该向量的状态更新方程由高斯过程噪声驱动，表示为x[k]＝x[k-1]+u[k-1](9)其中，u[k-1]是高斯过程噪声；分别表示收发端的到达角、发射角的角度变化方差；其次，给定滤波器初始状态向量和协方差矩阵：x0＝E[x0](10)P0＝Σu(11)。上述的步骤三所述计算容积点集，包括：步骤3.1：根据k-1时刻的估计状态向量x[k-1]，计算容积点集为：Xi[k-1]＝S[k-1]ξi+x[k-1]，i＝1,…,m(12)步骤3.2：根据k-1时刻的协方差矩阵计算S[k-1]，公式如下：步骤3.3：定义ξi，公式如下：i＝1,…,m(14)其中，m＝2n，n为状态向量的维数。上述的步骤四所述进行状态预测，包括：步骤4.1：根据系统模型传播容积点，公式如下：Xi[k|k-1]＝Xi[k-1]，i＝1,…,m(15)步骤4.2：计算状态预测均值，公式如下：步骤4.3：计算预测协方差，公式如下：步骤4.4：根据系统模型进行容积点非线性变换，公式如下：Zi[k|k-1]＝g[Xi[k|k-1]]＝y[Xi[k|k-1]]-n，i＝1,…,m(18)步骤4.5：计算测量预测均值，公式如下：步骤4.6：计算预测误差协方差，公式如下：步骤4.7：计算预测互协方差，公式如下：上述的步骤五所述状态更新，包括：步骤5.1：计算卡尔曼增益，公式如下：K[k]＝PZXP-1ZZ(22)步骤5.2：更新状态向量，公式如下：x[k]＝x[k|k-1]+K[k](y[k]-Zi[k|k-1])(23)步骤5.3：更新协方差矩阵，公式如下：P[k]＝P[k|k-1]-K[k]PZZK[k]T(24)。上述的步骤六所述阈值判断，检测跟踪的有效性，具体为：当波束成形向量f和波束合并向量w的角度与CKF求得的最优估计值偏差超过阈值γ，即波束宽度的一半时，对两者的角度进行迭代更新，返回步骤三，直到迭代次数达到跟踪时刻T。本专利技术具有以下有益效果：1)本专利技术在模拟波束成形系统下，利用容积变换解决均值和协方差的非线性传递问题，具有完全的数值稳定性，提高了波束跟踪精度。2)本专利技术提出的基于容积卡尔曼滤波的波本文档来自技高网...

【技术保护点】
1.一种基于容积卡尔曼滤波的毫米波波束跟踪方法，其特征在于，包括：/n步骤一：建立通信信道模型；/n步骤二：给定滤波器初始条件，完成初始化；/n步骤三：计算容积点集；/n步骤四：进行状态预测；/n步骤五：更新状态；/n步骤六：阈值判断，检测跟踪的有效性。/n

【技术特征摘要】
1.一种基于容积卡尔曼滤波的毫米波波束跟踪方法，其特征在于，包括：
步骤一：建立通信信道模型；
步骤二：给定滤波器初始条件，完成初始化；
步骤三：计算容积点集；
步骤四：进行状态预测；
步骤五：更新状态；
步骤六：阈值判断，检测跟踪的有效性。


2.根据权利要求1所述的一种基于容积卡尔曼滤波的毫米波波束跟踪方法，其特征在于，所述步骤一建立的通信信道模型，发射端和接收端均配置均匀线性阵列，其信道模型表示为：



其中，L表示信道路径的个数，αl[k]是第k个时刻第l条路径的信道复增益系数,服从一阶复高斯马尔科夫模型，表示为αl[k]＝ραl[k-1]+ζ[k-1]；
其中ρ表示信道相关系数；

和分别是第k个时刻第l条路径对应的接收端和发射端的阵列响应向量，其中和分别表示第l条路径的信号到达角和信号发射角，[·]H表示矩阵的共轭转置变换。


3.根据权利要求2所述的一种基于容积卡尔曼滤波的毫米波波束跟踪方法，其特征在于，所述发射端和接收端的接收信号表示为：
y[k]＝wHhHf+n(1)
其中，y[k]表示第k个时刻的接收信号，h表示信道矩阵，f表示发射端模拟波束成形系统的波束成形向量，w表示接收端模拟波束成形系统的波束合并向量，n为均值为0、方差为R的加性高斯白噪声。


4.根据权利要求3所述的一种基于容积卡尔曼滤波的毫米波波束跟踪方法，其特征在于，对于有Nr个阵元的接收端的均匀线性阵列，阵列响应向量表示为：



对于有Nt个阵元的发射端的均匀线性阵列，阵列响应向量表示为：



其中，表示阵元间距，λ为信号波长，[·]T表示矩阵的转置变换；
波束成形向量f和波束合并向量w分别表示为








5.根据权利要求4所述的一种基于容积卡尔曼滤波的毫米波波束跟踪方法，其特征在于，所述步骤二初始化中，首先定义容积卡尔曼滤波中的状态向量为：
x[k]＝[αR[k],αI[k],θr[k],θt[k]]T(8)
其中，θr[k],θt[k]分别表示信号到达角和信号发射角，α[k]＝αR[k]+jαI[k]，即αR[k],αI[k]别表示路径增益的实部和虚部；
该向量的状态更新方程由高斯过程噪声驱动，表示为
x[k]＝x[k-1]+u[k-1](9)
其中，u[k-1]是高斯过程噪声；

【专利技术属性】
技术研发人员：仲伟志，李鹏辉，张俊杰，朱秋明，陈小敏，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32