基于奇異值分解的液壓信號(hào)時(shí)域分解方法

張小明，唐 建，張梅軍

(中國(guó)人民解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，南京 210007)

基于奇異值分解的液壓信號(hào)時(shí)域分解方法

張小明，唐 建，張梅軍

(中國(guó)人民解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，南京 210007)

為了抑制模態(tài)混疊和降低分量中的噪聲含量，提出了一種基于奇異值分解的液壓信號(hào)時(shí)域分解方法。根據(jù)奇異值分解的兩點(diǎn)特性：①每個(gè)頻率成分對(duì)應(yīng)兩個(gè)大小相當(dāng)?shù)钠娈愔担虎诟黝l率對(duì)應(yīng)的奇異值的大小與該頻率的振幅呈正相關(guān)，該方法先選取原信號(hào)中的某一頻率，向其中疊加頻率相同、振幅已知的周期信號(hào)，使疊加信號(hào)中該頻率的振幅最大，這樣與其對(duì)應(yīng)的奇異值一定位于對(duì)角矩陣的前兩階，解決了原信號(hào)該頻率的奇異值階數(shù)無(wú)法確定的問(wèn)題，繼而選取前兩階奇異值重構(gòu)，再減去前步加入的周期信號(hào)，即還原出原信號(hào)中該頻率的時(shí)間序列。同樣，對(duì)于原信號(hào)中的其他頻率用相同方法處理，最終獲得一組分量。經(jīng)實(shí)驗(yàn)，該方法較EMD不僅能有效消除模態(tài)混疊，而且降低了分量中的噪聲含量。

奇異值分解；構(gòu)造信號(hào)；疊加信號(hào)；模態(tài)混疊；噪聲含量

液壓信號(hào)可以看成是由若干個(gè)不同頻率的頻率族分量組成，而每一個(gè)頻率族的分量所對(duì)應(yīng)的幅值可以看成是調(diào)幅信號(hào)[1]，將這些成分單獨(dú)地分開(kāi)是后續(xù)處理過(guò)程的關(guān)鍵步驟，在將液壓信號(hào)分解到不同的頻率通道中的研究中，主要有以下兩種方法：

(1) 小波分解及小波包分解。不同于小波分解僅對(duì)低頻部分繼續(xù)分解，小波包分解兼對(duì)低頻和高頻部分實(shí)施分解，以獲得更高的頻率分辨率。當(dāng)需要獲得信號(hào)在一個(gè)頻段的時(shí)間序列的時(shí)候，可以重構(gòu)該頻段的分解系數(shù)獲得。但是限于怎樣合適選取小波基函數(shù)及特定基函數(shù)對(duì)于特定信號(hào)是否適用，使得小波分解及小波包分解的應(yīng)用于實(shí)測(cè)信號(hào)存在困難。

(2) EMD(Empirical Mode Decomposition)及其改進(jìn)方法。EMD是由Huang等[2-3]提出的一種信號(hào)的時(shí)域分解方法，信號(hào)通過(guò)EMD被分解成一系列內(nèi)稟模態(tài)函數(shù) (Intrinsic Mode Function， IMF)，這些IMF含有特定的頻率成分。在處理實(shí)測(cè)信號(hào)時(shí)，EMD主要存在兩點(diǎn)缺陷：

①模態(tài)混疊 (Mode Aliasing)。EMD分解的最理想的效果是每個(gè)IMF分量中僅含有一種頻率成分。而實(shí)測(cè)信號(hào)在分解時(shí)常會(huì)出現(xiàn)單個(gè)IMF分量中存在多個(gè)頻率成分，或者單個(gè)頻率成分會(huì)發(fā)散到不同的IMF分量中的現(xiàn)象，即模態(tài)混疊。為了改善模態(tài)混疊，Wu等[4]提出的EEMD (Ensemble Empirical Mode Decomposition)通過(guò)向原信號(hào)中加入多組白噪聲分解后再求平均的方法來(lái)抑制該缺陷，但該方法面臨著如何合適選取噪聲組數(shù)和方差的困難；Ryan等[5]提出掩膜信號(hào)法，通過(guò)對(duì)原信號(hào)分別加上和減去一個(gè)掩膜信號(hào)分別進(jìn)行EMD分解，再對(duì)IMF分量做平均，趙玲等[6]對(duì)平均瞬時(shí)頻率進(jìn)行修正來(lái)改進(jìn)掩膜信號(hào)，但是針對(duì)不同的信號(hào)需要確定不同的掩膜信號(hào)，而且用希爾伯特變換估計(jì)出的瞬時(shí)頻率偏差比較大，該方法對(duì)分離振幅和頻段比較接近的信號(hào)效果并不好。

