張?zhí)祢U 高永升 杜曉華 侯瑞玲

(重慶郵電大學(xué)信號與信息處理重慶市重點實驗室,重慶 400065)

1.引 言

直擴超寬帶(DS-UWB)信號采用極窄沖激脈沖和直接序列擴頻(DS-SS)技術(shù)相結(jié)合以獲得超寬頻譜,其具有非常低的平均發(fā)射功率,這較大地增強了其安全性、低截獲/檢測概率(LPI/D),現(xiàn)在已廣泛應(yīng)用于民用和軍事通信中。但非合作通信方要想對DS-UWB信號進(jìn)行盲檢測與估計,并獲取所傳輸?shù)男畔⒕蜁兊卯惓＠щy,最好的辦法就是先對DSUWB信號參數(shù)、擴頻序列進(jìn)行盲估計,再和主動通信一樣獲取所傳輸?shù)男畔ⅰ?/p>

目前,已報道的關(guān)于低信噪比DS-UWB信號盲檢測與估計方面的研究較少,而且大多都是針對單用戶DS-UWB信號。關(guān)于單用戶DS-UWB信號盲檢測的研究,文獻(xiàn)[1]給出了希爾伯特—黃法、隨機共振法、基于特征值分解法和分段相關(guān)累積(PACA)的檢測方法。文獻(xiàn)[2]利用循環(huán)平穩(wěn)性進(jìn)行了單用戶DS-UWB信號的盲檢測和參數(shù)估計。文獻(xiàn)[3]通過小波的分解重構(gòu)算法對UWB信號進(jìn)行了不同分辨率級別上的時延估計。但是,針對多用戶DSUWB信號的檢測與盲估計研究卻很少。

針對多用戶DS-UWB信號的碼片寬度估計問題,主要研究了多用戶DS-UWB信號的特點及其循環(huán)譜密度函數(shù),利用多用戶DS-UWB信號循環(huán)譜的譜冗余特性,完成多用戶DS-UWB信號的碼片寬度估計。

2.多用戶DS-UWB信號模型

所研究的超寬帶(UWB)信號是指超寬帶沖激無線電信號,以低占空比的超短沖激脈沖作為載體,采用DS-SS技術(shù)實現(xiàn)信息傳輸。首先,用PN序列對待發(fā)送的二進(jìn)制符號序列進(jìn)行擴頻調(diào)制;其次,用窄脈沖進(jìn)行脈沖幅度調(diào)制(PAM)。

在多用戶DS-UWB系統(tǒng)中,第k個用戶的信號波形可以寫為

q(t)的傅立葉變換為

式中,tau稱為脈沖形成因子,單位為秒,隨著脈沖形成因子增大,二階高斯脈沖波形變寬,但其能量譜密度的帶寬變窄,選擇符合超寬帶信號的二階高斯脈沖可通過選擇合適的脈沖形成因子獲得。注意到高斯脈沖的持續(xù)時間很長,這將不可避免地導(dǎo)致脈沖混疊和符號間干擾(ISI),因此,有必要對式(2)的高斯脈沖進(jìn)行有效截斷,截斷長度的上限由不超過切普(chip)持續(xù)時間Tc決定,下限受產(chǎn)生極窄脈沖的技術(shù)水平限制,這里采用已經(jīng)證明[Giancola 2002]的截斷長度Tu,令Tu=1.1tau表示對高斯脈沖的有效截斷,由此占空比可定義為Tu/Tc.

當(dāng)系統(tǒng)中存在K個用戶時,K個用戶共用空間頻譜資源,就產(chǎn)生了多用戶干擾(MUI)。所有的分析都是基于以下假設(shè)：

1)所有信源產(chǎn)生二進(jìn)制b序列,該序列由獨立同分布且0和1等概出現(xiàn)的隨機變量構(gòu)成。

2)所有信源脈沖重復(fù)頻率均為1/Tb.

3)信息碼互相獨立且等概率,每個碼隨機產(chǎn)生,且對應(yīng)一個PN序列。

4)無線電信號在無多徑影響的信道上傳輸,第n個用戶與其參考接收機之間的無多徑的信道沖激響應(yīng)是路徑增益Ak和傳輸時延τk的函數(shù)。傳輸時延是相互獨立的,且在[0,Tb]上服從均勻分布。沖激響應(yīng)表達(dá)式為 ：h(k)(t)=Akδ(t-τk),其中 δ(?)為Dirac沖激函數(shù)。

