譚博, 鄭華, 裴承鳴

(西北工業(yè)大學(xué) 動力與能源學(xué)院, 陜西 西安 710072)

一種用于FTPVS數(shù)據(jù)處理的非平穩(wěn)度定義

譚博, 鄭華, 裴承鳴

(西北工業(yè)大學(xué) 動力與能源學(xué)院, 陜西 西安 710072)

為解決非平穩(wěn)程度量化問題及優(yōu)化非平穩(wěn)信號處理方法的精度,針對試驗(yàn)數(shù)據(jù)特點(diǎn),提出了一種基于時(shí)頻分析的非平穩(wěn)度定義方法,使得處理方法的各項(xiàng)參數(shù)能夠隨非平穩(wěn)度變化進(jìn)行實(shí)時(shí)調(diào)整,從而提高數(shù)據(jù)處理結(jié)果的精度。通過仿真試驗(yàn),對非平穩(wěn)度定義的合理性和實(shí)用性進(jìn)行了研究,并通過風(fēng)洞試驗(yàn)及飛行試驗(yàn)的數(shù)據(jù)處理證實(shí)了其工程實(shí)用性及對非平穩(wěn)數(shù)據(jù)處理方法精度的影響。試驗(yàn)結(jié)果表明,所提出的非平穩(wěn)度不僅可以可信地表征信號非平穩(wěn)程度,而且還可以輔助提高FTPVS數(shù)據(jù)處理結(jié)果的精度。

非平穩(wěn)度; 替代數(shù)據(jù); FTPVS

0 引言

在傳統(tǒng)的臺階式顫振試驗(yàn)中,往往要求試驗(yàn)的各項(xiàng)參數(shù)在達(dá)到指定水平后保持一段時(shí)間,以完成試驗(yàn)信號的采集工作,以此保證采集信號的分段平穩(wěn)性,使得傳統(tǒng)的平穩(wěn)信號分析方法可以用于后續(xù)處理工作,但也因此產(chǎn)生了諸如周期長、成本高、不符合真實(shí)飛行器的飛行狀態(tài)且在臨界速度附近的安全性較差等不足。為了克服這些缺點(diǎn),近年來出現(xiàn)了一類連續(xù)變速的顫振試驗(yàn)(FTPVS)方法。在這類試驗(yàn)中,各項(xiàng)參數(shù)不斷變化,數(shù)據(jù)采集工作同步進(jìn)行,通過對信號的實(shí)時(shí)處理以研究試驗(yàn)對象的當(dāng)前狀態(tài)。對這類試驗(yàn)而言,其采集信號多為非平穩(wěn)隨機(jī)過程,這類信號處理結(jié)果的性能會受到信號非平穩(wěn)程度的影響。當(dāng)非平穩(wěn)程度低時(shí),甚至可以采用平穩(wěn)信號處理方法進(jìn)行分析;而當(dāng)非平穩(wěn)程度高時(shí),若不能及時(shí)對處理方法做出相應(yīng)調(diào)整,將會很難取得理想的結(jié)果, 為分析處理帶來難度。這使得針對非平穩(wěn)程度的量化工作變得十分必要[1-3]。

至于非平穩(wěn)度的定義及計(jì)算方法,由于非平穩(wěn)信號在各個(gè)領(lǐng)域內(nèi)有不同的特點(diǎn),目前尚未有明確統(tǒng)一的標(biāo)準(zhǔn),現(xiàn)有的幾種定義多為針對特定應(yīng)用領(lǐng)域內(nèi)的實(shí)際問題而提出,往往只適用于一類或幾類信號,不能用于解決其他領(lǐng)域內(nèi)的問題,因此無法用于FTPVS信號非平穩(wěn)度的量化。而對于這類試驗(yàn)信號而言,其顯著特點(diǎn)為能量分布隨時(shí)間會發(fā)生變化,在試驗(yàn)臨界狀態(tài)下這一變化會明顯加劇。由此認(rèn)為，通過對能量分布變化劇烈程度的描述,可以得到這類信號的非平穩(wěn)度。針對這一特點(diǎn),本文提出了一種基于時(shí)頻分布的非平穩(wěn)度定義及計(jì)算方法,隨后通過仿真數(shù)據(jù)及試驗(yàn)數(shù)據(jù)對定義的合理性及工程實(shí)用性進(jìn)行了研究,最后結(jié)合工程試驗(yàn),研究了其對改善非平穩(wěn)信號處理方法精度的影響。結(jié)果表明,本文定義的非平穩(wěn)度對能量分布變化具有良好的線性響應(yīng),能夠可信地表示出信號能量分布變化的劇烈情況,依據(jù)此非平穩(wěn)度調(diào)整處理方法的參數(shù)可以取得更加理想的分析結(jié)果。

