基于天波重構(gòu)技術(shù)的eLORAN信號周期識別算法

劉時堯 華 宇 張首剛

①(中國科學(xué)院國家授時中心 西安 710600)

②(中國科學(xué)院精密導(dǎo)航定位與定時技術(shù)重點研究室 西安 710600)

③(中國科學(xué)院時間頻率基準重點研究室 西安 710600)

④(中國科學(xué)院大學(xué) 北京 100049)

1 引言

羅蘭(LOng RAnge Navigation, LORAN)C系統(tǒng)是美國海軍于20世紀50年代出于軍事目的而發(fā)展起來的陸基導(dǎo)航系統(tǒng)，是一種遠程雙曲線導(dǎo)航定位系統(tǒng)[1]。增強型羅蘭(enhanced LORAN, eLORAN)是由Loran-C系統(tǒng)演變而來，其采用Eurofix數(shù)據(jù)鏈技術(shù)進行三態(tài)脈位調(diào)制(Pulse Position Modulation,PPM)實現(xiàn)增強數(shù)據(jù)信息發(fā)播[2]，具有作用距離遠、穩(wěn)定性好、抗干擾能力強等優(yōu)點。近年來，美國、韓國、英國、俄羅斯等多國均重新啟動了新型eLORAN系統(tǒng)的建設(shè)研究，而我國也于2017年啟動了“高精度地基授時系統(tǒng)”項目，旨在成為北斗衛(wèi)星系統(tǒng)的可靠地基備份以完善我國的星地一體化高精度定位、導(dǎo)航、授時(Positioning, Navigation and Timing, PNT)系統(tǒng)。

周期識別是信號處理中的核心問題，其穩(wěn)定性極易受到天波、連續(xù)波、交叉干擾以及各種噪聲的影響，可能引起跳周(1載波周期=10 μs)。另外，對于添加了數(shù)據(jù)調(diào)制的eLORAN信號，很容易因為將天波信號誤捕獲為地波而造成解碼誤碼率增大，最終導(dǎo)致授時錯誤及導(dǎo)航定位的極大誤差。天波干擾指的是經(jīng)過電離層反射的天波輻射信號與地面波間接收時延差小于信號的有效時長時對地面波造成的重疊效應(yīng)[3]。根據(jù)國際海事無線電技術(shù)委員會(Radio Technical Commission for Maritime, RTCM)發(fā)布的長波接收機最低性能標準，天波一般位于地波信號后的37.5～1500 μs[4]。為了保證PNT精度，eLORAN接收機要求利用地波信號進行周期識別，因此一般將地波脈沖組第1個脈沖的第3載波周期正向過0點(30 μs處)作為標準過0點。

近年來出現(xiàn)了不少周期識別算法的研究[5—9]，包括優(yōu)化包絡(luò)相關(guān)法、多徑估計延遲鎖定環(huán)(Multipath Estimation Delay Lock Loop, MEDLL)算法等，但均是單純針對尋找標準過0點的研究，而忽略了討論如何判斷天波干擾和地波識別問題。另外這些均是時域方法，容易受到各種強干擾而導(dǎo)致信號畸變，因此在惡劣條件下的適用性不高。針對以上難點，本文通過比較結(jié)合現(xiàn)有方法并提出聯(lián)合算法，在低信噪比、強天波干擾環(huán)境下解決了天地波的識別問題，并在此基礎(chǔ)上有效構(gòu)造天波信號，利用去除天波后的偽地波成功實現(xiàn)周期識別，為后續(xù)解調(diào)解碼等過程提供了保障。

2 信號體制

Loran-C信號以脈沖組的形式發(fā)射，而eLORAN信號與傳統(tǒng)的Loran-C信號的區(qū)別主要是在第3～8個脈沖增加了數(shù)據(jù)調(diào)制，標準eLORAN脈沖時域波形用公式表示為[10]

實際中由于多跳天波離地波較遠，因此只考慮距離較近的一跳天波的影響。受到一跳天波干擾的信號時域波形可用公式表示為

其中，sground(t),ssky(t),τg,τs,A,B分別為地波、天波信號及其時延差和幅值；n(t)為噪聲干擾。

