王曉君,尚 燕,孫壽浩,邱 峻

(1.河北科技大學(xué)信息科學(xué)與工程學(xué)院,石家莊050018;2.石家莊市經(jīng)緯度科技有限公司,石家莊050000)

基于MDBZP的C/A碼信號(hào)高靈敏捕獲算法*

王曉君1,**,尚 燕1,孫壽浩1,邱 峻2

(1.河北科技大學(xué)信息科學(xué)與工程學(xué)院,石家莊050018;2.石家莊市經(jīng)緯度科技有限公司,石家莊050000)

衛(wèi)星信號(hào)的捕獲是GPS接收機(jī)信號(hào)處理的首要任務(wù)。針對(duì)在微弱信號(hào)的情況下傳統(tǒng)捕獲方法不能很好地捕獲到衛(wèi)星信號(hào)的問(wèn)題,提出一種改進(jìn)的二倍分組塊補(bǔ)零(MDBZP)方法為高靈敏度GPS接收機(jī)提供理論依據(jù)。該方法將連續(xù)相干積分分解成普通循環(huán)相關(guān)和隨后的傅里葉變換來(lái)產(chǎn)生不同的多普勒頻移搜索,并針對(duì)連續(xù)積分跨越多個(gè)導(dǎo)航電文數(shù)據(jù)位的情況給出解決方案,同時(shí)考慮了多普勒效應(yīng)對(duì)碼長(zhǎng)造成的影響。仿真結(jié)果表明,此方法對(duì)低至15dB-Hz的微弱信號(hào)C/A碼信號(hào)仍有較好的捕獲能力。最后給出了不同載噪比情況下所對(duì)應(yīng)的檢測(cè)信噪比。

GPS接收機(jī);微弱信號(hào)處理;粗捕獲碼;高靈敏捕獲;修改的二倍分組塊補(bǔ)零

1 引 言

全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System,GNSS)已經(jīng)滲透到各個(gè)領(lǐng)域,具有廣泛的應(yīng)用前景。目前,要想完成對(duì)目標(biāo)的定位,接收到的可見(jiàn)衛(wèi)星必須多達(dá)4顆[1]。然而處于特殊情況下,例如森林和室內(nèi)等環(huán)境,可見(jiàn)衛(wèi)星數(shù)量比較少,接收機(jī)一般接收到不超過(guò)3顆衛(wèi)星信號(hào)。究其原因,處在空曠的環(huán)境下,GPS信號(hào)較強(qiáng),載噪比(C/N0)可以達(dá)到44 dB-Hz以上;處在惡劣條件下,GPS信號(hào)衰減比較嚴(yán)重,C/N0往往低于36 dB-Hz,必須提高接收機(jī)的靈敏度。

弱GPS信號(hào)的處理和應(yīng)用正在成為一個(gè)新興的研究熱點(diǎn),弱信號(hào)處理關(guān)鍵技術(shù)之一是信號(hào)捕獲[2]。D.M,Lin.討論了采用GPS軟件接收機(jī)捕獲微弱信號(hào)的問(wèn)題,采用了20 ms的相干積分時(shí)間,接著又采用非相干積分,該算法利用199個(gè)非相干積分能夠捕獲到載噪比為21 dB-Hz的信號(hào)。D.M. Akos.對(duì)捕獲弱信號(hào)的相干積分性能分析結(jié)果說(shuō)明,捕獲一個(gè)32 dB-Hz的載噪比信號(hào)至少需要8 ms的相干積分時(shí)間,捕獲一個(gè)22 dB-Hz的載噪比信號(hào)至少需要200 ms的相干積分時(shí)間,捕獲一個(gè)17 dB-Hz的載噪比信號(hào)至少需要400 ms的相干積分時(shí)間。

面對(duì)微弱GPS信號(hào)捕獲的困難,本文提出了采用較少的步驟來(lái)處理長(zhǎng)的數(shù)據(jù)相關(guān)運(yùn)算,即二倍分組塊補(bǔ)零算法(Double Block Zero Padding,DBZP),同時(shí)對(duì)該算法進(jìn)行了修改,針對(duì)長(zhǎng)時(shí)間的相關(guān)積分會(huì)跨越數(shù)據(jù)位和碼長(zhǎng)的多普勒效應(yīng)問(wèn)題給出了具體解決方案,從而為后面的信號(hào)跟蹤部分提供了準(zhǔn)確的多普勒頻移和C/A碼相位。

