梁 芳， 趙林軍

(1.陜西理工學(xué)院 物理與電信工程學(xué)院， 陜西 漢中 723000;2.陜西理工學(xué)院 電氣工程學(xué)院， 陜西 漢中 723000)

0 引 言

直序列擴(kuò)頻系統(tǒng)中常使用延遲抖動(dòng)環(huán)(Delay-Locked Loop，DLL)實(shí)現(xiàn)收發(fā)偽碼相位同步[1]。根據(jù)偽碼長度與采樣頻率的關(guān)系，DLL又分為等量采樣環(huán)路(Commensurate Sampling Loop，CSL)與非等量采樣環(huán)路(Non-Commensurate Sampling Loop，NCSL)[2]。由于NCSL可在較低采樣條件下，實(shí)現(xiàn)較高的時(shí)間分別率而得到廣泛應(yīng)用。

2002年，QUIRK等[2]提出了NCS技術(shù)偽碼跟蹤環(huán)，并指出當(dāng)偽碼每chip的時(shí)寬Tc與采樣間隔Ts滿足kTs=nTc時(shí)，跟蹤環(huán)路可獲得Tc/k的跟蹤精度，而沒有對(duì)采樣頻率與偽碼環(huán)路相關(guān)函數(shù)的影響進(jìn)行分析。2006年，QUIRK等[3]對(duì)NCSL的累加時(shí)間進(jìn)行了分析，提出環(huán)路的累積時(shí)長應(yīng)是未采樣偽碼周期與其一個(gè)整周期內(nèi)非等量采樣樣點(diǎn)數(shù)的最小公倍數(shù)。這對(duì)實(shí)時(shí)性要求較高的系統(tǒng)而言，工程上難以實(shí)現(xiàn)。2010年，柯颋等[4]論述了多普勒頻移可顯著改善等量采樣偽碼跟蹤環(huán)的時(shí)間分辨力，同時(shí)指出偽碼跟蹤環(huán)路的時(shí)間鑒別力與采樣頻率之間不存在簡潔的邏輯關(guān)系。因此，在偽碼跟蹤環(huán)路的數(shù)字化實(shí)現(xiàn)時(shí)，往往需要花費(fèi)大量時(shí)間對(duì)其相關(guān)參數(shù)進(jìn)行調(diào)試。2014年，JIN等[5]提出當(dāng)η=Ts/Tc=I.d1d2…dv，累積相關(guān)時(shí)長L=10v(L是整數(shù)倍的偽碼周期長)，且η·10v與10v互素時(shí)，DLL中的收發(fā)偽碼互相關(guān)函數(shù)具有理想包絡(luò)。

綜上所述，當(dāng)DLL環(huán)路中的偽碼碼片速率給定，由于系統(tǒng)中多種因素的影響，實(shí)際的NCSL采樣頻率不能滿足文獻(xiàn)[4]的約束條件，因此，按其難以評(píng)估采樣NCS技術(shù)的偽碼跟蹤環(huán)路的性能，且文獻(xiàn)[2-5]中給出的相關(guān)約束條件本身就難以實(shí)現(xiàn)。故研究NCSL偽碼跟蹤環(huán)的MATLAB建模與分析，對(duì)設(shè)計(jì)者快速評(píng)判NCSL環(huán)路具有實(shí)際工程意義。

1 偽碼DLL跟蹤環(huán)路的結(jié)構(gòu)及其參數(shù)設(shè)計(jì)

由文獻(xiàn)[3,6]知，目前在直擴(kuò)通信系統(tǒng)中常使用DLL跟蹤環(huán)實(shí)現(xiàn)發(fā)信號(hào)偽碼與本地再生偽碼相位的同步，該環(huán)路的結(jié)構(gòu)如圖1所示。其中r(t)為接收到基帶直擴(kuò)信號(hào)，r(k)為采樣量化器(ADC)的輸出。為了實(shí)現(xiàn)偽碼相位的精確同步，三路相關(guān)累積分別完成r(k)與本地再生偽碼的超前相位PN(k+1)、當(dāng)前相位PN(k)和滯后相位PN(k-1)在一個(gè)或數(shù)個(gè)偽碼周期內(nèi)的相乘累加運(yùn)算[7]，再按照相應(yīng)的鑒相算法由上述三路相關(guān)累積輸出量給出本地偽碼相位與接收信號(hào)中的偽碼相位的誤差e(k)，最后由e(k)控制本地再生偽碼發(fā)生器更新一組新的偽碼相位輸出PN(i)，(i=k+1，k，k-1)。其中，環(huán)路濾波器F(z)可使得本地偽碼相位平穩(wěn)轉(zhuǎn)移。NCO為本地偽碼發(fā)生器提供工作時(shí)鐘，在文中所述的仿真分析中，NCO由一個(gè)門限ε代替。當(dāng)y(k)≥ε，本地偽碼再生器相位向左移動(dòng)更新；當(dāng)y(k)≤-ε時(shí)，本地偽碼再生器相位向右移動(dòng)更新；當(dāng)y(k)∈(-ε，ε)時(shí)，本地偽碼再生器相位不更新。

