一種基于權(quán)重法的Nakagami-m復衰落隨機序列生成方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于無線通信信道建模及模擬領(lǐng)域，具體來說本發(fā)明設及一種針對 化kagami-m復衰落隨機序列生成方法。
【背景技術(shù)】
[0002] 無線移動通信環(huán)境中，化kagami-m衰落相較于傳統(tǒng)的瑞利分布或萊斯分布可W更 好更廣泛地表征小尺度衰落特征?；痥agami-m衰落是通過變參數(shù)伽馬分布的密度函數(shù)擬合 實驗數(shù)據(jù)所得到的。它可W模擬從嚴重、適中、輕微到無衰落的信道環(huán)境，包含了瑞利分布 和萊斯分布。
[0003] 目前化kagami-m衰落的仿真主要是針對其包絡分布的仿真，主要有逆變換法、舍 棄法和組合法；針對化kagami-m衰落相位分布的仿真方法較少，一般采用舍棄法。組合 法是一種簡單的擬合模型，仿真速度快，但是誤差較大，特別是m值不在其適用范圍之內(nèi) 時，仿真結(jié)果與理論值幾乎完全不符。逆變換模型實用性較差精度也不高，原因在于求解 化kagami-m累積分布函數(shù)的反函數(shù)十分復雜且需要系數(shù)的擬合捜索。舍棄法精度高，但舍 棄帽子函數(shù)選擇復雜，仿真速度慢，效率較低。因此，有必要研究新的高效的復衰落信道建 模方法。

【發(fā)明內(nèi)容】

