CN113949439B - 一种面向无人机空对空通信的几何随机信道建模方法 - Google Patents

Abstract

本发明公开了一种面向无人机空对空通信的几何随机信道建模方法，该方法首先考虑了无人机收发端的三维任意变速移动，并通过通信场景计算通信系统中的有效散射区域、散射簇以及散射点；其次，引入了马尔可夫生灭过程来描述簇在时间轴上的消失和演化；最后，通过计算复杂信道冲激响应，用于计算分析信道统计特性；本发明综合考虑了无人机收发双端的变速任意轨迹移动，使用生灭过程模拟信道变化，通过随机分布生成角度参数，能够保证所建立的模型可以有效模拟信道非平稳特性，符合真实通信场景。

Description

一种面向无人机空对空通信的几何随机信道建模方法

技术领域

本发明涉及信道建模技术领域，特别是涉及一种面向无人机空对空通信的几何随机信道建模方法。

背景技术

无人机由于其成本较低，使用周期长，易于保养等优势在众多领域中得到了广泛的应用。在无人机通信方面，由于其可以灵活部署，操作简便等优势，被用于应急通信，临时基站等。由于无人机的快速发展，无人机独特属性也使无人机信道区别于旧有传统信道特性，因此需要建立一个合适的三维动态几何随机模型用于模拟仿真无人机信道，从而指导无人机通信的系统设计和优化。

在无人机通信场景中，由于无人机的快速随机移动，导致无人机的通信环境在不断的变化，不同时刻会有不同的信道状态且变化较大，因此无人机信道具有非平稳特性。一般将无人机信道根据收发端不同分为无人机空对空信道和无人机空对地信道，目前大部分研究集中分析无人机空对地信道，关于无人机空对空信道也是一个较为重要的应用场景，因此建立一个通用准确的无人机空对空通信几何随机信道模型是必要的。

发明内容

有鉴于此，本发明的目的在于提供一种面向无人机空对空通信的几何随机信道建模方法，用以建立一个通用准确的无人机空对空通信几何随机信道模型，准确分析描述无人机空对空通信链路的情况。

为了实现上述目的，本发明采用如下技术方案：

一种面向无人机空对空通信的几何随机信道建模方法，包括如下步骤：

步骤S1、获取无人机的接收端以及发送端在初始时刻的实际位置和实际速度，通过构建椭球模型，对所述实际速度进行变换，得到所述接收端以及发送端的相对速度；

步骤S2、使用马尔可夫生灭过程模拟由于收发端快速移动导致的簇的生灭过程，并且随机生成某一时刻的到达角，并通过到达角随机生成散射点位置坐标及其时变坐标；

步骤S3、根据所述接收端和发送端的坐标以及散射点的时变坐标计算出时变时延参数、时变功率参数以及时变角度参数，其中，所述时变时延参数包括视距成分时延和非视距成分时延，所述时变角度参数包括视距分量的到达角、视距分量的离开角、非视距分量的到达角、非视距分量的离开角；

步骤S4、根据步骤S3中的时变时延参数，时变功率参数以及时变角度参数计算信道冲激响应。

进一步的，所述步骤S1具体包括：

步骤S101、获取所述接收端以及发送端在初始时刻的实际位置和实际速度；

步骤S102、根据步骤S101中获取的实际位置和实际速度，并且以所述接收端为参考点，构建椭球模型，再确定所述接收端以及发送端在模型中的坐标，其中，所述接收端在该模型的坐标系中，x轴以及y轴的值均为0，z轴的值为其离地高度；

步骤S103、根据步骤S101中获取的实际位置和实际速度，将初始时刻的接收端以及发送端进行速度转换，使得所述接收端在模型中为竖直上下变速运动，所述发送端为相对于接收端的任意方向的变速移动，再获得所述接收端以及发送端的相对速度。

进一步的，所述步骤S2具体包括：

步骤S201、设定初始时刻的簇的数量，并且根据接收端和发送端之间的相对距离以及随机分布生成簇的时延参数，再通过该时延参数计算得到簇的几何参数，其中，

所述簇的时延参数的表达式为：