圖1 基帶偽碼跟蹤環(huán)原理

現(xiàn)就偽碼DLL跟蹤環(huán)中的相關(guān)算法簡述如下：設(shè)跟蹤環(huán)路的離散信號(hào)為

(1)

式中S為信號(hào)功率，d為符號(hào)，每符號(hào)1個(gè)完整的Gold碼長，k為采樣序號(hào)，kτ為信號(hào)傳輸時(shí)延τ經(jīng)離散后的值。DLL環(huán)路中的超前相關(guān)、當(dāng)前相關(guān)與滯后相關(guān)分別為

(2)

式中累積時(shí)長L=nfsTb，其中n為正整數(shù)，Tb為符號(hào)d的寬度，即累積可在一個(gè)符號(hào)內(nèi)或多個(gè)符號(hào)內(nèi)進(jìn)行。根據(jù)文獻(xiàn)[8]，環(huán)路鑒相算法采用歸一化點(diǎn)積功率型算法[9]，于是，鑒相器的輸出誤差e(k)為

(3)

環(huán)路中數(shù)字濾波F(z)采用理想二階數(shù)字濾波器，其傳遞函數(shù)為

(4)

濾波器F(z)系數(shù)C1與C2分別為

C1=4π·ξ·Bn·Tb,

C2=(2π·Bn·Tb)2,

(5)

式中ξ=0.707為阻尼系數(shù)，Bn為環(huán)路噪聲帶寬。環(huán)路濾波器的誤差輸出量為：

y(k)=(C1+C2)e(k)-C1e(k-1)+y(k-1)。

(6)

2 偽碼DLL環(huán)的非等量采樣

按照文獻(xiàn)[4]所述，圖1偽碼跟蹤環(huán)的實(shí)現(xiàn)有等量采樣與非等量采樣兩種技術(shù)。其中，采用等量采樣技術(shù)時(shí)，易引起時(shí)間分辨模糊現(xiàn)象，如圖2所示。其中的連續(xù)方波為實(shí)際信號(hào)的波形，時(shí)間軸上的“點(diǎn)”代表采樣位置，設(shè)偽碼每chip采2個(gè)樣點(diǎn)(Ts=1/fs)。由于采樣脈沖起始位置的不確定性，使得T1時(shí)刻DLL環(huán)路中本地再生偽碼采樣(如圖2(a)所示)與接收信號(hào)中的偽碼采樣(如圖2(b)所示)和T2時(shí)刻DLL環(huán)路中本地再生偽碼采樣與接收信號(hào)中偽碼采樣(如圖2(c)所示)相關(guān)時(shí)，結(jié)果將是完全一樣。而實(shí)際對(duì)應(yīng)的連續(xù)信號(hào)的相關(guān)顯然是存在差異。于是，等量采樣的偽碼跟蹤環(huán)就產(chǎn)生了時(shí)間跟蹤的模糊問題。

圖2 等量采樣時(shí)間分辨模糊示意圖

為此，QUIRK等[2-3]提出了非等量采樣技術(shù)，如圖3所示。圖中每碼片樣點(diǎn)數(shù)目不同，T1時(shí)刻偽碼跟蹤環(huán)路中收發(fā)偽碼互相關(guān)峰顯然是不同于T2時(shí)刻偽碼跟蹤環(huán)路中的偽碼互相關(guān)峰的，因此，二者是可分辨的。故采用NCS技術(shù)后，DLL環(huán)路可以區(qū)分圖2所示的時(shí)延τ1與時(shí)延τ2。

圖3 非等量采樣技術(shù)示意圖

3 仿真流程與代碼實(shí)現(xiàn)

