一种面向车联网的动态MIMO信道建模及参数计算方法与流程

本发明涉及车联网移动通信信道建模仿真，属于无线通信领域，具体是一种考虑复杂交通状况下车速随机变化的动态mimo信道建模方法，以及信道参数计算方法。

背景技术：

随着道路交通密度的快速增加，道路安全和交通流量管理问题亟待解决。车辆到车辆(vehicle-to-vehicle,v2v)通信系统能够有效提高道路交通的安全性和便利性，成为第5代移动通信所关注的关键传输技术之一。多输入多输出技术(multipleinputmultipleoutput，mimo)能够在不增加带宽和传输功率的前提下，实现高速数据传输，显著提高了通信系统容量与可靠性。因此，构建符合实际交通状况的v2v动态mimo信道模型，对v2v系统的仿真和性能评估以及系统的改进设计具有重要意义。

在v2v通信系统中，由于收发端和散射体均处于快速运动状态，导致其信道呈现明显的非平稳特性。目前，大部分v2v信道模型假设移动端以恒定速度沿特定方向移动，然而，在实际交通环境中，由于周围车辆和设施以及交通灯等的影响，车辆在移动过程中会经历加速或者减速的过程，并且在转弯和路面不平地段，车辆会改变其运动方向。运动速度大小和方向的变化会对信道特性产生重要影响，但考虑车辆速度变化的v2v信道模型还很少。另一方面，目前大部分非平稳信道建模方法没有考虑到不同场景段信道状态之间的连续性，导致输出多普勒频率和相位与实际不吻合。因此，有必要研究一种符合实际交通状况并能输出精确多普勒频率的v2v动态mimo信道建模及参数演进方法。

技术实现要素：

为了准确分析及研究复杂交通状况下的v2v动态mimo信道模型及特性，本发明提出一种面向车联网的动态mimo信道建模及参数计算方法。

本发明所采用的技术方案有：一种面向车联网的动态mimo信道建模及参数计算方法，步骤如下：

第一步：针对复杂交通状况下的v2v动态mimo通信场景，将收发车辆之间的mimo信道建模为

其中，s表示发送端天线数目，u表示接收端天线数目，发送端第s根天线与接收端第u根天线之间的信道hu,s(t,τ)表示为不同功率和时延的多条传播路径信号的叠加，将第n条路径的归一化信道衰落系数建模为

其中，m表示第n条路径的散射支路数，表示由于收发端和散射体位置变化引起的多普勒相移，表示由天线位置和散射环境变化引起的相移，表示随机初始相移；

第二步：实时计算第n条路径第m条散射支路的多普勒相移采用方法如下：

1)实时计算移动发射端即mt和移动接收端即mr的位置矢量，方法如下：

步骤一：假设初始时刻mt和mr的位置矢量分别为lmt(t0)，lmr(t0)；

步骤二：迭代计算t+δt时刻mt和mr的位置矢量，方法如下

其中，vmt(t),vmr(t)分别表示mt和mr的移动速度矢量。

步骤三：令t＝t+δt，返回步骤二，计算下一时刻mt和mr的位置矢量，直到仿真结束；

2)假设收发端之间第n条散射路径的第一个和最后一个散射体分别记为an,zn，实时计算an,zn的位置矢量，采用方法如下：

步骤一：假设初始时刻散射体an,zn的位置矢量分别为

步骤二：迭代计算t+δt时刻散射体an,zn的位置矢量，方法如下

其中，表示散射体an以及zn的移动速度矢量。

步骤三：令t＝t+δt，返回步骤二，计算下一时刻散射体an,zn的位置矢量，直到仿真结束；

3)计算mt与散射体an之间的信号角度均值以及散射体zn与mr之间的信号角度均值，方法如下：

其中，(·)^x，(·)^y，(·)^z分别表示(·)矢量的横、纵和竖坐标值，表示发送端与散射体an之间的距离，表示接收端与散射体zn之间的距离；

4)根据用户预先指定或实测的角度统计分布，动态计算并产生m组时变角度值θmt,n,m(t)，θmr,n,m(t)，φmt,n,m(t)，φmr,n,m(t)，并实时构建第n条路径第m条散射支路收发端方向矢量，方法如下：

