• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于能量衰減比的雙通道源數(shù)目估計(jì)方法*

    2016-08-03 01:13:28周以齊
    振動(dòng)、測(cè)試與診斷 2016年2期

    于 剛, 周以齊

    (山東大學(xué)高效潔凈機(jī)械制造教育部重點(diǎn)實(shí)驗(yàn)室 濟(jì)南,250061)

    ?

    基于能量衰減比的雙通道源數(shù)目估計(jì)方法*

    于剛,周以齊

    (山東大學(xué)高效潔凈機(jī)械制造教育部重點(diǎn)實(shí)驗(yàn)室濟(jì)南,250061)

    摘要利用源信號(hào)到達(dá)兩傳感器間具有不同能量衰減比的特性,提出了一種新的源數(shù)估計(jì)方法,解決了因傳感器數(shù)量不足而無(wú)法準(zhǔn)確估計(jì)源信號(hào)數(shù)目的問(wèn)題。首先,利用線性時(shí)頻變換方法得到兩觀測(cè)信號(hào)在頻域的能量分布,然后,計(jì)算能量散點(diǎn)圖中對(duì)應(yīng)角度上的能量總和;最后,通過(guò)峰值檢測(cè)法實(shí)現(xiàn)源數(shù)目的自動(dòng)檢測(cè)。通過(guò)理論分析、仿真和實(shí)驗(yàn),證明了該方法的有效性。本方法為盲源分離算法處理振動(dòng)、噪聲信號(hào),提供了可靠的先驗(yàn)信息。

    關(guān)鍵詞源數(shù)目估計(jì); 不相關(guān)源; 相關(guān)源; 欠定

    引言

    盲源分離(blind source separation, 簡(jiǎn)稱BSS)技術(shù)僅通過(guò)測(cè)量信號(hào)即可估計(jì)出源信號(hào),為提取振動(dòng)、噪聲信號(hào)特征提供了新的思路[1-3]。BSS算法的成功應(yīng)用需要滿足傳感器數(shù)目(m)大于或者等于源信號(hào)數(shù)目(n)的假設(shè)。但實(shí)際測(cè)量中,信號(hào)源個(gè)數(shù)往往是無(wú)法預(yù)知的,這一點(diǎn)限制了BSS方法在工程中的應(yīng)用。當(dāng)m

    傳統(tǒng)源數(shù)目估計(jì)有主分量分析和奇異值分解法,通過(guò)識(shí)別測(cè)量信號(hào)協(xié)方差矩陣的非零奇異值或者非零特征值來(lái)估計(jì)不相關(guān)源的個(gè)數(shù)[4]。文獻(xiàn)[5]結(jié)合奇異值分解和聚類方法解決了不相關(guān)源的數(shù)目估計(jì)問(wèn)題。文獻(xiàn)[6]將觀測(cè)信號(hào)變換到頻域中,利用非負(fù)矩陣分解的方法實(shí)現(xiàn)源數(shù)目估計(jì)。文獻(xiàn)[4-6]所用方法均需要滿足傳感器數(shù)目大于源信號(hào)。

    多個(gè)文獻(xiàn)針對(duì)欠定情況下源數(shù)估計(jì)問(wèn)題進(jìn)行了研究。文獻(xiàn)[7]通過(guò)觀測(cè)信號(hào)間的功率譜比值估計(jì)了不相關(guān)源的數(shù)目。文獻(xiàn)[8]將觀測(cè)信號(hào)從時(shí)域變換到小波域,利用重構(gòu)的小波系數(shù)信號(hào)進(jìn)行奇異值分解,得到源個(gè)數(shù)估計(jì)。文獻(xiàn)[9]利用經(jīng)驗(yàn)?zāi)B(tài)分解得到觀測(cè)信號(hào)的本征模函數(shù),同時(shí)利用奇異值分解得到多個(gè)虛擬觀測(cè)通道的特征分布,再利用貝葉斯信息準(zhǔn)則判斷源信號(hào)數(shù)目。文獻(xiàn)[10]針對(duì)機(jī)械運(yùn)行過(guò)程中的動(dòng)態(tài)故障源,采用拓展四階累積量矩陣自適應(yīng)地估計(jì)超定、正定及欠定情況下的源數(shù)目。文獻(xiàn)[8-10]通過(guò)將觀測(cè)信號(hào)進(jìn)行分解,把欠定問(wèn)題轉(zhuǎn)化為正定問(wèn)題后,再估計(jì)源信號(hào)數(shù)目。由于需要用到信號(hào)分解方法,增加了計(jì)算復(fù)雜度,同時(shí)易受到干擾噪聲信號(hào)的影響。

    本研究利用信號(hào)源在傳遞過(guò)程中到達(dá)兩傳感器具有不同能量衰減比的特性,依靠?jī)赏ǖ赖臏y(cè)量信號(hào)可估計(jì)任意數(shù)量的不相關(guān)信號(hào)源,同時(shí)對(duì)信號(hào)間相關(guān)成分不大的相關(guān)源可實(shí)現(xiàn)正確的數(shù)目估計(jì)。

    1算法原理

    信號(hào)源以能量的形式向外傳播,源信號(hào)到達(dá)傳感器時(shí)具有一定的衰減。當(dāng)機(jī)械設(shè)備穩(wěn)定運(yùn)行時(shí),單位時(shí)間內(nèi)不同傳感器間所獲得同一源的能量比值是固定的。僅存在單個(gè)源時(shí),其能量計(jì)算可用公式

    (1)

    其中:Ei(i=1,2,…,n)為第i個(gè)傳感器接收到的信號(hào)能量;xi為第i個(gè)傳感器接收到的時(shí)域信號(hào)。

    但機(jī)械設(shè)備運(yùn)行中存在的源信號(hào)有多個(gè),既有平穩(wěn)信號(hào)又有非平穩(wěn)信號(hào),同時(shí)還混有噪音信號(hào)。無(wú)法僅利用實(shí)際測(cè)量信號(hào)在時(shí)域中的公式計(jì)算多個(gè)源的能量。

    線性時(shí)頻變換是一種可以將時(shí)域信號(hào)展開到時(shí)頻域中的方法,較為常用的有短時(shí)傅里葉變換、連續(xù)小波變換等。定義觀測(cè)信號(hào)在時(shí)頻域中的頻率能量分布公式為

    (2)

    其中:Ei(f)為第i個(gè)傳感器接收到的信號(hào)在頻率上的能量分布;Re(xi(t,f))與Im(xi(t,f))為第i個(gè)傳感器信號(hào)在時(shí)頻變換后的實(shí)部與虛部。

    如果源信號(hào)之間不相關(guān)或者僅部分相關(guān)(即沒(méi)有頻率共存或者僅部分頻率共存),當(dāng)存在某一個(gè)頻率點(diǎn)fk僅屬于一個(gè)源時(shí),E1(fk)/E2(fk)即為此源信號(hào)到達(dá)兩傳感器之間的能量衰減比。系統(tǒng)穩(wěn)定運(yùn)行時(shí),對(duì)所有屬于同一源的頻率點(diǎn),兩通道間Ei(f)的比值是固定的。證明如下。

    假設(shè)有兩個(gè)源混入到兩觀測(cè)信號(hào)

    (3)

    其中:a11,a21,a12,a22分別為源s1和s2到達(dá)傳感器x1和x2的衰減;Δt1,Δt2為源信號(hào)之間的延遲。

    對(duì)式(3)做時(shí)頻變換之后的實(shí)、虛部分別為

    (4)

    當(dāng)某個(gè)頻率成分fk僅存在于s1時(shí),則有

    (5)

    (6)

    在頻點(diǎn)fk的兩通道能量比為

    (7)

    當(dāng)某個(gè)頻率成分fh僅存在于s2時(shí),則有

    (8)

    式(7)與式(8)表明,當(dāng)把所有E1(f)與E2(f)的值畫在散點(diǎn)圖上時(shí),若a11/a21≠a12/a22,則會(huì)形成兩條方向不同的聚類直線,其角度代表不同的源信號(hào)在兩傳感器間的能量衰減比,聚類直線的個(gè)數(shù)即為信號(hào)源的個(gè)數(shù)。以仿真信號(hào)為例說(shuō)明算法過(guò)程。

    2仿真驗(yàn)證

    構(gòu)造5個(gè)仿真信號(hào),S1,S3為非平穩(wěn)調(diào)制信號(hào),模擬旋轉(zhuǎn)機(jī)械的故障信號(hào);S2,S4為周期平穩(wěn)信號(hào),模擬設(shè)備正常運(yùn)行信號(hào);S5為振蕩衰減信號(hào),模擬振動(dòng)沖擊信號(hào),信號(hào)長(zhǎng)度為1 024。為了說(shuō)明算法對(duì)相關(guān)源的估計(jì)能力,利用系數(shù)c12與c34來(lái)調(diào)整S1,S2間和S3,S4間相關(guān)成分的大小。各源信號(hào)的時(shí)域波形如圖1,設(shè)定c12=0,c34=0。源信號(hào)具體參數(shù)為

    S1=sin(100πt)(1+cos(20πt))

    S2=sin(200πt)+c12sin(120πt)

    S3=sin(300πt)(1+cos(20πt))

    S4=sin(400πt)+c34sin(320πt)

    設(shè)混合矩陣為

    圖1 源信號(hào)S1~S5時(shí)域波形Fig.1 The sources S1~S5 in time-domain

    構(gòu)造觀測(cè)信號(hào),同時(shí)混入20%的噪聲信號(hào)

    x1(t)=a11S1(t)+a12S2(t+Δt2)+

    a13S3(t)+a14S4(t+Δt4)+a15S5(t)

    x2(t)=a21S1(t+Δt1)+a22S2(t)+

    a23S3(t+Δt3)+a24S4(t)+a25S5(t+Δt5)其中:Δti代表不同源到達(dá)兩傳感器間的延遲。

    2.1仿真信號(hào)源數(shù)目估計(jì)

    以不相關(guān)源為例,設(shè)定c12=0,c34=0。選擇連續(xù)小波變換對(duì)信號(hào)x1(t)、x2(t)進(jìn)行時(shí)頻變換,選取cmor3-3為小波基函數(shù),尺度為102 4。根據(jù)式(2),求得小波變換后兩觀測(cè)信號(hào)在頻域中的能量分布。在兩通道頻域能量值所形成的散點(diǎn)圖中,所有屬于同一個(gè)源的頻點(diǎn)能量值會(huì)聚類成不同角度的直線,圖2中共有5條明顯的聚類直線。為了實(shí)現(xiàn)自動(dòng)檢測(cè),將坐標(biāo)軸從橫軸到縱軸按角度分為90份,計(jì)算每個(gè)角度上散點(diǎn)值的總和,得到90個(gè)能量和值,見圖3。然后利用峰值檢測(cè)方法得到每條聚類直線對(duì)應(yīng)的角度和峰值的個(gè)數(shù)[11],如圖中“*”標(biāo)示,峰值個(gè)數(shù)即為源信號(hào)數(shù)目。通過(guò)設(shè)定一定的閥值(設(shè)為最大值的0.02),可以消除幅值較小的噪音峰值影響,以此來(lái)實(shí)現(xiàn)自動(dòng)判定源個(gè)數(shù)的目的。

    圖2 能量散點(diǎn)圖Fig.2 The energy scatter plot

    圖3 能量峰值圖Fig.3 The energy peak plot

    算法具體流程為:

    1) 對(duì)兩傳感器信號(hào)x1(t)和x2(t)進(jìn)行連續(xù)小波變換,得到兩組小波變換系數(shù)X1(t,f)和X2(t,f);

    2) 根據(jù)式(2),計(jì)算每個(gè)頻率下的能量值,即E1(f)與E2(f);

    3) 利用E1(f)與E2(f)得到兩通道信號(hào)的能量散點(diǎn)圖;

    4) 計(jì)算每個(gè)角度上對(duì)應(yīng)的能量和,將頻域能量值從散點(diǎn)圖變換到峰值圖中;

    5) 設(shè)置一定的閥值,利用峰值檢測(cè)法得到源數(shù)目的估計(jì)。

    2.2相關(guān)成分對(duì)源數(shù)估計(jì)的影響

    本節(jié)討論相關(guān)成分大小對(duì)源數(shù)目估計(jì)的影響。當(dāng)c12=0.1,c34=0.1;c12=0.2,c34=0.2;c12=0.4,c34=0.4時(shí),所對(duì)應(yīng)的能量峰值圖如圖4所示。從圖中可以看出,當(dāng)源信號(hào)之間的相關(guān)成分較小時(shí),可以通過(guò)設(shè)定一定的閥值,忽略掉相關(guān)成分形成的峰值,從而實(shí)現(xiàn)源數(shù)的正確估計(jì)。隨著相關(guān)成分比重的增加,相關(guān)成分所形成的虛假峰值會(huì)淹沒(méi)源信號(hào)的峰值,導(dǎo)致出現(xiàn)誤估計(jì)。

    圖4 c12,c34取三組不同數(shù)值時(shí)相關(guān)源的能量峰值圖Fig.4 The energy peak plot of correlated source signals corresponding to three different values of c12,c34

    2.3衰減比的接近程度對(duì)源數(shù)估計(jì)的影響

    從仿真信號(hào)來(lái)看,提出的算法可以正確地估計(jì)出任意數(shù)量的不相關(guān)源數(shù)目,但前提是各源信號(hào)之間具有不同的能量衰減比。實(shí)際測(cè)試中,如果存在兩源信號(hào)到達(dá)兩傳感器具有相近或者相同的衰減比時(shí),則兩源信號(hào)在能量散點(diǎn)圖中形成的聚類直線會(huì)較為接近甚至重合,會(huì)誤將兩個(gè)源信號(hào)估計(jì)為一個(gè)源信號(hào)。本節(jié)討論源信號(hào)傳遞到兩傳感器衰減比的接近程度對(duì)源數(shù)估計(jì)的影響。

    假設(shè)兩源信號(hào)在散點(diǎn)圖中所形成聚類直線的方向向量分別為A1,A2,θ為A1,A2的角距離。將坐標(biāo)軸等分為90份,即從1°到90°,然后計(jì)算每個(gè)劃分的角度空間所包含能量值的總和,得到90個(gè)能量和值,再利用峰值檢測(cè)法得到峰值的個(gè)數(shù),即為源信號(hào)數(shù)目。A1,A2所形成峰值的個(gè)數(shù)與θ的關(guān)系如圖5所示。

    圖5 A1,A2夾角與峰值個(gè)數(shù)的關(guān)系Fig.5 Relationship between the angle of A1,A2 and the number of peak

    圖5所示的幾種情況中,(a),(b)會(huì)形成1個(gè)峰值,(c),(d)會(huì)形成2個(gè)峰值。所以,當(dāng)0°≤θ≤1°時(shí),會(huì)將兩個(gè)源信號(hào)誤估計(jì)為一個(gè)源信號(hào);1°<θ≤2°時(shí),形成峰值的個(gè)數(shù)與A1,A2在坐標(biāo)軸中的相對(duì)位置有關(guān),具有一定的隨機(jī)性;2°<θ≤90°時(shí),會(huì)形成2個(gè)峰值,可以正確地估計(jì)出源信號(hào)數(shù)目,定義(2°,90°]為有效估計(jì)區(qū)間。有效估計(jì)區(qū)間的最小值為坐標(biāo)軸單個(gè)劃分區(qū)間的2倍。因此,增加坐標(biāo)軸的劃分份數(shù),將增大有效估計(jì)區(qū)間的范圍,有利于提高源數(shù)估計(jì)的精度。同時(shí),實(shí)際測(cè)試中變換位置多測(cè)幾組數(shù)據(jù),使得不同源信號(hào)到達(dá)傳感器具有不同的衰減比,一定程度上也會(huì)減少誤估計(jì)情況的出現(xiàn)。

    3實(shí)驗(yàn)驗(yàn)證

    為進(jìn)一步驗(yàn)證算法的有效性,將分別對(duì)實(shí)測(cè)的振動(dòng)、噪聲信號(hào)進(jìn)行源估計(jì)。每個(gè)實(shí)例用4張圖來(lái)說(shuō)明算法過(guò)程。前2個(gè)為兩測(cè)量信號(hào)經(jīng)小波變換后的時(shí)頻分布圖,第3個(gè)為兩通道信號(hào)的能量散點(diǎn)圖,第4個(gè)為能量峰值圖。

    3.1實(shí)例1

    測(cè)量信號(hào)來(lái)自于某泵水系統(tǒng)中齒輪箱的振動(dòng)加速度數(shù)據(jù)。已知齒輪箱傳動(dòng)為兩級(jí)傳動(dòng),主動(dòng)輪輪齒齒數(shù)為13,輸入轉(zhuǎn)速為997 r/min,采樣頻率為3 200 Hz。因此,通過(guò)計(jì)算得出此齒輪箱系統(tǒng)的1階嚙合頻率為13×997/60=216 Hz。

    從測(cè)量信號(hào)1和2的時(shí)頻圖(圖6(a)和(b))可以看出,有4個(gè)主要頻率成分,其頻率分別約為216,281,432和864 Hz。則216,432和864 Hz頻率成分應(yīng)為齒輪嚙合的1階頻率及其2階與4階倍頻。281 Hz成分應(yīng)為齒輪安裝存在的輕微偏心所造成的1階邊頻成分。從能量散點(diǎn)圖(圖6(c))中可以看到,有4條聚類直線。能量峰值圖(圖6(d))中,有4個(gè)峰值被檢測(cè)出來(lái)。

    圖6 齒輪箱振動(dòng)加速度信號(hào)Fig.6 Gearbox vibration acceleration signals

    3.2實(shí)例2

    測(cè)量信號(hào)為聲壓信號(hào),采集于某挖掘機(jī)消聲器排氣口附近。此挖掘機(jī)采用4缸4沖程內(nèi)燃機(jī),當(dāng)前轉(zhuǎn)速為2 400 r/min。消聲器排氣噪聲來(lái)自于內(nèi)燃機(jī),主要包括氣門開閉時(shí)拍擊引起的機(jī)械噪聲與燃燒爆發(fā)引起的燃燒噪聲。由內(nèi)燃機(jī)相關(guān)參數(shù)可計(jì)算得出:發(fā)火頻率為(2400×2)/60=80 Hz,氣門開合拍擊的主要頻率為(2 400×4)/30=320 Hz。

    數(shù)據(jù)采樣頻率為5 120 Hz,由于其主要頻率成分集中在0~500 Hz之間,為了更清晰地展示測(cè)量信號(hào)的時(shí)頻特征,截取信號(hào)這一部分的時(shí)頻分布圖。從觀測(cè)信號(hào)1和2的時(shí)頻圖(圖7(a)和(b))可以看到,有3個(gè)主要頻率成分,分別約為80,160 Hz和230~370 Hz的寬頻成分。80與160 Hz成分應(yīng)為發(fā)火頻率的1階頻率與2階頻率。230~370 Hz的寬頻成分應(yīng)為氣門開合拍擊噪聲在消聲器中產(chǎn)生的氣流再生噪聲,故其頻率呈現(xiàn)寬頻特性。從散點(diǎn)圖(圖7(c))可以看到,有3條聚類直線。在峰值圖中(圖7(d)),有3個(gè)明顯的峰值被估計(jì)出來(lái)。

    圖7 聲壓信號(hào)Fig.7 Sound pressure signals

    3.3實(shí)例3

    數(shù)據(jù)采集于某懸臂梁錘擊實(shí)驗(yàn)的兩個(gè)測(cè)點(diǎn)通道,從兩通道信號(hào)的時(shí)頻分布圖(圖8(a)和(b))中可以看出主要有3階模態(tài)被激發(fā)出來(lái),其頻率分別約為100,200和600 Hz。如果每一階模態(tài)響應(yīng)可當(dāng)作一個(gè)源信號(hào),則兩通道信號(hào)中存在有3個(gè)源。從能量散點(diǎn)圖(圖8(c))中可以看出,有3條明顯的聚類直線。能量峰值圖(圖8(d))中有3個(gè)峰值被檢測(cè)出來(lái)。

    圖8 模態(tài)信號(hào)Fig.8 Modal signals

    4結(jié)束語(yǔ)

    本研究利用源信號(hào)傳遞到兩個(gè)傳感器具有不同的能量衰減比來(lái)估計(jì)信號(hào)源數(shù)目,解決了因傳感器數(shù)量不足而無(wú)法準(zhǔn)確估計(jì)源信號(hào)數(shù)目的問(wèn)題。僅利用兩傳感器信號(hào)可估計(jì)出任意數(shù)量的信號(hào)源,但同時(shí)也存在局限性:a.針對(duì)相關(guān)源信號(hào)的數(shù)目估計(jì),當(dāng)信號(hào)間相關(guān)成分較小時(shí),可通過(guò)設(shè)置一定的閥值濾除虛假峰值的影響;當(dāng)相關(guān)成分較大時(shí),形成的虛假峰值會(huì)淹沒(méi)真實(shí)源信號(hào)的峰值,易造成誤估計(jì)。b.當(dāng)兩個(gè)源信號(hào)到達(dá)兩傳感器具有相近或者相同的衰減比時(shí),所提算法會(huì)將兩個(gè)源信號(hào)誤估計(jì)為一個(gè)源信號(hào)。針對(duì)后者局限性,可以通過(guò)增加坐標(biāo)軸的劃分份數(shù)來(lái)提高源數(shù)估計(jì)精度;同時(shí),通過(guò)調(diào)整測(cè)試位置,使得不同源信號(hào)到達(dá)傳感器具有不同的衰減比,以減少誤估計(jì)情況的出現(xiàn)。

    實(shí)驗(yàn)階段,通過(guò)機(jī)械工程領(lǐng)域中3種常用實(shí)測(cè)信號(hào)進(jìn)行驗(yàn)證,結(jié)果表明,該方法能夠準(zhǔn)確估計(jì)出源數(shù)目,可以為盲源分離方法提供可靠的先驗(yàn)信息,同樣也適用于其他工程實(shí)踐應(yīng)用對(duì)信號(hào)源個(gè)的數(shù)估計(jì)。

    參考文獻(xiàn)

    [1]李舜酩. 轉(zhuǎn)子振動(dòng)故障信號(hào)的盲分離[J]. 航空動(dòng)力學(xué)報(bào), 2006, 20(5): 751-756.

    Li Shunming. Blind source separation of rotor vibration faults[J]. Journal of Aerospace Power, 2006, 20(5): 751-756. (in Chinese)

    [2]艾延廷, 費(fèi)成巍, 張鳳玲, 等. ICA在航空發(fā)動(dòng)機(jī)振動(dòng)信號(hào)盲源分離中的應(yīng)用[J]. 振動(dòng)、測(cè)試與診斷, 2010, 30(6): 671-674.

    Ai Yanting, Fei Chengwei, Zhang Fengling, et al. Blind source separation for aero-engines vibration signal by independent component analysis[J]. Journal of Vibration, Measurement & Diagnosis, 2010, 30(6): 671-674. (in Chinese)

    [3]Kerschen G, Poncelet F, Golinval J C. Physical interpretation of independent component analysis in structural dynamics[J]. Mechanical Systems and Signal Processing, 2007, 21(4): 1561-1575.

    [4]Kompella M S, Davies P, Bernhard R J, et al. A technique to determine the number of incoherent sources contributing to the response of a system[J]. Mechanical Systems and Signal Processing, 1994, 8(4): 363-380.

    [5]焦衛(wèi)東, 楊世錫, 吳昭同. 基于源數(shù)估計(jì)的旋轉(zhuǎn)機(jī)械源盲分離[J]. 中國(guó)機(jī)械工程, 2003, 14(14): 1184-1187.

    Jiao Weidong, Yang Shixi, Wu Zhaoton. A method of blind source separation for rotating machinery based on estimation of the number of sources[J]. China Mechanical Engineering, 2003, 14(14): 1184-1187. (in Chinese)

    [6]李寧, 史鐵林. 基于非負(fù)矩陣分解的盲信號(hào)源數(shù)估計(jì)[J]. 中國(guó)機(jī)械工程, 2007, 18(22): 2734-2737.

    Li Ning, Shi Tielin. Estimation of the number of blind sources based on non-negative matrix factorization[J]. China Mechanical Engineering, 2007, 18(22): 2734-2737. (in Chinese)

    [7]李寧, 史鐵林. 一種新的盲信號(hào)源數(shù)估計(jì)方法研究[J]. 中國(guó)機(jī)械工程, 2007, 18(19): 2298-2302.

    Li Ning, Shi Tielin. A new method for estimating the number of blind sources[J]. China Mechanical Engineering, 2007, 18(19): 2298-2302. (in Chinese)

    [8]毋文峰, 陳小虎, 蘇勛家, 等. 機(jī)械振動(dòng)源數(shù)估計(jì)的小波方法[J]. 機(jī)械科學(xué)與技術(shù), 2011, 30(10): 1679-1682.

    Wu Wenfeng, Chen Xiaohu, Su Xunjia, et al. Wavelet decomposition algorithm for source number estimation of mechanical vibration[J]. Mechanical Science and Technology for Aerospace Engineering. 2011, 30(10): 1679-1682. (in Chinese)

    [9]葉紅仙, 楊世錫, 楊將新. 基于EMD-SVD-BIC的機(jī)械振動(dòng)源數(shù)估計(jì)方法[J]. 振動(dòng)、測(cè)試與診斷, 2010, 30(3): 330-334.

    Ye Hongxian, Yang Shixi, Yang Jianxin. Mechanical vibration source number estimation based on EMD-SVD-BIC[J]. Journal of Vibration, Measurement & Diagnosis, 2010, 30(3): 330-334. (in Chinese)

    [10]陽(yáng)小燕, 唐一科, 周雄. 動(dòng)態(tài)故障源數(shù)估計(jì)的自適應(yīng)盲源分離方法[J]. 振動(dòng)、測(cè)試與診斷, 2013,33 (2): 290-293.

    Yang Xiaoyan, Tan Yike, Zhou Xiong. Method of adaptive blind source separation based on dynamic fault source number estimation[J]. Journal of Vibration, Measurement & Diagnosis, 2013,33 (2): 290-293. (in Chinese)

    [11]Billauer E. peakdet: Peak detection using MATLAB.[EB/OL].(2008-08-01)[2012-07-20]. http://www.billauer.co.il/peakdet.html.

    E-mail:yugang2010@163.com

    doi:10.16450/j.cnki.issn.1004-6801.2016.02.016

    收稿日期:2014-06-12;修回日期:2014-08-27

    中圖分類號(hào)TH17

    第一作者簡(jiǎn)介:于剛,男,1987年12月生,博士生。主要研究方向?yàn)樾盘?hào)處理。曾發(fā)表《General linear chirplet transform》(《Mechanical Systems and Signal Processing》2016,Vol.70,No.1)等論文。

    *國(guó)家科技支撐計(jì)劃資助項(xiàng)目(2015BAF07B04)

    亚洲婷婷狠狠爱综合网| 最后的刺客免费高清国语| 国产美女午夜福利| 中文在线观看免费www的网站| 日日摸夜夜添夜夜添av毛片| 身体一侧抽搐| 国产精品久久久久久精品电影| 最近视频中文字幕2019在线8| 毛片女人毛片| 国产综合懂色| av.在线天堂| 少妇人妻一区二区三区视频| 偷拍熟女少妇极品色| 国产淫片久久久久久久久| 国产乱人偷精品视频| 搞女人的毛片| 国产一区二区三区av在线| 亚洲18禁久久av| 亚洲婷婷狠狠爱综合网| 中文在线观看免费www的网站| 一个人看的www免费观看视频| 欧美日韩在线观看h| 久久精品综合一区二区三区| 91在线精品国自产拍蜜月| 国产av国产精品国产| 18禁裸乳无遮挡免费网站照片| 波野结衣二区三区在线| 在线观看美女被高潮喷水网站| 亚洲熟妇中文字幕五十中出| 国产亚洲最大av| 精品少妇黑人巨大在线播放| 看十八女毛片水多多多| 伊人久久精品亚洲午夜| 你懂的网址亚洲精品在线观看| 亚洲精品亚洲一区二区| 亚洲成色77777| 精品国产露脸久久av麻豆 | 边亲边吃奶的免费视频| 午夜福利视频1000在线观看| 欧美不卡视频在线免费观看| 天堂av国产一区二区熟女人妻| 日本黄色片子视频| 91精品国产九色| 亚洲精品一区蜜桃| 青春草视频在线免费观看| 白带黄色成豆腐渣| 搞女人的毛片| 夜夜看夜夜爽夜夜摸| 国产免费福利视频在线观看| 午夜福利高清视频| 亚洲成人中文字幕在线播放| 日韩精品青青久久久久久| 欧美bdsm另类| 黄片wwwwww| xxx大片免费视频| 亚洲精品乱久久久久久| 国产av不卡久久| 最近视频中文字幕2019在线8| 在线观看美女被高潮喷水网站| 欧美激情在线99| 成人午夜高清在线视频| 在线免费观看的www视频| 国产成人精品一,二区| 在线免费观看不下载黄p国产| 真实男女啪啪啪动态图| 欧美三级亚洲精品| 插逼视频在线观看| 三级男女做爰猛烈吃奶摸视频| 夫妻性生交免费视频一级片| 成年女人看的毛片在线观看| 成人毛片a级毛片在线播放| 男女边摸边吃奶| 亚洲熟女精品中文字幕| 汤姆久久久久久久影院中文字幕 | 亚洲国产精品专区欧美| 亚洲一级一片aⅴ在线观看| 国产一区亚洲一区在线观看| 国产亚洲最大av| 一夜夜www| 亚洲高清免费不卡视频| 日日干狠狠操夜夜爽| 69av精品久久久久久| 午夜福利在线观看吧| 国语对白做爰xxxⅹ性视频网站| 免费观看的影片在线观看| 国产成人aa在线观看| 啦啦啦啦在线视频资源| 久热久热在线精品观看| 亚洲精品一区蜜桃| 亚洲三级黄色毛片| 国产伦理片在线播放av一区| 大香蕉久久网| 成人毛片60女人毛片免费| av免费在线看不卡| 国产精品久久视频播放| 一级av片app| 色综合站精品国产| 久久久久久伊人网av| 精品不卡国产一区二区三区| 99久久人妻综合| 国产黄色视频一区二区在线观看| 一区二区三区四区激情视频| 久久午夜福利片| 激情 狠狠 欧美| 日韩欧美三级三区| 国产成人精品婷婷| av国产免费在线观看| 1000部很黄的大片| 日韩人妻高清精品专区| 日韩中字成人| 亚洲欧美精品专区久久| 国产成人免费观看mmmm| 一级毛片黄色毛片免费观看视频| 免费黄网站久久成人精品| 日韩成人伦理影院| 日日摸夜夜添夜夜添av毛片| 欧美最新免费一区二区三区| 日韩人妻高清精品专区| 免费电影在线观看免费观看| 精品国产三级普通话版| 美女黄网站色视频| av在线天堂中文字幕| 午夜老司机福利剧场| 色视频www国产| 国产高清三级在线| 高清视频免费观看一区二区 | 国产成人a∨麻豆精品| 最近的中文字幕免费完整| 久久精品国产鲁丝片午夜精品| 精品国内亚洲2022精品成人| 九草在线视频观看| 久久精品人妻少妇| 久久这里只有精品中国| 欧美性感艳星| 成人av在线播放网站| 精品99又大又爽又粗少妇毛片| 51国产日韩欧美| 美女主播在线视频| 国产精品人妻久久久影院| 好男人视频免费观看在线| 亚洲av免费高清在线观看| 国产精品久久久久久久久免| 国产伦精品一区二区三区四那| av福利片在线观看| 国产精品无大码| 亚洲av福利一区| 在线观看免费高清a一片| 少妇熟女欧美另类| 国产大屁股一区二区在线视频| 国产探花极品一区二区| av网站免费在线观看视频 | 国产高清国产精品国产三级 | 国产极品天堂在线| 免费不卡的大黄色大毛片视频在线观看 | 国产中年淑女户外野战色| 人人妻人人看人人澡| 亚洲乱码一区二区免费版| 丝袜美腿在线中文| 久久久精品免费免费高清| 日本与韩国留学比较| 色播亚洲综合网| 在线观看人妻少妇| 国产三级在线视频| 久久国内精品自在自线图片| 久久久久久久久久成人| 国产av不卡久久| a级毛色黄片| 亚洲成色77777| 一级片'在线观看视频| 日韩一区二区三区影片| 欧美xxⅹ黑人| 色综合站精品国产| 午夜久久久久精精品| 中文字幕av成人在线电影| 少妇裸体淫交视频免费看高清| 国产91av在线免费观看| 男的添女的下面高潮视频| 免费大片18禁| 十八禁国产超污无遮挡网站| 免费看a级黄色片| 亚洲aⅴ乱码一区二区在线播放| 免费黄频网站在线观看国产| 人人妻人人澡欧美一区二区| 亚洲国产精品成人综合色| 真实男女啪啪啪动态图| 亚洲成人中文字幕在线播放| 麻豆成人午夜福利视频| 三级经典国产精品| 国产精品精品国产色婷婷| 中文字幕免费在线视频6| 日本-黄色视频高清免费观看| 综合色av麻豆| 国产色爽女视频免费观看| 久久精品人妻少妇| av卡一久久| 日韩电影二区| 日韩伦理黄色片| 国内精品宾馆在线| 熟妇人妻久久中文字幕3abv| 国产精品一区二区三区四区久久| 男人狂女人下面高潮的视频| 国产精品一区二区性色av| 青春草亚洲视频在线观看| 亚洲精品乱码久久久v下载方式| 国产精品一区二区性色av| 日本av手机在线免费观看| 国产一区有黄有色的免费视频 | 国产亚洲精品久久久com| 国产伦理片在线播放av一区| 亚洲18禁久久av| 又爽又黄无遮挡网站| 精品熟女少妇av免费看| 午夜亚洲福利在线播放| 99热这里只有是精品50| av国产免费在线观看| 在线 av 中文字幕| 国产精品.久久久| 国产亚洲午夜精品一区二区久久 | 天天一区二区日本电影三级| 国产成年人精品一区二区| 国产精品一二三区在线看| av天堂中文字幕网| 精品一区在线观看国产| 黄色欧美视频在线观看| 又大又黄又爽视频免费| 韩国av在线不卡| 夫妻性生交免费视频一级片| 一个人免费在线观看电影| 又黄又爽又刺激的免费视频.| 汤姆久久久久久久影院中文字幕 | 亚洲最大成人中文| 亚洲自偷自拍三级| 内地一区二区视频在线| 天天躁夜夜躁狠狠久久av| 18禁在线无遮挡免费观看视频| 国产伦精品一区二区三区四那| 亚洲美女视频黄频| 99久国产av精品| 国内精品美女久久久久久| 久久97久久精品| 联通29元200g的流量卡| 国产精品久久久久久精品电影| 三级国产精品片| 在线天堂最新版资源| 国产黄色小视频在线观看| 国产日韩欧美在线精品| 精品久久久久久久久亚洲| 久久国内精品自在自线图片| 成人亚洲精品一区在线观看 | 91精品一卡2卡3卡4卡| 亚洲国产高清在线一区二区三| 搞女人的毛片| 亚洲熟女精品中文字幕| 久久精品久久精品一区二区三区| 免费黄网站久久成人精品| 国产精品一区二区三区四区免费观看| 男女下面进入的视频免费午夜| 天天躁日日操中文字幕| 色尼玛亚洲综合影院| 十八禁网站网址无遮挡 | 中文字幕av成人在线电影| 777米奇影视久久| 国产成人aa在线观看| 亚洲无线观看免费| 黄片无遮挡物在线观看| 91在线精品国自产拍蜜月| 欧美高清性xxxxhd video| 插逼视频在线观看| 亚洲av福利一区| 51国产日韩欧美| 搞女人的毛片| av免费观看日本| 久久韩国三级中文字幕| 能在线免费观看的黄片| 欧美性猛交╳xxx乱大交人| 在线免费观看不下载黄p国产| 国产色爽女视频免费观看| 国产午夜福利久久久久久| av又黄又爽大尺度在线免费看| 最近手机中文字幕大全| 亚洲天堂国产精品一区在线| .国产精品久久| 亚洲怡红院男人天堂| av卡一久久| 在线观看人妻少妇| 欧美三级亚洲精品| 亚洲自偷自拍三级| 日韩亚洲欧美综合| 亚洲av成人精品一区久久| 一级毛片黄色毛片免费观看视频| 男女下面进入的视频免费午夜| 在线观看免费高清a一片| 免费黄频网站在线观看国产| 国产精品一区二区性色av| 精品一区在线观看国产| 一级av片app| 免费大片黄手机在线观看| av免费观看日本| 插阴视频在线观看视频| 国产激情偷乱视频一区二区| 日韩一本色道免费dvd| 亚洲最大成人av| 久久久久性生活片| 欧美+日韩+精品| 3wmmmm亚洲av在线观看| 国产高清有码在线观看视频| 亚洲精品国产av蜜桃| 久久久久国产网址| 国产一区二区在线观看日韩| 卡戴珊不雅视频在线播放| 亚洲精品国产成人久久av| 久久久久性生活片| 欧美人与善性xxx| 最新中文字幕久久久久| 性色avwww在线观看| 波多野结衣巨乳人妻| 日韩欧美精品v在线| 亚洲高清免费不卡视频| 中文字幕人妻熟人妻熟丝袜美| 成人毛片60女人毛片免费| 国产精品爽爽va在线观看网站| av女优亚洲男人天堂| 国产在线男女| 日本与韩国留学比较| 两个人视频免费观看高清| 亚洲图色成人| av.在线天堂| 噜噜噜噜噜久久久久久91| 亚洲精品国产av蜜桃| 寂寞人妻少妇视频99o| av在线天堂中文字幕| 别揉我奶头 嗯啊视频| 精品一区二区三区人妻视频| 丰满乱子伦码专区| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 激情五月婷婷亚洲| 别揉我奶头 嗯啊视频| 成人毛片60女人毛片免费| 久久精品久久精品一区二区三区| 久久久a久久爽久久v久久| 亚洲伊人久久精品综合| 99久久人妻综合| 日韩精品有码人妻一区| 狠狠精品人妻久久久久久综合| 自拍偷自拍亚洲精品老妇| 非洲黑人性xxxx精品又粗又长| 午夜福利在线观看免费完整高清在| 大片免费播放器 马上看| 岛国毛片在线播放| 免费观看无遮挡的男女| 大香蕉久久网| 亚洲婷婷狠狠爱综合网| 亚洲成色77777| 男女视频在线观看网站免费| 夫妻性生交免费视频一级片| 国语对白做爰xxxⅹ性视频网站| 三级男女做爰猛烈吃奶摸视频| 欧美成人a在线观看| 日韩,欧美,国产一区二区三区| 美女xxoo啪啪120秒动态图| 91精品一卡2卡3卡4卡| 国产视频首页在线观看| 丝袜喷水一区| 国产国拍精品亚洲av在线观看| 精品人妻一区二区三区麻豆| 97超视频在线观看视频| 青春草视频在线免费观看| 亚洲国产欧美在线一区| 青春草视频在线免费观看| 国产男人的电影天堂91| 青春草视频在线免费观看| 一级片'在线观看视频| 一级爰片在线观看| 有码 亚洲区| 国产激情偷乱视频一区二区| 国产高清国产精品国产三级 | 非洲黑人性xxxx精品又粗又长| 国产乱人视频| 亚洲精品第二区| 日韩一本色道免费dvd| 亚洲欧美一区二区三区国产| 国产亚洲5aaaaa淫片| 一级黄片播放器| 日产精品乱码卡一卡2卡三| 午夜激情欧美在线| 街头女战士在线观看网站| a级毛色黄片| 亚洲乱码一区二区免费版| 干丝袜人妻中文字幕| 午夜福利网站1000一区二区三区| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 久久精品夜色国产| 国产一区二区三区av在线| 极品教师在线视频| 久久久久精品久久久久真实原创| 日本三级黄在线观看| 精品国产三级普通话版| 精品久久久久久久久久久久久| 精品一区二区免费观看| 精品久久久噜噜| 少妇的逼好多水| a级一级毛片免费在线观看| 国产亚洲午夜精品一区二区久久 | 天天躁夜夜躁狠狠久久av| 最近的中文字幕免费完整| 日产精品乱码卡一卡2卡三| 麻豆av噜噜一区二区三区| 亚洲乱码一区二区免费版| 国产精品伦人一区二区| 国产精品av视频在线免费观看| 蜜桃久久精品国产亚洲av| 亚洲国产精品专区欧美| 欧美一级a爱片免费观看看| 亚洲美女视频黄频| 麻豆乱淫一区二区| 全区人妻精品视频| 精品久久国产蜜桃| 亚洲美女搞黄在线观看| 国产伦一二天堂av在线观看| 亚洲成人中文字幕在线播放| 亚洲av免费高清在线观看| 岛国毛片在线播放| 丰满乱子伦码专区| 国产大屁股一区二区在线视频| 深夜a级毛片| 精品欧美国产一区二区三| 岛国毛片在线播放| 搡老乐熟女国产| 亚洲在线观看片| 精品久久久久久久久亚洲| 在线观看免费高清a一片| 日韩视频在线欧美| 亚洲欧美日韩无卡精品| 国产av国产精品国产| 成人一区二区视频在线观看| 99热网站在线观看| 中文字幕免费在线视频6| 自拍偷自拍亚洲精品老妇| 乱系列少妇在线播放| 亚洲精品自拍成人| 亚洲精品中文字幕在线视频 | 久久久久网色| freevideosex欧美| 久久久久久久国产电影| 久久久久久久久久黄片| 99久久精品一区二区三区| 国产视频首页在线观看| 岛国毛片在线播放| 国产亚洲精品久久久com| 天天躁夜夜躁狠狠久久av| 国产亚洲一区二区精品| 亚洲怡红院男人天堂| 一级爰片在线观看| 我要看日韩黄色一级片| 九草在线视频观看| 久久久久国产网址| 亚洲高清免费不卡视频| 午夜福利网站1000一区二区三区| 成年免费大片在线观看| 亚洲怡红院男人天堂| 午夜精品一区二区三区免费看| 男女国产视频网站| 国产毛片a区久久久久| 亚洲va在线va天堂va国产| 日本黄大片高清| 麻豆乱淫一区二区| 国产精品av视频在线免费观看| 一级爰片在线观看| 成人漫画全彩无遮挡| 久久精品国产亚洲av涩爱| 乱人视频在线观看| 看免费成人av毛片| av国产久精品久网站免费入址| 亚洲第一区二区三区不卡| 熟女人妻精品中文字幕| 欧美日韩精品成人综合77777| 久久精品国产亚洲av天美| 国产精品av视频在线免费观看| 国产日韩欧美在线精品| 日韩强制内射视频| 国产成人精品福利久久| 午夜福利视频1000在线观看| freevideosex欧美| 亚洲不卡免费看| 亚洲精品色激情综合| 18+在线观看网站| 中国国产av一级| 欧美激情久久久久久爽电影| 最近手机中文字幕大全| 免费高清在线观看视频在线观看| 精品久久久久久成人av| 久久久国产一区二区| 伦精品一区二区三区| 丰满少妇做爰视频| 在线观看av片永久免费下载| 欧美极品一区二区三区四区| 久久人人爽人人爽人人片va| 内射极品少妇av片p| 一区二区三区乱码不卡18| av网站免费在线观看视频 | ponron亚洲| 91精品国产九色| 免费不卡的大黄色大毛片视频在线观看 | 五月伊人婷婷丁香| 国产日韩欧美在线精品| 国产综合精华液| .国产精品久久| 精品人妻熟女av久视频| 日韩伦理黄色片| 国产色婷婷99| 国产黄片视频在线免费观看| 你懂的网址亚洲精品在线观看| 最近2019中文字幕mv第一页| 国内揄拍国产精品人妻在线| 欧美97在线视频| 免费大片18禁| 亚洲国产高清在线一区二区三| 麻豆精品久久久久久蜜桃| 日韩不卡一区二区三区视频在线| 秋霞伦理黄片| 卡戴珊不雅视频在线播放| 边亲边吃奶的免费视频| 蜜桃亚洲精品一区二区三区| 亚洲国产欧美在线一区| 伊人久久精品亚洲午夜| 精品欧美国产一区二区三| 少妇人妻一区二区三区视频| 精品久久久噜噜| 国产黄色小视频在线观看| 精品熟女少妇av免费看| 国产单亲对白刺激| 中文精品一卡2卡3卡4更新| 中国美白少妇内射xxxbb| 六月丁香七月| 婷婷色麻豆天堂久久| 欧美高清性xxxxhd video| 国产亚洲午夜精品一区二区久久 | 3wmmmm亚洲av在线观看| 啦啦啦韩国在线观看视频| 一区二区三区四区激情视频| 国产一级毛片七仙女欲春2| 亚洲精品中文字幕在线视频 | 蜜桃亚洲精品一区二区三区| 国产男人的电影天堂91| 国内精品一区二区在线观看| 简卡轻食公司| 亚洲欧美成人综合另类久久久| 久久久精品免费免费高清| 嫩草影院入口| 亚洲精品久久久久久婷婷小说| av黄色大香蕉| 亚洲最大成人av| 美女脱内裤让男人舔精品视频| a级毛色黄片| 美女黄网站色视频| 亚洲在久久综合| 中文乱码字字幕精品一区二区三区 | 国产老妇伦熟女老妇高清| 亚洲欧美日韩无卡精品| 精品一区二区三卡| 91久久精品国产一区二区成人| 国产永久视频网站| 亚洲图色成人| 中国国产av一级| 肉色欧美久久久久久久蜜桃 | 人妻一区二区av| 最近手机中文字幕大全| 亚洲精品日韩av片在线观看| 国产淫语在线视频| 嫩草影院入口| 亚洲精品,欧美精品| 又黄又爽又刺激的免费视频.| 国产黄色小视频在线观看| 午夜福利网站1000一区二区三区| 99热这里只有精品一区| 综合色av麻豆| 你懂的网址亚洲精品在线观看| 免费看日本二区| 日韩av免费高清视频| 国产精品一区二区在线观看99 | 人妻夜夜爽99麻豆av| 极品教师在线视频| 国国产精品蜜臀av免费| 人人妻人人澡人人爽人人夜夜 | 精品一区二区三区人妻视频| 国产91av在线免费观看| 国产精品爽爽va在线观看网站| 夜夜看夜夜爽夜夜摸| 国产精品蜜桃在线观看| 久久久精品94久久精品| 亚洲成人av在线免费| 女人被狂操c到高潮| 久久人人爽人人爽人人片va| 亚洲经典国产精华液单| 少妇人妻精品综合一区二区| 美女黄网站色视频| 黄色配什么色好看| 国产美女午夜福利| 免费电影在线观看免费观看| 少妇熟女aⅴ在线视频| 三级毛片av免费| 精品一区二区免费观看| 熟妇人妻久久中文字幕3abv| 能在线免费观看的黄片| 国产伦精品一区二区三区视频9| av在线亚洲专区| 精品一区二区三区视频在线| 男女下面进入的视频免费午夜| 国产日韩欧美在线精品| 乱人视频在线观看| 最近中文字幕2019免费版| 成人鲁丝片一二三区免费| 亚洲一级一片aⅴ在线观看|