②IMF分量中噪聲含量較高。分解得到的IMF分量中有的分量的幅值譜中的波峰向兩邊發(fā)散較大，分量中參雜了不必要的噪聲成分。實(shí)際操作時(shí)，為了有效降低IMF分量中的噪聲含量，往往對(duì)實(shí)測(cè)信號(hào)先降噪，再分解。

奇異值分解 (Singular Value Decomposition， SVD)是一種通過(guò)矩陣運(yùn)算對(duì)信號(hào)進(jìn)行處理的方法，根據(jù)Hankel矩陣構(gòu)造形式將一維信號(hào)構(gòu)造成信號(hào)矩陣，該矩陣通過(guò)奇異值分解被分解成兩個(gè)酉矩陣和一個(gè)對(duì)角陣的乘積，該方法在信號(hào)處理尤其是信號(hào)降噪中有著廣泛的使用[7-12]。其中最具代表性的是趙學(xué)智等人提出的奇異值差分譜降噪方法，該方法合理地選擇了有效重構(gòu)奇異值的階數(shù)，很好地從含噪信號(hào)中還原出理想信號(hào)。

根據(jù)奇異值分解的相關(guān)特性，本文提出了一種液壓信號(hào)的時(shí)域分解方法。該方法根據(jù)原信號(hào)中某一頻率的構(gòu)造一相同頻率、振幅已知的周期信號(hào)，與原信號(hào)疊加，使得疊加信號(hào)中該頻率對(duì)應(yīng)的奇異值所在階數(shù)可以確定，繼而選取這些階數(shù)上奇異值重構(gòu)出疊加信號(hào)中該頻率的時(shí)間序列，再減去加入的周期信號(hào)，最終還原出原信號(hào)中該頻率的時(shí)間序列。在與EMD進(jìn)行比較后，發(fā)現(xiàn)該方法不僅有效消除了模態(tài)混疊，而且降低了分量中的噪聲含量。

1 奇異值分解及其差分譜理論

1.1 奇異值分解理論

實(shí)測(cè)信號(hào)x可以根據(jù)Hankel矩陣形式被構(gòu)造成一個(gè)m×n的信號(hào)矩陣Hx，矩陣構(gòu)造形式為

(1)

式中：N為x的長(zhǎng)度；1

對(duì)Hx進(jìn)行奇異值分解

Hx=USVT

(2)

式中：U為m×m的酉矩陣；S為m×n的對(duì)角矩陣；V為n×n的酉矩陣。其中S可表示為

(3)

式中：∑=diag(σ1,σ2,…,σr) (σ1>σ2>…>σr>0)。

這種變換的本質(zhì)是將原始信號(hào)分解為一系列分量信號(hào)的簡(jiǎn)單線性疊加，而這一系列分量可以通過(guò)單獨(dú)選取S中的對(duì)角線上的單個(gè)奇異值重構(gòu)獲得，而且重構(gòu)出的信號(hào)的相位和它在原信號(hào)中疊加的相位是不變的，即具有零相位偏移特性。

1.2 奇異值差分譜理論

奇異值分解被有效用于信號(hào)的降噪中，在文獻(xiàn)[13]中，定義差分值bi來(lái)描述σi的變化情況