1 FTPVS信號特點(diǎn)

對于試驗(yàn)中的采集信號而言,可以通過時(shí)頻表示的方法將頻譜隨時(shí)間變化的一維時(shí)域信號x(t)或頻域信號X(ω)映射成為時(shí)間-頻率平面上的二維信號Sx(t,ω),以Hermite函數(shù)作為窗函數(shù)的時(shí)頻分布為例,得到信號x(t)的Hermite窗函數(shù)時(shí)頻分布如下:

式中:hk(t)為k階Hermite函數(shù)。

對于傳統(tǒng)顫振試驗(yàn)而言,通常認(rèn)為其采集信號為分段平穩(wěn)信號。在每一段采集時(shí)間內(nèi)，信號能量分布不會隨時(shí)間發(fā)生變化,因此其在任意時(shí)間點(diǎn)的局部能量分布應(yīng)等于其全局能量分布的均值,即:

式中:〈·〉表示分布的期望。

對于FTPVS試驗(yàn)而言,由于在試驗(yàn)過程中各項(xiàng)因素不斷變化,其能量分布Sx,K(t,f)在采集時(shí)間內(nèi)也不停地變化,這一變化在接近臨界狀態(tài)時(shí)的幅度會大幅加劇。此時(shí)其局部能量分布與全局能量分布之間不再存在上式中的關(guān)系,為了刻畫局部能量分布變化的劇烈程度,提出了一種基于時(shí)頻分布的非平穩(wěn)度定義及計(jì)算方式。

2 非平穩(wěn)度定義

借助Theiler等[4]提出的替代數(shù)據(jù)法,可以得到原始信號的具有平穩(wěn)性質(zhì)的替代數(shù)據(jù),這些替代數(shù)據(jù)具有與原始信號相同的幅頻特性,且保持了其二階統(tǒng)計(jì)矩不變。依據(jù)第1節(jié)的討論可知,替代數(shù)據(jù)局部時(shí)頻分布應(yīng)與全局時(shí)頻分布的均值相等。而對于原始信號,其局部時(shí)頻分布應(yīng)為一個(gè)彼此不同的分布集合,{S(t1,f),S(t2,f), …,S(tn,f)},計(jì)算這個(gè)集合的方差,就可以得到一個(gè)表征原始信號x(t)頻域能量分布隨時(shí)間變化的“起伏”程度的量θ1[5]:

θ1=E(cn-E(c)) (n=1,…,N)

其中:

cn=D(S(tn,f),〈S(tn,f)〉n)

式中:D(…)為兩個(gè)分布之間距離的計(jì)算方法。為了靈敏地反映信號的幅度以及頻率變化,這里選用Basseville[6]給出的距離定義。

由平穩(wěn)信號的特點(diǎn)可知,在理想的條件下替代數(shù)據(jù)的θ值應(yīng)為零,但實(shí)際應(yīng)用中,由于采用相位的隨機(jī)性,θ值往往會分散在一個(gè)相對較小的范圍內(nèi),將該范圍記為平穩(wěn)過程的θ值區(qū)間。對多個(gè)替代數(shù)據(jù)的θj值進(jìn)行分布統(tǒng)計(jì)并設(shè)置閾值,給出非平穩(wěn)度的定義為:

式中:θ1為原始信號的θ值;θ0為替代數(shù)據(jù)θ值得到的平穩(wěn)區(qū)間的閾值。易知信號的DNS值恒為正且無上限,平穩(wěn)信號的DNS值小于或等于1,而非平穩(wěn)信號的DNS值則大于1。

3 仿真試驗(yàn)驗(yàn)證

為了對定義DNS的合理性進(jìn)行驗(yàn)證,參考文獻(xiàn)[7-8],通過為一個(gè)頻率為2 Hz、采樣率為128、采樣長度為128的正弦信號引入幅值系數(shù)和常系數(shù)來考察二者對信號DNS的影響,最終計(jì)算得到的DNS結(jié)果如表1所示。二者的計(jì)算公式如下:

x(t)′=kC+kAx(t)