基于MATLAB實(shí)現(xiàn)NCSL偽碼DLL環(huán)路的建模流程如圖4所示。其中，fso是基帶采樣頻率的中心頻率，Δ為采樣頻率變化區(qū)間最大距離，仿真采樣速率變化范圍為[fso-Δ/2，fso+Δ/2](單位為MHz)。v為采樣頻率變化步長。偽碼DLL跟蹤環(huán)路假設(shè)已實(shí)現(xiàn)±1/2 chip捕獲。為此程序設(shè)計(jì)每次相關(guān)后，本地偽碼的相位移動(dòng)1 Sample。為了降低仿真實(shí)現(xiàn)難度，圖1中的信號(hào)樣點(diǎn)量化采用浮點(diǎn)數(shù)，不考慮量化誤差與運(yùn)算誤差。

圖4 偽碼DLL跟蹤 環(huán)建模流程圖

從第二小節(jié)知，等量采樣偽碼DLL環(huán)路與非等量采樣偽碼DLL環(huán)路的最大區(qū)別是偽碼每碼片的采樣數(shù)不同。且由圖3可知，在采樣頻率給定的條件下，實(shí)際系統(tǒng)中信號(hào)的首樣點(diǎn)落在[0，Ts]區(qū)間具有不確定性。因此，基于MATLAB實(shí)現(xiàn)非等量采樣偽碼DLL環(huán)路的關(guān)鍵是實(shí)現(xiàn)對(duì)r(t)與本地再生偽碼PN(·)的非等量采樣。為了偽碼便于實(shí)現(xiàn)，r(t)中的每個(gè)符號(hào)(待傳數(shù)據(jù))采用一個(gè)完整的偽碼周期進(jìn)行調(diào)制。故對(duì)r(t)與本地再生偽碼PN(·)的采樣具有相似性。

下文著重介紹偽碼的非等量采樣的實(shí)現(xiàn)。為論述方便，首先給出偽碼的非等量采樣函數(shù)的MATLAB代碼，具體如下：

function [out]=NCS(fs，fc，n，seq)； %n為1采樣間隔等分?jǐn)?shù)

Ts=1/fs； %采樣間隔

Tc=1/fc； %偽碼每碼片時(shí)寬

len=length(seq)； %偽碼長

Sample_Delay=n-1； %設(shè)置首樣隨機(jī)時(shí)延

Start=randint(1，1，[1，Sample_Delay])；

kk=Start/(n*fs)； %計(jì)算首樣時(shí)延時(shí)間

i=1；

j=0；

while j

while kk

temp=mod(j，1023)+1；

out(i)=seq(temp)；

i=i+1；

kk=kk+Ts；

end

kk=kk-Tc；

j=j+1；

end

end

在NCS子函數(shù)中，seq是由為1 023 chip的Gold序列；Sample_Delay用于設(shè)置偽碼離散化首樣點(diǎn)時(shí)延范圍，該時(shí)延范圍的最大值為n，再由MATLAB系統(tǒng)自帶函數(shù)randint()在[0，Sample_Delay]范圍內(nèi)產(chǎn)生隨機(jī)時(shí)延，以模擬實(shí)際系統(tǒng)中的采樣脈沖的不確定性。這是因?yàn)橛?jì)算機(jī)內(nèi)的數(shù)據(jù)均是離散的，而時(shí)間系統(tǒng)中，在[0，Ts]內(nèi)有無窮個(gè)時(shí)間點(diǎn)，為了簡化模擬，代碼中預(yù)設(shè)了一個(gè)變量n為[0，Ts]內(nèi)的預(yù)設(shè)等分?jǐn)?shù)。如此以來，結(jié)合雙while嵌套語句，最終完成偽碼或信號(hào)r(t)的非等量采樣。

4 仿真結(jié)果與分析

以下是NCSL的仿真分析。

實(shí)驗(yàn)1 Gold碼非等量采樣方法實(shí)現(xiàn)。為便于觀察，實(shí)驗(yàn)參數(shù)設(shè)為：fc=1 023，fs=2 240，delay=9，代碼運(yùn)行后，Start=7。在上述參數(shù)下，對(duì)應(yīng)離散后的結(jié)果，可視為[0，Ts]內(nèi)又分了10等分非等量采樣，從第7個(gè)點(diǎn)開始，以Ts為間隔分布。表1列出了NCS函數(shù)運(yùn)行結(jié)果，其中所用符號(hào)與第3節(jié)NCS()函數(shù)中變量相對(duì)應(yīng)。

表1 偽碼NCS采樣結(jié)果表

