潘微宇，趙知?jiǎng)?，2

（1.杭州電子科技大學(xué) 通信工程學(xué)院，浙江 杭州 310018；2.中國(guó)電子科技集團(tuán)第36 研究所 通信系統(tǒng)信息控制技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，浙江 嘉興 314001）

0 引言

直接擴(kuò)頻序列碼分多址（Direct Sequence Code Division Multiple Access，DS-CDMA）采用高速率的偽碼對(duì)信息碼進(jìn)行擴(kuò)頻調(diào)制，具有較強(qiáng)的隱蔽性和抗干擾能力[1]。DS-CDMA 信號(hào)可根據(jù)擴(kuò)頻方式分為短碼（Short Code，SC）DS-CDMA 和長(zhǎng)碼（Long Code，LC）DS-CDMA 兩種信號(hào)，為了進(jìn)一步提高保密性，利用長(zhǎng)碼對(duì)短碼擴(kuò)頻信號(hào)做加擾處理，得到長(zhǎng)短碼（Long and Short Code，LSC）DS-CDMA 信號(hào)。SC-DS-CDMA 信號(hào)用一周期擴(kuò)頻碼調(diào)制一個(gè)信息碼，文獻(xiàn)[2]-[3]分別采用多重信號(hào)分類(lèi)法[2]和平行因子法[3]實(shí)現(xiàn)了各用戶偽碼盲估計(jì)；文獻(xiàn)[4]利用特征值分解法和Givens 矩陣旋轉(zhuǎn)法來(lái)估計(jì)SC-DS-CDMA 信號(hào)各用戶的偽碼，但在尋找Givens 矩陣旋轉(zhuǎn)角度時(shí)采用極值法，難以得到最佳旋轉(zhuǎn)角度，估計(jì)性能不佳；LC-DS-CDMA 信號(hào)用一周期擴(kuò)頻碼調(diào)制多個(gè)信息碼，文獻(xiàn)[5-6]通過(guò)構(gòu)建張量模型，分別利用TUCKER 分解[5]和CP 分解[6]估計(jì)LCDS-CDMA 信號(hào)的偽碼；文獻(xiàn)[7]采用迭代最小二乘投影算法和庫(kù)搜索估計(jì)異步LC-DS-CDMA 信號(hào)各用戶的偽碼；文獻(xiàn)[8]采用矩陣填充理論、奇異值分解和FastICA 算法估計(jì)LSC-DS-CDMA 信號(hào)各用戶復(fù)合碼，并利用m 序列三階相關(guān)特性估計(jì)出各用戶長(zhǎng)短碼；文獻(xiàn)[9]利用FastICA 算法分離LSC-DS-CDMA 信號(hào)各用戶復(fù)合碼，并基于分圓陪集理論和三階相關(guān)特性估計(jì)長(zhǎng)短碼。上述研究都只針對(duì)理想信號(hào)，實(shí)際接收信號(hào)總會(huì)帶有小的殘余頻偏[10]，使得上述方法無(wú)法直接使用。文獻(xiàn)[11]-[12]分別利用線性調(diào)頻Z 變換法和無(wú)跡卡爾曼濾波法估計(jì)殘余頻偏，但無(wú)法完全消除頻偏；文獻(xiàn)[13]利用特征值分解法得到帶殘余頻偏的偽碼，并通過(guò)全數(shù)字鎖相環(huán)（Digital Phase Locked Loop，DPLL）對(duì)頻偏進(jìn)行跟蹤消除，實(shí)現(xiàn)了長(zhǎng)碼直擴(kuò)信號(hào)的偽碼估計(jì)，但并未對(duì)LSC-DS-CDMA 信號(hào)進(jìn)行研究。