公式中，τ_LOS(t)＝D_LOS(t)/c，τ_LOS(t)为视距成分时延，D_LOS(t)为接收端以及发送端之间的距离，c为光速，服从指数分布；

所述簇的几何参数的表达式为：

a′＝τ_l*c/2

步骤S202、为每一个簇随机生成若干对水平到达角和垂直到达角，其中通过每一对到达角均获得一个散射点的位置坐标；

步骤S203、根据所述接收端以及发送端的坐标，计算无人机通信场景的地面有效散射区域，并根据所述地面有效散射区域以及随机生成的水平到达角和垂直到达角，计算对应的散射点的位置坐标；设定所述散射点为相对地面静止，根据接收端的实际速度将所述散射点转换为相对接收端的移动散射点，并且计算该移动散射点的时变坐标；其中，

所述地面有效散射区域的表达式为：

公式中，b'²＝a'²-f'²，(x,y)为地面有效散射区域上的任意一点坐标，(x_c,y_c)表示地面有效散射区域的中心位置坐标，a,b为地面有效散射区域建模为椭圆的长半轴和短半轴长度，θ为收发端水平面投影点的夹角，以接收端为水平面坐标原点，则发送端的坐标为(x_T,y_T)，则/>h_T表示为发送端在坐标系中z轴的值，h_R表示为接收端在坐标系中z轴的值，f'表示为接收端与发送端之间距离的一半；

所述散射点的时变位置坐标S_l,n(t)在每个任意时刻t的表达式为：

公式中，(x_l,n,y_l,n,z_l,n)为散射点的坐标；

表示为水平到达角，/>表示为垂直到达角；

所述时变坐标S_l,n(t)的表达式为：

S_l,n(t)＝S_l,n(t-Δt)-v_RΔt

公式中，v_R表示实际速度，Δt指一段极小的时间，在这段时刻假定速度不变，散射簇固定不变；

步骤S204、首先根据所述接收端以及发送端的相对速度，计算其对应的移动场景参数，并且设定生成概率和消亡概率，然后通过马尔可夫生灭过程模拟簇的演进过程，并且计算得到新生成簇的数学期望，最后将经过马尔可夫生灭过程模拟后，新生成的簇代入到步骤S201中且继续执行S201-步骤S204，计算相应的时延参数和相应的近地散射点坐标以及该近地散射点对应的移动散射点的时变坐标；其中，

所述新生成簇的数学期望的表达式为：

公式中，D_c为基于不同场景的相关因子，χ_p(t,Δt)为所述的移动场景参数，λ_R为所述消亡概率，λ_G为所述生成概率；

所述移动场景参数的的表达式为：

χ_p(t,Δt)＝χ_T(t,Δt)+χ_R(t,Δt)

公式中，v'_R(t),v'_T(t)表示转换后的相对速度。

进一步的，所述步骤S3具体包括：

步骤S301、根据接收端以及发送端和散射点的距离计算得到视距成分时延和非视距成分时延，其中，视距成分时延由接收端和发送端的距离获得，表达式为：

τ_LOS(t)＝D_LOS(t)/c

公式中，D_LOS(t)为接收端以及发送端之间的距离，c为光速；

非视距成分时延通过选择每一个簇上任意一个散射点，计算其到接收端以及发送端的距离和求得，表达式为：

τ_l(t)＝(D_l,R(t)+D_l,T(t))/c

公式中，D_l,R/T(t)为散射点到R或者T的距离，其中，R表示为接收端，T表示为发送端；

步骤S302、根据步骤S301中求得的非视距成分时延，计算得到时变功率参数，表达式为：

公式中，τ_l(t)表示非视距成分时延，r_τ在LoS分量中为3，在NLoS分量中为2.1，σ为延迟扩展，其对数值lg(σ)服从正态分布，其值服从(0,ζ²)的高斯分布，ζ＝3；

为了让簇的功率和为1，所以对簇的功率进行归一化处理，表达式为：

公式中，L(t)为时变的有效散射簇数量；

步骤S303、根据不同时刻的接收端和发送端的坐标以及散射点的时变坐标，计算角度参数，其中，

视距分量的到达角的表达式为：