bi=σi-σi+1(1

(4)

式中，min(m，n)為Hankel矩陣行數(shù)和列數(shù)中的較小值。當(dāng)bi達(dá)到最大值時(shí)，此時(shí)奇異值σi所代表的分量實(shí)現(xiàn)了從有用信號(hào)分量向噪聲分量的轉(zhuǎn)變，接下來(lái)選取該點(diǎn)及其之前的奇異值為有效奇異值，其他奇異值置零，然后通過(guò)式(2)重構(gòu)出降噪后的信號(hào)矩陣，再重排矩陣元素獲得降噪信號(hào)。

2 基于奇異值分解的信號(hào)時(shí)域分解方法

2.1 方法的提出

本文提出的方法旨在消除分量中模態(tài)混疊，這就要求待分解的原信號(hào)必須是平穩(wěn)信號(hào)，因?yàn)槠椒€(wěn)信號(hào)中的頻率具有可數(shù)性和時(shí)不變性的特點(diǎn)，這就給原信號(hào)中不同頻率能夠被單獨(dú)分開(kāi)提供了可能；同時(shí)原信號(hào)的截取應(yīng)滿足采樣定理和足夠的采樣長(zhǎng)度，這樣生成的序列才能重現(xiàn)振動(dòng)波形中的所有頻率成分及獲得較高的頻率分辨率。文章開(kāi)始提到液壓信號(hào)可以看成由不同頻率族分量的調(diào)幅信號(hào)組成，表示液壓信號(hào)同樣是平穩(wěn)信號(hào)，適合用本文提出的方法進(jìn)行分解。

通過(guò)錢征文等[14]的研究和相關(guān)實(shí)驗(yàn)，總結(jié)出奇異值分解在處理平穩(wěn)信號(hào)時(shí)存在以下兩點(diǎn)性質(zhì)：

(1) 當(dāng)原信號(hào)中除噪聲外僅含有一個(gè)頻率成分，該信號(hào)的Hankel矩陣經(jīng)分解后得到的對(duì)角陣S的前兩階的奇異值較大，而其他階奇異值很??；當(dāng)除噪聲外含有兩個(gè)頻率成分，分解后S的前4階的奇異值較大，而其他階奇異值很小。依此類推：當(dāng)原信號(hào)除噪聲外含有k個(gè)頻率成分，經(jīng)分解后S的前2k階奇異值較大，其他階奇異值較小。而對(duì)于信噪比較高的信號(hào)，仍然滿足上述規(guī)律。

(2) 成對(duì)出現(xiàn)的奇異值的大小與該頻率的振幅呈正相關(guān)，即某頻率的振幅越大，經(jīng)分解后在S中得到的奇異值越大，它們?cè)赟中所處的階數(shù)越靠前。而奇異值大小與頻率大小無(wú)明顯關(guān)系。

現(xiàn)構(gòu)造y1對(duì)以上性質(zhì)加以說(shuō)明，y1中含有4個(gè)分量，分量1是頻率10 Hz，振幅為1的正弦分量；分量2是頻率為20 Hz，振幅為2的正弦分量；分量3是頻率為30 Hz，振幅為3的正弦分量；分量4是均值為0，方差為1的白噪聲。y1的采樣頻率為1 000 Hz，含有1 024個(gè)采樣點(diǎn)。圖1是y1及其各分量的時(shí)域圖。

y1=sin(20πt)+2sin(40πt)+3sin(60πt)+wn

(5)

圖1 y1及各分量的時(shí)域圖Fig.1 y1’s components’ time-domain plots

構(gòu)造y1的Hankel矩陣并對(duì)其進(jìn)行奇異值分解，獲得U1、S1和V1。畫出S1對(duì)角線上奇異值隨其所在階數(shù)的變化火柴梗圖，如圖2所示。