針對(duì)帶殘余頻偏的LSC-DS-CDMA 信號(hào)，本文按長(zhǎng)碼周期和短碼周期進(jìn)行二次分段構(gòu)成信號(hào)矩陣，利用特征值分解、全數(shù)字鎖相環(huán)和Givens 旋轉(zhuǎn)消除殘余頻偏和估計(jì)得到各用戶復(fù)合碼，提出了結(jié)合隨機(jī)差分和模擬退火思想的鯨魚(yú)優(yōu)化算法來(lái)求解總方差最小值的Givens矩陣的最佳旋轉(zhuǎn)角，最后根據(jù)梅西算法和相關(guān)運(yùn)算估計(jì)各用戶的長(zhǎng)碼和短碼。

1 信號(hào)模型

假設(shè)LSC-DS-CDMA 信號(hào)中各用戶擴(kuò)頻短碼和長(zhǎng)擾碼的碼片速率相等，且已完成盲同步[8]，以擴(kuò)頻短碼碼片速率對(duì)接收到的信號(hào)進(jìn)行采樣，則包含R 個(gè)用戶的帶殘余載波的LSC-DS-CDMA 信號(hào)模型為：

其中，f△為殘余頻偏；φ ∈（0，2π）表示均勻分布的隨機(jī)相位；Ar、dr（n）、sr（n）分別表示第r 個(gè)用戶的信號(hào)幅度、信息碼和復(fù)合碼，sr（n）=br（n）cr（n），br（n）表示周期為Nb的擴(kuò)頻短碼，選用Walsh碼，cr（n）表示周期為Nc的長(zhǎng)擾碼，選用m 序列；各個(gè)用戶的信息碼、短碼和長(zhǎng)擾碼均相互獨(dú)立；Nh=PNc為接收信號(hào)長(zhǎng)度，P 表示信號(hào)周期個(gè)數(shù)，M=「Nc/Nb表示一周期復(fù)合碼所調(diào)制的信息碼個(gè)數(shù)，“「”為向上取整；w（n）表示均值為零、方差為σ2的復(fù)高斯白噪聲。

先以長(zhǎng)擾碼周期Nc對(duì)接收信號(hào)進(jìn)行分段，yp=[y（（p-1）Nc+1）y（（p-1）Nc+2）…y（pNc）]H，p=1，2，…，P。由于一周期的長(zhǎng)擾碼包含M 個(gè)擴(kuò)頻短碼，因此可利用短碼周期Nb對(duì)yp再進(jìn)行分段，ypm=[y（（p-1）Nc+（m-1）Nb+1）y（（p-1）Nc+（m-1）Nb+2）…y（（p-1）Nc+mNb）]H，將不同接收信號(hào)周期p中的第m 段信號(hào)依次進(jìn)行排列，可得共M 個(gè)信號(hào)矩陣，表示為：

由式（2）可知，每個(gè)信號(hào)矩陣中均包含了各用戶的復(fù)合碼片段信息和載波頻偏，因此首先需要消除載波頻偏影響，依次估計(jì)出各用戶完整的復(fù)合碼，最后可由復(fù)合碼估計(jì)得到各用戶的長(zhǎng)短碼。

2 基于DPLL 和Givens 旋轉(zhuǎn)的復(fù)合碼估計(jì)

令z（n）=ej（2πf△n+φ），矩陣ypm的向量形式可表示為：

其中，wpm表示第周期p 中的第m 段信號(hào)的噪聲，srm=[sr（（m-1）Nb+1）…sr（mNb）]H，dr[（p-1）M+m]表示srm對(duì)應(yīng)的信息 碼。Zpm=diag{z[（p-1）Nc+（m-1）Nb+1]，z[（p-1）Nc+（m-1）Nb+2]，…，z[（p-1）Nc+mNb]}，還可表示為：

其中，Z11=diag{ej（2πf△nb）}，由于各用戶間擴(kuò)頻短碼和長(zhǎng)擾碼相互獨(dú)立，且不同用戶的信息碼與噪聲統(tǒng)計(jì)獨(dú)立[14]，結(jié)合式（4），可得ypm的自相關(guān)矩陣為：

2.1 殘余頻偏消除