上述的三類傳統(tǒng)算法主要基于圖像整體信息的變化來選定關(guān)鍵幀，容易造成關(guān)鍵幀選取錯(cuò)誤、計(jì)算量大、實(shí)時(shí)性差等問題。因此，本文在此基礎(chǔ)上使用一種基于深度學(xué)習(xí)的目標(biāo)檢測方法，通過建立卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Network，CNN)模型，分類提取視頻中列車頭部、尾部及車身所在關(guān)鍵幀，使得基于深度學(xué)習(xí)的目標(biāo)檢測在關(guān)鍵幀提取的應(yīng)用中成為可能。

式中:kC為常系數(shù);kA為幅值系數(shù)。

表1 不同幅值系數(shù)及常系數(shù)下DNS的值Table 1 DNS of signals with different parameters

由計(jì)算結(jié)果可以看到,由于本文定義的DNS表征的是信號能量分布隨時(shí)間變化的快慢程度,而幅值系數(shù)不會改變信號能量分布變化,所以它對信號的DNS值沒有影響;但是常數(shù)系數(shù)的加入相當(dāng)于為信號添加了低頻的趨勢信號,會影響信號的能量分布,因此計(jì)算得到的DNS值有小幅度波動,但是總體而言,仍處于前文討論的平穩(wěn)信號DNS范圍內(nèi)。

為了考察所定義的DNS對信號非平穩(wěn)程度表示的可信度,利用仿真生成變幅和變頻兩類正弦信號以觀察DNS值隨頻率或幅值變化速率的改變情況。由前文可知,為信號引入整體幅值系數(shù)不影響信號DNS值,因此對幅值變化速率的考量以幅值變化倍率的速度為自變量,而對頻率變化則包括起始頻率和終止頻率兩個(gè)方面。所生成仿真信號采樣率為128,采樣長度為1 024。最終得到對應(yīng)DNS值的結(jié)果如圖1和圖2所示。

圖1 變幅信號DNS值計(jì)算結(jié)果Fig.1 DNS of variable amplitude signals

圖2 變頻信號DNS值計(jì)算結(jié)果Fig.2 DNS of variable frequency signals

圖1中,x軸為信號幅度變化的倍率ABL,曲線的最低點(diǎn)對應(yīng)平穩(wěn)信號的DNS值,位于最低點(diǎn)左側(cè)的各點(diǎn)對應(yīng)信號幅值隨時(shí)間減小的信號對應(yīng)的DNS值,DNS最大值對應(yīng)的信號其幅值在終止時(shí)刻為初始時(shí)刻的0.1倍;位于最低點(diǎn)右側(cè)的各點(diǎn)表示信號幅值隨時(shí)間增大時(shí)的DNS值,最大值對應(yīng)的信號幅值在終止時(shí)刻為初始時(shí)刻的10倍,曲線兩側(cè)均近似線性且較為對稱。這一現(xiàn)象表明,信號幅值產(chǎn)生等倍率的縮小或增大時(shí)對應(yīng)的DNS值一致,且與DNS值近似呈現(xiàn)正比關(guān)系。在圖2中,x軸為掃頻信號初始時(shí)刻的頻率f0,y軸為掃頻信號終止時(shí)刻的頻率ff,z軸為信號DNS值。由圖中可以看到,在起始頻率和終止頻率相等的xy平面對角線上的DNS值相對最小,這是因?yàn)檫@些DNS值為不同頻率正弦信號的計(jì)算結(jié)果。DNS值在平面兩側(cè)逐漸上升,且表現(xiàn)出較好的線性與對稱特性,這表明DNS值僅與頻率變化速率相關(guān),具有較好的指示意義。

4 顫振試驗(yàn)驗(yàn)證

FTPVS是一類新興的飛行器顫振試驗(yàn)方法,與傳統(tǒng)的試驗(yàn)方法相比,其根本區(qū)別在于其采集數(shù)據(jù)非平穩(wěn),這使得傳統(tǒng)的顫振數(shù)據(jù)處理方法難以應(yīng)用。由于非平穩(wěn)數(shù)據(jù)處理方法的數(shù)值性能很大程度上受到觀測信號非平穩(wěn)程度的影響,因此,對被測信號的非平穩(wěn)度進(jìn)行測量變得十分必要。下面以實(shí)際的風(fēng)洞試驗(yàn)及飛行試驗(yàn)數(shù)據(jù)為例,驗(yàn)證本文方法在實(shí)際試驗(yàn)數(shù)據(jù)處理中的應(yīng)用效果。

4.1 顫振風(fēng)洞試驗(yàn)

