陳雪峰 ，劉 洋（.西安航空電子科技有限公司 陜西 西安 70075；.西安愛生技術(shù)集團(tuán) 陜西 西安 70075）

衛(wèi)星導(dǎo)航抗干擾天線，一般由陣列天線、方向圖形成網(wǎng)絡(luò)和自適應(yīng)處理單元3部分組成。陣列天線接收到多路衛(wèi)星信號(hào)后，經(jīng)過通道電路同時(shí)輸入到方向圖形成網(wǎng)絡(luò)和自適應(yīng)處理單元，方向圖形成網(wǎng)絡(luò)對(duì)各路衛(wèi)星信號(hào)加權(quán)并疊加，作為抗干擾天線的輸出，自適應(yīng)處理單元對(duì)陣列天線送來的信號(hào)進(jìn)行抽樣處理并進(jìn)行自適應(yīng)權(quán)計(jì)算后，給方向圖形成網(wǎng)絡(luò)提供相應(yīng)的加權(quán)量，再反過來對(duì)各路衛(wèi)星信號(hào)進(jìn)行加權(quán)調(diào)整，相當(dāng)于使各陣元的增益或相位發(fā)生改變，從而在天線陣的方向圖中產(chǎn)生對(duì)著干擾源方向的零點(diǎn)，將干擾信號(hào)抵消，從而達(dá)到抗干擾的目的[1-2]。

在實(shí)際使用中發(fā)現(xiàn)抗干擾天線的抗干擾性能時(shí)好時(shí)壞，抖動(dòng)較大，很不穩(wěn)定，魯棒性不強(qiáng)。主要原因是一般都直接把上面自適應(yīng)處理單元計(jì)算的權(quán)值更新給方向圖形成網(wǎng)絡(luò)，而工程中自適應(yīng)處理單元輸出結(jié)果不是很平滑，有時(shí)還會(huì)出現(xiàn)跳變，所以就會(huì)影響抗干擾天線的性能，為了克服此問題，本論文提出自適應(yīng)處理單元在計(jì)算出權(quán)值后不直接更新，而是按多項(xiàng)式平滑濾波算法對(duì)濾波器權(quán)值進(jìn)行修正，然后更新自適應(yīng)處理單元，保證抗干擾天線性能的穩(wěn)定性；

1 平滑濾波的數(shù)學(xué)模型

用矩陣形式表示為：X=Pβ+ε

根據(jù)最小二乘法求βj的最小二乘估計(jì)，即通過求

設(shè)X(n)中的一組數(shù)據(jù)為X(i),i=-M,…,0, …,M,構(gòu)造計(jì)算X(i)在N+а處的L階（L≥0的整數(shù)）微商的最優(yōu)線性無偏估計(jì)值[5]，得:

權(quán)系數(shù)WN-1的無偏最優(yōu)估計(jì)，與P，L采樣點(diǎn)的個(gè)數(shù)N以及а有關(guān)，當(dāng)P，L，N和а給定的情況下，WN-1是完全確定的量，可看作一個(gè)“窗”，常常取固定數(shù)量的一段觀測值設(shè)計(jì)“窗”，然后采用移動(dòng)“窗”的方法來對(duì)整段數(shù)據(jù)進(jìn)行處理.直觀上來看，即是對(duì)選取的N個(gè)觀測數(shù)據(jù)進(jìn)行了P階多項(xiàng)式擬合，獲得N個(gè)權(quán)系數(shù)WN-1，把它們存貯起來，下一次對(duì)一個(gè)“新”的觀測數(shù)據(jù)和相應(yīng)的觀測序列相乘相加，即可獲得下一時(shí)刻的估值，逐步后移，直到全部擬合成為止，這種算法簡單易行，同時(shí)也提高了計(jì)算速度。據(jù)此，這種算法稱作“移窗法”，“移動(dòng)”的擬合值為：

式中s為移動(dòng)步長，為整數(shù)，當(dāng)s=1時(shí)，為逐點(diǎn)移動(dòng)法。

2 抗干擾權(quán)系數(shù)的平滑應(yīng)用

2.1 平滑濾波原始數(shù)據(jù)的錄取

這里主要完成對(duì)自適應(yīng)抗干擾權(quán)值系數(shù)原始數(shù)據(jù)的積累，在抗干擾天線中, 自適應(yīng)處理單元按線性約束最小均方誤差(LCMV)準(zhǔn)則做自適應(yīng)權(quán)系數(shù)的計(jì)算,求解空域的自適應(yīng)權(quán)值[6]：