記Gm=Z11Hm=Z11[h1m，h2m，…，hRm]，hrm=[hrm（1），hrm（2），…，hrm（Nb）]H表示Hm中第r個(gè)特征向量。由于Z11=diag{ej（2πf△nb）}，則Grm=Z11hrm可以展開(kāi)為：

利用DPLL 消除殘余載波[13]，得到第r 個(gè)特征向量的估計(jì)值rm。將Gm中R 個(gè)特征向量依次送入DPLL 進(jìn)行殘余載波的跟蹤消除，得到Hm的估計(jì)值

2.2 復(fù)合碼估計(jì)

信號(hào)子空間Hm中的R 個(gè)特征向量無(wú)法保證與R個(gè)用戶的復(fù)合碼一致，從空間幾何的角度看，Hm和真正的復(fù)合碼序列子空間Sm=[s1m，s2m，…，sRm]存在一定的夾角[4]，利用Givens 旋轉(zhuǎn)可實(shí)現(xiàn)信號(hào)子空間對(duì)應(yīng)的角度旋轉(zhuǎn)，以得到真正的復(fù)合碼序列子空間估計(jì)。Givens 旋轉(zhuǎn)矩陣表示為：

其中，θr1，r2∈（0，2π），0＜r1，r2≤R。將視為荷載矩陣，利用（R-1）個(gè)Givens 旋轉(zhuǎn)矩陣T1，2，T2，3，…，TR-1，R對(duì)其進(jìn)行旋轉(zhuǎn)，即：

3 鯨魚(yú)優(yōu)化算法的改進(jìn)及偽碼估計(jì)

3.1 鯨魚(yú)優(yōu)化算法（WOA）

WOA 算法具體分為包圍獵物、螺旋泡沫網(wǎng)狩獵和搜索獵物三部分。利用適應(yīng)度函數(shù)f（X（t+1））來(lái)衡量鯨魚(yú)，令包圍獵物和螺旋泡沫網(wǎng)狩獵兩種方式的概率p 相同[17]，當(dāng)p＜0.5 時(shí)采取包圍獵物方式，p≥0.5 時(shí)采取螺旋泡沫網(wǎng)狩獵方式，更新公式為：

其中，t 為當(dāng)前迭代次數(shù)；X（t）為當(dāng)前鯨魚(yú)個(gè)體位置，X*（t）為當(dāng)前最優(yōu)個(gè)體位置；b為常數(shù)，ε 為[-1，1]間的隨機(jī)數(shù)，p 為[0，1]間的隨機(jī)數(shù)；I1和I2為系數(shù)變量，分別表示為：

式中，i1和i2為 [0，1]間的隨機(jī) 向量，tmax為最大迭代次數(shù)，α 為由2 線性遞減至0的收斂因子。

鯨魚(yú)除了圍捕獵物外，還可以隨機(jī)尋找新的獵物，當(dāng)|I1|＜1 時(shí)利用式（10）進(jìn)行位置更新，當(dāng)|I1|≥1時(shí)隨機(jī)選擇一個(gè)鯨魚(yú)個(gè)體作為當(dāng)前最優(yōu)個(gè)體，其他個(gè)體向其靠攏，表達(dá)式為：

其中，Xr（t）為隨機(jī)選擇的鯨魚(yú)個(gè)體位置。

3.2 結(jié)合模擬退火思想的隨機(jī)差分變異鯨魚(yú)優(yōu)化算法（SA-IWOA）

為了提升鯨魚(yú)算法的全局搜索能力，首先對(duì)收斂因子α 改進(jìn)如下：

改進(jìn)后的α 在算法初期數(shù)值較大且緩慢減小，可以擴(kuò)大搜索范圍，加強(qiáng)全局搜索能力；在后期α 減小快，加強(qiáng)局部搜索能力。

然后利用隨機(jī)差分變異策略對(duì)位置更新后的鯨魚(yú)種群進(jìn)行如式（16）所示的擾動(dòng)更新：

