一种基于权重法的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)大小