RST為接收到的數(shù)據(jù)協(xié)方差陣，根據(jù)各陣列天線接收到的數(shù)據(jù)估計(jì)協(xié)方差陣RST；

式中，M時(shí)間的快拍數(shù)，X(t) 為t時(shí)刻的數(shù)據(jù)向量，上標(biāo)H表示共軛操作；

可見在原始權(quán)系統(tǒng)的錄取中，一次自適應(yīng)抗干擾權(quán)系數(shù)計(jì)算包括：協(xié)方差矩陣計(jì)算，估耗時(shí)約0.64 us,方差矩陣的求逆計(jì)算，耗時(shí)約2.27 ms,自適應(yīng)權(quán)Wa計(jì)算，耗時(shí)約1.45 ms,方向圖形成網(wǎng)絡(luò)計(jì)算，耗時(shí)約1.2 ms,一次自適應(yīng)抗干擾權(quán)系數(shù)計(jì)算時(shí)間為5.56 Ms，考慮到實(shí)際情況，本論文中自適應(yīng)權(quán)值計(jì)算間隔為10 ms。

2.2 抗干擾權(quán)平滑濾波函數(shù)的確定

從本論文的第二節(jié)可知，在P，L，N和а確定后，抗干擾權(quán)平滑濾波函數(shù)就確定了。因?yàn)樽赃m應(yīng)權(quán)值數(shù)據(jù)在此只需要對(duì)其做位置平滑濾波, 所以取L=0 (L=1 為速度的平滑濾波;L= 2為加速度的平滑濾波)。在實(shí)時(shí)處理中,為了對(duì)下一個(gè)時(shí)刻進(jìn)行預(yù)測,與實(shí)測值比較,取a＞0,而在事后處理中,常采用中心平滑技術(shù)a=-(N- 1) /2。為了保證自適應(yīng)權(quán)值的更新時(shí)間不超過0.1 s,本論文中自適應(yīng)權(quán)值采樣點(diǎn)的個(gè)數(shù)N=10。由于數(shù)據(jù)的多樣性,綜合考慮偶然誤差和方法誤差, 在工程中擬合多項(xiàng)式一般選4～6階，在考慮到本應(yīng)用的實(shí)時(shí)性,擬合多項(xiàng)式定為4階,即p=4。

當(dāng)確定了窗的大小N,p,L及a后, 濾波系數(shù)就確定了,但濾波系數(shù)會(huì)帶來對(duì)實(shí)測數(shù)據(jù)的放大或縮小,因此需對(duì)權(quán)系數(shù)做歸一化處理。而后采用移窗技術(shù),對(duì)整段遙測數(shù)據(jù)進(jìn)行平滑濾波處理。

2.3 抗干擾權(quán)系數(shù)的平滑

論文前面已經(jīng)對(duì)北斗天線前端接收的數(shù)據(jù)按自適應(yīng)濾波準(zhǔn)則計(jì)算出了抗干擾權(quán)系數(shù)，并確定了抗干擾權(quán)平滑濾波函數(shù)，這里將對(duì)這些原始的抗干擾權(quán)系數(shù)按照平滑函數(shù)進(jìn)行平滑，抗干擾權(quán)系數(shù)平滑濾波可按處理可按以下算法流程進(jìn)行，如圖1所示。

圖1 抗干擾權(quán)系數(shù)的平滑流程Fig.1 The smooth flow of anti-jamming coefficient

3 應(yīng)用效果

圖2 測試場景Fig.2 The test scenario

圖3 靜態(tài)沒有平滑的抗干擾權(quán)系數(shù)Fig.3 Static not smooth anti-jamming coefficient

圖5 動(dòng)態(tài)沒有平滑的抗干擾權(quán)系數(shù)Fig.5 dynamic not smooth anti-jamming coefficient

我們對(duì)上述平滑算法在抗干擾天線樣機(jī)上進(jìn)行了效果測試，測試的場景圖如圖2所示，測試在暗室內(nèi)進(jìn)行，衛(wèi)星信號(hào)通由衛(wèi)星信號(hào)源通過天線向抗干擾天線輻射，同時(shí)通過3個(gè)輻射線向抗干擾天線輻射3個(gè)干擾，此時(shí)抗干擾天線接收信號(hào)并進(jìn)行抗干擾處理后，輸出到衛(wèi)星接收機(jī)進(jìn)行衛(wèi)星信號(hào)的解算。為了分析平滑算法在北斗抗干擾天線中的應(yīng)用效果，我們在同一個(gè)場景下，分別記錄了算法應(yīng)用前后的數(shù)據(jù)，測試分靜態(tài)和動(dòng)態(tài)，具體情況如下。

