基于模糊函數(shù)頻率軸投影的CDMA信號(hào)時(shí)/頻差估計(jì)

楊林森 張子敬 郭付陽(yáng)

(西安電子科技大學(xué) 西安710071)

0 引言

隨著衛(wèi)星通信技術(shù)的飛速發(fā)展，地面干擾輻射源對(duì)于衛(wèi)星正常通信的威脅越來(lái)越受到廣泛的關(guān)注，同步軌道雙星時(shí)差和頻差定位技術(shù)是一種常見(jiàn)的對(duì)地面輻射源定位的技術(shù)，已經(jīng)得到了成功的應(yīng)用［1］。

同步軌道雙星定位中，通常將與目標(biāo)輻射源直接通信的衛(wèi)星稱(chēng)為主星，在主星附近輔助定位的衛(wèi)星稱(chēng)為鄰星。由于天線(xiàn)輻射一般具有指向性，主星往往處于地面輻射源的波束主瓣方向內(nèi)，而鄰星常處于波束的某個(gè)旁瓣之內(nèi)。地面輻射源天線(xiàn)主瓣和旁瓣發(fā)出的同源信號(hào)有且僅有功率的差別，但是由于主星和鄰星的位置以及速度不同，地面輻射源發(fā)射信號(hào)經(jīng)過(guò)兩顆衛(wèi)星分別轉(zhuǎn)發(fā)后就會(huì)產(chǎn)生時(shí)差和頻差。雙星時(shí)/頻差定位原理正是通過(guò)測(cè)量雙星接收信號(hào)間的時(shí)差和頻差，再聯(lián)合地球表面方程實(shí)現(xiàn)對(duì)未知輻射源的定位［1-3］。

CDMA 信號(hào)是一種常用的通信信號(hào)，在CDMA信號(hào)通信體制內(nèi)，各用戶(hù)通過(guò)使用相互正交的用戶(hù)特征碼對(duì)數(shù)字通信內(nèi)容進(jìn)行調(diào)制，實(shí)現(xiàn)了多個(gè)用戶(hù)在相同頻段上的同時(shí)通信。由于CDMA 各用戶(hù)信號(hào)在時(shí)間和頻率上相互重疊，各用戶(hù)信號(hào)的信干比很低，導(dǎo)致各用戶(hù)時(shí)/頻差的精確估計(jì)具有很大困難。另外，由于CDMA 信號(hào)帶寬較大，相比于一般信號(hào)在相同的信號(hào)積累時(shí)間長(zhǎng)度內(nèi)，采樣的數(shù)據(jù)長(zhǎng)度明顯增加，導(dǎo)致運(yùn)算量顯著提高。模糊函數(shù)是時(shí)/頻差估計(jì)最常用的工具，然而直接計(jì)算CDMA 信號(hào)的模糊函數(shù)所需的運(yùn)算量過(guò)大。為解決CDMA 信號(hào)中各個(gè)用戶(hù)相互干擾的問(wèn)題，在文獻(xiàn)［4］中提出利用主星信號(hào)逆向分析得到的各用戶(hù)的重構(gòu)信號(hào)，用重構(gòu)信號(hào)進(jìn)行估計(jì)時(shí)/頻差。為了解決運(yùn)算量大的問(wèn)題，文獻(xiàn)［4］根據(jù)實(shí)際問(wèn)題中兩路信號(hào)頻差遠(yuǎn)小于采樣頻率的特點(diǎn)，對(duì)混合積信號(hào)進(jìn)行積分抽取，顯著提升了計(jì)算模糊函數(shù)的速度。

本文在文獻(xiàn)［4］提出的利用主星重構(gòu)信號(hào)進(jìn)行時(shí)/頻差估計(jì)的基礎(chǔ)上，提出一種基于模糊函數(shù)頻率軸投影的CDMA 信號(hào)時(shí)/頻差方法。該方法只需若干次FFT 即可實(shí)現(xiàn)，與現(xiàn)有方法對(duì)比，在運(yùn)算速度方面具有顯著優(yōu)勢(shì)。

1 信號(hào)模型