設(shè)定天地波時延差Δτ為100 μs，天地波幅度比(Sky-Ground Ratio, SGR)為6 dB并加以10 dB的白噪聲，根據(jù)式(1)、式(2)可以得到eLORAN標準脈沖及混合信號波形如圖1所示。

根據(jù)《船用羅蘭C接收設(shè)備技術(shù)標準》[10]，羅蘭C信號的信噪比(Signal Noise Ratio, SNR)指信號接收電平與大氣噪聲均方根電平值之比，可用無量綱的數(shù)或分貝數(shù)表示。其中信號電平定義為脈沖包絡(luò)起點之后25 μs處連續(xù)波有效值電平，本文之后的仿真SNR均是根據(jù)以上定義為標準。另外，信號進入接收機后，首先會經(jīng)過梳狀濾波器(40次左右線性累加)，在保證信號不畸變的同時提高信噪比(15 dB左右)，故文中的仿真SNR如無特別說明均是在天線接收端信噪比基礎(chǔ)上增加15 dB后的結(jié)果。

3 算法描述

eLORAN接收機進行信號處理，實現(xiàn)高精度PNT功能是依賴地波信號的，而天波信號的混入會影響周期識別效果，造成以下兩類致命的錯誤：(1)將天波誤捕獲為地波信號造成的地波識別錯誤；(2)天波影響地波的位置估計而導(dǎo)致的跳周錯誤。為了解決以上問題，提高尋找地波第3周過0點的準確度，本文對天波特性進行了深入分析，并提出了一種針對天波干擾的周期識別聯(lián)合算法。該算法可以分為以下3個部分，如圖2所示。

(1) 天地波識別模塊。當天地波位置較近或天波信號較強時，加之各種噪聲干擾等因素的影響，會造成地波識別錯誤，導(dǎo)致數(shù)十微秒甚至上百微秒級別的時延估計(Time Delay Estimation, TDE)誤差并最終引起授時定位偏差。而該模塊的目的正是在于將識別信號首先定位在地波信號上。

(2) 天地波分離模塊。經(jīng)研究發(fā)現(xiàn)，高強度的天波干擾會對地波的位置估計產(chǎn)生較大影響。因此，經(jīng)天地波識別、決策后，需要進行天地波分離的信號進入此模塊，通過自適應(yīng)匹配算法得到天波的準確位置及強度并消除天波，得到較為干凈的偽地波信號。

(3) 周期識別模塊。該模塊將對以上兩個模塊的輸出結(jié)果進行進一步處理，利用基本無天波干擾的偽地波信號尋找到標準過零點位置，實現(xiàn)周期識別。

3.1 天地波識別模塊

天地波識別模塊的核心功能是判斷天波的有無，而該過程的實現(xiàn)主要是利用信號時延差、幅度比估計，以及二者的相互關(guān)系進行分析判斷。目前已經(jīng)存在一些較為成熟的天地波時延差估計方法，包括多重信號分類(MUltiple SIgnal Classification,MUSIC), (Auto Regressive Moving Average,ARMA)及ESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques)算法等[11—14]，都可以在不是很極端的信噪比環(huán)境下達到TDE需求。其中MUSIC算法分辨率更高，峰值更尖銳，但信噪比要求也最高；ARMA模型及ESPRIT算法抗噪能力稍強，但分辨率及估計準確度較差。但以上提到的算法均存在一個致命問題，即均是在“天波存在”假設(shè)下提出的技術(shù)，研究重點主要為天、地波的時延差估計，不能體現(xiàn)信號幅度特征，無法在盲測時利用幅度信息輔助判斷是否存在天波。近年來的一些新的研究成果，如EXIP-WRELAX(EXtended Invariance Principle Weighted fourier transform and RELAXation)算法、神經(jīng)網(wǎng)絡(luò)算法等[15—17]也都存在同樣的問題，因此在實際應(yīng)用時容易受到天波或交叉干擾等影響而導(dǎo)致錯捕或誤識別且無法糾正，這也是目前各種接收終端中周期識別一直無法較好實現(xiàn)的關(guān)鍵所在。

