• <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ù)碼雷管減震爆破試驗研究
    成年女人毛片免费观看观看9 | 午夜福利,免费看| 伦理电影免费视频| 免费观看性生交大片5| 亚洲欧美一区二区三区久久| 啦啦啦在线免费观看视频4| a级毛片黄视频| 国产精品蜜桃在线观看| 亚洲av在线观看美女高潮| 亚洲欧美成人精品一区二区| 国产精品国产三级专区第一集| 精品午夜福利在线看| 久久97久久精品| 国产 精品1| 新久久久久国产一级毛片| 成年人免费黄色播放视频| 男女免费视频国产| 可以免费在线观看a视频的电影网站 | 人人妻人人添人人爽欧美一区卜| 国产日韩欧美视频二区| 免费看不卡的av| 午夜福利免费观看在线| 别揉我奶头~嗯~啊~动态视频 | 中国三级夫妇交换| 热re99久久国产66热| 夜夜骑夜夜射夜夜干| 老司机亚洲免费影院| 久久久久久免费高清国产稀缺| 18禁观看日本| 色94色欧美一区二区| 国产一区二区三区av在线| 欧美日韩亚洲国产一区二区在线观看 | 久久精品熟女亚洲av麻豆精品| 国产精品一区二区在线观看99| 这个男人来自地球电影免费观看 | 亚洲,欧美精品.| 满18在线观看网站| 男女高潮啪啪啪动态图| 亚洲国产精品成人久久小说| 一区二区日韩欧美中文字幕| 午夜福利在线免费观看网站| 亚洲精品美女久久久久99蜜臀 | 国产 精品1| 中文字幕av电影在线播放| av天堂久久9| 成人影院久久| 久久ye,这里只有精品| netflix在线观看网站| 一本大道久久a久久精品| 国产老妇伦熟女老妇高清| 热99国产精品久久久久久7| 欧美人与性动交α欧美软件| 中文字幕精品免费在线观看视频| 亚洲综合精品二区| 免费日韩欧美在线观看| 侵犯人妻中文字幕一二三四区| 久久精品久久久久久久性| 久久精品aⅴ一区二区三区四区| 精品国产乱码久久久久久男人| 日本一区二区免费在线视频| netflix在线观看网站| 熟妇人妻不卡中文字幕| 国产精品一区二区在线观看99| 最黄视频免费看| 9色porny在线观看| 悠悠久久av| 日韩,欧美,国产一区二区三区| 两性夫妻黄色片| 亚洲精品国产色婷婷电影| 欧美精品亚洲一区二区| 亚洲第一青青草原| 日韩精品有码人妻一区| 一级毛片我不卡| 精品亚洲乱码少妇综合久久| 国产免费福利视频在线观看| 精品一区在线观看国产| 国产一级毛片在线| 成人亚洲欧美一区二区av| 日本wwww免费看| 精品少妇内射三级| 国产乱人偷精品视频| 大陆偷拍与自拍| 视频区图区小说| 国产免费视频播放在线视频| 无遮挡黄片免费观看| 一级毛片黄色毛片免费观看视频| av又黄又爽大尺度在线免费看| 不卡av一区二区三区| videos熟女内射| 老司机影院毛片| 午夜老司机福利片| 五月天丁香电影| 男人添女人高潮全过程视频| 这个男人来自地球电影免费观看 | 午夜福利免费观看在线| 国产一区亚洲一区在线观看| 最近2019中文字幕mv第一页| 久久精品国产亚洲av涩爱| 国产黄色免费在线视频| 看免费av毛片| 女人高潮潮喷娇喘18禁视频| 男女下面插进去视频免费观看| 国产成人啪精品午夜网站| 色婷婷久久久亚洲欧美| 国产伦理片在线播放av一区| 日本一区二区免费在线视频| 老司机亚洲免费影院| 青春草视频在线免费观看| 51午夜福利影视在线观看| videosex国产| av电影中文网址| 亚洲成人手机| 少妇猛男粗大的猛烈进出视频| 久久99热这里只频精品6学生| 精品国产超薄肉色丝袜足j| a级片在线免费高清观看视频| 久久精品亚洲av国产电影网| 好男人视频免费观看在线| 97在线人人人人妻| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区| 国产日韩欧美视频二区| 国产精品一区二区在线观看99| 少妇人妻久久综合中文| 欧美日韩精品网址| 久久精品国产亚洲av涩爱| 亚洲欧洲国产日韩| 老司机影院成人| 久热这里只有精品99| av网站免费在线观看视频| 日本猛色少妇xxxxx猛交久久| 99热全是精品| www.av在线官网国产| 99热国产这里只有精品6| 亚洲国产欧美网| 精品亚洲乱码少妇综合久久| 大话2 男鬼变身卡| 久久久久久久久免费视频了| av在线app专区| 国产成人精品福利久久| 国产免费福利视频在线观看| 狠狠婷婷综合久久久久久88av| 又黄又粗又硬又大视频| 精品久久久精品久久久| 丝袜脚勾引网站| 国产亚洲一区二区精品| 国产亚洲av片在线观看秒播厂| 国产成人一区二区在线| 永久免费av网站大全| 亚洲美女视频黄频| 久久热在线av| 精品国产一区二区久久| 黑人猛操日本美女一级片| av电影中文网址| 日韩制服丝袜自拍偷拍| 亚洲男人天堂网一区| av网站免费在线观看视频| 久久久久国产一级毛片高清牌| 国产1区2区3区精品| 青春草国产在线视频| 国产精品香港三级国产av潘金莲 | 久久久久精品人妻al黑| 这个男人来自地球电影免费观看 | 国产淫语在线视频| 中文字幕人妻丝袜制服| 丝袜人妻中文字幕| 国产精品 欧美亚洲| 美国免费a级毛片| 18禁动态无遮挡网站| 在线观看免费视频网站a站| 国产av精品麻豆| 老鸭窝网址在线观看| 男女边摸边吃奶| 十分钟在线观看高清视频www| 欧美久久黑人一区二区| 日本欧美国产在线视频| 国产爽快片一区二区三区| 亚洲视频免费观看视频| 精品免费久久久久久久清纯 | 最近最新中文字幕大全免费视频 | svipshipincom国产片| 热99国产精品久久久久久7| 久久久久久人妻| 亚洲成色77777| 少妇的丰满在线观看| 亚洲三区欧美一区| 最近中文字幕2019免费版| 男女边吃奶边做爰视频| 亚洲精品一区蜜桃| 男人操女人黄网站| 久久久久精品性色| 人人妻人人澡人人爽人人夜夜| 成人三级做爰电影| 性色av一级| 精品国产一区二区三区久久久樱花| 国产亚洲午夜精品一区二区久久| 成年av动漫网址| 制服丝袜香蕉在线| 少妇猛男粗大的猛烈进出视频| 国产日韩欧美在线精品| av免费观看日本| 日韩欧美精品免费久久| 19禁男女啪啪无遮挡网站| 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区国产| 精品国产露脸久久av麻豆| 精品第一国产精品| 午夜影院在线不卡| 国产精品国产三级专区第一集| 日韩欧美精品免费久久| 人妻一区二区av| 亚洲精品乱久久久久久| 男人舔女人的私密视频| 超色免费av| 99久久99久久久精品蜜桃| 午夜福利一区二区在线看| 在线观看免费午夜福利视频| 天堂俺去俺来也www色官网| 免费高清在线观看日韩| 久久99热这里只频精品6学生| 国产色婷婷99| 日韩一卡2卡3卡4卡2021年| 久久久久国产精品人妻一区二区| 成人亚洲精品一区在线观看| 久久ye,这里只有精品| 亚洲av在线观看美女高潮| 亚洲激情五月婷婷啪啪| 亚洲精品日本国产第一区| 九草在线视频观看| 91老司机精品| 国产不卡av网站在线观看| 亚洲欧美成人综合另类久久久| 久久久久网色| 亚洲成色77777| 丰满迷人的少妇在线观看| 日本91视频免费播放| 免费人妻精品一区二区三区视频| 美女国产高潮福利片在线看| 99精国产麻豆久久婷婷| 啦啦啦在线免费观看视频4| 考比视频在线观看| av又黄又爽大尺度在线免费看| 晚上一个人看的免费电影| 在线看a的网站| 一区二区三区四区激情视频| 欧美老熟妇乱子伦牲交| 亚洲一级一片aⅴ在线观看| 国产乱来视频区| 亚洲中文av在线| 欧美中文综合在线视频| 亚洲美女搞黄在线观看| 久久99精品国语久久久| 国产老妇伦熟女老妇高清| 不卡视频在线观看欧美| 国产高清国产精品国产三级| 男女无遮挡免费网站观看| 日韩精品有码人妻一区| 国产精品一二三区在线看| 狠狠婷婷综合久久久久久88av| 啦啦啦中文免费视频观看日本| 欧美变态另类bdsm刘玥| av网站在线播放免费| 99久久人妻综合| 亚洲激情五月婷婷啪啪| 久久久久网色| 中文字幕人妻丝袜制服| 色播在线永久视频| 久久鲁丝午夜福利片| 亚洲综合色网址| 久久精品国产a三级三级三级| 欧美精品一区二区免费开放| 亚洲三区欧美一区| 久久免费观看电影| 久久久久网色| 满18在线观看网站| 热99国产精品久久久久久7| videosex国产| 卡戴珊不雅视频在线播放| 国产又色又爽无遮挡免| 日韩av免费高清视频| av线在线观看网站| 日本wwww免费看| 一区二区三区激情视频| 精品福利永久在线观看| 中文字幕精品免费在线观看视频| 国产欧美亚洲国产| 我的亚洲天堂| 亚洲第一青青草原| 成人18禁高潮啪啪吃奶动态图| 免费少妇av软件| 又大又黄又爽视频免费| 悠悠久久av| 亚洲一区二区三区欧美精品| 色婷婷久久久亚洲欧美| 在线看a的网站| 人妻人人澡人人爽人人| 国产一区二区激情短视频 | 熟女少妇亚洲综合色aaa.| 精品国产一区二区三区四区第35| 午夜老司机福利片| 天堂俺去俺来也www色官网| 男女免费视频国产| 在线 av 中文字幕| 精品免费久久久久久久清纯 | 国产又爽黄色视频| 国产精品亚洲av一区麻豆 | xxx大片免费视频| 韩国精品一区二区三区| 日本午夜av视频| 超碰成人久久| 天天添夜夜摸| 久久 成人 亚洲| 久久精品国产亚洲av涩爱| 黑人猛操日本美女一级片| 精品一区二区三卡| 亚洲精品国产一区二区精华液| 亚洲伊人色综图| 国产精品成人在线| 少妇人妻精品综合一区二区| 纵有疾风起免费观看全集完整版| 一本大道久久a久久精品| 欧美最新免费一区二区三区| 午夜91福利影院| 欧美激情极品国产一区二区三区| 日本爱情动作片www.在线观看| 国产不卡av网站在线观看| 男女边吃奶边做爰视频| 搡老岳熟女国产| 99热国产这里只有精品6| 婷婷成人精品国产| 国产精品国产三级专区第一集| 亚洲欧美一区二区三区国产| 视频在线观看一区二区三区| 97精品久久久久久久久久精品| 亚洲国产精品一区二区三区在线| 男人爽女人下面视频在线观看| 亚洲第一区二区三区不卡| 亚洲国产欧美网| 国产探花极品一区二区| 日韩成人av中文字幕在线观看| 久久久精品免费免费高清| 国产男人的电影天堂91| 国产免费又黄又爽又色| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲伊人色综图| 一区二区av电影网| 久久精品国产a三级三级三级| 成人国产麻豆网| 久久久久久免费高清国产稀缺| 免费高清在线观看视频在线观看| 久久精品久久久久久久性| 超色免费av| 韩国av在线不卡| 久久免费观看电影| 麻豆精品久久久久久蜜桃| 久久天躁狠狠躁夜夜2o2o | 国产午夜精品一二区理论片| www.精华液| 美女午夜性视频免费| 操美女的视频在线观看| 国产精品麻豆人妻色哟哟久久| e午夜精品久久久久久久| 九九爱精品视频在线观看| 成人国产麻豆网| 久久人人97超碰香蕉20202| 国产熟女午夜一区二区三区| 女的被弄到高潮叫床怎么办| 国产精品秋霞免费鲁丝片| 国产欧美日韩一区二区三区在线| av免费观看日本| 晚上一个人看的免费电影| 如日韩欧美国产精品一区二区三区| 少妇 在线观看| 一级片免费观看大全| 国产老妇伦熟女老妇高清| 亚洲欧美精品自产自拍| 成年人午夜在线观看视频| 亚洲欧美一区二区三区黑人| 国产老妇伦熟女老妇高清| 国产免费福利视频在线观看| 美女国产高潮福利片在线看| 电影成人av| 老司机在亚洲福利影院| 麻豆av在线久日| 狂野欧美激情性bbbbbb| 国产精品成人在线| 老汉色∧v一级毛片| 国产黄频视频在线观看| 国产精品嫩草影院av在线观看| 欧美97在线视频| 国产色婷婷99| 日韩人妻精品一区2区三区| 丁香六月欧美| 午夜福利一区二区在线看| 99精品久久久久人妻精品| 中文乱码字字幕精品一区二区三区| 欧美老熟妇乱子伦牲交| 少妇人妻久久综合中文| 亚洲精品成人av观看孕妇| 亚洲七黄色美女视频| 精品国产国语对白av| www.av在线官网国产| 亚洲精品久久午夜乱码| 黄片小视频在线播放| 久久久久精品国产欧美久久久 | 波野结衣二区三区在线| 一级毛片电影观看| 午夜福利视频在线观看免费| 婷婷色综合大香蕉| 亚洲国产精品成人久久小说| 国产精品一区二区在线观看99| 亚洲情色 制服丝袜| 国产又爽黄色视频| 久久久精品国产亚洲av高清涩受| 亚洲精品第二区| 欧美xxⅹ黑人| 丁香六月欧美| 亚洲欧美成人精品一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 嫩草影视91久久| 国产精品一国产av| 国产精品99久久99久久久不卡 | 高清欧美精品videossex| 国产精品二区激情视频| 少妇的丰满在线观看| 国产一级毛片在线| 中文天堂在线官网| 色婷婷久久久亚洲欧美| 最近最新中文字幕大全免费视频 | 水蜜桃什么品种好| 精品福利永久在线观看| 国产成人a∨麻豆精品| 考比视频在线观看| 操美女的视频在线观看| 婷婷成人精品国产| 欧美日韩福利视频一区二区| 国产日韩欧美亚洲二区| 亚洲欧洲日产国产| 亚洲国产精品999| 免费观看性生交大片5| 亚洲av电影在线进入| 2021少妇久久久久久久久久久| 国产女主播在线喷水免费视频网站| 国产精品一区二区精品视频观看| 男人操女人黄网站| 精品亚洲成a人片在线观看| 国产一区亚洲一区在线观看| av视频免费观看在线观看| 亚洲av电影在线观看一区二区三区| 99re6热这里在线精品视频| 国产深夜福利视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 精品国产超薄肉色丝袜足j| 操出白浆在线播放| 在线观看人妻少妇| 激情五月婷婷亚洲| 不卡视频在线观看欧美| 在线免费观看不下载黄p国产| 永久免费av网站大全| 久久人人爽人人片av| 中文字幕人妻熟女乱码| videos熟女内射| 国产精品三级大全| 欧美日本中文国产一区发布| 午夜91福利影院| 在线看a的网站| 麻豆精品久久久久久蜜桃| 国产欧美日韩综合在线一区二区| 午夜影院在线不卡| 欧美久久黑人一区二区| 成人亚洲欧美一区二区av| 久久久久精品性色| 日本色播在线视频| 欧美精品高潮呻吟av久久| 熟女少妇亚洲综合色aaa.| 考比视频在线观看| 一二三四中文在线观看免费高清| av.在线天堂| 一级黄片播放器| 美女主播在线视频| 无遮挡黄片免费观看| 国产精品三级大全| 亚洲在久久综合| 一区二区三区乱码不卡18| 亚洲国产精品一区二区三区在线| 在线观看www视频免费| bbb黄色大片| 亚洲欧美色中文字幕在线| av天堂久久9| 亚洲成色77777| 亚洲欧美中文字幕日韩二区| 制服诱惑二区| 少妇精品久久久久久久| 国产精品女同一区二区软件| 美女脱内裤让男人舔精品视频| a 毛片基地| 日本av免费视频播放| 精品人妻一区二区三区麻豆| 国产视频首页在线观看| 国产成人精品在线电影| 亚洲欧美精品综合一区二区三区| 国产欧美日韩综合在线一区二区| 青春草视频在线免费观看| 午夜免费鲁丝| 在线观看免费午夜福利视频| 国产日韩欧美亚洲二区| 亚洲欧洲精品一区二区精品久久久 | 成人亚洲精品一区在线观看| 欧美日韩精品网址| 日本一区二区免费在线视频| h视频一区二区三区| 精品少妇久久久久久888优播| 中文字幕人妻丝袜制服| 啦啦啦啦在线视频资源| 亚洲av男天堂| 尾随美女入室| 午夜福利视频精品| 一级,二级,三级黄色视频| 免费黄网站久久成人精品| 日本猛色少妇xxxxx猛交久久| 九草在线视频观看| 欧美日韩av久久| 一区在线观看完整版| 亚洲第一av免费看| 久久青草综合色| 久久人妻熟女aⅴ| 亚洲成人一二三区av| 国产成人免费观看mmmm| 免费在线观看完整版高清| 国产成人a∨麻豆精品| 亚洲国产欧美一区二区综合| www.自偷自拍.com| 伦理电影免费视频| 国产精品av久久久久免费| 久久久久国产精品人妻一区二区| 菩萨蛮人人尽说江南好唐韦庄| 午夜免费男女啪啪视频观看| 亚洲成人手机| 日韩av不卡免费在线播放| 亚洲精品国产av蜜桃| 国产精品三级大全| 波野结衣二区三区在线| 国产日韩欧美视频二区| 18禁裸乳无遮挡动漫免费视频| 免费看不卡的av| 日韩欧美精品免费久久| 街头女战士在线观看网站| 亚洲国产毛片av蜜桃av| 久久精品久久精品一区二区三区| 天天操日日干夜夜撸| 国产爽快片一区二区三区| 黑人猛操日本美女一级片| 国产男女超爽视频在线观看| 韩国高清视频一区二区三区| 色吧在线观看| 一本一本久久a久久精品综合妖精| 美国免费a级毛片| 香蕉国产在线看| 国产精品久久久av美女十八| 观看美女的网站| 欧美xxⅹ黑人| 中文字幕亚洲精品专区| 国产激情久久老熟女| 久久女婷五月综合色啪小说| 麻豆精品久久久久久蜜桃| 各种免费的搞黄视频| 欧美人与善性xxx| 国产片特级美女逼逼视频| 又大又黄又爽视频免费| 蜜桃在线观看..| 亚洲自偷自拍图片 自拍| 日本av免费视频播放| 亚洲国产看品久久| av卡一久久| 欧美日韩成人在线一区二区| 日本91视频免费播放| 久久99热这里只频精品6学生| 中文天堂在线官网| 日日啪夜夜爽| 观看美女的网站| 亚洲av欧美aⅴ国产| 伊人久久大香线蕉亚洲五| av天堂久久9| 制服丝袜香蕉在线| 国产爽快片一区二区三区| av天堂久久9| 亚洲七黄色美女视频| 国产爽快片一区二区三区| 国产精品一区二区在线不卡| 国产成人精品久久久久久| 熟妇人妻不卡中文字幕| 亚洲国产精品一区三区| 满18在线观看网站| 亚洲图色成人| 亚洲av电影在线观看一区二区三区| 国产免费福利视频在线观看| 自线自在国产av| 黑人欧美特级aaaaaa片| 亚洲七黄色美女视频| 自线自在国产av| 国产精品女同一区二区软件| 天天操日日干夜夜撸| 日日啪夜夜爽| 最近中文字幕高清免费大全6| 日韩一本色道免费dvd| 波多野结衣av一区二区av| 丁香六月欧美| 久久99精品国语久久久| 美女中出高潮动态图| 黄色毛片三级朝国网站| 超色免费av| 国产精品香港三级国产av潘金莲 | 国产一区有黄有色的免费视频| 国产99久久九九免费精品| 国产精品 国内视频|