假設(shè)主星和鄰星接收到的CDMA 基帶信號(hào)為

其中sk(t)為主星接收到的第k個(gè)用戶(hù)信號(hào)，K為用戶(hù)個(gè)數(shù)，Δtk和Δfk為鄰星接收到相對(duì)于主星的第k個(gè)用戶(hù)信號(hào)的時(shí)差和多普勒頻差，n1(t)和n2(t)分別為主星和鄰星的噪聲。可以看出，鄰星接收到的信號(hào)為各用戶(hù)信號(hào)經(jīng)過(guò)不同時(shí)延和多普勒頻移后的線(xiàn)性疊加，且各用戶(hù)的時(shí)差和多普勒頻差不同。

對(duì)于主星接收信號(hào)中的第k個(gè)用戶(hù)信號(hào)sk(t)，有

其中dk(m)為用戶(hù)傳輸?shù)牡趍個(gè)符號(hào)，代表著用戶(hù)的通信內(nèi)容，M為用戶(hù)傳輸?shù)姆?hào)總數(shù)，Tb為單個(gè)符號(hào)持續(xù)的時(shí)間，CDMA 信號(hào)中各用戶(hù)通過(guò)特征波形區(qū)分，通信中傳輸?shù)拿恳晃环?hào)都需要用用戶(hù)特征波形調(diào)整，第k個(gè)用戶(hù)的特征波形為

其中bk(l)為特征波形的第l個(gè)碼片，Tc為單個(gè)碼片持續(xù)時(shí)間，且有Tb= LTc，p(t)為持續(xù)時(shí)間為T(mén)c的矩形脈沖。為了抑制通信中各用戶(hù)間的相互干擾，CDMA 信號(hào)常用正交偽隨機(jī)碼對(duì)各用戶(hù)的通信內(nèi)容進(jìn)行編碼。

2 模糊函數(shù)在頻率軸的投影

兩路信號(hào)f(t)和g(t)的互模糊函數(shù)(Ambiguity function，AF)定義為:

其中τ，v分別為時(shí)延和多普勒頻率。模糊函數(shù)可以認(rèn)為是信號(hào)在時(shí)域和頻域的二維相關(guān)，其峰值位置對(duì)應(yīng)兩路信號(hào)的時(shí)差和頻差。然而，當(dāng)信號(hào)較長(zhǎng)時(shí)，二維搜索模糊函數(shù)峰值的計(jì)算量過(guò)大，難以實(shí)時(shí)處理。

我們定義Pfg(v)為信號(hào)f(t)和g(t)的模糊函數(shù)χfg(τ，v)在頻率軸上的投影，文獻(xiàn)［5］給出了將模糊函數(shù)投影到頻率軸的快速算法:

其中F(v)和G(v)分別為信號(hào)f(t)和g(t)的傅立葉變換。式(5)將計(jì)算模糊函數(shù)在頻率軸上的投影簡(jiǎn)化為了兩路信號(hào)自相關(guān)乘積的傅立葉變換，因此，計(jì)算模糊函數(shù)在頻率軸上的投影只需要若干次FFT 即可實(shí)現(xiàn)。

3 CDMA 信號(hào)時(shí)/頻差估計(jì)

由于CDMA 信號(hào)中包含有多個(gè)用戶(hù)，且各用戶(hù)的時(shí)/頻差不同，直接對(duì)主星和鄰星接收信號(hào)做互模糊函數(shù)會(huì)產(chǎn)生多峰，當(dāng)用戶(hù)時(shí)/頻差接近時(shí)，峰值難以分辨，無(wú)法精確估計(jì)時(shí)/頻差。在文獻(xiàn)［4］中，通過(guò)對(duì)主星CDMA 信號(hào)逆向分析獲取各用戶(hù)的特征波形和通信內(nèi)容，將各用戶(hù)信號(hào)進(jìn)行重構(gòu)。若將重構(gòu)的主星各用戶(hù)信號(hào)分別與鄰星信號(hào)做互模糊函數(shù)，能夠有效去除用戶(hù)間的干擾，得到只有單峰的互模糊函數(shù)，最終能夠提高時(shí)/頻差估計(jì)精度。然而，利用重構(gòu)的主星各用戶(hù)信號(hào)計(jì)算完整的模糊函數(shù)，運(yùn)算量大的問(wèn)題仍然沒(méi)有解決。為了減少計(jì)算模糊函數(shù)的運(yùn)算量，文獻(xiàn)［4］利用現(xiàn)實(shí)中多普勒頻差遠(yuǎn)小于采樣頻率的特點(diǎn)，提出對(duì)混合積信號(hào)積分抽取的方法來(lái)實(shí)現(xiàn)CDMA 信號(hào)各用戶(hù)時(shí)/頻差快速估計(jì)，其實(shí)質(zhì)是通過(guò)對(duì)混合積信號(hào)降采樣來(lái)降低計(jì)算模糊函數(shù)的運(yùn)算量。