2 傳統(tǒng)的基于FFT的GPS信號(hào)的捕獲基本原理

傳統(tǒng)的FFT碼相位捕獲算法是輸入信號(hào)經(jīng)過(guò)傅里葉變換后與經(jīng)過(guò)傅里葉變換(FFT)的本地偽碼相乘,輸出結(jié)果經(jīng)過(guò)傅里葉逆變換(IFFT)轉(zhuǎn)換為時(shí)域信號(hào),傅里葉逆變換輸出的模值表示輸入信號(hào)與本地偽碼的相關(guān)結(jié)果,從而完成捕獲GPS信號(hào)的算法[2]。

GPS接收機(jī)RF前端信號(hào)經(jīng)過(guò)下變頻的數(shù)字中頻信號(hào)模型

其中,r(tk)是tk采樣時(shí)刻中頻信號(hào)的輸出,采樣間隔Ts=1/fs,d為導(dǎo)航數(shù)據(jù);c為接收到的PRN碼,其碼長(zhǎng)為1 023 chip,碼片速率為1.023 MHz;fIF為中頻載波頻率,fd為多普勒頻移,τ為碼相位,θ0為初始相位,n為零均值、單位方差的高斯白噪聲;A為信號(hào)幅度,它從噪聲方差σ2n歸一化為1中推導(dǎo)出, C/N0為載噪比(單位dB-Hz),有

假設(shè)^τ、^fd分別為碼相位和多普勒的頻移估計(jì)值,輸入的信號(hào)和本地的載波相乘累加后得到的同相分量I、正交分量Q分別為

圖1 FFT捕獲結(jié)構(gòu)框圖Fig.1 FFT acquisition block diagram

3 改進(jìn)的DBZP微弱C/A碼信號(hào)的捕獲算法

上面介紹的捕獲算法通常在一個(gè)C/A碼的周期里面進(jìn)行運(yùn)算就可以捕獲到信號(hào),但在信號(hào)微弱環(huán)境中卻很難實(shí)現(xiàn),想要檢測(cè)出微弱信號(hào)就要進(jìn)行長(zhǎng)時(shí)間的相關(guān)積分來(lái)提高信噪比,這就需要采用微弱信號(hào)的捕獲算法。

3.1 算法原理及實(shí)現(xiàn)步驟

修改的二倍分組塊補(bǔ)零(Modify Double Block Zero Padding,MDBZP)方法是將連續(xù)相干積分分解成普通循環(huán)相關(guān)和隨后的傅里葉變換來(lái)產(chǎn)生不同多普勒頻移的搜索,同時(shí)考慮了碼長(zhǎng)多普勒效應(yīng)所引起的問(wèn)題。圖2是基于MDBZP的捕獲算法框圖。接收信號(hào)變換成基帶信號(hào)并產(chǎn)生

通過(guò)采樣時(shí)間Ts=1/fs來(lái)設(shè)置碼時(shí)延之間的間隔,在一個(gè)C/A碼周期內(nèi)獲得Nt個(gè)碼時(shí)延采樣點(diǎn)。相干積分時(shí)間為TI,用頻率間隔fres=1/TI來(lái)產(chǎn)生多普勒倉(cāng),每個(gè)TI時(shí)間內(nèi)的數(shù)據(jù)劃分成分組的塊數(shù)等于多普勒倉(cāng)數(shù)Nfd=2fdcovTI,其中fdcov為最大的多普勒覆蓋范圍。每個(gè)分組塊包含S=fSTI/Nfd個(gè)采樣點(diǎn)。

圖2 基于MDBZP的捕獲算法框圖Fig.2 Block diagram of the MDBZP-based acquisition algorithm

