一种基于权重法的Nakagami-m复衰落随机序列生成方法

一种基于权重法的Nakagami-m复衰落随机序列生成方法
【技术领域】
[0001] 本发明属于无线通信信道建模及模拟领域，具体来说本发明设及一种针对 化kagami-m复衰落随机序列生成方法。
【背景技术】
[0002] 无线移动通信环境中，化kagami-m衰落相较于传统的瑞利分布或莱斯分布可W更 好更广泛地表征小尺度衰落特征。化kagami-m衰落是通过变参数伽马分布的密度函数拟合 实验数据所得到的。它可W模拟从严重、适中、轻微到无衰落的信道环境，包含了瑞利分布 和莱斯分布。
[0003] 目前化kagami-m衰落的仿真主要是针对其包络分布的仿真，主要有逆变换法、舍 弃法和组合法；针对化kagami-m衰落相位分布的仿真方法较少，一般采用舍弃法。组合 法是一种简单的拟合模型，仿真速度快，但是误差较大，特别是m值不在其适用范围之内 时，仿真结果与理论值几乎完全不符。逆变换模型实用性较差精度也不高，原因在于求解 化kagami-m累积分布函数的反函数十分复杂且需要系数的拟合捜索。舍弃法精度高，但舍 弃帽子函数选择复杂，仿真速度慢，效率较低。因此，有必要研究新的高效的复衰落信道建 模方法。

【发明内容】