1）靜態(tài)測試

在衛(wèi)星導(dǎo)航抗干擾天線和干擾源都靜止不動(dòng)的情況下進(jìn)行抗干擾測試,在測試中把自適應(yīng)抗干擾權(quán)系數(shù)從抗干擾天線中導(dǎo)出，使用平滑濾波前后數(shù)據(jù)分別如圖3和圖4所示，很明顯平滑效果很好，并且對(duì)干擾的突變，濾波算法也能很好還原變化中的權(quán)值；另外從衛(wèi)星接收機(jī)上觀察導(dǎo)航解算在加入平滑濾波后變的十分穩(wěn)定。

2）動(dòng)態(tài)測試

在衛(wèi)星導(dǎo)航抗干擾天線在轉(zhuǎn)臺(tái)上以每秒5度的速度轉(zhuǎn)動(dòng)的情況下進(jìn)行動(dòng)態(tài)抗干擾性能測試,在測試中把自適應(yīng)抗干擾權(quán)系數(shù)從抗干擾天線中導(dǎo)出，使用平滑濾波前后數(shù)據(jù)分別如圖5和6所示，平滑效果也很好，整個(gè)濾波算法也能很好還原權(quán)值變化，沒有失真；另外從衛(wèi)星接收機(jī)上觀察導(dǎo)航解算在加入平滑濾波后比沒加前抗干擾性能穩(wěn)定了很多。

圖4 靜態(tài)加平滑后的抗干擾權(quán)系數(shù)Fig.4 Static smooth anti-jamming coefficient

圖6 動(dòng)態(tài)加平滑后的抗干擾權(quán)系數(shù)Fig.6 dynamic smooth anti-jamming coefficient

4 結(jié) 論

試驗(yàn)證明本論文所確定的抗干擾權(quán)平滑濾波函數(shù)及濾波處理過程是種有效, 對(duì)抗干擾權(quán)值有良好的平滑效果, 本論文成功利用正交多項(xiàng)式平滑濾波算法解決了，北斗抗干擾天線在實(shí)際使用中抗干擾性能時(shí)好時(shí)壞，魯棒性不強(qiáng)的問題。

[1]費(fèi)華連. 導(dǎo)航戰(zhàn)中的GPS干擾與抗干擾技術(shù)[J].航空電子技術(shù)2001, 32(1): 18-22.

FEI Hua-lian .GPS interference and anti-interference technology in Navigation War[J].Avionics Technology,2001,32(1): 18-22.

[2]仉智權(quán)，朱暢.功率倒置法在北斗抗干擾天線中的應(yīng)用[ J].電子設(shè)計(jì)工程, 2012, 20(4):42-46.

ZHANG Zhi-quan,ZHU Chang. Application of power inversion algorithm in Compass anti-jamming antenna system[J].ElectronicDesign Engineering, 2012, 20(4):42-46.

[3]胡峰,孫國基. 基于多項(xiàng)式擬合的容錯(cuò)平滑與容錯(cuò)微分平滑[ J].工程數(shù)學(xué)學(xué)報(bào), 2000(17):28-31.

HU Feng,SUNGuo-ji. Fault-tolerant smoothing and tolerant differential smoothing based on polynomial fitting[J].Journal of Engineering Thematics, 2000 (17):28-31.

[4]Thornley DJ.Anisotropic multidimensional Savitzhy Golay kernals for smoothing,differentiation and reconstruction[J].Departmental Technical Reports, 2006(7):45-49

[5]WIKIPEDIA.SavitzkyGolay smoothing filter[EB/OL](2010-02-15). http://en.wikipedia.org/wiki/savitzky-Golay_smoothing_filter.

[6]張琳,初海彬.衛(wèi)星導(dǎo)航用戶機(jī)空時(shí)濾波抗干擾技術(shù)[J].北京航空航天大學(xué)學(xué)報(bào), 2007, 33(1):76-80.

ZHANG Lin CHU Hai-Bin.Space-time interference suppress iontechnologyOfsatellite navigation receiver[J]. Journal of Beijing University of Aeronautics and Astronautics, 2007,33(1): 76-80.