經(jīng)典頻譜相除算法的最大優(yōu)勢在于可以較好地反映信號幅度特征，這是判斷信號存在與否，實現(xiàn)天地波信號的識別功能的關(guān)鍵。因此本文對此算法進行了深入研究以改進其在強干擾環(huán)境下的分辨率及抗干擾性能要求。其基本原理推導(dǎo)為[18,19]

eLORAN信號的能量主要集中在85～115 kHz，因此在頻譜相除算法中通常會使用窗函數(shù)以降低帶外噪聲及干擾。文獻[18,19]均使用了帶寬50 kHz的漢寧窗進行頻域截斷，但是此種窗在天地波時延差較小時由于分辨率不足而無法達到效果，因此本文對窗函數(shù)進行了改進。通過對各類窗函數(shù)的研究比較可以得到一致的結(jié)論：窗口越窄，去噪效果越好，IFFT變換后峰越平坦(分辨率低)，估計越準確；反之，窗口越寬，峰越尖銳，分辨率越高，而去噪效果也越差。另外分析表明，窗函數(shù)的下降梯度越大，兩側(cè)越平坦，則干擾抑制能力越強。因此為了兼顧高分辨率和干擾抑制這兩方面需求，可以在更寬窗口的基礎(chǔ)上對其梯度進行優(yōu)化。經(jīng)過大量對比驗證，本文選擇了在式(4)中取140 kHz寬的Parzen窗及h=2.3，以提高窗口下降梯度并使帶外部分更加平緩。根據(jù)以上改進，通過式(3)、式(4)步驟提取到了更惡劣條件下(SNR=5 dB(天線接收端—10 dB)，SGR=18 dB, Δτ=38 μs)的天地波信息。

通過圖3(a)可以看出，改進的Parzen窗在保證窗寬度的同時對梯度有了明顯改善，圖3(b)的峰信息對比也顯示出改進窗的高分辨率可以成功將天波和地波區(qū)分開來，并維持了較好的去噪效果。但超過強天波閾值( Δτ≤60 μs, SGR＞5 dB)的高強天波會嚴重影響TDE的準確度，加之窗函數(shù)引起δ峰左右兩側(cè)的旁瓣效應(yīng)(Δτ＜35.5 μs，分別小于—22 dB,—33 dB)及噪聲干擾等會使TDE誤差遠大于5 μs。這對周期識別來說極其致命，同時也是以往許多研究中都無法克服的問題。另外，強噪聲及其他臺站的交叉干擾信號等造成的尖峰也可能會使真實信號峰被淹沒，導(dǎo)致信號識別錯誤，故僅利用1次頻譜相除的峰信息分辨天波和地波是缺乏可信度的。為此，本文基于大數(shù)法則的思想設(shè)計了天地波識別算法，其步驟如下：

(1) 利用捕獲框?qū)⑿盘柨刂圃诓蓸訑?shù)據(jù)中心附近，連續(xù)采集N組峰信息(建議不少于50組以保證樣本數(shù)充足，每組1 ms)，累加并取極大值點Pm，則由于信號已捕獲，Pm必為天波或地波信號標識；在其左右各取1個次大峰Pl, Pr，計算并預(yù)存峰值比及時延差信息；

(2) 在每組峰信息中尋找最接近Pm位置的峰Qm(i)，并在其左右各統(tǒng)計1個次大峰Ql(i)和Qr(i)共N組。若某側(cè)存在信號，則應(yīng)有大量次大峰集中在Pl或Pr附近，因此可分別統(tǒng)計此兩點±10 μs范圍內(nèi)Ql(i)和Qr(i)的個數(shù)，利用經(jīng)驗個數(shù)閾值(可取N的50%～75%)分別判斷左右兩側(cè)是否存在eLORAN信號，并根據(jù)旁瓣的幅度、位置特性排除將旁瓣誤判為信號的可能性；

(3) 根據(jù)天波判斷結(jié)果進行信號輸出抉擇：