5)信道的輸出受到雙邊功率譜密度為N0/2(單位為W/Hz)的白噪聲n(t)的污染。

因此,接收到的K個用戶DS-UWB信號的波形為

當(dāng) τ1=τ2= … =τk時,式(4)表示同步多用戶DS-UWB系統(tǒng);反之,稱其為異步多用戶DS-UWB系統(tǒng)。

3.算法原理

3.1 信號的循環(huán)平穩(wěn)特性

周期平穩(wěn)信號是一種特殊的非平穩(wěn)信號,其具有周期平穩(wěn)性,即其統(tǒng)計參數(shù),如均值、相關(guān)函數(shù)等,都隨著時間作周期性的變化。對于這種信號更適合建模為循環(huán)平穩(wěn)隨機過程,循環(huán)平穩(wěn)隨機過程可用循環(huán)譜密度函數(shù)來表征,該理論被稱為循環(huán)譜理論或譜相關(guān)理論[5]。

假設(shè)x(t)是以T為周期的平穩(wěn)過程,現(xiàn)將其自相關(guān)函數(shù)展開成傅氏級數(shù),系數(shù)記為循環(huán)自相關(guān)函數(shù)(τ),得

再對(τ)做傅氏變換,得到x(t)的循環(huán)譜密度函數(shù)(也稱循環(huán)譜)為

式中：α稱為信號x(t)的循環(huán)頻率,α=m/T(m是所有整數(shù));f稱為x(t)頻譜頻率。式(5)中(τ)是x(t)的兩個復(fù)頻移函數(shù)u(t)和v(t)的互相關(guān)函數(shù),即

其中：

所以

式中：

式(11)中(t,f)稱為循環(huán)周期圖。式(12)中XT(t,f)稱為x(t)的短時傅立葉變換。從式(11)可知,循環(huán)譜密度函數(shù)(f)可表示為x(t)在頻率f+α/2與f-α/2處譜分量間的互相關(guān)密度,因此,循環(huán)譜密度函數(shù)(f)又稱為譜相關(guān)函數(shù)。

在循環(huán)相關(guān)檢測和循環(huán)平穩(wěn)信號的參數(shù)估計中,循環(huán)譜密度函數(shù)是一個重要的概念,能夠完全表征循環(huán)平穩(wěn)信號的很多特征。注意：實際上,循環(huán)平穩(wěn)信號x(t)在頻率f+α/2與f-α/2處譜分量是完全相關(guān)的,這種特殊的譜相關(guān)特性稱為譜冗余,一般的平穩(wěn)信號(包括高斯白噪聲)就沒有這種性質(zhì);利用譜冗余,可以從平穩(wěn)的干擾中將循環(huán)平穩(wěn)信號提取出來;不同的循環(huán)平穩(wěn)信號具有不一樣的循環(huán)頻率α,可以利用不同的循環(huán)譜密度函數(shù)檢測出不同的循環(huán)平穩(wěn)信號。也就是說：平穩(wěn)噪聲和干擾在α≠0處不呈現(xiàn)譜相關(guān),即(f)=0;而對于多用戶DS-UWB信號,有α=m/Tc(m∈Z),Tc就是多用戶DS-UWB信號的碼片寬度,因此,利用譜冗余就可以完成碼片寬度的估計。為方便計算,采用FFT Accumulation Method(FAM)算法[5]來對循環(huán)譜密度進(jìn)行估計。

3.2 循環(huán)自相關(guān)和循環(huán)譜的疊加特性

在一個比較寬的帶寬內(nèi),通常有許多信號和背景噪聲,設(shè)接收信號c(t)由信號rl(t)(l=1,2,…,L)和噪聲n(t)組成,即

假設(shè)L個信號和噪聲彼此統(tǒng)計獨立,那么其循環(huán)自相關(guān)和循環(huán)譜滿足疊加特性[7]