圖2中可以發(fā)現(xiàn)，前6階奇異值(772.34、762.55、520.89、506.16、256.98、250.05)的數(shù)值比較大，為有用信號(hào)的奇異值，而之后的奇異值很小，為噪聲成分的奇異值。而且前6階奇異值可以很明顯地分成3組：772.34和762.55為一組；520.89和506.16為一組；256.98和250.05為一組。由前述的奇異值分解的性質(zhì)可得：這3組奇異值由3個(gè)頻率成分分解獲得，而且這3個(gè)頻率成分的振幅依次遞減。再參照的的頻率組成可以確定：772.34和762.55是30Hz分量的奇異值，520.89和506.16時(shí)20Hz分量的奇異值，256.98和250.05為10Hz分量的奇異值。

圖2 奇異值-階數(shù)變化圖Fig. 2 Singular value-order plot

分別選取前3對(duì)奇異值重構(gòu)出3個(gè)時(shí)間序列，圖3為3個(gè)序列時(shí)域圖及幅值譜。從圖3中可以明顯看出3個(gè)序列各自對(duì)應(yīng)著y1中的30Hz、20Hz和10Hz的頻率成分，和前面的結(jié)論相吻合。

從上述過(guò)程可以看出，奇異值大小對(duì)振幅的分辨是敏感的，假設(shè)原信號(hào)中某個(gè)頻率成分的振幅最大，且在數(shù)值上較其他頻率的振幅占有絕對(duì)優(yōu)勢(shì)，那么該成分對(duì)應(yīng)的奇異值一定位于對(duì)角陣S的前兩階，此時(shí)只要選取S的前兩階就可以重構(gòu)出該頻率的時(shí)間序列。但是實(shí)測(cè)信號(hào)中也可能存在振幅相當(dāng)?shù)念l率成分，這樣它們對(duì)應(yīng)的奇異值大小也就難以區(qū)分。

圖3 3組奇異值重構(gòu)出的時(shí)間序列的時(shí)域圖及幅值譜Fig.3 Time-domain plots and amplitude spectrums of three pairs singular values

擬定任意平穩(wěn)信號(hào)x的分解步驟

