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

    基于支持向量機(jī)的壓電微操作平臺非線性特性描述

    2016-12-13 04:33:47胡俊峰鄭昌虎蔡建陽
    中國機(jī)械工程 2016年22期
    關(guān)鍵詞:分點(diǎn)驅(qū)動器柵格

    胡俊峰 鄭昌虎 蔡建陽

    江西理工大學(xué),贛州,341000

    ?

    基于支持向量機(jī)的壓電微操作平臺非線性特性描述

    胡俊峰 鄭昌虎 蔡建陽

    江西理工大學(xué),贛州,341000

    壓電驅(qū)動器的遲滯現(xiàn)象會使微操作平臺出現(xiàn)非線性問題,嚴(yán)重影響了其運(yùn)動精度和重復(fù)定位精度。為了解決該問題,提出了一種基于Preisach模型與支持向量機(jī)的聯(lián)合建模方法來描述微操作平臺的非線性特性。以一種一維微操作平臺為對象,以壓電驅(qū)動電壓和所對應(yīng)的平臺輸出位移為樣本點(diǎn),采用支持向量機(jī)理論建立反映平臺非線性回歸模型,利用該模型預(yù)測非樣本點(diǎn)所對應(yīng)的平臺輸出位移,結(jié)合Preisach模型可精確預(yù)測平臺在任意電壓序列作用下的輸出位移。為了驗(yàn)證所建立模型的有效性,進(jìn)行了實(shí)驗(yàn)研究,任意選取2組不同的輸入電壓序列,利用支持向量機(jī)回歸模型和Preisach模型分別得到所選取的電壓序列對應(yīng)的輸出位移的預(yù)測值,在相同的電壓序列作用下進(jìn)行實(shí)驗(yàn)得到其實(shí)測值,將實(shí)測值與預(yù)測值進(jìn)行比較分析,結(jié)果表明,2組實(shí)測值與預(yù)測值之間的相對誤差范圍分別為0.6%~2.1%、0.02%~2.1%,預(yù)測位移與實(shí)測位移非常接近,說明所建立的模型能精確描述微操作平臺的非線性特性,以實(shí)現(xiàn)其精確運(yùn)動。

    微操作平臺;壓電驅(qū)動器;Preisach模型;支持向量機(jī);非線性特性

    0 引言

    目前,高精度、高分辨率的精密微操作平臺在精密定位、微機(jī)械操作、醫(yī)療器械、超精密加工、航空航天等前沿領(lǐng)域應(yīng)用越來越廣泛[1-2]。柔順機(jī)構(gòu)是利用柔性鉸鏈的彈性變形傳遞運(yùn)動和力的一種新型免裝配機(jī)構(gòu)。柔順機(jī)構(gòu)具有無摩擦、無間隙、分辨率高和加工簡單等優(yōu)點(diǎn),特別適宜作為微操作平臺機(jī)構(gòu)[2]。壓電陶瓷驅(qū)動器(PZT)具有剛度高、響應(yīng)速度快、功耗低、輸出力大等優(yōu)點(diǎn)[3-6],PZT和柔順機(jī)構(gòu)結(jié)合而成的微操作平臺具有廣泛應(yīng)用前景。但是,壓電驅(qū)動器在電場作用下會出現(xiàn)遲滯現(xiàn)象,使其具有多值映射性和記憶性等非線性特性,使得微操作平臺出現(xiàn)非線性問題,致使其重復(fù)精度降低,嚴(yán)重影響平臺的控制精度。目前,微操作平臺的模型僅考慮了柔順機(jī)構(gòu)的力學(xué)模型,沒有考慮壓電驅(qū)動器的非線性特性對平臺性能的影響,因此,需要建立能反映壓電驅(qū)動器遲滯非線性特性對平臺影響的模型,以便基于該模型對平臺進(jìn)行精確位移控制。

    為了建立反映遲滯非線性特性的模型,國內(nèi)外學(xué)者提出了Maxwell滑動模型[3]、Duhem模型[4]、Prandtl-Ishlinskii模型[5]和Preisach模型[6-9]等。但是這些模型僅用于反映壓電驅(qū)動器的非線性,沒有應(yīng)用于壓電微操作平臺。其中,Preisach模型是目前應(yīng)用廣泛的描述非線性現(xiàn)象的模型之一,它具有原理清晰、結(jié)構(gòu)簡單、精度高、易于和控制方法相結(jié)合等優(yōu)點(diǎn)。但是,建立能精確反映遲滯非線性的Preisach模型的關(guān)鍵是對遲滯環(huán)進(jìn)行回歸擬合,目前建立回歸模型的方法主要有神經(jīng)網(wǎng)絡(luò)[7-9]、模糊算法[10]、滑模算法[11]和插值算法[12-14]等。但是,采用以上建模方法均需要較大的樣本數(shù)據(jù)量,由于小樣本數(shù)據(jù)量會造成神經(jīng)網(wǎng)絡(luò)建模過程中出現(xiàn)“過擬合”現(xiàn)象和模糊算法的模糊規(guī)則不清等問題,滑模算法增加了Preisach模型的復(fù)雜性。支持向量機(jī)可以解決采用小樣本數(shù)據(jù)回歸非線性擬合問題[15],所以支持向量機(jī)回歸理論適合于建立壓電微操作平臺的非線性模型。同時(shí),支持向量機(jī)從結(jié)構(gòu)風(fēng)險(xiǎn)最小化的角度保證了模型的最大泛化能力,根據(jù)該理論建立的模型能精確預(yù)測驅(qū)動器在任意輸入電壓所對應(yīng)的平臺輸出位移。

    本文以一種一維微操作平臺為對象,采用Preisach模型和支持向量機(jī)聯(lián)合建模的方法描述由壓電陶瓷驅(qū)動器遲滯現(xiàn)象所引起的微操作平臺的非線性特性。首先,建立反映微操作平臺非線性特性的離散化形式的Preisach模型;然后,采用支持向量機(jī)理論建立反映壓電驅(qū)動器電壓與所對應(yīng)的平臺輸出位移之間非線性關(guān)系的回歸模型,根據(jù)該模型可預(yù)測任一電壓所對應(yīng)的輸出位移,從而根據(jù)Preisach模型計(jì)算出壓電驅(qū)動器任意輸入電壓系列的對應(yīng)的平臺輸出位移;最后,進(jìn)行實(shí)驗(yàn)驗(yàn)證分析,以證明所提出方法的可行性和有效性。

    1 壓電微操作平臺的Preisach模型

    圖1所示為一維微操作平臺[2],由于壓電驅(qū)動器的遲滯現(xiàn)象會使平臺出現(xiàn)圖2所示的非線性特性,圖2中(αi,βi)為施加在壓電驅(qū)動器第i個電壓極值點(diǎn)(i=1,2,…,n;αi≥βi),當(dāng)αi>βi時(shí),(αi,βi)表示電壓從0上升αi后再下降到βi,x(αi,βi)為此時(shí)平臺的輸出位移,由圖2可知,在相同電壓驅(qū)動下平臺的輸出位移不相同,取決于施加電壓順序。f(αi,βi)為驅(qū)動器電壓從0上升至αi后下降至βi時(shí)的位移變化量,它可表示為

    f(αi,βi)=x(αi,αi)-x(αi,βi)

    (1)

    圖1 微操作平臺

    圖2 微操作平臺的非線性特性

    圖2所示的非線性可采用Preisach模型描述,圖2中箭頭朝上的曲線為電壓上升階段,箭頭朝下的曲線為電壓下降階段,離散化Preisach模型[12]可表示為

    (2)

    式中,u(t)、f(t)分別為壓電驅(qū)動器在t時(shí)刻的輸入電壓和平臺的位移輸出值;n為飽和電壓的等分個數(shù)。

    由式(1)和式(2)可知,欲求出在t時(shí)刻的輸出位移f(t),需要獲得x(αi,βi)。通過實(shí)驗(yàn)得到x(αi,βi)的過程如下:首先將壓電陶瓷驅(qū)動器的飽和電壓U平均分為n等分,每個等分點(diǎn)的電壓值為kU/n(k=1,2,…,n),將壓電陶瓷驅(qū)動器的電壓從0上升至任一等分點(diǎn)kU/n,記下此時(shí)的電壓(kU/n,kU/n)和平臺輸出位移x(kU/n,kU/n);然后,將電壓下降到0,記下下降過程中經(jīng)過的等分點(diǎn)上的電壓(kU/n,mU/n)和平臺輸出位移x(kU/n,mU/n),m=k-1,k-2,…,1。遍歷所有等分點(diǎn),可以得到壓電陶瓷驅(qū)動器在等分點(diǎn)上的電壓與所對應(yīng)的位移。

    現(xiàn)通過一實(shí)例說明上述模型,設(shè)壓電驅(qū)動器的飽和電壓為60 V,將其均分為3等分,等分點(diǎn)電壓分別為0、20 V、40 V、60 V,如圖3所示。將壓電驅(qū)動器的電壓從0上升至60 V,記下此時(shí)極值電壓(60,60)和輸出位移x(60,60),然后將電壓從60 V下降至0,記下下降過程經(jīng)過的等分點(diǎn)電壓極值點(diǎn)(60,40)、(60,20)、(60,0)和相應(yīng)的平臺輸出位移x(60,40)、x(60,20)、x(60,0)。同理,將電壓從0上升至40 V再降到0,得到x(40,40)、x(40,20)和x(40,0)。最后將電壓上升至20 V再減為0,得到x(20,20)與x(20,0)。

    圖3 說明Preisach模型的實(shí)例

    由圖3可知,對壓電陶瓷驅(qū)動器施加的電壓都在等分點(diǎn)上,可以得到各電壓對應(yīng)的輸出位移。設(shè)壓電驅(qū)動器輸入0-60 V-20 V-40 V的序列電壓,其極值點(diǎn)電壓均在等分點(diǎn)上,根據(jù)式(2)計(jì)算平臺輸出位移的步驟如下:

    (1)根據(jù)施加電壓的順序可確定α1= 60 V,α2= 40 V,β0= 0,β1= 20 V,β2= 0,u(t)= 20 V。

    (2)當(dāng)i=1時(shí),由式(2)可得

    f(α1,β0)-f(α1,β1)=f(60,0)-f(60,20)=

    [x(60,60)-x(60,0)]-[x(60,60)-x(60,20)]

    (3)

    (3)當(dāng)i=2時(shí),由式(2)可得

    f(α2,β1)-f(α1,β2)=f(40,20)-f(40,0)=

    [x(40,40)-x(40,20)]-[x(40,40)-x(40,0)]

    (4)

    (4)由式(2)可計(jì)算

    f(u(t),β2)=f(20,0)=x(20,20)-x(20,0)

    (5)

    (5)微操作平臺的輸出位移為

    x(60,20)-x(60,0)+x(40,0)-

    x(40,20)+x(20,20)-x(20,0)

    (6)

    根據(jù)上述求解過程,當(dāng)輸入電壓序列均在等分點(diǎn)上時(shí),可以通過直接利用在等分點(diǎn)上的x(α,β)計(jì)算壓電驅(qū)動器所對應(yīng)的平臺輸出位移。但是,當(dāng)電壓不在等分點(diǎn)上時(shí),例如當(dāng)壓電驅(qū)動器輸入電壓順序?yàn)?-50 V-30 V-45 V-0-10 V,利用式(2)計(jì)算其輸出位移為x(50,30)-x(50,0)+x(45,0)-x(45,30)+x(10,10)-x(10,0)。但各x(α,β)無法得到,可采用支持向量機(jī)建立反映圖3所示的輸出位移x(α,β)與輸入電壓(α,β)之間關(guān)系的遲滯環(huán)Ⅰ、Ⅱ、Ⅲ的回歸模型,然后根據(jù)該模型預(yù)測任一點(diǎn)的x(α,β)值,從而根據(jù)式(2)計(jì)算出平臺的輸出位移。

    2 基于支持向量機(jī)的非線性模型

    為了建立反映微操作平臺非線性回歸模型,等分點(diǎn)上驅(qū)動器的輸入電壓與平臺輸出位移構(gòu)成支持向量機(jī)的樣本點(diǎn){[(α1,β1),x1(α1,β1)],…,[(αm,βm),xm(αm,βm)]},m為樣本點(diǎn)數(shù)。由于輸出位移x(α,β)與輸入電壓(α,β)之間關(guān)系為非線性,采用支持向量機(jī)理論通過非線性映射函數(shù)將由不同輸入電壓順序的(α,β)形成的輸入樣本空間映射到高維特征空間,在高維特征空間構(gòu)造反映微操作平臺非線性的回歸函數(shù)為

    f(z)=ω·φ(z)+b

    (7)

    式中,ω為超平面系數(shù);z為某遲滯環(huán)對應(yīng)的電壓(α,β);b為待定偏置量;φ()為未知的高維函數(shù);ω·φ(z)為ω與φ(z)的內(nèi)積;f(z)為x(α,β)。

    (8)

    (9)

    其中,R(zi,zk)表示核函數(shù),R(zi,zk)=φ(zi)·φ(zk),即φ(zi)與φ(zk)的內(nèi)積,選取徑向基(RBF)核函數(shù),它可表示為

    R(z,zi)=exp(-g‖z-zi‖2)

    (10)

    式中,g為核函數(shù)的寬度參數(shù)。

    (11)

    偏置量b可通過Karush-Kuhn-Tucker條件求得

    (12)

    聯(lián)立式(7)、式(11)和式(12)可得反映微操作平臺非線性的回歸函數(shù):

    (13)

    為了衡量式(13)所示的回歸模型的精度,采用平方相關(guān)系數(shù)進(jìn)行評價(jià),它可表示為

    r2=

    (14)

    式中,k為測試點(diǎn)數(shù);f(zi)、yi分別為第i個測試點(diǎn)的預(yù)測位移和實(shí)測位移。

    由支持向量機(jī)理論[16]可知,懲罰系數(shù)c和核函數(shù)寬度參數(shù)g對能否建立精確反映微操作平臺非線性模型具有重要作用。懲罰系數(shù)c用于權(quán)衡所建立的模型復(fù)雜度和精度,由式(8)可知,系數(shù)c越大,對樣本點(diǎn)的擬合程度越高,模型復(fù)雜度高,所建立的反映平臺非線性的模型泛化能力變差;系數(shù)c越小,對樣本點(diǎn)數(shù)據(jù)中超出模型與實(shí)際輸出位移的誤差要求ε的樣本懲罰越小,使訓(xùn)練誤差變大,建立模型的精度不高[15]。核函數(shù)寬度g與樣本點(diǎn)的輸入空間范圍有關(guān),樣本點(diǎn)越多取值越大,反之取值越小。所以,需要根據(jù)實(shí)驗(yàn)樣本點(diǎn)優(yōu)化選擇合適的系數(shù)c和g,以建立能精確預(yù)測在不同驅(qū)動電壓所對應(yīng)的輸出位移。

    目前對參數(shù)c和g優(yōu)化的方法主要有梯度法、粒子群法等,相比這些優(yōu)化方法,柵格化搜索方法的優(yōu)點(diǎn)是迭代搜索簡單,能同時(shí)搜索多個參數(shù)并保證其在劃分網(wǎng)格中的全局最優(yōu)解,便于并行計(jì)算,運(yùn)行效率高,所以采用柵格搜索方法尋找能精確建立平臺回歸模型的最優(yōu)參數(shù)c和g。

    由于通過位移傳感器采集的樣本點(diǎn)數(shù)據(jù)具有隨機(jī)誤差,該誤差會影響參數(shù)(c,g)的性能評價(jià),對于同一組(c,g),擬合函數(shù)的精確性會隨樣本數(shù)據(jù)變化而變化,尤其對于小樣本數(shù)據(jù),樣本數(shù)據(jù)的隨機(jī)性對參數(shù)優(yōu)選的影響較大,不利于模型的推廣。因此,在柵格搜索過程中應(yīng)用交叉驗(yàn)證法對每組參數(shù)組合的性能進(jìn)行綜合評價(jià)。將交叉驗(yàn)證與柵格搜索方法相結(jié)合,以均方誤差的均值最小為目標(biāo)對參數(shù)c與g進(jìn)行優(yōu)選,以提高參數(shù)優(yōu)選的效率和準(zhǔn)確性,消除樣本數(shù)據(jù)的隨機(jī)誤差對平臺回歸模型精確性的影響。

    采用柵格搜索與交叉驗(yàn)證相結(jié)合的方法對參數(shù)(c,g)的尋優(yōu)過程如圖4所示,具體步驟如下。

    圖4 參數(shù)(c,g)的尋優(yōu)流程

    (15)

    (2)設(shè)置(c,g)的搜索范圍和步長,由于初始無法估計(jì)c與g的大致范圍,搜索范圍應(yīng)盡可能大,所以采用對數(shù)坐標(biāo)形式(log2c,log2g)對其進(jìn)行柵格化。

    (3)將訓(xùn)練數(shù)據(jù)等分為K組(K≥2),對于柵格中每一組(c,g),以任意一組作為測試組,其余K-1組作為訓(xùn)練組,訓(xùn)練模型后對測試組進(jìn)行預(yù)測,并計(jì)算這一組(c,g)對應(yīng)的均方誤差,其定義為

    (16)

    式中,f(xi)為預(yù)測值;yi為樣本實(shí)際值。

    (4)對于與步驟(3)相同的(c,g),更換另一測試組,取剩余K-1組作為訓(xùn)練集,根據(jù)式(16)計(jì)算均方誤差,直至對K組都進(jìn)行一次預(yù)測后,取K組均方誤差的均值ΔMSE作為該組(c,g)的預(yù)測誤差。

    (5)更換參數(shù)組合(c,g),重復(fù)步驟(3)和(4),逐次計(jì)算網(wǎng)格中各參數(shù)組合(c,g)下的均方誤差,對所有的組合的ΔMSE進(jìn)行比較,ΔMSE最小的參數(shù)組合即網(wǎng)格區(qū)間內(nèi)的最佳(c,g)值。

    (6)判斷ΔMSE是否小于Δmin,若是則結(jié)束流程;否則根據(jù)ΔMSE的分布情況,縮小參數(shù)搜索范圍和步長,重復(fù)步驟(3)~(5)直至滿足要求。

    3 實(shí)驗(yàn)驗(yàn)證

    3.1 實(shí)驗(yàn)裝置

    圖5 實(shí)驗(yàn)裝置

    圖6 實(shí)驗(yàn)數(shù)據(jù)采集過程

    如圖5所示,實(shí)驗(yàn)裝置由微操作平臺、壓電陶瓷驅(qū)動器、驅(qū)動電源、位移傳感器、信號調(diào)理器、D/A卡、工控機(jī)等組成。壓電驅(qū)動器為哈爾濱芯明天公司生產(chǎn)的80VS12,電壓取值范圍為0~150 V;驅(qū)動電源為芯明天公司生產(chǎn)的XE-501,功率放大倍數(shù)為15;位移傳感器為PI公司生產(chǎn)的D-E20.200,分辨率為0.2 μm,線性度為0.3 %;D/A卡為NI公司生產(chǎn)的PCI-6713。實(shí)驗(yàn)樣本數(shù)據(jù)點(diǎn)采集過程如圖6所示,根據(jù)第1節(jié)所述的施加電壓順序,采用LabView軟件編寫電壓輸入程序,通過D/A卡將工控機(jī)輸出的電壓經(jīng)功率放大作為壓電驅(qū)動器的輸入電壓,采用位移傳感器測量所對應(yīng)的輸出位移,位移信號經(jīng)信號調(diào)理器通過USB接口輸入至工控機(jī),采用LabView編寫位移采集程序,以獲得對應(yīng)各個電壓等分點(diǎn)的輸出位移。

    3.2 微操作平臺的非線性模型的實(shí)驗(yàn)驗(yàn)證

    由于合適的參數(shù)(c,g)對于能否建立精確反映微操作平臺非線性模型非常關(guān)鍵,下面通過擬合平臺外環(huán)曲線來說明柵格搜索和交叉驗(yàn)證方法找尋(c,g)值的有效性。由于VS80型壓電驅(qū)動器的輸入電壓超過120 V會影響其使用壽命,故將驅(qū)動器的飽和電壓設(shè)為120 V。為了擬合平臺的外環(huán)曲線,將飽和電壓等分為40個點(diǎn),并測量該等分點(diǎn)對應(yīng)的輸出位移得到樣本數(shù)據(jù)點(diǎn),見表1。

    表1 用于擬合微操作平臺外環(huán)曲線的樣本點(diǎn)

    c與g的初始尋優(yōu)范圍為[2-10,210],尋優(yōu)步長設(shè)為2-1,K=3,ε=0.001,Δmin=10-4μm2。圖7所示為c和g尋優(yōu)范圍內(nèi)每個柵格對應(yīng)的ΔMSE值,圖8所示為對應(yīng)的ΔMSE值等高線。由圖7和圖8可看出ΔMSE的變化趨勢,ΔMSE接近于0的c與g的范圍如圖8中虛線部分所示,在該范圍內(nèi)ΔMSE值均小于Δmin,在該范圍內(nèi)尋優(yōu)選取最小的c值所對應(yīng)的一組(c,g),可得到c=32,g=1.4,ΔMSE=2.3625×10-6μm2。根據(jù)該(c,g)值可擬合出描述平臺外環(huán)的回歸曲線,如圖9所示。

    圖7 均方差均值的分布

    圖8 均方差均值的等高線

    圖9 微操作平臺外環(huán)的擬合曲線

    為了說明回歸曲線的精確性,選取12個測試樣本點(diǎn)進(jìn)行驗(yàn)證分析,測試樣本點(diǎn)包括外環(huán)上升曲線的6組點(diǎn)與下降曲線的6組點(diǎn),見表2。

    圖10所示為12個測試點(diǎn)的實(shí)測位移與由模型得到的預(yù)測位移之間的誤差,由圖10可知,誤差值范圍為[-0.15,0.15]μm。同時(shí),將12個測試點(diǎn)的實(shí)測位移與預(yù)測位移分別代入式(14)和式(16),可分別得到衡量回歸曲線精度的指標(biāo)r2=0.999 99和EMS=9.504×10-7μm2,由該指標(biāo)可知平臺的外環(huán)曲線擬合度較好,說明交叉驗(yàn)證的柵格搜索方法能確定合適的(c,g)值以建立精確反映平臺非線性模型。

    表2 測試點(diǎn)和對應(yīng)的位移預(yù)測值與實(shí)測值

    圖10 誤差分布

    下面通過實(shí)驗(yàn)驗(yàn)證分析采用支持向量機(jī)和Preisach模型描述微操作平臺非線性的有效性。將壓電驅(qū)動器的輸入電壓分成16等分,得到等分點(diǎn)上的極值點(diǎn)電壓(α,β)和平臺對應(yīng)的輸出位移x(α,β)構(gòu)成的樣本點(diǎn),見表3。設(shè)c與g的尋優(yōu)范圍為[2-10,210],步長設(shè)為2-1,K = 4,ε= 0.001,Δmin=10-5μm2。根據(jù)圖3所示流程可得到最優(yōu)(c,g)值組合,即c = 32,g = 4,該組合對應(yīng)的均方誤差的均值ΔMSE=1.485×10-6μm2,說明(c,g)值能滿足設(shè)定的要求。根據(jù)該參數(shù)組合采用支持向量機(jī)對表3中的樣本點(diǎn)進(jìn)行訓(xùn)練,得到圖11所示的回歸曲線,該曲線反映了平臺的非線性特性。根據(jù)式(14)可計(jì)算出反映曲線擬合精度的指標(biāo)r2=0.999 986,說明所擬合的曲線能精確地反映平臺的非線性特性。

    圖11 微操作平臺的擬合曲線

    為了驗(yàn)證分析運(yùn)用Preisach模型和支持向量機(jī)所建立的模型能預(yù)測不同輸入電壓序列的輸出位移,可通過比較分析當(dāng)壓電驅(qū)動器在不同輸入電壓序列作用下的平臺輸出位移的實(shí)際值與預(yù)測值來實(shí)現(xiàn)。選取第1測試組為極值點(diǎn)位于等分點(diǎn)上的電壓序列0-120 V-60 V-75 V-15 V-105 V,第2測試組為極值點(diǎn)不在等分點(diǎn)上的電壓序列0-112.5 V-19.5 V-69 V-36 V-93 V-68 V。由式(1)和式(2)表示的Preisach模型可知,要計(jì)算出電壓序列的輸出位移,需要得到不同輸入電壓(α,β)所對應(yīng)的位移x(α,β),可利用圖11所示的回歸模型得到x(α,β),第1和第2測試組各輸入電壓對應(yīng)的位移x(α,β)分別見表4和表5。

    表3 用于擬合微操作平臺遲滯環(huán)的樣本點(diǎn) μm

    表4 第1測試組電壓對應(yīng)的x(α,β)預(yù)測值

    表5 第2測試組電壓對應(yīng)的x(α,β)預(yù)測值

    根據(jù)表4和表5中的x(α,β),將其代入式(1)和式(2)表示的Preisach模型,則可以分別計(jì)算出對應(yīng)于兩組測試組各電壓序列的輸出位移f(z)的預(yù)測值。根據(jù)圖5所示的實(shí)驗(yàn)過程,對壓電陶瓷驅(qū)動器分別施加上述兩組序列電壓,可以分別得到不同輸入電壓所對應(yīng)的輸出位移的實(shí)測值。第1和第2測試組的各輸入電壓對應(yīng)的輸出位移的實(shí)測值、預(yù)測值、相對誤差分別見表6和表7,可知,第1和第2測試組的相對誤差范圍分別為0.6%~2.1%和0.02%~2.1%,說明預(yù)測位移與實(shí)測位移非常接近,也就是說,聯(lián)合Preisach模型和支持向量機(jī)回歸模型能精確預(yù)測平臺在驅(qū)動器輸入任意電壓序列所對應(yīng)的輸出位移,從而實(shí)時(shí)預(yù)測平臺的輸出位移。

    表6 第1測試組的實(shí)測位移和預(yù)測位移

    表7 第2測試組的實(shí)測位移和預(yù)測位移

    4 結(jié)語

    為了精確描述壓電微操作平臺非線性特性,本文提出了一種Preisach模型與支持向量機(jī)回歸方法相結(jié)合的建模方法。提出了一種基于交叉驗(yàn)證和柵格化搜索方法,對非線性模型精度有較大影響的懲罰系數(shù)c與核函數(shù)參數(shù)g進(jìn)行尋優(yōu)。以一種一維微操作平臺為對象,利用支持向量機(jī)回歸理論擬合的非線性回歸曲線,并結(jié)合傳統(tǒng)的Preisach模型,可對任意輸入電壓系列的平臺輸出位移進(jìn)行精確的預(yù)測。實(shí)驗(yàn)驗(yàn)證分析結(jié)果表明,基于交叉驗(yàn)證的柵格搜索方法能確定合適的(c,g)值以建立反映平臺非線性的精確回歸模型,聯(lián)合Preisach模型能精確預(yù)測任意輸入電壓序列的輸出位移,說明所建立的模型能精確描述微操作平臺的非線性特性,以實(shí)現(xiàn)其精確運(yùn)動。

    [1]YueYi,GaoFeng,ZhaoXianchao,etal.RelationshipamongInput-force,Payload,StiffnessandDisplacementofa3-DOFPerpendicularParallelMicro-Manipulator[J].MechanismandMachineTheory, 2010, 45(5):756-771.

    [2] 胡俊峰, 徐貴陽, 郝亞洲.基于響應(yīng)面法的微操作平臺多目標(biāo)優(yōu)化[J].光學(xué)精密工程, 2015, 23(4):1096-1104.HuJunfeng,XuGuiyang,HaoYazhou.Multi-objectiveOptimizationofaNovelMicro-manipulationStageBasedonResponseSurfaceMethod[J].OpticsandPrecisionEngineering, 2015, 23(4): 1096-1104.

    [3]VominhT,TjahjowidodoT,RamonH,etal.ANewApproachtoModelingHysteresisinaPneumaticArtificialMuscleUsingtheMaxwell-slipModel[J].IEEE/ASMETransactionsonMechatronics, 2011, 16(1): 177-186.

    [4]LinCJ,LinPT.TrackingControlofaBiaxialPiezo-actuatedPositioningStageUsingGeneralizedDuhemModel[J].Computers&MathematicswithApplications, 2012, 64(5): 766-787.

    [5]JiangH,JiH,QiuJ,etal.AModifiedPrandtl-ishlinskiiModelforModelingAsymmetricHysteresisofPiezoelectricActuators[J].IEEETransactionsonUltrasonicsFerroelectrics&FrequencyControl, 2010, 57(5): 1200-1210.

    [6]XiaoShunli,LiYangmin.ModelingandHighDynamicCompensatingtheRate-dependentHysteresisofPiezoelectricActuatorsviaaNovelModifiedInversePreisachModel[J].IEEETransactionsonControlSystemsTechnology, 2013, 21(5): 1549-1557.

    [7]MohammadRZ,MohsenF,HassanS,etal.HysteresisNonlinearityIdentificationUsingNewPreisachModel-basedArtificialNeuralNetworkApproach[J].JournalofAppliedMathematics, 2011, 2011(22):1-21.

    [8] 劉向東, 修春波, 李黎,等.遲滯非線性系統(tǒng)的神經(jīng)網(wǎng)絡(luò)建模[J].壓電與聲光, 2007, 29(1): 106-108.LiuXiangdong,XiuChunbo,LiLi,etal.HysteresisModelingUsingNeuralNetworks[J].Piezoelectectrics&Acoustooptics, 2007, 29(1): 106-108.

    [9] 趙新龍, 譚永紅.對Preisach類的遲滯非線性神經(jīng)網(wǎng)絡(luò)建模[J].控制理論與應(yīng)用, 2006, 23(4): 581-585.ZhouXinlong,TanYonghong.ModelingPreisach-typeHysteresisNonlinearityUsingNeuralNetworks[J].ControlTheory&Applications, 2006, 23(4): 581-585.

    [10] 張建華, 葛紅宇, 李宏勝,等.壓電致動器Preisach模型的模糊插值算法及實(shí)現(xiàn)[J].南京理工大學(xué)學(xué)報(bào),2011, 35(6): 780-785.ZhangJianhua,GeHongyu,LiHongsheng,etal.FuzzyInterpolationAlgorithmandRealizationofPreisachModelforPiezoActuator[J].JournalofNanjingUniversityofScienceandTechnology,2011, 35(6):780-785.

    [11] 賴志林, 劉向東, 耿潔,等.壓電陶瓷執(zhí)行器遲滯的滑模逆補(bǔ)償控制[J].光學(xué)精密工程,2011,19(6):1281-1290.LaiZhilin,LiuXiangdong,GengJie,etal.SlidingModeControlofHysteresisofPiezoceramicActuatorBasedonInversePreisachCompensation[J].OpticsandPrecisionEngineering, 2011, 19(6): 1281-1290.

    [12] 魏燕定, 陶惠峰.壓電驅(qū)動器遲滯特性的Preisach模型研究[J].壓電與聲光, 2004, 26(5): 364-367.WeiYanding,TaoHuifeng.StudythePreisachModelofHysteresisinPiezoelectricActuator[J].Piezoelectrics&Acoustooptics, 2004, 26(5): 364-367.

    [13] 李黎, 劉向東,王偉,等.壓電陶瓷執(zhí)行器遲滯特性的廣義非線性Preisach模型及其數(shù)值實(shí)現(xiàn)[J].光學(xué)精密工程, 2007, 15(5): 706-712.LiLi,LiuXiangdong,WangWei,etal.GeneralizedNonlinearPresachModelforHysteresisNonlinearityofPiezoceramicActuatorandItsNumericalImplementation[J].OpticsandPrecisionEngineering,2007,15(5):706-712.

    [14] 陳輝, 譚永紅, 周杏鵬,等.壓電陶瓷執(zhí)行器的動態(tài)模型辨識與控制[J].光學(xué)精密工程, 2012, 20(1):88-95.ChenHui,TanYonghong,ZhouXingpeng,etal.IdentificationandControlofDynamicModelingforPiezoceramicActuator[J].OpticsandPrecisionEngineering, 2012, 20(1): 88-95.

    [15]SmolaAJ,Sch?lkopfB.ATutorialonSupportVectorRegression[J].StatisticsandComputing, 2004, 14(3): 199-222.

    [16] 韓彥彬, 白廣忱, 李曉穎,等.基于支持向量機(jī)柔性機(jī)構(gòu)動態(tài)可靠性分析[J].機(jī)械工程學(xué)報(bào),2014,50(11):86-92.HanYanbin,BaiGuangchen,LiXiaoying,etal.DynamicReliabilityAnalysisofFlexibleMechanismBasedonSupportVectorMachine[J].JournalofMechanicalEngineering, 2014, 50(11): 86-92.

    (編輯 陳 勇)

    Description of Nonlinear Characteristics of Piezoelectric Micro-manipulation Stage Based on SVM

    Hu Junfeng Zheng Changhu Cai Jianyang

    Jiangxi University of Science and Technology,Ganzhou,Jiangxi,341000

    The nonlinear characteristics of micro-manipulation stage would be occurred due to the hysteresis of piezoelectric actuator, which might seriously affect the movement precision and repeated positioning accuracy of the stage. In order to solve the issue, a hybrid modeling method was proposed to describe the hysteresis characteristics based on Preisach and SVM. As for a one dimensional micro-manipulation stage, the sample points were chosen as the piezoelectric driving voltages and corresponding output displacements of the stage. The regression model of the hysteresis loop of the PEA was established by using SVM theory. The model might be used to predict the displacements corresponding to the voltages not at sample points. Combined with the Preisach model, the output displacement corresponding to any voltage sequence might be obtained. The experimental studies were carried out in order to verify the validity of the presented model. Arbitrarily selecting two sets of input voltage sequences, and the prediction value of the output displacements corresponding to the selected voltage sequences might be obtained by using the regression model and Preisach model. The measured values were obtained under the same voltage sequences by experimental tests. Comparing the measured values and predicted values, the results show that the relative error ranges of the two sets of the measured and predicted values are as 0.6%~2.1% and 0.02%~2.1%, respectively. Therefore, the predicted displacements are very close to the measured ones. It illustrates that the proposed model may accurately describe the nonlinear characteristics of the micro-manipulation stage in order to realize its precise movements.

    micro-manipulation stage; piezoelectric actuator; Preisach model; support vector machine(SVM); nonlinear characteristics

    2016-01-11

    國家自然科學(xué)基金資助項(xiàng)目(51265016,51565016)

    TH703;TP274

    10.3969/j.issn.1004-132X.2016.22.006

    胡俊峰,男,1978年生。江西理工大學(xué)機(jī)電工程學(xué)院副教授。主要研究方向?yàn)槿犴槞C(jī)構(gòu)及智能控制。發(fā)表論文30余篇。鄭昌虎,男,1991年生。江西理工大學(xué)機(jī)電工程學(xué)院碩士研究生。蔡建陽,男,1991年生。江西理工大學(xué)機(jī)電工程學(xué)院碩士研究生。

    猜你喜歡
    分點(diǎn)驅(qū)動器柵格
    基于鄰域柵格篩選的點(diǎn)云邊緣點(diǎn)提取方法*
    藏起驅(qū)動器號確保數(shù)據(jù)安全
    電腦愛好者(2020年6期)2020-05-26 09:27:33
    來自低谷的你
    青年生活(2020年13期)2020-05-26 01:51:33
    定比分點(diǎn)之換底分點(diǎn)伸縮法
    五禽戲“動作節(jié)分點(diǎn)”劃分與學(xué)練建議(三)
    健身氣功(2018年2期)2018-06-04 06:51:02
    壓電陶瓷驅(qū)動器的遲滯特性
    不同剖面形狀的柵格壁對柵格翼氣動特性的影響
    SiC基和Si基永磁同步電動機(jī)驅(qū)動器的比較
    基于CVT排布的非周期柵格密度加權(quán)陣設(shè)計(jì)
    新型低電壓大變形微驅(qū)動器數(shù)值求解及仿真
    成人国产综合亚洲| 久久久久久大精品| 亚洲国产欧美人成| 一级作爱视频免费观看| 男女做爰动态图高潮gif福利片| 欧美日韩瑟瑟在线播放| 国产欧美日韩精品一区二区| 特大巨黑吊av在线直播| 久久久久九九精品影院| 亚洲aⅴ乱码一区二区在线播放| 国产野战对白在线观看| 小蜜桃在线观看免费完整版高清| 啪啪无遮挡十八禁网站| 国产精品一区二区性色av| 蜜桃久久精品国产亚洲av| 国产精品日韩av在线免费观看| 免费在线观看影片大全网站| av专区在线播放| 欧美黄色淫秽网站| 午夜福利在线观看吧| 天堂√8在线中文| 久久人人精品亚洲av| 欧美最新免费一区二区三区 | 国模一区二区三区四区视频| www.999成人在线观看| 日本黄大片高清| 国产精品永久免费网站| 久久久久久久精品吃奶| 女生性感内裤真人,穿戴方法视频| 一二三四社区在线视频社区8| 亚洲综合色惰| 午夜a级毛片| 久久99热这里只有精品18| 好男人电影高清在线观看| 国产 一区 欧美 日韩| 美女cb高潮喷水在线观看| 免费av毛片视频| 午夜福利在线在线| 久久久久国产精品人妻aⅴ院| 国产aⅴ精品一区二区三区波| 欧美丝袜亚洲另类 | 麻豆国产97在线/欧美| 婷婷六月久久综合丁香| 亚洲精品日韩av片在线观看| 国产av在哪里看| 欧美区成人在线视频| 亚洲欧美精品综合久久99| or卡值多少钱| 神马国产精品三级电影在线观看| 国产单亲对白刺激| 搡老妇女老女人老熟妇| 欧美在线黄色| 国产高潮美女av| 中文资源天堂在线| 99久久99久久久精品蜜桃| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 在线看三级毛片| 美女高潮的动态| 高清毛片免费观看视频网站| 亚洲精品粉嫩美女一区| 成人三级黄色视频| 国产又黄又爽又无遮挡在线| 97超视频在线观看视频| 精品国产亚洲在线| 美女被艹到高潮喷水动态| 一级av片app| 成年人黄色毛片网站| 深爱激情五月婷婷| 亚州av有码| 免费看a级黄色片| 久久精品91蜜桃| 一级黄色大片毛片| 天堂网av新在线| 观看免费一级毛片| 99久国产av精品| 日本 欧美在线| 国产av不卡久久| 日韩欧美 国产精品| 三级男女做爰猛烈吃奶摸视频| 黄色一级大片看看| 白带黄色成豆腐渣| 欧美最新免费一区二区三区 | 在线看三级毛片| 欧洲精品卡2卡3卡4卡5卡区| 免费看美女性在线毛片视频| 国产精品自产拍在线观看55亚洲| 一二三四社区在线视频社区8| 久久国产乱子免费精品| 好男人电影高清在线观看| 国产三级中文精品| 一区二区三区高清视频在线| 成人国产一区最新在线观看| 亚洲精品一区av在线观看| 日韩成人在线观看一区二区三区| 村上凉子中文字幕在线| 国产精品久久久久久人妻精品电影| 欧美性猛交╳xxx乱大交人| 午夜福利18| 黄色一级大片看看| 国产精品人妻久久久久久| 内地一区二区视频在线| 18美女黄网站色大片免费观看| 超碰av人人做人人爽久久| 亚洲国产欧洲综合997久久,| 美女大奶头视频| 中文字幕人成人乱码亚洲影| 男女下面进入的视频免费午夜| 久久精品国产自在天天线| 男女做爰动态图高潮gif福利片| 午夜a级毛片| 亚洲精品一区av在线观看| 成人特级黄色片久久久久久久| 男人和女人高潮做爰伦理| 男女床上黄色一级片免费看| 看十八女毛片水多多多| 免费观看精品视频网站| 欧美国产日韩亚洲一区| 久久久国产成人精品二区| 一级黄片播放器| 日韩 亚洲 欧美在线| 男插女下体视频免费在线播放| 男人和女人高潮做爰伦理| 毛片一级片免费看久久久久 | 国产精品野战在线观看| 国产综合懂色| 亚洲在线观看片| 亚洲在线观看片| 中文字幕人妻熟人妻熟丝袜美| www.熟女人妻精品国产| 国内久久婷婷六月综合欲色啪| 日韩欧美国产在线观看| 成人欧美大片| 又爽又黄a免费视频| 亚洲成a人片在线一区二区| 午夜福利视频1000在线观看| 日日干狠狠操夜夜爽| 免费看美女性在线毛片视频| 亚洲无线在线观看| 国产成年人精品一区二区| 国产三级黄色录像| 免费大片18禁| 免费搜索国产男女视频| 久久久色成人| 51国产日韩欧美| 老司机福利观看| 午夜久久久久精精品| 国产伦精品一区二区三区四那| 日韩精品中文字幕看吧| av福利片在线观看| 国产黄片美女视频| 99久久精品一区二区三区| 99热这里只有是精品50| 99热这里只有是精品在线观看 | 天堂网av新在线| 99热这里只有是精品在线观看 | 91狼人影院| 日韩欧美国产在线观看| 亚洲一区高清亚洲精品| 三级毛片av免费| 久久久国产成人精品二区| 色综合欧美亚洲国产小说| 中文字幕久久专区| 在线观看舔阴道视频| 757午夜福利合集在线观看| 1000部很黄的大片| 亚洲成人久久性| 欧美高清成人免费视频www| 嫩草影院新地址| 波多野结衣高清作品| 久久国产乱子免费精品| 亚洲精品一区av在线观看| 永久网站在线| 国内精品美女久久久久久| 免费在线观看影片大全网站| 悠悠久久av| 啦啦啦观看免费观看视频高清| 国产精品久久久久久久久免 | 可以在线观看的亚洲视频| 最近中文字幕高清免费大全6 | 国产在线男女| 欧美+日韩+精品| 一级黄色大片毛片| 欧美又色又爽又黄视频| 一a级毛片在线观看| 琪琪午夜伦伦电影理论片6080| 欧美性感艳星| 美女高潮喷水抽搐中文字幕| 亚洲片人在线观看| 久久国产乱子免费精品| 天堂动漫精品| 免费大片18禁| 国产精品精品国产色婷婷| 一个人免费在线观看电影| 日本一二三区视频观看| 国内精品久久久久久久电影| 欧洲精品卡2卡3卡4卡5卡区| 免费人成在线观看视频色| 免费看光身美女| 精品国产亚洲在线| 亚洲国产精品sss在线观看| 97超视频在线观看视频| 午夜亚洲福利在线播放| 亚洲av美国av| 成人特级av手机在线观看| 在线免费观看的www视频| 国产不卡一卡二| 九九久久精品国产亚洲av麻豆| 精品久久久久久,| 在线观看av片永久免费下载| 身体一侧抽搐| 在线观看av片永久免费下载| 人妻夜夜爽99麻豆av| 搡老岳熟女国产| 国产成年人精品一区二区| 久久6这里有精品| 国产成人av教育| 国产aⅴ精品一区二区三区波| 国产免费男女视频| 精品人妻偷拍中文字幕| 国产一区二区在线观看日韩| 99国产精品一区二区三区| 国产久久久一区二区三区| 欧美色欧美亚洲另类二区| 可以在线观看的亚洲视频| 国产中年淑女户外野战色| 蜜桃亚洲精品一区二区三区| 内地一区二区视频在线| 啦啦啦观看免费观看视频高清| 在线观看舔阴道视频| av天堂中文字幕网| 日韩欧美 国产精品| 精品久久久久久久久亚洲 | 国内少妇人妻偷人精品xxx网站| 少妇丰满av| 高清日韩中文字幕在线| 亚洲av免费在线观看| 亚洲七黄色美女视频| 丝袜美腿在线中文| 女人被狂操c到高潮| 欧美一级a爱片免费观看看| 国产精品1区2区在线观看.| 色在线成人网| 99在线人妻在线中文字幕| 中国美女看黄片| 成年女人毛片免费观看观看9| 小蜜桃在线观看免费完整版高清| 亚洲经典国产精华液单 | 精品午夜福利在线看| 亚洲七黄色美女视频| 毛片女人毛片| 国产v大片淫在线免费观看| 国产高清三级在线| 成人精品一区二区免费| 精品国产亚洲在线| av欧美777| 性色av乱码一区二区三区2| 99久久精品国产亚洲精品| 在线播放国产精品三级| 久久久久国产精品人妻aⅴ院| 99国产极品粉嫩在线观看| 九色成人免费人妻av| 久久午夜福利片| 人妻久久中文字幕网| 午夜视频国产福利| 国产成人a区在线观看| 久久久久久久久中文| 一本综合久久免费| 一个人看视频在线观看www免费| 亚洲成人免费电影在线观看| 精品一区二区三区视频在线观看免费| 91久久精品电影网| 日韩欧美精品v在线| 国产一区二区激情短视频| 嫩草影院精品99| 亚洲精品久久国产高清桃花| 欧美极品一区二区三区四区| 国产亚洲精品av在线| 97碰自拍视频| 欧美成人一区二区免费高清观看| 简卡轻食公司| 99热精品在线国产| 免费在线观看影片大全网站| 久久99热这里只有精品18| av视频在线观看入口| 欧美黄色淫秽网站| 精品久久久久久久久久久久久| 亚洲avbb在线观看| 国产精品电影一区二区三区| 99精品在免费线老司机午夜| 国产探花极品一区二区| 亚洲一区二区三区色噜噜| 日韩中字成人| 国产一区二区激情短视频| 国产一区二区三区在线臀色熟女| 午夜福利视频1000在线观看| 亚洲avbb在线观看| 免费在线观看日本一区| 国内精品一区二区在线观看| 嫩草影视91久久| 一二三四社区在线视频社区8| 天天躁日日操中文字幕| 欧美3d第一页| 可以在线观看毛片的网站| 少妇高潮的动态图| 国产精品综合久久久久久久免费| 久久久久久久久大av| 男插女下体视频免费在线播放| 最近最新中文字幕大全电影3| 日韩中文字幕欧美一区二区| 我的女老师完整版在线观看| www.www免费av| 老女人水多毛片| 老司机深夜福利视频在线观看| 精品一区二区三区视频在线观看免费| 欧美又色又爽又黄视频| 国产欧美日韩一区二区三| 欧美成人性av电影在线观看| 午夜福利成人在线免费观看| 白带黄色成豆腐渣| 深夜a级毛片| 免费电影在线观看免费观看| 成人无遮挡网站| 丰满人妻熟妇乱又伦精品不卡| 在线观看66精品国产| 国产私拍福利视频在线观看| 18禁裸乳无遮挡免费网站照片| 又黄又爽又刺激的免费视频.| 亚洲精品在线美女| 亚洲av免费在线观看| 国产亚洲精品久久久com| 男女那种视频在线观看| 人妻夜夜爽99麻豆av| 五月玫瑰六月丁香| 国产精品爽爽va在线观看网站| 三级毛片av免费| 久久精品国产亚洲av天美| 白带黄色成豆腐渣| 久久久久久九九精品二区国产| 免费av毛片视频| 乱人视频在线观看| 成人国产一区最新在线观看| 久久精品影院6| 高清毛片免费观看视频网站| 精品99又大又爽又粗少妇毛片 | 亚洲精品456在线播放app | 午夜免费男女啪啪视频观看 | 亚洲人成伊人成综合网2020| 淫妇啪啪啪对白视频| 欧美三级亚洲精品| 五月玫瑰六月丁香| 69人妻影院| 老司机午夜十八禁免费视频| 午夜免费男女啪啪视频观看 | 少妇的逼好多水| 男插女下体视频免费在线播放| 久久久久免费精品人妻一区二区| 午夜免费成人在线视频| 国产成人福利小说| 俺也久久电影网| 黄色丝袜av网址大全| 婷婷亚洲欧美| avwww免费| 99久久无色码亚洲精品果冻| 亚洲第一欧美日韩一区二区三区| 99国产精品一区二区三区| 亚洲国产精品sss在线观看| 免费观看人在逋| 一卡2卡三卡四卡精品乱码亚洲| 国产精品乱码一区二三区的特点| 美女大奶头视频| 天美传媒精品一区二区| 午夜福利成人在线免费观看| 国语自产精品视频在线第100页| 国产毛片a区久久久久| 91av网一区二区| 99在线人妻在线中文字幕| 免费观看的影片在线观看| 91麻豆av在线| 亚洲成人免费电影在线观看| 丰满人妻一区二区三区视频av| 99在线视频只有这里精品首页| 午夜日韩欧美国产| 久久人人精品亚洲av| 毛片一级片免费看久久久久 | 亚洲精品色激情综合| 嫩草影院精品99| 久9热在线精品视频| 国产综合懂色| 亚洲专区中文字幕在线| or卡值多少钱| 美女高潮的动态| 国产亚洲精品综合一区在线观看| 中文字幕av成人在线电影| 听说在线观看完整版免费高清| 欧美成人性av电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区三区视频了| 亚洲无线在线观看| 欧美成人免费av一区二区三区| 成年版毛片免费区| 国产精品一及| 精品人妻熟女av久视频| 一个人看视频在线观看www免费| 88av欧美| 嫁个100分男人电影在线观看| 国产精品久久电影中文字幕| 九九在线视频观看精品| 亚洲美女视频黄频| 午夜激情欧美在线| 99精品在免费线老司机午夜| 国内少妇人妻偷人精品xxx网站| 哪里可以看免费的av片| 男人和女人高潮做爰伦理| 看免费av毛片| 亚洲国产高清在线一区二区三| 亚洲av免费在线观看| 免费在线观看成人毛片| 好看av亚洲va欧美ⅴa在| 午夜福利欧美成人| 亚洲欧美清纯卡通| 日日摸夜夜添夜夜添av毛片 | 欧美+亚洲+日韩+国产| 午夜亚洲福利在线播放| 国产在视频线在精品| 国产精品久久久久久亚洲av鲁大| 老鸭窝网址在线观看| 九色国产91popny在线| 亚洲乱码一区二区免费版| 搡老妇女老女人老熟妇| 亚洲在线观看片| h日本视频在线播放| 99热只有精品国产| 18美女黄网站色大片免费观看| 日本一本二区三区精品| 亚洲一区二区三区色噜噜| 国产精品日韩av在线免费观看| 首页视频小说图片口味搜索| 亚洲av二区三区四区| 在线免费观看的www视频| 日韩 亚洲 欧美在线| 国产精品乱码一区二三区的特点| 他把我摸到了高潮在线观看| 成人av在线播放网站| 亚洲一区二区三区不卡视频| 久久99热这里只有精品18| 日韩亚洲欧美综合| 91狼人影院| 黄片小视频在线播放| 麻豆成人午夜福利视频| 在线观看舔阴道视频| 亚洲美女视频黄频| 搡老妇女老女人老熟妇| 亚洲美女搞黄在线观看 | 久久久久亚洲av毛片大全| 亚洲最大成人av| 婷婷丁香在线五月| 久久久久免费精品人妻一区二区| 久久精品夜夜夜夜夜久久蜜豆| 久9热在线精品视频| 黄色配什么色好看| 99在线人妻在线中文字幕| 蜜桃亚洲精品一区二区三区| 欧美三级亚洲精品| 少妇裸体淫交视频免费看高清| 高潮久久久久久久久久久不卡| 国产三级在线视频| 高清毛片免费观看视频网站| 一本综合久久免费| 亚洲av免费在线观看| 亚洲人成网站在线播| 欧美最新免费一区二区三区 | 色综合站精品国产| 国产69精品久久久久777片| 免费在线观看影片大全网站| 日韩欧美精品免费久久 | 亚洲最大成人av| 男人狂女人下面高潮的视频| 精品福利观看| 九九久久精品国产亚洲av麻豆| 成年人黄色毛片网站| 久久久久久久久中文| 精品人妻视频免费看| 18禁裸乳无遮挡免费网站照片| 亚洲成av人片在线播放无| 亚洲av美国av| 亚洲美女搞黄在线观看 | 在现免费观看毛片| 亚洲国产日韩欧美精品在线观看| 欧美bdsm另类| aaaaa片日本免费| 国产精品98久久久久久宅男小说| 乱人视频在线观看| 国产综合懂色| 日本五十路高清| 最后的刺客免费高清国语| 人人妻人人看人人澡| 亚洲国产欧洲综合997久久,| 亚洲av不卡在线观看| 身体一侧抽搐| 一区二区三区免费毛片| 欧美日韩国产亚洲二区| 九色成人免费人妻av| www.熟女人妻精品国产| 亚洲精品在线观看二区| 嫁个100分男人电影在线观看| 熟妇人妻久久中文字幕3abv| 亚洲av美国av| 特大巨黑吊av在线直播| 久久久久久九九精品二区国产| 日本 欧美在线| 中文资源天堂在线| 天堂动漫精品| 日本在线视频免费播放| 村上凉子中文字幕在线| 国产欧美日韩一区二区精品| 国产乱人伦免费视频| 俄罗斯特黄特色一大片| 观看免费一级毛片| 日韩成人在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 日本三级黄在线观看| 国产免费一级a男人的天堂| 日本a在线网址| 国产高清三级在线| 在线免费观看的www视频| 久久久久国内视频| 黄色配什么色好看| 在线观看美女被高潮喷水网站 | 亚洲午夜理论影院| 天堂动漫精品| 高清在线国产一区| 在线观看av片永久免费下载| 亚洲五月婷婷丁香| 国产视频一区二区在线看| 久久午夜福利片| 免费无遮挡裸体视频| 国产精品女同一区二区软件 | 九色国产91popny在线| 九九久久精品国产亚洲av麻豆| 欧美高清性xxxxhd video| 欧美在线黄色| 中文字幕精品亚洲无线码一区| 不卡一级毛片| 色播亚洲综合网| 国产伦在线观看视频一区| 麻豆国产97在线/欧美| 高潮久久久久久久久久久不卡| 黄色女人牲交| 国产成人啪精品午夜网站| www日本黄色视频网| 性欧美人与动物交配| 亚洲精品一区av在线观看| 18禁黄网站禁片午夜丰满| 国产一区二区三区视频了| 午夜免费激情av| 女同久久另类99精品国产91| 嫩草影院新地址| 老司机深夜福利视频在线观看| av在线天堂中文字幕| 18+在线观看网站| 女生性感内裤真人,穿戴方法视频| 国产在线精品亚洲第一网站| 在线观看一区二区三区| 如何舔出高潮| 色在线成人网| 成人鲁丝片一二三区免费| 村上凉子中文字幕在线| 亚洲成a人片在线一区二区| 男人舔奶头视频| 熟女人妻精品中文字幕| 国产视频一区二区在线看| 两个人视频免费观看高清| 婷婷亚洲欧美| 深爱激情五月婷婷| a级一级毛片免费在线观看| 国产精品日韩av在线免费观看| 十八禁国产超污无遮挡网站| 搡女人真爽免费视频火全软件 | 高潮久久久久久久久久久不卡| 精品一区二区免费观看| 校园春色视频在线观看| 免费看a级黄色片| 国产亚洲欧美98| 自拍偷自拍亚洲精品老妇| 亚洲成a人片在线一区二区| 亚洲精品一区av在线观看| 美女高潮的动态| 日韩成人在线观看一区二区三区| 亚洲第一欧美日韩一区二区三区| 国产野战对白在线观看| 国产单亲对白刺激| 美女免费视频网站| 精品午夜福利视频在线观看一区| 美女大奶头视频| 欧美极品一区二区三区四区| 少妇丰满av| 丰满人妻熟妇乱又伦精品不卡| 99久久精品热视频| 真实男女啪啪啪动态图| 亚洲av免费在线观看| 国产精品电影一区二区三区| 露出奶头的视频| 日韩国内少妇激情av| 夜夜爽天天搞| 精品不卡国产一区二区三区| 桃红色精品国产亚洲av| 亚洲中文字幕日韩| 午夜精品在线福利| 人妻夜夜爽99麻豆av| 能在线免费观看的黄片| 91字幕亚洲| 国内精品久久久久久久电影| 自拍偷自拍亚洲精品老妇| 性欧美人与动物交配| eeuss影院久久|