第一步:將所接收的一個(gè)相干積分時(shí)間TI內(nèi)的基帶信號(hào)采樣劃分成若干個(gè)分組子塊,每?jī)蓚€(gè)相鄰的分組子塊組合成一個(gè)分組塊,并產(chǎn)生Nfd個(gè)重疊分組塊,每個(gè)重疊分組塊的數(shù)據(jù)長(zhǎng)度為2S,最后一個(gè)分組子塊與下一個(gè)相干積分時(shí)間的第一個(gè)數(shù)據(jù)子塊合成。

第二步:考慮了多普勒補(bǔ)償?shù)谋镜谻/A碼也用相同的fS采樣產(chǎn)生,并且劃分成Nfd個(gè)分組子塊,每個(gè)分組子塊用S個(gè)零進(jìn)行填充,每個(gè)分組塊的大小為2S。為了更好的對(duì)C/A碼進(jìn)行多普勒補(bǔ)償,在整個(gè)多普勒頻移范圍內(nèi)等間隔設(shè)置補(bǔ)償頻率來(lái)生成C/A碼。第i種補(bǔ)償頻率fdi為

其中,W為補(bǔ)償?shù)亩嗥绽疹l率種類。

第三步:將一個(gè)相干周期內(nèi)所有的基帶信號(hào)數(shù)據(jù)與某一種特定移位情況C/A碼完成部分相關(guān)運(yùn)算,稱之為一個(gè)捕獲步。在每個(gè)捕獲步上,在接收基帶信號(hào)的采樣和本地產(chǎn)生C/A碼采樣的每?jī)蓚€(gè)對(duì)應(yīng)的分組塊之間進(jìn)行循環(huán)相關(guān)計(jì)算,保留每個(gè)合并的分組塊的最先S個(gè)循環(huán)相關(guān)值。為了敘述清楚,用二維數(shù)據(jù)索引標(biāo)示每個(gè)分組塊的相關(guān)結(jié)果。以第一捕獲步為例,分別用(1,1),(2,1),…,(Nfd,1)標(biāo)示分組塊的相關(guān)結(jié)果。每個(gè)捕獲步從所有的分組塊中保留的點(diǎn)將組成大小為Nfd×S的數(shù)據(jù)矩陣。

第四步:每個(gè)捕獲步將產(chǎn)生S個(gè)不確定的碼時(shí)延的部分相干積分,對(duì)于一個(gè)C/A碼所有的碼時(shí)延,計(jì)算相關(guān)運(yùn)算的總捕獲步數(shù)為Nep,其中Nep為一個(gè)C/A碼周期內(nèi)計(jì)算相關(guān)運(yùn)算的捕獲步數(shù),其值為Nt/S。每個(gè)捕獲步中,參與相關(guān)運(yùn)算的C/A碼分組子塊分別向左循環(huán)移動(dòng)一次,每一次移動(dòng),第一個(gè)分組塊就會(huì)變成最后一個(gè)分組塊。在Nep個(gè)捕獲步結(jié)束時(shí),將產(chǎn)生一個(gè)大小為Nfd×Nt的數(shù)據(jù)矩陣,這個(gè)矩陣是Nep個(gè)Nfd×S的數(shù)據(jù)矩陣的行擴(kuò)展。

第五步:對(duì)矩陣的每列做Nfd點(diǎn)的FFT運(yùn)算得到相干積分,每個(gè)單元對(duì)應(yīng)于多普勒倉(cāng)和碼時(shí)延的相關(guān)運(yùn)算,得到矩陣Ci,i為多普勒補(bǔ)償頻率的序號(hào)。在矩陣Ci中對(duì)應(yīng)于第i個(gè)多普勒范圍的Ni行的相干積分才添加到總相干積分C中,Ci矩陣的其余行將被舍棄。

第六步:在進(jìn)行L次總非相干積分后,如果矩陣中的某個(gè)單元有一個(gè)值超過(guò)門限γ,那么就可以推斷出已捕獲到一顆衛(wèi)星,包含最大能量值的單元將對(duì)應(yīng)于所估計(jì)的碼時(shí)延和多普勒頻移;若沒(méi)有單元幅值超過(guò)設(shè)定的門限值,則未捕獲該衛(wèi)星信號(hào)。

3.2 多普勒效應(yīng)的補(bǔ)償