由于CDMA 信號(hào)中各用戶(hù)的特征碼采用相互正交的偽隨機(jī)編碼序列，因而CDMA 各用戶(hù)發(fā)射信號(hào)具有尖銳的自相關(guān)函數(shù)，而各用戶(hù)間的互相關(guān)函數(shù)值接近零。對(duì)于偽隨機(jī)序列，其模糊函數(shù)為圖釘型［6］，存在一個(gè)很尖的主瓣，主瓣峰值遠(yuǎn)高于旁瓣高度。若將主星重構(gòu)的各用戶(hù)信號(hào)與鄰星接收信號(hào)的互模糊函數(shù)在頻率軸上進(jìn)行投影將得到一個(gè)具有很高尖峰的曲線(xiàn)，投影的峰值的位置對(duì)應(yīng)原于模糊函數(shù)峰值的頻率坐標(biāo)，由投影的峰值位置可以估計(jì)出用戶(hù)雙星接收信號(hào)的頻差，對(duì)雙星信號(hào)補(bǔ)償頻差后通過(guò)一次相關(guān)，由相關(guān)峰位置即可估計(jì)出用戶(hù)雙星接收信號(hào)的時(shí)差。根據(jù)上面的分析，利用式(5)可以實(shí)現(xiàn)模糊函數(shù)在頻率軸上投影的快速計(jì)算，因此，本文提出了基于模糊函數(shù)在頻率軸上投影的CDMA 信號(hào)時(shí)/頻差估計(jì)快速方法，其具體流程為:

(1)，對(duì)于主星重構(gòu)的第k個(gè)用戶(hù)信號(hào)sk(t)和鄰星接收信號(hào)r2(t)，按式(5)計(jì)算互模糊函數(shù)在頻率軸上的投影Pskr2(v)。

(2)，搜索Pskr2(v)的峰值，記峰值位置為v0，則第k個(gè)用戶(hù)的多普勒頻差估計(jì)為。

步驟(4)，對(duì)重構(gòu)的各用戶(hù)信號(hào)依次按上述步驟處理，直到估計(jì)出所有用戶(hù)的時(shí)/頻差。

4 仿真實(shí)驗(yàn)

假設(shè)仿真的CDMA 通信系統(tǒng)中包含8 個(gè)用戶(hù)，各用戶(hù)使用長(zhǎng)度為31 的gold 偽隨機(jī)序列作為用戶(hù)特征碼。在T=5ms 的信號(hào)持續(xù)時(shí)間內(nèi)，8 個(gè)用戶(hù)同時(shí)進(jìn)行通信，數(shù)字通信符號(hào)傳輸速率為322.58kbps，單個(gè)符號(hào)的傳輸時(shí)間Tb=3.1μs，每傳輸一位符號(hào)都需要用31 位的用戶(hù)特征碼進(jìn)行調(diào)制，因此用戶(hù)特征碼碼片持續(xù)時(shí)間為T(mén)c=0.1μs，接收信號(hào)的采樣率為20MHz。雙星接收的8 用戶(hù)CDMA 信號(hào)中各用戶(hù)對(duì)應(yīng)的時(shí)延差為［0，0.2，1，-1，2，-4，-0.3，5］(μs)，多普勒頻差為［0，1，50，74，12，-4，6.8，-8］(KHz)，為了驗(yàn)證本方法的性能，特別選取第一個(gè)用戶(hù)和第二個(gè)用戶(hù)的時(shí)/頻差參數(shù)差距很小。

