任春雷,周小陳,張炳誠,薛小龍,武艷奎,楊懿
(1.國防科技大學空天科學學院,湖南 長沙 410000;2.中國航天科工飛航技術(shù)研究院,北京 100074;3.中國航發(fā)北京航科發(fā)動機控制系統(tǒng)科技有限公司,北京 102200;4.北京控制與電子技術(shù)研究所,北京 100038;5.龐巴迪(天津)航空服務(wù)有限公司,天津 300300;6.北京航天試驗技術(shù)研究所,北京 100074)
在液體火箭發(fā)動機試驗中,根據(jù)所測壓力信號隨時間的變化可分為穩(wěn)態(tài)壓力和脈動壓力兩大類。穩(wěn)態(tài)壓力的測量結(jié)果是壓力平均值,用于研究試驗臺和發(fā)動機的性能。脈動壓力則主要針對發(fā)動機推力室、泵、燃氣發(fā)生器、推進劑管路等高壓、高溫、強振動測點進行脈動壓力參數(shù)的測量與分析,其主要目的是研究發(fā)動機啟動特性、不穩(wěn)定燃燒、管路特性以及分析發(fā)動機的故障原因[1]。由于脈動壓力測量與穩(wěn)態(tài)壓力測量在測量方法、測量系統(tǒng)構(gòu)成和數(shù)據(jù)結(jié)構(gòu)等方面有較大的差別。因此采用準確、有效的方法對脈動壓力數(shù)據(jù)進行分析是非常重要的。
從時間域內(nèi)壓力變化趨勢的角度進行分析,脈動壓力與穩(wěn)態(tài)壓力最大的區(qū)別在于其變化時間序列具有非線性、隨機性的復雜特征[2]。業(yè)內(nèi)一般采用計算一段時間內(nèi)的平均值或者觀察某個時間點的壓力值的方法對穩(wěn)態(tài)壓力數(shù)據(jù)進行分析。而對于高頻采樣下的高精度脈動壓力數(shù)據(jù),則不僅需要關(guān)注其在時域內(nèi)的幅值信息,而且還需要根據(jù)信號時域和頻域的特征,結(jié)合實際研究的需要,采用時頻分析技術(shù)深度挖掘數(shù)據(jù)在時域和頻域的特性信息。目前,在信號時頻域分析技術(shù)領(lǐng)域,較為常用的方法有傅里葉變換、短時傅里葉變換、快速傅里葉變換和小波分析方法等等。這幾種時頻分析方法各具特點。
傅里葉變換和快速傅里葉變換方法建立了信號在時、頻域的相互轉(zhuǎn)換關(guān)系。通過繪制信號在頻域內(nèi)的全局圖譜,可以清晰觀察所分析的信號在時域內(nèi)的頻率-幅值特征。該方法的不足之處在于缺乏局域分析能力,即無法定位所得到的頻率-幅值特性在時域內(nèi)的發(fā)生時刻[3-4]。短時傅里葉變換的核心思想是將一個非線性信號看作若干個短的線性信號的集合,通過設(shè)定一個尺度不變的窗函數(shù),將窗函數(shù)在時間軸上進行平移,以此來覆蓋整個非線性信號,分析多個短信號局部的時間-頻率信息。該方法的不足之處在于實際應(yīng)用中需要根據(jù)信號的特征選擇對應(yīng)的窗函數(shù),且分析效果取決于所選定的窗函數(shù),分析的結(jié)果具有一定的不確定度性[5]。小波分析方法的實質(zhì)是通過分析特定時間域內(nèi)原始信號的特征,人為選定與之相匹配的基函數(shù)。將基函數(shù)在尺度上進行相應(yīng)的伸縮,建立一個函數(shù)簇。將函數(shù)簇在指定的間域內(nèi)進行多次平移,并對分解得到信息進行濾波等處理后進行重構(gòu),能夠得到原始信號在時間域內(nèi)局部特征信息。小波分析克服傅里葉變換不能兼顧時間和頻率分辨力的缺點,具有局部自適應(yīng)分析能力,在工程應(yīng)用中得到了廣泛的應(yīng)用。其局限性體現(xiàn)在以下兩點:首先,小波分析以傅里葉變換為基礎(chǔ),雖然提供了一個尺度可變的時間-頻率窗,但是過度調(diào)整窗口的尺度會產(chǎn)生信息泄露的問題。其次,由于已經(jīng)根據(jù)原始信號的特征人為選定了與之相匹配的基函數(shù),且在后續(xù)的分析過程中無法進行更改,不同的小波基函數(shù)的分析結(jié)果往往是不同的[6-7]。
為滿足精確描述頻率隨時間變化的需要,美籍華人科學家Norden E.Huang提出了一種新的信號處理方法:經(jīng)驗?zāi)B(tài)分解(empirical mode decomposition,EMD)[8]。該方法對于非平穩(wěn)線性信號的分析具有非常好的效果,其核心思想是對非平穩(wěn)線性信號進行平穩(wěn)化處理,逐級分解信號中所包含的不同尺度下的波動或者變化趨勢,將分解得到的一系列不同特征尺度的數(shù)據(jù)序列稱之為本征模態(tài)函數(shù)(Intrinsic Mode Function,IMF)。對每一個IMF進行希爾伯特變換(Hilbert Transform,HT),從而得到不同特征尺度下的瞬時振幅和瞬時頻率。該理論的創(chuàng)新之處在于定義了瞬時頻率的概念,引進了EMD和IMF的概念。通過對IMF分量進行希爾伯特變換得到時頻域的Hilbert能量分布頻譜圖,建立了譜分析與時域的對應(yīng)關(guān)系,能夠精確分析原始信號時頻域的分布規(guī)律[9-10]。
汪偉等[11]針對跑步機用戶進行行走速度估計過程中出現(xiàn)的加速度漂移問題,根據(jù)跑步機用戶行走速度數(shù)據(jù)的特征,利用EMD方法在不同頻率段內(nèi)將擺動腳速度進行分解。通過選擇與用戶步頻相近的頻段,摒棄無關(guān)頻段,并對本證模態(tài)函數(shù)進行重構(gòu),從而解決了加速度飄逸的問題。采用集成經(jīng)驗?zāi)B(tài)分解法對速度積分后得到的擺動腳位移進行分解與重構(gòu),有效解決了積分累積誤差引起的模態(tài)混疊現(xiàn)象。徐可、李俊等[12-13]充分利用EMD方法在自適應(yīng)分析方面的優(yōu)勢,采用將EMD分解后的本征模態(tài)函數(shù)所包含的能量特征和支持向量機(SVM)遺傳算法相結(jié)合的方法對非平穩(wěn)軸承振動信號處理進行分析處理,分析的結(jié)果表明該方法在識別和定位滾動軸承故障方面具有良好的效果。汪朝海等[14]利用EMD方法在主成分分析、信息提取和自適應(yīng)分析方面的優(yōu)勢,對數(shù)據(jù)分析EMD分解。將IMF函數(shù)和殘差函數(shù)組成的特征向量作為滾動軸承故障識別的主要特征。對振動數(shù)據(jù)的分析結(jié)果顯示,與支持向量機、BP神經(jīng)網(wǎng)絡(luò)等算法相比,該方法在故障識別率、運行時間和操作便捷性方面更具有優(yōu)勢。楊建華等[15]對降噪后的滾動軸承故障數(shù)據(jù)進行EMD分解。通過對比分析軸承故障仿真信號和滾動軸承實驗信號的之間的特性信息,能夠?qū)Π瑥娫肼暤臐L動軸承數(shù)據(jù)進行故障分析和定位。晏紅艷等[16]針對傳統(tǒng)相干噪聲壓縮去噪法對含噪信號進行去噪過程中容易消除有效噪聲的問題,結(jié)合多通道聯(lián)合時頻分析法,將含噪地震信號分解為多個頻率段下的分量。利用EMD方法在特征信號提取和自適應(yīng)分解方面的優(yōu)勢,在消除了含噪地震信號中噪聲的同時也有效保存了原始信號中的有用信息。
本文針對發(fā)動機試驗脈動壓力數(shù)據(jù)的非平穩(wěn)非線性特征,結(jié)合經(jīng)驗?zāi)B(tài)分解在自適應(yīng)特征分析、瞬時頻率的精確定位能力、局部瞬時表達能力以及提取信號分量等方面的優(yōu)勢,對某型號發(fā)動機試驗的脈動壓力數(shù)據(jù)進行分析。
工程應(yīng)用中,脈動壓力測量系統(tǒng)具有頻率響應(yīng)快,采集系統(tǒng)采樣率高的特征,能夠精確測量推力室高溫、高壓和劇烈振動環(huán)境下的壓力變化信息。測量系統(tǒng)一般由數(shù)采和分析系統(tǒng)、供電系統(tǒng)、脈動壓力傳感器、信號轉(zhuǎn)換裝置、循環(huán)冷卻系統(tǒng)和信號轉(zhuǎn)換裝置等幾個部分組成[8],測量系統(tǒng)組成原理如圖1所示。
圖1 脈動壓力測量系統(tǒng)組成原理圖Fig.1 Schematic diagram of fluctuating pressure measurement system composition
脈動壓力傳感器的功能主要是完成脈動壓力測點的測量,實現(xiàn)測點壓力值到電信號的轉(zhuǎn)換。傳感器必須具有耐強振動、高溫和高壓的工作特性。脈動壓力傳感器的安裝以及測壓導管的長度等都有一定的要求[8]。
數(shù)采和分析系統(tǒng)主要完成數(shù)據(jù)的采集和分析,對其采樣率和頻響也有相應(yīng)的要求,采集系統(tǒng)采樣率一般大于等于10 kHz。
循環(huán)冷卻系統(tǒng)主要功能是在傳感器工作過程中通過水冷或者氣冷的方式降低傳感器內(nèi)部環(huán)境的溫度,保障傳感器正常工作溫度環(huán)境,延長工作壽命,提高系統(tǒng)可靠性。
對于一個實信號x(t),可以構(gòu)建一個復信號z(t),定義為[10]
式中:a(t),θ(t)分別為實信號x(t)的瞬時幅值和瞬時相位。
瞬時頻率等于瞬時相位的導數(shù),即
令EMD分解得到的IMF分量為hi(t),對其進行希爾伯特變換為
根據(jù)公式構(gòu)造解析信號
則幅值函數(shù)為
式中:ai(t)為解析信號的瞬時幅值;φi(t)為解析信號的瞬時相位。
瞬時頻率為
希爾伯特譜為
希爾伯特邊際譜為
式中:T為所分析原始信號的總時長。從式(9)可知,希爾伯特譜能夠在時、頻域內(nèi)精確描述信號的幅值的變化規(guī)律。
工程應(yīng)用中對EMD分解后的征模態(tài)函數(shù)(IMF)分量進行HT變換時,IMF分量在數(shù)學上需要滿足時間歷程中極值點個數(shù)和過零點個數(shù)相等或者相差1個和局部極大值、極小值相加的均值為0的要求[17]。雖然IMF函數(shù)能夠合理地反映信號內(nèi)部在每一個周期上的波動性,但是并非在所有情況下滿足上述要求的IMF函數(shù)都能夠計算出瞬時頻率。其主要的原因在于工程應(yīng)用中原始信號往往較為復雜。真實信號中包含大量不相關(guān)的噪聲分量使得原始信號中包含多個振動模態(tài)分量。針對該問題,Huang和相關(guān)研究人員對IMF函數(shù)及相應(yīng)原始信號的特征進行了研究和修訂,提出符合對信號進行EMD分解及對IMF函數(shù)進行HT變化的原始信號需要滿足以下假定條件[17]:①原始信號至少包含一個極大值點和一個極小值點;②相鄰兩個極值點之間的時間間隔為特征時間尺度;③對于信號中無極值點,僅有拐點的情況,采用對信號進行微分的方法獲得極值點。
EMD分解的步驟主要分為:
1)定位x(t)中所有極大值點和極小值點。采用三次樣條插值法繪制x(t)的上下包絡(luò)線(極大值點形成上包絡(luò)線,極小值點形成的下包絡(luò)線)。計算上、下包絡(luò)線的均值m1。x(t)與m1的差記為h1,即
將h1作為新信號x(t),重復式(11)的步驟。直至所得到的分量hi滿足IMF函數(shù)的基本條件[17]時,將滿足條件的第1個一階IMF分量稱之為C1,C1包含信號中的最高頻成分。
2)將一階IMF分量C1從x(t)中分離出來,令剩余的原始信號為r1,令
將r1作為新信號,重復步驟1),直到第n階的殘余信號成為單調(diào)函數(shù)。
3)數(shù)學上x(t)可視為n個IMF分量和一個殘余分量的和,即
式中:rn(t)為信號中的殘量。各IMF分量Cj(t)則代表了信號中所包含的所有不同高、低頻率的成分。
在實際操作過程中需利用計算機編程完成信號的EMD分解過程。分解流程一般如圖2所示。
2.2.1 自適應(yīng)時頻分析
從2.1章節(jié)的基本理論可知,EMD理論建立了瞬時頻率的物理意義。通過將原始信號分解為多個IMF分量,對各IMF分量進行HT變換得到原始信號在時、頻域的希爾伯特譜。瞬時特征量精確定位了希爾伯特譜在時頻域的特征。相較傳統(tǒng)的時頻分析方法,EMD方法的優(yōu)點在于定義了一組根據(jù)原始信號特征而制定的具有自適應(yīng)分解特性的基函數(shù),從而使該方法在實際應(yīng)用中具有良好的自適應(yīng)分析優(yōu)勢。
2.2.2 信號局部瞬時特征表征
根據(jù)2.1章節(jié)的理論可知,通過給出瞬時頻率的定義,即對EMD分解過程中所得到的多個IMF分量的相位進行求導,能夠從物理意義層面對瞬時頻率在局部瞬時特性方面的優(yōu)勢進行清晰、明確地說明。相較傳統(tǒng)傅里葉變換在頻率混疊、虛假諧波方面的缺陷,EMD分解中多個IMF分量從高階到低階的排序也有效解決了上述問題。多個IMF分量瞬時頻率相位求導的方法則體現(xiàn)了該方法具有對原始信號的局部特征進行精確描述的能力。
2.2.3 主成分分析
在工程應(yīng)用中,復雜的非線性非平穩(wěn)隨機信號往往包含了大量真實信號和各種不相關(guān)的噪聲分量。對復雜信號進行分析的重要目的是將復雜原始信號中有價值的真實信號提取出來。從2.1章節(jié)的介紹可知,EMD方法能夠?qū)⒃夹盘栔械念l率分量按照從高到低的順序分離提取出來以滿足信號分析的需要。數(shù)據(jù)分解的結(jié)果能夠較好地體現(xiàn)各頻率分量在原始信號中的比重。因此,該方法也是一種能夠?qū)碗s信號進行主成分分析、提取的重要工具。
基于脈動壓力測量系統(tǒng)高頻響、高采樣率的特性,脈動壓力數(shù)據(jù)具有采樣率高、數(shù)據(jù)容量大的顯著特征[18]。某型號發(fā)動機試車中,按照圖1組建脈動壓力參數(shù)測量系統(tǒng)。傳感器采用Kulite-CT-375水冷型脈動壓力傳感器。傳感器的性能參數(shù)和計量機構(gòu)檢定結(jié)論分別如表1、表2所示。
表1 CT-375型傳感器性能參數(shù)Tab.1 Performance parameters of CT-375 sensor
表2 CT-375型傳感器檢定結(jié)果Tab.2 Verification results of CT-375sensor
推力室脈動壓力參數(shù)名為Pis1,脈動參數(shù)采集系統(tǒng)的采樣率為10 kHz。選取試車中穩(wěn)定段0.2 s(2000個樣本點)的數(shù)據(jù)進行分析。試驗全程脈動壓力數(shù)據(jù)如圖3所示。相鄰測量部位穩(wěn)態(tài)壓力測量數(shù)據(jù)如圖4所示。測點名:Pcio1,穩(wěn)態(tài)參數(shù)采集系統(tǒng)的采樣率為1 kHz。由于脈動壓力測點和穩(wěn)態(tài)壓力測點分屬不同的測量系統(tǒng),所用數(shù)據(jù)采集設(shè)備不同,故二者在坐標軸上顯示的時間有所區(qū)別,但發(fā)動機總工作時間是相同的。
從圖3和圖4可知,脈動壓力測點數(shù)據(jù)和穩(wěn)態(tài)壓力測點數(shù)據(jù)在波形變化趨勢上基本一致,印證了測量數(shù)據(jù)的合理性。相較高采樣率、高精度脈動壓力測量系統(tǒng),穩(wěn)態(tài)壓力測量系統(tǒng)的測量數(shù)據(jù)在波形上較為光滑。而脈動壓力測量系統(tǒng)憑借高采樣率、高頻響和高測量精度的優(yōu)勢,能夠非常精確地捕捉到推力室內(nèi)壓力的波動變化情況。
圖3 某型號發(fā)動機試驗脈動壓力全程數(shù)據(jù)圖Fig.3 Full range data diagram of fluctuating pressure of a certain type of engine test
圖4 某型號發(fā)動機試驗穩(wěn)態(tài)壓力全程數(shù)據(jù)圖Fig.4 Full range data diagram of steady-state pressure of a certain type of engine test
脈動壓力數(shù)據(jù)穩(wěn)定段0.2 s和0.02 s局部數(shù)據(jù)分別如圖5和圖6所示。圖5中,在發(fā)動機試車穩(wěn)定段26.10~26.30 s內(nèi),推力室壓力在8.614~8.911 MPa之間震蕩波動,變化范圍達到0.297 MPa。圖6清晰地顯示了推力室內(nèi)壓力突變的情況,具有高采樣率特征的脈動壓力測量數(shù)據(jù)在精確描述推力室內(nèi)脈動壓力在時域內(nèi)波動情況的同時,也為工程技術(shù)人員采用多種時頻分析方法對其進行分析奠定了重要的基礎(chǔ)。
圖5 發(fā)動機試驗脈動壓力局部數(shù)據(jù)圖1Fig.5 Local data diagram 1 of engine test fluctuating pressure
圖6 發(fā)動機試驗脈動壓力局部數(shù)據(jù)圖2Fig.6 Local data diagram 2 of engine test fluctuating pressure
根據(jù)2.1章節(jié)的理論基礎(chǔ)和方法,對穩(wěn)定段26.10~26.30 s內(nèi)的脈動壓力數(shù)據(jù)進行經(jīng)驗?zāi)B(tài)分解,分解結(jié)果如圖7所示。
圖7 EMD分解結(jié)果圖Fig.7 EMD decomposition result
從圖7中EMD分解結(jié)果可知,采用EMD分解方法能夠?qū)⒚}動壓力數(shù)據(jù)中的頻率分量分離出來。IMF1~IMF4階分量分別顯示了脈動壓力數(shù)據(jù)所在時域內(nèi)所包含的頻率成分以及幅值相對大小。其中IMF1階分量在原始信號中所占的比重最大,頻率分量最高。IMF2~IMF4階分量在原始信號中所占比例逐漸降低。IMF2,IMF3和IMF4分量的幅值波動趨勢要大于IMF1分量的幅值波動。R主要為系統(tǒng)中的一些其他不相關(guān)低頻分量,變化過程非常緩慢且無規(guī)律,對數(shù)據(jù)的分析基本無影響。
對IMF函數(shù)進行HT變換的希爾伯特譜如圖8所示。
圖8 IMF1~IMF4時間頻率曲線圖Fig.8 IMF1~IMF4 time frequency curve
由圖8可知第1層~第4層時間-頻率曲線非常清晰準確地顯示了全時間域內(nèi)頻率的變化情況。在第1,2層中,3500~5000 Hz左右的中頻帶是數(shù)據(jù)的主要組成成分,1500 Hz以下的低頻帶所占的比重較少。第2,3和4層時間-頻率曲線頻率值變化情況顯示在該時域內(nèi)0.03,0.048,0.155和0.188 s存在多個頻率值突變。結(jié)合推力室上下游相應(yīng)的溫度、壓力測點、試驗時序和試驗工況等信息,對于分析判斷推力室內(nèi)燃燒狀態(tài)具有非常重要的作用。
對分解得到的IMF分量進行重構(gòu)后的時間-幅值特性曲線如圖9所示。脈動壓力數(shù)據(jù)的FFT頻譜分析圖如圖10所示。
結(jié)合圖5、圖8、圖9和圖10進行分析可得到以下幾點結(jié)論:
圖9 EMD重構(gòu)時間-幅值圖Fig.9 EMD reconstruction time-amplitude diagram
圖10 FFT頻譜分析圖Fig.10 FFT spectrum analysis diagram
1)通過對比圖5與圖9可知,圖9中重構(gòu)后的時間-幅值特性曲線相比原始數(shù)據(jù)具有更高的分辨力,能夠更加清晰地分辨時間域內(nèi)的幅值波動情況。不僅可以清晰分辨脈動壓力周期性波動的情況,而且可以精準定位在0.12,0.18 s的幅值振蕩。說明在該試車的穩(wěn)定段內(nèi),隨著推進劑燃燒對壓強振蕩響應(yīng),燃燒產(chǎn)生的能量持續(xù)注入燃燒室工作系統(tǒng)引起脈動壓力有規(guī)律地振蕩。通過對比同時段推力室對應(yīng)的振動測點的振動數(shù)據(jù),脈動壓力有規(guī)律的振蕩在振動數(shù)據(jù)中也得到了體現(xiàn)。此外,根據(jù)EMD分解的特點,針對渦輪泵、燃發(fā)器等振動劇烈的測點,還可參考文獻[8],[18]中的方法,結(jié)合相應(yīng)振動測點的數(shù)據(jù)進行分析。
2)圖8中第1層和第2層的時間頻率特性曲線與圖10中FFT分析結(jié)果相吻合。圖8中第3層~第4層時間-頻率特性曲線非常清晰地顯示了低頻分量的變化相對較緩慢的特性,也能非常清晰定位對應(yīng)的變化時間點。
3)由圖10中FFT分析結(jié)果可知,全時域內(nèi)幅值峰值多集中在低頻區(qū)(0~2000 Hz)和中頻區(qū)(3500~5500 Hz)。這與圖8中EMD分解時間-頻率曲線的結(jié)果相吻合。但FFT的不足之處是分析結(jié)果只能顯示全時域數(shù)據(jù)在頻域內(nèi)的信息,無法定位某個幅值峰值所發(fā)生的時間信息,也無法分析信號中的各頻率成分。
分別采用Haar小波和Db9小波對脈動壓力數(shù)據(jù)進行7層分解的結(jié)果如圖11至圖14所示。
圖11 Haar小波5~7層分解細節(jié)系數(shù)圖Fig.11 Detail coefficient diagram of Haar wavelet 5~7 level decomposition
圖12 Haar小波1~4層分解細節(jié)系數(shù)圖Fig.12 Detail coefficient diagram of Haar wavelet 1~4 level decomposition
圖13 Db9小波5~7層分解細節(jié)系數(shù)圖Fig.13 Detail coefficient diagram of Db9 wavelet 5~7 level decomposition
圖14 Db9小波1~4層分解細節(jié)系數(shù)圖Fig.14 Detail coefficient diagram of Db9 wavelet 1~4 level decomposition
由小波分析圖11至圖14可知,小波分析的結(jié)果與所選小波基函數(shù)有關(guān)(文獻[18]介紹了多種小波基函數(shù)的特征)。由于Haar小波具有緊支撐、單個矩形波、矩形狀階梯變化等特征,分解得到的多層細節(jié)分量具有矩形特征,尤其是d4~d7層細節(jié)信息尤為明顯。而Db系小波具有對稱性好、頻域的局部化能力強和消失矩光滑性好等特點,分解得到的多層細節(jié)分量的光滑性較好。d1~d5層分解信號中均能定位原始信號中幅值產(chǎn)生突變的時間信息。d6~d7層的分解結(jié)果均因分解層數(shù)高,過濾掉了過多的有用信息而存在一定程度的失真。小波分析的結(jié)果可以結(jié)合相應(yīng)的振動測點數(shù)據(jù),進一步分析所測部件的動態(tài)特性[8]。小波基函數(shù)的選擇與分解層數(shù)的差異導致小波分析所得到的結(jié)果存在一定的差異,也給準確分析數(shù)據(jù)特征帶來了困難。
本文介紹了型號發(fā)動機試驗中脈動壓力測量系統(tǒng)的組成和數(shù)據(jù)特征。分析了采用傳統(tǒng)傅里葉變換、小波分析等方法對非線性非平穩(wěn)隨機數(shù)據(jù)進行時頻分析的缺點和不足。重點介紹了EMD的基本理論。從自適應(yīng)時頻分析、局部瞬時定位分析和信號主成分特征分析和提取三個方面介紹了該理論的優(yōu)點以及實際工程應(yīng)用的詳細步驟。分別采用FFT方法、基于不同小波基函數(shù)和分解層數(shù)的小波分析方法和EMD分解方法對某型號發(fā)動機試驗脈動壓力數(shù)據(jù)進行分析。對比分析了FFT方法和小波分析方法對非平穩(wěn)非線性數(shù)據(jù)的不足。脈動壓力數(shù)據(jù)的EMD分解結(jié)果顯示該方法能夠非常清晰地自適應(yīng)分離原始信號中的高、低頻分量,顯示數(shù)據(jù)在時間-頻率域的多維特征。該方法可以與其他數(shù)據(jù)分析方法相結(jié)合,為產(chǎn)品型號設(shè)計和試驗人員分析非線性非平穩(wěn)型數(shù)據(jù)提供重要的參考依據(jù),具有重要的參考和推廣價值。