其中，“·”為點(diǎn)乘運(yùn)算，I3為[0，1]間的隨機(jī)向量，Xp（·）表示擾動(dòng)位置，（t）為隨機(jī)選擇的鯨魚(yú)個(gè)體位置。增強(qiáng)種群個(gè)體的多樣性，防止算法早熟而陷入局部最優(yōu)。

最后利用模擬退火的思想，在每次擾動(dòng)更新后，如果新個(gè)體的適應(yīng)度函數(shù)f（Xp（t+1））比更新前的f（X（t））小，則接受更新；否則通過(guò)接受概率判斷是否更新，這樣使種群個(gè)體往最優(yōu)解靠近的同時(shí)，使算法具有跳出局部最優(yōu)的能力。接受概率γ 為：

其中，β表示溫度。

3.3 基于SA-IWOA 的最優(yōu)角度求解

令X=[θ1，2，θ2，3，…，θR-1，R]為鯨魚(yú)個(gè)體，將矩陣的總方差值視作適應(yīng)度函數(shù)，即，根據(jù)式（7）～式（9），可由X 計(jì)算得到。利用SA-IWOA 求解的具體步驟如下：

（1）初始化算法參數(shù)：種群數(shù)目J、最大迭代次數(shù)tmax、初始溫度β、i1和i2；

（3）根據(jù)式（10）和式（14）更新鯨魚(yú)位置，利用式（7）計(jì)算Givens 旋轉(zhuǎn)矩陣?yán)檬剑?）計(jì)算m，計(jì)算各鯨魚(yú)個(gè)體的f（X（t））；

（4）根據(jù)式（16）對(duì)位置更新后的鯨魚(yú)個(gè)體進(jìn)行擾動(dòng)更新，利用式（7）計(jì)算Givens 旋轉(zhuǎn)矩陣?yán)檬剑?）計(jì)算m，計(jì)算每個(gè)鯨魚(yú)個(gè)體的f（Xp（t+1））；

（5）若f（Xp（t+1））＜f（X（t）），則X（t+1）=Xp（t+1）；否則，產(chǎn)生一個(gè)隨機(jī)數(shù)I4∈[0，1]，利用式（17）計(jì)算γ，當(dāng)I4落在區(qū)間[γ，1]內(nèi)時(shí)，X（t+1）=Xp（t+1），反之X（t+1）=X（t+1），然后進(jìn)行退溫操作：β′=0.95β；

（6）記錄此時(shí)總方差值最小的鯨魚(yú)個(gè)體即最優(yōu)個(gè)體Xg，判斷迭代次數(shù)是否為最大值，若是，則輸出最優(yōu)個(gè)體，否則分別根據(jù)式（11）、式（12）和式（15）更新參數(shù)I1、I2和α；返回步驟（3）繼續(xù)執(zhí)行。

3.4 長(zhǎng)擾碼和擴(kuò)頻短碼的估計(jì)

由式（2）可知，依次估計(jì)出M 個(gè)信號(hào)矩陣Ym中的各用戶復(fù)合碼片段，然后按順序進(jìn)行拼接，即可得到共RM個(gè)長(zhǎng)度為Nc的模糊序列其 中。R 個(gè)用戶的復(fù)合碼就包含在這RM個(gè)模糊序列中，將其記為。但由于每段復(fù)合碼都可以獨(dú)立取正負(fù)號(hào)，導(dǎo)致存在幅度模糊的問(wèn)題，且還存在短碼的影響，無(wú)法估計(jì)各用戶的長(zhǎng)擾碼，通過(guò)兩次延遲相乘的方式消除幅度模糊和短碼影響：

其中，1≤n≤Nb，1≤λ1≤R。

所以可得本文提出的帶殘余載波的LSC-DS-CDMA信號(hào)偽碼估計(jì)算法（簡(jiǎn)記為L(zhǎng)SC-SA-IWOA）主要步驟如下：