圖1 為主星和鄰星接收信號(hào)在信噪比為0dB 的情況下互模糊函數(shù)的局部仿真圖。從圖1 中可以明顯看出，CDMA 信號(hào)的模糊函數(shù)為圖釘型，存在若干個(gè)又高又窄的峰值，每一個(gè)峰值對(duì)應(yīng)于一個(gè)用戶(hù)的時(shí)/頻差。由于第一個(gè)用戶(hù)和第二個(gè)用戶(hù)時(shí)/頻差相距較近，對(duì)應(yīng)的模糊函數(shù)峰值部分重疊而無(wú)法清晰分辨，因此圖1 中的模糊函數(shù)只能看到7 個(gè)明顯峰值。

圖1 雙星接收信號(hào)的互模糊函數(shù)

為了避免用戶(hù)間的干擾，本文利用從主星接收信號(hào)中逆向分析得到的各用戶(hù)重構(gòu)信號(hào)進(jìn)行時(shí)/頻差估計(jì)。圖2 為第二個(gè)用戶(hù)的主星重構(gòu)信號(hào)與鄰星接收信號(hào)之間的局部互模糊函數(shù)。對(duì)比圖1 和圖2可以看出，由于主星重構(gòu)信號(hào)中只包含單個(gè)用戶(hù)的信號(hào)，利用主星重構(gòu)信號(hào)得到的模糊函數(shù)只存在單峰且模糊函數(shù)旁瓣更低，因此利用主星重構(gòu)信號(hào)能夠大大降低各用戶(hù)間的交叉干擾，進(jìn)而更準(zhǔn)確地估計(jì)出各個(gè)用戶(hù)的時(shí)/頻差。

圖2 用戶(hù)2 主星重構(gòu)信號(hào)與鄰星信號(hào)的互模糊函數(shù)

從圖1 和圖2 中可以看出，由于CDMA 信號(hào)中用戶(hù)特征碼具有偽隨機(jī)性，CDMA 信號(hào)的模糊函數(shù)為圖釘型，其具有又高又陡的主瓣且主瓣峰值遠(yuǎn)高于旁瓣高度，因此CDMA 信號(hào)模糊函數(shù)在頻率軸上投影的峰值點(diǎn)將對(duì)應(yīng)其模糊函數(shù)峰值點(diǎn)的頻率坐標(biāo)。圖3(a)為用戶(hù)2 主星重構(gòu)信號(hào)與鄰星接收信號(hào)間的模糊函數(shù)在頻率軸上投影的曲線(xiàn)。由圖3(a)可以看出CDMA 信號(hào)在頻率軸上的投影將產(chǎn)生一個(gè)明顯的峰值，根據(jù)圖3 中峰值的位置，可以估計(jì)出第二個(gè)用戶(hù)的頻差為0.999 kHz。對(duì)雙星接收信號(hào)補(bǔ)償頻差后再做相關(guān)，利用相關(guān)峰位置可以估計(jì)出該用戶(hù)的時(shí)差。圖3(b)為兩路信號(hào)補(bǔ)償頻率后相關(guān)運(yùn)算結(jié)果的仿真圖，由相關(guān)峰位置可以估計(jì)出雙星接收信號(hào)的時(shí)差為0.2μs。利用上述方法，可以分別估計(jì)出雙星接收的CDMA 信號(hào)中各用戶(hù)的時(shí)差為［0，0.2，1，-1，2，-4，-0.3，5］(μs)，頻差為［0，0.999，49.997，73.996，11.999，-4，6.8，-8］(kHz)。文獻(xiàn)［4］中利用各用戶(hù)主星重構(gòu)信號(hào)，提出了基于積分抽取的CDMA 時(shí)/頻差估計(jì)方法，而本文在其基礎(chǔ)上，利用模糊函數(shù)在頻率軸上的投影實(shí)現(xiàn)CDMA 信號(hào)時(shí)/頻差的快速估計(jì)。下面通過(guò)對(duì)比這兩種方法的復(fù)乘次數(shù)，說(shuō)明本方法在運(yùn)算速度上的優(yōu)勢(shì)。