實(shí)現(xiàn)對(duì)弱信號(hào)進(jìn)行捕獲,需要進(jìn)行長(zhǎng)時(shí)間相干積分和非相干積分提高信噪比,多普勒效應(yīng)將會(huì)對(duì)接收信號(hào)的碼片長(zhǎng)度產(chǎn)生影響,會(huì)導(dǎo)致與本地的復(fù)制碼碼延遲的不同,這會(huì)造成積分增益損失[4]。為了避免該情況,將整個(gè)多普勒范圍劃分成M個(gè)小范圍來(lái)處理由碼長(zhǎng)多普勒效應(yīng)引起的問(wèn)題,定義fdi為第i個(gè)范圍的中間頻率,如公式(6)所示,以及定義Ni為在i個(gè)范圍內(nèi)的若干個(gè)多普勒倉(cāng)數(shù),并產(chǎn)生W個(gè)復(fù)制碼類型,第i個(gè)多普勒范圍的相干積分計(jì)算將采用第i個(gè)碼類型,補(bǔ)償?shù)膹?fù)制碼可表示為

3.3 處理未知位邊沿位置

信號(hào)捕獲的過(guò)程中,相對(duì)于相干積分,非相干積分的取模運(yùn)算使結(jié)果存在平方損耗,增加相干積分時(shí)間可以使平方損失減小[5]。但加長(zhǎng)相干積分運(yùn)算會(huì)出現(xiàn)跨越數(shù)據(jù)位的問(wèn)題,為解決該情況,避免使相關(guān)運(yùn)算結(jié)果跨越導(dǎo)航電文,就要從數(shù)據(jù)開(kāi)始跳變的邊沿開(kāi)始對(duì)數(shù)據(jù)進(jìn)行處理。由于一個(gè)數(shù)據(jù)位是20 ms,將存在20個(gè)不確定位邊沿的位置,每個(gè)位邊沿將對(duì)應(yīng)1 ms碼周期的起始點(diǎn),如果相干積分的起始點(diǎn)沒(méi)有與數(shù)據(jù)位的起始點(diǎn)對(duì)準(zhǔn),會(huì)在積分過(guò)程中存在一些損失。將采用Nb個(gè)并行的非相干積分?jǐn)?shù)組來(lái)保存數(shù)據(jù),Nb的選擇不要超過(guò)20,并在20個(gè)不確定邊沿范圍內(nèi),均勻地選擇其間隔。因此,可以算出Nb的選擇總數(shù)不應(yīng)該小于4。雖然4≤Nb≤20將引起一些積分損失,但是這可以減少總的處理時(shí)間,積分損失也不會(huì)影響最后的捕獲結(jié)果。

3.4 處理未知數(shù)據(jù)位

當(dāng)相干積分跨越導(dǎo)航數(shù)據(jù)位時(shí),數(shù)據(jù)位的跳變將會(huì)影響捕獲能力以及捕獲的準(zhǔn)確性[6]。算法利用一個(gè)長(zhǎng)的相干積分時(shí)間TI,其相干積分時(shí)間為一個(gè)數(shù)據(jù)位長(zhǎng)度Tdms的N倍,需要對(duì)N個(gè)數(shù)據(jù)位的最可靠組合進(jìn)行估計(jì),因此將有2N-1個(gè)矩陣對(duì)應(yīng)于某個(gè)不確定的數(shù)據(jù)位組合,從中選擇最可靠數(shù)據(jù)位的組合加到先前的總非相干積分中。由于運(yùn)算量將隨著N的增加呈級(jí)數(shù)增長(zhǎng),所以N的選擇不要過(guò)大。

3.5 信號(hào)捕獲性能分析

在信號(hào)捕獲過(guò)程中虛警概率和檢測(cè)概率是兩個(gè)重要的指標(biāo),通常依據(jù)虛警概率要求確定捕獲的門限,然后再根據(jù)檢測(cè)概率來(lái)分析信號(hào)的捕獲性能[7]。當(dāng)信號(hào)捕獲門限為γ時(shí),對(duì)應(yīng)虛警概率為

檢測(cè)概率為

其中,β=(C/N0)MT,T=0.001 s,IL-1(x)為一階修正貝塞爾函數(shù),M為預(yù)檢測(cè)積分時(shí)間TI的碼周期數(shù),L為整個(gè)捕獲過(guò)程中的非相干積分總數(shù)。