[0004] 为了达到上述目的，本发明提出一种基于权重法的化kagami-m复衰落分布仿真 方法，所提权重法根据概率密度函数的解析式生成服从化kagami-m衰落分布的大量随机 样本数据，然后对样本数据进行均匀随机采样，得到仿真精度所需要的目标随机序列。一种 基于权重法的化kagami-m复衰落随机序列生成方法，根据化kagami-m衰落分布函数，并采 用权重法生成幅度随机序列和相位随机序列，并使用所述幅度随机序列和相位随机序列产 生复随机序列。
[0005] 在上述技术方案的基础上，其包括W下步骤：
[0006] S1根据化kagami-m衰落幅度和相位衰落分布函数，采用权重法生成幅度随机序 列RDC和相位随机序列0怔]，其中K为目标随机序列长度；
[0007]S2由上述生成的化kagami-m幅度随机序列和相位随机序列，产生复随机序列 ZDC=R〇(]*e9w，K为目标随机序列长度。
[000引在上述技术方案的基础上，所述步骤S1中生成幅度随机序列RDC包括，
[0009] S1. 1对给定的化kagami-m幅度概率密度函数fijC)在定义域[a,b]上均匀采样得 至IJN点定义域采样序列而[闲，其中第i个定义域采样点而(1) =a+化-a/N)*i，i= 1，2，" N，a为定义域下限，b为定义域上限；
[0010] S1. 2将定义域采样序列而[闲代入幅度概率密度函数fE(.)中，得到大小为N值 域序列Ye[闲，其中第i个值域采样点YeQ] =fE狂kQ));
[0011] SI. 3计算幅度权重系数Z,[i]:权重系数定义为归一化的值域序列值乘w随机样 本数据序列大小M，并向下取整数；即 其中M〉10*K;
，
[0012]S1. 4按照定义域采样值XeQ)大小顺序，根据随机样本权重系数与[1]，依次将个 与山个X^a)加入到幅度随机样本数据序列Se[M]中，序列Se[M]即为满足幅度概率密度 分布fK(.)但无随机性的本数据；
[0013]S1.5生成所需长度K的输出随机序列RDC，其包括，首先生成[1，M]上的整数随 机变量的均匀随机分布序列Vk怔]，即Vk怔]=l+rand(l，K)*(M-l)，rancKl，K)生成1行 K列的均匀随机序列，rand函数为标准的均匀分布生成函数；其次W均匀随机分布序列值 Vc[时的排序，在样本数据序列Sc[M]中进行取样得到目标随机序列R[K]，其中所取的第j 个元素R(j) =Sk(VkU))，j= 1，2--K，此时输出的目标序列R[K]是满足服从r^Jakagami-m 幅度概率密度函数的幅度目标随机序列。
[0014] 在上述技术方案的基础上，所述步骤S1中生成相位随机序列0怔]包括，
[0015]S1. 6对给定的化kagami-m相位概率密度函数f0(.)在定义域[C，d]上均匀采 样得到N点定义域采样序列Xe[闲，其中第i个定义域采样点Xe(i) =c+(d-c/N)*i，i= 1，2…N，c为定义域下限，d为定义域上限；
[0016]S1. 7将定义域采样序列Xe[闲代入相位概率密度函数fe(.)中，得到大小为N值 域序列Ye[闲，其中第i个值域采样点Ye(i] =fe狂e(i));
[0017] SI. 8计算相位权重系数Ze[i]:权重系数定义为归一化的值域序列值乘w随机 样本数据序列大小M，并向下取整数；良F
[0018]S1.9按照定义域采样值Xe(i)大小顺序，根据随机样本权重系数Ze[i]，依次将个 Zg[i]个Xe(i)加入到相位随机样本数据序列Se[M]中，序列Se[M]即为满足相位概率密 度分布fe(.)的但无随机性的样本数据；
[0019]S1. 10生成所需长度K的输出随机序列0怔]。首先生成[1，M]上的整数随机变量 的均匀随机分布序列Ve怔]，即Ve怔]=l+rand(l，K)*(M-1)，rand(l，K)生成1行K列的 均匀随机序列，rand函数为标准的均匀分布生成函数；其次W均匀随机分布序列值Ve怔] 的排序，在样本数据序列Se[M]中进行取样得到目标随机序列0 [时，其中所取的第j个元 素0 (j) =Sg(Vg(j))，j= 1，2--K，此时输出的目标序列0 [K]是满足服从r^Jakagami-m 相位概率密度函数的相位目标随机序列。
[0020] 与现有技术相比较，本发明的有益效果在于：
[0021] (1)本发明提出了一种高效基于权重法的复衰落分布随机序列生成方法，该方法 W少量的空间复杂度换取时间效率和仿真精度；
[0022] (2)该方法的仿真精度可根据实际仿真需要，通过调整样值采样点数和目标序列 点数来实现，可避免其他方法高精度需求时插值操作带来的误差。
[0023](3)本发明中提出的权重法实际上可模拟任意具有解析表达式的分布函数的随机 序列生成，具有很大的通用性。
【附图说明】
[0024] 图1为本发明一种基于权重法的化kagami-m复衰落随机序列生成方法流程图；
[0025] 图2为本发明权重法生成任意具有解析式表达的概率密度函数分布随机序列流 程图；
[0026] 图3为本发明m=4, Q=1的化kagami-m幅度包络序列仿真结果与其他方法对 比图；
[0027] 图4为本发明Q=1的化kagami-m幅度包络序列仿真精度与其他方法对比图； [002引图5为本发明Q=1的化kagami-m幅度包络序列仿真归一化时间复杂度与其他 方法对比图；
[0029] 图6为本发明m=4, Q=1的化kagami-m相位包络序列仿真结果与其他方法对 比图；
[0030] 图7为本发明Q=1的化kagami-m相位包络序列仿真精度与其他方法对比图；
[003U 图8为本发明Q=1的化kagami-m相位包络序列仿真归一化时间复杂度与其他 方法对比图。 具体实施例
[0032] 下面结合具体实施例对本发明进行详细说明。现针对给定m=4, Q=1时的 化kagami-m衰落复衰落分布仿真，本发明过程具体实施方案如下；
[0033] 请参考图1与图2。
[0034]S1根据化kagami-m衰落幅度衰落分布函数，采用权重法生成幅度随机序列R[K] 和相位随机序列0 [K]，K= 1〇5;其中生成幅度随机序列R[K]和相位随机序列0怔]为两 并列的步骤，首先生成幅度随机序列或先生成相位随机序列均可。
[003引所述步骤S1包括，
[0036]S1. 1对给定的化kagami-m幅度概率密度函数f,(r)在定义域[a,b]上均匀采样 得到N= 1〇5点定义域采样序列Xe[闲，其中而(1) =a+化-a/N)*i，i= 1，2"，N ;
[0037]
[00測 SI. 2将定义域采样序列而[闲代入幅度概率密度函数fK(r)中，得到N=105点值 域序列Yk[闲，YeU] =fK狂K(i));
[0039] SI. 3计算幅度权重系数Z,[i]:权重系数定义为归一化的值域序列值乘W随机样 本数据序列大小M= 1〇6,并向下取整数；即 - >
[0040]S1. 4按照定义域采样值XeQ)大小 顺序，根据随机样本权重系数与[1]，依次将个 与山个X^a)加入到幅度随机样本数据序列Se[M]中，序列Se[M]即为满足幅度概率密度 分布fc(r)的随机样本数据；
[0041]S1.5生成所需长度K= 1〇5的输出随机序列RDC。首先生成[1，M]上的整数随 机变量的均匀随机分布序列Vk怔]，即Vk怔]=ramlna，M，K);其次W均匀随机分布序列值 Vc怔]的排序，在随机样本数据序列中Sc[M]进行取样得到目标随机序列RDC，即R(j)= Sk(VkO))，j= 1，2. . . .K。此时输出的目标序列RDC是满足服从Nakagami-m幅度概率密 度函数的随机序列。
[0042] 在上述技术方案的基础上，所述步骤S1包括，
[0043] S1. 11对给定的化kagami-m相位概率密度函数fg( 0 )在定义域[c，d]上均匀采 样得到N= 1〇5点定义域采样序列X0[闲，其中Xg(i) =c+(d-c/N)*i,i= 1，;
[0044]
[0045] SI. 12将定义域采样序列Xe[闲代入相位概率密度函数f( 0 )中，得到大小为N值 域序列Ye[闲，Ye(i]=fe狂e(i));
[0046] SI. 13计算相位权重系数Ze[i]:权重系数定义为归一化的值域序列值乘W随机 样本数据序列大小M= 1〇6,并向下取整数；印 - >
[0047] S1. 14按照定义域采样值Xea)大小顺序，根据随机样本权重系数Ze山，依次将 个Ze山个X0a)加入到相位随机样本数据序列Se[M]中，序列Se[M]即为满足概率密度 分布fe( 0 )的随机样本数据；
[0048] S1. 15生成所需长度K= 105的输出随机序列0怔]。首先生成[1，M]上的整数随 机变量的均匀随机分布序列Ve怔]，即Ve怔]=l+rand(l，K)*(M-1);其次W均匀随机分布 序列值Ve怔]的排序，在随机样本数据序列中Se[M]进行取样得到目标随机序列0怔]，即 0 (j) =Se(Ve0))，j= 1，2. . . .K。此时输出的目标序列0怔]是满足服从化kagami-m 相位概率密度函数的随机序列。
[0049] S2由上述生成的化kagami-m幅度随机序列和相位随机序列，产生复随机序列 2怔]=R〇(]*eew，K= 1〇5。完成对复衰落随机序列的生成。
[0050] 下面结合仿真图对本发明相对现有技术的优点进行进一步说明。图4和图5分别 给出了不同m参数下（Q= 1)的幅度分布仿真精度和算法复杂度对比情况。权重法仿真 精度与逆变换法、舍弃法基本相同，明显优于组合法；权重法算法速度最快，其仿真时长相 比最慢的舍弃法提高了 400倍到1400倍，较组合法提高1倍；较逆变换法提高7倍左右。 图6和图7分别给出了不同m参数下（Q= 1)的相位分布仿真精度和算法复杂度对比情 况。权重法和舍弃法的均方误差基本相同，但权重法仿真速度比舍弃法快100-600倍。
[0051] 本发明提出了一种高效基于权重法的复衰落分布随机序列生成方法，该方法W少 量的空间复杂度换取时间效率和仿真精度；并且该方法的仿真精度可根据实际仿真需要， 通过调整样值采样点数和目标序列点数来实现，可避免其他方法高精度需求时插值操作带 来的误差。同时，本发明中提出的权重法实际上可模拟任意具有解析表达式的分布函数的 随机序列生成，具有很大的通用性。
【主权项】
1. 一种基于权重法的Nakagami-m复衰落随机序列生成方法，其特征在于： 根据Nakagami-m衰落衰落分布函数，并采用权重法生成幅度随机序列和相位随机序 列，并使用所述幅度随机序列和相位随机序列产生复随机序列。2. 如权利要求1所述的一种基于权重法的Nakagami-m复衰落随机序列生成方法，其特 征在于：其包括以下步骤： Sl根据Nakagami-m衰落幅度衰落分布函数，采用权重法生成幅度随机序列R[K]和相 位随机序列Θ[Κ]，其中K为目标随机序列长度； S2由上述生成的Nakagami-m幅度随机序列和相位随机序列，产生复随机序列Ζ[Κ]= R[K] *e0 [K]，e为指数函数的底数，K为目标随机序列长度。3. 如权利要求2所述的一种基于权重法的Nakagami-m衰落复衰落分布仿真方法，其特 征在于：所述步骤Sl中生成幅度随机序列R[K]包括， SI. 1对给定的Nakagami-m幅度概率密度函数fK(.)在定义域[a, b]上均勾采样得到 N点定义域采样序列XK[N],其中第i个定义域采样点XK(i) = a+(b-a/N)*i, i = 1，2…N， a为定义域下限，b为定义域上限； S1.2将定义域采样序列XK[N]代入幅度概率密度函数&(.)中，得到大小为N值域序 列YK[N]，其中第i个值域采样点YK(i] = fK(XK(i)); SI. 3计算幅度权重系数ZK[i]:权重系数定义为归一化的值域序列值乘以随机样本数 据序列大小M，并向下取整数实中M>10*K; SI. 4按照定义域采样值Xk (i)大小顺序，根据随机样本权重系数Zk [i]，依次将个Zk[i] 个Xk⑴加入到幅度随机样本数据序列SK[M]中，序列S K[M]即为满足幅度概率密度分布 fK(.)但无随机性的样本数据； SI. 5生成所需长度K的输出随机序列R[K]，其包括，首先生成[1，M]上的整数随机变 量的均匀随机分布序= l+rand(l，K)*(M-l)，rand(l，K)生成1行K列 的均匀随机序列，rand函数为标准的均匀分布生成函数；其次以均匀随机分布序列值V K[K] 的排序，在样本数据序列Sk [Μ]中进行取样得到目标随机序列R[Κ]，其中所取的第j个元素 R(j) = SK(VK(j))，j = 1，2. ... K，此时输出的目标序列R[K]是满足服从Nakagami-m幅度 概率密度函数的幅度目标随机序列。4. 如权利要求2所述的一种基于权重法的Nakagami-m衰落复衰落分布仿真方法，其特 征在于：所述步骤Sl中生成相位随机序列Θ [K]包括， SI. 6对给定的Nakagami-m相位概率密度函数f 0 (.)在定义域[c，d]上均勾采样得到 N点定义域采样序列X0 [N],其中第i个定义域采样点X0⑴=c+(d-c/N)*i, i = 1，2···Ν， c为定义域下限，d为定义域上限； SI. 7将定义域采样序列X0 [N]代入相位概率密度函数f0 (.)中，得到大小为N值域序 列Y0 [N]，其中第i个值域采样点Ye⑴=f e (Xe⑴）； SI. 9计算相位权重系数Z0 [i]:权重系数定义为归一化的值域序列值乘以随机样本数 据序列大小M，并向下取整数；@SI. 9按照定义域采样值X0 (i)大小顺序，根据随机样本权重系数Z0 [i]，依次将个 Z0 [i]个X0 (i)加入到相位随机样本数据序列S0 [M]中，序列S0 [M]即为满足相位概率密 度分布(.)的但无随机性的样本数据； SI. 10生成所需长度K的输出随机序列θ [K]，首先生成[1，M]上的整数随机变量的 均匀随机分布序列V0 〇(]，即V0 [K] = l+rand(l，K)*(M-l)，rand(l，K)生成1行K列的均 匀随机序列，rand函数为标准的均匀分布生成函数；其次以均匀随机分布序列值V 0[K]的 排序，在样本数据序列S0[M]中进行取样得到目标随机序列θ [K]，其中所取的第j个元素 9 (j) = Se (V0 (j)), j = 1, 2--K，此时输出的目标序列θ [κ]是满足服从Nakagamii^ 位概率密度函数的相位目标随机序列。
【专利摘要】一种基于权重法的Nakagami-m复衰落随机序列生成方法,根据Nakagami-m衰落幅度和相位衰落分布函数，并采用权重法分别生成幅度随机序列和相位随机序列，并使用所述幅度随机序列和相位随机序列产生复随机序列。本发明以少量的空间复杂度换取时间效率和仿真精度；可根据实际仿真需要，通过调整样值采样点数和目标序列点数来实现，可避免其他方法高精度需求时插值操作带来的误差。并且本发明实际上可模拟任意具有解析表达式的分布函数的随机序列生成，具有很大的通用性。
【IPC分类】G06F7/58
【公开号】CN104881267
【申请号】CN201510096288
【发明人】石磊, 赵蕾, 郭振, 刘彦明, 李小平
【申请人】西安电子科技大学
【公开日】2015年9月2日
【申请日】2015年3月4日