某次風(fēng)洞試驗(yàn)的觀測通道試驗(yàn)數(shù)據(jù)及對應(yīng)的風(fēng)速隨時(shí)間的變化趨勢如圖3和圖4所示,被測信號DNS值變化如圖5所示。

圖3 風(fēng)洞試驗(yàn)實(shí)測信號Fig.3 Measured signal in wind tunnel test

圖5 信號的DNS值Fig.5 DNS of signal

4.2 顫振飛行試驗(yàn)

本次試驗(yàn)數(shù)據(jù)來自某型飛機(jī)跨聲速顫振飛行試驗(yàn)[9]。所選用數(shù)據(jù)中包含了跨聲速階段的采集數(shù)據(jù),采集數(shù)據(jù)如圖6所示。對飛行高度和飛行速度進(jìn)行歸一化處理后,得到的隨時(shí)間變化情況如圖7所示。計(jì)算被測信號每一秒的DNS值,得到其隨時(shí)間的變化如圖8所示。

圖6 飛行試驗(yàn)實(shí)測信號Fig.6 Measured signals of flight test

圖7 飛行參數(shù)Fig.7 Flight parameters

圖8 信號的DNS值Fig.8 DNS of signal

與風(fēng)洞試驗(yàn)相比,飛行實(shí)測信號受到的影響更多,因此信號的DNS變化也更加多樣。由圖7可以看到,在19～30 s時(shí)間段內(nèi),飛行器的高度表讀數(shù)及速度表讀數(shù)同時(shí)發(fā)生了跳變,這段時(shí)間正對應(yīng)飛行器突破音障的時(shí)間,飛行器受到外部條件的影響,其固有模態(tài)發(fā)生了變化,致使信號能量分布在短時(shí)間內(nèi)大幅變動。由圖8可以看到,信號的DNS值在19～30 s時(shí)間段發(fā)生了突變,而突破音障后這一影響便不再存在,信號的DNS值再次回落至較低水平,這說明本文DNS定義對結(jié)構(gòu)模態(tài)的改變十分敏感,也證實(shí)了本文方法的有效性。

4.3 非平穩(wěn)度在FTPVS數(shù)據(jù)處理中的應(yīng)用

在4.1節(jié)和4.2節(jié)的試驗(yàn)數(shù)據(jù)分析中,可以看出本文定義的非平穩(wěn)度可以客觀地反映出監(jiān)測信號的非平穩(wěn)程度變化,從而為判斷測試對象的當(dāng)前狀態(tài)提供可靠的依據(jù)。同時(shí),應(yīng)用信號非平穩(wěn)度還可以提高數(shù)據(jù)處理方法的精度。下面參考文獻(xiàn)[10]，以非平穩(wěn)信號處理中,常見的應(yīng)用粒子濾波建立系統(tǒng)的時(shí)變自適應(yīng)模型方法為例,說明如何應(yīng)用非平穩(wěn)度提高建立模型的精度。

選用4.1節(jié)中的風(fēng)洞試驗(yàn)信號作為分析對象。首先采用粒子數(shù)設(shè)為100的一般粒子濾波方法建立被測對象的四階系統(tǒng)模型，所得模型輸出與實(shí)際采集信號的相對誤差δ如圖9所示。

圖9 一般粒子濾波結(jié)果相對誤差Fig.9 Relative error of traditional particle filtering results

由圖9可見,在信號非平穩(wěn)度低的時(shí)段(0～200 s),采用一般粒子濾波方法經(jīng)過短時(shí)間的自適應(yīng)調(diào)整,就可建立相對誤差較小的系統(tǒng)模型;而在非平穩(wěn)度上升的時(shí)段中(200～290 s),隨著信號的非平穩(wěn)度上升,相對誤差逐漸增加,直至試驗(yàn)結(jié)束時(shí)仍存在較大的相對誤差。此時(shí)應(yīng)用傳統(tǒng)的方法已經(jīng)難以建立對象的精確系統(tǒng)模型。接下來利用信號的非平穩(wěn)度對粒子濾波方法中的粒子數(shù)目進(jìn)行調(diào)整,其數(shù)目變化與非平穩(wěn)度變化曲線一致,圖10為調(diào)整后得到的模型輸出與實(shí)際采集信號的相對誤差。

圖10 調(diào)整參數(shù)后粒子濾波結(jié)果Fig.10 Particle filtering results after parameter adjustment