[0004] 為了達到上述目的，本發(fā)明提出一種基于權(quán)重法的化kagami-m復衰落分布仿真 方法，所提權(quán)重法根據(jù)概率密度函數(shù)的解析式生成服從化kagami-m衰落分布的大量隨機 樣本數(shù)據(jù)，然后對樣本數(shù)據(jù)進行均勻隨機采樣，得到仿真精度所需要的目標隨機序列。一種 基于權(quán)重法的化kagami-m復衰落隨機序列生成方法，根據(jù)化kagami-m衰落分布函數(shù)，并采 用權(quán)重法生成幅度隨機序列和相位隨機序列，并使用所述幅度隨機序列和相位隨機序列產(chǎn) 生復隨機序列。
[0005] 在上述技術(shù)方案的基礎(chǔ)上，其包括W下步驟：
[0006] S1根據(jù)化kagami-m衰落幅度和相位衰落分布函數(shù)，采用權(quán)重法生成幅度隨機序 列RDC和相位隨機序列0怔]，其中K為目標隨機序列長度；
[0007]S2由上述生成的化kagami-m幅度隨機序列和相位隨機序列，產(chǎn)生復隨機序列 ZDC=R〇(]*e9w，K為目標隨機序列長度。
[000引在上述技術(shù)方案的基礎(chǔ)上，所述步驟S1中生成幅度隨機序列RDC包括，
[0009] S1. 1對給定的化kagami-m幅度概率密度函數(shù)fijC)在定義域[a,b]上均勻采樣得 至IJN點定義域采樣序列而[閑，其中第i個定義域采樣點而(1) =a+化-a/N)*i，i= 1，2，" N，a為定義域下限，b為定義域上限；
[0010] S1. 2將定義域采樣序列而[閑代入幅度概率密度函數(shù)fE(.)中，得到大小為N值 域序列Ye[閑，其中第i個值域采樣點YeQ] =fE狂kQ));
[0011] SI. 3計算幅度權(quán)重系數(shù)Z,[i]:權(quán)重系數(shù)定義為歸一化的值域序列值乘w隨機樣 本數(shù)據(jù)序列大小M，并向下取整數(shù)；即 其中M〉10*K;
，
[0012]S1. 4按照定義域采樣值XeQ)大小順序，根據(jù)隨機樣本權(quán)重系數(shù)與[1]，依次將個 與山個X^a)加入到幅度隨機樣本數(shù)據(jù)序列Se[M]中，序列Se[M]即為滿足幅度概率密度 分布fK(.)但無隨機性的本數(shù)據(jù)；
[0013]S1.5生成所需長度K的輸出隨機序列RDC，其包括，首先生成[1，M]上的整數(shù)隨 機變量的均勻隨機分布序列Vk怔]，即Vk怔]=l+rand(l，K)*(M-l)，rancKl，K)生成1行 K列的均勻隨機序列，rand函數(shù)為標準的均勻分布生成函數(shù)；其次W均勻隨機分布序列值 Vc[時的排序，在樣本數(shù)據(jù)序列Sc[M]中進行取樣得到目標隨機序列R[K]，其中所取的第j 個元素R(j) =Sk(VkU))，j= 1，2--K，此時輸出的目標序列R[K]是滿足服從r^Jakagami-m 幅度概率密度函數(shù)的幅度目標隨機序列。
[0014] 在上述技術(shù)方案的基礎(chǔ)上，所述步驟S1中生成相位隨機序列0怔]包括，
[0015]S1. 6對給定的化kagami-m相位概率密度函數(shù)f0(.)在定義域[C，d]上均勻采 樣得到N點定義域采樣序列Xe[閑，其中第i個定義域采樣點Xe(i) =c+(d-c/N)*i，i= 1，2…N，c為定義域下限，d為定義域上限；
[0016]S1. 7將定義域采樣序列Xe[閑代入相位概率密度函數(shù)fe(.)中，得到大小為N值 域序列Ye[閑，其中第i個值域采樣點Ye(i] =fe狂e(i));
[0017] SI. 8計算相位權(quán)重系數(shù)Ze[i]:權(quán)重系數(shù)定義為歸一化的值域序列值乘w隨機 樣本數(shù)據(jù)序列大小M，并向下取整數(shù)；良F
[0018]S1.9按照定義域采樣值Xe(i)大小順序，根據(jù)隨機樣本權(quán)重系數(shù)Ze[i]，依次將個 Zg[i]個Xe(i)加入到相位隨機樣本數(shù)據(jù)序列Se[M]中，序列Se[M]即為滿足相位概率密 度分布fe(.)的但無隨機性的樣本數(shù)據(jù)；
[0019]S1. 10生成所需長度K的輸出隨機序列0怔]。首先生成[1，M]上的整數(shù)隨機變量 的均勻隨機分布序列Ve怔]，即Ve怔]=l+rand(l，K)*(M-1)，rand(l，K)生成1行K列的 均勻隨機序列，rand函數(shù)為標準的均勻分布生成函數(shù)；其次W均勻隨機分布序列值Ve怔] 的排序，在樣本數(shù)據(jù)序列Se[M]中進行取樣得到目標隨機序列0 [時，其中所取的第j個元 素0 (j) =Sg(Vg(j))，j= 1，2--K，此時輸出的目標序列0 [K]是滿足服從r^Jakagami-m 相位概率密度函數(shù)的相位目標隨機序列。
[0020] 與現(xiàn)有技術(shù)相比較，本發(fā)明的有益效果在于：
[0021] (1)本發(fā)明提出了一種高效基于權(quán)重法的復衰落分布隨機序列生成方法，該方法 W少量的空間復雜度換取時間效率和仿真精度；
[0022] (2)該方法的仿真精度可根據(jù)實際仿真需要，通過調(diào)整樣值采樣點數(shù)和目標序列 點數(shù)來實現(xiàn)，可避免其他方法高精度需求時插值操作帶來的誤差。
[0023](3)本發(fā)明中提出的權(quán)重法實際上可模擬任意具有解析表達式的分布函數(shù)的隨機 序列生成，具有很大的通用性。
【附圖說明】
[0024] 圖1為本發(fā)明一種基于權(quán)重法的化kagami-m復衰落隨機序列生成方法流程圖；
[0025] 圖2為本發(fā)明權(quán)重法生成任意具有解析式表達的概率密度函數(shù)分布隨機序列流 程圖；
[0026] 圖3為本發(fā)明m=4, Q=1的化kagami-m幅度包絡序列仿真結(jié)果與其他方法對 比圖；
[0027] 圖4為本發(fā)明Q=1的化kagami-m幅度包絡序列仿真精度與其他方法對比圖； [002引圖5為本發(fā)明Q=1的化kagami-m幅度包絡序列仿真歸一化時間復雜度與其他 方法對比圖；
[0029] 圖6為本發(fā)明m=4, Q=1的化kagami-m相位包絡序列仿真結(jié)果與其他方法對 比圖；
[0030] 圖7為本發(fā)明Q=1的化kagami-m相位包絡序列仿真精度與其他方法對比圖；
[003U 圖8為本發(fā)明Q=1的化kagami-m相位包絡序列仿真歸一化時間復雜度與其他 方法對比圖。 具體實施例
[0032] 下面結(jié)合具體實施例對本發(fā)明進行詳細說明?，F(xiàn)針對給定m=4, Q=1時的 化kagami-m衰落復衰落分布仿真，本發(fā)明過程具體實施方案如下；
[0033] 請參考圖1與圖2。
[0034]S1根據(jù)化kagami-m衰落幅度衰落分布函數(shù)，采用權(quán)重法生成幅度隨機序列R[K] 和相位隨機序列0 [K]，K= 1〇5;其中生成幅度隨機序列R[K]和相位隨機序列0怔]為兩 并列的步驟，首先生成幅度隨機序列或先生成相位隨機序列均可。
[003引所述步驟S1包括，
[0036]S1. 1對給定的化kagami-m幅度概率密度函數(shù)f,(r)在定義域[a,b]上均勻采樣 得到N= 1〇5點定義域采樣序列Xe[閑，其中而(1) =a+化-a/N)*i，i= 1，2"，N ;
[0037]
[00測 SI. 2將定義域采樣序列而[閑代入幅度概率密度函數(shù)fK(r)中，得到N=105點值 域序列Yk[閑，YeU] =fK狂K(i));
[0039] SI. 3計算幅度權(quán)重系數(shù)Z,[i]:權(quán)重系數(shù)定義為歸一化的值域序列值乘W隨機樣 本數(shù)據(jù)序列大小M= 1〇6,并向下取整數(shù)；即 - >
[0040]S1. 4按照定義域采樣值XeQ)大小