如果多個目標(biāo)信號的循環(huán)頻率彼此不同,背景噪聲不是循環(huán)平穩(wěn)的,就可以根據(jù)各自的循環(huán)頻率進(jìn)行檢測、信號分離、參數(shù)(載頻、符號率等)估計。特別地,如果只有一個信號,則c(t)=r(t)+n(t),其中n(t)是高斯白噪聲,那么(τ)=(τ),(f)=(f).

3.3 多用戶DS-UWB信號譜相關(guān)分析

為方便分析多用戶DS-UWB信號的循環(huán)譜,首先研究單用戶DS-UWB信號。由于噪聲不具有循環(huán)平穩(wěn)特性,為了簡化運算,以下推導(dǎo)暫不考慮噪聲的影響。單用戶DS-UWB信號可視為一個特殊的脈沖幅度調(diào)制(PAM)信號,其采用寬度遠(yuǎn)小于碼片(chip)寬度的二階高斯脈沖進(jìn)行調(diào)制。文獻(xiàn)[7]已經(jīng)給出了各種調(diào)制信號的循環(huán)譜密度函數(shù)。現(xiàn)在考慮單個用戶的情況,系統(tǒng)接收到的第k個信號波形為

常見的調(diào)制系統(tǒng)可被模型化為一個以平穩(wěn)過程為激勵的多輸入線性周期時變變換(LPTV)系統(tǒng)。所以,很容易用LPTV變化參數(shù)來研究信號的循環(huán)平穩(wěn)特性。循環(huán)譜密度為的輸入信號通過LPTV系統(tǒng)的輸出循環(huán)譜密度為[7]

如果LPTV是時不變系統(tǒng),上式簡化為

式中：α=m/T,T是LPTV系統(tǒng)的周期且T=Tc;Gn(f)和Gm(f)是沖激響應(yīng)為h(t,u)的傅立葉級數(shù)的系數(shù)gn(τ)和gm(τ)的傅立葉變換。

首先考慮一般的幅度調(diào)制(AM),傳統(tǒng)AM信號的模型為

其中p(t)為載波,其傅立葉級數(shù)為

所以式(19)就可看成一個輸入為a(t)的LPTV變換,其沖激響應(yīng)的傅立葉級數(shù)系數(shù)為

式中,δ(τ)是Dirac沖激函數(shù),且

直接應(yīng)用式(17)可以得到循環(huán)譜密度函數(shù)

當(dāng)

對所有的β=m/Tc,有pβ=1/Tc;利用式(23)得到循環(huán)譜密度函數(shù)為

如果時間序列x(t)通過一個沖激響應(yīng)為q(t)的濾波器,則y(t)=x(t)?q(t),y(t)就是PAM信號。將替代式(19)中的a(t),則式(16)可看作x(t)通過一個沖激響應(yīng)為q(t)的LPTV系統(tǒng)的處理過程,根據(jù)式(17)就得到系統(tǒng)接收到的第k個用戶信號的循環(huán)譜密度函數(shù)

式中：

如果把信息序列和PN碼序列調(diào)制的序列看作具有一般平穩(wěn)特性[8],則有

根據(jù)循環(huán)譜的疊加特性,多用戶DS-UWB信號的循環(huán)譜密度函數(shù)為

從式(30)可以看出：理論上多用戶DS-UWB信號只是在單用戶的基礎(chǔ)上改變了譜相關(guān)函數(shù)的相位和幅度信息,并沒有改變信號的循環(huán)頻率和循環(huán)譜密度函數(shù)的譜相關(guān)特性。多用戶DS-UWB信號(同步或者異步)仍然具有循環(huán)平穩(wěn)性,循環(huán)頻率就是碼片寬度Tc的倒數(shù)的整數(shù)倍。多用戶DS-UWB信號的循環(huán)譜密度將在α=m/Tc處產(chǎn)生譜線。因此,在某一頻譜頻率f處沿α軸搜索譜線間隔,碼片寬度就能準(zhǔn)確估計出來。

4.仿真實驗

仿真實驗1：