公式中，D_LOS(t)表示为接收端以及发送端之间的距离，为视距分量的水平到达角，/>为视距分量的垂直到达角；

视距分量的离开角的表达式为：

公式中，为视距分量的水平离开角，/>为视距分量的垂直离开角；

非视距分量的到达角的表达式为：

公式中，为第l个散射簇上的第n个散射点散射分量到接收端的水平到达角，为第l个散射簇上的第n个散射点散射分量到接收端的垂直到达角，(x_l,n,y_l,n,z_l,n)为第l个散射簇上的第n个散射点坐标，D_l,R(t)为第l个散射簇到接收端的距离；

非视距分量的离开角的表达式为：

公式中，为第l个散射簇上的第n个散射点散射分量到发射端的水平离开角，为第l个散射簇上的第n个散射点散射分量到发射端的垂直离开角，D_l,T(t)为第l个散射簇到发射端的距离。

进一步的，步骤S4具体步骤如下：

步骤S401、计算接收端第p根天线和发送端第q根天线间的衰落信道的冲激响应，表达式为：

公式中，P_l(t)为时变功率参数，h_l,pq(t)为第l个抽头的系数，τ_l为第l个抽头的时延参数，δ(τ-τ_l)表示时延的冲激函数，τ为时延变量，L(t)表示随着时间变化的有效散射簇数量；

步骤S402、计算视距信道和非视距信道的抽头系数，其中，

视距分量系数的表达式为：

非视距分量系数的表达式为：

公式中，l>2，λ为传输波长，r_p为接收端R的第p根天线相对接收端R天线阵列中心的相对距离，r_p为发送端T的第q根天线相对发送端T天线阵列中心的相对距离，θ_l,n为随机相移，是指每个散射射线自己随机生成的相移，是在(0，2π]上均匀分布的随机变量；Φ_LOS(t)和/>表示发射角，Ψ_LOS(t)和/>代表到达角，f_Tm表示发送端的最大多普勒频移，f_Rm表示接收端最大多普勒频移。

本发明的有益效果是：

本发明综合考虑了无人机收发双端的变速任意轨迹移动，使用生灭过程模拟信道变化，通过随机分布生成角度参数，能够保证所建立的模型可以有效模拟信道非平稳特性，符合真实通信场景。本发明适用于多个频段场景，具有一定程度的普适性。

附图说明

图1为实施例1中一种面向无人机空对空通信的几何随机信道建模方法的流程示意图。

图2为实施例1中椭球模型的示意图，其中，R表示为接收端，T表示为发送端，h_R和h_T为R、T离地高度，v_R和v_T依次表示R、T的实际速度，sca象征任意一个散射分量为scatter的前三个字母，los表示视距分量。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1

参见图1和图2，本实施例提供一种面向无人机空对空通信的几何随机信道建模方法，该方法考虑了无人机收发双端的三维变速任意移动，模型使用马尔可夫生灭过程模拟了簇的演进过程，从而保证可以有效模拟无人机空对空信道的时域非平稳性。

该方法具体包括：

步骤S1、获取无人机的接收端以及发送端在初始时刻的实际位置和实际速度，通过构建椭球模型，对实际速度进行变换，得到接收端以及发送端的相对速度；

步骤S1具体包括：

步骤S101、获取接收端R以及发送端T在初始时刻的实际位置和实际速度；

步骤S102、根据步骤S101中获取的实际位置和实际速度，并且以接收端为参考点，构建椭球模型，再确定接收端R以及发送端T在模型中的坐标，其中，接收端R在该模型的坐标系中，x轴以及y轴的值均为0，z轴的值为其离地高度；

更具体的说，在本实施例中，设定R端的坐标为(0,0,h_R)，那么T端的坐标则为(x_T,y_T,h_T)，h_R和h_T为R、T离地高度，x_T和y_T为T到R的距离。

步骤S103、根据步骤S101中获取的实际位置和实际速度，将初始时刻的接收端以及发送端进行速度转换，使得接收端在模型中为竖直上下变速运动，发送端为相对于接收端的任意方向的变速移动，再获得接收端以及发送端的相对速度。