(a) 若存在天波則根據(jù)信號時延關(guān)系確定天波、地波，并保存預(yù)估幅度比D、時延差 Δτ及等信息。由于弱天波干擾對地波的位置判斷的影響不大，因此可根據(jù)前段提到的強天波閾值進行判決，將參數(shù)D和 Δτ在閾值范圍內(nèi)的信號輸入天地波分離模塊。

(b) 若判斷無天波信號干擾或參數(shù)值在強天波閾值范圍之外，則直接將當前信號及預(yù)估的地波位置信息τg0輸入周期識別模塊。

該模塊是后續(xù)模塊的基礎(chǔ)，其準確率主要受交叉干擾、突發(fā)干擾及噪聲的強度及頻度的影響，其中窗函數(shù)及帶通濾波器(Band-Pass Filter, BPF)可以高效抑制大量的帶外干擾及噪聲，而梳狀濾波器可壓制絕大多數(shù)的交叉干擾信號。另外，在接收設(shè)備非常接近某干擾臺站時，交叉干擾信號強度也有可能高于遠處臺站的跟蹤信號。此種情況下可同時捕獲該干擾信號并利用其脈沖組重復(fù)周期關(guān)系對消大部分交叉干擾能量。

接收信號經(jīng)過各種前端處理后，仍可能有少量無法預(yù)測的同頻干擾或突發(fā)干擾等影響信號峰，因此在工程應(yīng)用時應(yīng)根據(jù)不同地點測試數(shù)據(jù)綜合考慮并給出經(jīng)驗個數(shù)閾值以供接收設(shè)備設(shè)定。

該算法的計算復(fù)雜度主要取決于每組數(shù)據(jù)采樣點數(shù)M，以及頻譜相除和峰值檢測等步驟。其中：傅里葉變換(Fast Fourier Transform, FFT)、IFFT的時間復(fù)雜度均為O(M ·log2M)；除法及窗函數(shù)運算復(fù)雜度均為O(M)；檢測最大峰、次大峰的復(fù)雜度為O(M)。故利用N組峰信息進行天地波識別的時間復(fù)雜度為：O(N ·(2·M ·log2M+M)+N ·M+2·M+N ·M)=O(N ·M ·log2M)。相比于天地波分離模塊所需的10 MHz采樣率，此模塊對采樣率要求不高，在實際應(yīng)用中可使用現(xiàn)有接收機通常使用的2 MHz采樣率以節(jié)省80%的計算量。以上算法利用了大數(shù)法則的理論，利用多組峰信息并選擇合適的個數(shù)閾值進行聯(lián)合判決，極好地解決了交叉干擾、同頻干擾以及其他突發(fā)干擾的影響，在利用模擬源及新型接收終端進行測試時幾乎沒有出現(xiàn)過錯誤，較好地實現(xiàn)了不同天波干擾環(huán)境下的天地波識別，準確地將捕獲框鎖定在地波信號上。

3.2 天地波分離模塊

天地波分離的核心是準確地重構(gòu)天波，而強天波干擾會極大影響地波的位置估計精度，使TDE誤差大于5 μs并導(dǎo)致周期識別跳周。針對此問題，本文設(shè)計了自適應(yīng)去天波算法，包含了兩階搜索過程以及用于匹配的混合信號幅度2維網(wǎng)格，以保證準確估計天波時延τs及幅度B的同時盡可能節(jié)省計算量。其中網(wǎng)格中映射的數(shù)值為無干擾、無噪聲環(huán)境下不同幅度(范圍0.1～18倍，步長0.1)、時延差(范圍37.5～200 μs，步長0.1 μs)天波干擾下的混合信號幅值，由仿真軟件生成后存儲于DSP存儲器中。該算法其流程圖如圖4所示，并包括如圖4的步驟。

(1) 將預(yù)估的時延差Δτ及幅度比D與預(yù)存的混合信號幅度網(wǎng)格C進行匹配，根據(jù)其數(shù)值與經(jīng)過帶通濾波器的實際信號幅度間的比例關(guān)系以及 Δτ的搜索范圍，通過式(5)預(yù)估出第1階搜索中時延循環(huán)對應(yīng)的天波幅度搜索中心B0(k)。由于不同時延差對應(yīng)的天波幅度估計初值差別可能較大，因此，在第1階段時延差循環(huán)中得到不同的B0(k)可以達到縮小幅度循環(huán)范圍，提高效率的目的