步驟1 從x的時(shí)域圖上得到它的振幅A，從幅值譜上準(zhǔn)確得到其頻率成分f1、f2、…、fk(f1

(6)

(7)

步驟4 構(gòu)造xnew的Hankel矩陣Hnew，對(duì)其進(jìn)行奇異值分解，獲得兩個(gè)酉矩陣Unew和Vnew，一個(gè)對(duì)角矩陣Snew。

步驟5 選取Snew中的前兩階奇異值重構(gòu)出一個(gè)信號(hào)xb1，xb1是疊加信號(hào)xnew中f1對(duì)應(yīng)的時(shí)間序列。

(8)

步驟7 針對(duì)f2，重復(fù)步驟1～步驟6，最終得到原信號(hào)中fi(i=1,2,…,k)對(duì)應(yīng)的時(shí)間序列xi(i=1,2,…,k)。

2.2 基于仿真信號(hào)的性能檢驗(yàn)

實(shí)驗(yàn)發(fā)現(xiàn)，EMD在分解含有頻率成分較接近的信號(hào)時(shí)易產(chǎn)生模態(tài)混疊，因此構(gòu)造信號(hào)y2(見(jiàn)圖4)。

圖4 y2的時(shí)域圖Fig. 4 y2’s time-domain plots

y2中含有10Hz、12Hz和14Hz的頻率成分，且各成分振幅均為1。該信號(hào)含1 024個(gè)采樣點(diǎn)，采樣頻率為1 000Hz。

y2=sin(20πt)+sin(24πt)+sin(28πt)

(9)

先對(duì)y2進(jìn)行EMD分解，圖5為IMF分量的時(shí)域圖和幅值譜。從圖5(b)中明顯地看出IMF1和IMF2中出現(xiàn)嚴(yán)重的模態(tài)混疊。

再對(duì)y2用提出方法進(jìn)行分解，先從時(shí)域圖上得到y(tǒng)2的振幅為A2=2.93。據(jù)式(6)構(gòu)造一個(gè)頻率為10Hz

(10)

圖6為xi、x2、x3的時(shí)域圖和幅值譜。從圖6(b)中可以看出每個(gè)幅值譜明顯只含有一種頻率成分，不存在模態(tài)混疊。且各頻率依次為10Hz、12Hz、14Hz，和原信號(hào)中的頻率成分吻合。且圖6(a)中xi的振幅都穩(wěn)定在1處，和原信號(hào)中的各頻率振幅一致。

圖5 IMF分量的時(shí)域圖和幅值譜Fig.5 IMF’s time-domain plots and amplitude spectrums

圖6 xi的時(shí)域圖和幅值譜Fig.6 xi’s time-domain plots and amplitude spectrums

2.3 基于液壓信號(hào)的性能檢驗(yàn)

實(shí)驗(yàn)所用的液壓綜合試驗(yàn)臺(tái)上分布著泄漏控制閥、調(diào)速閥、氣穴控制閥和堵塞控制閥等。通過(guò)不同閥門的是否開(kāi)合，開(kāi)合大小的控制，來(lái)實(shí)現(xiàn)液壓系統(tǒng)不同工況的控制。圖7為實(shí)驗(yàn)現(xiàn)場(chǎng)圖和傳感器布置圖。實(shí)驗(yàn)采用121A100型壓電式IEPE加速度傳感器通過(guò)EMT690D信號(hào)采集儀與電腦連接。如圖7(b)中A點(diǎn)所示，傳感器通過(guò)磁力吸附在液壓缸的水平方向，來(lái)測(cè)取物塊提升過(guò)程中液壓缸的振動(dòng)信號(hào)。

圖7 實(shí)驗(yàn)相關(guān)圖片F(xiàn)ig.7 Experiment related pictures

分別測(cè)得液壓缸在堵塞，泄漏，氣穴工況的振動(dòng)信號(hào)，y3，y4，y5是從3個(gè)信號(hào)中截取出的1 024個(gè)連續(xù)信號(hào)點(diǎn)，采樣頻率為1 000Hz，它們的時(shí)域圖如圖8所示。

(a) 堵塞信號(hào)y3

(b) 泄漏信號(hào)y4

(c) 氣穴信號(hào)y5圖8 不同狀態(tài)下信號(hào)的時(shí)域圖Fig.8 Signals of different status in time domain

2.3.1 堵塞信號(hào)y3的分解

(1) 本文提出的方法

作y3的幅值譜如圖9所示。

圖9 y3的幅值譜Fig.9 y3’s amplitude spectrum

從圖9中精確地得出幅值譜中的明顯頻率成分，分別為19.53Hz、38.09Hz、57.62Hz和173.80Hz，對(duì)應(yīng)圖9中的f1、f2、f3、f4。再?gòu)膱D8(a)中得到y(tǒng)3的振幅A3=2.28。

(11)

圖10中可以看出y3被單獨(dú)分解成了4個(gè)分量，不存在模態(tài)混疊。4個(gè)頻率分別是19.53 Hz、38.09 Hz、57.62 Hz和173.80 Hz。和圖9中y3的幅值譜中的主要頻率成分吻合，而且每個(gè)頻率的幅值和圖9中各頻率幅值相當(dāng)，每個(gè)頻率被極大地提取出來(lái)。同時(shí)各分量的頻域譜線比較平滑，噪聲含量比較小。

(2) EMD分解

對(duì)y3先用小波閾值降噪后再進(jìn)行EMD分解，得到一組IMF分量。圖11為IMF分量的時(shí)域圖和幅值譜。從圖11(b)中可以看出：①IMF1和IMF2發(fā)生了模態(tài)混疊；②IMF分量的頻域圖比較散亂，分量中噪聲含量較高。

圖10 xi的時(shí)域圖和幅值譜Fig.10 xi’s time-domain plots and amplitude spectrums

圖11 IMF分量的時(shí)域圖和幅值譜Fig.11 IMF’s time-domain plots and amplitude spectrums

2.3.2 泄漏信號(hào)y4的分解

(1) 本文提出的方法進(jìn)行處理。

作y4的幅值譜如圖12所示。