更具体的说，在本实施例中，可以通过无人机收发端的移动的水平角和俯仰角来描述无人机收发端的移动，γ_R和γ_T来表示R、T的水平角，δ_R和δ_T来表示R、T的俯仰角，具体是：

v_R＝v_R[cosγ_Rcosδ_R sinγ_Rcosδ_R sinδ_R]

v_T＝v_T[cosγ_Tcosδ_T sinγ_Tcosδ_T sinδ_T]

再将无人机收发端的移动转换为相对R的移动，R保留垂直方向移动，表达式为：

v'_R＝v_R[0 0 sinδ_R]

v'_T＝[v_Tcosγ_Tcosδ_T-v_Rcosγ_Rcosδ_R v_Tsinγ_Tcosδ_T-v_Rsinγ_Rcosδ_R v_Tsinδ_T]

步骤S2、使用马尔可夫生灭过程模拟由于收发端快速移动导致的簇的生灭过程，并且随机生成某一时刻的到达角，并通过到达角随机生成散射点位置坐标及其时变坐标；

步骤S2具体包括：

步骤S201、设定初始时刻的簇的数量，并且根据接收端和发送端之间的相对距离以及随机分布生成簇的时延参数，再通过该时延参数计算得到簇的几何参数，其中，

所述簇的时延参数的表达式为：

公式中，τ_LOS(t)＝D_LOS(t)/c，τ_LOS(t)为视距成分时延，D_LOS(t)为接收端以及发送端之间的距离，c为光速，服从指数分布；

所述簇的几何参数的表达式为：

a′＝τ_l*c/2

步骤S202、使用von-Mises分布和cosine分布，为每一个簇随机生成若干对水平到达角和垂直到达角，其中通过每一对到达角均获得一个散射点的位置坐标；

步骤S203、根据所述接收端以及发送端的坐标，计算无人机通信场景的地面有效散射区域，并根据所述地面有效散射区域以及随机生成的水平到达角和垂直到达角计算对应的散射点的位置坐标S_l,n(t)；

设定所述散射点为相对地面静止，根据接收端的实际速度将所述散射点转换为相对接收端的移动散射点，并且计算该移动散射点的时变坐标；其中，

所述地面有效散射区域的表达式为：

公式中，b'²＝a'²-f'²，(x,y)为地面有效散射区域上的任意一点坐标，(x_c,y_c)表示地面有效散射区域的中心位置坐标，a,b为地面有效散射区域建模为椭圆的长半轴和短半轴长度，θ为收发端水平面投影点的夹角，以接收端为水平面坐标原点，则发送端的坐标为(x_T,y_T)，则/>h_T表示为发送端在坐标系中z轴的值，h_R表示为接收端在坐标系中z轴的值，f'表示为接收端与发送端之间距离的一半；

所述散射点的位置坐标S_l,n(t)在每个任意时刻t的表达式为：

公式中，(x_l,n,y_l,n,z_l,n)为散射点的坐标；

表示为水平到达角，/>表示为垂直到达角；

所述时变坐标S_l,n(t)的表达式为：

S_l,n(t)＝S_l,n(t-Δt)-v_RΔt

公式中，v_R表示实际速度，Δt指一段极小的时间，在这段时刻假定速度不变，散射簇固定不变；

步骤S204、首先根据所述接收端以及发送端的相对速度，计算其对应的移动场景参数，并且设定生成概率和消亡概率，然后通过马尔可夫生灭过程模拟簇的演进过程，并且计算得到新生成簇的数学期望，最后将经过马尔可夫生灭过程模拟后，新生成的簇代入到步骤S201中且继续执行S201-步骤S204，计算相应的时延参数和相应的近地散射点坐标以及该近地散射点对应的移动散射点的时变坐标；其中，

所述新生成簇的数学期望的表达式为：

公式中，D_c为基于不同场景的相关因子，χ_p(t,Δt)为所述的移动场景参数，λ_R为所述消亡概率，λ_G为所述生成概率；

所述移动场景参数的的表达式为：

χ_p(t,Δt)＝χ_T(t,Δt)+χ_R(t,Δt)

公式中，v'_R(t),v'_T(t)表示转换后的相对速度。