4 仿真分析

通過(guò)對(duì)真實(shí)衛(wèi)星發(fā)射的信號(hào)進(jìn)行仿真,然后再把噪聲疊加到模擬的數(shù)據(jù)上來(lái)產(chǎn)生接收機(jī)收到的信號(hào),根據(jù)不同的信噪比來(lái)產(chǎn)生不同的信號(hào)幅值[8]。首先對(duì)接收機(jī)射頻前端接收的信號(hào)進(jìn)行下變頻得到中頻載波數(shù)字信號(hào),然后把C/A碼和導(dǎo)航電文信息加到下變頻的中頻信號(hào)上,在設(shè)定噪聲為零均值、單位方差的高斯白噪聲的情況下,根據(jù)載噪比計(jì)算公式推出信號(hào)的幅度A。本文在算法仿真時(shí),數(shù)字中頻信號(hào)的載波頻率fIF=1 405 kHz,采樣頻率fS= 5 700 kHz,一個(gè)C/A碼周期的采樣點(diǎn)數(shù)為Nt= 5 700,相干積分時(shí)間TI=80 ms。將多普勒頻移搜索范圍為10 kHz,即fdcov=5 kHz。由此可知,分組塊數(shù),即多普勒倉(cāng)數(shù)Nfd=800個(gè),多普勒頻率分辨力fres=6.25 Hz,每個(gè)分組塊包含點(diǎn)數(shù)S=570個(gè)。多普勒補(bǔ)償頻率種類數(shù)W取為5,對(duì)應(yīng)于補(bǔ)償頻率為(-4 kHz-2 kHz 0 Hz 2 kHz 4 kHz),非相干積分次數(shù)L=30。在積分矩陣模塊中最大值元素出現(xiàn)位置所對(duì)應(yīng)的列坐標(biāo)就是基帶信號(hào)的碼相位,對(duì)應(yīng)的行坐標(biāo)標(biāo)定了多普勒頻率的大小。

4.1 信號(hào)的捕獲結(jié)果

采用相干積分時(shí)間TI=40 ms,相當(dāng)于2個(gè)數(shù)據(jù)位,非連續(xù)積分時(shí)間L=30,共計(jì)60個(gè)數(shù)據(jù)位時(shí)長(zhǎng)的捕獲算法。圖3顯示了在不同載噪比的情況下,歸一化信號(hào)幅值時(shí),信號(hào)的捕獲結(jié)果。隨著噪聲的增加,捕獲結(jié)果就會(huì)產(chǎn)生影響,在載噪比為19 dB-Hz的情況下,噪聲峰值就比較大,虛警概率比較大,信號(hào)不能被捕獲。

圖3 2個(gè)數(shù)據(jù)位長(zhǎng)信號(hào)的捕獲結(jié)果仿真Fig.3 Two data bits signal capture simulation

采用相干積分時(shí)間TI=80 ms,相當(dāng)于4個(gè)數(shù)據(jù)位,非連續(xù)積分時(shí)間L=30,共計(jì)120個(gè)數(shù)據(jù)位時(shí)長(zhǎng)的捕獲算法。圖4顯示了在不同載噪比的情況下,歸一化信號(hào)幅值時(shí),信號(hào)的捕獲結(jié)果。在增加連續(xù)積分時(shí)間的情況下,最大相關(guān)峰值明顯大于噪聲峰值,信號(hào)可以被捕獲。仿真結(jié)果顯示,該算法可以捕獲載噪比低至15 dB-Hz的微弱信號(hào)。

圖4 4個(gè)數(shù)據(jù)位長(zhǎng)信號(hào)的捕獲結(jié)果仿真Fig.4 Four data bits signal capture simulation

4.2 信號(hào)的檢測(cè)信噪比