圖12 y4的幅值譜Fig.12 y4’s amplitude spectrum

從圖12中精確地得出幅值譜中的明顯頻率成分，分別為19.53 Hz、39.06 Hz、57.62 Hz、77.15 Hz、154.30 Hz和173.80 Hz，對(duì)應(yīng)圖12中的f1、f2、f3、f4、f5、f6。再?gòu)膱D8(a)中得到y(tǒng)4的振幅A4=1.63。用提出的方法對(duì)y4進(jìn)行分解，得到分量xi(i=1,2,…,6)。圖13為xi的時(shí)域圖和幅值譜。從圖中看出，經(jīng)分解后y4中的6個(gè)頻率成分被單獨(dú)分開(kāi)，不存在模態(tài)混疊，且各頻率成分保持了原信號(hào)中的幅值，分量的頻域譜線仍比較光滑，噪聲含量較小。

(2) EMD分解。

對(duì)y3先用小波閾值降噪后再進(jìn)行EMD分解，得到一組IMF分量。圖14為IMF分量的時(shí)域圖和幅值譜。從圖14(b)中可以看出IMF1和IMF2中存在著模態(tài)混疊，同時(shí)分量的頻域譜線仍比較雜亂，噪聲含量較高。

圖13 xi的時(shí)域圖和幅值譜Fig.13 xi’s time-domain plots and amplitude spectrums

圖14 IMF分量的時(shí)域圖和幅值譜Fig.14 IMF’s time-domain plots and amplitude spectrums

2.3.3 氣穴信號(hào)的分解

(1) 本文提出的方法進(jìn)行處理。

作y5的幅值譜如圖15所示。

圖15 y5的幅值譜Fig.15 y5’s amplitude spectrum

從圖15中精確地得出幅值譜中的明顯頻率成分，分別為3.91 Hz、15.93 Hz、57.62 Hz、77.15 Hz和99.68 Hz，

對(duì)應(yīng)圖15中的f1、f2、f3、f4、f5。再?gòu)膱D8(a)中得到y(tǒng)5的振幅A4=1.94。用提出的方法對(duì)y5進(jìn)行分解，得到分量xi(i=1,2,…,5)。圖16為xi的時(shí)域圖和幅值譜。經(jīng)分解，y5中的5個(gè)頻率成分被單獨(dú)分開(kāi)，不存在模態(tài)混疊，且各頻率成分保持了原信號(hào)中的振動(dòng)幅值，同樣各分量中噪聲含量較小。

(2) EMD分解。

對(duì)y5先用小波閾值降噪后再進(jìn)行EMD分解，得到一組IMF分量。圖17為IMF分量的時(shí)域圖和幅值譜。從圖17(b)中可以看出IMF1和IMF2中存在著模態(tài)混疊，同時(shí)分量的頻域譜線仍比較雜亂，噪聲含量較高。

圖16 xi的時(shí)域圖和幅值譜Fig.16 xi’s time-domain plots and amplitude spectrums

圖17 IMF分量的時(shí)域圖和幅值譜Fig.17 IMF’s time-domain plots and amplitude spectrums

綜合液壓系統(tǒng)的三種工況的分解結(jié)果，將各工況下的原信號(hào)和分量中頻率及其幅值匯總至表1。

從表1中可以看出：

(1) 原信號(hào)經(jīng)分解后，分量中的頻率幅值都略有減小(用下降箭頭標(biāo)出)，這是因?yàn)榉纸鈺r(shí)每次僅選取前兩個(gè)較大的奇異值重構(gòu)，而代表噪聲的小奇異值沒(méi)有參與重構(gòu)，噪聲一直被壓制在剩余信號(hào)中，從而減小了噪聲在幅值譜中的疊加。

