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

    基于稀疏表征的降噪方法及其在振動激勵控制中的應(yīng)用

    2016-04-01 03:53:16孟利波秦毅合燁郭磊
    中國測試 2016年11期
    關(guān)鍵詞:振動優(yōu)化信號

    孟利波,秦毅,合燁,郭磊

    (國家山區(qū)公路工程技術(shù)研究中心,重慶400044)

    基于稀疏表征的降噪方法及其在振動激勵控制中的應(yīng)用

    孟利波,秦毅,合燁,郭磊

    (國家山區(qū)公路工程技術(shù)研究中心,重慶400044)

    針對振動控制中因混有噪聲而導(dǎo)致控制精度降低,影響振動控制研究效果的現(xiàn)實問題,提出基于稀疏表征的降噪方法,并通過引入峭度作為測度指標改進基于分離變量的增廣拉格朗日收斂算法(SALSA)的基追蹤降噪方法。然后,提出獲取優(yōu)化的拉格朗日乘數(shù)的方法,從而使得基追蹤降噪可以獲得更佳的降噪性能。最后將降噪信號反饋給前端控制器,以達到預(yù)期的控制目的。通過仿真實例,驗證基于SALSA的基追蹤降噪方法的有效性。結(jié)果表明:將該方法用于水潤滑軸承實驗的振動控制,能夠有效地去除所測力矩信號中不同強度的高斯白噪聲,因而可以較好地應(yīng)用在振動控制中的反饋環(huán)節(jié)。

    稀疏表征;降噪;SALSA算法;峭度;振動控制

    0 引言

    在研究旋轉(zhuǎn)機械中的軸承、齒輪等零部件的動態(tài)特性時,往往需要對該部件進行動態(tài)加載。這一過程就是一個振動控制過程。但在振動控制的過程中,扭矩傳感器所采集的力信號往往會受到噪聲的影響,使得加載力帶有隨機誤差,并最終影響加載精度。因此,必須抑制這些噪聲以提高加載器的加載精度。

    目前,在振動控制的降噪處理中,廣泛是采用小波或小波包進行閾值降噪處理[1-2]。Li和Dang[3]提出改進的小波閾值降噪算法對無線通信信號進行降噪。Xu和Wang[4]將小波閾值和解調(diào)技術(shù)相結(jié)合用于激光信號降噪。Chen等[5]利用斯坦無偏風(fēng)險估計組提出重疊多小波組閾值降噪方法,并應(yīng)用于軋機傳動系統(tǒng)狀態(tài)監(jiān)測。但是以上方法均是采用二進小波變換,這雖然能達到無冗余存貯與重建信號的目的,但隨著分解層數(shù)的增加,各層、各頻段序列的數(shù)據(jù)點數(shù)也減半、采樣頻率也減半,當數(shù)據(jù)點數(shù)太少時,信號細節(jié)會丟失,并且小波分解的結(jié)果存在著各頻帶間能量的交疊問題[6-7]。于是,學(xué)者們提出了基于經(jīng)驗?zāi)J椒纸?empirical mode decomposition,EMD)的信號降噪方法。例如,Yang等[8]提出了EMD區(qū)間閾值降噪方法;Tian等[9]提出了清晰首區(qū)間閾值EMD降噪方法(clear first interval thresholded empirical mode decomposition,EMD-CIIT)。但是EMD存在模態(tài)混疊問題,且容易受到強噪聲干擾。

    自然環(huán)境和工程中大部分噪聲都是高斯白噪聲,目前信號降噪方法也主要針對該類噪聲。為了更好地實現(xiàn)自適應(yīng)信號降噪,本文基于基追蹤(BP)問題,采用了一種新的線性凸優(yōu)化算法——SALSA算法[10-11],得到了信號的一種稀疏表示,并通過引入譜峭度作為目標函數(shù)來改進SALSA算法,從而使改進后的SALSA算法在降噪過程中具備更好的適應(yīng)性,將其運用在軸承的振動控制中,降噪效果好,計算效率高。

    1 稀疏信號

    稀疏信號就是指信號可以用少數(shù)幾個特征向量的線性組合進行表示。信號的稀疏表示,就是找到一種簡潔的方式來表示信號,使得絕大部分變換系數(shù)的值接近于零或等于零,從而使得到的變換信號是稀疏或者近似稀疏的。

    對于某一個信號,均可以用欠定方程表示為

    其中,A為M×N矩陣,y為長度為M的向量,x為長度為N的向量,且N>M。

    該系統(tǒng)未知量的個數(shù)多于方程的個數(shù),同時矩陣A的寬度大于其長度,當假定AA*可逆時,則方程組有無窮解。

    為解決該方程,常用的方法是基于最小二乘的方法。對于本文,為使得信號更加稀疏,采用的是基于基追蹤的降噪方法(BPD)。首先需引入范數(shù)l1定義為

    然后求解式(1)的方法就是使得x的絕對值之和最小,即基追蹤(BP)問題,如下式:

    當y含有噪聲時,這種狀況下,需尋求一種近似的目標函數(shù)式,該式即為基追蹤降噪(BPD)問題:

    對于傳統(tǒng)的最小二乘法,它是求取平方和的最小值,相比求取絕對值和的最小值,其對信號中的較大值更加敏感,如圖1所示。因此,當采用最小二乘法時,為保證平方和最小,需要取得少量的較大值,因為相比于小值,它們的影響更大,因此,最后求取的信號中有更多的小值,也就造成了信號的不稀疏。相反,基于基追蹤降噪(BPD)的方法就不會包含很多較小值,從而會獲得更加稀疏的信號。針對高斯白噪聲信號,本文選用冗余傅里葉變換基進行降噪,該變換基A定義為

    式中,0≤m≤M-1和0≤n≤N-1。

    圖1 和x2函數(shù)圖像

    2 稀疏表征算法——SALSA

    由式(4)發(fā)現(xiàn),由于‖x‖1不可微,給計算增加一定難度。并且容易發(fā)現(xiàn),該問題屬于凸優(yōu)化問題,基于這種性質(zhì),就有內(nèi)部局部最小值。此外,對于凸優(yōu)化問題,目前已有很多算法能夠進行求解。例如ISTA[12]和Split Bregman迭代算法[13],其能夠保證在每次迭代后成本函數(shù)值的減小,但是這種算法具有收斂慢的缺點。最近興起的分離變量的增廣拉格朗日收斂算法(SALSA),其與Split Bregman迭代算法均先采用了分離變量,但SALSA的解決方式是基于增廣的拉格朗日模型,而該模型是解決優(yōu)化問題中更基礎(chǔ)更標準的工具。在實踐的過程中也證明了該算法有很好的收斂性質(zhì)。

    在RN上定義兩個凸能量函數(shù)f1(x)和f2(x),其中x為自變量,給出極小化模型:

    為分離變量,引入中間變量v,并將v作為函數(shù)f2的自變量,這樣上式就轉(zhuǎn)化為以下約束問題:

    該問題顯然等價于問題(6)。欲解決該問題,可以利用增廣型拉格朗日模型解決,模型如下:

    式中,λ為拉格朗日乘數(shù),μ≥0被稱為懲罰參數(shù)。

    對于該增廣型拉格朗日模型,可以引入序列dk,進行迭代,從而達到不斷收斂的效果,迭代算法過程如下:

    要想很好地利用上述模型,必須能夠解決式(9),可以采用分裂的方法,分別按x和v迭代求極小,從而有效地求得方程的極小值。具體步驟如下:

    步驟1:

    步驟2:

    為了有效解決本文討論的線性凸優(yōu)化問題公式(4),分別令

    由此可以建立完整的SALSA迭代計算過程,如圖2所示。

    相比傳統(tǒng)的懲罰函數(shù)方法,SALSA迭代算法有很多優(yōu)點。首先,它基于l1正則項的問題,使得其收斂速度很快。當其收斂很快時,就可以直接利用無約束問題。其次,λ是一個固定常數(shù),只需選擇一個合適的值使得條件數(shù)達到最少,從而使迭代優(yōu)化算法快速收斂。通過計算出的x就可以估算出原始信號^y=Ax。

    圖2 SALSA算法流程圖

    3 利用SALSA算法進行去噪仿真

    現(xiàn)有一組含有500個采樣點的音頻信號,其采樣率為16 000,持續(xù)時間約為31 ms,用s(m)表示,圖形如圖3(a)所示,可以看出該音頻信號主要集中在低頻部分,但高頻部分仍有部分微弱的信號。用w(m)表示采樣點數(shù)為500的高斯噪聲信號,由于其具有隨機性,因此會造成信號的不稀疏??捎脃(m)表示含有噪聲的信號如式(12),其圖形如圖3(d)所示,頻域圖涵蓋整個范圍,造成信號的不稀疏。

    圖3 原始信號、加噪信號、降噪信號的時域波形以及對應(yīng)的頻域圖

    然后利用SALSA算法對其進行稀疏表征。采用50次迭代,得到降噪信號,如圖3(e)所示??梢钥闯?,噪聲信號大大削減,而波形的基本形狀幾乎沒發(fā)生變化,并且與低通濾波不同,高頻部分并沒有完全衰減,仍然保留了原始信號中的部分高頻成分。

    4 基于譜峭度的算法優(yōu)化

    前面已經(jīng)討論過,在使用SALSA算法時,需選取合適的λ。經(jīng)過多次試驗發(fā)現(xiàn):若選擇較大的λ,則降噪后的信號,不僅噪聲大大衰減;有用信號也大大衰減;反之,若選擇較小的λ,將不能完全把噪聲濾除。

    峭度(kur)表示故障形成的大幅值脈沖出現(xiàn)的概率,按下式定義為

    式中,Xi為原信號,為信號均值。為了將脈沖響應(yīng)與背景噪聲的差距拉大以提高信噪比,峭度系數(shù)取脈沖響應(yīng)幅值的4次方為判斷依據(jù),準確度大大提高。于是在尋求算法優(yōu)化階段同時引入了信噪比和峭度來尋找最佳λ。

    在仿真實驗中引入余弦信號作為模擬仿真信號,將其與隨機信號進行疊加得到一個模擬試驗中的信號。在仿真實驗中引入一個比例參數(shù)rate:

    從圖中可見,結(jié)論1)是成立的。在同樣的噪聲和λ取值范圍下,令余弦信號幅值取值范圍為(1.0,0.9,…,0.1),此時所得近似最佳λ(1 030),如表1所示。從表中可知,在隨機信號一定時,rate值不會發(fā)生很大的變化,由rate值求出的近似最佳λ也是非常接近在信噪比與λ關(guān)系圖上得出的最佳λ(1 030),于是,在真實信號未知的情況下可以參考模擬仿真信號下得到的rate值反求出一個十分接近最佳真實值的λ,從而獲得優(yōu)化的降噪結(jié)果。

    圖4 余弦信號取不同幅值時降噪結(jié)果信噪比隨λ變化的示意圖

    表1 不同余弦信號幅值對應(yīng)的近似最佳λ

    5 實驗驗證

    將上述降噪方法運用在水潤滑軸承激勵控制中,實驗平臺如圖5所示。該實驗臺能夠通過二通道電液伺服協(xié)調(diào)加載系統(tǒng)對實驗軸進行靜態(tài)加載和動態(tài)激勵。在進行動態(tài)激勵時,可根據(jù)需要,選擇相應(yīng)的波形,如正弦波、三角波、方波等進行激勵。

    圖5 軸承加載實驗臺

    對于本實驗,可采用15 Hz的正弦波進行激勵。通過扭矩傳感器采集到的力矩信號如圖6(a)所示,通過選定一系列的λ值進行降噪,預(yù)估噪聲強度為0.7。根據(jù)前述方法求得真實信號近似最佳λ值為110。通過扭矩傳感器采集到的力矩信號和采用本文方法獲得的降噪信號如圖6(b)所示??梢钥闯?,所加載的力矩由于受到噪聲污染,含有許多毛刺,但經(jīng)過本文的方法處理后,具有明顯的降噪效果。將降噪后的信號反饋給前端控制器,用來不斷修正被控量與輸入量之間的偏差,從而能夠達到更準確的控制效果。

    圖6 低噪聲時實測信號與降噪信號

    為了進一步驗證本文所提方法的有效性,通過調(diào)整信號調(diào)理裝置,增加噪聲強度,所測得信號如圖7(a)所示。采用本文方法所得降噪結(jié)果如圖7(b)所示。同樣地,可以將降噪后的信號反饋給前端控制器,以提高控制精度。

    圖7 高噪聲時實測信號與降噪信號

    6 結(jié)束語

    對于混有噪聲的振動控制信號,利用本文敘述的方法可以達到良好的降噪效果,可見引入峭度和反求最佳λ值的方法,使得降噪效果更加明顯,更加具有適應(yīng)性,最后將降噪信號反饋給前端控制器,用來不斷修正被控量與輸入量之間的偏差,從而能夠提高控制精度。該方法的不足之處在于,為了獲取最佳λ,通常要使λ取值范圍較大,這就會顯著地增加計算量。但是就目前計算機性能而言,該方法仍然具有很快的計算速度。

    [1]劉愛平,鄒立華,戴素亮,等.振動控制研究中的小波包降噪處理[J].振動與沖擊,2008,27(2):95-97.

    [2]馮永勝,董云峰.小波降噪在雷達制導(dǎo)空空導(dǎo)彈導(dǎo)引頭上的應(yīng)用[J].戰(zhàn)術(shù)與導(dǎo)彈技術(shù),2007(6):76-80.

    [3]LI R J,DANG A H.A blind detection scheme based on modified wavelet denoising algorithm for wireless optical communications[J].Optics Communications,2015,353 (15):165-170.

    [4]XU F,WANG Y Q.Signal enhancement of a novel multiaddress codinglidar backscatters based on a combined technique of demodulation and wavelet de-noising[J].Optics and Lasers in Engineering,2015(74):122-129.

    [5]CHEN J L,WAN Z G,PAN J,et al.Customized maximaloverlapmultiwavelet denoising with data-driven group threshold for condition monitoring of rolling mill drivetrain[J].Mechanical Systems and Signal Processing,2016,(68-69):44-67.

    [6]INOUE T,SUEOKA A,KANEMOTO H,et al.Detection of minute signs of a small fault in a periodic or aquasiperiodic signal by the harmonic wavelet transform[J].Mechanical System and Signal Processing,2007,21(5): 2041-2055.

    [7]張文斌,周曉軍,楊先勇,等.基于諧波窗方法的轉(zhuǎn)子軸心軌跡提純[J].振動與沖擊,2009,28(8):74-77.

    [8]YANG G L,LIU Y Y,WANG Y Y,et al.EMD interval thresholding denoising based on similarity measure to select relevant modes[J].Signal Processing,2015(109):95-109

    [9]TIAN P F,ZHANG L,CAO X J,et al.The Application of EMD-CIIT Lidar Signal Denoising Method in Aerosol Detection[J].Procedia Engineering,2015(102):1233-1237.

    [10]AFONSO M V,BIOUCAS D J M,F(xiàn)IGUEIREDO M A T.Fast image recovery using variable splitting and constrained optimization[J].IEEE TransImage Process,2010,19(9): 2345-2356.

    [11]AFONSO M V,BIOUCAS-D J M,F(xiàn)IGUEIREDO M A T.An augmented Lagrangian approach to the constrained optimization formulation of imaging inverse problems[J].IEEE Trans Image Process,2011,20(3):681-695.

    [12]BECKA,TEBOULLE M.A fast iterative shrinkage thresholding algorithm for linear inverse problems[J].SIAM Journal on Imaging Science,2009,2(1):183–202.

    [13]GOLDSTEIN T,OSHER S.The split Bregman algorithm method for L1 regularized problems[J].SIAM Journal on Imaging Sciences,2009,2(2):323-343.

    (編輯:李妮)

    De-noising method based on sparse representation and its application invibration control

    MENG Libo,QIN Yi,HE Ye,GUO Lei
    (National Engineering and Research Center for Mountainous Highways,Chongqing 400044,China)

    Aimed at the fact that the control accuracy declines due to the noise effect on the vibration control,a de-noising method based on sparse representation is put forward to improve the control signal.To obtain a better noise-reduction result,the basis pursuit denoising method based on split variable augmented Lagrangian shrinkage algorithm(SALSA)is improved by using kurtosis as measurement index.Then,a method for obtaining the optimized Lagrange multiplier parameter is proposed,which can improve the denoising performance of basis pursuit denoising.Finally the de-noising signal is fed back into the front controller in order to obtain the expected control effect.The simulation result validates the effectiveness of basis pursuit denoising based on SALSA.Finally,the proposed method is applied to the vibration control of a water water-lubricated bearing test,and the results show that this method can effectively remove the noises with different intensity from the sampled torque signals,thus it can be well used in the feedback element of the vibration control.

    sparse representation;de-noising;SALSA algorithm;kurtosis;vibration control

    A

    1674-5124(2016)11-0094-06

    10.11857/j.issn.1674-5124.2016.11.020

    2016-01-07;

    2016-03-04

    國家山區(qū)公路工程技術(shù)研究中心(GSGZJ-2014-04)

    孟利波(1977-),男,重慶市人,研究員,博士,研究方向為結(jié)構(gòu)監(jiān)測與控制。

    猜你喜歡
    振動優(yōu)化信號
    振動的思考
    超限高層建筑結(jié)構(gòu)設(shè)計與優(yōu)化思考
    民用建筑防煙排煙設(shè)計優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    振動與頻率
    基于FPGA的多功能信號發(fā)生器的設(shè)計
    電子制作(2018年11期)2018-08-04 03:25:42
    中立型Emden-Fowler微分方程的振動性
    日本熟妇午夜| 国产成年人精品一区二区| 国产精品av视频在线免费观看| 色哟哟哟哟哟哟| 国产成人精品久久久久久| 三级男女做爰猛烈吃奶摸视频| 成人亚洲欧美一区二区av| 在线天堂最新版资源| 成人亚洲欧美一区二区av| 国产一级毛片七仙女欲春2| 国产高清不卡午夜福利| 精品不卡国产一区二区三区| 欧美最新免费一区二区三区| 亚洲第一区二区三区不卡| 国产真实伦视频高清在线观看| 哪里可以看免费的av片| 国产真实乱freesex| 国产亚洲av片在线观看秒播厂 | av卡一久久| 在现免费观看毛片| 国产日韩欧美在线精品| 成人亚洲欧美一区二区av| 男人舔女人下体高潮全视频| 国产黄色小视频在线观看| 国产麻豆成人av免费视频| 国产色婷婷99| 不卡一级毛片| 在线免费观看的www视频| 亚洲av熟女| 免费看av在线观看网站| 一个人看视频在线观看www免费| 日日啪夜夜撸| 免费看美女性在线毛片视频| 久久精品夜色国产| 尾随美女入室| 哪个播放器可以免费观看大片| 乱人视频在线观看| 久久人人爽人人片av| 亚洲丝袜综合中文字幕| 日本免费一区二区三区高清不卡| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 一本精品99久久精品77| av福利片在线观看| 人人妻人人澡欧美一区二区| 国产 一区精品| 久久久久久久久久成人| 国内久久婷婷六月综合欲色啪| 中文字幕久久专区| 精品99又大又爽又粗少妇毛片| 精品日产1卡2卡| 男插女下体视频免费在线播放| 1024手机看黄色片| 毛片女人毛片| 成人性生交大片免费视频hd| 久久精品国产亚洲网站| 欧美成人免费av一区二区三区| 久久久久久久久久久丰满| 日韩大尺度精品在线看网址| 国产日韩欧美在线精品| 伦精品一区二区三区| 国产老妇女一区| 日本免费一区二区三区高清不卡| 久久精品国产自在天天线| 亚洲美女视频黄频| 2022亚洲国产成人精品| 国产高清激情床上av| 国产精品日韩av在线免费观看| 18禁在线无遮挡免费观看视频| 亚洲成人久久爱视频| 级片在线观看| 久久亚洲国产成人精品v| 中文亚洲av片在线观看爽| 网址你懂的国产日韩在线| 夫妻性生交免费视频一级片| 亚洲自拍偷在线| 91麻豆精品激情在线观看国产| 欧美一区二区亚洲| 国产亚洲av嫩草精品影院| 蜜桃久久精品国产亚洲av| 国产伦在线观看视频一区| 亚洲成人av在线免费| 国产午夜福利久久久久久| 国产不卡一卡二| 乱人视频在线观看| 在线免费十八禁| 亚洲无线在线观看| 日韩高清综合在线| 国产91av在线免费观看| 欧美最黄视频在线播放免费| 亚洲,欧美,日韩| 中文亚洲av片在线观看爽| 我的老师免费观看完整版| 插逼视频在线观看| 欧美+亚洲+日韩+国产| 久久99精品国语久久久| 有码 亚洲区| 成人午夜高清在线视频| 1024手机看黄色片| 男人狂女人下面高潮的视频| 亚洲国产精品久久男人天堂| a级一级毛片免费在线观看| 久久久久久久久久久免费av| 两个人的视频大全免费| 久久久久久伊人网av| 亚洲国产欧美人成| 国产麻豆成人av免费视频| 国产精品人妻久久久久久| 亚洲欧美精品专区久久| 18禁在线无遮挡免费观看视频| av黄色大香蕉| 91午夜精品亚洲一区二区三区| 色综合站精品国产| 久久精品国产亚洲网站| 亚洲欧美日韩东京热| 久久亚洲精品不卡| 久久国内精品自在自线图片| 91精品一卡2卡3卡4卡| 只有这里有精品99| 成人无遮挡网站| 国产亚洲av片在线观看秒播厂 | 非洲黑人性xxxx精品又粗又长| 国产午夜精品论理片| 中文欧美无线码| 成人午夜高清在线视频| 成人美女网站在线观看视频| 欧美精品一区二区大全| 国产精品日韩av在线免费观看| 久久99热6这里只有精品| 搞女人的毛片| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕av成人在线电影| 中国国产av一级| 国产精品伦人一区二区| 国产精品久久久久久久电影| 国产一区二区亚洲精品在线观看| 欧美丝袜亚洲另类| 99久久人妻综合| 99视频精品全部免费 在线| 久久精品国产清高在天天线| 热99re8久久精品国产| 蜜桃久久精品国产亚洲av| 看免费成人av毛片| 99久久精品热视频| 两个人视频免费观看高清| 免费看日本二区| 99久久九九国产精品国产免费| 亚洲av.av天堂| 精品欧美国产一区二区三| 免费一级毛片在线播放高清视频| 在线观看一区二区三区| 麻豆成人午夜福利视频| 不卡视频在线观看欧美| 亚洲自偷自拍三级| 免费看日本二区| 国产高清不卡午夜福利| 乱人视频在线观看| 最近最新中文字幕大全电影3| 1000部很黄的大片| 亚洲国产精品成人综合色| 一本精品99久久精品77| 男的添女的下面高潮视频| 精品久久久久久久久av| 久久精品国产亚洲网站| 在线天堂最新版资源| 麻豆国产97在线/欧美| 成人鲁丝片一二三区免费| 黄色配什么色好看| 少妇的逼好多水| 日本黄色片子视频| 国产白丝娇喘喷水9色精品| 在线免费观看的www视频| 欧美xxxx性猛交bbbb| 熟妇人妻久久中文字幕3abv| 亚洲国产欧美在线一区| 极品教师在线视频| 久久久久久久久中文| 美女 人体艺术 gogo| 激情 狠狠 欧美| 一边摸一边抽搐一进一小说| 欧美极品一区二区三区四区| 欧美成人免费av一区二区三区| 国产色爽女视频免费观看| 亚洲,欧美,日韩| 婷婷精品国产亚洲av| 久久人人精品亚洲av| 观看美女的网站| 99riav亚洲国产免费| 国产精品一区www在线观看| 直男gayav资源| 亚洲av中文字字幕乱码综合| 插逼视频在线观看| 成人三级黄色视频| 丰满乱子伦码专区| 性色avwww在线观看| 蜜桃久久精品国产亚洲av| 久久亚洲国产成人精品v| 国产爱豆传媒在线观看| 少妇人妻一区二区三区视频| 18禁在线播放成人免费| 三级毛片av免费| 国产精品国产高清国产av| 在线免费观看的www视频| 少妇熟女欧美另类| 在线天堂最新版资源| 日韩强制内射视频| 99热这里只有是精品50| 久久国产乱子免费精品| 免费观看的影片在线观看| 国产成人午夜福利电影在线观看| 成人国产麻豆网| 欧美色视频一区免费| 欧美人与善性xxx| 国产精品国产三级国产av玫瑰| 一级黄片播放器| 国产白丝娇喘喷水9色精品| 长腿黑丝高跟| 亚洲国产精品久久男人天堂| 只有这里有精品99| 日韩高清综合在线| 人妻制服诱惑在线中文字幕| 亚洲欧洲日产国产| 精品人妻一区二区三区麻豆| 精品人妻视频免费看| 男人的好看免费观看在线视频| 日韩在线高清观看一区二区三区| 久久婷婷人人爽人人干人人爱| 亚洲欧美中文字幕日韩二区| h日本视频在线播放| 床上黄色一级片| 国产精品久久久久久久久免| 又粗又爽又猛毛片免费看| 69av精品久久久久久| 免费人成在线观看视频色| 美女大奶头视频| 免费观看的影片在线观看| 久久久久久久久久成人| 久久久久久国产a免费观看| 女的被弄到高潮叫床怎么办| 国产精品久久久久久亚洲av鲁大| 久久国内精品自在自线图片| 一个人看的www免费观看视频| 青春草国产在线视频 | 亚洲精品日韩av片在线观看| 国产美女午夜福利| 全区人妻精品视频| 色播亚洲综合网| 99久久精品一区二区三区| 日韩精品有码人妻一区| 久久久国产成人精品二区| 国产成人精品久久久久久| av视频在线观看入口| 日韩大尺度精品在线看网址| 久久精品国产清高在天天线| 亚洲人成网站在线播| 在线观看av片永久免费下载| 嫩草影院入口| 久久久久久久久久成人| 国产视频首页在线观看| 老女人水多毛片| 美女大奶头视频| 国产精品久久久久久精品电影小说 | 欧美高清成人免费视频www| 高清日韩中文字幕在线| 国产精品蜜桃在线观看 | 在线观看午夜福利视频| 亚洲欧美日韩东京热| 亚洲欧洲日产国产| 日韩欧美国产在线观看| 哪个播放器可以免费观看大片| 亚洲国产色片| 你懂的网址亚洲精品在线观看 | 麻豆av噜噜一区二区三区| 日本黄色视频三级网站网址| 成人午夜高清在线视频| 小蜜桃在线观看免费完整版高清| av又黄又爽大尺度在线免费看 | 十八禁国产超污无遮挡网站| 只有这里有精品99| 最近视频中文字幕2019在线8| 亚洲精品国产成人久久av| 亚洲自拍偷在线| 免费av毛片视频| 中国国产av一级| 日韩一区二区三区影片| 91狼人影院| 久久精品人妻少妇| 丰满的人妻完整版| 欧美日韩乱码在线| 欧美色视频一区免费| 晚上一个人看的免费电影| 91精品一卡2卡3卡4卡| 永久网站在线| 久久久久久伊人网av| 中文字幕熟女人妻在线| 最近2019中文字幕mv第一页| 国产在线精品亚洲第一网站| 欧美日韩综合久久久久久| 3wmmmm亚洲av在线观看| 美女xxoo啪啪120秒动态图| 日本撒尿小便嘘嘘汇集6| 国产精品人妻久久久久久| 久久久午夜欧美精品| 九色成人免费人妻av| 99热6这里只有精品| 小蜜桃在线观看免费完整版高清| 男的添女的下面高潮视频| 久久久精品欧美日韩精品| 青春草国产在线视频 | 欧美zozozo另类| 欧美xxxx黑人xx丫x性爽| 午夜精品一区二区三区免费看| 神马国产精品三级电影在线观看| 日韩av不卡免费在线播放| 狂野欧美激情性xxxx在线观看| 99热全是精品| 国产成人一区二区在线| 久久国产乱子免费精品| 亚洲成人av在线免费| 日本与韩国留学比较| 一本久久精品| 搞女人的毛片| 精品久久久久久久久av| 校园春色视频在线观看| 91久久精品电影网| 日韩av在线大香蕉| 美女被艹到高潮喷水动态| 亚洲精品乱码久久久久久按摩| 一个人免费在线观看电影| 51国产日韩欧美| 人人妻人人澡人人爽人人夜夜 | 欧美性猛交黑人性爽| 桃色一区二区三区在线观看| 日韩在线高清观看一区二区三区| 国产精品无大码| 欧美zozozo另类| 国产又黄又爽又无遮挡在线| 91午夜精品亚洲一区二区三区| av视频在线观看入口| 婷婷色av中文字幕| 老司机影院成人| 午夜免费男女啪啪视频观看| 日产精品乱码卡一卡2卡三| or卡值多少钱| 少妇熟女欧美另类| 日本熟妇午夜| 久久99热这里只有精品18| 亚洲五月天丁香| 国产成人精品一,二区 | 欧美性猛交╳xxx乱大交人| 欧美成人精品欧美一级黄| 国产毛片a区久久久久| 亚洲欧美成人综合另类久久久 | 精品人妻熟女av久视频| 在线播放无遮挡| 国产精品无大码| 国模一区二区三区四区视频| 国产av麻豆久久久久久久| 欧美日韩在线观看h| 精品久久久久久久末码| 麻豆一二三区av精品| 中文欧美无线码| 国产成人a∨麻豆精品| 人妻久久中文字幕网| 国产精品一二三区在线看| 尾随美女入室| 国产精品久久久久久精品电影| 亚洲欧美中文字幕日韩二区| 欧美一级a爱片免费观看看| 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| 日本av手机在线免费观看| 日本黄色片子视频| 国产精品.久久久| 国产一区二区在线观看日韩| 国产精品久久久久久av不卡| 国产亚洲av片在线观看秒播厂 | 国产成人精品婷婷| 嘟嘟电影网在线观看| 国产极品精品免费视频能看的| 国产一区亚洲一区在线观看| 成人av在线播放网站| 国产在线男女| 久久精品久久久久久噜噜老黄 | 天美传媒精品一区二区| 国产久久久一区二区三区| 久久热精品热| 天天躁日日操中文字幕| 国产av麻豆久久久久久久| 久久精品久久久久久噜噜老黄 | 国产亚洲欧美98| 成人特级黄色片久久久久久久| 超碰av人人做人人爽久久| 国产成人精品一,二区 | 免费看日本二区| 1000部很黄的大片| 欧美在线一区亚洲| 少妇的逼水好多| 亚洲一级一片aⅴ在线观看| 中出人妻视频一区二区| 国产精品久久久久久亚洲av鲁大| 麻豆av噜噜一区二区三区| 国产午夜精品久久久久久一区二区三区| 久久精品久久久久久久性| 日本与韩国留学比较| 成人二区视频| 搡女人真爽免费视频火全软件| 日本在线视频免费播放| 亚洲精品亚洲一区二区| 国产精品无大码| 草草在线视频免费看| 狠狠狠狠99中文字幕| 国产 一区 欧美 日韩| 日日啪夜夜撸| 国产精品久久视频播放| 在线国产一区二区在线| 久久中文看片网| 又爽又黄无遮挡网站| 欧美精品一区二区大全| 欧美xxxx黑人xx丫x性爽| 亚洲在久久综合| 久久久久久久午夜电影| 国产色爽女视频免费观看| 男人舔女人下体高潮全视频| 久久精品国产亚洲网站| 岛国毛片在线播放| 波野结衣二区三区在线| 搡女人真爽免费视频火全软件| 国产精品美女特级片免费视频播放器| 成人高潮视频无遮挡免费网站| 波多野结衣高清作品| 亚洲久久久久久中文字幕| 亚洲精品自拍成人| 一区二区三区免费毛片| 亚洲国产日韩欧美精品在线观看| 亚洲熟妇中文字幕五十中出| 在线播放国产精品三级| 黄片无遮挡物在线观看| 成人性生交大片免费视频hd| 国产熟女欧美一区二区| 两个人视频免费观看高清| kizo精华| 一边亲一边摸免费视频| 日韩欧美国产在线观看| www日本黄色视频网| 久久精品国产亚洲网站| 国产高潮美女av| 午夜免费激情av| 国产精品1区2区在线观看.| 午夜免费激情av| 国产女主播在线喷水免费视频网站 | 寂寞人妻少妇视频99o| 麻豆国产av国片精品| 亚洲精品自拍成人| 女的被弄到高潮叫床怎么办| 国产日本99.免费观看| 人人妻人人澡欧美一区二区| 欧美高清性xxxxhd video| av免费在线看不卡| 国产精品电影一区二区三区| 乱人视频在线观看| 国产黄片美女视频| 老司机影院成人| 久久欧美精品欧美久久欧美| 亚洲av不卡在线观看| 亚洲av男天堂| 天堂av国产一区二区熟女人妻| 国产精品一区二区性色av| 超碰av人人做人人爽久久| 大又大粗又爽又黄少妇毛片口| 国产精品无大码| 成人鲁丝片一二三区免费| 亚洲欧美精品专区久久| 中文字幕精品亚洲无线码一区| 亚洲成人久久性| 熟女人妻精品中文字幕| 三级国产精品欧美在线观看| 男人的好看免费观看在线视频| 全区人妻精品视频| 亚洲最大成人手机在线| 亚洲电影在线观看av| 高清毛片免费观看视频网站| 久久久久国产网址| 天美传媒精品一区二区| 熟妇人妻久久中文字幕3abv| 深爱激情五月婷婷| 免费av观看视频| 国国产精品蜜臀av免费| 99热这里只有是精品50| 免费av毛片视频| 亚洲不卡免费看| 99热这里只有是精品在线观看| 国产一级毛片七仙女欲春2| 免费搜索国产男女视频| 人人妻人人看人人澡| 变态另类丝袜制服| 最近最新中文字幕大全电影3| 热99在线观看视频| 午夜福利成人在线免费观看| 亚洲精品日韩av片在线观看| 国产亚洲欧美98| 日产精品乱码卡一卡2卡三| 最近的中文字幕免费完整| 久久热精品热| 午夜爱爱视频在线播放| 麻豆久久精品国产亚洲av| 美女xxoo啪啪120秒动态图| 国产精品电影一区二区三区| 在线免费观看的www视频| 别揉我奶头 嗯啊视频| 黄片wwwwww| 精品日产1卡2卡| 97人妻精品一区二区三区麻豆| a级毛色黄片| 免费看a级黄色片| 国产日本99.免费观看| 国产精品久久久久久亚洲av鲁大| 最近2019中文字幕mv第一页| 99久久久亚洲精品蜜臀av| 1024手机看黄色片| 1000部很黄的大片| 久久九九热精品免费| 亚洲国产高清在线一区二区三| 日本一二三区视频观看| 国产精品久久久久久久电影| 国产蜜桃级精品一区二区三区| 亚洲第一电影网av| 搡女人真爽免费视频火全软件| 极品教师在线视频| 国产伦精品一区二区三区视频9| 亚洲人成网站高清观看| 免费看美女性在线毛片视频| 全区人妻精品视频| 国产视频内射| 久久精品国产清高在天天线| 97超碰精品成人国产| videossex国产| 91久久精品电影网| 午夜福利成人在线免费观看| 哪里可以看免费的av片| 岛国毛片在线播放| 精品人妻熟女av久视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲丝袜综合中文字幕| 色播亚洲综合网| 国产成人精品久久久久久| 女人被狂操c到高潮| 日韩三级伦理在线观看| 熟妇人妻久久中文字幕3abv| 一本一本综合久久| eeuss影院久久| 日韩视频在线欧美| 欧美zozozo另类| 欧美一区二区精品小视频在线| 亚洲av电影不卡..在线观看| 国产国拍精品亚洲av在线观看| 亚洲国产精品合色在线| 成人性生交大片免费视频hd| 成人国产麻豆网| 赤兔流量卡办理| 亚洲精品影视一区二区三区av| 国产精品综合久久久久久久免费| 免费黄网站久久成人精品| 色综合站精品国产| 日本成人三级电影网站| 日本黄色片子视频| 国产高清视频在线观看网站| 成人高潮视频无遮挡免费网站| av女优亚洲男人天堂| 亚洲精品成人久久久久久| 联通29元200g的流量卡| 99久久无色码亚洲精品果冻| 亚洲精品亚洲一区二区| 一本精品99久久精品77| 亚洲av一区综合| 国产免费一级a男人的天堂| 特大巨黑吊av在线直播| 成年女人永久免费观看视频| 夫妻性生交免费视频一级片| 桃色一区二区三区在线观看| 国产高清视频在线观看网站| 国产精品国产高清国产av| 国产麻豆成人av免费视频| 丰满的人妻完整版| 国产成人一区二区在线| 在现免费观看毛片| 国内揄拍国产精品人妻在线| 中文字幕av在线有码专区| 18禁裸乳无遮挡免费网站照片| 免费看a级黄色片| 成人美女网站在线观看视频| 99国产精品一区二区蜜桃av| 国产三级中文精品| 精品人妻一区二区三区麻豆| 在线观看美女被高潮喷水网站| 成人毛片60女人毛片免费| 亚洲久久久久久中文字幕| 中文字幕熟女人妻在线| 国产午夜福利久久久久久| 老司机福利观看| 久久精品综合一区二区三区| 成人亚洲欧美一区二区av| 女人被狂操c到高潮| 久久精品影院6| 免费av不卡在线播放| 精品无人区乱码1区二区| 亚洲精品成人久久久久久| 精品人妻一区二区三区麻豆| 一区福利在线观看| 狂野欧美激情性xxxx在线观看| 国产精品人妻久久久久久| 日韩成人av中文字幕在线观看| 特级一级黄色大片| av卡一久久| 国产一级毛片七仙女欲春2|