由圖10可見,當(dāng)信號非平穩(wěn)度上升時(shí),由于粒子數(shù)目也隨之調(diào)整,因此在相對誤差達(dá)到一定程度后開始逐漸減小,至試驗(yàn)結(jié)束時(shí),這一誤差已經(jīng)減小至一個(gè)可以接受的范圍內(nèi)。由此可見,依據(jù)非平穩(wěn)度調(diào)節(jié)算法，可以在很大程度上提高相關(guān)參數(shù)分析結(jié)果的精度,具有良好的工程應(yīng)用前景。

5 結(jié)束語

本文借助替代數(shù)據(jù)生成方法,給出了一種適用于FTPVS數(shù)據(jù)處理的非平穩(wěn)度的定義及計(jì)算方法。仿真結(jié)果證實(shí)了其對于信號的非平穩(wěn)性具有良好的分辨能力,且對于非平穩(wěn)度的變化具有較強(qiáng)的指示能力。通過對實(shí)際連續(xù)變速風(fēng)洞試驗(yàn)以及飛行試驗(yàn)數(shù)據(jù)的非平穩(wěn)度進(jìn)行計(jì)算,驗(yàn)證了這一定義在實(shí)際數(shù)據(jù)處理的應(yīng)用中,可以明確地反映出被測信號能量的變化,從而反映出試驗(yàn)對象結(jié)構(gòu)模態(tài)的變化,具有較強(qiáng)的理論研究價(jià)值和潛在的工程應(yīng)用前景。

[1] 王宏禹,邱天爽,陳喆.非平穩(wěn)隨機(jī)信號分析與處理[M].北京:國防工業(yè)出版社,2008:1-52.

[2] 向馗,蔣靜坪.時(shí)間序列的非平穩(wěn)度分析[J].科技通報(bào),2007,23(1):1-5.

[3] 張海勇.基于局域波法的非平穩(wěn)隨機(jī)信號分析中若干問題的研究[D].大連: 大連理工大學(xué),2001.

[4] Theiler J,Eubank S,Longtin A.Testing for nonlinearity in time series:the method of surrogate data[J].Physica D,1992,58(1):77-94.

[5] Borgnat P,Flandrin P,Honeine P,et al.Testing stationarity with surrogates:a time-frequency approach[J].IEEE,Transactions on Signal Processing,2010,58(7):3459-3470.

[6] Basseville M.Distances measures for signal processing and pattern recognition[J].Signal Process,1989,18(4):349-369.

[7] 趙善紅.顫振邊界預(yù)測的時(shí)頻共振方法研究[D].西安:西北工業(yè)大學(xué),2002.

[8] 張文,劉玉田.自適應(yīng)粒子群優(yōu)化算法及其在無功優(yōu)化中的應(yīng)用[J].電網(wǎng)技術(shù),2006,30(8):19-24.

[9] 陳志敏,薄煜明,吳盤龍.基于自適應(yīng)粒子群優(yōu)化的新型粒子濾波在目標(biāo)跟蹤中的應(yīng)用[J].控制與決策,2013,28(2):193-200.

[10] 鄭華,譚博,裴承鳴.粒子濾波算法在多傳感器測量中的應(yīng)用[J].現(xiàn)代電子技術(shù),2014,37(1):24-26,30.

(編輯：姚妙慧)

A methodology of non-stationary degree suitable for FTPVS data processing

TAN Bo, ZHENG Hua, PEI Cheng-ming

(School of Power and Energy, NWPU, Xi’an 710072, China)

Aiming at qualifying of non-stationary and better analysis results, a methodology of non-stationary degree based on time-frequency analysis is proposed in this paper according to the characteristics of actual collected data, which makes it possible to modify the related parameters of processing method according to non-stationary degree. Non-stationary degree proposed in this paper is tested through both simulation data and actual experimental data. The results show the non-stationary degree in this paper could credibly indicate the non-stationary degree and could be used to improve the accuracy of processing methods for FTPVS.

non-stationary degree; surrogate data; FTPVS

2015-03-05;

2015-06-02;

時(shí)間:2015-08-17 11:04

國家自然科學(xué)基金資助(11302175)

譚博(1987-)，男，陜西西安人，博士研究生，主要研究領(lǐng)域?yàn)闀r(shí)變參數(shù)建模、非平穩(wěn)信號處理； 鄭華(1983-)，男，河南南陽人，講師，主要研究領(lǐng)域?yàn)樽赃m應(yīng)控制算法、非平穩(wěn)信號處理。

V215.34

A

1002-0853(2015)06-0560-05