5)实时计算第n条路径第m条散射支路的多普勒相移，方法如下：

其中，k＝2πf0/c表示波数，f0和c分别为载波频率和光速；

第三步：实时计算第n条路径第m条散射支路的相移采用方法如下：

1)计算收发端天线位置旋转矩阵，方法如下：

其中，分别表示发送端移动速度的方位角和俯仰角，分别表示接收端移动速度的方位角和俯仰角；

2)计算第n条路径第m条散射支路的相移方法如下：

其中，dmt,s(t0)，dmr,u(t0)分别表示发送端第s根天线和接收端第u根天线初始时刻的位置；

第四步：实时产生第n条路径第m条散射支路的随机初始相移产生方法如下，在[0,2π)范围内产生m个服从均匀分布的随机变量

本发明具有如下有益效果：

(1)、本发明提出的v2v动态mimo信道建模方法，综合考虑了收发车辆移动、散射体移动、运动速度及轨迹变化等多种影响电波传播的因素。

(2)、本发明提出的v2v动态mimo信道建模方法和参数计算方法，支持时变的信道模型参数，适用于任意三维天线阵列及任意车辆运动轨迹的传播场景，并能够保证输出信道衰落相位的连续性。

附图说明：

图1为本发明方案中复杂交通状况下的v2v动态mimo通信场景。

图2为采用本发明方案产生的二维角度分布演进过程。

图3为采用本发明方案产生的信道时变衰落包络及相位。

具体实施方式：

下面结合附图对本发明作进一步的说明。

本发明一种面向车联网的动态mimo信道建模及参数计算方法，具体步骤如下：

第一步：如图1所示，针对复杂交通状况下的v2v动态mimo通信场景，本发明将收发车辆之间的mimo信道建模为

其中，s表示发送端天线数目，u表示接收端天线数目，发送端第s根天线与接收端第u根天线之间的信道hu,s(t,τ)可表示为不同功率和时延的多条传播路径信号的叠加，本发明将第n条路径的归一化信道衰落系数建模为

其中，m表示第n条路径的散射支路数，表示由于收发端和散射体位置变化引起的多普勒相移，表示由天线位置和散射环境变化引起的相移，表示随机初始相移。

第二步：实时计算第n条路径第m条散射支路的多普勒相移本发明采用方法如下：

1)实时计算移动发射端(mobiletransmitter，mt)和移动接收端(mobilereceiver，mr)的位置矢量，方法如下：

步骤一：假设初始时刻mt和mr的位置矢量分别为lmt(t0)，lmr(t0)；

步骤二：迭代计算t+δt时刻mt和mr的位置矢量，方法如下

其中，vmt(t),vmr(t)分别表示mt和mr的移动速度矢量。

步骤三：令t＝t+δt，返回步骤二，计算下一时刻mt和mr的位置矢量，直到仿真结束。

2)假设收发端之间第n条散射路径的第一个和最后一个散射体分别记为an,zn，实时计算an,zn的位置矢量，本发明采用方法如下：

步骤一：假设初始时刻散射体an,zn的位置矢量分别为

步骤二：迭代计算t+δt时刻散射体an,zn的位置矢量，方法如下

其中，表示散射体an以及zn的移动速度矢量。

步骤三：令t＝t+δt，返回步骤二，计算下一时刻散射体an,zn的位置矢量，直到仿真结束。

3)计算mt与散射体an之间的信号角度均值以及散射体zn与mr之间的信号角度均值，方法如下：

其中，(·)^x，(·)^y，(·)^z分别表示(·)矢量的横、纵和竖坐标值，表示发送端与散射体an之间的距离，表示接收端与散射体zn之间的距离。

4)根据用户预先指定或实测的角度统计分布，动态计算并产生m组时变角度值θmt,n,m(t)，θmr,n,m(t)，φmt,n,m(t)，φmr,n,m(t)，并实时构建第n条路径第m条散射支路收发端方向矢量，方法如下：

5)实时计算第n条路径第m条散射支路的多普勒相移，方法如下：

其中，k＝2πf0/c表示波数，f0和c分别为载波频率和光速。

第三步：实时计算第n条路径第m条散射支路的相移本发明采用方法如下：

