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

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

    2017-08-31 11:55:45張小明張梅軍
    振動(dòng)與沖擊 2017年16期
    關(guān)鍵詞:時(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

    猜你喜歡
    時(shí)域振幅幅值
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    基于S變換的交流電網(wǎng)幅值檢測(cè)系統(tǒng)計(jì)算機(jī)仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    滬市十大振幅
    正序電壓幅值檢測(cè)及諧波抑制的改進(jìn)
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    低壓電力線信道脈沖噪聲的幅值與寬度特征
    精品不卡国产一区二区三区| 国产亚洲精品综合一区在线观看| 国产伦在线观看视频一区| 内射极品少妇av片p| 欧美日本视频| 成人亚洲精品av一区二区| 欧美高清性xxxxhd video| 韩国av在线不卡| 麻豆av噜噜一区二区三区| 成人性生交大片免费视频hd| avwww免费| 久久热精品热| 国产在视频线在精品| 18禁裸乳无遮挡免费网站照片| 久久人人精品亚洲av| 国产高清视频在线播放一区| av在线蜜桃| 波多野结衣高清无吗| 神马国产精品三级电影在线观看| 伦理电影大哥的女人| 中文资源天堂在线| 欧美人与善性xxx| 国语自产精品视频在线第100页| 日本黄色视频三级网站网址| 全区人妻精品视频| 国产色婷婷99| 国产在线男女| 看免费成人av毛片| 人妻丰满熟妇av一区二区三区| 精品无人区乱码1区二区| 女人十人毛片免费观看3o分钟| 99热网站在线观看| 亚洲 国产 在线| 欧美激情国产日韩精品一区| 在线播放国产精品三级| 欧美xxxx性猛交bbbb| av视频在线观看入口| 丰满人妻一区二区三区视频av| 日本三级黄在线观看| 日韩高清综合在线| 国产一区二区在线av高清观看| av天堂在线播放| 看十八女毛片水多多多| 免费观看精品视频网站| 欧美丝袜亚洲另类 | 18禁裸乳无遮挡免费网站照片| 床上黄色一级片| 国产 一区精品| 亚洲专区国产一区二区| 欧美最黄视频在线播放免费| 又爽又黄无遮挡网站| 嫩草影院入口| 久久久久久国产a免费观看| 亚洲成人中文字幕在线播放| 在线a可以看的网站| 99精品在免费线老司机午夜| 性欧美人与动物交配| 99久久精品国产国产毛片| 观看美女的网站| 国产国拍精品亚洲av在线观看| 少妇被粗大猛烈的视频| 男女啪啪激烈高潮av片| 久久香蕉精品热| 特级一级黄色大片| 日日夜夜操网爽| 99热这里只有精品一区| aaaaa片日本免费| 日本一二三区视频观看| 亚洲av熟女| 亚洲狠狠婷婷综合久久图片| 88av欧美| 国产av不卡久久| 国产精品亚洲一级av第二区| 中国美白少妇内射xxxbb| 国产 一区 欧美 日韩| 又爽又黄无遮挡网站| 人妻少妇偷人精品九色| 97热精品久久久久久| 嫩草影院新地址| 成人特级av手机在线观看| 99久久中文字幕三级久久日本| 国产真实伦视频高清在线观看 | 69人妻影院| 91久久精品国产一区二区三区| 丰满人妻一区二区三区视频av| 国产不卡一卡二| 此物有八面人人有两片| 色哟哟哟哟哟哟| 精品不卡国产一区二区三区| 亚洲最大成人中文| 午夜免费男女啪啪视频观看 | 干丝袜人妻中文字幕| 乱系列少妇在线播放| 久久久久九九精品影院| 免费大片18禁| 亚洲精品国产成人久久av| 国产精品国产三级国产av玫瑰| www日本黄色视频网| 日本免费一区二区三区高清不卡| 亚洲经典国产精华液单| 国产探花极品一区二区| 天堂av国产一区二区熟女人妻| 国产av一区在线观看免费| 成熟少妇高潮喷水视频| 欧美xxxx性猛交bbbb| 久9热在线精品视频| 日本一二三区视频观看| 男女视频在线观看网站免费| 制服丝袜大香蕉在线| 久久久久久国产a免费观看| 88av欧美| 亚洲综合色惰| 97人妻精品一区二区三区麻豆| 亚洲 国产 在线| 久久精品91蜜桃| 婷婷亚洲欧美| 91在线观看av| 日韩大尺度精品在线看网址| 亚洲人成网站在线播| 大型黄色视频在线免费观看| h日本视频在线播放| 黄色配什么色好看| 亚洲午夜理论影院| 黄色欧美视频在线观看| 99热这里只有精品一区| 免费在线观看成人毛片| 国产精品日韩av在线免费观看| 国产三级中文精品| 国产探花极品一区二区| 中文字幕高清在线视频| 国产色爽女视频免费观看| 成人亚洲精品av一区二区| 欧美3d第一页| 国产日本99.免费观看| 免费观看的影片在线观看| videossex国产| 69av精品久久久久久| 国产伦人伦偷精品视频| 淫秽高清视频在线观看| 亚洲精品粉嫩美女一区| 免费高清视频大片| 亚州av有码| 丰满人妻一区二区三区视频av| 精品一区二区免费观看| 国内揄拍国产精品人妻在线| 欧美一区二区精品小视频在线| 欧美区成人在线视频| 网址你懂的国产日韩在线| 国产精品亚洲美女久久久| 午夜视频国产福利| 中文字幕久久专区| 国产在线精品亚洲第一网站| 久久精品国产鲁丝片午夜精品 | 国产一区二区三区视频了| 九九在线视频观看精品| 亚洲色图av天堂| 国产精品一区二区三区四区免费观看 | 亚洲精品日韩av片在线观看| 日日啪夜夜撸| 久久精品人妻少妇| 国产一区二区三区视频了| 亚洲成人免费电影在线观看| 一个人看的www免费观看视频| 欧美人与善性xxx| 熟妇人妻久久中文字幕3abv| 国产亚洲精品av在线| 日韩欧美在线乱码| 人人妻人人看人人澡| 波多野结衣巨乳人妻| 可以在线观看毛片的网站| 亚洲欧美日韩高清在线视频| 色尼玛亚洲综合影院| 内地一区二区视频在线| 精品99又大又爽又粗少妇毛片 | 亚洲专区中文字幕在线| 国产爱豆传媒在线观看| 国产高清视频在线播放一区| 国产一级毛片七仙女欲春2| av在线蜜桃| 悠悠久久av| 12—13女人毛片做爰片一| 三级国产精品欧美在线观看| 色哟哟哟哟哟哟| 亚洲第一区二区三区不卡| 免费观看在线日韩| 国产午夜福利久久久久久| 国产精品综合久久久久久久免费| 悠悠久久av| 悠悠久久av| 22中文网久久字幕| 一进一出抽搐gif免费好疼| 国产精品综合久久久久久久免费| 日韩一本色道免费dvd| 又粗又爽又猛毛片免费看| 大型黄色视频在线免费观看| 亚洲精品亚洲一区二区| av天堂在线播放| 精品久久久久久,| 午夜福利欧美成人| 精品国产三级普通话版| 久久这里只有精品中国| 久久99热这里只有精品18| 一进一出好大好爽视频| 97超视频在线观看视频| 国产高清视频在线观看网站| 少妇的逼好多水| 国产国拍精品亚洲av在线观看| 国产精品三级大全| 中文字幕熟女人妻在线| or卡值多少钱| 一个人免费在线观看电影| 日本爱情动作片www.在线观看 | 久久久久久久午夜电影| 欧美性猛交╳xxx乱大交人| ponron亚洲| 午夜影院日韩av| 日韩中字成人| 一进一出好大好爽视频| 久久人人爽人人爽人人片va| 丰满乱子伦码专区| 一级毛片久久久久久久久女| 欧美激情在线99| 一区二区三区免费毛片| 春色校园在线视频观看| 麻豆国产av国片精品| 午夜激情欧美在线| 美女高潮的动态| 99热这里只有是精品在线观看| 免费看日本二区| 亚洲第一电影网av| 久久国产精品人妻蜜桃| 99精品在免费线老司机午夜| 成人精品一区二区免费| 91狼人影院| 国产乱人伦免费视频| 亚洲电影在线观看av| 日韩人妻高清精品专区| 无遮挡黄片免费观看| 日本 av在线| 亚洲av美国av| 久久久精品欧美日韩精品| 91在线精品国自产拍蜜月| 午夜老司机福利剧场| 国产淫片久久久久久久久| 亚洲人成伊人成综合网2020| 91精品国产九色| 欧美高清性xxxxhd video| 99热这里只有是精品50| 免费一级毛片在线播放高清视频| 老司机深夜福利视频在线观看| 色综合色国产| 亚洲国产高清在线一区二区三| 99久久久亚洲精品蜜臀av| 国产精品福利在线免费观看| 久久婷婷人人爽人人干人人爱| 波多野结衣高清作品| 综合色av麻豆| 99久久无色码亚洲精品果冻| 久久国内精品自在自线图片| 亚洲欧美激情综合另类| 亚洲avbb在线观看| 成人美女网站在线观看视频| 韩国av一区二区三区四区| 国产乱人伦免费视频| 日本黄大片高清| 精品人妻1区二区| 婷婷六月久久综合丁香| 国产一区二区在线av高清观看| 校园人妻丝袜中文字幕| 俺也久久电影网| 美女高潮的动态| 中文在线观看免费www的网站| a级一级毛片免费在线观看| 尾随美女入室| 全区人妻精品视频| 国产男靠女视频免费网站| 草草在线视频免费看| 国产av不卡久久| 99国产精品一区二区蜜桃av| 国产淫片久久久久久久久| 久久草成人影院| 黄色欧美视频在线观看| 欧美激情久久久久久爽电影| 不卡一级毛片| 久久99热这里只有精品18| 狠狠狠狠99中文字幕| 熟女人妻精品中文字幕| 国产在线男女| 亚洲美女黄片视频| 国产高清有码在线观看视频| 国产精品久久电影中文字幕| 国产熟女欧美一区二区| 亚洲18禁久久av| 中文字幕精品亚洲无线码一区| 在线观看舔阴道视频| 国产高清激情床上av| 在线免费十八禁| 中文字幕av在线有码专区| 婷婷色综合大香蕉| 毛片女人毛片| 国产激情偷乱视频一区二区| 亚洲自拍偷在线| 成人无遮挡网站| 精品国内亚洲2022精品成人| 亚洲真实伦在线观看| 国产精品久久久久久av不卡| 国产色婷婷99| 午夜影院日韩av| 性色avwww在线观看| 99久久精品一区二区三区| 国产精品无大码| 一边摸一边抽搐一进一小说| 精品久久久久久久久久免费视频| 国产女主播在线喷水免费视频网站 | 999久久久精品免费观看国产| 久久精品综合一区二区三区| 91麻豆av在线| 婷婷精品国产亚洲av| 久久精品影院6| x7x7x7水蜜桃| 亚洲三级黄色毛片| 91久久精品电影网| 亚洲熟妇熟女久久| 91在线精品国自产拍蜜月| 精品一区二区三区视频在线| 欧洲精品卡2卡3卡4卡5卡区| 日韩在线高清观看一区二区三区 | 五月伊人婷婷丁香| 亚洲av日韩精品久久久久久密| 97热精品久久久久久| 91麻豆av在线| 在线播放无遮挡| av在线蜜桃| 国产精品一及| 免费搜索国产男女视频| 久久99热6这里只有精品| 可以在线观看毛片的网站| 亚洲成人免费电影在线观看| 一a级毛片在线观看| 淫秽高清视频在线观看| 亚洲欧美日韩高清专用| 精品久久久久久成人av| 日日摸夜夜添夜夜添av毛片 | 久久天躁狠狠躁夜夜2o2o| 久久久久性生活片| 精品久久久久久久久久免费视频| 精品福利观看| 18禁黄网站禁片免费观看直播| 久久久久久大精品| 国产精品嫩草影院av在线观看 | 国产成人一区二区在线| 亚洲专区中文字幕在线| 亚洲国产欧洲综合997久久,| 18禁裸乳无遮挡免费网站照片| 九九爱精品视频在线观看| 日韩国内少妇激情av| av天堂在线播放| 综合色av麻豆| bbb黄色大片| 国产av麻豆久久久久久久| 我的老师免费观看完整版| 国产精品国产三级国产av玫瑰| 日本 欧美在线| 99热这里只有精品一区| 国产精品国产高清国产av| 97人妻精品一区二区三区麻豆| 俄罗斯特黄特色一大片| x7x7x7水蜜桃| 午夜日韩欧美国产| 男女啪啪激烈高潮av片| 色哟哟·www| 可以在线观看的亚洲视频| 成人一区二区视频在线观看| 午夜福利在线观看免费完整高清在 | 午夜亚洲福利在线播放| 两人在一起打扑克的视频| 亚洲国产日韩欧美精品在线观看| 男女啪啪激烈高潮av片| 久久久久久大精品| 一级黄片播放器| 免费av观看视频| 男人舔奶头视频| 亚洲内射少妇av| 99在线视频只有这里精品首页| 久久人妻av系列| 亚洲人成网站在线播放欧美日韩| 日本 av在线| 99国产精品一区二区蜜桃av| 老司机午夜福利在线观看视频| 久久久久免费精品人妻一区二区| 亚洲专区中文字幕在线| 国产精品一区二区三区四区久久| 男人的好看免费观看在线视频| 黄色丝袜av网址大全| av天堂中文字幕网| 久久久久久久久大av| 97热精品久久久久久| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 国产主播在线观看一区二区| 91久久精品国产一区二区成人| 亚洲av电影不卡..在线观看| 欧美激情久久久久久爽电影| 久久精品久久久久久噜噜老黄 | 国产精品日韩av在线免费观看| 亚洲成人久久爱视频| 亚洲中文日韩欧美视频| 日韩人妻高清精品专区| 亚洲精品粉嫩美女一区| 可以在线观看的亚洲视频| 午夜日韩欧美国产| 2021天堂中文幕一二区在线观| 亚洲av一区综合| 亚洲人成网站在线播| 国产高清不卡午夜福利| 亚洲av美国av| 能在线免费观看的黄片| 欧美日本视频| 一级a爱片免费观看的视频| 国产精品久久久久久av不卡| 中文字幕人妻熟人妻熟丝袜美| 午夜福利在线观看免费完整高清在 | 久久久久久久精品吃奶| 高清在线国产一区| 国内精品美女久久久久久| 色综合站精品国产| 久久精品国产99精品国产亚洲性色| 亚洲国产精品合色在线| 色综合亚洲欧美另类图片| 亚洲图色成人| av.在线天堂| 久久精品综合一区二区三区| 国产成人av教育| 人妻制服诱惑在线中文字幕| 亚洲综合色惰| 国产精品美女特级片免费视频播放器| 欧美日韩国产亚洲二区| 色综合亚洲欧美另类图片| 精品人妻熟女av久视频| 日韩国内少妇激情av| 高清日韩中文字幕在线| 少妇人妻精品综合一区二区 | 婷婷色综合大香蕉| 桃红色精品国产亚洲av| 在线国产一区二区在线| 国产主播在线观看一区二区| 真人一进一出gif抽搐免费| 亚洲无线观看免费| 国产精品久久久久久精品电影| av女优亚洲男人天堂| 国产午夜精品久久久久久一区二区三区 | 人人妻人人看人人澡| 欧美bdsm另类| 99精品久久久久人妻精品| 国产精品乱码一区二三区的特点| 亚洲一级一片aⅴ在线观看| 日本撒尿小便嘘嘘汇集6| .国产精品久久| 亚洲精品国产成人久久av| 欧美高清性xxxxhd video| 悠悠久久av| 亚洲七黄色美女视频| 亚洲在线观看片| 亚洲国产精品sss在线观看| 亚洲欧美清纯卡通| 搡老熟女国产l中国老女人| 国产精品人妻久久久影院| 18禁在线播放成人免费| 午夜久久久久精精品| 一区二区三区高清视频在线| 免费看美女性在线毛片视频| 真实男女啪啪啪动态图| av中文乱码字幕在线| 夜夜爽天天搞| 校园春色视频在线观看| 免费看光身美女| 五月玫瑰六月丁香| 免费看av在线观看网站| 免费av不卡在线播放| 22中文网久久字幕| 精品午夜福利在线看| 久久久久性生活片| 日本成人三级电影网站| 在线观看一区二区三区| 麻豆成人午夜福利视频| 久久久久久大精品| 亚洲精品乱码久久久v下载方式| 日本熟妇午夜| 国产午夜精品久久久久久一区二区三区 | 亚洲成人久久爱视频| 亚洲色图av天堂| 国产蜜桃级精品一区二区三区| 无人区码免费观看不卡| 中文字幕av在线有码专区| 国产高潮美女av| 免费人成视频x8x8入口观看| 在线播放国产精品三级| 日韩精品青青久久久久久| 九九热线精品视视频播放| 国产精品永久免费网站| 小说图片视频综合网站| 高清日韩中文字幕在线| 欧美丝袜亚洲另类 | www.色视频.com| 高清毛片免费观看视频网站| 99久久成人亚洲精品观看| 麻豆国产av国片精品| 国产伦精品一区二区三区四那| 久久亚洲精品不卡| 亚洲欧美清纯卡通| 国产精品99久久久久久久久| 精华霜和精华液先用哪个| 狠狠狠狠99中文字幕| 直男gayav资源| 欧美人与善性xxx| 日韩欧美精品v在线| 啪啪无遮挡十八禁网站| 禁无遮挡网站| 亚洲av五月六月丁香网| or卡值多少钱| 色av中文字幕| 久久欧美精品欧美久久欧美| 免费看美女性在线毛片视频| 日本在线视频免费播放| av在线观看视频网站免费| 欧美性猛交╳xxx乱大交人| 国产精品一区二区免费欧美| 大型黄色视频在线免费观看| 国模一区二区三区四区视频| 日本一本二区三区精品| ponron亚洲| 午夜老司机福利剧场| 精品一区二区三区av网在线观看| 国产不卡一卡二| 日本 欧美在线| 国产精品亚洲美女久久久| 草草在线视频免费看| 最新在线观看一区二区三区| 又黄又爽又免费观看的视频| 露出奶头的视频| 91麻豆av在线| 日日啪夜夜撸| 亚洲黑人精品在线| 亚洲男人的天堂狠狠| 久久香蕉精品热| avwww免费| xxxwww97欧美| av福利片在线观看| 欧美成人性av电影在线观看| 真实男女啪啪啪动态图| 午夜a级毛片| 男女视频在线观看网站免费| 一个人看视频在线观看www免费| 成年女人永久免费观看视频| 国产v大片淫在线免费观看| 日韩欧美国产一区二区入口| 国产精品女同一区二区软件 | 亚洲成人精品中文字幕电影| 欧美xxxx黑人xx丫x性爽| 人妻制服诱惑在线中文字幕| 久久久成人免费电影| 男人和女人高潮做爰伦理| 国产熟女欧美一区二区| 99久久无色码亚洲精品果冻| 99热这里只有是精品50| 国产精品久久久久久亚洲av鲁大| 麻豆成人午夜福利视频| 亚洲午夜理论影院| 精品久久久久久久久av| 成熟少妇高潮喷水视频| .国产精品久久| 女生性感内裤真人,穿戴方法视频| 老司机福利观看| 婷婷精品国产亚洲av在线| 日韩国内少妇激情av| 我要看日韩黄色一级片| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品日韩av片在线观看| 亚洲精品影视一区二区三区av| 成人综合一区亚洲| 村上凉子中文字幕在线| 女人十人毛片免费观看3o分钟| 伦精品一区二区三区| 久久久成人免费电影| 日韩 亚洲 欧美在线| 亚洲电影在线观看av| 听说在线观看完整版免费高清| 亚洲不卡免费看| 此物有八面人人有两片| 久久久国产成人免费| 欧美+日韩+精品| www.色视频.com| 97热精品久久久久久| 白带黄色成豆腐渣| 国产精品免费一区二区三区在线| 色视频www国产| 国产三级中文精品| 亚洲精品久久国产高清桃花| av福利片在线观看| 国产精品精品国产色婷婷| 日韩在线高清观看一区二区三区 | 国产精品无大码| 亚洲精品久久国产高清桃花| 免费av不卡在线播放| 久久久久久大精品| 99久久成人亚洲精品观看| 中文字幕av成人在线电影| 国产成人aa在线观看| 日韩中文字幕欧美一区二区| 高清毛片免费观看视频网站| 久久久久免费精品人妻一区二区| 黄色视频,在线免费观看| 99热这里只有是精品50| 国产精品人妻久久久影院| 一卡2卡三卡四卡精品乱码亚洲|