結(jié)果分析：由表1可知，在實(shí)驗(yàn)1參數(shù)條件下，文中所述的非等量采樣代碼運(yùn)行后，在Gold碼的第4個(gè)碼片與第9個(gè)碼片內(nèi)采得3個(gè)樣點(diǎn)，而在第1、2、3、5、6、7、8、10碼片內(nèi)各采得2個(gè)樣點(diǎn)，可見文中所述的分等量采樣實(shí)現(xiàn)方法是正確的。

實(shí)驗(yàn)2采樣頻率對(duì)非等量采樣DLL環(huán)路偽碼相關(guān)函數(shù)的影響分析。選擇采樣頻率，應(yīng)使得環(huán)路中的收發(fā)偽碼互相關(guān)函數(shù)R(τ)具有良好的對(duì)稱性，這是影響環(huán)路時(shí)延分辨力的關(guān)鍵。由于文獻(xiàn)[2-5]給出的約束關(guān)系在實(shí)際系統(tǒng)中難以長期滿足，故需要通過大量實(shí)驗(yàn)對(duì)偽碼跟蹤環(huán)所使用的采樣頻率、相關(guān)累積時(shí)間進(jìn)行評(píng)估。仿真參數(shù)設(shè)置為：BPSK偽碼采用1 023平衡Gold碼調(diào)制的基帶信號(hào)，符號(hào)噪聲比均為-10 dB。其中，偽碼速率為1.023 Mchip/s，采樣速率為(4.092+Δ) MHz/s，(Δ∈[-0.005，0.005])，頻率變化步長為1 kHz。偽碼DLL環(huán)中的累積時(shí)間L為一個(gè)完整的偽碼周期(具體點(diǎn)數(shù)由程序根據(jù)偽碼周期與所使用的采樣頻率自動(dòng)計(jì)算獲得)。圖5畫出了DLL中偽碼互相關(guān)函數(shù)的5條曲線，為了便于觀察，它們分別取其各自互相關(guān)函數(shù)峰值左右各3個(gè)點(diǎn)所繪制。由圖5可以看出：(1)fs∈[4.089，4.091 9] MHz時(shí)，R(τ)在峰值左右均具有良好的對(duì)稱性；(2)fs≤4.087 MHz時(shí)，R(τ)的右斜邊下降較快；(3)fs>4.092 MHz時(shí)，R(τ)的左斜邊下降較快。

結(jié)果分析：若假設(shè)實(shí)際的非等量采樣偽碼DLL環(huán)路使用的采樣頻率的中心為4.090 MHz，信號(hào)r(t)中的多普勒頻移或采樣時(shí)鐘在(4.090±0.001) MHz內(nèi)變化時(shí)，環(huán)路中的偽碼相關(guān)函數(shù)仍具有良好的對(duì)稱性，即環(huán)路的時(shí)間分辨力仍保持良好。這一結(jié)論在更高采樣頻率時(shí)，具有類似情況出現(xiàn)。同時(shí)，采樣頻率對(duì)偽碼DLL中的偽碼互相關(guān)函數(shù)的峰值有一定影響，但其影響在該頻率窗內(nèi)大約為0.32 dB以內(nèi)?？紤]到這一因素，故非等量采樣偽碼跟蹤環(huán)所使用的采樣頻率小于4倍過采樣時(shí)，其捕獲門限的選擇應(yīng)以4倍等量采樣門限為準(zhǔn)(由圖5看出，fs=4fc=4.092 MHz時(shí)，其偽碼互相關(guān)峰應(yīng)小于fs=4.091 9 MHz的偽碼互相關(guān)峰)。

同時(shí)，由圖5又可以得出：若在4倍偽碼速率附近實(shí)現(xiàn)非等量采樣偽碼跟蹤環(huán)路時(shí)，應(yīng)該選擇稍小于4倍的偽碼速率作為系統(tǒng)的采樣頻率，而不能選取稍高于4倍的偽碼速率為系統(tǒng)的采樣時(shí)鐘，這是由于采用前者可獲得偽碼捕獲門限的下界。而選擇后者，顯然難以確定這一門限值。

實(shí)驗(yàn)3 非等量采樣偽碼DLL環(huán)路濾波器響應(yīng)分析。仿真參數(shù)設(shè)置同實(shí)驗(yàn)2，環(huán)路噪聲帶寬Bn=2 Hz。圖6畫出了不同采樣頻率下，非等量采樣偽碼DLL環(huán)路濾波器輸出曲線。

圖5 NCSL環(huán)路自相關(guān)函數(shù)曲線 圖6 偽碼跟蹤環(huán)環(huán)路濾波器輸出曲線