(2) 以τs0,B0(k)為中心進行第1階段搜索，選擇相對較寬的搜索范圍以及較大的步長循環(huán)構(gòu)造并去除偽天波，得到偽地波組；利用偽地波組的峰值信息與基準信號組合并與偽地波組對消得到殘差組并計算均方根(Root Mean Square, RMS)，利用最小均方原則得到粗略的天波時延差及幅度參數(shù)值,B′。

其中，信號輸入接收機后會進行電平調(diào)整，故本文后續(xù)仿真統(tǒng)一將采樣信號幅度標準化為10，并選擇10 MHz采樣率以滿足天波幅度及時延估計精度需求。另外，進入此模塊的信號中天波功率相對較大，因此利用頻譜相除技術(shù)時對天波參數(shù)的估計誤差相對較小，時延誤差范圍一般不超過5 μs，而幅度誤差一般不超過25%。據(jù)此，兩階段搜索范圍及步長推薦設(shè)置如下：

根據(jù)上文中天地波分離模塊的詳細步驟說明以及圖4的具體實現(xiàn)流程，式(6)—式(9)結(jié)合各過程變量給出了幾處關(guān)鍵節(jié)點的計算方法，其中式(6)、式(7)表示兩階搜索中偽天波對消、偽地波對消及殘差RMS的循環(huán)計算方法；式(8)、式(9)為構(gòu)造最優(yōu)偽天波、產(chǎn)生最優(yōu)偽地波的計算方法

根據(jù)圖5(a)所示，去天波前的天、地波時延估計分別為138.9 μs, 95.0 μs，與真實值偏差分別為1.1 μs, 5.0 μs，天地波幅度比D估計為2.1562。而運用本節(jié)提出的天波分離算法后，根據(jù)仿真過程量可知天波時延估計誤差為0，幅度比估計為2.0662，偏差為3.31%，即基本上去除了天波，如圖5(b)、圖5(c)所示。

3.3 周期識別模塊

以上兩個模塊成功地實現(xiàn)了天地波的識別和分離，極大降低了天波干擾的影響并準確定位地波信號。而周期識別模塊的目的是準確尋找偽地波信號的標準過零點，其中的關(guān)鍵正是識別正確的載波周期。目前對此問題的方法研究主要分為以下兩類：

(1) 時域方法：周期累加值法、波形匹配法、優(yōu)化包絡(luò)相關(guān)算法、MEDLL算法等。

時域算法對信噪比要求較高，其中周期累加值法、周期和值法、波形匹配法的信噪比下限為23～25 dB；峰值比算法有了一定的性能提升，在信噪比達到15 dB 以上時能夠有效檢測標準過0點[9]，因此有一定的工程應(yīng)用價值；包絡(luò)相關(guān)算法分辨率較低，容易受到波形畸變影響且計算量相對較大；MEDLL算法在包絡(luò)相關(guān)算法的基礎(chǔ)上結(jié)合了包絡(luò)求導(dǎo)技術(shù)[6]，一定程度提高了分辨率及信噪比性能，但當信噪比低于14 dB或天波過強時甚至無法分辨。

(2) 變換域方法：MUSIC算法、ARMA模型、ESPRIT算法、IFFT頻譜相除算法等。

相對于時域算法，各種變換域方法復(fù)雜度更高，但在信噪比條件及穩(wěn)定性方面有著明顯優(yōu)勢，3.1節(jié)總結(jié)了幾種方法的性能特點。其中，ARMA模型及ESPRIT算法抗噪能力稍強，但分辨率及估計準確度較差，容易導(dǎo)致跳周；MUSIC算法及頻譜相除算法相對穩(wěn)定，在天波干擾較小時估計精確度極高，但MUSIC算法的信噪比要求相對較高[14]且計算量更大。