假設(shè)系統(tǒng)中有5個用戶,并具有理想功率控制,且系統(tǒng)中各個用戶的傳輸功率相同,Pow=-30 d Bm,采樣頻率fc=5×1010Hz,信息符號numbits=10 bit,碼片寬度Tc=2.0×10-9s,脈沖成形因子tau=0.25×10-9s,每比特映射的脈沖數(shù)Ns=63位,PN碼長Np=63位,高斯單周期脈沖沖擊響應(yīng)持續(xù)時間Tm=0.5×10-9s.加入零均值高斯白噪聲,信噪比SNR=-10 d B.利用FAM算法計算多用戶DS-UWB信號的循環(huán)譜密度,采樣點數(shù)N′=12755,1/Δf的數(shù)據(jù)段內(nèi)的離散傅立葉變換的點數(shù)N=32,Δt=12755Ts.通過計算機仿真可得到5個用戶DS-UWB信號(同步和異步)的循環(huán)譜密度函數(shù),其“循環(huán)頻率α-頻譜頻率f-循環(huán)密度函數(shù)幅值”的局部放大圖如圖1和圖2所示。

由于PAM調(diào)制脈沖q(t)不同于傳統(tǒng)的矩形脈沖,所以多用戶DS-UWB信號的循環(huán)譜密度函數(shù)與傳統(tǒng)的PAM信號有所不同。從仿真結(jié)果可以看出：在f=0兩邊,多用戶DS-UWB信號的循環(huán)譜密度函數(shù)出現(xiàn)兩個對稱的峰值,而不是在f=0處有單個峰值。

仿真實驗2：

由式(30)可看出,多用戶DS-UWB信號具有循環(huán)平穩(wěn)性,其循環(huán)譜密度將在α=m/Tc處產(chǎn)生譜線。從圖1和圖2中發(fā)現(xiàn)：在某個頻譜頻率f的截面上,會出現(xiàn)一個幅度最大的譜線,在產(chǎn)生最大譜線的頻譜頻率fmax附近搜索次大譜線,兩譜線之間的寬度就是估計的碼片速率。圖3和圖4是fmax=7.815×109Hz的截面局部放大圖。在仿真實驗中,頻譜頻率fmax截面圖只是截取了循環(huán)頻率α=0正負(fù)5倍1/Tc的寬度,并且對所有的坐標(biāo)進(jìn)行了歸一化處理,歸一化頻率f′max=0.1379。從圖3和圖4中可以看到,由于多用戶之間干擾的存在,使得α≠k/Tc時出現(xiàn)類似噪聲的干擾譜線,但是并沒有改變信號的循環(huán)平穩(wěn)性。比較圖3和圖4,同步情況下截面圖中的譜線幅度大于異步的情況,說明異步多用戶DS-UWB信號在α=m/Tc處的譜線更容易被噪聲淹沒。

仿真實驗中,經(jīng)過搜索可得譜線之間的間隔ΔN=127,Δα=1/12755Ts=3.92 ×106Hz,碼片速率1/Tc=ΔN/12755×Ts=4.9784×108bit/s,和真實值相差2.16×106bit/s,差值小于 Δα;估計值和真值之比為0.9957;估計的碼片寬度T′c=2.0087×10-9s.由此可見,用此方法可以很好地估計多用戶DS-UWB信號的碼片寬度。

仿真實驗3：

為進(jìn)一步驗證本算法的有效性及信噪比對估計性能的影響,在仿真實驗2的基礎(chǔ)上做M=400次的MonteCarlo仿真,得到(同步/異步情況下)要準(zhǔn)確估計偽碼寬度所需的累加次數(shù)隨信噪比變化的曲線,如圖5所示。

圖5 兩種算法下信噪比和平均累加次數(shù)的關(guān)系

從圖5可以看出,同步多用戶DS-UWB信號在-24 dB的信噪比下仍可準(zhǔn)確估計碼片寬度,異步多用戶DS-UWB信號在相同條件下跟同步情況下的差別并不明顯,同時所得到的最低信噪比容限比文獻(xiàn)[2]中所得的單用戶結(jié)果要小。同時,在同一信噪比下,異步多用戶DS-UWB信號比同步多用戶DS-UWB信號需要稍微多一點的累加次數(shù)才能準(zhǔn)確地估計偽碼寬度。由此可見,用此方法可在較低的信噪比下準(zhǔn)確地估計多用戶DS-UWB信號的碼片寬度,且同步多用戶DS-UWB信號的估計性能要略優(yōu)于異步多用戶DS-UWB信號。