在MDBZP捕獲算法中,捕獲門限的設(shè)定是非常重要的,通常是根據(jù)虛警概率和檢測(cè)概率來(lái)得到的[9]。通過(guò)運(yùn)算在矩陣模塊中有最大值元素出現(xiàn),那么根據(jù)其位置就可以確定碼相位和多普勒頻移的值,將所有單元中的最大值當(dāng)作信號(hào),在最大值周圍的值也看作是處于一個(gè)碼片內(nèi),其余所有的相關(guān)值看作噪聲。接收機(jī)前端的采樣頻率是5 700 kHz, 1 ms內(nèi)大約是5個(gè)采樣點(diǎn),所以在最大值左右各4個(gè)采樣點(diǎn)是一個(gè)碼片內(nèi)的相關(guān)值,剩余單元的所有值都當(dāng)成噪聲來(lái)看待。信號(hào)功率和噪聲功率是信號(hào)檢測(cè)過(guò)程中的兩個(gè)重要參數(shù)Pds和Pdn分別為

式中,M表示相關(guān)值中的最大值,b表示噪聲的平均值, ni表示噪聲的相關(guān)值,K表示噪聲相關(guān)值的個(gè)數(shù)。

檢測(cè)信噪比定義[10]為

從檢測(cè)信噪比中就可以分析出不同C/N0下信號(hào)的捕獲效果,如圖5所示。

圖5 信號(hào)的檢測(cè)信噪比Fig.5 Detection signal-to-noise ratio of the signal

上述的仿真結(jié)果表明,采用新的處理捕獲算法可以更好地捕獲處于低信噪比GPS信號(hào),很大程度上提高了檢測(cè)信噪比。很明顯,新的捕獲方法在捕獲低信噪比衛(wèi)星信號(hào)時(shí)比傳統(tǒng)FFT捕獲方法有了很大的改善。

5 結(jié) 語(yǔ)

本文將P(Y)碼捕獲技術(shù)應(yīng)用于微弱GPS信號(hào)的C/A碼捕獲,給出了一種改進(jìn)的基于二倍分組塊補(bǔ)零原理的微弱GPS C/A碼信號(hào)的捕獲算法及其實(shí)現(xiàn)步驟。算法充分考慮了C/A的多普勒頻率補(bǔ)償。文中討論了如何處理未知數(shù)據(jù)位和未知位邊沿位置等問(wèn)題,使信號(hào)捕獲不受數(shù)據(jù)位的影響。本文給出相干積分時(shí)間TI=80 ms、非相干積累次數(shù)L= 30時(shí)的仿真情況,結(jié)果表明,本算法對(duì)多普勒頻率的分辨力達(dá)到6.25 Hz,可以較好地捕獲到載噪比低至15 dB-Hz的微弱信號(hào),還可以確定數(shù)據(jù)位跳變邊沿的位置,可有效實(shí)現(xiàn)微弱GPS C/A碼的捕獲問(wèn)題,從而實(shí)現(xiàn)導(dǎo)航系統(tǒng)的高靈敏定位。

[1] Xie Gang.Principles of GPS and Receiver Design[M]. Beijing:Publishing House of Electronics Industry,2009.

[2] Ziedan N I.GNSS Receivers for Weak Signals[M].New York:Artech House,2006.

[3] Morton Y T,Miller M,Tsui J,et al.GPS civil signal self -interference mitigation during weak signal acquisition [J].IEEE Transactions on Signal Processing,2007,55 (12):5859-5863.

[4] Shanmugam S K,Nielsen J,Lachapelle G.Enhanced differential detection scheme for weak GPS signal acquisition[C]//Proceedings of 2007 ION GNSS.Fort Worth:ION,2007:1-14.

[5] 胡從瑋,李曉玲,安雷.一種改進(jìn)的GPS微弱信號(hào)捕獲方法[J].武漢大學(xué)學(xué)報(bào),2008,33(8):821-823. HU Cong-wei,LI Xiao-ling,AN Lei.An improved methodfor GPS weak signal acquisition[J].Journal of Wuhan University,2008,33(8):821-823.(in Chinese)

[6] Wei Yu,Bo Zheng,Watson R,et al.Differential combing for acquiring weak GPS signals[J].Signal Processing,2007(87):824-840.

[7] 翟傳潤(rùn),張婧.GPS軟件接收機(jī)微弱信號(hào)捕獲算法研究[J].電子測(cè)量技術(shù),2008,31(5):137-140. ZHAI Chuan-run,ZHANG Jing.GPS receiver weak signal acquisition algorithm[J].Electronic Measurement Technology,2008,31(5):137-140.(in Chinese)