(2) 對(duì)比不同工況下的頻率組成，堵塞、泄漏信號(hào)含有四個(gè)相同的頻率(38.09Hz和39.06Hz被視為同一頻率)，但是泄漏信號(hào)多出了77.15Hz的中高頻分量和154.3Hz、173.8Hz的高頻分量，這就導(dǎo)致泄漏狀態(tài)下的執(zhí)行裝置的提升聲音更加尖銳。而氣穴信號(hào)相比其他兩種工況出現(xiàn)了幅值較大的3.91的低頻分量，導(dǎo)致氣穴狀態(tài)的工作聲音比較低沉。

表1 不同狀態(tài)下原信號(hào)和分量中頻率及其幅值

3 結(jié) 論

本文利用奇異值分解中奇異值對(duì)不同頻率振幅大小的分辨能力，提出了一種基于奇異值分解的液壓信號(hào)的時(shí)域分解方法，方法根據(jù)原信號(hào)的頻率組成，依次向其中疊加頻率相同、振幅已知的周期信號(hào)，再選取特定的奇異值重構(gòu)，最終還原出每個(gè)頻率的時(shí)間序列。

分別通過(guò)仿真信號(hào)和液壓信號(hào)對(duì)該方法性能進(jìn)行驗(yàn)證，并與EMD的分解結(jié)果就模態(tài)混疊和噪聲含量做了比較。本文提出的方法不僅消除了分量中模態(tài)混疊，而且抑制了分量中的噪聲含量。

[ 1 ] 唐宏賓， 吳運(yùn)新， 滑廣軍，等. 基于EMD包絡(luò)譜分析的液壓泵故障診斷方法[J]. 振動(dòng)與沖擊， 2012， 31(9)： 44-48.TANG Hongbin， WU Yunxin， HUA Guangjun， et al. Fault diagnosis of pump using EMD and envelope spectrum analysis[J]. Journal of Vibration and Shock， 2012， 31(9)： 44-48.

[ 2 ] HUANG N E， SHEN Z， LONG S R， et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[C]// Proceedings of the Royal Society of London A： Mathematical， Physical and Engineering Sciences. London： The Royal Society， 1998：903-995.

[ 3 ] HUANG N E， WU Manli， QU Wendong， et al. Applications of Hilbert-Huang transform to non-stationary financial time series analysis[J]. Applied Stochastic Models in Business and Industry， 2003， 19(3)： 245-268.

[ 4 ] WU Zhaohua， HUANG N E. Ensemble empirical mode decomposition： a noise-assisted data analysis method[J]. Advances in Adaptive Data Analysis， 2009， 1(1)： 1-41.

[ 5 ] RYAN D， KAISER J F. The use of a masking signal to improve empirical mode decomposition[C]// International Conference on Acousic， Speech and Signal Processing 2005 IEEE. Cambridge： MIT Press， 2005： 485-488.

[ 6 ] 趙玲， 劉小峰， 秦樹(shù)人，等. 消除經(jīng)驗(yàn)?zāi)B(tài)分解中混疊現(xiàn)象的改進(jìn)掩膜信號(hào)法[J]. 振動(dòng)與沖擊， 2010， 29(9)： 13-17. ZHAO Ling， LIU Xiaofeng， QIN Shuren， et al. Use of masking signal to improve empirical modede composition[J]. Journal of Vibration and Shock， 2010， 29(9)： 13-17.

[ 7 ] HASSANPOUR H， ZEHTABIAN A. Time domain signal enhancement based on an optimized singular vector denoising algorithm[J]. Digital Signal Processing， 2012， 22(5)： 786-794.

[ 8 ] 劉鎏， 閆云聚， 李鵬博. 奇異譜分解在超聲速無(wú)人機(jī)聲振試驗(yàn)數(shù)據(jù)處理中的應(yīng)用[J]. 振動(dòng)與沖擊， 2015， 34(3)： 28-34. LIU Liu， YAN Yunju， LI Pengbo. Singular value spectral decomposition and its application in acoustic vibration test data processing of a supersonic aircraft[J]. Journal of Vibration and Shock， 2015， 34(3)： 28-34.

