IDMA通信系統(tǒng)中的粒子群交織算法

2 IDMA通信系統(tǒng)中的檢測(cè)原理與方法
圖1為IDMA通信系統(tǒng)發(fā)送和接收部分的結(jié)構(gòu)圖。該系統(tǒng)發(fā)送端包括K個(gè)用戶，第k個(gè)用戶發(fā)送的碼元序列為：dk=[dk(1)，…，dk(i)，…，dk(I)](I為發(fā)送序列碼元長(zhǎng)度)，經(jīng)碼長(zhǎng)為S的重復(fù)編碼擴(kuò)頻得到序列為：Gk=[Gk(1)，…，Gk(j)，…，Gk(J)]，J為擴(kuò)頻后的碼長(zhǎng)，再經(jīng)交織器打亂順序重排后成為發(fā)送序列為：Xk=[xk(1)，…，xk(j)，…xk(J)]。

系統(tǒng)接收端采用Turbo型迭代譯碼結(jié)構(gòu)，主要由基本信號(hào)檢測(cè)器ESE、解交織器和K個(gè)用戶譯碼器DEC組成。系統(tǒng)中考慮完全同步無(wú)記憶信道，接收信號(hào)r(j)表示為：

式中，n(j)為高斯白噪聲采樣，hk為第k個(gè)用戶的信道衰落系數(shù)，ζk(k)為第后個(gè)用戶第j個(gè)碼片的失真，根據(jù)中心極限定理，ζk(j)滿足近似高斯分布，可用均值和方差函數(shù)表示。
基本信號(hào)檢測(cè)器模塊ESE的檢測(cè)算法概括如下：

式中E(．)表示均值，V(．)表示方差。
譯碼器DEC反饋的先驗(yàn)信息經(jīng)交織后得到eDEC[xk(j)]，更新ESE中接收信號(hào)的均值和方差，并產(chǎn)生輸入到譯碼器DEC的外信息eESE(xk(j))，外信息經(jīng)過(guò)解交織，作為譯碼器DEC輸入端的先驗(yàn)信息，如此循環(huán)迭代直到規(guī)定的次數(shù)以后，K個(gè)用戶的譯碼器分別產(chǎn)生相應(yīng)信息序列的硬判決值dk。

3 基于粒子群的交織算法
3．1 粒子群算法的相關(guān)定義與操作
粒子群算法是一種基于迭代的優(yōu)化方法，具有易于實(shí)現(xiàn)、需要調(diào)整的參數(shù)少等優(yōu)點(diǎn)，并且在較少的迭代次數(shù)的情況下就可獲得全局最優(yōu)解。設(shè)粒子群粒子個(gè)數(shù)為M，在一個(gè)D維的搜索空間中，粒子i(1≤i≤M)在第k次迭代時(shí)的位置信息為速度信息表示為搜索空間的任意數(shù)，即粒子到目前為止所經(jīng)歷的最好位置為群體中所有粒子到目前為止所經(jīng)歷過(guò)的最好位置為其中b為具有最優(yōu)位置粒子的索引。為使粒子群算法能正確解決交織問(wèn)題，這里定義粒子的位置和速度的含義及相關(guān)操作如下：
(1)粒子定義粒子的位置定義為一序列。對(duì)于搜索空間為D維的種群，序列長(zhǎng)度為D。假設(shè)某個(gè)粒子j的位置為xj，D=5，則可表示為xj=(1，2，3，4，5)；速度定義為粒子位置的變換集，即一組置換序列的有序列表，表示為：v={(ik,jk)，ik,jk∈D}，k∈{1，2，…，|v|，其中|v|表示該速度所含置換序列的個(gè)數(shù)。
(2)置換操作 假設(shè)某粒子i的位置為xi，定義置換序列(mi,ni)，置換操作用以交換xi中第mi和ni個(gè)值的位置，則x'i=xi+(mi,ni)，其中x'i為經(jīng)過(guò)置換操作后得到的新位置。
(3)加法操作包括粒子速度與速度的加法操作及粒子速度和位置的加法操作。設(shè)vi、vj、vk分別表示第i、第j及第k個(gè)粒子速度，xk為第k個(gè)粒子位置。vi+vj表示兩個(gè)速度相加的操作，其結(jié)果為兩個(gè)置換序列合并，產(chǎn)生一個(gè)新置換序列串；vk+xk表示速度和位置的加法操作，即將一組置換序列依次作用于某個(gè)粒子位置。其結(jié)果為一個(gè)新位置。
(4)減法操作主要指粒子位置與位置的減法操作。該操作相減后結(jié)果為一組置換序列，即速度。設(shè)xk為第k個(gè)粒子位置，xi為第i個(gè)粒子位置，則xi-xk為一個(gè)置換序列。例如：xi=(1，2，3，4，5)，xk=(2，3，1，4，5)，由于xi(1)=1，xk(3)=1，第1個(gè)交換序列為(1，3)，xj=xk(1，3)=(1，3，2，4，5)；又xi(2)=2，xi(3)=2，第2個(gè)交換序列為(2，3)，xi=xj+(2，3)，因此經(jīng)上述操作得到：xi=xk+{(1，3)，(2，3)}，所以xi-xk={(1，3)，(2，3)}。
(5)乘法操作 指實(shí)數(shù)與粒子速度的乘法操作。對(duì)于在(0，1)任意實(shí)數(shù)c，設(shè)速度v有i個(gè)置換序列，則乘法操作截取速度置換序列，使新的速度置換序列個(gè)數(shù)為|cxi|(cxi取整)。
根據(jù)以上定義，粒子群的更新公式可描述為：

式中，c1與c2是兩個(gè)正的常數(shù)，稱為加速因子，r1和r2為分布于[0，1]間的隨機(jī)數(shù)。
3．2 粒子群交織算法
對(duì)于碼片長(zhǎng)度為J的序列，其交織方式有J!種。當(dāng)J較小，而用戶數(shù)較大時(shí)，隨機(jī)產(chǎn)生的交織序列之間的互相關(guān)系數(shù)接近1，交織區(qū)分用戶時(shí)，嚴(yán)重影響兩個(gè)用戶間通信的性能。粒子群算法能搜索到全局最優(yōu)解，可選擇互相關(guān)性最弱的交織序列區(qū)分用戶，從而提高通信性能。
在數(shù)據(jù)傳輸中，碼片長(zhǎng)度J=IxS，其中I為傳輸序列碼元長(zhǎng)度。S為擴(kuò)頻碼長(zhǎng)度。粒子群交織是從J!個(gè)交織序列中選擇K個(gè)互相關(guān)性弱的序列區(qū)分用戶，但J!個(gè)解在實(shí)際操作中運(yùn)算量較大，故選擇N個(gè)作為初始解。粒子群交織算法中選擇互相關(guān)矩陣作為適應(yīng)度函數(shù)。X為N×J的數(shù)據(jù)矩陣，是N個(gè)可能的解，每個(gè)解是長(zhǎng)度為J的交織序列，其元素為X(n,J)∈｛－1，1}(n=1，…，N，j=1，…，J)；INDEX表示一個(gè)，N×J的矩陣，其值是對(duì)應(yīng)數(shù)據(jù)矩陣X的交織序列的索引值；R是X的互相關(guān)系數(shù)矩陣，為N×N維，R的元素R(i，j)由下式計(jì)算得到：

式中，Xi是矩陣X的第i行向量，μi=E(Xi)，E表示數(shù)學(xué)期望。