（1）利用長(zhǎng)碼和短碼周期對(duì)接收信號(hào)進(jìn)行2 次分段，構(gòu)成帶殘余載波的接收信號(hào)矩陣Ym；

（2）對(duì)信號(hào)矩陣進(jìn)行特征值分解，得到帶殘余載波的信號(hào)子空間Gm；

（4）利用SA-IWOA 算法求解最優(yōu)旋轉(zhuǎn)角，并利用Givens旋轉(zhuǎn)矩陣對(duì)信號(hào)子空間進(jìn)行旋轉(zhuǎn)，得到各用戶的復(fù)合碼估計(jì)m；

（5）通過(guò)兩次延遲相乘的方式消除幅度模糊和短碼影響，并根據(jù)梅西算法和相關(guān)運(yùn)算分別獲得長(zhǎng)碼和短碼估計(jì)。

4 算法仿真與性能分析

在仿真實(shí)驗(yàn)中，各用戶的信號(hào)幅度相同，即Ar=A；信息碼是隨機(jī)生成的序列{+1，-1}；定義信噪比（Signal Noise Ratio，SNR）為表示接收信號(hào)方差；長(zhǎng)擾碼使用8 階m 序列，擴(kuò)頻短碼使用周期長(zhǎng)度為Nb=64 的OVSF 碼；種群數(shù)目J=50，最大迭代次數(shù)tmax=400，初始溫度β=1。長(zhǎng)碼和短碼的誤碼率的平均值稱(chēng)為長(zhǎng)短碼誤碼率，用其衡量算法的偽碼估計(jì)性能。以下仿真結(jié)果均為400 次蒙特卡洛實(shí)驗(yàn)的平均。

4.1 實(shí)驗(yàn)1：不同頻偏對(duì)長(zhǎng)短碼估計(jì)性能的影響

接收信號(hào)周期數(shù)P=100，用戶數(shù)R=2，當(dāng)頻偏分別為f△=0、f△=±0.02、f△=±0.03、f△=±0.04時(shí)，本文算法對(duì)長(zhǎng)短碼估計(jì)的長(zhǎng)短碼誤碼率曲線如圖1 所示。

圖1 不同頻偏對(duì)長(zhǎng)短碼估計(jì)性能的影響

由圖1 可知，當(dāng)|f△|相同時(shí)，誤碼率基本相同；當(dāng)|f△|逐漸增大時(shí)，誤碼率也逐漸升高，當(dāng)頻偏分別為f△=0、f△=0.02、f△=0.03 和f△=0.04時(shí)，誤碼率低于1%時(shí)本文算法所需的信噪比分別為-7.1 dB、-5.5 dB、-4.8 dB、-3.9 dB，說(shuō)明DPLL 能有效去除殘余頻偏。

4.2 實(shí)驗(yàn)2：不同周期數(shù)對(duì)長(zhǎng)短碼估計(jì)性能的影響

殘余頻偏f△=0.02，用戶數(shù)R=2，當(dāng)接收信號(hào)周期數(shù)P 分別為50、100 和150時(shí)，本文算法的長(zhǎng)短碼誤碼率曲線如圖2 所示。

圖2 不同周期數(shù)對(duì)長(zhǎng)短碼估計(jì)性能的影響

由圖2 可知，接收信號(hào)周期數(shù)越多，誤碼率越低。當(dāng)接收信號(hào)周期數(shù)P 分別為50、100 和150，誤碼率低于1%時(shí)本文算法所需的信噪比分別為-2.5 dB、-5.5 dB、-6.9 dB。因?yàn)榻邮招盘?hào)周期數(shù)越多，其包含的有用信息更多，特征值分解法的結(jié)果越準(zhǔn)確，利用Givens 矩陣旋轉(zhuǎn)和SA-IWOA 算法估計(jì)的復(fù)合碼也就越準(zhǔn)確，最終利用梅西算法和相關(guān)運(yùn)算估計(jì)長(zhǎng)短碼的誤碼率就越低。

4.3 實(shí)驗(yàn)3：不同用戶數(shù)對(duì)長(zhǎng)短碼估計(jì)性能的影響

