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

    模式自適應(yīng)小波時能密度法及其在微差爆破振動信號分析中的應(yīng)用

    2014-06-07 05:55:20凌同華劉浩然
    煤炭學(xué)報 2014年10期
    關(guān)鍵詞:小波基微差小波

    張 勝,凌同華,劉浩然,曹 峰

    (1.長沙理工大學(xué)土木與建筑學(xué)院,湖南長沙 410114;2.長沙理工大學(xué)橋梁工程安全控制省部共建教育部重點實驗室,湖南長沙 410114)

    模式自適應(yīng)小波時能密度法及其在微差爆破振動信號分析中的應(yīng)用

    張 勝1,2,凌同華1,劉浩然1,曹 峰1

    (1.長沙理工大學(xué)土木與建筑學(xué)院,湖南長沙 410114;2.長沙理工大學(xué)橋梁工程安全控制省部共建教育部重點實驗室,湖南長沙 410114)

    如何解決適合爆破振動信號特征的爆破振動信號小波基構(gòu)造、添加及其實現(xiàn)等問題,始終是困擾廣大研究者的難題。為準(zhǔn)確確定微差爆破振動信號實際微差延遲時間,在模式自適應(yīng)小波原理的基礎(chǔ)上,構(gòu)造了從實測微差爆破振動信號中分離出子信號的模式自適應(yīng)小波基,然后將該模式自適應(yīng)小波基添加到Matlab小波分析工具箱中,并提出了一種新的基于爆破振動信號分析的模式自適應(yīng)小波時能密度法。以實驗設(shè)計的多段微差爆破振動信號進(jìn)行模式自適應(yīng)小波時能密度法分析為例,驗證了該方法的可行性。該方法亦可用于分析實測微差爆破振動信號在實際微差延遲時間的確定,并與db8小波法和模式自適應(yīng)小波法進(jìn)行比較,結(jié)果表明模式自適應(yīng)小波時能密度法具有更好的應(yīng)用效果,從而為爆破振動信號分析在實際工程中的應(yīng)用提供了參考。

    爆破振動;微差時間;模式自適應(yīng)小波時能密度法;小波基構(gòu)造;小波基添加

    微差爆破實踐中,普通雷管起爆都有延期問題,即便是所謂的高精度雷管,也有±10 ms的誤差,由于延期雷管不能準(zhǔn)確按設(shè)計延期時間起爆,導(dǎo)致微差爆破起爆順序混亂,影響了實際微差爆破效果與安全[1-2]。因此,研究與確定微差爆破振動信號實際微差延期時間對改進(jìn)爆破方式、優(yōu)化爆破參數(shù)和控制爆破效果具有重要的指導(dǎo)意義。

    從微差爆破振動信號實測資料分析,微差爆破振動信號具有持時短、衰減快以及時頻局部化的特點,屬于典型的非平穩(wěn)信號。長期以來限于理論的發(fā)展,只好將微差爆破振動信號視為簡單的平穩(wěn)信號,分析和處理平穩(wěn)信號最常用的方法是Fourier分析。近年來,隨著科學(xué)技術(shù)的發(fā)展,小波時頻局部化分析已廣泛應(yīng)用于工程技術(shù)領(lǐng)域,用小波時頻局部化分析和處理非平穩(wěn)信號已激起人們很高的熱忱[3]。小波時頻局部化分析方法是一種窗口大小固定但其形狀、時頻窗都可改變的時頻局部化分析方法[4]。小波時頻局部化分析技術(shù)在微差爆破振動信號特征提取、頻譜分析、微差延期時間識別、波形預(yù)測和反應(yīng)譜分析等方面具有良好的應(yīng)用效果[5]。然而,采用經(jīng)典小波進(jìn)行小波時頻局部化分析時,往往是根據(jù)研究對象的特點從已有的Matlab小波分析工具箱中選擇一個比較合適的小波基,雖然可以滿足微差爆破振動信號分析處理的一般要求,但由于小波基以及尺度選取的不確定性事實,已嚴(yán)重影響了小波時頻局部化分析的實際應(yīng)用效果[6-8]。因此,如何根據(jù)微差爆破振動信號的特征選擇或構(gòu)造與特征信號更好匹配的最優(yōu)小波基、提高信號處理的效率與質(zhì)量是小波理論在爆破振動信號分析中亟待解決的問題[9-13]。

    本文以此為切入點,在模式自適應(yīng)小波原理的基礎(chǔ)上,運(yùn)用時頻轉(zhuǎn)換技術(shù),將實測微差爆破振動信號分離出子信號并對其進(jìn)行模式自適應(yīng)波形匹配,構(gòu)造出與子信號相似度高的模式自適應(yīng)爆破振動信號小波基(簡稱模式自適應(yīng)小波,下同),然后將該模式自適應(yīng)小波添加到Matlab小波分析工具箱中,并提出一種新的基于爆破振動信號小波時頻局部化分析的模式自適應(yīng)小波時能密度法,將其用于實驗設(shè)計信號和實測微差爆破振動信號的微差延期時間識別,取得了滿意的結(jié)果,從而為爆破振動信號分析在實際工程中的應(yīng)用提供參考。

    1 模式自適應(yīng)小波時能密度法

    1.1 連續(xù)小波變換

    函數(shù)ψ(t)∈L2(R)(L2(R)為平方可積的信號空間)的連續(xù)傅里葉變換為ψ^(ω)[14-16]。當(dāng)ψ^(ω)滿足小波允許條件:

    時,稱ψ(t)為一個小波基。若將該小波基ψ(t)經(jīng)伸縮和平移后,就得到一個小波基序列。

    對于實數(shù)a,b且a≠0的情況下,小波基序列為

    式中,a為尺度因子;b為平移因子。

    對于微差爆破振動信號f(t)∈L2(R)的連續(xù)小波變換為

    由式(2)可以看出,小波變換的實質(zhì)是將微差爆破振動信號f(t)與小波基序列ψa,b(t)做卷積,求得微差爆破振動信號f(t)在小波基序列ψa,b(t)上的投影。若要微差爆破振動信號中的特征成分在小波系數(shù)時頻譜上更好地顯示,則需要將微差爆破振動信號f(t)在小波基序列ψa,b(t)上的投影系數(shù)盡可能地最大化,也就是說小波基的波形應(yīng)與微差爆破振動信號特征波形具有較高的相似度。

    1.2 模式自適應(yīng)小波

    模式自適應(yīng)小波是利用最小平方優(yōu)化方法來構(gòu)建一個與給定模式近似的小波[17]。具體方法有多項式逼近法(polynomial)和常數(shù)正交函數(shù)空間投影法(orth.a(chǎn)nd constants)2種。如果信號特征模式比較簡單,則可以選擇多項式逼近法;如果信號特征模式比較復(fù)雜,則可以選擇常數(shù)正交函數(shù)空間投影法。參數(shù)規(guī)則度(regularity)定義了模式自適應(yīng)小波在[0, 1]上的邊界約束,可以是“無”或“連續(xù)”或“可微”。

    基于上述構(gòu)造模式自適應(yīng)小波的思想[18-19],提出了以下步驟:

    (1)根據(jù)給定的檢測模式f,產(chǎn)生模式自適應(yīng)小波ψf(t)。不失一般性,假設(shè)supp(f)=[0,1],則有

    (2)檢測信號上所有的報警,對任意b和a>0的問題,搜索信號小波能量的局部極大值。

    (3)檢測和丟棄所有錯誤的報警,這一規(guī)則必須應(yīng)用于確定每個報警是否錯誤。

    1.3 模式自適應(yīng)小波時能密度法

    若模式自適應(yīng)小波滿足式(1)的容許性條件,則說明模式自適應(yīng)小波變換是完備的并且保持能量守恒,根據(jù)Moyal定理,有下式成立

    式(3)表明模式自適應(yīng)小波變換系數(shù)平方的積分同被分析的微差爆破振動信號的能量成正比。因此,可以把看作是模式自適應(yīng)小波時間-尺度平面上的模式自適應(yīng)小波能量密度函數(shù),即給出了以尺度因子a和平移因子b為中心的、尺度間隔為Δa、時間間隔為Δb的能量。根據(jù)模式自適應(yīng)小波能量密度的概念,式(3)也可以寫成

    在模式自適應(yīng)小波變換中,尺度因子a在某種程度上對應(yīng)于實際頻率f,因此,式(5)給出了微差爆破振動信號所有頻帶的能量隨平移因子b的分布情況,稱為模式自適應(yīng)小波時能密度函數(shù)。

    2 模式自適應(yīng)小波時能密度法的構(gòu)造、添加與實現(xiàn)

    2.1 可行性分析

    從理論上講,任何函數(shù)f(x)∈L2(R)只要滿足小波基允許條件就能作為一個新的小波基進(jìn)行選用,但選擇具有緊支撐性、一定正則性、曲線外形與微差爆破振動信號有較好相似性的函數(shù)作為小波基函數(shù),能在時頻域上都取得良好的時頻局部化分析效果。從某一實測微差爆破振動信號中分離出的子信號速度-時間曲線如圖1所示及其相應(yīng)的功率譜曲線見圖2[5]。從圖1,2可以看出,爆破振動子信號在時頻域上能量都比較集中、且衰減速度較快,也就是說爆破振動子信號在時頻域上均具有良好的時頻局部化特性。因此,可以根據(jù)微差爆破振動信號的特點,由爆破振動子信號來進(jìn)行模式自適應(yīng)波形匹配,構(gòu)造出與爆破振動子信號相似度高的模式自適應(yīng)小波。

    圖1 子信號的速度-時間曲線Fig.1 The velocity vs time curve of blast vibration sub-signal

    圖2 子信號的頻譜圖Fig.2 Spectrogram of blast vibration sub-signal

    2.2 模式自適應(yīng)小波的構(gòu)造

    在Matlab語言環(huán)境下,通過利用Wavelet Toolbox工具箱中的New Wavelet for CWT界面對如圖1所示的爆破振動子信號進(jìn)行模式自適應(yīng)波形匹配,構(gòu)造出與爆破振動子信號相似度高的sub_signal.mat函數(shù)(其中與原始信號的相關(guān)系數(shù)高達(dá)0.99)。構(gòu)造好的模式自適應(yīng)小波如圖3所示。

    圖3 子信號的模式自適應(yīng)小波設(shè)計Fig.3 Pattern adapted admissible wavelet design of blast vibration sub-signal

    2.3 模式自適應(yīng)小波的添加

    在Matlab編程語言工作平臺上,可以使用wavemngr命令添加該模式自適應(yīng)小波[20]。添加好的模式自適應(yīng)小波顯示如圖4所示。

    圖4 子信號的模式自適應(yīng)小波顯示Fig.4 Pattern adapted admissible wavelet display of blast vibration sub-signal

    2.4 模式自適應(yīng)小波時能密度法的實現(xiàn)

    2.4.1 模式自適應(yīng)小波的實現(xiàn)

    爆破振動信號延期時間識別在時間上觀察得愈細(xì)致,就愈需要壓縮觀察范圍,即用高頻小波作細(xì)致的觀察和提取,故應(yīng)選用較小的尺度。模式自適應(yīng)小波基函數(shù)支撐區(qū)間為[0,1],Matlab中的CWT算法無法實現(xiàn)模式自適應(yīng)小波基函數(shù)尺度小于1的情況,因此,有必要采用數(shù)值積分的方法。

    將式(2)中的時間t和平移因子b作離散處理,按照數(shù)字信號處理的習(xí)慣把時間t離散為jTs,把平移因子b離散為iTs,其中Ts是采樣時間,則式(2)可變?yōu)?/p>

    其中,dt為采樣時間間隔,取dt=Ts,式(6)可變?yōu)?/p>

    其中,f(jdt)為微差爆破振動信號按照采樣時間間隔dt連續(xù)采集的數(shù)據(jù)。式(7)小波變換的結(jié)果可存放在計算機(jī)設(shè)定數(shù)組中,即實現(xiàn)了模式自適應(yīng)小波變換。

    2.4.2 模式自適應(yīng)小波時能密度法的實現(xiàn)

    采用模式自適應(yīng)小波作為基函數(shù),選取合適的起始尺度以及尺度間隔,對實驗設(shè)計信號和實測微差爆破振動信號進(jìn)行模式自適應(yīng)小波時能密度法分析,根據(jù)模式自適應(yīng)小波時能密度曲線中出現(xiàn)的奇異點位置即可得到各段雷管的實際起爆時刻,從而可以識別出微差爆破的實際延遲時間。模式自適應(yīng)小波時能密度法的流程如圖5所示。

    圖5 模式自適應(yīng)小波時能密度法的流程Fig.5 Flowchart of pattern adapted wavelet time-energy density method

    3 實驗設(shè)計與分析

    3.1 實驗設(shè)計

    圖6為爆破振動子信號在0~80 ms內(nèi)每隔20 ms進(jìn)行一段波形疊加而構(gòu)成的多段微差爆破振動信號,各分段微差爆破振動信號出現(xiàn)的時刻分別為0,20,40,60,80 ms,若多段微差爆破振動信號的采樣頻率為5 000 Hz,則多段微差爆破振動信號的采用時間間隔為1/5 000 s=0.2 ms,這樣各分段微差爆破振動信號出現(xiàn)的時刻又可以用相應(yīng)的時間序列位置來表示,也就是說各分段微差爆破振動信號出現(xiàn)的時間序列位置分別為0,100,200,300,400。若將最低段次雷管產(chǎn)生的多段微差爆破振動信號表示為x(n),則其他各分段微差爆破振動信號可以分別表示為x(n-100),x(n-200),x(n-300),x(n-400),疊加后的多段微差爆破振動信號則可以表示為

    式中,y(n)為疊加后的多段微差爆破振動信號;x(nnm)為不同微差間隔下的分段微差爆破振動信號;m為爆破振動子信號的條數(shù)。

    圖6 信號y(n)的速度-時間曲線Fig.6 The velocity vs time curve of signal y(n)

    3.2 結(jié)果分析

    由于多段微差爆破振動信號y(n)是在爆破振動子信號x(n)的基礎(chǔ)上進(jìn)行4次信號疊加而成,每次爆破振動子信號的疊加都會引起多段微差爆破振動信號y(n)的局部能量的突變,因此,可采用小波變換模極大值法識別出每次爆破振動子信號疊加的時刻,從而確定多段微差爆破振動信號y(n)中各相鄰爆破振動子信號間的實際微差延期時間間隔。目前,在微差爆破振動信號小波變換模極大值法中用得最多的是db8小波,本文也選用db8小波[21]。在Matlab小波分析工具箱中分別選用db8小波、模式自適應(yīng)小波對多段微差爆破振動信號y(n)進(jìn)行連續(xù)小波變換并取模值,結(jié)果見表1和圖7(a),(b)。圖7(c)為采用模式自適應(yīng)小波作為基函數(shù),尺度因子下限為0.1、上限為10和間隔為0.1時多段微差爆破振動信號y(n)的模式自適應(yīng)小波時能密度曲線。

    表1 采用3種方法對實驗信號的識別效果比較Table 1 Comparison of identification effect for experiment signals by three different methods

    圖7 信號y(n)的db8小波、模式自適應(yīng)小波變換模值以及模式自適應(yīng)小波時能密度曲線Fig.7 db8 wavelet,pattern adapted wavelet transform modulus and pattern adapted wavelet time-energy density curve of signal y(n)

    從表1、圖7可以清楚地看出,多段微差爆破振動信號y(n)的db8小波、模式自適應(yīng)小波變換模極大值以及模式自適應(yīng)小波時能密度曲線中均出現(xiàn)了5個奇異點,位置分別為0.044 4,0.064 4,0.084 4, 0.104 4,0.124 4和0.042 4,0.062 4,0.082 4, 0.102 4,0.122 4以及0.042 4,0.062 4,0.082 4, 0.102 4,0.122 4 s,表明該多段微差爆破振動信號是由5段爆破振動子信號疊加而成。由于微差爆破延期時間間隔是前后兩段雷管起爆時刻間的時間間隔,若將第1個奇異點位置作為最低段次雷管的起爆時刻,則該批次雷管的實際起爆時刻均為0,20,40,60, 80 ms,段間微差延期時間間隔均為20 ms。由此可以表明,db8小波法、模式自適應(yīng)小波法以及模式自適應(yīng)小波時能密度法對多段微差爆破振動信號y(n)中各分段爆破振動子信號的延期時間間隔識別效果較理想,可以有效識別多段微差爆破振動信號中的各奇異點位置。同時,方法2和方法3得到的多段微差爆破振動信號奇異點位置是完全相同的,由此表明模式自適應(yīng)小波時能密度法分析微差爆破振動信號的設(shè)計時間間隔是完全可靠的。

    4 實例分析

    采用爆破振動分析儀和速度傳感器對某地下礦進(jìn)行多次微差爆破振動信號測試,圖8為從中抽取的一條爆破振動速度-時間曲線,分別采用db8小波法、模式自適應(yīng)小波法以及模式自適應(yīng)小波時能密度法(其中尺度下限為0.1、上限為10和間隔為0.1)對其進(jìn)行分析處理,結(jié)果如圖9所示。

    圖8 爆破振動信號的速度-時間曲線Fig.8 The velocity vs time curve of blast vibration signal

    從圖9可以清楚地看出,微差爆破振動信號的db8小波、模式自適應(yīng)小波變換模值以及模式自適應(yīng)小波時能密度曲線中均出現(xiàn)了6個奇異點,位置分別為0.044 8,0.082 4,0.208 4,0.245 6,0.306 0,0.344 0和0.042 8,0.082 4,0.198 8,0.248 8,0.296 4,0.346 0以及0.042 8,0.082 4,0.198 4,0.250 0,0.296 0, 0.346 0 s,段間微差延期時間間隔分別為37.6, 126.0,37.2,60.4,38.0和39.6,116.4,50.0,47.6, 49.6以及39.6,116.0,51.6,46.0,50.0 ms。分別將db8小波法、模式自適應(yīng)小波法以及模式自適應(yīng)小波時能密度法得到的實際微差延期時間間隔與雷管的設(shè)計延期時間間隔進(jìn)行比較(表2)。

    圖9 微差爆破振動信號的db8小波、模式自適應(yīng)小波變換模值以及模式自適應(yīng)小波時能密度曲線Fig.9 db8 wavelet,pattern adapted wavelet transform modulus and pattern adapted wavelet time-energy density curve of measured millisecond blast vibration signal

    表2 采用3種方法對微差爆破振動信號的識別效果比較Table 2 Comparison of identification effect for millisecond blast vibration signals by three different methods

    從表2可以看出,采用3種不同小波分析法得到的2~4段次、5~6段次、6~7段次雷管起爆的實際微差延期時間間隔均在設(shè)計微差延期時間間隔范圍內(nèi),屬正常起爆;而4~5段次、7~9段次雷管起爆的實際微差延期時間間隔均超出了設(shè)計微差延期時間間隔范圍,屬非正常起爆。由此表明4~5段次、7~9段次雷管進(jìn)行微差爆破時需謹(jǐn)慎使用。雖然3種不同小波分析法對同一微差爆破振動信號進(jìn)行微差延期時間間隔的識別,效果均較理想,但通過圖9的比較不難看出,模式自適應(yīng)小波法的分辨率要比db8小波法的高,表明在微差爆破振動信號的奇異性檢測方面,基于爆破振動子信號的模式自適應(yīng)小波法比Matlab工具箱已有小波法的效果好,具有較好的突出微差爆破振動信號奇異點的位置,從而驗證了將爆破振動子信號作為爆破振動信號分析中的小波基是可行的,解決了適合爆破振動信號特征的小波基構(gòu)造、添加及其實現(xiàn)等問題,為分析微差爆破振動信號時選擇以及優(yōu)化小波基提出了一種新的方法。同時,模式自適應(yīng)小波時能密度法比模式自適應(yīng)小波法分辨率更高,壓制了隨機(jī)噪聲的干擾,突出了爆破振動信號奇異點位置,提高了爆破振動信號處理的效率與質(zhì)量。

    5 結(jié) 論

    (1)在模式自適應(yīng)連續(xù)小波原理的基礎(chǔ)上,構(gòu)造了從實測微差爆破振動信號中分離出子信號的爆破振動信號小波基,然后將該爆破振動信號小波基添加到Matlab小波分析工具箱中,從而解決了適合爆破振動信號特征的小波基構(gòu)造與添加及實現(xiàn)等問題,并提出了一種新的基于爆破振動信號小波分析的模式自適應(yīng)小波時能密度法,并且該方法能夠自適應(yīng)產(chǎn)生“基”,即不需要任何先驗基函數(shù)。

    (2)通過微差延期時間間隔下的多段微差爆破振動信號驗證了模式自適應(yīng)小波時能密度法在微差爆破振動信號分析中的可行性,然后在工程實例中比較了db8小波法、模式自適應(yīng)小波法以及模式自適應(yīng)小波時能密度法在微差起爆延時間隔的實際應(yīng)用效果,結(jié)果表明基于爆破振動信號特征的模式自適應(yīng)小波法、模式自適應(yīng)小波時能密度法在確定爆破振動信號微差延時方面比db8小波法效果更理想,且模式自適應(yīng)小波時能密度法比模式自適應(yīng)小波法分辨率更高,更好的突出了爆破振動信號奇異點位置,提高了爆破振動信號處理的效率與質(zhì)量,從而為小波理論在爆破工程實際中的應(yīng)用研究奠定了基礎(chǔ)。

    [1] 韓 博,馬芹永.煤礦巖巷毫秒延期爆破振動測試與控制技術(shù)研究[J].煤炭學(xué)報,2013,38(2):209-214.

    Han Bo,Ma Qinyong.Research on testing and controlling techniques of millisecond delay blasting vibration in coal mine roadway[J].Journal of China Coal Society,2013,38(2):209-214.

    [2] 凌同華,李夕兵.基于小波變換的時-能分布確定微差爆破的實際延遲時間[J].巖石力學(xué)與工程學(xué)報,2004,23(13):2266-2270.

    Ling Tonghua,Li Xibing.Time-energy analysis based on wavelet transform for identifying real delay time in millisecond blasting[J].Chinese Journal of Rock Mechanics and Engineering,2004,23 (13):2266-2270.

    [3] 凌同華,廖艷程,張 勝.沖擊荷載下巖石聲發(fā)射信號能量特征的小波包分析[J].振動與沖擊,2010,29(10):127-130.

    Ling Tonghua,Liao Yancheng,Zhang Sheng.Application of wavelet packet method in frequency band energy distribution of rock acoustic emission signals under impact loading[J].Journal of Vibration and Shock,2010,29(10):127-130.

    [4] 趙明生,梁開水,余德運(yùn),等.段數(shù)對爆破振動信號的時頻特性影響分析[J].煤炭學(xué)報,2012,37(1):55-61.

    Zhao Mingsheng,Liang Kaishui,Yu Deyun,et al.Effect of segments on time frequency characteristics of blasting vibration signals[J].Journal of China Coal Society,2012,37(1):55-61.

    [5] 李夕兵,凌同華,張義平.爆破震動信號分析理論與技術(shù)[M].北京:科學(xué)出版社,2009:66-81.

    [6] 李 臻.自適應(yīng)提升小波的構(gòu)造及其在齒輪箱早期故障診斷中的應(yīng)用[J].煤炭學(xué)報,2010,35(S):228-231.

    Li Zhen.The construction of adaptive lifting wavelet and its application in incipient fault diagnosis of gearbox[J].Journal of China Coal Society,2010,35(S):228-231.

    [7] 凌同華.爆破震動效應(yīng)及其災(zāi)害的主動控制[D].長沙:中南大學(xué),2004:34-52.

    [8] 劉 高.深部大規(guī)模爆破振動信號分析與災(zāi)害控制研究[D].長沙:中南大學(xué),2011:44-57.

    [9] 耿艷峰,馮叔初.小波構(gòu)造綜述[J].石油大學(xué)學(xué)報(自然科學(xué)版),2004,28(1):127-131.

    Geng Yanfeng,Feng Shuchu.Overview of wavelet construction[J]. Journal of the University of Petroleum,China(Edition of Natural Science),2004,28(1):127-131.

    [10] 宋光明.爆破振動小波包分析理論與應(yīng)用研究[D].長沙:中南大學(xué),2001:18-30.

    [11] 李春庚.自適應(yīng)小波的構(gòu)造及其在信號處理中的應(yīng)用[D].大連:大連海事大學(xué),2007:1-16.

    [12] 陳倩倩.基于巖土工程信號特征的小波基構(gòu)造及其算法與實現(xiàn)研究[D].長沙:長沙理工大學(xué),2013:1-9.

    [13] 丁愛玲,石光明,張 寧,等.波形匹配小波設(shè)計及其在信號數(shù)據(jù)壓縮中的應(yīng)用[J].電子與信息學(xué)報,2007,29(4):804-807.

    Ding Ailing,Shi Guangming,Zhang Ning,et al.Signal compression and design of wavelet based on waveform matching[J].Journal of Electronics&Information Technology,2007,29(4):804-807.

    [14] 程正興,楊守志,馮曉霞.小波分析的理論、算法、進(jìn)展和應(yīng)用[M].北京:國防工業(yè)出版社,2007:1-11.

    [15] 張銳戈,譚永紅.基于最優(yōu)Morlet小波和隱馬爾可夫模型的軸承故障診斷[J].振動與沖擊,2012,31(12):5-8,27.

    Zhang Ruige,Tan Yonghong.Fault diagnosis of rolling element bearings based on optimal morlet wavelet and hidden markov model [J].Journal of Vibration and Shock,2012,31(12):5-8,27.

    [16] 栗茂林,梁 霖,王孫安,等.基于連續(xù)小波系數(shù)非線性流形學(xué)習(xí)的沖擊特征提取方法[J].振動與沖擊,2012,31(1):106-111,126.

    Li Maolin,Liang Lin,Wang Sun’an,et al.Mechanical impact feature extraction method based on nonlinear manifold learning of continuous wavelet coefficients[J].Journal of Vibration and Shock, 2012,31(1):106-111,126.

    [17] 張煜東,朱 慶,陳賢卿,等.基于模式的自制小波在EBPSK信號檢測中的應(yīng)用[J].東南大學(xué)學(xué)報(自然科學(xué)版),2011, 41(4):691-694.

    Zhang Yudong,Zhu Qing,Chen Xianqing,et al.Pattern-based custom wavelet used for EBPSK signal detection[J].Journal of Southeast University(Natural Science Edition),2011,41(4):691-694.

    [18] Mesa H.Adapted wavelets for pattern detection[A].Progress in Pattern Recognition,Image Analysis and Applications[C].Springer Berlin Heidelberg,2005:933-944.

    [19] Chapa J O,Rao R M.Algorithms for designing wavelets to match a specified signal[J].IEEE Transactions on Signal Processing, 2000,48(12):3395-3406.

    [20] 胡昌華,李國華,周 濤.基于MATLAB 7.x的系統(tǒng)分析與設(shè)計—小波分析(第三版)[M].西安:西安電子科技大學(xué)出版社,2008:236-245.

    [21] 嚴(yán) 鵬,盧文波,羅 憶,等.基于小波變換時-能密度分析的爆破開挖過程中地應(yīng)力動態(tài)卸載振動到達(dá)時刻識別[J].巖石力學(xué)與工程學(xué)報,2009,28(S1):2836-2844.

    Yan Peng,Lu Wenbo,Luo Yi,et al.Identification of arriving time of vibration induced by geostress dynamic unloading during blasting-excavation employing method of time-energy analysis based on wavelet transform[J].Chinese Journal of Rock Mechanics and Engineering,2009,28(S1):2836-2844.

    Pattern adapted wavelet time-energy density method and its application in millisecond blast vibration signal analysis

    ZHANG Sheng1,2,LING Tong-hua1,LIU Hao-ran1,CAO Feng1

    (1.School of Civil Engineering and Architecture,Changsha University of Science&Technology,Changsha 410114,China;2.Key Laboratory for Safety Control of Bridge Engineering,Ministry of Education and Hunan Province,Changsha University of Science&Technology,Changsha 410114,China)

    How to solve wavelet basis construction,addition,realization and other issues that suitable for characterizing blast vibration signal is a challenge for many researchers.To determine real millisecond delay time of millisecond blast vibration signal accurately,a pattern adapted wavelet basis separated from measured millisecond blast vibration signal was constructed based on the principle of pattern adapted wavelet.The pattern adapted wavelet basis was added to Wavelet Toolbox in Matlab,and a new pattern adapted wavelet time-energy density method was proposed based on blast vibration signal analysis.The feasibility of the method was verified with the pattern adapted wavelet time-energy density method of millisecond blast vibration signal for experimental design.The method can also be used to analyze the measured millisecond blast vibration signal to determine real millisecond delay time.Compared with db8 wavelet method and with pattern adapted wavelet method,the results indicate that the pattern adapted wavelet time-energy density method has a better application effect.Thereby,the method provides a reference for blast vibration signal analysisin practical engineering application.

    blast vibration;millisecond time;pattern adapted wavelet time-energy density method;wavelet basis construction;wavelet basis addition

    TD235

    A

    0253-9993(2014)10-2007-07

    2014-03-14 責(zé)任編輯:常 琛

    國家自然科學(xué)基金資助項目(51278071);湖南省研究生科研創(chuàng)新資助項目(CX2013B368);長沙理工大學(xué)橋梁工程安全控制省部共建教育部重點實驗室開放基金資助項目(13KB03)

    張 勝(1984—),男,湖南益陽人,博士研究生。E-mail:zhangsheng0403311@163.com。通訊作者:凌同華(1968—),男,湖南雙峰人,教授,博士生導(dǎo)師,博士。E-mail:lingtonghua@163.com

    張 勝,凌同華,劉浩然,等.模式自適應(yīng)小波時能密度法及其在微差爆破振動信號分析中的應(yīng)用[J].煤炭學(xué)報,2014,39(10): 2007-2013.

    10.13225/j.cnki.jccs.2014.0325

    Zhang Sheng,Ling Tonghua,Liu Haoran,et al.Pattern adapted wavelet time-energy density method and its application in millisecond blast vibration signal analysis[J].Journal of China Coal Society,2014,39(10):2007-2013.doi:10.13225/j.cnki.jccs.2014.0325

    猜你喜歡
    小波基微差小波
    車軸軸徑微差測量儀測量不確定度評定
    構(gòu)造Daubechies小波的一些注記
    市政基坑砼支撐梁微差延期爆破拆除技術(shù)
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    利用小波變換分析電能質(zhì)量擾動問題中的電壓驟升影響
    小波閾值圖像去噪中小波基選擇
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    小波非參數(shù)回歸分析方法的實現(xiàn)及比較研究*
    基于FPGA小波變換核的設(shè)計
    電測與儀表(2014年8期)2014-04-04 09:19:38
    金堆城南露天數(shù)碼雷管減震爆破試驗研究
    一个人观看的视频www高清免费观看| 国产亚洲5aaaaa淫片| 久久精品国产亚洲av天美| 三级毛片av免费| 亚洲精品中文字幕在线视频 | 99久久中文字幕三级久久日本| 国产成人精品婷婷| 国产黄色小视频在线观看| 久久久久网色| av线在线观看网站| 国产在线男女| 国产乱人偷精品视频| 国产一区亚洲一区在线观看| 精品久久久久久久人妻蜜臀av| 午夜福利视频精品| 亚洲精品影视一区二区三区av| 欧美性感艳星| 国产高清国产精品国产三级 | 国产高清有码在线观看视频| 色吧在线观看| 国产精品久久久久久精品电影| 亚洲美女视频黄频| 日本色播在线视频| 2021天堂中文幕一二区在线观| 亚洲精品久久久久久婷婷小说| 久久久久久久亚洲中文字幕| 成年女人在线观看亚洲视频 | 久久精品国产亚洲av天美| 国产一级毛片在线| 国产精品伦人一区二区| 人妻系列 视频| 成人av在线播放网站| 亚洲精品日本国产第一区| 日本爱情动作片www.在线观看| 国产精品国产三级国产av玫瑰| 26uuu在线亚洲综合色| 国产精品福利在线免费观看| 男女啪啪激烈高潮av片| 午夜福利高清视频| 精品一区二区三区人妻视频| 黄色日韩在线| 人人妻人人澡人人爽人人夜夜 | 又大又黄又爽视频免费| 久久国产乱子免费精品| 精品久久久久久久久亚洲| 床上黄色一级片| 欧美激情在线99| 秋霞在线观看毛片| 午夜视频国产福利| 久久99精品国语久久久| 最近视频中文字幕2019在线8| 久久午夜福利片| 99热这里只有是精品在线观看| av黄色大香蕉| 在线播放无遮挡| 婷婷色麻豆天堂久久| 欧美bdsm另类| 热99在线观看视频| 亚洲精品aⅴ在线观看| 精品欧美国产一区二区三| 国产精品熟女久久久久浪| 欧美xxxx性猛交bbbb| 国产69精品久久久久777片| 麻豆成人午夜福利视频| 国产69精品久久久久777片| 99久国产av精品| 久久久久久久大尺度免费视频| 国产在线一区二区三区精| 免费少妇av软件| 国产激情偷乱视频一区二区| 精品久久久噜噜| 啦啦啦韩国在线观看视频| 18+在线观看网站| 日韩欧美国产在线观看| 亚洲在线自拍视频| 成人一区二区视频在线观看| 少妇的逼水好多| 午夜免费观看性视频| 亚洲精品一二三| 91精品一卡2卡3卡4卡| 十八禁网站网址无遮挡 | av一本久久久久| 深爱激情五月婷婷| 亚洲精品国产成人久久av| 天美传媒精品一区二区| av在线播放精品| 性插视频无遮挡在线免费观看| 日日摸夜夜添夜夜爱| 丰满少妇做爰视频| 在线观看免费高清a一片| 禁无遮挡网站| 美女cb高潮喷水在线观看| 简卡轻食公司| 最近的中文字幕免费完整| 国内精品美女久久久久久| 免费av观看视频| 亚洲性久久影院| 亚洲色图av天堂| 免费观看在线日韩| 午夜福利在线在线| av国产免费在线观看| 免费电影在线观看免费观看| 麻豆国产97在线/欧美| 国产成人免费观看mmmm| www.色视频.com| 中文字幕av成人在线电影| av天堂中文字幕网| 我的老师免费观看完整版| 男人舔奶头视频| 国产午夜精品久久久久久一区二区三区| 欧美日韩精品成人综合77777| 欧美bdsm另类| 日韩欧美三级三区| 最新中文字幕久久久久| 日韩av不卡免费在线播放| 精品熟女少妇av免费看| 婷婷色麻豆天堂久久| 乱人视频在线观看| 国产在线一区二区三区精| 免费不卡的大黄色大毛片视频在线观看 | 日本免费a在线| 天天一区二区日本电影三级| 免费无遮挡裸体视频| h日本视频在线播放| h日本视频在线播放| 久久久亚洲精品成人影院| 日日摸夜夜添夜夜添av毛片| 日日摸夜夜添夜夜添av毛片| 熟女电影av网| 97人妻精品一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 久久久久精品久久久久真实原创| 国产永久视频网站| 亚洲自拍偷在线| 日本与韩国留学比较| 97人妻精品一区二区三区麻豆| 美女主播在线视频| av在线观看视频网站免费| 日韩不卡一区二区三区视频在线| 久久久久久久久大av| 最近最新中文字幕免费大全7| 久久精品夜夜夜夜夜久久蜜豆| 国产 一区精品| 99热网站在线观看| 丰满乱子伦码专区| 久久精品久久精品一区二区三区| 最近手机中文字幕大全| 秋霞在线观看毛片| 精品亚洲乱码少妇综合久久| 精品久久久久久久人妻蜜臀av| 国产精品.久久久| 国产有黄有色有爽视频| av福利片在线观看| 中文字幕av成人在线电影| 少妇被粗大猛烈的视频| 美女黄网站色视频| 天美传媒精品一区二区| 人妻系列 视频| 欧美极品一区二区三区四区| 国产高清不卡午夜福利| 国产成人精品婷婷| 一级黄片播放器| 天堂av国产一区二区熟女人妻| 永久网站在线| 搡老乐熟女国产| 久久久久性生活片| 99re6热这里在线精品视频| 久久国产乱子免费精品| 丝袜美腿在线中文| 国产免费又黄又爽又色| 欧美日韩一区二区视频在线观看视频在线 | 又黄又爽又刺激的免费视频.| av天堂中文字幕网| 九色成人免费人妻av| 人妻少妇偷人精品九色| 亚洲精品一二三| 只有这里有精品99| 美女内射精品一级片tv| av黄色大香蕉| 男女边摸边吃奶| 在线观看美女被高潮喷水网站| 免费在线观看成人毛片| 人妻制服诱惑在线中文字幕| 熟女电影av网| 成年女人看的毛片在线观看| 性色avwww在线观看| 97超视频在线观看视频| 18禁在线播放成人免费| 欧美另类一区| 国产黄片视频在线免费观看| 午夜免费激情av| 晚上一个人看的免费电影| 午夜久久久久精精品| 18+在线观看网站| av天堂中文字幕网| 国产免费一级a男人的天堂| 亚洲精品亚洲一区二区| 极品少妇高潮喷水抽搐| 91精品一卡2卡3卡4卡| 久久久午夜欧美精品| 国产精品国产三级专区第一集| 黄片wwwwww| 草草在线视频免费看| 只有这里有精品99| av免费观看日本| 国产精品一区二区三区四区久久| 国精品久久久久久国模美| 国产精品一区www在线观看| 内地一区二区视频在线| 欧美日本视频| av.在线天堂| 日本午夜av视频| 有码 亚洲区| 嫩草影院入口| 超碰97精品在线观看| 尤物成人国产欧美一区二区三区| 精品久久久久久久人妻蜜臀av| 2021天堂中文幕一二区在线观| 2018国产大陆天天弄谢| 久久久久久久久久黄片| 人人妻人人澡人人爽人人夜夜 | 国产欧美日韩精品一区二区| 国产成人精品一,二区| 精品国内亚洲2022精品成人| 人妻一区二区av| 欧美区成人在线视频| 中文欧美无线码| 高清av免费在线| 一区二区三区乱码不卡18| 99热这里只有是精品在线观看| 国产精品久久久久久av不卡| 18禁裸乳无遮挡免费网站照片| 亚洲国产精品国产精品| 国产激情偷乱视频一区二区| 免费观看在线日韩| 国内精品美女久久久久久| 99热网站在线观看| 午夜老司机福利剧场| 午夜久久久久精精品| 三级男女做爰猛烈吃奶摸视频| 亚洲美女视频黄频| 国产精品久久视频播放| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久久成人| 男女下面进入的视频免费午夜| 男女边吃奶边做爰视频| 亚洲欧美清纯卡通| 成人二区视频| 国产在线男女| 亚洲欧美精品自产自拍| 中文欧美无线码| 久久综合国产亚洲精品| 亚洲国产精品成人久久小说| 国产免费又黄又爽又色| 大陆偷拍与自拍| 精品久久久久久电影网| 卡戴珊不雅视频在线播放| 人妻系列 视频| 精品国产一区二区三区久久久樱花 | 国产又色又爽无遮挡免| 亚洲精品视频女| 少妇的逼好多水| 亚洲va在线va天堂va国产| 国产色爽女视频免费观看| 肉色欧美久久久久久久蜜桃 | 亚洲图色成人| 国产亚洲一区二区精品| 亚洲美女搞黄在线观看| 午夜视频国产福利| 亚洲av二区三区四区| 日韩一本色道免费dvd| 免费播放大片免费观看视频在线观看| 我的女老师完整版在线观看| 热99在线观看视频| 天堂√8在线中文| 99久久精品国产国产毛片| 国产在线一区二区三区精| 不卡视频在线观看欧美| 国产永久视频网站| 又粗又硬又长又爽又黄的视频| 国产精品一及| 校园人妻丝袜中文字幕| 建设人人有责人人尽责人人享有的 | 插阴视频在线观看视频| 免费黄色在线免费观看| 日韩国内少妇激情av| 亚洲精品乱久久久久久| 少妇人妻精品综合一区二区| 免费少妇av软件| 国产一区二区亚洲精品在线观看| 一级毛片电影观看| 丰满乱子伦码专区| 久久久久网色| 午夜精品一区二区三区免费看| 一级毛片黄色毛片免费观看视频| 久久99热这里只有精品18| 免费看美女性在线毛片视频| 亚洲天堂国产精品一区在线| 色网站视频免费| 乱人视频在线观看| 日韩欧美国产在线观看| 九九在线视频观看精品| 舔av片在线| 欧美极品一区二区三区四区| 91午夜精品亚洲一区二区三区| 国产精品久久久久久av不卡| 在线播放无遮挡| 免费不卡的大黄色大毛片视频在线观看 | 一区二区三区乱码不卡18| 色5月婷婷丁香| 精品久久国产蜜桃| 精品一区二区三卡| 欧美精品国产亚洲| 久热久热在线精品观看| 精品人妻视频免费看| 欧美xxxx性猛交bbbb| 日本-黄色视频高清免费观看| 亚洲国产精品成人综合色| 免费人成在线观看视频色| 久久久成人免费电影| 99热6这里只有精品| 亚洲色图av天堂| 亚洲熟妇中文字幕五十中出| 99热全是精品| 91午夜精品亚洲一区二区三区| 国产精品久久久久久av不卡| 欧美变态另类bdsm刘玥| 97超视频在线观看视频| 高清毛片免费看| 久久精品人妻少妇| 国语对白做爰xxxⅹ性视频网站| 久久亚洲国产成人精品v| xxx大片免费视频| 天堂俺去俺来也www色官网 | 在线观看av片永久免费下载| av卡一久久| 免费av观看视频| 国产成人91sexporn| 九九爱精品视频在线观看| 神马国产精品三级电影在线观看| 色播亚洲综合网| 亚洲欧美日韩东京热| 日韩av免费高清视频| 欧美 日韩 精品 国产| 91久久精品电影网| 高清欧美精品videossex| 搞女人的毛片| 嫩草影院新地址| 日韩电影二区| 久久99热这里只有精品18| 日本欧美国产在线视频| 最近中文字幕2019免费版| 国产精品爽爽va在线观看网站| 久久这里有精品视频免费| 国产精品久久久久久精品电影| 国产69精品久久久久777片| 国产极品天堂在线| 一区二区三区乱码不卡18| 欧美人与善性xxx| 天堂√8在线中文| 国产一级毛片在线| av在线观看视频网站免费| 好男人在线观看高清免费视频| kizo精华| 最近视频中文字幕2019在线8| 成人亚洲精品av一区二区| 久久久亚洲精品成人影院| 天美传媒精品一区二区| 97超碰精品成人国产| xxx大片免费视频| 亚洲在线自拍视频| 亚洲精品成人久久久久久| 久久久色成人| 亚洲av二区三区四区| 久久精品综合一区二区三区| 中文在线观看免费www的网站| av国产免费在线观看| av女优亚洲男人天堂| 91精品一卡2卡3卡4卡| 国产高清国产精品国产三级 | 大片免费播放器 马上看| 在线播放无遮挡| 午夜激情久久久久久久| 午夜激情福利司机影院| 汤姆久久久久久久影院中文字幕 | 久久精品熟女亚洲av麻豆精品 | 人妻少妇偷人精品九色| 亚洲国产高清在线一区二区三| 卡戴珊不雅视频在线播放| 国产黄片美女视频| 免费观看性生交大片5| 国产又色又爽无遮挡免| 高清午夜精品一区二区三区| 搡老乐熟女国产| 中文字幕久久专区| 又大又黄又爽视频免费| 老师上课跳d突然被开到最大视频| 婷婷色综合大香蕉| 亚洲精品国产av蜜桃| 国产午夜精品久久久久久一区二区三区| 在线 av 中文字幕| 午夜福利高清视频| 午夜爱爱视频在线播放| 夜夜看夜夜爽夜夜摸| 久久国内精品自在自线图片| .国产精品久久| 欧美一区二区亚洲| 久久久久久久久久成人| 18禁在线播放成人免费| 熟妇人妻久久中文字幕3abv| 超碰av人人做人人爽久久| 免费黄色在线免费观看| 亚洲欧美一区二区三区黑人 | 丝袜美腿在线中文| 在线a可以看的网站| 亚洲aⅴ乱码一区二区在线播放| 日日摸夜夜添夜夜添av毛片| 高清av免费在线| 国产午夜精品一二区理论片| 一级a做视频免费观看| 丝袜喷水一区| 最近最新中文字幕免费大全7| 色尼玛亚洲综合影院| eeuss影院久久| 亚洲av电影不卡..在线观看| 日韩伦理黄色片| 欧美 日韩 精品 国产| 成人国产麻豆网| 欧美成人精品欧美一级黄| 国产一级毛片七仙女欲春2| 最近手机中文字幕大全| 国产伦在线观看视频一区| 亚洲18禁久久av| 亚洲成人一二三区av| 三级男女做爰猛烈吃奶摸视频| 亚州av有码| 少妇人妻精品综合一区二区| 97超碰精品成人国产| 美女大奶头视频| 简卡轻食公司| 久久精品国产亚洲av涩爱| 日日撸夜夜添| 国产精品美女特级片免费视频播放器| 十八禁网站网址无遮挡 | 国产免费一级a男人的天堂| 欧美日韩视频高清一区二区三区二| 久久这里有精品视频免费| 3wmmmm亚洲av在线观看| 久久久久精品性色| 不卡视频在线观看欧美| 777米奇影视久久| 免费黄频网站在线观看国产| 成人国产麻豆网| 麻豆精品久久久久久蜜桃| 少妇人妻一区二区三区视频| 亚洲美女搞黄在线观看| 日日撸夜夜添| 日韩一区二区视频免费看| 97超视频在线观看视频| 国产精品人妻久久久久久| 嫩草影院新地址| 国产老妇女一区| 成人av在线播放网站| 国产黄频视频在线观看| 亚洲不卡免费看| 亚洲精品乱久久久久久| 天堂av国产一区二区熟女人妻| 久久久久久久午夜电影| 中文字幕人妻熟人妻熟丝袜美| 色综合亚洲欧美另类图片| 日日啪夜夜爽| 国产伦精品一区二区三区视频9| 免费少妇av软件| 亚洲av成人av| 日韩成人伦理影院| 美女xxoo啪啪120秒动态图| kizo精华| 青春草亚洲视频在线观看| 欧美高清成人免费视频www| 亚洲精品色激情综合| 国产午夜精品一二区理论片| 青春草亚洲视频在线观看| 国产精品人妻久久久久久| 午夜亚洲福利在线播放| 国产精品一区二区在线观看99 | 久久久久网色| 在线天堂最新版资源| 亚洲精品国产av成人精品| 夫妻性生交免费视频一级片| 69人妻影院| 欧美+日韩+精品| 成人午夜高清在线视频| 国产高清三级在线| 嫩草影院精品99| 色5月婷婷丁香| 最近视频中文字幕2019在线8| 亚洲无线观看免费| 亚洲婷婷狠狠爱综合网| 国产熟女欧美一区二区| 五月天丁香电影| 肉色欧美久久久久久久蜜桃 | 美女cb高潮喷水在线观看| 国产乱人视频| 能在线免费观看的黄片| 欧美日韩精品成人综合77777| 欧美日韩亚洲高清精品| 国产一区二区在线观看日韩| 国产在视频线在精品| 国产一区亚洲一区在线观看| 亚洲人成网站在线观看播放| 成人漫画全彩无遮挡| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| 久久久欧美国产精品| 日日摸夜夜添夜夜爱| 欧美日韩综合久久久久久| 好男人在线观看高清免费视频| 精品人妻偷拍中文字幕| 免费黄色在线免费观看| 亚洲怡红院男人天堂| 丝袜喷水一区| 中国国产av一级| 不卡视频在线观看欧美| 91久久精品国产一区二区三区| 久久久久久久久中文| 啦啦啦中文免费视频观看日本| 国产精品一区二区在线观看99 | 亚洲天堂国产精品一区在线| or卡值多少钱| 国产亚洲精品av在线| 欧美日韩国产mv在线观看视频 | 18+在线观看网站| 成人漫画全彩无遮挡| 特级一级黄色大片| 日本欧美国产在线视频| 1000部很黄的大片| 中文字幕制服av| 亚洲国产精品专区欧美| 中文字幕av在线有码专区| 深爱激情五月婷婷| 国内精品一区二区在线观看| 人妻少妇偷人精品九色| 亚洲电影在线观看av| 又爽又黄无遮挡网站| 亚洲精品色激情综合| av黄色大香蕉| 亚洲av成人av| 久久久久免费精品人妻一区二区| 18禁在线播放成人免费| 一个人看的www免费观看视频| 一级毛片我不卡| 纵有疾风起免费观看全集完整版 | 伦理电影大哥的女人| 国产精品久久久久久av不卡| 亚洲av中文av极速乱| 美女高潮的动态| 亚洲国产欧美在线一区| 国产成人午夜福利电影在线观看| 国产成人福利小说| 成年版毛片免费区| 大香蕉97超碰在线| 日本色播在线视频| 亚洲天堂国产精品一区在线| 五月玫瑰六月丁香| 成人高潮视频无遮挡免费网站| 精品久久久久久久末码| 人人妻人人澡欧美一区二区| 一级黄片播放器| av免费在线看不卡| 成人欧美大片| 一个人免费在线观看电影| 成人毛片60女人毛片免费| 久久草成人影院| 亚洲美女视频黄频| 熟妇人妻不卡中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 精品一区二区三卡| 少妇被粗大猛烈的视频| 久久精品夜色国产| 国产黄频视频在线观看| 99热全是精品| 国产黄色免费在线视频| 成年免费大片在线观看| 日韩精品青青久久久久久| 亚洲av免费高清在线观看| 成人亚洲欧美一区二区av| 久久久精品94久久精品| 国产精品不卡视频一区二区| 99久久精品热视频| 成年av动漫网址| 黑人高潮一二区| 国产 一区 欧美 日韩| 亚洲av成人精品一二三区| 成人av在线播放网站| 国产精品一区二区性色av| 亚洲va在线va天堂va国产| 国产精品1区2区在线观看.| 男女视频在线观看网站免费| 成年版毛片免费区| 日本欧美国产在线视频| 免费观看在线日韩| 欧美日韩国产mv在线观看视频 | 亚洲精品乱码久久久久久按摩| av在线播放精品| 国产午夜精品一二区理论片| 九色成人免费人妻av| 人妻一区二区av| 亚洲av电影在线观看一区二区三区 | 日本黄大片高清| 免费观看性生交大片5| 久久精品国产亚洲av涩爱| 一级毛片久久久久久久久女| 日日摸夜夜添夜夜添av毛片| 亚洲精品乱久久久久久| 亚洲人成网站在线观看播放|