另外,圖5表明在進(jìn)行多用戶DS-UWB信號碼片寬度估計時,循環(huán)譜算法所能達(dá)到的信噪比容限高于分段相關(guān)累積法所能達(dá)到的信噪比容限,估計性能比較好。但是在程序運行中,循環(huán)譜算法估計參數(shù)時需要更多的時間來完成,計算量比較大。

5.結(jié) 論

針對多用戶DS-UWB信號模型,通過理論分析得到了多用戶DS-UWB信號的循環(huán)譜密度函數(shù)。理論推導(dǎo)發(fā)現(xiàn)系統(tǒng)中改變用戶數(shù)目的情況下,接收信號的循環(huán)譜密度函數(shù)僅僅在相位和幅度信息上改變,但信號的循環(huán)頻率并不改變。所以,利用循環(huán)相關(guān)理論也可以完成多用戶DS-UWB信號檢測及其碼片寬度估計。最后給出了仿真結(jié)果,驗證了理論分析的結(jié)果且和分段相關(guān)累積法進(jìn)行比較。從仿真結(jié)果看,用循環(huán)譜分析法估計多用戶DS-UWB信號的碼片寬度,可以在較低的信噪比下得到較好的估計效果,并且同步多用戶DS-UWB信號的估計性能要略優(yōu)于異步多用戶DS-UWB信號;此算法的估計性能明顯比分段相關(guān)累積法的估計性能好。

[1] 王明陽.非合作超寬帶沖激無線電信號檢測技術(shù)研究[D].2006,國防科學(xué)技術(shù)大學(xué)：信息與通信工程.

[2] 虞路勤,張?zhí)祢U,李越雷,等.DS-UWB信號的譜相關(guān)分析及碼片寬度估計[J].電視技術(shù),2010,34(4)：56-59.YU Luqin,ZHANG Tianqi,LI Yuelei,et al.The analysis of spectral correlation of DS-UWB signal and estimation of chip interval[J].Video Engineering,2010,34(4)：56-59.(in Chinese)

[3] 陳雋永,徐 慶,徐繼麟.超寬帶雷達(dá)中基于小波變化的時延估計方法[J].電波科學(xué)學(xué)報,2001,16(2)：245-248 CHEN Juanyong,XU Qing,XU Jilin.A new timedelay estimation for the ultra wide-band radar[J].Chinese Journal of Radio Science,2001,16(2)：245-248.(in Chinese)

[4] 張?zhí)祢U,苗 圃,馬國寧,等.多徑環(huán)境下的直擴信號偽碼周期估計[J].電波科學(xué)學(xué)報,2009,24(5)：973-978.ZHANG Tianqi,MIAO Pu,MA Guoning,et al.Period estimation of the PN sequence for direct sequence spread spectrum in multipath environment[J].Chinese Journal of Radio Science,2009,24(5)：973-978.(in Chinese)

[5] GARDNER W A.Spectral correlation of modulated signals：part I-analog modulation[J].IEEE Transactions on Communication,1987,35(6)：584-594.

[6] 喻 凱.基于循環(huán)譜的MPSK信號盲參數(shù)估計[D].2007,電子科技大學(xué)：通信與信息工程.

[7] ROBERTS RS,BROWN W A,LOOMISH H,et al.Computationally efficient algorithms for cyclic spectral analysis[J].IEEE Signal Processing Magazine,1991,8(2)：38-49.

[8] 張?zhí)祢U.DS及DS/CDMA信號的估計方法研究[D].2002,電子科技大學(xué)：電路與系統(tǒng).

[9] 張?zhí)祢U,周正中.直擴信號偽碼周期的譜檢測[J].電波科學(xué)學(xué)報,2001,16(4)：518-521.ZHANG Tianqi,ZHOU Zhengzhong.A new spectral method of periodic detection of PC sequence in lower SNRDS/SS signals[J].Chinese Journal of Radio Science,2001,16(4)：518-521.(in Chinese)

[10] GA RDNERW A.Exploitation of spestral redundancy in cyclostationary signals[J].IEEE Signal Processing Magazine.1991,8(2)：14-36.

[11] GARDNER W A.Measurement of spestral correlation[J].IEEE Transactions on Acoustics,Speech and Signal Processing,1986,34(5)：1111-1123.