摘 要： 提供一種干擾數據的拋物線波形擬合方法，利用該方法在對接收的數據進行拋物線波形擬合的過程中，無需進行復雜的逆矩陣計算，因而可以提高信號處理的效率，滿足探測器對數據處理的實時性要求。

　　關鍵詞： 計算機期刊投稿,隨機干擾,擬合算法,拋物線波形擬合,逆矩陣

　　Fast parabola fitting algorithm of stochastic interference data

　　AN Kai

　　(Shandong Aerospace Electro?technology Institute， Yantai 264670， China)

　　Abstract： A fast parabola fitting algorithm of the stochastic interference data is put forward in this paper. In the process of the parabola waveform fitting for the received data， it is unnecessary for the method to perform the inverse matrix calculation. Therefore， the efficiency of the signal processing can be improved. The method can meet the detector′s requirement of real?time data processing because it is fast enough.

　　Keywords： random interference; fitting algorithm; parabola waveform fitting; inverse matrix

　　0 引 言

　　在光電、電磁、激光和紅外等探測過程中，接收信號往往帶有隨機干擾[1?3]，只有剔除這些隨機干擾，才能提取出正確的接收數據，從而達到探測的目的。通常在無干擾情況下這些數據總是具有某些特性，例如核爆信號經若干過程處理，到達數據采集單元前端信號的波形近似于拋物線[4]。依據這些特性，利用曲線擬合方法就可以剔除隨機干擾，獲得更加精確的測量數據[5]。目前盡管曲線擬合的方法很多，如最小二乘法、偽逆矩陣法等，但計算過程比較復雜，以最小二乘法和偽逆矩陣法為例，二者都涉及復雜的逆矩陣運算，需要耗費大量的運算時間和資源，不僅嚴重影響了探測的實時性和精確性，也增加了探測器的研制成本。

　　本文將以核爆信號處理為背景，提供一種無需進行逆矩陣計算的拋物線波形擬合方法，可提高信號處理效率，從而滿足探測器對數據處理的實時性要求。

　　1 初始時刻為0的情形

　　假定在第[i]個采樣時刻接收到的帶干擾數據為[y(i)]，[i=1，2，…，N]，而拋物線的方程為：

　　[y=at2+bt+c]

　　因此有：

　　[y(1)y(2)?y(N)=1111222???1NN2?cba]

　　記：

　　[P=1111222???1NN2]

　　并稱之為擬合矩陣。于是有：

　　[PTy(1)y(2)?y(N)=PTPcba]

　　[cba=(PTP)-1PTy(1)y(2)?y(N)]

　　式中：[PT]表示矩陣[P]的轉置;[(PTP)-1]表示矩陣[(PTP)]的逆矩陣。因此：

　　記：

　　[A=k=1Ny(k)B=k=1Nky(k)C=k=1Nk2y(k)]

　　[D=3(3N2+3N+2)N(N-1)(N-2)E=12(2N+1)(8N+11)N(N2-1)(N2-4)F=180N(N2-1)(N2-4)G=-18(2N+1)N(N-1)(N-2)H=-180N(N-1)(N2-4)I=30N(N-1)(N-2)]

　　則：

　　[cba=DGIGEHIHFABC=DGIGEHIHF=AD+BG+CIAG+BE+CHAI+BH+CF]

　　即：

　　[a=AI+BH+CFb=AG+BE+CHc=AD+BG+CI]

　　2 初始時刻為[t0≠0]的情形

　　這一情形擬合算法的實現步驟包括采集數據、平移數據、中間結果計算和擬合系數計算，如圖1所示。

　　假定在第[i]個采樣時刻接收到的帶干擾數據為[y(i)，][i=t0+1，t0+2，…，t0+N，]待擬合的拋物線為[y=at2+bt+c。] 此時上述擬合矩陣中元素的數值較大，因此容易引起許多不良后果，在計算中容易導致溢出，運算更加復雜。為此，可以將采樣數據沿橫坐標左移[t0]，得到的數據為：

　　[y(i)=y(t0+i)，i=1，2，…，N]

　　利用上述逆合方法可以求出相應的拋物線[y=at2+bt+c]中的系數[a，b，c。]由坐標平移可得：

　　[at2+bt+c=a(t-t0)2+b(t-t0)+c]

　　即：

　　[at2+bt+c=at2+(b-2at0)t+at02-bt0+c]

　　因此：

　　[a=a，b=(b-2at0)c=at02-bt0+c]

　　擬合步驟如圖1所示。

　　圖1 拋物線擬合算法實現步驟

　　3 采樣周期[Δt≠1]的情形

　　以上討論針對的是采樣周期[Δt=1]的情形，其實上述結果也可以推廣到[Δt≠1]的情形，此時拋物線擬合方程組可表示為：

　　[y(1)y(2)?y(N)=1Δt(Δt)212Δt22(Δt)2???1NΔtN2(Δt)2cba]

　　或：

　　[y(1)y(2)?y(N)=1111222???1NN2cbΔta(Δt)2]

　　因此，如果按照采樣周期[Δt=1]的擬合方法，得到的擬合拋物線為：

　　[y=at2+bt+c]

　　則擬合曲線可表示為：

　　[y=a(Δt)2t2+b?Δt?t+c]

　　4 仿 真

　　對于波形為[y=-t2+t]的拋物線信號，考慮時段[t∈[0，1]]內的波形擬合問題。假定隨機干擾信號為白噪聲[N(0， 0.05)，]以采樣周期為0.01 s為例，可以得到100個采樣數據。拋物線及其擬合曲線如圖2所示。

　　采用上述擬合算法，得到的擬合拋物線為：

　　[y=-1.009 4t2+1.055 8t-0.009 4]

　　需要指出的是當白噪聲為[N(0，0.01)]時，圖2中的兩條曲線完全重合，從圖形上已無法分辨，仿真時為了從圖形上區分兩條曲線，采用方差較大的白噪聲[N(0， 0.05)。]

　　5 結 論

　　為了提高信號探測的精度，數據的采樣周期必須足夠短，對實時擬合的情形，這意味著擬合的速度必須足夠高。而在擬合過程中最復雜的運算就是逆矩陣的計算。本文中提出的方法避免了逆矩陣的運算，因此運算速度極大地提高，縮短了數據的采樣周期，從而進一步提高了信號探測的精度。

　　圖2 拋物線及其擬合曲線

　　參考文獻

　　[1] 劉香翠，程翔，張良，等.納米TiN煙幕干擾激光和紅外性能研究[J].激光與紅外，2011，41(8)：920?924.

　　[2] 劉佐成.高頻開關電源中的電磁干擾問題及電磁兼容[J].現代電子技術，2007，30(10)：166?167.

　　[3] 劉衛東，劉尚合，魏明.抑制周期性窄帶干擾的時域重構優化算法性能分析[J].高電壓技術，2010，36(5)：1211?1216.

　　[4] 劉代志，王仁明，慕曉冬.天基核爆信號識別處理系統技術方案研究[J].核電子學與探測技術，2003，23(4)：292?295.