CN105978647B - 一种三维双移动mimo传播信道的建模及仿真方法 - Google Patents
一种三维双移动mimo传播信道的建模及仿真方法 Download PDFInfo
- Publication number
- CN105978647B CN105978647B CN201610237783.4A CN201610237783A CN105978647B CN 105978647 B CN105978647 B CN 105978647B CN 201610237783 A CN201610237783 A CN 201610237783A CN 105978647 B CN105978647 B CN 105978647B
- Authority
- CN
- China
- Prior art keywords
- angle
- receiving
- transmitting
- dimensional
- arrival
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000004088 simulation Methods 0.000 claims abstract description 29
- 239000013598 vector Substances 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000009826 distribution Methods 0.000 claims abstract description 16
- 230000033001 locomotion Effects 0.000 claims abstract description 10
- 238000004891 communication Methods 0.000 claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000005562 fading Methods 0.000 claims description 10
- 230000010363 phase shift Effects 0.000 claims description 6
- 238000005303 weighing Methods 0.000 claims description 6
- 238000009827 uniform distribution Methods 0.000 claims description 5
- 238000005315 distribution function Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 3
- 230000001186 cumulative effect Effects 0.000 abstract 2
- 239000011159 matrix material Substances 0.000 abstract 2
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3911—Fading models or fading generators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/0082—Monitoring; Testing using service channels; using auxiliary channels
- H04B17/0087—Monitoring; Testing using service channels; using auxiliary channels using auxiliary channels or channel simulators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Mobile Radio Communication Systems (AREA)
- Radio Transmission System (AREA)
Abstract
本发明提出一种三维双移动MIMO传播信道的建模及仿真方法，供通信系统优化研究和系统性能评估使用。在三维双移动MIMO传播系统中，所述建模及仿真实现方法包括如下步骤：对三维双移动MIMO传播信道进行理论建模，获得其理论数学模型；将数学理论模型转换为便于仿真实现的累加模型；利用用户指定或实测的方位角/俯仰角联合分布特性分别计算获得三维离散到达角/离开角的取值；利用收、发端移动速度大小、方向及三维离开角/到达角计算收、发端的多普勒频移；利用收、发端天线阵列流型，获得其导引矢量矩阵；将三维离散到达角/离开角、收/发端多普勒频移、收/发端导引矢量矩阵代入三维双移动MIMO传播信道的累加形式表达式，计算获得三维双移动MIMO传播信道。
Description
技术领域：
本发明涉及无线信道的建模及仿真方法，属于无线通信领域，具体是一种综合考虑收发端双移动、方位角/俯仰角联合分布、多普勒频移、天线阵列等因素的三维双移动MIMO传播信道的建模及仿真方法。
背景技术：
多输入多输出技术(Multiple Input Multiple Output，MIMO)是当前无线移动通信领域的关键技术，它能够充分利用空间位置的多天线，不增加带宽或发送功率的前提下，有效对抗无线信道衰落，从而大幅度提高系统的频谱利用率和容量，实现高速数据传输。
建立准确而有效的MIMO信道模型是实现MIMO移动通信系统及确立传输策略的理论基础，而现有成熟的信道模型及仿真方法大多局限于二维入射形式。在实际传播环境中，由于不同高度的散射体随机分布和收发端天线自身的高度，导致电磁波经过反射、绕射和折射等方式后从不同三维方向入射到接收天线，即除了水平维的方位角外，还有垂直维的俯仰角。同时，越来越多的通信系统中的收发双方均存在移动情况，因此，研究三维双移动MIMO传播信道的模型具有重要理论意义。另一方面，为了评估MIMO通信系统在实际MIMO通信环境下的性能，也需要实时产生MIMO信道，故该信道的仿真实现具有重要工程应用价值。
发明内容：
为了准确分析及研究三维双移动MIMO传播信道的模型及特性，本发明提出一种综合考虑收发端双移动、方位角/俯仰角联合分布、多普勒频移、天线阵列等因素的三维双移动MIMO传播信道的建模及仿真方法。
本发明采用如下技术方案：一种三维双移动MIMO传播信道的建模及仿真方法，包括如下步骤：
第一步：考虑发射和接收信号三维传播的情况，建立双移动MIMO信道模型为如下形式
其中，ar(αr,βr)、at(αt,βt)分别为收、发端天线阵列的导引矢量；αr/αt为收/发端方位角，βr/βt为收/发端俯仰角；p(αt,βt)为发射信号离开角的联合概率密度函数，p(αr,βr)为接收信号到达角的联合概率密度函数；为支路附加相移，服从[0,2π]的均匀分布；fr、ft分别表示收发端移动导致的多普勒频移分量；j为虚部单位；
第二步：为了快速实现三维双移动MIMO传播信道的仿真，将第一步中的理论模型改写为如下仿真模型形式
其中，为发射端天线阵列的导引矢量，为接收端天线阵列的导引矢量；为收/发端方位角，为收/发端俯仰角；为收发端各支路的多普勒频移；N2、N1为接收端、发射端的散射支路数目；为支路附加相移，服从[0,2π]的均匀分布，j为虚部单位；
第三步：根据用户预先指定或实测的方位角/俯仰角联合分布，动态计算并产生N2、N1组到达角和离开角的方位角/俯仰角值；
第四步：利用离开角/到达角取值以及收/发端移动速度矢量计算N1N2支路的多普勒频移；
第五步：利用收/发端天线阵列的三维布局策略，获得收/发端天线阵列的导引矢量进而获得
第六步：将上述步骤中计算获得的多普勒频移天线阵列导引矢量代入三维双移动MIMO传播信道的仿真模型中，最后实时输出MIMO信道衰落。
进一步地，第三步中具体包括如下方法：
1)假设到达角或离开角的归一化二维联合分布函数为其中离开角情况令n＝n1，到达角情况令n＝n2，利用等体积法将该函数曲线围成的体积等分为N份，其中离开角情况令N＝N1-1，到达角情况令N＝N2-1，即
同时，令对的空间采样步长相等，即
2)以为中心点，在左半区间取N/4个αt/r值，并在右半区间对称取N/4个αt/r值；以为中心点，在左半区间取N/4个βt/r值，并在右半区间对称取N/4个βt/r值；
3)结合步骤1)和2)可获得(N+1)组方位角和俯仰角
进一步地，第四步中具体包括如下方法：
令为收发支路的多普勒频移，vt、vr、分别为发射端以及接收端的速度大小、方向，fc为通信载频，分别为到达角、离开角，各支路的多普勒频移包含发射端频移和接收端频移两部分因素，且可利用下式进行计算
本发明具有如下有益效果：
(1)、本发明提出的三维双移动MIMO传播信道模型，综合考虑了三维信号传播以及收发端二者都存在运动的情况；
(2)、本发明提出的三维双移动MIMO传播信道仿真方法，适用于方位角/俯仰角非独立情况，且二者的联合概率分布可服从任意分布；
(3)、本发明提出的三维双移动MIMO传播信道仿真方法，便于计算机仿真和FPGA硬件实现，同时输出信道具备时间选择性衰落、频率选择性衰落和空间选择性衰落等特征。
附图说明：
图1为基于收发端双移动散射环境下的三维MIMO信道模型。
图2为本发明所列举实施例中Von Mises Fisher二维联合概率理论分布。
图3为采用本发明产生的方位角和俯仰角的二维统计分布。
图4是本发明所列举实施例中产生的N1N2条支路的多普勒频移。
图5是本发明所列举实施例中产生的MIMO信道衰落的包络统计分布。
图6是本发明所列举实施例中产生的MIMO信道衰落的归一化自相关函数模值。
具体实施方式：
本发明三维双移动MIMO传播信道的建模及仿真方法，包括如下步骤：
第一步：考虑发射和接收信号三维传播的情况，建立双移动MIMO信道模型为如下形式
其中，ar(αr,βr)、at(αt,βt)分别为收、发端天线阵列的导引矢量；αr/αt为收/发端方位角(入射信号投影与x轴夹角)，βr/βt为收/发端俯仰角(入射信号与xoy平面夹角)；p(αt,βt)为发射信号离开角的联合概率密度函数，p(αr,βr)为接收信号到达角的联合概率密度函数；为支路附加相移，服从[0,2π]的均匀分布；fr、ft分别表示收发端移动导致的多普勒频移分量；j为虚部单位；
第二步：为了快速实现三维双移动MIMO传播信道的仿真，将第一步中的理论模型改写为如下仿真模型形式
其中，为发射端天线阵列的导引矢量，为接收端天线阵列的导引矢量；为收/发端方位角(入射信号投影与x轴夹角)，为收/发端俯仰角(入射信号与xoy平面夹角)；为收发端各支路的多普勒频移；N2、N1为接收端、发射端的散射支路数目；为支路附加相移，服从[0,2π]的均匀分布；j为虚部单位；
第三步：根据用户预先指定或实测的方位角/俯仰角联合分布，动态计算并产生N2、N1组到达角和离开角的方位角/俯仰角值；
第四步：利用离开角/到达角取值以及收/发端移动速度矢量计算N1N2支路的多普勒频移；
第五步：利用收/发端天线阵列的三维布局策略，获得收/发端天线阵列的导引矢量进而获得
第六步：将上述步骤中计算获得的多普勒频移天线阵列导引矢量代入三维双移动MIMO传播信道的仿真模型中，最后实时输出MIMO信道衰落。
其中，第三步中具体包括如下方法：
1)假设到达角或离开角的归一化二维联合分布函数为(离开角情况令n＝n1，到达角情况令n＝n2)，利用等体积法将该函数曲线围成的体积等分为N份(离开角情况令N＝N1-1，到达角情况令N＝N2-1)，即
同时，令对的空间采样步长相等，即
2)以为中心点，在左半区间取N/4个αt/r值，并在右半区间对称取N/4个αt/r值；以为中心点，在左半区间取N/4个βt/r值，并在右半区间对称取N/4个βt/r值；
3)结合步骤1)和2)可获得(N+1)组方位角和俯仰角
其中，第四步中具体包括如下方法：
令为收发支路的多普勒频移，vt、vr、分别为发射端以及接收端的速度大小、方向，fc为通信载频，分别为到达角、离开角，各支路的多普勒频移包含发射端频移和接收端频移两部分因素，且可利用下式进行计算
下面具体通过一个实施例来说明本发明三维双移动MIMO传播信道的建模及仿真方法，本实施例采用表1所示的仿真场景参数，MIMO信道的产生步骤如下所示：
表1：仿真场景参数设置
1、计算产生N+1＝41组服从Von Mises Fisher二维分布的离开角/到达角取值，根据定义及角度参数，可得
将上述分布函数代入具体实施方式部分中的方程组式(3)和(4)，求解获得的收发端方位角及俯仰角的N组取值(中心点除外)如表2所示。
表2：收发端方位角及俯仰角取值
2、将收发端的移动速度vt/vr、方向及方位角俯仰角参数代入多普勒频移表达式，计算获得各散射支路的多普勒频移如图3所示。
3、根据仿真参数，计算获得收发天线阵列的导引矢量分别为
4、将上述步骤获得的仿真参数代入仿真模型式(2)，实时计算输出MIMO信道衰落。
本实施例所获得的效果可以通过图1-图5仿真实验中所获得的具体数据作进一步的说明。我们可以看到：1)从图1、2可以看出利用本发明提出的到达角/离开角的产生方法产生出的方位角/俯仰角的联合统计分布与与理论分布吻合；2)从图4、5可以看出双移动信道的幅值统计分布、包络自相关特性与传统的单移动瑞利信道的不同。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下还可以作出若干改进，这些改进也应视为本发明的保护范围。
Claims (2)
1.一种三维双移动MIMO传播信道的建模及仿真方法，其特征在于：包括如下步骤
第一步：考虑发射和接收信号三维传播的情况，建立双移动MIMO信道模型为如下形式
其中，ar(αr,βr)、at(αt,βt)分别为收、发端天线阵列的导引矢量；αr/αt为收/发端方位角，βr/βt为收/发端俯仰角；p(αt,βt)为发射信号离开角的联合概率密度函数，p(αr,βr)为接收信号到达角的联合概率密度函数；为支路附加相移，服从[0,2π]的均匀分布；fr、ft分别表示收发端移动导致的多普勒频移分量；j为虚部单位；
第二步：为了快速实现三维双移动MIMO传播信道的仿真，将第一步中的理论模型改写为如下仿真模型形式
其中，为发射端天线阵列的导引矢量，为接收端天线阵列的导引矢量；为收/发端方位角，为收/发端俯仰角；为收发端各支路的多普勒频移；N2、N1为接收端、发射端的散射支路数目；为支路附加相移，服从[0,2π]的均匀分布；j为虚部单位；
第三步：根据用户预先指定或实测的方位角/俯仰角联合分布，动态计算并产生N2、N1组到达角和离开角的方位角/俯仰角值
1)假设到达角或离开角的归一化二维联合分布函数为其中离开角情况令n＝n1，到达角情况令n＝n2，利用等体积法将该函数曲线围成的体积等分为N份，其中离开角情况令N＝N1-1，到达角情况令N＝N2-1，即
同时，令对的空间采样步长相等，即
2)以为中心点，在左半区间取N/4个αt/r值，并在右半区间对称取N/4个αt/r值；以为中心点，在左半区间取N/4个βt/r值，并在右半区间对称取N/4个βt/r值；
3)结合步骤1)和2)可获得(N+1)组方位角和俯仰角
第四步：利用离开角/到达角取值以及收/发端移动速度矢量计算N1N2支路的多普勒频移；
第五步：利用收/发端天线阵列的三维布局策略，获得收/发端天线阵列的导引矢量进而获得
第六步：将上述步骤中计算获得的多普勒频移天线阵列导引矢量代入三维双移动MIMO传播信道的仿真模型中，最后实时输出MIMO信道衰落。
2.如权利要求1所述的三维双移动MIMO传播信道的建模及仿真方法，其特征在于：第四步中具体包括如下方法：
令为收发支路的多普勒频移，vt、vr、分别为发射端以及接收端的速度大小、方向，fc为通信载频，分别为到达角、离开角，各支路的多普勒频移包含发射端频移和接收端频移两部分因素，且可利用下式进行计算
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610237783.4A CN105978647B (zh) | 2016-04-15 | 2016-04-15 | 一种三维双移动mimo传播信道的建模及仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610237783.4A CN105978647B (zh) | 2016-04-15 | 2016-04-15 | 一种三维双移动mimo传播信道的建模及仿真方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105978647A CN105978647A (zh) | 2016-09-28 |
| CN105978647B true CN105978647B (zh) | 2018-06-22 |
Family
ID=56989537
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610237783.4A Expired - Fee Related CN105978647B (zh) | 2016-04-15 | 2016-04-15 | 一种三维双移动mimo传播信道的建模及仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105978647B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108347287B (zh) * | 2017-01-24 | 2021-04-27 | 上海大唐移动通信设备有限公司 | 一种无线环境仿真方法及装置 |
| CN107425895B (zh) * | 2017-06-21 | 2020-07-03 | 西安电子科技大学 | 一种基于实测的3d mimo统计信道建模方法 |
| CN108418645B (zh) * | 2018-01-26 | 2020-11-06 | 南京航空航天大学 | 一种非平稳移动通信信道建模及参数平滑演进方法 |
| CN108768565B (zh) * | 2018-04-17 | 2020-10-20 | 南京航空航天大学 | 一种面向车联网的动态mimo信道建模及参数计算方法 |
| CN110958069A (zh) * | 2019-12-24 | 2020-04-03 | 南京信息工程大学 | 一种基于隧道环境的三维非平稳宽带双簇信道建模方法 |
| CN111181673B (zh) * | 2020-01-08 | 2021-11-26 | 重庆邮电大学 | 一种基于双移动场景的3d无线信道建模方法 |
| CN112511199B (zh) * | 2020-11-10 | 2022-04-08 | 河海大学 | 一种低复杂度3d_mimo相关信道快速生成方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103297988A (zh) * | 2013-05-10 | 2013-09-11 | 北京邮电大学 | 适用于三维传播信道的建模方法及装置 |
| CN103747456A (zh) * | 2014-01-20 | 2014-04-23 | 南京信息工程大学 | 基于三维空间域多天线mimo统计信道的建模方法 |
| CN104994517A (zh) * | 2015-06-10 | 2015-10-21 | 哈尔滨工业大学 | 一种高空平台mimo通信系统三维信道建模方法 |
-
2016
- 2016-04-15 CN CN201610237783.4A patent/CN105978647B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103297988A (zh) * | 2013-05-10 | 2013-09-11 | 北京邮电大学 | 适用于三维传播信道的建模方法及装置 |
| CN103747456A (zh) * | 2014-01-20 | 2014-04-23 | 南京信息工程大学 | 基于三维空间域多天线mimo统计信道的建模方法 |
| CN104994517A (zh) * | 2015-06-10 | 2015-10-21 | 哈尔滨工业大学 | 一种高空平台mimo通信系统三维信道建模方法 |
Non-Patent Citations (2)
| Title |
|---|
| 一种基于MIMO系统的三维空时无线信道仿真模型;鲁凌云等;《中国铁道科学》;20060731;第27卷(第4期);第91-96段 * |
| 基于空间信道模型的MIMO系统建模与仿真;戴幻尧等;《中国电子科学研究院学报》;20130430(第2期);第125-129段 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105978647A (zh) | 2016-09-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105978647B (zh) | 一种三维双移动mimo传播信道的建模及仿真方法 | |
| CN109639378B (zh) | 矩形隧道无线传播信道建模方法 | |
| CN111245480B (zh) | 一种基于三维空间的双簇几何信道建模方法 | |
| CN103747456B (zh) | 基于三维空间域多天线mimo统计信道的建模方法 | |
| CN108418645B (zh) | 一种非平稳移动通信信道建模及参数平滑演进方法 | |
| CN109412673B (zh) | 面向无人机通信的几何随机信道模型的实时仿真方法 | |
| WO2023169590A1 (zh) | 一种适用于全频段全场景的6g普适信道建模方法 | |
| CN110417498A (zh) | 一种用于三维空间的虚拟宽带几何信道的建立方法 | |
| CN108768565B (zh) | 一种面向车联网的动态mimo信道建模及参数计算方法 | |
| CN110212955B (zh) | 一种基于射线的3d mimo信道建模的方法 | |
| CN110958069A (zh) | 一种基于隧道环境的三维非平稳宽带双簇信道建模方法 | |
| CN103018732B (zh) | 基于空时联合优化的mimo雷达波形合成方法 | |
| CN108494512B (zh) | 一种毫米波大规模mimo信道模拟系统及方法 | |
| CN114665998B (zh) | 空时一致性下的三重非平稳无线通信信道建模方法 | |
| Jiang et al. | Novel 3-D irregular-shaped geometry-based channel modeling for semi-ellipsoid vehicle-to-vehicle scattering environments | |
| CN113992290A (zh) | 一种面向轨道角动量无线通信的几何随机信道建模方法 | |
| Guan et al. | Channel sounding and ray tracing for thz channel characterization | |
| CN103984836A (zh) | 一种3d-mimo天线极化增益的联合存储方法 | |
| CN114499724A (zh) | 一种面向低轨卫星通信的空时频非平稳传输特性分析方法 | |
| Liao et al. | 3-D massive MIMO channel model for high-speed railway wireless communication | |
| Bai et al. | A 3-D wideband multi-confocal ellipsoid model for wireless MIMO communication channels | |
| CN113162712B (zh) | 一种基于传播图理论的多波束信道建模方法 | |
| CN112994814B (zh) | 一种大规模信道仿真的fpga系统 | |
| CN104219009A (zh) | 基于VonMisesFisher分布的MIMO信号衰落相关性算法 | |
| CN114696932B (zh) | 用于车联网通信的空时频非平稳无线通信信道建模方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180622 |