接收信號(hào)周期數(shù)P=100，殘余頻偏f△=0.02，當(dāng)用戶數(shù)R 分別為2、3 和4時(shí)，本文算法的長(zhǎng)短碼誤碼率曲線如圖3 所示。

圖3 不同用戶數(shù)對(duì)長(zhǎng)短碼估計(jì)性能的影響

由圖3 可知，用戶數(shù)越多，長(zhǎng)短碼估計(jì)性能越差。因?yàn)橛脩魯?shù)越多，用戶間信號(hào)干擾越大，導(dǎo)致特征值分解法的結(jié)果誤差越大，使得Givens 矩陣旋轉(zhuǎn)的誤差越大，導(dǎo)致最終的長(zhǎng)短碼估計(jì)性能越差。

4.4 實(shí)驗(yàn)4：鯨魚(yú)算法對(duì)比

接收信號(hào)周期數(shù)P=100，殘余頻偏f△=0.02，用戶數(shù)R=2，當(dāng)信噪比為-9 dB時(shí)，應(yīng)用本文所提出的SA-IWOA算法、混沌鯨魚(yú)算法（COWA）[17]和隨機(jī)差分鯨魚(yú)算法（IWOA）[20]得到的適應(yīng)度函數(shù)即總方差值隨迭代次數(shù)的變化如圖4 所示。

圖4 總方差值隨迭代次數(shù)的變化圖

由圖4 可知，本文所提的SA-IWOA 算法收斂速度最快，求得的總方差值最小，即估計(jì)的旋轉(zhuǎn)角度優(yōu)于IWOA和COWA 算法的。

4.5 實(shí)驗(yàn)5：偽碼估計(jì)算法對(duì)比

實(shí)驗(yàn)條件同實(shí)驗(yàn)4，將CWOA 和IWOA 用于本文算法中分別簡(jiǎn)記為L(zhǎng)SC-CWOA 和LSC-IWOA，本文算法LSC-SA-IWOA 和LSC-CWOA、LSC-IWOA、文獻(xiàn)[4]算法的長(zhǎng)短碼誤碼率如圖5 所示。

圖5 長(zhǎng)短碼誤碼率對(duì)比

由圖5 可知，本文算法的誤碼率最低，其次是LSCIWOA 算法，然后是LSC-CWOA 算法，最后是文獻(xiàn)[4]的極值法。極值法在估計(jì)總方差的最小值時(shí)，極易得到局部最優(yōu)解，缺乏全局搜索的能力；CWOA 利用混沌序列增加了全局搜索的能力，但局部搜索的能力不強(qiáng)；IWOA 利用隨機(jī)差分變異的方式增加全局搜索的能力，但變異的結(jié)果隨機(jī)，差的變異結(jié)果在一定程度上影響了最優(yōu)角度的求解。

5 結(jié)論

針對(duì)帶殘余載波的LSC-DS-CDMA 信號(hào)的偽碼估計(jì)，先根據(jù)長(zhǎng)碼周期和短碼周期對(duì)接收信號(hào)進(jìn)行分段并構(gòu)建矩陣；接著利用特征值分解和DPLL 得到消除殘余頻偏的用戶特征向量子空間；再根據(jù)Givens 矩陣旋轉(zhuǎn)，利用結(jié)合SA-IWOA 算法求出最佳旋轉(zhuǎn)角，估計(jì)出各用戶的復(fù)合碼序列并依次拼接；最后利用梅西算法和相關(guān)運(yùn)算估計(jì)各用戶的長(zhǎng)碼和短碼。仿真實(shí)驗(yàn)表明，本文算法能有效估計(jì)帶殘余頻偏的LSC-DS-CDMA 信號(hào)偽碼，估計(jì)性能優(yōu)于文獻(xiàn)[4]的極值法，所提出的SA-IWOA 算法的尋優(yōu)性能優(yōu)于IWOA 和CWOA。