曾 鑫，熊 波，于德海，邢永強(qiáng)

(1.海軍某工程建設(shè)指揮部，北京 100036;2.海軍航空工程學(xué)院，山東 煙臺 264001;3.92514 部隊，山東 煙臺 264001)

彈載高度表是控制導(dǎo)彈飛行的關(guān)鍵部件，其工作性能的好壞直接決定著導(dǎo)彈飛行軌跡的穩(wěn)定程度和命中概率［1］。掠海飛行的導(dǎo)彈容易受海面各種雜波干擾，影響導(dǎo)彈正常的掠海飛行，嚴(yán)重時可能導(dǎo)致導(dǎo)彈墜毀。因此高度表的狀態(tài)好壞必須在使用前進(jìn)行測試。本文提出了采用數(shù)字?jǐn)M合和濾波相結(jié)合的方法對高度表進(jìn)行高度信號峰值與脈沖個數(shù)進(jìn)行測試。該方法比直接進(jìn)行數(shù)字濾波的數(shù)字處理過程在幅頻特性上衰減變小。因此，在高度表對高度信號進(jìn)行峰值測試時，該方法具有較高的測量精度。此外，對干擾序列進(jìn)行濾波時可以選擇較低截止頻率，相對于直接濾波輸出的信號具有較好的平滑特性。進(jìn)行脈沖計數(shù)時，平滑性較高的波形計數(shù)精度較高。

1 某型彈載高度表測試系統(tǒng)設(shè)計

1.1 高度表工作原理

高度表常見的調(diào)制方式有:正弦波調(diào)制、三角波調(diào)制和鋸齒波調(diào)制。反艦導(dǎo)彈通常使用鋸齒波調(diào)制的無線電高度表。以鋸齒波調(diào)制方式為例，設(shè)發(fā)射信號為鋸齒波調(diào)制的連續(xù)波，Tm為調(diào)制周期(遠(yuǎn)大于高度表最大工作高度對應(yīng)的回波延遲t2-t1)，Δf 為調(diào)頻帶寬(遠(yuǎn)大于高度表最大工作高度對應(yīng)的差拍信號fb)，調(diào)制斜率為K =Δf/Tm，目標(biāo)的徑向速度可忽略(反艦導(dǎo)彈的速度遠(yuǎn)大于目標(biāo)的速度)。t1時刻高度表發(fā)射的線性調(diào)頻波的頻率為f1，在t2時刻收到回波信號，此時發(fā)射信號的頻率變?yōu)閒2?；夭ㄐ盘柡桶l(fā)射信號在混頻器中混頻、檢波后得到差拍信號fb。高度表工作的信號調(diào)制如圖1 所示［2］。

差拍信號得到高度信號的計算公式為

其中:H 為所需測量的高度;c 為電磁波傳播速度。

處于動態(tài)平衡時Δf 和c 是常數(shù)，fb在瞬態(tài)平衡下也是常數(shù)，此時高度H 與調(diào)制周期成正比。采用直接進(jìn)行周期測試，難以保證高度表測量的精度。鋸齒波周期是固定的，不同高度產(chǎn)生的調(diào)制周期對應(yīng)不同的鋸齒波個數(shù)。因此，可通過鋸齒波計數(shù)的方式對調(diào)制周期的長度進(jìn)行測量:

其中:K1，K2表示恒定差拍下的不同調(diào)制斜率。

1.2 高度表測試系統(tǒng)

當(dāng)導(dǎo)彈掠海飛行時回波信號的干擾來自整個海平面。海面雜波影響回波信號的頻率，從而影響調(diào)制周期的長度。對海面雜波的處理，由高度表調(diào)制體制完成。本文的重點(diǎn)是將調(diào)制信號進(jìn)行轉(zhuǎn)化，然后對調(diào)制信號的峰值和長度進(jìn)行測試。高度表測試系統(tǒng)的結(jié)構(gòu)［2］如圖2 所示。

適配器的作用是將高度表產(chǎn)生的調(diào)制信號轉(zhuǎn)化為數(shù)據(jù)采集卡可測試的數(shù)據(jù)，即將信號進(jìn)行電壓轉(zhuǎn)換。在信號通過分壓轉(zhuǎn)換電路過程中，產(chǎn)生了一定的干擾信號;測試系統(tǒng)采用DAQ(數(shù)據(jù)采集)技術(shù)進(jìn)行數(shù)據(jù)采集時，會帶來一定的干擾;此外，高度表內(nèi)部產(chǎn)生的調(diào)制信號也帶有噪聲。因此，需要對采集信號進(jìn)行去噪處理。

圖2 測試系統(tǒng)結(jié)構(gòu)

2 信號去噪方法

2.1 數(shù)字低通濾波算法

調(diào)制信號是一個帶噪聲的鋸齒波信號。將信號進(jìn)行FFT 變換，作出信號的幅頻特性如圖3 所示。

圖3 調(diào)制信號頻譜圖

由幅頻特性圖可知，調(diào)制信號中含有不同頻率的干擾信號，需要對該信號進(jìn)行低通濾波。信號頻譜圖中，不同的頻率分量通過截止頻率為ωc濾波器后產(chǎn)生的幅值衰減為