[ 9 ] 劉敏， 張英堂， 李志寧，等. 基于自適應(yīng)奇異值標(biāo)準(zhǔn)譜和 EMD 的柴油機(jī)故障診斷[J]. 車用發(fā)動(dòng)機(jī)， 2015(2)： 77-82. LIU Min， ZHANG Yingtang， LI Zhining， et al. Diesel engine fault diagnosis based on adaptive sigular value stantard spectrum and empirical mode decomposition[J]. Vehicle Engine， 2015(2)： 77-82.

[10] 王超， 孔凡讓， 黃偉國(guó)，等. 改進(jìn)的奇異值分解在軸承故障診斷中的應(yīng)用[J]. 振動(dòng)工程學(xué)報(bào)， 2014， 27(2)： 296-303. WANG Chao， KONG Fanrang， HUANG Weiguo， et al. Application of improved singular value decompostion in bearing fault diagnosis[J]. Journal of Vibration Engineering， 2014，27 (2)： 296-303.

[11] 胥永剛， 謝志聰， 孟志鵬，等. 基于奇異值分解的磁記憶信號(hào)特征提取方法[J]. 振動(dòng)、測(cè)試與診斷， 2014， 34(6)： 1105-1109. XU Yonggang， XIE Zhicong， MENG Zhipeng， et al. Feature extraction method of magnetic memory signal based on SVD[J]. Journal of Vibration， Measurement & Diagnosis， 2014， 34(6)： 1105-1109.

[12] 徐彥凱， 雙凱. 自適應(yīng)奇異值分解瞬變信號(hào)檢測(cè)研究[J]. 電子與信息學(xué)報(bào)， 2014， 36(3)： 583-588. XU Yankai， SHUANG Kai. Detection of transient signal based on adaptive singular value decomposition[J]. Journal of Electronics & Information Technology，2014， 36(3)： 583-588.

[13] 趙學(xué)智， 葉邦彥， 陳統(tǒng)堅(jiān). 奇異值差分譜理論及其在車床主軸箱故障診斷中的應(yīng)用[J]. 機(jī)械工程學(xué)報(bào)， 2010， 46(1)： 100-108. ZHAO Xuezhi， YE Bangyan， CHEN Tongjian. Difference spectrum theory of singular value and its application to the fault diagnosis of headstock of lathe[J]. Journal of Mechanical Engineering， 2010， 46(1)： 100-108.

[14] 錢征文， 程禮， 李應(yīng)紅. 利用奇異值分解的信號(hào)降噪方法[J]. 振動(dòng)、測(cè)試與診斷， 2011， 31(4)： 459-463. QIAN Zhengwen， CHENG Li， LI Yinghong. Noise reduction method based on singular value decomposition[J]. Journal of Vibration， Measurement & Diagnosis， 2011， 31(4)： 459-463.

A hydraulic signal decomposition method based on singular value decomposition

ZHANG Xiaoming， TANG Jian， ZHANG Meijun

(Engineering Institute of Engineering Corps， PLA University of Science and Technology， Nanjing 210007， China)

In order to suppress the mode aliasing phenomenon and reduce noise composition， a signal decomposition method in time domain based on singular value decomposition (SVD) was proposed. Based on two features of SVD， firstly each frequency corresponds to two sizeable singular values. Secondly singular values are positively related to the amplitude of its corresponded frequency. The method was conducted by adding a known simulation sine signal with appropriate amplitude to make the location of singular values easier to be identified. Then， time series were reconstructed by choosing related singular values. Finally， the time series of a certain frequency could be achieved by subtracting added simulation signal. By comparing with the EMD， it is effectively confirmed that the method can both eliminate mode aliasing and reduce noise composition.

singular value decomposition； construct signals； superposition of signals； mode aliasing； noise composition

國(guó)家自然科學(xué)基金 (51175511；61472392)；省青年基金(BK20150724)

2016-03-21 修改稿收到日期： 2016-06-28

張小明 男，碩士生，1991年12月生

唐建 女， 碩士生導(dǎo)師， 1977年1月生

TH137

A

10.13465/j.cnki.jvs.2017.16.015