一种可快速实现复合衰落信道累积分布性能评估的方法
【专利摘要】本发明公开了一种可快速实现复合衰落信道累积分布性能评估的方法，包括以下几个步骤：(1)采用Gamma分布近似Lognormal分布用于模拟阴影效应，进而构建Gamma?Gamma分布用于近似原Gamma?Lognormal复合衰落信道模型；(2)在此近似模型的基础上，推导出以Meijer?G函数表示的接收信噪比累积分布函数的闭合表达式；(3)由Meijer?G函数计算出原复合衰落信道累积分布的性能。本发明所提出的Gamma?Gamma分布近似方案可以简化原精确模型的复杂无穷积分的表达，从而得出接收信噪比累积分布函数的闭合形式，进而降低CDF函数公式计算的复杂度，有利于快速分析、评估通信系统。
【专利说明】
一种可快速实现复合衰落信道累积分布性能评估的方法
技术领域
[0001] 本发明涉及一种可快速实现复合衰落信道累积分布性能评估的方法，尤其适用于 综合考虑路径损耗、阴影效应、小尺度衰落以及加性高斯白噪声的复杂移动通信系统。
【背景技术】
[0002] 无线通信系统的通信性能在很大程度上受制于系统工作的无线通信信道。在无线 通信系统中，衰落是指接收端接收到的信号幅度和相位随着传输的不同路径和时间变化而 发生随机改变的现象，这种现象将极大地影响接收机的性能。由于衰落反映的一个方面是 信号功率问题，所以根据功率减小的快慢和程度，人们将衰落分成大尺度衰落和小尺度衰 落。小尺度衰落是指电磁波在信道传播过程中受到反射、散射、绕射等因素的影响而导致接 收端所接收的信号是来自各方向的电磁波的叠加，该信号在小范围内能够引起剧烈的波动 即多径衰落;所谓的大尺度衰落是指信号在传播过程中由发射机和接收机之间的障碍物如 建筑物、高山和丛林等的阻挡造成的能量损耗和波动。在实际的无线通信系统中，天线收发 两端之间的传播路径的多样性和各种复杂的地形会导致信道的小尺度衰落具有非常强的 随机性;此外，现代无线通信系统中收发端的天线在小区内间隔较远，例如由于接入距离的 不同使得每个分布式天线端口（基站)与移动台之间的信号所经历的路径损耗不同，所以天 线端口（基站)与移动台之间的无线信号在传输过程中经历的衰落就不仅仅是小尺度衰落， 还需考虑诸如阴影衰落、路径损耗等大尺度衰落因素的不利影响。
[0003] 常见的大尺度传播模型一般包括对数距离路径损耗和对数正态阴影模型，而小尺 度衰落模型则包括Nakagami衰落、瑞利衰落和莱斯衰落等。由于在实际的无线通信环境中， Nakagami分布能很好地对包括城市在内的衰落信道中传输的信号包络进行建模，并综合了 常用的瑞利和莱斯分布等纯散射和叠加视距传输的情形，故该模型一经提出就得到业界的 广泛认可与应用。
【发明内容】
[0004] 针对现有技术存在的不足，本发明目的是提供一种可快速实现复合衰落信道累积 分布性能评估的方法，本发明所提出的Gamma-Gamma分布近似方案可以简化原精确模型的 复杂无穷积分的表达，从而得出接收信噪比累积分布函数的闭合形式，进而降低CDF函数公 式计算的复杂度，有利于快速分析、评估通信系统诸如中断概率、信道容量等性能指标。
[0005] 为了实现上述目的，本发明是通过如下的技术方案来实现：
[0006] 本发明的一种可快速实现复合衰落信道累积分布性能评估的方法，包括以下几个 步骤：
[0007] (1)设置信道参数，采用Gamma分布近似(近似指的是Gamma分布的概率密度函数 (TOF)和Lognormal分布的概率密度函数形状相似，也即两者的概率密度函数曲线能较好地 相互重合。为本领域通用术语)Lognorma 1分布用于模拟阴影效应，进而构建Gamma-Gamma分 布用于近似（该近似指的就是Gamma-Gamma分布的FOF和原Gamma-Lognormal分布的FOF相 似。为本领域通用术语)原Ga_a-Lognormal复合衰落信道模型；
[0008] (2)在复合衰落信道模型（此模型指：用来近似的Gamma-Gamma分布模型）的基础 上，推导出以Meijer-G函数表示的接收信噪比累积分布函数的闭合表达式；
[0009] ( 3)由Mei jer-G函数通过查阅公式值列表或者数值计算软件（如Matlab或 Mathematics)计算出原复合衰落信道累积分布的性能。
[0010] 步骤(1)中，所述复合衰落信道模型的构建方法如下：
[0011]在小尺度Nakagami衰落信道下，无线通信系统传输信号的包络α服从Nakagami分 布，其H)F为：
[0013] 上式中m和ω是Nakagami分布的两个重要参数，表达式分别为：
[0015] 其中，E[ ·]表示求均值，Var[ ·]表示求方差，Γ( ·)表示伽马函数，ω是衰落幅 度α的均方值，m被称为形状因子或衰落指数，表示此时小尺度衰落的严重程度，其取值满足 m 彡 1/2;
[0016] 衰落指数m的不同取值存在几种特殊情况:当m= 1/2时，其退化为单边高斯分布;当m= 1 时，正好是瑞利分布;当m>l时，Nakagami分布可以等效为莱斯因子为 的莱斯分布；
[0017] 在考虑该Nakagami衰落信道中存在加性高斯白噪声的情况下，接收端每个符号对 应的平均接收信噪比；Γ和瞬时接收信噪比γ存在如下关系：
[0019]其中，No和Es分别为高斯白噪声的功率谱密度和信号发送功率；
[0020]依据上式可知，单个符号的瞬时接收信噪比γ和接收信号包络α两者概率密度函 数之间存在如下关系：
[0022]根据随机变量PDF与其函数所得新随机变量TOF之间的雅克比行列式变换规则，可 得接收单个符号瞬时信噪比γ的roF为
[0024] 该表达式清楚地说明随机变量γ服从Gamma分布；
[0025] 如果信道中同时存在大尺度路径损耗与阴影衰落，则其平均接收信噪比^服从对 数正态分布，其H)F为
[0027]上式中，ξ = 10/1η10为一固定常数;以和〇分别为障碍物对接收信号包络功率产生 的平均路径损耗和随机波动标准差；
[0028]在综合考虑Nakagami衰落、路径损耗和阴影衰落的情况下，可得此时复合衰落信 道模型的接收信噪比γ的roF为：
[0030] 其中i指的是平均接收信噪比，同时也是该无穷积分的自变量，从上式(7)的表达 形式上可以看出，模拟实际复杂通信环境中的复合衰落信道模型其接收信噪比服从Ga_a-Lognormal 分布。
[0031] 步骤(2)中，所述接收信噪比累积分布函数的闭合表达式如下：
[0032]对(7)式得到的PDF进行积分即可得接收信噪比γ的CDF为：
[0034]用Gamma分布代替表达式(7)中的对数正态分布以对对数阴影衰落进行建模，也即 以Ga_a*布表示的平均接收信噪比的PDF为：
[0036] 上式中，η为Gamma分布的阶数;X表示平均功率；由Gamma分布得到的近似公式(9) 与原精确Lognormal分布得到的公式(6)之间，核心参数间的变换关系是：
[0038] 上式中，Φ( ·)和V (·)分别是digamma和trigamma函数;所以可得到近似后的复 合衰落信道中接收信噪比γ的FOF为：
[0040]从上式（11)的表达形式可以看出，构建的近似复合衰落信道模型其接收信噪比服 [0042]其中，K(m_n)( ·)是(m-n)阶第二类修正Bessel函数;再对上式进行积分可得接收 信噪比γ的⑶F为
[0044] 其中(4)是Mei jer-G函数且0<k<q，(XKp<q;k，1，p，q为整数；
[0045] 公式（13)给出了一种可以用工程数学上常用的表列函数（是专业术语，英文名为 "Tabulated function"，翻译过来就是"表列函数"）来表示的，复合衰落信道接收信噪比 CDF的闭合表达式通过数值仿真软件计算。
[0046] 上述数值仿真软件具体采用的是Mat lab或者Mathematic。
[0047] 本发明提出的一种快速实现复合衰落信道累积分布函数性能评估的方法，使之便 于分析系统诸如中断概率、信道容量等性能指标。该方案是在构建复合衰落模型时，仍以 Nakagami分布反映传输信号的小尺度衰落，而采用Gamma分布近似对数正态(Lognormal)分 布来描述路径损耗和阴影效应的大尺度衰落特性，得出经过近似后的接收端每个符号信噪 比的概率密度函数服从Ga_a-Gamma分布。基于此，该复合衰落信道的累积分布函数即可化 简为一个以所谓"表列函数" Meijer-G函数表示的闭合表达式，这将有利于进行计算机快速 数值仿真以分析复合衰落信道的特性。
【附图说明】
[0048] 图1为基于蜂窝小区结构的无线通信系统通用模型图；
[0049] 图2为本发明的一种可快速实现复合衰落信道累积分布性能评估的方法工作流程 图；
[0050] 图3为以衰落指数m作变量，采用Gamma-Gamma近似复合衰落信道接收信噪比累积 分布特性与精确Ga_a-lognormal累积分布函数曲线对比图；
[0051] 图4为以阴影衰落程度〇作变量，采用Gamma-Gamma近似复合衰落信道接收信噪比 累积分布特性与精确累积分布函数曲线对比图；
[0052] 图5为以平均路径损耗μ作变量，采用Gamma-Gamma近似复合衰落信道接收信噪比 累积分布特性与精确累积分布函数曲线对比图。
【具体实施方式】
[0053] 为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合
【具体实施方式】，进一步阐述本发明。
[0054]本发明的实施例可以广泛应用于存在复杂通信环境的现代无线通信系统，即通过 提出一种快速评估信道累积分布性能的方法为诸如无线通信系统中断概率及信道容量的 分析与评估提供有力依据。包含复合衰落信道的现代无线通信系统的一般模型示意图如图 1所示:在一个类似蜂窝结构的无线通信环境中，基站发射的无线信号在区域范围内进行远 距离传输时，由于信号在传输过程中会遇到诸如高大建筑、丛林等障碍物的阻挡而发生大 尺度衰落;而信号在到达接收端附近时，由于无线信号受到周围散射、反射环境等因素的影 响而产生多径小尺度衰落。所以类似手机等移动设备的接收端接收的信号就是发送端发送 的原始信号通过具有很强随机性的信道后得到的，也因此，对无线通信信道的建模研究一 直都是移动通信研究过程中的重点。
[0055]本发明在综合考虑路径损耗与阴影效应等大尺度衰落、Nakagami小尺度衰落以及 白噪声的情况下构建了一个符合实际通信传输条件的复合衰落信道模型，基于该模型得到 的接收信号信噪比的概率密度函数(Probability Density Function JDF)是一种精确的 信号功率统计特性表达式，在此基础上依据该TOF表达式得到的累积分布函数(Cumulat i ve Distribution Function，⑶F)则可以进一步评估现代通信系统中诸如中断概率、信道容量 等重要系统性能指标，故对该CDF函数的计算与评估就显得尤为重要。由于上述复合衰落信 道接收信号信噪比的TOF表达式通常是一个复杂的无穷积分形式，故无法进行化简得到闭 型的累积分布函数表达式，所以不利于进一步开展相关研究，其已成为现代通信系统性能 分析的一个难点。本发明即给出一种针对上述复合衰落信道的简化近似模型，从而能够快 速实现对复合衰落信道累积分布函数特性的性能分析和研究。
[0056]参见图2,在小尺度Nakagami衰落信道下，无线通信系统传输信号的包络α服从 Nakagami 分布，其]为：
[0058] 上式中m和ω是Nakagami分布的两个重要参数，表达式分别为：
[0060]其中，E[ ·]表示求均值，Var[ ·]表示求方差。Γ ( ·)表示伽马函数，ω是衰落幅度α 的均方值，m被称为形状因子或衰落指数，表示此时小尺度衰落的严重程度，其取值满足1/2。 衰落指数m的不同取值存在几种特殊情况:当m=l/2时，其退化为单边高斯分布;当m=l时，正好 是瑞利分布；当m>l时，Nakagami分布可以等效为莱斯因子为
.的 莱斯分布。在考虑该Nakagami衰落信道中存在加性高斯白噪声的情况下，接收端每个符号 对应的平均接收信噪比f和瞬时接收信噪比γ存在如下关系：
[0062]其中，No和Es分别为高斯白噪声的功率谱密度和信号发送功率。依据上式可知，单 个符号的瞬时接收信噪比γ和接收信号包络α两者概率密度函数之间存在如下关系：
[0064]根据随机变量PDF与其函数所得新随机变量TOF之间的雅克比行列式变换规则，可 得接收单个符号瞬时信噪比γ的roF为
[0066]该表达式清楚地说明随机变量γ服从Gamma分布。
[0067]如果信道中同时存在大尺度路径损耗与阴影衰落，则其平均接收信噪比f服从对 数正态分布，其roF为
[0069]上式中，ξ = 10/1η10，μ和〇(均以dB为单位），分别为障碍物对接收信号包络功率产 生的平均路径损耗和随机波动标准差。在综合考虑Nakagami衰落、路径损耗和阴影衰落的 情况下，可得此时复合衰落信道模型的接收信噪比γ的TOF为：
[0071]从上式(7)的表达形式上可以看出，模拟实际复杂通信环境中的复合衰落信道模 型其接收信噪比服从Gamma-Lognormal分布。再对上式得到的TOF进行积分即可得接收信噪 比γ的CDF为：
[0073]由公式(8)可知，复合衰落信道的接收信噪比γ的CDF为一复杂的无穷积分式，无 法写成工程数学中通用的闭型表达式的形式，这将非常不利于开展基于该累积分布函数表 达式而建立的现代无线通信系统中诸如中断概率、信道容量等性能指标的评估和研究工 作。所以本发明拟提供一种简化近似方法来解决此问题，即用Gamma分布代替表达式(8)中 的对数正态分布以对对数阴影衰落进行建模，也即以Gamma分布表示的平均接收信噪比的 roF为：
[0075] 上式中，η为Gamma分布的阶数;X表示平均功率。由Gamma分布得到的近似公式(9) 与原精确Lognormal分布得到的公式(6)之间，核心参数间的变换关系是：
[0077] 上式中，Φ( ·)和V (·)分别是digamma和trigamma函数。所以可得到近似后的复 合衰落信道中接收信噪比γ的FOF为：
[0079]从上式（11)的表达形式可以看出，构建的近似复合衰落信道模型其接收信噪比服 从6&1111]1&-6&1111]1&分布。令七=8/^和11 = 111/乂，经推导可得
[0081]其中，Ku-n)( ·)是(m-n)阶第二类修正Bessel函数。再对上式进行积分可得接收 信噪比γ的⑶F为
[0083] 其中，Gg.(·1')是Mei jer-G函数且0<k<q，0<Kp<q;k，1，p，q为整数。公式（13) 给出了一种可以用工程数学上常用的"表列函数"来表示的，有关复合衰落信道接收信噪比 CDF的闭合表达式，该表达式可以通过诸如数值仿真软件Matlab、Mathematic等快速、精确 计算。
[0084] 图3是以衰落指数m作变量，采用Gamma-Gamma近似得到的复合衰落信道接收信噪 比累积分布函数与原精确累积分布函数曲线对比图。从该图可以看出，在给定典型大尺度 衰落影响因子如路径损耗μ = 2(ΜΒ和阴影衰落程度〇 = 8dB的前提下，改变小尺度衰落形状 因子m所得的近似累积分布函数CDF随接收信噪比X变化的曲线和精确累积分布函数性能曲 线之间有比较好的近似效果，即Gamma-Gamma分布能较好地反映真实复合衰落信道的累积 分布特性。
[0085] 图4是以阴影衰落程度〇作变量，采用Gamma-Gamma近似复合衰落信道接收信噪比 累积分布函数与精确累积分布函数曲线对比图。在给定典型小尺度衰落影响因子如衰落指 数m=l和大尺度衰落的影响因子路径损耗μ = 2(ΜΒ的情况下，改变阴影衰落随机波动标准 差σ所得的近似和精确累积分布函数性能曲线之间的近似效果图。从该图可以看出，不论随 机波动标准差σ如何取值，近似处理后的累积分布函数曲线图与原精确模型曲线图之间的 误差较小，且两类曲线发展趋势基本一致，即该近似处理具有很高的精确性。
[0086] 图5是以平均路径损耗μ作变量，采用Gamma-Gamma近似复合衰落信道接收信噪比 累积分布函数与精确累积分布函数曲线对比图。在给定典型大尺度衰落影响因子即衰落随 机波动标准差〇 = 8dB和小尺度衰落影响因子即衰落指数m=l的情形下，改变平均路径损耗 μ所得的近似累积分布函数性能曲线亦能很好的近似精确累积分布函数性能曲线。
[0087] 由以上本发明给出的具体实施过程可以看出，该对复合衰落信道累积分布性能进 行快速评估的近似处理方法在不同的衰落影响因子条件下皆具有非常强的适用性和准确 性:也即对于决定小尺度衰落严重程度的衰落指数m和决定大尺度衰落严重程度的随机波 动标准差〇、平均路径损耗μ的典型取值，该方法均能比较准确地近似精确接收信噪比累积 分布函数曲线，进而反映复合信道的累积分布特性。
[0088] 总之，本发明所提出的Gamma-Gamma分布近似方案可以简化原精确模型的复杂无 穷积分的表达从而得出接收信噪比累积分布函数的闭合形式，进而降低CDF函数公式计算 的复杂度，有利于快速分析、评估通信系统诸如中断概率、信道容量等性能指标。
[0089]以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术 人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本 发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变 化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其 等效物界定。
【主权项】
1. 一种可快速实现复合衰落信道累积分布性能评估的方法，其特征在于，包括以下几 个步骤： (1) 设置信道参数，采用Gamma分布近似Lognormal分布用于模拟阴影效应，进而构建 Ga_a_Ga_a*布用于近似原Ga_a_Lognormal复合衰落信道模型； (2) 在复合衰落信道模型的基础上，推导出以Meijer-G函数表示的接收信噪比累积分 布函数的闭合表达式； (3) 由Meijer-G函数通过查阅公式值列表或者数值计算软件计算出原复合衰落信道累 积分布的性能。2. 根据权利要求1所述的可快速实现复合衰落信道累积分布性能评估的方法，其特征 在于，步骤(1)中，所述复合衰落信道模型的构建方法如下： 在小尺度Nakagami衰落信道下，无线通信系统传输信号的包络α服从Nakagami分布，其 I3DF 为：(1) 上式中m和ω是Nakagami分布的两个重要参数，表达式分别为：(2) 其中，Ε[ ·]表示求均值，Var[ ·]表示求方差，Γ( ·)表示伽马函数，ω是衰落幅度α的 均方值，m被称为形状因子或衰落指数，表示此时小尺度衰落的严重程度，其取值满足1/ 2；衰落指数m的不同取值存在几种特殊情况:当m= 1/2时，其退化为单边高斯分布;当m= 1时， 正好是瑞利分布;当m>l时，Nakagami分布可以等效为莱斯因子） 的莱斯分布； 在考虑该Nakagami衰落信道中存在加性高斯白噪声的情况下，接收端每个符号对应的 平均接收信噪比f和瞬时接收信噪比γ存在如下关系：(3) 其中，No和Es*别为高斯白噪声的功率谱密度和信号发送功率； 依据上式可知，单个符号的瞬时接收信噪比γ和接收信号包络α两者概率密度函数之 间存在如下关系：⑷ 根据随机变量PDF与其函数所得新随机变量PDF之间的雅克比行列式变换规则，可得接 收单个符号瞬时信噪比γ的TOF为 (D) 该表达式清楚地说明随机变量γ服从Gamma分布； 如果信道中同时存在大尺度路径损耗与阴影衰落，则其平均接收信噪比^服从对数正 态分布，其I3DF为C6) 上式中，ξ = 10/1η10为一固定常数;以和〇分别为障碍物对接收信号包络功率产生的平 均路径损耗和随机波动标准差； 在综合考虑Nakagami衰落、路径损耗和阴影衰落的情况下，可得此时复合衰落信道模 型的接收信噪比γ的F1DF为：其中，^指的是平均接收信噪比，同时也是该无穷积分的自变量，从上式(7)的表达形式 上可以看出，模拟实际复杂通信环境中的复合衰落信道模型其接收信噪比服从Gamma-Lognormal 分布。3.根据权利要求2所述的可快速实现复合衰落信道累积分布性能评估的方法，其特征 在于，步骤(2)中，所述接收信噪比累积分布函数的闭合表达式如下： 对(7)式得到的PDF进行积分即可得接收信噪比γ的CDF为：用Gamma分布代替表达式（7)中的对数正态分布以对对数阴影衰落进行建模，也即以 Ga_a*布表示的平均接收信噪比的TOF为：(9) 上式中，η为Gamma分布的阶数;X表示平均功率；由Gamma分布得到的近似公式(9)与原 精确Lognorma 1分布得到的公式(6)之间，核心参数间的变换关系是：(10) 上式中，Φ( ·)和V (·)分别是digamma和trigamma函数;所以可得到近似后的复合衰 落信道中接收信噪比γ的F1DF为： (11) 从上式（11)的表达形式可以看出，构建的近似复合衰落信道模型其接收信噪比服从 6已臟已-6已臟已分布；令七=8/乂和11=111/乂，经推导可得(12) 其中，Ku-n)( ·)是(m-n)阶第二类修正Bessel函数;再对上式进行积分可得接收信噪比 γ的CDF为其中，丨(.丨·)是Mei jer-G函数且(Xk彡q，(XKp彡q;k，l，p，q为整数； 公式（13)给出了一种可以用工程数学上常用的表列函数来表示的，复合衰落信道接收 信噪比CDF的闭合表达式通过数值仿真软件计算。4.根据权利要求3所述的可快速实现复合衰落信道累积分布性能评估的方法，其特征 在于，所述数值仿真软件具体采用的是Matlab或者Mathematic。
【文档编号】H04W24/06GK106027183SQ201610371412
【公开日】2016年10月12日
【申请日】2016年5月30日
【发明人】孙蔓, 李岳衡, 奉凤, 薛团结, 郭臣, 徐荣蓉, 孙得娣, 潘进勇, 居美艳, 黄平
【申请人】河海大学