当通过马尔可夫生灭过程完成t到t+Δt时间段的簇的演进过程后，再将新生成的簇代入到步骤S201中，继续进行，即可生成相应的时延以及散射点参数。

步骤S3、根据所述接收端和发送端的坐标以及散射点的时变坐标计算出时变时延参数、时变功率参数以及时变角度参数，其中，所述时变时延参数包括视距成分时延和非视距成分时延，所述时变角度参数包括视距分量的到达角、视距分量的离开角、非视距分量的到达角、非视距分量的离开角；

所述步骤S3具体包括：

步骤S301、根据接收端以及发送端和散射点的距离计算得到视距成分时延和非视距成分时延，其中，视距成分时延由接收端和发送端的距离获得，表达式为：

τ_LOS(t)＝D_LOS(t)/c

公式中，D_LOS(t)为接收端以及发送端之间的距离，c为光速；

非视距成分时延通过选择每一个簇上任意一个散射点，计算其到接收端以及发送端的距离和求得，表达式为：

τ_l(t)＝(D_l,R(t)+D_l,T(t))/c

公式中，D_l,R/T(t)为散射点到R或者T的距离，其中，R表示为接收端，T表示为发送端；

步骤S302、根据步骤S301中求得的非视距成分时延，计算得到时变功率参数，表达式为：

公式中，τ_l(t)表示非视距成分时延，r_τ在LoS分量中为3，在NLoS分量中为2.1，σ为延迟扩展，其对数值lg(σ)服从正态分布，在LoS分量中均值为-0.24log₁₀(1+f_c)-7.14，方差为0.38，在NLoS分量中的均值为-0.24log₁₀(1+f_c)-6.83，方差为0.16log₁₀(1+f_c)+0.28；Y_l其值服从(0,ζ²)的高斯分布，ζ＝3；

为了让簇的功率和为1，所以对簇的功率进行归一化处理，表达式为：

公式中，L(t)为时变的有效散射簇数量；

步骤S303、根据不同时刻的接收端和发送端的坐标以及散射点的时变坐标，计算角度参数，其中，

视距分量的到达角的表达式为：

公式中，D_LOS(t)表示为接收端以及发送端之间的距离，为视距分量的水平到达角，/>为视距分量的垂直到达角；

视距分量的离开角的表达式为：

公式中，为视距分量的水平离开角，/>为视距分量的垂直离开角；

非视距分量的到达角的表达式为：

公式中，为第l个散射簇上的第n个散射点散射分量到接收端的水平到达角，为第l个散射簇上的第n个散射点散射分量到接收端的垂直到达角，(x_l,n,y_l,n,z_l,n)为第l个散射簇上的第n个散射点坐标，D_l,R(t)为第l个散射簇到接收端的距离；

非视距分量的离开角的表达式为：

公式中，为第l个散射簇上的第n个散射点散射分量到发射端的水平离开角，为第l个散射簇上的第n个散射点散射分量到发射端的垂直离开角，D_l,T(t)为第l个散射簇到发射端的距离。

步骤S4、根据步骤S3中的时变时延参数，时变功率参数以及时变角度参数计算信道冲激响应。

步骤S4具体步骤如下：

步骤S401、计算接收端第p根天线和发送端第q根天线间的衰落信道的冲激响应，表达式为：

公式中，P_l(t)为时变功率参数，h_l,pq(t)为第l个抽头的系数，τ_l为第l个抽头的时延参数，δ(τ-τ_l)表示时延的冲激函数，τ为时延变量，L(t)表示随着时间变化的有效散射簇数量；

步骤S402、计算视距信道和非视距信道的抽头系数，其中，

视距分量系数的表达式为：

非视距分量系数的表达式为：

公式中，l>2，λ为传输波长，r_p为接收端R的第p根天线相对接收端R天线阵列中心的相对距离，r_q为发送端T的第q根天线相对发送端T天线阵列中心的相对距离，θ_l,n为随机相移，是指每个散射射线自己随机生成的相移，是在(0，2π]上均匀分布的随机变量；Φ_LOS(t)和/>表示发射角，Ψ_LOS(t)和/>代表到达角，f_Tm表示发送端的最大多普勒频移，f_Rm表示接收端最大多普勒频移。