經(jīng)過對上述算法理論的深入分析可以看出，各種時域算法較為簡單實用，但對信號接收環(huán)境有較高的要求。如峰值比算法，在噪聲及干擾環(huán)境適宜時有較高的穩(wěn)定性，但也容易受到天波殘差以及各種噪聲干擾引起的波形畸變的影響，導(dǎo)致跳周錯誤，在工程應(yīng)用中可結(jié)合環(huán)境估計結(jié)果決策使用。而本文研究重點為惡劣環(huán)境下的周期識別，噪聲環(huán)境較差，因此更適合變換域方法。下文選擇了其中穩(wěn)定性高，計算量相對較小且信噪比要求更低的頻譜相除算法進行周期識別過程(經(jīng)仿真章節(jié)驗證，SNR下限接近—2 dB，遠低于上述方法)。

利用頻譜相除方法推算標準過零點位置的基本流程與其他各種算法大致相同，即首先識別正確的載波周期，再檢測標準過零點位置，其具體步驟描述如下：

(1) 不需要天地波分離的信號在天地波識別后直接輸入此模塊，并將預(yù)存的地波位置τg0記為PA；對天地波分離模塊輸出的偽地波則利用頻譜相除技術(shù)，預(yù)估信號的起始位置點PA。由于偽地波中的天波殘差較小，對IFFT計算的分辨率要求不高，故可以在計算中選擇較窄的窗函數(shù)以實現(xiàn)更大程度的噪聲抑制，突出信號尖峰。

(2) 根據(jù)估計所得地波起始位置點PA推后30 μs預(yù)估第3載波周期過0點PB。雖然基本無天波干擾，但PA點的估計誤差會導(dǎo)致PB點與準確的標準過0點PC之間存在一個小的偏差。

(3) 根據(jù)PB點位置尋找最接近PB點的載波正向近0點PC并作為推算的標準過0點。

根據(jù)以上方式繼續(xù)對3.2節(jié)示例中輸出的偽地波進行處理，可得信號位置信息如圖6所示。

從幅度歸一化位置估計對比圖6(a)中可以看出，去除天波干擾后，偽地波信號的起始位置估計誤差已控制在半個載波周期范圍內(nèi)(本例為0.3 μs)，隨后如圖6(b)所示，點PC準確地標識了標準過0點的位置。該方法利用了頻譜相除技術(shù)的高穩(wěn)定性，后續(xù)仿真結(jié)果可以驗證，在去除大部分天波干擾或無天波干擾時，利用偽地波可以得到精確的標準過0點位置以及地波的起始位置。

4 仿真實驗及結(jié)果分析

第3節(jié)分解討論了該自適應(yīng)周期識別算法的各個模塊。本節(jié)對整體算法進行了充分的仿真實驗及結(jié)果展示，并對主要過程量進行了對比分析。其中，由于天地波識別算法在仿真實驗及不同地點的拷機測試中均無出錯，故沒有給出實驗數(shù)據(jù)。

根據(jù)eLORAN接收機的最低性能標準要求[4]，當天線接收端信號SNR≥—10 dB, SGR≤12～26 dB,Δτ在37.5～60 μs范圍時接收設(shè)備應(yīng)能進入鎖定狀態(tài)，而新型eLORAN接收機更要求在此范圍內(nèi)能夠成功解碼。本文算法與3.1節(jié)提到的各種類似的變換域方法均有其信噪比極限，在噪聲功率過大時噪聲峰會淹沒信號峰，算法有效性也會急劇下降。其中在利用不同的窗及BPF處理且保證信號不失真時，頻譜相除算法信噪比下限不低于—2 dB(天線接收端—17 dB左右)。而在實際應(yīng)用中，一般認為天線接收端信噪比小于0 dB時為低信噪比環(huán)境。基于以上分析，本文針對各種惡劣環(huán)境，對仿真中的可變參數(shù)作如下設(shè)計：

(1) 輸入信號：標準eLORAN信號+天波信號+隨機白噪聲，采樣率為10 MHz；

(2) 天地波時延差Δτ：38～60 μs，步長為1 μs；天地波幅度比D：SGR=0～24 dB，步長為6 dB；

(3) 信噪比SNR：0～10 dB(天線接收端—15～—5 dB)，步長為1 dB。

