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

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

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

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

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

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

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

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

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

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

    0 引言

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

    目前,在振動控制的降噪處理中,廣泛是采用小波或小波包進行閾值降噪處理[1-2]。Li和Dang[3]提出改進的小波閾值降噪算法對無線通信信號進行降噪。Xu和Wang[4]將小波閾值和解調(diào)技術相結合用于激光信號降噪。Chen等[5]利用斯坦無偏風險估計組提出重疊多小波組閾值降噪方法,并應用于軋機傳動系統(tǒng)狀態(tài)監(jiān)測。但是以上方法均是采用二進小波變換,這雖然能達到無冗余存貯與重建信號的目的,但隨著分解層數(shù)的增加,各層、各頻段序列的數(shù)據(jù)點數(shù)也減半、采樣頻率也減半,當數(shù)據(jù)點數(shù)太少時,信號細節(jié)會丟失,并且小波分解的結果存在著各頻帶間能量的交疊問題[6-7]。于是,學者們提出了基于經(jīng)驗模式分解(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)自適應信號降噪,本文基于基追蹤(BP)問題,采用了一種新的線性凸優(yōu)化算法——SALSA算法[10-11],得到了信號的一種稀疏表示,并通過引入譜峭度作為目標函數(shù)來改進SALSA算法,從而使改進后的SALSA算法在降噪過程中具備更好的適應性,將其運用在軸承的振動控制中,降噪效果好,計算效率高。

    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)化問題中更基礎更標準的工具。在實踐的過程中也證明了該算法有很好的收斂性質(zhì)。

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

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

    該問題顯然等價于問題(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的高斯噪聲信號,由于其具有隨機性,因此會造成信號的不稀疏。可用y(m)表示含有噪聲的信號如式(12),其圖形如圖3(d)所示,頻域圖涵蓋整個范圍,造成信號的不稀疏。

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

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

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

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

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

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

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

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

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

    表1 不同余弦信號幅值對應的近似最佳λ

    5 實驗驗證

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

    圖5 軸承加載實驗臺

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

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

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

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

    6 結束語

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

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

    [2]馮永勝,董云峰.小波降噪在雷達制導空空導彈導引頭上的應用[J].戰(zhàn)術與導彈技術,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]張文斌,周曉軍,楊先勇,等.基于諧波窗方法的轉子軸心軌跡提純[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ū)公路工程技術研究中心(GSGZJ-2014-04)

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

    猜你喜歡
    振動優(yōu)化信號
    振動的思考
    科學大眾(2023年17期)2023-10-26 07:39:14
    超限高層建筑結構設計與優(yōu)化思考
    民用建筑防煙排煙設計優(yōu)化探討
    關于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    振動與頻率
    天天愛科學(2020年6期)2020-09-10 07:22:44
    基于FPGA的多功能信號發(fā)生器的設計
    電子制作(2018年11期)2018-08-04 03:25:42
    中立型Emden-Fowler微分方程的振動性
    国产日韩一区二区三区精品不卡 | 成人国语在线视频| 国产亚洲最大av| 日韩免费高清中文字幕av| 纯流量卡能插随身wifi吗| 人妻一区二区av| 日韩一本色道免费dvd| 一本—道久久a久久精品蜜桃钙片| 亚洲精品av麻豆狂野| 建设人人有责人人尽责人人享有的| 精品亚洲成a人片在线观看| 看免费成人av毛片| 精品国产乱码久久久久久小说| 人人澡人人妻人| 免费黄色在线免费观看| av专区在线播放| 天堂中文最新版在线下载| 亚洲国产精品专区欧美| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久久免| 国产一区二区三区综合在线观看 | 国产精品女同一区二区软件| 韩国高清视频一区二区三区| 亚洲欧美成人精品一区二区| 国产极品天堂在线| 欧美日韩视频高清一区二区三区二| 久久久久国产精品人妻一区二区| 18+在线观看网站| 亚洲av欧美aⅴ国产| 国产一区亚洲一区在线观看| 久久久久网色| 赤兔流量卡办理| 男的添女的下面高潮视频| 亚洲精品色激情综合| 人妻人人澡人人爽人人| 日韩视频在线欧美| 婷婷色麻豆天堂久久| 国产片内射在线| 九色亚洲精品在线播放| 亚洲无线观看免费| 日韩中字成人| 午夜老司机福利剧场| 亚洲不卡免费看| www.av在线官网国产| 成人黄色视频免费在线看| 免费大片18禁| 国产欧美日韩一区二区三区在线 | 免费观看的影片在线观看| 极品少妇高潮喷水抽搐| 国产免费福利视频在线观看| av一本久久久久| 狂野欧美白嫩少妇大欣赏| av专区在线播放| 亚洲av欧美aⅴ国产| 色哟哟·www| 国产精品久久久久久精品电影小说| 午夜影院在线不卡| 一级,二级,三级黄色视频| 色哟哟·www| 国产精品一区www在线观看| 国产淫语在线视频| 久久99蜜桃精品久久| 亚洲情色 制服丝袜| 天堂俺去俺来也www色官网| 另类亚洲欧美激情| 亚洲精品日韩av片在线观看| 日本av手机在线免费观看| 亚洲精品456在线播放app| 不卡视频在线观看欧美| 日日撸夜夜添| 九色亚洲精品在线播放| 国产一区二区三区综合在线观看 | 91久久精品国产一区二区成人| 中文字幕亚洲精品专区| 你懂的网址亚洲精品在线观看| 十分钟在线观看高清视频www| 一级毛片我不卡| tube8黄色片| 永久网站在线| 国产亚洲最大av| 久久99热6这里只有精品| 尾随美女入室| 九九久久精品国产亚洲av麻豆| videossex国产| 一个人免费看片子| 日韩av不卡免费在线播放| 免费av中文字幕在线| 最近的中文字幕免费完整| 国产又色又爽无遮挡免| 婷婷色综合大香蕉| 一区二区三区乱码不卡18| 国产高清三级在线| 亚洲精品乱码久久久久久按摩| 亚洲精品国产色婷婷电影| 99热这里只有是精品在线观看| 精品少妇黑人巨大在线播放| 日本wwww免费看| 欧美日韩国产mv在线观看视频| 久久久亚洲精品成人影院| 激情五月婷婷亚洲| 亚洲av不卡在线观看| 女人久久www免费人成看片| 国产伦精品一区二区三区视频9| 王馨瑶露胸无遮挡在线观看| 国产淫语在线视频| 一区在线观看完整版| av卡一久久| 新久久久久国产一级毛片| 欧美日韩一区二区视频在线观看视频在线| 精品一区在线观看国产| 亚洲人成网站在线播| 免费av不卡在线播放| 亚洲av.av天堂| 2021少妇久久久久久久久久久| 成人综合一区亚洲| 国产 一区精品| 啦啦啦视频在线资源免费观看| 伊人久久国产一区二区| 国产免费一级a男人的天堂| 久久久精品94久久精品| 亚洲人成网站在线播| 99国产综合亚洲精品| 国产男人的电影天堂91| 日本黄色片子视频| 99视频精品全部免费 在线| 国产欧美另类精品又又久久亚洲欧美| 麻豆精品久久久久久蜜桃| 国产探花极品一区二区| 亚洲av.av天堂| 美女国产高潮福利片在线看| 久久av网站| 久久精品久久久久久噜噜老黄| 一级毛片我不卡| 精品视频人人做人人爽| 日本与韩国留学比较| 亚洲欧美日韩卡通动漫| 美女中出高潮动态图| 亚洲精品乱码久久久久久按摩| 视频中文字幕在线观看| 日本91视频免费播放| 国产女主播在线喷水免费视频网站| 伦理电影免费视频| 亚洲国产精品国产精品| 亚洲av中文av极速乱| 校园人妻丝袜中文字幕| 国产亚洲av片在线观看秒播厂| 亚洲av日韩在线播放| 国产视频首页在线观看| 99久久人妻综合| 高清黄色对白视频在线免费看| 国产片内射在线| 国产精品蜜桃在线观看| 人妻夜夜爽99麻豆av| 国产精品一区二区在线不卡| videos熟女内射| 久久热精品热| 91精品伊人久久大香线蕉| 一级爰片在线观看| 成人黄色视频免费在线看| 午夜老司机福利剧场| 久久久精品免费免费高清| 中文欧美无线码| 午夜91福利影院| 男人添女人高潮全过程视频| 欧美人与性动交α欧美精品济南到 | 国产在线视频一区二区| 777米奇影视久久| 国产视频内射| 狠狠婷婷综合久久久久久88av| 大码成人一级视频| 中文字幕制服av| 成年美女黄网站色视频大全免费 | 日韩熟女老妇一区二区性免费视频| 各种免费的搞黄视频| 亚洲五月色婷婷综合| 欧美激情极品国产一区二区三区 | 在现免费观看毛片| 久久鲁丝午夜福利片| 七月丁香在线播放| 欧美成人精品欧美一级黄| 一级毛片黄色毛片免费观看视频| 免费大片18禁| 日韩精品免费视频一区二区三区 | 国产欧美另类精品又又久久亚洲欧美| 特大巨黑吊av在线直播| 母亲3免费完整高清在线观看 | 久久精品国产自在天天线| 最近手机中文字幕大全| 十八禁网站网址无遮挡| 欧美另类一区| 狠狠精品人妻久久久久久综合| 欧美三级亚洲精品| 国精品久久久久久国模美| 精品久久久噜噜| 成人手机av| 精品久久久精品久久久| 一个人看视频在线观看www免费| 一级毛片 在线播放| 97精品久久久久久久久久精品| 国产精品一区二区三区四区免费观看| 搡老乐熟女国产| 久久久久久人妻| 人妻夜夜爽99麻豆av| 日本爱情动作片www.在线观看| 又黄又爽又刺激的免费视频.| 2022亚洲国产成人精品| 亚洲综合精品二区| 九草在线视频观看| 久久ye,这里只有精品| 十八禁高潮呻吟视频| 午夜日本视频在线| 亚洲成人手机| 欧美+日韩+精品| 中国三级夫妇交换| 亚洲av不卡在线观看| 边亲边吃奶的免费视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品一二三| 亚洲伊人久久精品综合| 超碰97精品在线观看| 免费观看在线日韩| 日日摸夜夜添夜夜爱| 久久婷婷青草| 色吧在线观看| 晚上一个人看的免费电影| 成人毛片60女人毛片免费| 午夜免费鲁丝| 日韩精品免费视频一区二区三区 | 亚洲综合色网址| 在线亚洲精品国产二区图片欧美 | 欧美成人精品欧美一级黄| 纵有疾风起免费观看全集完整版| 一级毛片我不卡| 高清视频免费观看一区二区| 一级a做视频免费观看| 99热这里只有精品一区| 欧美精品人与动牲交sv欧美| 婷婷色综合www| 十八禁网站网址无遮挡| 国产日韩一区二区三区精品不卡 | 亚洲国产av影院在线观看| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕| 久久久久精品性色| 国产成人精品福利久久| 欧美性感艳星| 久久这里有精品视频免费| 精品一区二区免费观看| 飞空精品影院首页| 在线观看免费日韩欧美大片 | 国产黄频视频在线观看| 水蜜桃什么品种好| 大话2 男鬼变身卡| 亚洲综合精品二区| 久久国产精品大桥未久av| 亚洲情色 制服丝袜| 成人综合一区亚洲| 亚洲精品色激情综合| 狠狠精品人妻久久久久久综合| 亚洲av中文av极速乱| 免费观看的影片在线观看| 日本黄色日本黄色录像| 国产日韩一区二区三区精品不卡 | 我的女老师完整版在线观看| 18禁观看日本| 制服诱惑二区| 久久久精品区二区三区| 2022亚洲国产成人精品| 久久久久国产精品人妻一区二区| 国产免费一级a男人的天堂| 热re99久久国产66热| 亚洲国产精品999| 国产成人精品婷婷| 亚洲第一av免费看| 我的女老师完整版在线观看| 欧美 日韩 精品 国产| 午夜日本视频在线| 色婷婷久久久亚洲欧美| 亚洲成人av在线免费| 欧美日韩视频高清一区二区三区二| 精品一品国产午夜福利视频| 午夜日本视频在线| 精品人妻熟女av久视频| 两个人的视频大全免费| 99久国产av精品国产电影| 日韩熟女老妇一区二区性免费视频| 亚洲av电影在线观看一区二区三区| 国产一区有黄有色的免费视频| 成人国产麻豆网| 欧美日韩av久久| 你懂的网址亚洲精品在线观看| 高清不卡的av网站| www.av在线官网国产| 日本猛色少妇xxxxx猛交久久| 精品亚洲成国产av| 精品一区二区免费观看| av线在线观看网站| 99热国产这里只有精品6| 下体分泌物呈黄色| www.av在线官网国产| 成人影院久久| 久久久精品免费免费高清| 免费观看的影片在线观看| 国产欧美亚洲国产| 校园人妻丝袜中文字幕| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 亚洲精华国产精华液的使用体验| 一级毛片电影观看| 午夜视频国产福利| 欧美少妇被猛烈插入视频| 天天影视国产精品| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜免费资源| 国产国语露脸激情在线看| 在线观看www视频免费| 母亲3免费完整高清在线观看 | 大片电影免费在线观看免费| 中文字幕久久专区| 大片电影免费在线观看免费| 在线观看一区二区三区激情| 成人国产av品久久久| 国产男女超爽视频在线观看| 一二三四中文在线观看免费高清| 少妇高潮的动态图| 啦啦啦视频在线资源免费观看| 国产精品99久久久久久久久| 综合色丁香网| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 国产毛片在线视频| 免费大片18禁| 国产在线视频一区二区| 日韩成人av中文字幕在线观看| 日韩成人伦理影院| 黑人欧美特级aaaaaa片| 国产探花极品一区二区| 一区二区三区四区激情视频| 欧美激情极品国产一区二区三区 | 亚洲av不卡在线观看| 欧美日韩在线观看h| 久久99热6这里只有精品| 这个男人来自地球电影免费观看 | 看非洲黑人一级黄片| 999精品在线视频| 国产精品偷伦视频观看了| 日韩一区二区三区影片| 亚洲精品久久久久久婷婷小说| 午夜激情av网站| 国产成人精品无人区| 我要看黄色一级片免费的| 国产有黄有色有爽视频| 久久午夜福利片| 精品酒店卫生间| 久久影院123| 国产av一区二区精品久久| 日韩av免费高清视频| 日本欧美视频一区| 国产精品三级大全| 18+在线观看网站| 国产亚洲最大av| 热99久久久久精品小说推荐| 美女xxoo啪啪120秒动态图| 久久毛片免费看一区二区三区| 免费av中文字幕在线| 两个人免费观看高清视频| 国产精品三级大全| 日韩电影二区| 91久久精品国产一区二区三区| 80岁老熟妇乱子伦牲交| 搡老乐熟女国产| 97超碰精品成人国产| 激情五月婷婷亚洲| 伦精品一区二区三区| 国产成人av激情在线播放 | 国产片内射在线| 国产精品人妻久久久久久| 成人影院久久| 国产探花极品一区二区| 一二三四中文在线观看免费高清| 亚洲经典国产精华液单| 插阴视频在线观看视频| 国产免费一级a男人的天堂| 亚洲国产色片| 亚洲精品视频女| 51国产日韩欧美| 日韩av不卡免费在线播放| 亚洲五月色婷婷综合| 蜜臀久久99精品久久宅男| 99热这里只有是精品在线观看| 成人手机av| 久热久热在线精品观看| 免费播放大片免费观看视频在线观看| 亚洲国产成人一精品久久久| 99热网站在线观看| 久久久久久久亚洲中文字幕| 一边亲一边摸免费视频| 另类亚洲欧美激情| 日韩成人av中文字幕在线观看| 人妻制服诱惑在线中文字幕| 国产精品久久久久久av不卡| 九九在线视频观看精品| 曰老女人黄片| 欧美三级亚洲精品| 精品少妇内射三级| 午夜福利视频在线观看免费| 国产成人精品一,二区| 亚洲精品乱码久久久久久按摩| 好男人视频免费观看在线| 亚洲精品久久久久久婷婷小说| 黄色毛片三级朝国网站| 亚洲人成网站在线观看播放| 亚州av有码| 蜜桃久久精品国产亚洲av| 国产黄色免费在线视频| 国产精品久久久久成人av| 日韩,欧美,国产一区二区三区| 亚洲综合色网址| 久久精品国产鲁丝片午夜精品| 午夜av观看不卡| 欧美xxxx性猛交bbbb| 免费av不卡在线播放| 在线观看美女被高潮喷水网站| 久久婷婷青草| 最近最新中文字幕免费大全7| 丝袜在线中文字幕| 大话2 男鬼变身卡| 王馨瑶露胸无遮挡在线观看| 99re6热这里在线精品视频| √禁漫天堂资源中文www| 午夜福利视频在线观看免费| 欧美日韩在线观看h| 亚洲精品日本国产第一区| 久久久久国产网址| 国产精品一区二区在线不卡| 久久狼人影院| 精品久久久久久电影网| 建设人人有责人人尽责人人享有的| 久久99一区二区三区| av电影中文网址| 夫妻性生交免费视频一级片| 欧美日韩视频高清一区二区三区二| 考比视频在线观看| 亚洲国产欧美在线一区| 蜜桃久久精品国产亚洲av| 街头女战士在线观看网站| 五月天丁香电影| 熟女电影av网| 99久久综合免费| 日本黄色日本黄色录像| 久久午夜综合久久蜜桃| 九九久久精品国产亚洲av麻豆| 亚洲精品成人av观看孕妇| 三级国产精品欧美在线观看| 国模一区二区三区四区视频| 日韩,欧美,国产一区二区三区| 亚洲国产精品一区三区| av专区在线播放| 亚洲国产欧美在线一区| 涩涩av久久男人的天堂| 搡女人真爽免费视频火全软件| 免费看av在线观看网站| 亚洲精品自拍成人| 国产精品一区www在线观看| 免费观看在线日韩| 一级毛片我不卡| 波野结衣二区三区在线| 亚洲精品久久午夜乱码| 91精品一卡2卡3卡4卡| 高清av免费在线| 国产高清三级在线| 亚洲av.av天堂| 久久精品国产a三级三级三级| 3wmmmm亚洲av在线观看| 婷婷色av中文字幕| 蜜臀久久99精品久久宅男| 久久精品夜色国产| 美女大奶头黄色视频| 老司机影院成人| 夜夜看夜夜爽夜夜摸| 能在线免费看毛片的网站| 免费大片黄手机在线观看| 中国国产av一级| 久久久欧美国产精品| 欧美国产精品一级二级三级| 久久精品熟女亚洲av麻豆精品| 一本—道久久a久久精品蜜桃钙片| 老熟女久久久| 视频区图区小说| 精品99又大又爽又粗少妇毛片| 91aial.com中文字幕在线观看| 人妻系列 视频| 一区二区三区精品91| 97精品久久久久久久久久精品| 久久免费观看电影| 国产精品99久久99久久久不卡 | 亚洲精品日韩av片在线观看| 热re99久久精品国产66热6| 一区二区三区精品91| 日韩一区二区三区影片| 亚洲av中文av极速乱| 国产精品一区www在线观看| 国产女主播在线喷水免费视频网站| 亚洲四区av| 高清黄色对白视频在线免费看| 女性被躁到高潮视频| 满18在线观看网站| 精品一区二区三区视频在线| 精品国产露脸久久av麻豆| 国产av国产精品国产| 亚洲国产日韩一区二区| 91在线精品国自产拍蜜月| 人妻制服诱惑在线中文字幕| 一区二区三区四区激情视频| 国产成人精品福利久久| 91精品国产国语对白视频| 性高湖久久久久久久久免费观看| 国产av国产精品国产| 女性生殖器流出的白浆| 精品少妇内射三级| .国产精品久久| 亚洲精品国产av成人精品| 永久网站在线| 亚洲欧美中文字幕日韩二区| 日韩av不卡免费在线播放| 亚洲精品第二区| 中文欧美无线码| 女的被弄到高潮叫床怎么办| av福利片在线| 18禁观看日本| 亚洲伊人久久精品综合| 一个人免费看片子| 久久人人爽人人爽人人片va| 久久这里有精品视频免费| 91久久精品国产一区二区三区| 国产又色又爽无遮挡免| 丰满乱子伦码专区| 久久久欧美国产精品| 亚洲欧洲精品一区二区精品久久久 | 在线免费观看不下载黄p国产| 我的老师免费观看完整版| 简卡轻食公司| 亚洲激情五月婷婷啪啪| 国产视频内射| 国产熟女午夜一区二区三区 | 日本wwww免费看| 中文字幕人妻熟人妻熟丝袜美| 我的女老师完整版在线观看| 青春草视频在线免费观看| 国产精品欧美亚洲77777| 欧美变态另类bdsm刘玥| 一级二级三级毛片免费看| 十八禁网站网址无遮挡| 午夜福利视频精品| 九九在线视频观看精品| 午夜福利影视在线免费观看| 久久人妻熟女aⅴ| 黄色视频在线播放观看不卡| 这个男人来自地球电影免费观看 | 国产精品人妻久久久久久| 国产综合精华液| 日韩,欧美,国产一区二区三区| 卡戴珊不雅视频在线播放| 欧美最新免费一区二区三区| 久久国内精品自在自线图片| 人妻夜夜爽99麻豆av| 最近2019中文字幕mv第一页| freevideosex欧美| 久久久久久久大尺度免费视频| 一级片'在线观看视频| 久久鲁丝午夜福利片| 亚洲欧美一区二区三区国产| 美女国产视频在线观看| 欧美激情 高清一区二区三区| 免费观看的影片在线观看| 婷婷色综合大香蕉| 久久久精品免费免费高清| 亚洲情色 制服丝袜| 国产精品一区二区三区四区免费观看| 国模一区二区三区四区视频| 天堂中文最新版在线下载| av.在线天堂| 亚洲国产精品成人久久小说| 99热国产这里只有精品6| 热99久久久久精品小说推荐| 高清欧美精品videossex| 美女主播在线视频| 久久毛片免费看一区二区三区| 国产在视频线精品| 欧美日韩精品成人综合77777| 午夜视频国产福利| 亚洲欧美清纯卡通| 日本午夜av视频| 日韩一区二区视频免费看| videossex国产| 亚洲熟女精品中文字幕| 人妻 亚洲 视频| 欧美日韩国产mv在线观看视频| 国产精品一国产av| 男女无遮挡免费网站观看| 99久久综合免费| 如日韩欧美国产精品一区二区三区 | 国产精品久久久久久久电影| 特大巨黑吊av在线直播| 91久久精品国产一区二区三区| 精品熟女少妇av免费看| 18禁动态无遮挡网站| 搡老乐熟女国产| 我的老师免费观看完整版| 久久精品国产鲁丝片午夜精品| 人人妻人人添人人爽欧美一区卜| 两个人的视频大全免费| 亚洲少妇的诱惑av| 亚洲色图综合在线观看| 亚洲精品久久久久久婷婷小说|