本发明未详述之处，均为本领域技术人员的公知技术。

以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。

Claims (1)

1.一种面向无人机空对空通信的几何随机信道建模方法，其特征在于，包括如下步骤：

步骤S1、获取无人机的接收端以及发送端在初始时刻的实际位置和实际速度，通过构建椭球模型，对所述实际速度进行变换，得到所述接收端以及发送端的相对速度；

步骤S2、使用马尔可夫生灭过程模拟由于收发端快速移动导致的簇的生灭过程，并且随机生成某一时刻的到达角，并通过到达角随机生成散射点位置坐标及其时变坐标；

步骤S3、根据所述接收端和发送端的坐标以及散射点的时变坐标计算出时变时延参数、时变功率参数以及时变角度参数，其中，所述时变时延参数包括视距成分时延和非视距成分时延，所述时变角度参数包括视距分量的到达角、视距分量的离开角、非视距分量的到达角、非视距分量的离开角；

步骤S4、根据步骤S3中的时变时延参数，时变功率参数以及时变角度参数计算信道冲激响应；

所述步骤S1具体包括：

步骤S101、获取所述接收端以及发送端在初始时刻的实际位置和实际速度；

步骤S102、根据步骤S101中获取的实际位置和实际速度，并且以所述接收端为参考点，构建椭球模型，再确定所述接收端以及发送端在模型中的坐标，其中，所述接收端在该模型的坐标系中，x轴以及y轴的值均为0，z轴的值为其离地高度；

步骤S103、根据步骤S101中获取的实际位置和实际速度，将初始时刻的接收端以及发送端进行速度转换，使得所述接收端在模型中为竖直上下变速运动，所述发送端为相对于接收端的任意方向的变速移动，再获得所述接收端以及发送端的相对速度；

所述步骤S2具体包括：

步骤S201、设定初始时刻的簇的数量，并且根据接收端和发送端之间的相对距离以及随机分布生成簇的时延参数，再通过该时延参数计算得到簇的几何参数，其中，

所述簇的时延参数的表达式为：

公式中，τ_LOS(t)＝D_LOS(t)/c，τ_LOS(t)为视距成分时延，D_LOS(t)为接收端以及发送端之间的距离，c为光速，服从指数分布；

所述簇的几何参数的表达式为：

a′＝τ_l*c/2

步骤S202、为每一个簇随机生成若干对水平到达角和垂直到达角，其中通过每一对到达角均获得一个散射点的位置坐标；

步骤S203、根据所述接收端以及发送端的坐标，计算无人机通信场景的地面有效散射区域，并根据所述地面有效散射区域以及随机生成的水平到达角和垂直到达角，计算对应的散射点的位置坐标；设定所述散射点为相对地面静止，根据接收端的实际速度将所述散射点转换为相对接收端的移动散射点，并且计算该移动散射点的时变坐标；其中，

所述地面有效散射区域的表达式为：

公式中，b'²＝a'²-f'²，(x,y)为地面有效散射区域上的任意一点坐标，(x_c,y_c)表示地面有效散射区域的中心位置坐标，a,b为地面有效散射区域建模为椭圆的长半轴和短半轴长度，θ为收发端水平面投影点的夹角，以接收端为水平面坐标原点，则发送端的坐标为(x_T,y_T)，则/>h_T表示为发送端在坐标系中z轴的值，h_R表示为接收端在坐标系中z轴的值，f'表示为接收端与发送端之间距离的一半；

所述散射点的时变位置坐标S_l,n(t)在每个任意时刻t的表达式为：

公式中，(x_l,n,y_l,n,z_l,n)为散射点的坐标；

表示为水平到达角，/>表示为垂直到达角；

所述时变坐标S_l,n(t)的表达式为：

S_l,n(t)＝S_l,n(t-Δt)-v_RΔt

公式中，v_R表示实际速度，Δt指一段极小的时间，在这段时刻假定速度不变，散射簇固定不变；