[8] Tsui J B Y.Fundamentals of Global Positioning System Receivers a Software Approach[M].Hoboken,New Jersey:John Wiley&Sons,Inc,2005:230-235.

[9] Mark P.GNSS Solutions:Carrier-to-Noise algorithms [J].Inside GNSS,2010,5(1):20-27.

[10] Kokkonen M,Pietial S.A New Bit Synchronization Method for a GPS Receiver[C]//Proceedings of 2002 IEEE Positioning,Location,and Navigation Symposium.Palm Springs,CA:IEEE,2002:85-90.

WANG Xiao-jun was born in Chicheng, Hebei Province,in 1973.He received the Ph.D. degree from Beijing Institute of Technology in 2007.He is now a professor of Hebei University of Science and Technology and postdoctor of the 54th Research Institute of CETC.His research concerns base band signal processing in satellite navigation.

Email:eicont@139.com

尚 燕(1982—),女,河北邢臺(tái)人,碩士,河北科技大學(xué)信息學(xué)院講師,主要研究方向?yàn)樾盘?hào)處理、圖像的稀疏表示;

SHANG Yan was born in Xingtai,Hebei Province,in 1982.She is now a lecturer with the M.S.degree.Her research concerns digital image processing.

Email:shangyan_sy@126.com

孫壽浩(1988—),男,山東德州人,河北科技大學(xué)信息學(xué)院碩士研究生,目前主要研究方向高靈敏度GNSS接收機(jī)信號(hào)處理;

SUN Shou-hao was born in Dezhou,Shandong Province,in 1988.He is now a graduate student.His research concerns signal processing in high-sensitivity GNSS receiver.

Emain:sunshouhao2010@126.com

邱 峻(1972—),男,石家莊市經(jīng)緯度科技有限公司總工程師,主要從事“北斗”一代、二代系列應(yīng)用終端產(chǎn)品研發(fā)。

QIU Jun was born in 1972.He is now a chief engineer of Shijiazhang LBH Technology Co.,Ltd..He is engaged in product development of Beidou-1 and Beidou-2 application terminals.

Email:qiujun0088@163.com

High-Sensitivity Capture Algorithm of C/A Code Signal Based on Modify Double Block Zero Padding

WANG Xiao-jun1,SHANG Yan1,SUN Shou-hao1,QIU Jun2
(1.College of Information and Science,Hebei University of Science and Technology,Shijiazhuang 050018,China; 2.Shijiazhuang LBH Technology Co.,Ltd.,Shijiazhuang 050000,China)

Signal capture is the first task of signal processing in GPS receiver.Sometimes,GPS signals are too weak to be

or processed by traditional FFT methods.A method of Modify Double Block Zero Padding (MDBZP)is presented to solve the problem of weak signal acquisition for high-sensitivity GPS receivers.At the same time,the solution with respect to continuous integration time crossing multiple data bits of navigation message is also provided and the influence of Doppler effect on code length is taken into account.The simulation results show that the method is better to capture weak GPS C/A-code signal whose carrier noise ratio is even low to 15 dB-Hz.Finally,the detection SNRs are given corresponding to different CNRs.

GPS receiver;weak signal processing;C/A code;high-sensitivity capture;MDBZP

Foundation of Hebei University of Science and Technology(No.XL201064)

TN967.1

A

1001-893X(2013)10-1299-06

王曉君(1973—),男,河北赤城人,2007年于北京理工大學(xué)獲微電子學(xué)與固體電子學(xué)專業(yè)博士學(xué)位,現(xiàn)為河北科技大學(xué)教授、中國(guó)電子科技集團(tuán)公司第五十四所博士后,主要研究方向?yàn)樾l(wèi)星導(dǎo)航基帶信號(hào)處理技術(shù);

10.3969/j.issn.1001-893x.2013.10.009

2013-04-19;

2013-06-24 Received date:2013-04-19;Revised date:2013-06-24

河北科技大學(xué)校立基金項(xiàng)目(XL201064)

**通訊作者:eicont@139.com Corresponding author:eicont@139.com