圖3 CDMA 信號(hào)時(shí)/頻差估計(jì)仿真

對(duì)于文獻(xiàn)［4］中提出的積分抽取算法，假設(shè)雙星接收信號(hào)長(zhǎng)度為N且抽取倍數(shù)為D，則計(jì)算一次混合積信號(hào)需要N次復(fù)乘，抽取后混合積長(zhǎng)度降為N/D，對(duì)抽取的混合積做一次FFT 的運(yùn)算量為N D，計(jì)算模糊函數(shù)需要重復(fù)N次上述過(guò)程，則用文獻(xiàn)［4］中提出的算法估計(jì)時(shí)/頻差需要的運(yùn)算量約為。

對(duì)于本文提出的算法，同樣假設(shè)信號(hào)長(zhǎng)度為N，由式(5)可以看出，計(jì)算模糊函數(shù)在頻率軸上投影只需要計(jì)算5 次長(zhǎng)度為N的FFT，得到了頻差估計(jì)后，對(duì)信號(hào)補(bǔ)償頻率，最后通過(guò)相關(guān)實(shí)現(xiàn)時(shí)差估計(jì)，本文提出的算法需要的運(yùn)算量約為5Nlog2N+N+3(2N-1)log2(2N-1)。

對(duì)于上面仿真的參數(shù)，得到信號(hào)長(zhǎng)度N=100000，假設(shè)文獻(xiàn)［4］中的方法使用的抽取倍數(shù)為D=100 。表1 為本文提出的方法和文獻(xiàn)［4］中的方法在運(yùn)算量上的對(duì)比，可以看出，本文提出的算法在運(yùn)算量上具有顯著的優(yōu)勢(shì)，即使當(dāng)利用文獻(xiàn)［4］中提出的算法抽取倍數(shù)達(dá)到D=100 倍的情況下，本文算法在運(yùn)算量上仍然具有明顯的優(yōu)勢(shì)。

表1 積分抽取與本文方法運(yùn)算量對(duì)比

5 結(jié)論

針對(duì)雙星定位問(wèn)題中，現(xiàn)有CDMA 信號(hào)時(shí)/頻差估計(jì)方法運(yùn)算量較大的問(wèn)題，本文提出了基于模糊函數(shù)在頻率軸投影的CDMA 信號(hào)時(shí)/頻差估計(jì)方法。針對(duì)CDMA 信號(hào)模糊函數(shù)為圖釘型的特點(diǎn)，本文通過(guò)求模糊函數(shù)在頻率軸上的投影即可實(shí)現(xiàn)CDMA 信號(hào)間時(shí)/頻差的估計(jì)。仿真結(jié)果說(shuō)明本文提出的方法能夠快速估計(jì)CDMA 信號(hào)中各用戶(hù)的時(shí)/頻差，相比已有方法，運(yùn)算量顯著減少。

［1］Haworth D P，Smith N G，Bardelli R.Interference Iocalization for EUTELSAT satellites the first European transmitter location system［J］.International Journal of Satellite Communiactions，1997，15(3):155-183.

［2］龍寧，曹廣平，王勤果.雙星時(shí)差頻差定位系統(tǒng)中的多信號(hào)定位技術(shù)［J］.電訊技術(shù)，2011，51(2):16-20.

［3］嚴(yán)航，姚山峰.基于參考站的低軌雙星定位誤差校正分析［J］.計(jì)算機(jī)仿真，2012，29(10):114-118.

［4］嚴(yán)航，朱珍珍.基于積分抽取的時(shí)/頻差參數(shù)估計(jì)方法［J］.宇航學(xué)報(bào)，2013，34(1):99-105.

［5］Sharif M R，Abeysekera S S.Efficient wideband signal parameter estimation using a Radon-Ambiguity transform slice［J］.IEEE Trans.on AES，2007，43(2):673-688.

［6］Nadav Levanon，Eli Mozeson.Radar signals［M］.New Jersey:John Wiley and Sons，2004:100-164.