步骤S204、首先根据所述接收端以及发送端的相对速度，计算其对应的移动场景参数，并且设定生成概率和消亡概率，然后通过马尔可夫生灭过程模拟簇的演进过程，并且计算得到新生成簇的数学期望，最后将经过马尔可夫生灭过程模拟后，新生成的簇代入到步骤S201中且继续执行S201-步骤S204，计算相应的时延参数和相应的近地散射点坐标以及该近地散射点对应的移动散射点的时变坐标；其中，

所述新生成簇的数学期望的表达式为：

公式中，D_c为基于不同场景的相关因子，χ_p(t,Δt)为所述的移动场景参数，λ_R为所述消亡概率，λ_G为所述生成概率；

所述移动场景参数的的表达式为：

χ_p(t,Δt)＝χ_T(t,Δt)+χ_R(t,Δt)

公式中，v'_R(t),v'_T(t)表示转换后的相对速度；

所述步骤S3具体包括：

步骤S301、根据接收端以及发送端和散射点的距离计算得到视距成分时延和非视距成分时延，其中，视距成分时延由接收端和发送端的距离获得，表达式为：

τ_LOS(T)＝D_LOS(t)/c

公式中，D_LOS(T)为接收端以及发送端之间的距离，c为光速；

非视距成分时延通过选择每一个簇上任意一个散射点，计算其到接收端以及发送端的距离和求得，表达式为：

τ_l(t)＝(D_l,R(t)+D_l,T(t))/c

公式中，D_l,R/T(t)为散射点到R或者T的距离，其中，R表示为接收端，T表示为发送端；

步骤S302、根据步骤S301中求得的非视距成分时延，计算得到时变功率参数，表达式为：

公式中，τ_l(t)表示非视距成分时延，r_τ在LoS分量中为3，在NLoS分量中为2.1，σ为延迟扩展，其对数值lg(σ)服从正态分布，其值服从(0,ζ²)的高斯分布，ζ＝3；

为了让簇的功率和为1，所以对簇的功率进行归一化处理，表达式为：

公式中，L(t)为时变的有效散射簇数量；

步骤S303、根据不同时刻的接收端和发送端的坐标以及散射点的时变坐标，计算角度参数，其中，

视距分量的到达角的表达式为：

公式中，D_LOS(t)表示为接收端以及发送端之间的距离，为视距分量的水平到达角，为视距分量的垂直到达角；

视距分量的离开角的表达式为：

公式中，为视距分量的水平离开角，/>为视距分量的垂直离开角；

非视距分量的到达角的表达式为：

公式中，为第l个散射簇上的第n个散射点散射分量到接收端的水平到达角，/>为第l个散射簇上的第n个散射点散射分量到接收端的垂直到达角，(x_l,n,y_l,n,z_l,n)为第l个散射簇上的第n个散射点坐标，D_l,R(t)为第l个散射簇到接收端的距离；

非视距分量的离开角的表达式为：

公式中，为第l个散射簇上的第n个散射点散射分量到发射端的水平离开角，/>为第l个散射簇上的第n个散射点散射分量到发射端的垂直离开角，D_l,T(t)为第l个散射簇到发射端的距离；

步骤S4具体步骤如下：

步骤S401、计算接收端第p根天线和发送端第q根天线间的衰落信道的冲激响应，表达式为：

公式中，P_l(t)为时变功率参数，h_l,pq(t)为第l个抽头的系数，τ_l为第l个抽头的时延参数，δ(τ-τ_l)表示时延的冲激函数，τ为时延变量，L(t)表示随着时间变化的有效散射簇数量；

步骤S402：计算视距信道和非视距信道的抽头系数，其中，

视距分量系数的表达式为：

非视距分量系数的表达式为：

公式中，l>2，λ为传输波长，r_p为接收端R的第p根天线相对接收端R天线阵列中心的相对距离，r_q为发送端T的第q根天线相对发送端T天线阵列中心的相对距离，θ_l,n为随机相移，是指每个散射射线自己随机生成的相移，是在(0，2π]上均匀分布的随机变量；Φ_LOS(t)和/>表示发射角，Ψ_LOS(t)和/>代表到达角，f_Tm表示发送端的最大多普勒频移，f_Rm表示接收端最大多普勒频移。