1)计算收发端天线位置旋转矩阵，方法如下：

其中，分别表示发送端移动速度的方位角和俯仰角，分别表示接收端移动速度的方位角和俯仰角。

2)计算第n条路径第m条散射支路的相移方法如下：

其中，dmt,s(t0)，dmr,u(t0)分别表示发送端第s根天线和接收端第u根天线初始时刻的位置。

第四步：实时产生第n条路径第m条散射支路的随机初始相移产生方法如下，在[0,2π)范围内产生m个服从均匀分布的随机变量

下面具体通过附图和实施例来说明本发明一种面向车联网的动态mimo信道建模及参数计算方法，本案例选用2x2mimo系统，仿真时长为5s，采样率fs＝10khz，信道状态更新间隔δt＝0.1s。本实例中，无线信道传播场景采用winner+信道标准中的城市微蜂窝场景，具体仿真参数如表1所示。

表1仿真参数

本实施例包括以下步骤：

1、计算第n条路径第m条散射支路的多普勒相移过程如下：

1)实时计算mt和mr的位置矢量，方法如下：

步骤一：初始时刻mt和mr的位置矢量分别表示为lmt(t0)＝[0,0,0]，lmr(t0)＝[0,0,0]；

步骤二：计算t+δt时刻mt和mr的位置矢量lmt(t+δt)，lmr(t+δt)；

步骤三：令t＝t+δt且返回步骤二，计算下一时刻mt和mr的位置矢量，直到仿真结束。

2)令收发端之间第n条散射路径的第一个和最后一个散射体为an,zn，并实时计算an,zn的位置矢量，方法如下：

步骤一：初始时刻散射体an,zn的位置矢量分别表示为

步骤二：计算t+δt时刻散射体an,zn的位置矢量

步骤三：令t＝t+δt且返回步骤二，计算下一时刻散射体an,zn的位置矢量，直到仿真结束。

3)计算mt与散射体an的信号角度均值以及散射体zn与mr的信号角度均值，本案例中，初始时刻第1条路径的角度均值为

4)本发明假设角度服从vonmisesfisher(vmf)统计分布，动态计算并产生1000组集中程度因子κ＝100的时变vmf角度值θmt,n,m(t)，θmr,n,m(t)，φmt,n,m(t)，φmr,n,m(t)。本案例中，初始时刻第1条路径第1条散射支路的角度为θmt,1,1(t0)＝0.2763，θmr,1,1(t0)＝0.3742，φmt,1,1(t0)＝0.6199，φmr,1,1(t0)＝-0.5773；实时计算第n条路径第m条散射支路相对收发端的信号单位矢量，本案例中，初始时刻第1条路径第1条散射支路相对收发端的信号单位矢量为smt,1,1(t0)＝[0.7830；0.5589；0.2728]，smr,1,1(t0)＝[0.7799；-0.5080；0.3655]。

5)实时计算第n条路径第m条散射支路的多普勒相移本案例中，第1条路径第1条散射支路为例，令初始时刻利用本发明方案仿真结束时刻(t＝5s)，

2、实时计算第n条路径第m条散射支路的相移过程如下：

1)计算收发端天线位置旋转矩阵本案例中，初始时刻收发端位置旋转矩阵为

2)计算第n条路径第m条散射支路的相移其中，四个子信道在初始时刻第1条路径第1条散射支路的相移分别为

3、实时产生第n条路径第m条散射支路的随机初始相移产生方法如下，在[0,2π)范围内产生m个服从均匀分布的随机变量本案例中，第1条路径第1条散射支路的随机初始相移为

4、将上述步骤获得的时变相位参数代入式(2)的信道模型中，最终获得该场景下v2v动态mimo信道。

本实施例得到的效果可以通过图2-图3仿真实验中所获得的具体数据进行进一步说明。我们看到：1)图2(a)和图2(b)分别给出了mt端和mr端的二维角度参数演进过程，从图中可以看出，角度均值随时间变化，而角度分布不变；2)图3给出了四条子信道衰落的包络和相位，可以看出本发明产生的信道衰落是时变的，且输出信道衰落相位连续。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下还可以作出若干改进，这些改进也应视为本发明的保护范围。