CN107171703B - 一种能够模拟室内多天线通信系统中衰落信号传播特性的方法 - Google Patents

一种能够模拟室内多天线通信系统中衰落信号传播特性的方法 Download PDF

Info

Publication number
CN107171703B
CN107171703B CN201710573679.7A CN201710573679A CN107171703B CN 107171703 B CN107171703 B CN 107171703B CN 201710573679 A CN201710573679 A CN 201710573679A CN 107171703 B CN107171703 B CN 107171703B
Authority
CN
China
Prior art keywords
cluster
angle
sub
fading
path
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
Application number
CN201710573679.7A
Other languages
English (en)
Other versions
CN107171703A (zh
Inventor
李岳衡
胡明亮
付婉
杨茂辉
汪成
居美艳
黄平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201710573679.7A priority Critical patent/CN107171703B/zh
Publication of CN107171703A publication Critical patent/CN107171703A/zh
Application granted granted Critical
Publication of CN107171703B publication Critical patent/CN107171703B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3911Fading models or fading generators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3912Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Electromagnetism (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种能够模拟室内多天线通信系统中衰落信号传播特性的方法,包括以下几个步骤:生成发射端第k簇的发射角;生成第k簇发射角的偏移角,然后通过偏移角与发射角相加得到发射端第k簇第l子径的发射角;使用同样的方法生成接收端第k簇第l子径的到达角;生成相位角;生成衰落指数m的Nakagami‑m随机数,Nakagami‑m随机数与相位角相乘组成复包络随机数;利用相位差异表达式与导向向量表达式生成发射端导向向量和接收端导向向量;将发射端导向矢量、接收端导向向量和复包络随机数代入信道冲击响应公式得到M*N信道矩阵。本发明能够准确地描述室内复杂传输环境中、尤其是多天线通信系统信号衰落的传播特性。

Description

一种能够模拟室内多天线通信系统中衰落信号传播特性的 方法
技术领域
本发明涉及一种能够模拟室内多天线通信系统中衰落信号传播特性的方法,适用于综合考虑阴影效应及小尺度衰落等复合衰落的复杂室内多天线无线通信场景。
背景技术
无线信号在传播过程中,由于受传播环境复杂性的影响,电磁波在遇到或通过障碍物时容易产生各种反射、折射、衍射等物理现象,从而导致接收端接收到的电信号实际上是幅度和相位随传输的不同路径和时间变化而发生随机改变的衰落信号。根据通过无线信道传输后信号的特点,无线衰落可以分为大尺度衰落(Large-scale Fading)和小尺度衰落(Small-scall Fading)。对于本发明所涉及的应用环境,如大型办公楼、大型超市、商场等复杂室内环境,总体而言,信号传输的距离较短,信号经历的衰落通常为小尺度衰落(小尺度衰落是指当移动台在一个较小的范围内移动,接收信号在短期内的快速波动,反映的是较短的距离和时间之内接收信号所呈现的快速起伏特性)。
小尺度衰落通常包括Rayleigh衰落、Rice衰落和Nakagami衰落:当发射端与接收端之间不存在直射路径,接收信号仅是通过周围障碍物的衍射、散射等途径到达接收端,此时的信号即包络服从Rayleigh分布;而当接收信号中存在直射信号分量时,此时的信号包络服从Rice分布;通过大量的数据分析,对于实际复杂的传输环境,使用Nakagami衰落来描述小尺度衰落更为准确和便利,这是因为Nakagami衰落可以通过调整衰落指数m的值,模拟包括单边高斯、Rayleigh、Rice、近似高斯衰落等,非常符合复杂通信环境。
近年来通过对室内通信系统的研究还发现,接收机接收到的信号包含非常多的射线簇(cluster),每一簇(cluster)又含有很多子径,每簇和每簇中的子径的到达时间基本上都服从独立泊松过程,从而构成丰富的多径信道传播环境。
在已有的传统描述室内衰落信道的S-V模型中,使用Rayleigh分布来描述簇中子径的幅度的方法不能很好地模拟复杂室内通信环境,故而存在一定的局限性。
发明内容
针对现有技术存在的不足,本发明目的是提供一种能够模拟室内多天线通信系统中衰落信号传播特性的方法,该方法综合考虑阴影效和小尺度衰落应给接收信号带来的不利影响,能够准确地描述室内复杂传输环境中、尤其是多天线通信系统信号衰落的传播特性。
为了实现上述目的,本发明是通过如下的技术方案来实现:
本发明的一种能够模拟室内多天线通信系统中衰落信号传播特性的方法,包括以下几个步骤:
(1)生成在(0,2π)范围内随机分布的发射端第k簇的发射角angle_T;
(2)利用公式(9)生成服从拉普拉斯变换的相对于第k簇发射角angle_T的偏移角w_kl_T,然后通过偏移角w_kl_T与发射角angle_T相加得到发射端第k簇第l子径的发射角theta_T;
(3)重复步骤(1)与步骤(2),使用同样的方法生成接收端第k簇第l子径的到达角theta_R;
(4)生成在(0,2π)范围内均匀分布的相位角;
(5)利用近似逆变换法生成衰落指数为m的Nakagami-m随机数,Nakagami-m随机数与相位角相乘组成复包络随机数;
(6)利用相位差异表达式与导向向量表达式即公式(12)生成发射端导向向量W_T_m和接收端导向向量W_R_n;
(7)将所述发射端导向矢量W_T_m、接收端导向向量W_R_n和复包络随机数代入信道冲击响应公式得到M*N信道矩阵。(通过该信道矩阵可以用现有的软件仿真得到信道冲激响应图以及信道容量PDF图,从而进一步验证该方法的特性)
本发明的步骤(4)、(5)中,记第k簇的第一个子径到达时间为Tk、平均到达率为A的泊松过程,每簇中的子径到达时间是服从平均到达率为a的泊松过程,则接收端信号簇与簇中子径到达时间分布分别服从以下指数分布:
Figure BDA0001350345180000021
Figure BDA0001350345180000022
τl,k表示第k簇中第l子径的到达时间,τ0,k是第k簇的第一条子径的到达时间,并设定为第k簇的到达时间,即τ0,k=Tk
设βl,k和θl,k分别指的是第k簇第l子径的幅度和相位,则K簇及每簇L条子径的时间和空间上的信道冲激响应表达式为:
Figure BDA0001350345180000023
其中,K为簇的数量,L为每簇中的子径数量,θT为发射机的发射角,θR为接收机的到达角,
Figure BDA0001350345180000024
为第k簇的平均发射角,
Figure BDA0001350345180000025
为第k簇的平均到达角,
Figure BDA0001350345180000026
为第k簇第l子径相对于
Figure BDA0001350345180000027
的发射角,
Figure BDA0001350345180000028
为第k簇第l子径相对于
Figure BDA0001350345180000029
的到达角,τk为第k簇的到达时间,τl,k为第k簇第l子径相对于τk的到达时间,βl,k为第k簇第l子径的复增益系数,
Figure BDA0001350345180000031
为第k簇第l子径的相位θl,k的复数表达形式,δ为一物理字符,如d(t-t0)表示时间在t0时刻有值,Xσ和βl,k分别代表服从对数正态分布的阴影效应和服从Nakagami分布的小尺度信道衰落接收变量值,Nakagami分布的概率密度函数表达式为:
Figure BDA0001350345180000032
对于Nakagami分布而言,上式(4)中
Figure BDA0001350345180000033
为Gamma函数,m称为衰落因子或衰落指数,表示小尺度衰落的严重程度,ω指第k簇第l子径的平均功率,表达式分别为:
Figure BDA0001350345180000034
Figure BDA0001350345180000035
Γ和γ分别是簇和子径中功率衰减的时间变量,
Figure BDA0001350345180000036
是第一簇第一子径的平均功率,E(·)表示求期望,在式(5)中,簇平均功率以指数
Figure BDA0001350345180000037
衰减,簇中子径平均功率以指数
Figure BDA0001350345180000038
衰减,当第一簇第一子径的平均功率
Figure BDA0001350345180000039
确定,簇中其他子径的平均功率通过上述公式得到,通常情况下不考虑路径损耗,
Figure BDA00013503451800000310
被归一化为1。
步骤(1)、(2)、(3)中,当阴影衰落由阻挡衰减决定时,阴影衰落用下面的模型来分析,其衰减表达式为:
s(d)=e-αd (7)
式中,d是障碍物物体的厚度,α是表征传输信道中各种障碍物的综合衰减因子;若第i个障碍物的衰减常数是αi,宽度的随机值di,那么信号经历的衰减表达式服从公式:
Figure BDA00013503451800000311
如果信号的传输路径中存在多个的障碍物,根据中心极限定理,∑iαidi可看作是服从高斯分布的随机序列;如此,logs(di)就是一个均值为μ、方差为σ的高斯随机变量;因此在公式(3)中引入对数正态分布
Figure BDA00013503451800000312
以反映阴影衰落的影响,σx表示随机变量X的方差;
对于发射端与接收端结构相同的系统,AOA与AOD的分布是相同的;根据室内测试数据以及发射端与接收端的对称性,AOA/AOD服从双边拉普拉斯分布:
Figure BDA0001350345180000041
式中,ω是相对于簇的平均值的发射角和到达角,σ(R,T)是以弧度表示的角度标准差,σP表示到达角或发射角的标准差。
步骤(6)中,由不同天线上接收或发射信号波程差导致的相位差异表达式可通过坐标轴旋转方法计算得到:
Figure BDA0001350345180000042
式中,k0=2π/λ0是自由空间波数,λ0是中心频率对应的波长,
Figure BDA0001350345180000043
Figure BDA0001350345180000044
分别是第m个接收天线与第n个发射天线的坐标,θ是接收天线或发射天线的相对到达角或发射角。
步骤(6)中,所述导向向量表达式为:
Figure BDA0001350345180000045
其中,
Figure BDA0001350345180000046
Figure BDA0001350345180000047
代表的是第m根接收天线和第n根发射天线的增益方向图,j为在复数中虚数的标志,如复数a+jb,则a是实部,b是虚部。
步骤(7)中,对于上述任意一次信道实现,假设天线均为全向辐射,在不考虑天线之间耦合效应的情形下,所述信道冲击响应公式如下:
Figure BDA0001350345180000048
其中,t表示信号到达时间。
衰落指数m分别取0.65、1、4。
本发明以Nakagami分布来表征室内接收信号射线簇中各子径幅度的分布,通过调整Nakagami衰落指数m的值,进而模拟不同衰落程度下的室内衰落环境;同时采用对数正态分布来近似大型障碍物所产生的阴影效应给信号带来的影响,并通过引入无线电波与天线之间的发射角(Angle of Departure,AOD)与到达角(Angle of Arrive,AOA)、以及天线辐射方向图和信号相程差等的概念来综合表征多天线传输的效果,进而得到描述该复杂室内衰落特性的信道冲激响应数学表达式及实现过程。
附图说明
图1(a)为发射端第k簇收发信号信道传输参数图;
图1(b)为接收端第k簇收发信号信道传输参数图;
图2为本发明中所设计的各簇子径复包络基于Nakagami分布时的信道容量PDF图;
图3(a)为本发明中所设计的各簇子径复包络基于Nakagami分布时的信道冲激响应图(衰落指数m=0.65);
图3(b)为本发明中所设计的各簇子径复包络基于Nakagami分布时的信道冲激响应图(衰落指数m=1);
图3(c)为本发明中所设计的各簇子径复包络基于Nakagami分布时的信道冲激响应图(衰落指数m=4);
图4为本发明所提出的模拟室内复杂多天线衰落信号传播特性方法的具体实现流程图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
本发明的实施例是在室内复杂多天线传输环境中,构造能够最大程度上逼近实际室内应用环境信道特性的一种有效方法,以为通信系统的软硬件设计、性能仿真与评估、以及后续性能的优化提供服务。
在室内复杂多天线通信系统中,通过大量的实验测试数据表明,接收机接收到的信号包往往括很多簇(cluster),每一簇(cluster)中包含一组子径,每簇和每簇中子径的到达时间都服从独立泊松过程,每一路径的时延可以是任意值大小。记第k簇的第一个子径是到达时间为Tk、平均到达率为A的泊松过程,每簇中的子径到达时间是服从平均到达率为a的泊松过程,则接收端信号簇与簇中子径到达时间分布分别服从以下指数分布:
Figure BDA0001350345180000051
Figure BDA0001350345180000052
这里τl,k表示第k簇中第l子径的到达时间,τ0,k是第k簇的第一条子径的到达时间,并设定为第k簇的到达时间(即τ0,k=Tk)。
在多天线通信系统中,当我们从时间和空间两个维度来考虑接收信号时,必然需要引入无线电波与天线之间的发射角AOD和到达角AOA这两个参量,单簇中(比如第k簇)的子径传输信号的发射参数和接收参数如图1(a)和图1(b)所示。设βl,k和θl,k分别指的是第k簇第l子径的幅度和相位,则K簇及每簇L条子径的时间和空间上的信道冲激响应表达式可表达为:
Figure BDA0001350345180000061
上式中各参数如表1所示;其中,Xσ和βl,k分别代表服从对数正态分布的阴影效应(后面会进一步解释)和服从Nakagami分布的小尺度信道衰落接收变量值,Nakagami分布的概率密度函数表达式为:
Figure BDA0001350345180000062
表1信道表达式中涉及的参数
Figure BDA0001350345180000063
对于Nakagami分布而言,上式(4)中
Figure BDA0001350345180000064
为Gamma函数,m和ω是Nakagami分布的两个重要参数,m称为衰落因子或衰落指数,表示小尺度衰落的严重程度,ω指第k簇第l子径的平均功率,表达式分别为:
Figure BDA0001350345180000065
Figure BDA0001350345180000066
这里Γ和γ分别是簇和子径中功率衰减的时间变量,
Figure BDA0001350345180000067
是第一簇第一子径的平均功率。在式(5)中,簇平均功率以指数
Figure BDA0001350345180000068
衰减,簇中子径平均功率以指数
Figure BDA0001350345180000069
衰减。一旦第一簇第一子径的平均功率
Figure BDA00013503451800000610
确定,簇中其他子径的平均功率通过上述公式就都能够得到,这里
Figure BDA00013503451800000611
由指定场景的路径决定,通常情况下不考虑路径损耗,
Figure BDA0001350345180000071
被归一化为1。
当阴影衰落主要由阻挡衰减决定时,阴影衰落可以用下面简单的模型来分析,其衰减表达式为:
s(d)=e-αd (19)
式中,d是障碍物物体的厚度,α是表征传输信道中各种障碍物的综合衰减因子。若第i个障碍物的衰减常数是αi,宽度的随机值di,那么信号经历的衰减表达式服从公式:
Figure BDA0001350345180000072
如果信号的传输路径中存在较多的障碍物,根据中心极限定理,∑iαidi可看作是服从高斯分布的随机序列。如此,logs(di)就是一个均值为μ、方差为σ的高斯随机变量。因此可以在公式(3)中引入对数正态分布
Figure BDA0001350345180000073
以反映阴影衰落的影响。
对于发射端与接收端结构相同的系统,AOA与AOD的分布应该是相同的。根据室内测试数据以及发射端与接收端的对称性,AOA/AOD服从双边拉普拉斯分布:
Figure BDA0001350345180000074
式中,ω是相对于簇的平均值的发射角和到达角,σ(R,T)是以弧度表示的角度标准差。对于上述任意一次信道实现,假设天线均为全向辐射,在不考虑天线之间耦合效应的情形下,信道冲击响应可写成:
Figure BDA0001350345180000075
其中导向向量的表达式为:
Figure BDA0001350345180000076
上述导向向量表达式中的
Figure BDA0001350345180000077
Figure BDA0001350345180000078
代表的是第m根接收天线和第n根发射天线的增益方向图,而由不同天线上接收或发射信号波程差导致的相位差异表达式可通过坐标轴旋转方法计算如下:
Figure BDA0001350345180000081
式中,k0=2π/λ0是自由空间波数,λ0是中心频率对应的波长,
Figure BDA0001350345180000082
Figure BDA0001350345180000083
分别是第m个接收天线与第n个发射天线的坐标。
图2、图3(a)、图3(b)、图3(c)分别给出了簇内子径包络服从Nakagami衰落时的室内多天线传输信道特性的计算机模拟效果图。在模拟过程中,主要信道参数取值如下:平均簇到达率A=0.023,平均子径到达率a=2.5,簇平均功率衰减指数Γ=7.4,子径平均功率衰减指数γ=4.3,对数正态标准差σX=3dB,簇的数量K=10,簇中子径的数量L=15,发射天线个数M=2,接收天线个数N=2。图2为本发明中所设计的各簇子径复包络基于Nakagami分布时的信道容量PDF图,如图所示,当衰落指数m分别取值为0.65和1时,PDF最大值所对应的信道容量值无明显差别,然而衰落指数m取值为4时,PDF最大值所对应的信道容量值远大于衰落指数m取值为0.65和1时的信道容量值,这是因为随着衰落指数m取值变大,信道衰落程度减小的缘故。如图3(a)、图3(b)、图3(c)所示,当衰落指数m较大时(比如m=4),相同时延情况下簇中子径的幅度值比m较小时(比如m=0.65)簇中子径的幅度值要大得多,这也与Nakagami衰落过程中,信号经历的衰落程度随衰落指数m的增大而减小这一基本物理特性相符合。
参见图4,以下为本发明所提出的模拟室内复杂多天线衰落信号传播特性方法的具体实现步骤及流程:
(1)生成在(0,2π)范围内随机分布的发射端第k簇的发射角angle_T;
(2)生成服从拉普拉斯变换的相对于第k簇发射角angle_T的偏移角w_kl_T,然后通过偏移角w_kl_T与发射角angle_T相加得到发射端第k簇第l子径的发射角theta_T;
(3)重复步骤1与步骤2,可以使用同样的方法生成接收端第k簇第l子径的到达角theta_R;
(4)生成在(0,2π)范围内均匀分布的相位角;
(5)利用近似逆变换法生成衰落指数m分别取0.65、1、4时的Nakagami-m随机数,与相位角相乘组成复包络随机数;
(6)不失一般性,发射天线与接收天线均为全向辐射,不考虑天线之间的耦合效应,利用相位函数表达式(12)与表达式(11)生成发射端导向向量W_T_m和接收端导向向量W_R_n;
(7)将发射端导向矢量、接收端导向向量和复包络随机数代入公式(10)得到M*N信道矩阵;
以下为本发明所提出的模拟室内复杂多天线衰落信号传播特性方法的具体实现算法:
Figure BDA0001350345180000091
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (7)

1.一种能够模拟室内多天线通信系统中衰落信号传播特性的方法,其特征在于,包括以下几个步骤:
(1)生成在(0,2π)范围内随机分布的发射端第k簇的发射角angle_T;
(2)生成服从拉普拉斯变换的相对于第k簇发射角angle_T的偏移角w_kl_T,然后通过偏移角w_kl_T与发射角angle_T相加得到发射端第k簇第l子径的发射角theta_T;
(3)重复步骤(1)与步骤(2),使用同样的方法生成接收端第k簇第l子径的到达角theta_R;
(4)生成在(0,2π)范围内均匀分布的相位角;
(5)利用近似逆变换法生成衰落指数为m的Nakagami-m随机数,Nakagami-m随机数与相位角相乘组成复包络随机数;
(6)利用相位差异表达式与导向向量表达式生成发射端导向向量W_T_m和接收端导向向量W_R_n;
(7)将所述发射端导向向量W_T_m、接收端导向向量W_R_n和复包络随机数代入信道冲击响应公式得到M*N信道矩阵。
2.根据权利要求1所述的能够模拟室内多天线通信系统中衰落信号传播特性的方法,其特征在于,步骤(4)、(5)中,记第k簇的第一个子径到达时间为Tk、平均到达率为A的泊松过程,每簇中的子径到达时间是服从平均到达率为a的泊松过程,则接收端信号簇与簇中子径到达时间分布分别服从以下指数分布:
fTk(Tk|Tk-1)=Aexp[-A(Tk-Tk-1)],k=1,2,... (1)
Figure FDA0002402292120000011
τl,k表示第k簇中第l子径的到达时间,τ0,k是第k簇的第一条子径的到达时间,并设定为第k簇的到达时间,即τ0,k=Tk
设βl,k和θl,k分别指的是第k簇第l子径的幅度和相位,则K簇及每簇L条子径的时间和空间上的信道冲激响应表达式为:
Figure FDA0002402292120000012
其中,K为簇的数量,L为每簇中的子径数量,θT为发射机的发射角,θR为接收机的到达角,
Figure FDA0002402292120000013
为第k簇的平均发射角,
Figure FDA0002402292120000014
为第k簇的平均到达角,
Figure FDA0002402292120000015
为第k簇第l子径相对于
Figure FDA0002402292120000016
的发射角,
Figure FDA0002402292120000017
为第k簇第l子径相对于
Figure FDA0002402292120000018
的到达角,τk为第k簇的到达时间,τl,k为第k簇第l子径相对于τk的到达时间,βl,k为第k簇第l子径的复增益系数,
Figure FDA00024022921200000213
为第k簇第l子径的相位θl,k的复数表达形式,δ为一物理字符,Xσ和βl,k分别代表服从对数正态分布的阴影效应和服从Nakagami分布的小尺度信道衰落接收变量值,Nakagami分布的概率密度函数表达式为:
Figure FDA0002402292120000021
对于Nakagami分布而言,上式(4)中
Figure FDA0002402292120000022
为Gamma函数,m称为衰落因子或衰落指数,表示小尺度衰落的严重程度,ω指第k簇第l子径的平均功率,表达式分别为:
Figure FDA0002402292120000023
Figure FDA0002402292120000024
Γ和γ分别是簇和子径中功率衰减的时间变量,
Figure FDA0002402292120000025
是第一簇第一子径的平均功率,E(·)表示求期望,在式(5)中,簇平均功率以指数
Figure FDA0002402292120000026
衰减,簇中子径平均功率以指数
Figure FDA0002402292120000027
衰减,当第一簇第一子径的平均功率
Figure FDA0002402292120000028
确定,簇中其他子径的平均功率通过公式(5)、(6)得到,不考虑路径损耗,
Figure FDA0002402292120000029
被归一化为1。
3.根据权利要求2所述的能够模拟室内多天线通信系统中衰落信号传播特性的方法,其特征在于,步骤(1)、(2)、(3)中,当阴影衰落由阻挡衰减决定时,阴影衰落用下面的模型来分析,其衰减表达式为:
s(d)=e-αd (7)
式中,d是障碍物物体的厚度,α是表征传输信道中各种障碍物的综合衰减因子;若第i个障碍物的衰减常数是αi,宽度的随机值di,那么信号经历的衰减表达式服从公式:
Figure FDA00024022921200000210
如果信号的传输路径中存在多个的障碍物,根据中心极限定理,∑iαidi可看作是服从高斯分布的随机序列;如此,logs(di)就是一个均值为μ、方差为σ的高斯随机变量;因此在公式(3)中引入对数正态分布
Figure FDA00024022921200000211
以反映阴影衰落的影响,
Figure FDA00024022921200000212
表示随机变量X的方差;
对于发射端与接收端结构相同的系统,AOA与AOD的分布是相同的;根据室内测试数据以及发射端与接收端的对称性,AOA/AOD服从双边拉普拉斯分布:
Figure FDA0002402292120000031
式中,ω是相对于簇的平均值的发射角和到达角,σ(R,T)是以弧度表示的角度标准差,σP表示到达角或发射角的标准差。
4.根据权利要求1所述的能够模拟室内多天线通信系统中衰落信号传播特性的方法,其特征在于,步骤(6)中,由不同天线上接收或发射信号波程差导致的相位差异表达式可通过坐标轴旋转方法计算得到:
Figure FDA0002402292120000032
式中,k0=2π/λ0是自由空间波数,λ0是中心频率对应的波长,
Figure FDA0002402292120000033
Figure FDA0002402292120000034
分别是第m个接收天线与第n个发射天线的坐标,θ是接收天线或发射天线的相对到达角或发射角。
5.根据权利要求4所述的能够模拟室内多天线通信系统中衰落信号传播特性的方法,其特征在于,步骤(6)中,所述导向向量表达式为:
Figure FDA0002402292120000035
其中,
Figure FDA0002402292120000036
Figure FDA0002402292120000037
代表的是第m根接收天线和第n根发射天线的增益方向图,j为在复数中虚数的标志。
6.根据权利要求5所述的能够模拟室内多天线通信系统中衰落信号传播特性的方法,其特征在于,步骤(7)中,对于任意一次信道实现,假设天线均为全向辐射,在不考虑天线之间耦合效应的情形下,所述信道冲击响应公式如下:
Figure FDA0002402292120000038
其中,t表示信号到达时间。
7.根据权利要求1所述的能够模拟室内多天线通信系统中衰落信号传播特性的方法,其特征在于,衰落指数m分别取0.65、1、4。
CN201710573679.7A 2017-07-14 2017-07-14 一种能够模拟室内多天线通信系统中衰落信号传播特性的方法 Expired - Fee Related CN107171703B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710573679.7A CN107171703B (zh) 2017-07-14 2017-07-14 一种能够模拟室内多天线通信系统中衰落信号传播特性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710573679.7A CN107171703B (zh) 2017-07-14 2017-07-14 一种能够模拟室内多天线通信系统中衰落信号传播特性的方法

Publications (2)

Publication Number Publication Date
CN107171703A CN107171703A (zh) 2017-09-15
CN107171703B true CN107171703B (zh) 2020-05-22

Family

ID=59824069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710573679.7A Expired - Fee Related CN107171703B (zh) 2017-07-14 2017-07-14 一种能够模拟室内多天线通信系统中衰落信号传播特性的方法

Country Status (1)

Country Link
CN (1) CN107171703B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109309536B (zh) * 2018-10-23 2020-07-17 河海大学 一种可降低Nakagami逆CDF函数近似复杂度的方法
CN109714120B (zh) * 2018-12-19 2020-05-22 河海大学 一种模拟耦合多天线室内空间衰落信道传播特性的方法
CN109861776A (zh) * 2019-01-23 2019-06-07 河海大学 一种模拟多天线室外空间耦合衰落传播特性的方法
CN109660308B (zh) * 2019-01-30 2020-09-04 江南大学 一种多墙嵌入损耗模型的建立方法及其应用
CN112311489A (zh) * 2019-08-01 2021-02-02 中移(苏州)软件技术有限公司 一种信道传输矩阵的确定方法、装置以及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0940931A2 (de) * 1998-03-03 1999-09-08 Deutsches Zentrum für Luft- und Raumfahrt e.V Verfahren zur Emulation von Mehrwegeschwund(Multipath-Fading)-Funkkanälen
CN102664669A (zh) * 2012-04-26 2012-09-12 厦门大学 一种提高室内分布式多输入多输出系统信道容量的方法
CN104660528A (zh) * 2013-11-25 2015-05-27 上海益尚信息科技有限公司 新型基于pn序列的ofdm系统导频信道估计方法及装置
CN106027183A (zh) * 2016-05-30 2016-10-12 河海大学 一种可快速实现复合衰落信道累积分布性能评估的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0940931A2 (de) * 1998-03-03 1999-09-08 Deutsches Zentrum für Luft- und Raumfahrt e.V Verfahren zur Emulation von Mehrwegeschwund(Multipath-Fading)-Funkkanälen
CN102664669A (zh) * 2012-04-26 2012-09-12 厦门大学 一种提高室内分布式多输入多输出系统信道容量的方法
CN104660528A (zh) * 2013-11-25 2015-05-27 上海益尚信息科技有限公司 新型基于pn序列的ofdm系统导频信道估计方法及装置
CN106027183A (zh) * 2016-05-30 2016-10-12 河海大学 一种可快速实现复合衰落信道累积分布性能评估的方法

Also Published As

Publication number Publication date
CN107171703A (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
CN107171703B (zh) 一种能够模拟室内多天线通信系统中衰落信号传播特性的方法
Tian et al. Semi-deterministic radio channel modeling based on graph theory and ray-tracing
CN108365903B (zh) 一种基于随机散射簇的三维Massive MIMO信道建模方法
Jiang et al. Physics-based 3D end-to-end modeling for double-RIS assisted non-stationary UAV-to-ground communication channels
CN109714120B (zh) 一种模拟耦合多天线室内空间衰落信道传播特性的方法
US9763216B2 (en) Radiator localization
Jiang et al. Three-dimensional geometry-based stochastic channel modeling for intelligent reflecting surface-assisted UAV MIMO communications
Cheng et al. Modeling and simulation for UAV air-to-ground mmWave channels
JP2014512725A (ja) オーバー・ザ・エア試験
CN104994517B (zh) 一种高空平台mimo通信系统三维信道建模方法
WO2021012833A1 (zh) 基于波束功率分配的雷达通信一体化协同探测方法及装置
JP2010226713A (ja) セルラー無線ネットワーク内におけるマルチアンテナの配備によって得られる容量および有効範囲の利得を近似し最適化する方法
Huang et al. A novel 3D GBSM for mmWave MIMO channels
Mehrnia et al. Novel maritime channel models for millimeter radiowaves
CN110212955B (zh) 一种基于射线的3d mimo信道建模的方法
Weng et al. Indoor Massive MIMO Channel Modelling Using Ray‐Launching Simulation
CN108957387A (zh) 一种卫星信号二维到达角估计方法及系统
CN108736995A (zh) 一种毫米波无线信道建模方法及装置
Zhou et al. Path loss model based on cluster at 28 GHz in the indoor and outdoor environments
CN111770527A (zh) 基于二维特征的可视与非可视信道识别方法及装置
JP5199935B2 (ja) 時空間チャネルシミュレータ
Blandino et al. Markov multi-beamtracking on 60 GHz mobile channel measurements
Jain et al. Efficient time domain HF geolocation using multiple distributed receivers
Ju et al. Statistical channel model of Wideband sub-THz radio propagation in indoor factories at 142 GHz: Towards 6G industrial wireless networks
Li et al. Sum rate analysis of MU‐MIMO with a 3D MIMO base station exploiting elevation features

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for 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

Granted publication date: 20200522

CF01 Termination of patent right due to non-payment of annual fee