因此，將帶噪聲的鋸齒波信號通過巴特沃思低通濾波器后產(chǎn)生一定程度的衰減和失真，影響信號測試幅值的正確性。

2.2 基于擬合的低通濾波

本文采用擬合［3-4］和數(shù)字濾波的方法，對信號進(jìn)行數(shù)據(jù)處理，減小信號在濾波過程中產(chǎn)生的幅度失真，提高信號幅值測試和脈沖計數(shù)的精度。

設(shè)被測信號在一個調(diào)制周期的采樣序列為f(1)，f(2)，…，f(n)，…，f(N)，稱為觀測值。對采樣后的信號進(jìn)行擬合，擬合的函數(shù)為

其中:k 為鋸齒波的斜率，是一個未知量;s 為鋸齒波的寬度，為已知量;m 為對應(yīng)的鋸齒波個數(shù)。

擬合函數(shù)f^(n)的求解過程比較簡單，可通過最小二乘法計算函數(shù)中的未知參數(shù)k 即可。首先計算擬合函數(shù)各點(diǎn)對應(yīng)的模型值序列，然后與采樣序列對應(yīng)的觀測值求均方差

圖4 基于擬合的數(shù)字低通濾波器信號處理流程

再對差值序列求方差最小時對應(yīng)參數(shù)值k0。將k0代入擬合函數(shù)得到不含參數(shù)的模型序列f^(1)=k0，f^(2)=2k0，…，然后對采樣點(diǎn)序列與擬合后的函數(shù)求殘差，獲得干擾序列Δf(1)，Δf(2)，…。對該序列進(jìn)行巴特沃思低通濾波，獲得濾波后的序列Δf '(1)，Δf '(2)，…，最后將濾波后的干擾序列與模型序列疊加，獲得新的觀測序列f '(1)，f '(2)，…。

根據(jù)上述信號處理原理，作出被測信號的數(shù)據(jù)處理流程如圖4 所示。

3 計算機(jī)仿真

原信號是帶有噪聲的鋸齒波信號，且信噪比較高，如圖5所示。

圖5 原始信號

采用巴特沃思低通濾波器進(jìn)行濾波后的信號產(chǎn)生了一定程度的幅值衰減，如圖6 所示。

圖6 直接數(shù)字濾波后的信號

經(jīng)過擬合處理的殘差序列是一個穩(wěn)定的帶有噪聲的平穩(wěn)信號，如圖7 所示。

圖7 數(shù)字?jǐn)M合后的殘差信號

對殘差序列進(jìn)行巴特沃思低通濾波，并在擬合信號的基礎(chǔ)上還原得到高度表如圖8 所示。

圖8 消除噪聲后的信號

從圖8 可以看出信號幅值衰減較小，而且信號失真程度較小。

4 實驗結(jié)果對比

通過Labwindows/CVI 軟件［5］控制高度表高度信號的測試和信號處理流程，將直接經(jīng)過數(shù)字濾波與基于擬合的數(shù)字濾波后的測試結(jié)果進(jìn)行比較。以100 m 高度測試為例，高度表測試系統(tǒng)通過數(shù)據(jù)采集卡采集5 組的數(shù)據(jù)，每組對應(yīng)的理論脈沖個數(shù)為1000。通過上述2 種方法分別對數(shù)據(jù)進(jìn)行處理，結(jié)果如表1 所示。

表1 測試結(jié)果對比

從表1 可以看出，基于擬合的數(shù)字濾波技術(shù)具有較高的測試精度。

5 結(jié)束語

本文采用數(shù)字?jǐn)M合與濾波相結(jié)合的方法對高度表測試信號進(jìn)行去噪處理。該方法是在數(shù)字濾波和函數(shù)擬合的基礎(chǔ)上發(fā)展起來的一種針對鋸齒波信號的濾波算法，較好地完成了高度表測試過程中高度測量的任務(wù)。該方法還可以對其他類型的確定信號進(jìn)行數(shù)據(jù)處理，提高信號的測試精度。對于現(xiàn)代化的武器裝備，各種測試指標(biāo)信號類型基本上是已知的，可以通過該方法提高系統(tǒng)的測試精度，提高對武器系統(tǒng)進(jìn)行射前檢查判斷的準(zhǔn)確程度。因此，該方法具有較為廣泛的應(yīng)用前景。

［1］王少云.KGW 無線電高度表測試儀的研制［J］.儀器儀表用戶，2006，13(5):25.

［2］常秀豐，李國林，廖輝榮.基于虛擬儀器的某型高度表測試系統(tǒng)軟件研制［J］.現(xiàn)在電子技術(shù)，2010，34(10):175-178.

［3］閆蓓，王斌，李媛.基于最小二乘法的橢圓擬合改進(jìn)算法［J］.北京航空航天大學(xué)學(xué)報，2008，34(3):295-298.

［4］趙世田，趙東標(biāo)，付瑩瑩.測量數(shù)據(jù)點(diǎn)的高精度B 樣條曲線擬合算法［J］. 計算機(jī)集成制造系統(tǒng)，2010，16(8):1708-1713.

［5］王建新，楊世鳳，隋美麗.LabWindows/CVI 測試技術(shù)及工程應(yīng)用［M］.北京:化學(xué)工業(yè)出版社，2006.