本文在MATLAB2017環(huán)境下，獨立調(diào)整各參數(shù)，并在不同情況下分別運行了1000次。由于篇幅所限，文中只展示了SGR=0 dB, 18 dB和24 dB幾種情況下的部分仿真結(jié)果及數(shù)據(jù)統(tǒng)計，如表1所示。

表1體現(xiàn)了在SGR=0 dB時該算法極高的性能。由于算法運行過程中幅度比估計值超過強天波門限范圍，因此跳過了天地波分離而直接進入周期識別模塊。結(jié)果表明，此方式在節(jié)省了效率的同時非常有效地實現(xiàn)了周期識別，驗證了算法的合理性。

表1 周期識別錯誤次數(shù)及準確率(SGR=0 dB)

圖7、圖8分別展示了SGR=18 dB和24 dB時的天波位置及幅度的估計誤差。可以看出，位置和幅度的估計精度隨著SNR的增加而提升。經(jīng)統(tǒng)計，位置誤差RMS最大值分別為0.0126 μs和0.0105 μs，幅值誤差率RMS最大值分別為1.09%和0.53%。通過對兩個參數(shù)非常準確的估計，可以有效地重構(gòu)并去除天波。

通過3.1節(jié)討論可知，在強天波干擾環(huán)境下，改進前的頻譜相除技術(shù)甚至無法區(qū)分天地波，即使改進了窗函數(shù)提高分辨率也仍存在極大的位置誤差，導(dǎo)致周期識別跳周頻發(fā)；而如表2、表3所示，經(jīng)過了天地波分離過程，極大地提高了在不同條件下周期識別的正確率，基本高于99%。

表2 天地波分離后周期識別錯誤次數(shù)及準確率(SGR=18 dB)

表3 天地波分離后周期識別錯誤次數(shù)及準確率(SGR=24 dB)

5 結(jié)論

本文所提聯(lián)合算法主要是針對各種惡劣環(huán)境所設(shè)計，這也是目前各種周期識別算法研究的難點所在。根據(jù)上文的算法討論和部分仿真結(jié)果展示以及更多的仿真數(shù)據(jù)統(tǒng)計，可以得出以下結(jié)論：

(1) 對窗函數(shù)的改進效果明顯，有效地將高強度天波干擾下的eLORAN天、地波信號區(qū)分開來，并可適配提出的天地波識別算法，成功實現(xiàn)對天、地波信號的準確定位。另外通過決策過程準確判斷是否需要進行天波分離，大幅節(jié)省了后續(xù)過程的計算量；

(2) 對天波的參數(shù)估計準確度指標均隨SNR及天地波時延差的增加而提高，在極端惡劣的噪聲及天波干擾環(huán)境下仍然能夠達到極高的水平，且本文所提2階搜索方式在大幅提高效率的同時可以有效地重構(gòu)并去除天波信號，輸出較為穩(wěn)定的偽地波信號；

(3) 本聯(lián)合算法可以保證在接收天線輸入端SNR=—15 dB的極低信噪比下也能夠保持極高的水準，周期識別準確率基本超過99%，完全可以滿足新型eLORAN接收機的需求。

eLORAN接收終端的PNT精度很大程度上取決于周期識別的準確度，但高強度的天波干擾一直是信號處理中的一個關(guān)鍵問題，嚴重限制了周期識別算法的完整性及有效性。針對這一問題，本文提出一種自適應(yīng)天波識別、分離和周期識別的聯(lián)合算法。經(jīng)過理論分析和各種不同接收環(huán)境下的仿真實驗，驗證了算法的可行性。本算法首先克服了經(jīng)典IFFT頻譜相除算法無法同時兼顧高分辨率和時延估計準確度需求的弊端，并結(jié)合其能夠體現(xiàn)信號幅度特征的特點，有效解決了以往研究中常被忽略的天地波信號識別問題。其次，能夠在惡劣的條件下成功將天地波的識別、估計、分離及尋找標準過零點等過程有機結(jié)合，在提高效率的同時保證了周期識別的正確率，具有良好的魯棒性和可實現(xiàn)性。同時，精確的天波估計，也為接收機處理交叉干擾，提高解調(diào)解碼的成功率提出了新的思路和理論參考。