結(jié)果分析：(1)按照實(shí)驗(yàn)2條件，若偽碼跟蹤環(huán)路采用的采樣頻率fs<4fc時(shí)，環(huán)路可迅速收斂；當(dāng)采樣頻率fs>4fc時(shí)，環(huán)路收斂速度較慢；(2)按照實(shí)驗(yàn)2與實(shí)驗(yàn)3的結(jié)論(1)選取恰當(dāng)?shù)牟蓸宇l率后，與fs=4fc的等量采樣偽碼跟蹤環(huán)路的收斂性相比較，二者無明顯差別。

5 結(jié)束語

在直接序列擴(kuò)頻系統(tǒng)中，由于信號(hào)中存在多普勒頻移以及本地時(shí)鐘的抖動(dòng)等因素的制約，實(shí)際的偽碼跟蹤環(huán)路的參數(shù)設(shè)置難以滿足相關(guān)文獻(xiàn)報(bào)道的約束條件，為此，筆者采用MATLAB對(duì)非等量采樣偽碼DLL環(huán)進(jìn)行了建模仿真分析。仿真結(jié)果表明：偽碼速率給定時(shí)，在近4倍(小于4倍)的過采樣偽碼跟蹤環(huán)中，存在一個(gè)大約2 kHz采樣頻率窗口，使得系統(tǒng)中所使用的采樣頻率落在該窗口內(nèi)時(shí)，偽碼DLL環(huán)路仍具有良好的時(shí)間分辨力。同時(shí)也發(fā)現(xiàn)，在相同條件下，當(dāng)采樣頻率小于該頻率窗的下限時(shí)，偽碼DLL環(huán)中的偽碼互相關(guān)函數(shù)右斜邊下降快于其左斜邊；當(dāng)采樣頻率大于該頻率窗的上限時(shí)，偽碼DLL環(huán)中的偽碼互相關(guān)函數(shù)左斜邊下降快于其右斜邊。同時(shí)，當(dāng)fs4fc(“”表示稍小)時(shí)，NCSL與CSL無明顯差別。文中所述的結(jié)論顯然對(duì)非等量采樣環(huán)路的硬件實(shí)現(xiàn)具有工程意義。

文中主要以4倍偽碼速率采樣時(shí)鐘論述，但是其結(jié)論對(duì)更高的采樣頻率有類似結(jié)論。讀者在具體設(shè)計(jì)實(shí)現(xiàn)中，可依據(jù)相關(guān)參數(shù)，按照文中方法，進(jìn)行環(huán)路參數(shù)分析，以確保偽碼跟蹤環(huán)路能良好工作。

[參考文獻(xiàn)]

[1] MILEANT A，MILLION S，HINEDI S，et al. The performance of the all-digital data transition tracking loop using nonlinear analysis[J]. Communications，IEEE Transactions on，1995，43(234)：1202-1215.

[2] QUIRK K J，SRINIVASAN M. Analysis of sampling and quantization effects on the performance of PN code tracking loops[C]//Communications，ICC 2002. IEEE International，2002：1480-1484.

[3] QUIRK K J，SRINIVASAN M. PN code tracking using noncommensurate sampling[J]. Communications， IEEE Transactions on，2006，54(7)：1349-1349.

[4] 柯颋，胡修林，劉禹圻，等．GNSS接收機(jī)中采樣率對(duì)時(shí)間鑒別力的影響[J]．宇航學(xué)報(bào)，2010，31(12)：2722-2729．

[5] JIN Xiao-jun，XU Zhao-bin，ZHANG Chao-jie，et al. Simple approach to determining parameters of noncommensurate sampling for optimal pseudo-noise code phase delay discrimination[J]. Electronics Letters，2014，50(4)：283-284.

[6] 張大衛(wèi)，胡修林，李晨．衛(wèi)星定位接收機(jī)高精度偽碼跟蹤的設(shè)計(jì)與實(shí)現(xiàn)[J]．電子技術(shù)應(yīng)用，2009，35(9)：108-111．

[7] SUWON K，YONG-HWAN L. Rapid acquisition of PN signals for DS/SS systems using a phase estimator[J]. Selected Areas in Communications，IEEE Journal on，2001，19(6)：1128-1137.

[8] 趙林軍．等量采樣偽碼跟蹤環(huán)的MATLAB實(shí)現(xiàn)[J]．陜西理工學(xué)院學(xué)報(bào)：自然科學(xué)版，2012，28(6)：15-18．

[9] 黃波，胡修林．一種高精度系統(tǒng)時(shí)間同步方法[J]．信息與控制，2011，40(3)：347-351．