999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

雙關節剛性機器人自適應BP神經網絡算法

2018-01-12 04:28:30楊航劉凌倪駿康張誠
西安交通大學學報 2018年1期

楊航, 劉凌, 倪駿康, 張誠

(西安交通大學電氣工程學院, 710049, 西安)

雙關節剛性機器人自適應BP神經網絡算法

楊航, 劉凌, 倪駿康, 張誠

(西安交通大學電氣工程學院, 710049, 西安)

為解決反向傳播(BP)神經網絡學習速度慢、泛化能力弱以及易陷入局部極小值等問題,提出了一種雙關節剛性機器人自適應BP神經網絡算法。首先,設計了一種新結構的參數可調激活函數,其映射范圍、斜率因子、水平和豎直位置等參數皆可自適應調整,以使BP神經網絡具有更強的非線性映射能力;然后,設計了一個模糊推理器來整定激活函數的斜率因子,以使斜率因子保持參數最優;最后,設計了一套控制系統并應用于雙關節剛性機器人系統的位置跟蹤控制,采用自適應BP神經網絡算法對該系統位置控制器的比例增益、積分增益和微分增益進行了實時整定。仿真結果表明:與使用經典固定參數S型激活函數的BP神經網絡算法相比,所提出的自適應BP神經網絡算法考慮了激活函數的自適應性,提高了BP神經網絡的學習速度和泛化能力,并抑制了假飽和現象,對位置誤差的收斂速度可提高20倍以上,且可使位置和速度誤差減小并接近于0。

剛性機器人;反向傳播神經網絡;自適應;激活函數;模糊推理

近年來,人們越來越關注機器人的跟蹤控制問題[1],所要實現的目標是減小機械臂實際移動軌跡和設定軌跡之間的偏差,并達到所期望的精度。然而,機器人系統是一個非線性、強耦合、高階和多變量的系統,其數學模型是一個具有時變結構的非線性模型[2]。由于機器人的動作常常可通過不同的方式和路徑來完成,手臂運動解不唯一,故需要在一定約束條件下考慮優化決策與控制問題,但是基于傳統算法的控制器已無法滿足控制要求。隨著機器人技術的發展,與其相關的控制算法逐漸豐富起來,如神經網絡控制、自適應控制、自校正控制、魯棒控制、變結構控制和預測控制等。由于神經網絡需要相對較少的系統動力學信息,且反向傳播(BP)神經網絡已被證明能夠在一定條件下以任何所需的精度逼近各種非線性函數,所以神經網絡算法已廣泛應用于機器人系統的控制[3]。

文獻[3-6]提出了一些使用神經網絡技術解決機器人系統控制問題的方案。雖然大多數神經網絡采用BP神經網絡或其改進形式,但經典的BP神經網絡有缺陷,這使它并不總能滿足實際應用的需求。激活函數是神經網絡解決非線性問題的關鍵,而經典的S型激活函數的映射范圍、斜率因子及位置參數是固定不變的,故其在學習及訓練過程中無法自動改變函數結構和尋找最佳參數,由此會導致BP神經網絡易陷入局部極小值、收斂速度慢和泛化能力弱等問題的出現。

BP神經網絡的學習能力與網絡復雜度之間的矛盾制約了其非線性映射能力,導致學習速度和收斂精度不理想。近十幾年來,針對BP神經網絡算法的不足,研究人員做了大量深入的研究,并提出了許多改進算法,如使用動量項加快離線訓練速度[7]、快速傳播算法[8]、二階優化[9]以及最優濾波法[10]等。文獻[11]提出了參數可調的激活函數模型;文獻[12]采用免疫遺傳算法離線優化和神經網絡自學習在線調整隸屬函數參數的方法,增強了水下機器人對環境變化的反應能力;文獻[13]提出了一種自適應T-S (Takagi-Sugeno)模糊神經網絡控制器,以實現欠驅動無人艇直線航跡的跟蹤控制。這些改進雖然改善了神經網絡算法的性能,但只是簡單調節了激活函數的陡峭程度和映射范圍,而映射區間的長度并不能改變,無法使神經網絡達到最佳的非線性映射能力。

為了克服由于固定的S型激活函數參數所導致的問題,本文提出了一種基于新型參數可調激活函數的自適應BP神經網絡算法,以提高閉環系統的魯棒性;隨后,提出了一種新結構的激活函數,其映射范圍、斜率因子、水平位置和豎直位置等參數皆可自適應調整,以使神經網絡具有更強的非線性映射能力,并提高學習速度和泛化能力;最后,通過設計的模糊推理器來整定激活函數的斜率因子。仿真結果證實了本文算法的有效性。

1 剛性機器人系統的數學模型

對于一個n關節剛性機器人系統,考慮摩擦力的動態方程[14]可以表示為

(1)

式中:M(q)∈n×n為對稱正定的慣性矩陣n×n為向心力和哥氏力矩陣n為摩擦力矢量;G(q)∈n為重力矢量n分別為機械臂的位置、速度和加速度;τ∈n為力矩矢量。

雙關節剛性機器人機械臂示意圖如圖1所示,其中:mi、li是第i個支架的質量和長度;q1、q2為關節1和關節2的位置角(rad);lci是從第i-1個關節點到第i個支架質心的長度,假設連桿材質均勻,則lc1=0.5l1,lc2=0.5l2;Ii為第i個支架的轉動慣量。

圖1 雙關節剛性機器人機械臂示意圖

雙關節剛性機器人的2個機械臂由2個永磁無刷電機驅動,而電機由電流伺服放大器驅動。伺服放大器的增益矩陣Ksa=kI(k>0),其中Ksa∈n×n,I是一個具有相應維數的單位矩陣,k為正的常數。永磁無刷電機的輸出轉矩[14]為

τ(t)=Kmim(t)=Kmid(t)=KmKsau(t)

(2)

式中:Km∈n×n為電機轉矩常數矩陣;im∈n為驅動電機的電流矢量;id∈n為電機的理想電流矢量;u∈n為伺服放大器的輸入電壓矢量。

本文的控制目標是驅動雙關節剛性機器人系統的機械臂以參考速度跟蹤參考位置軌跡,并實現位置跟蹤誤差和速度跟蹤誤差的最小化。

2 自適應BP神經網絡算法

2.1 參數可調的激活函數

激活函數是神經網絡解決非線性問題的關鍵。激活函數的種類繁多,如線性函數、閾值函數和S型函數。這里將基于經典固定參數S型激活函數的BP神經網絡定義為FPBP神經網絡,相應的算法稱為FPBP神經網絡算法。本文提出一種新型的參數可調激活函數,結構為

f=f(x,a,b,c,d)=

(3)

式中:a為區間因子(正實數);b為斜率因子(正實數);c和d為位置參數。

顯然,只要給出相應的參數,傳統的S型激活函數(1+e-x)-1只是式(3)的一個特例。因此,本文提出的激活函數具有更豐富的結構形式和非線性表達能力。

2.2 參數更新算法

相關文獻已經證明,只要隱藏層中節點的數量足夠,3層BP神經網絡具有模擬任何復雜非線性映射的能力[15]。本文提出的自適應BP神經網絡拓撲圖如圖2所示。

圖2 自適應BP神經網絡拓撲圖

在圖2中:自適應BP神經網絡的輸入層有m個輸入量xi(i=1,2,…,m);隱藏層有s個神經元,神經元閾值為γj(j=1,2,…,s),激活函數為本文提出的參數可調激活函數fj(j=1,2,…,s),其輸入量為Gj(j=1,2,…,s),輸出量為gj(j=1,2,…,s);自適應BP神經網絡的輸出層有n個神經元,層神經元閾值為θk(k=1,2,…,n),激活函數為本文的參數可調激活函數fok(k=1,2,…,n), 其輸入量為Yk(k=1,2,…,n),輸出量為yk(k=1,2,…,n);ωij表示輸入層第i個輸入量與隱藏層第j個神經元連接的權值;vjk表示隱藏層第j個輸入量與輸出層第k個神經元連接的權值。于是,有

(4)

gj=fj(Gj)

(5)

(6)

yk=fok(Yk)

(7)

方程(4)~(7)完成了自適應BP神經網絡從輸入到輸出的映射。將fj(j=1,2,…,s)和fok(k=1,2,…,n)簡寫為f,并對參數x、a、b、c和d分別求偏導,可得

(8)

(9)

(10)

(11)

(12)

(13)

則總誤差

(14)

2.2.1 神經網絡權值更新 輸出層權值vjk的更新公式為

(15)

考慮學習率和動量因子,輸出層權值更新為

(16)

(17)

式中:ηv為v的學習率;αv為v的動量因子。

隱藏層權值ωij的更新公式為

(18)

考慮學習率和動量因子,隱藏層權值更新為

(19)

(20)

式中:ηω為ω的學習率;αω為ω的動量因子。

2.2.2 激活函數參數更新

(1)區間因子a的更新。輸出層區間因子ak的更新為

(21)

隱藏層區間因子aj的更新為

(22)

考慮學習率和動量因子,區間因子更新為

(23)

(24)

式中:ηa為a的學習率;αa為a的動量因子。

(2)斜率因子b的更新。模糊推理不依賴于系統的數學模型,能獲得更好的動態控制性能,實現高精度、強魯棒性控制,提高系統的智能化程度,因此采用模糊推理算法設計斜率因子b的更新算法。將模糊子集劃分為{Z,PS,PM,PB},分別代表{零,正小,正中,正大}。選取高斯函數和三角形函數相結合的隸屬度函數。

首先,當位置誤差較大時,選用高斯隸屬度函數,且以快速減小誤差為主要目的。此時應選用較大的斜率因子,但當誤差迅速減小時,為防止神經網絡發生飽和,斜率因子可保持不變。

其次,當位置誤差減小到設定的區段之后,選用三角形隸屬度函數,且以防止BP神經網絡出現假飽和現象為主要目的。此時若誤差緩慢地減小,應選用較大的斜率因子,使神經網絡跳出平坦區域。

最后,當位置誤差繼續減小到模糊子集為Z的區間時,應選用較小的斜率因子,以防止BP神經網絡過飽和,提高系統的穩定性。根據大量的實驗,制定了斜率因子變化量Δb的模糊規則表,見表1。

表1 Δb的模糊規則表

模糊推理選用馬丹尼(Mardani)型,整定好的斜率因子b表示為

b(t+1)=b(0)+Δb(t)

(25)

式中:b(0)>0為斜率因子初值;Δb(t)為斜率因子變化量。為避免依賴精確的數學模型和長時間參數整定的困難,選用重心法作為解模糊的方法,并導出斜率因子變化量的查詢表,以實現斜率因子的在線整定。

(3)位置參數c的更新。輸出層位置參數ck的更新為

(26)

隱藏層位置參數cj的更新為

(27)

考慮學習率和動量因子,位置參數更新為

(28)

(29)

式中:ηc為c的學習率;αc為c的動量因子。

(4)位置參數d的更新。輸出層位置參數dk的更新為

(30)

隱藏層位置參數dj的更新為

(31)

考慮學習率和動量因子,位置參數更新為

(32)

(33)

式中:ηd為d的學習率;αd為d的動量因子。

學習率決定了參數移動到最優值的速度快慢,所以學習率對于算法的性能至關重要。附加動量因子使網絡在修正其權值時,不僅考慮誤差在梯度上的作用,而且考慮在誤差曲面上變化趨勢的影響。在沒有附加動量的作用時,網絡可能陷入淺的局部極小值,而利用附加動量的作用則有可能滑過這些極小值。

對于不同的輸入樣本,最佳的激活函數參數可能是不一致的。然而,由上述分析可以得出如下結論:本文所提出的新激活函數的參數可以自適應調整,使其具有更豐富的非線性表達能力,特別地,斜率因子可以根據位置誤差和它們的變化率,通過模糊推理器進行推理,因此斜率因子也具有自適應性。

3 仿真研究

3.1 系統模型

考慮圖1所示的雙關節剛性機器人系統機械臂。機器人系統的參數為:m1=2.0 kg,m2=1.0 kg,l1=0.2 m,l2=0.15 m,I1=61.25 g·m2,I2=20.42 g·m2,g=9.8 m/s2。為了便于比較FPBP神經網絡算法和本文提出的自適應BP神經網絡算法,除了激活函數的結構和激活函數的參數更新方法外,其余所有的條件和參數設置相同。為了更加客觀地進行對比,2種BP神經網絡算法的學習率和動量因子設置相同的值,即η=0.03,α=0.05。

對于式(1)所表示的2自由度機器人系統動態方程,各參數[16]為

(34)

(35)

(36)

(37)

qd=[qd1,qd2]T=[sin(πt),cos(πt)]T

(38)

其中t∈[0,tf]且tf=20 s。

FPBP神經網絡算法使用經典固定參數的S型激活函數,即

(39)

自適應BP神經網絡算法使用新型激活函數,初始值設定為:a(0)=b(0)=d(0)=0.5 rad,c(0)=0。

(40)

自適應模糊神經網絡控制系統框圖見圖3。

分別為機械臂的給定角位置、給定角速度和給定角加速度;qi(i=1,2)為關節的實際位置角;ei=qdi-qi(i=1,2)為位置角誤差為位置角誤差的變化率。圖3 自適應模糊神經網絡控制系統框圖

設u(k)為第k次采樣時刻控制器的輸出值,可得離散的位置控制器輸出

u(k)=u(k-1)+KpΔe(k)+Kie(k)+

Kd[Δe(k)-Δe(k-1)]

(41)

式中:Kp、Ki、Kd分別為位置控制器的比例增益、積分增益和微分增益;Δe(k)=e(k)-e(k-1)。

3.2 結果及分析

雙關節剛性機器人系統2個關節的位置跟蹤軌跡如圖4和圖5所示,從中可以看出:機器人系統的所有機械臂關節都可以通過自適應BP神經網絡算法成功地跟蹤所設定的軌跡,且位置誤差可以收斂并接近于0,而使用FPBP神經網絡算法不能準確地跟蹤期望軌跡,且穩態位置角誤差較大。

圖4 關節1的位置跟蹤軌跡

圖5 關節2的位置跟蹤軌跡

從圖4和圖5還可以看出,在機器人系統的參數相同的前提下,采用FPBP神經網絡算法的位置誤差比采用自適應BP神經網絡算法的位置誤差大,即FPBP神經網絡算法的跟蹤性能不如自適應BP神經網絡算法好。為了更清晰地比較2種算法的性能,對每個控制周期內采樣的一組位置跟蹤誤差求取均方根誤差(Erms),結果如圖6、圖7所示。

圖6 位置跟蹤誤差的均方根誤差曲線

圖7 速度跟蹤誤差的均方根誤差曲線

由圖6可以看出:采用FPBP神經網絡算法的機器人機械臂關節不能成功地跟蹤期望軌跡,位置跟蹤誤差下降緩慢;相反,采用自適應BP神經網絡算法的機器人機械臂關節的位置跟蹤誤差能在很短的時間內收斂并接近于0。從圖7可知,采用自適應BP神經網絡算法比采用FPBP神經網絡算法的速度跟蹤誤差小,說明自適應BP神經網絡算法極大地提高了BP神經網絡權值的收斂速度。自適應BP神經網絡算法的上述優勢的體現正是由于它能對激活函數的映射范圍、斜率及位置等參數進行調節,特別是能通過模糊推理器對激活函數的斜率因子進行自適應調整,從而保證了機械臂關節的高精度軌跡跟蹤性能。

4 結束語

本文針對BP神經網絡學習速度慢、泛化能力弱和易陷入局部極小值等問題,研究了自適應神經網絡控制策略,提出了一種新結構的參數可調激活函數,其映射范圍、斜率因子、水平位置和豎直位置等參數皆可自適應地調整,特別是激活函數的斜率因子可以通過模糊推理器進行整定,以使神經網絡具有更強的非線性映射能力。仿真結果表明,除激活函數的結構和參數更新方法不同之外,在相同參數條件下,自適應BP神經網絡算法比FPBP神經網絡算法具有更強的非線性表達能力、更快的收斂速度和更高的收斂精度,可以成功地驅動機械臂以參考速度跟蹤參考位置軌跡,且具有更好的動態性能。

[1] WANG L, LIU Z, CHEN C L, et al. Energy-efficient SVM learning control system for biped walking robots [J]. IEEE Transactions on Neural Networks & Learning Systems, 2013, 24(5): 831-837.

[2] LI Z, LI J, KANG Y. Adaptive robust coordinated control of multiple mobile manipulators interacting with rigid environments [J]. Automatica, 2010, 46(12): 2028-2034.

[3] CHEN C L P, WEN G X, LIU Y J, et al. Adaptive consensus control for a class of nonlinear multiagent time-delay systems using neural networks [J]. IEEE Transactions on Neural Networks & Learning Systems, 2014, 25(6): 1217-1226.

[4] CHEN C L P, LIU Y J, WEN G X. Fuzzy neural network-based adaptive control for a class of uncertain nonlinear stochastic systems [J]. IEEE Transactions on Cybernetics, 2014, 44(5): 583-593.

[5] LIU Y J, TONG S. Adaptive NN tracking control of uncertain nonlinear discrete-time systems with nonaffine dead-zone input [J]. IEEE Transactions on Cybernetics, 2014, 45(3): 497-505.

[6] LIU Y J, CHEN C L P, WEN G X, et al. Adaptive neural output feedback tracking control for a class of uncertain discrete-time nonlinear systems [J]. IEEE Transactions on Neural Networks, 2011, 22(7): 1162-1167.

[7] RUMELHART D E, HINTON G E, WILLIAMS R J. Learning internal representations by error propagation [M]∥Parallel Distributed Processing: Exploration of the Microstructure of Cognition. Cambridge, MA, USA: MIT Press, 1986: 318-362.

[8] 林軍, 倪宏, 孫鵬, 等. 一種采用神經網絡PID控制的自適應資源分配方法 [J]. 西安交通大學學報, 2013, 47(4): 112-117.

LIN Jun, NI Hong, SUN Peng, et al. Adaptive resource allocation based on neural network PID control [J]. Journal of Xi’an Jiaotong University, 2013, 47(4): 112-117.

[9] TSAI C, CHIH Y, WONG W H, et al. A hardware-efficient sigmoid function with adjustable precision for a neural network system [J]. IEEE Transactions on Circuits & Systems: II Express Briefs, 2015, 62(11): 1073-1077.

[10] 柏子游, 張勇, 虞烈. 基于模糊數學和神經網絡的數學形態學方法 [J]. 西安交通大學學報, 1999, 33(4): 104-105.

BAI Ziyou, ZHANG Yong, YU Lie. Mathematical morphological operations based on fuzzy set and neural network [J]. Journal of Xi’an Jiaotong University, 1999, 33(4): 104-105.

[11] CARRASCO-ROBLES M, SERRANO L. A novel minimum-size activation function and its derivative [J]. IEEE Transactions on Circuits & Systems: II Express Briefs, 2009, 56(4): 280-284.

[12] 梁霄, 張均東, 李巍, 等. 水下機器人T-S型模糊神經網絡控制 [J]. 電機與控制學報, 2010, 14(7): 99-104.

LIANG Xiao, ZHANG Jundong, LI Wei, et al. T-S fuzzy neural network control for autonomous underwater vehicles [J]. Electric Machines and Control, 2010, 14(7): 99-104.

[13] 董早鵬, 劉濤, 萬磊, 等. 基于Takagi-Sugeno模糊神經網絡的欠驅動無人艇直線航跡跟蹤控制 [J]. 儀器儀表學報, 2015, 36(4): 863-870.

DONG Zaopeng, LIU Tao, WAN Lei, et al. Straight-path tracking control of underactuated USV based on Takagi-Sugeno fuzzy neural network [J]. Chinese Journal of Scientific Instrument, 2015, 36(4): 863-870.

[14] LEWIS F L, YESILDIREK A, LIU K. Multilayer neural-net robot controller with guaranteed tracking performance [J]. IEEE Transactions on Neural Networks, 1996, 6(3): 703-715.

[15] 鄭南寧, 王龍, 胡超, 等. BP神經網絡的改進及其用于手寫數字識別的研究 [J]. 西安交通大學學報, 1992, 26(1): 1-12.

ZHENG Nanning, WANG Long, HU Chao, et al. Improved BP neural net and its application to handwritten numeral recognition [J]. Journal of Xi’an Jiaotong University, 1992, 26(1): 1-12.

[16] HE W, CHEN Y, YIN Z. Adaptive neural network control of an uncertain robot with full-state constraints [J]. IEEE Transactions on Cybernetics, 2016, 46(3): 620-629.

[本刊相關文獻鏈接]

程傳奇,郝向陽,李建勝,等.融合改進A*算法和動態窗口法的全局動態路徑規劃.2017,51(11):137-143.[doi:10.7652/xjtuxb201711019]

尹貴,張小棟,陳江城,等.下肢康復機器人按需輔助自適應控制方法.2017,51(10):39-46.[doi:10.7652/xjtuxb201710 007]

欒玉亮,榮偉彬,吳方勇,等.3-PPSR柔性并聯機器人力控制研究.2017,51(4):85-90.[doi:10.7652/xjtuxb201704013]

黃昔光,黃旭,李啟才.利用共形幾何代數的串聯機器人位置逆解求解方法.2017,51(1):9-12.[doi:10.7652/xjtuxb2017 01002]

周雪峰,孫廣彬,劉曉光,等.應用擴展零力矩點預觀控制和分解動量控制的仿人機器人全身運動規劃方法.2016,50(12):58-63.[doi:10.7652/xjtuxb201612010]

劉貴杰,劉鵬,穆為磊,等.采用能耗最優改進蟻群算法的自治水下機器人路徑優化.2016,50(10):93-98.[doi:10.7652/xjtuxb201610014]

韓冬,劉俊.采用代數建模方法的自適應觀測器設計.2015,49(7):88-91.[doi:10.7652/xjtuxb201507015]

閆志遠,杜志江,吳冬梅.手術機器人運動再現及最優布置系統.2013,47(12):116-122.[doi:10.7652/xjtuxb201312020]

AnAdaptiveBPNeuralNetworkAlgorithmfor2-JointRigidRobots

YANG Hang, LIU Ling, NI Junkang, ZHANG Cheng

(School of Electrical Engineering, Xi’an Jiaotong University, Xi’an 710049, China)

According to the existing problems that back-propagation (BP) neural network with fixed activation function parameters has the shortcomings of slow learning speed, weak generalization ability and ease of falling into local minimum, an adaptive BP neural network algorithm for 2-joint rigid robot is proposed. Firstly, a novel structure of activation function with adjustable parameters is designed. In order to make the BP neural network have better nonlinear mapping capability, its mapping range, steep degree, and horizontal and vertical position parameters can be adjusted adaptively. Then, a fuzzy controller is designed to adjust the slope factor. Thus, the optimality of the slope factor can be guaranteed. Finally, a control system is designed and applied to the tracking control of a 2-joint rigid robotic system. The adaptive BP neural network algorithm is adopted to tune the proportional, integral and differential gains of the position controller. Simulation results show that in comparison with the classical BP neural network algorithm based on Sigmoid activation function with fixed parameters, the adaptive BP neural network algorithm takes the adaptability of activation function into account, improves the learning speed and generalization ability, and restrains false saturation. The convergence rate of position errors can be increased by 20 times, and the position and velocity tracking errors can be reduced to a small value close to zero with the proposed algorithm.

rigid robot; back propagation neural network; adaptability; activation function; fuzzy inference

2017-04-16。

楊航(1992—),男,碩士生;劉凌(通信作者),男,副教授。 基金項目: 國家自然科學基金資助項目(51307130);中央高校基本科研業務費專項資金資助項目(1191320056)。

時間:2017-10-20

網絡出版地址: http:∥kns.cnki.net/kcms/detail/61.1069.T.20171020.1622.008.html

10.7652/xjtuxb201801019

TP183

A

0253-987X(2018)01-0129-07

(編輯 葛趙青)

主站蜘蛛池模板: 日本www色视频| 日韩欧美视频第一区在线观看| 国产精品无码制服丝袜| 中文字幕无码电影| 国产网站免费观看| 国产日韩欧美黄色片免费观看| 国产成人8x视频一区二区| 国产香蕉在线| 日本高清有码人妻| 伊人91在线| 国产丝袜第一页| 在线欧美一区| 国产精品不卡片视频免费观看| 国产成人免费高清AⅤ| 色综合热无码热国产| 精品丝袜美腿国产一区| AV老司机AV天堂| 国产Av无码精品色午夜| 国产jizz| 亚洲V日韩V无码一区二区 | 茄子视频毛片免费观看| 免费播放毛片| 亚洲第一成网站| 国产成人久久777777| 精品亚洲国产成人AV| 91最新精品视频发布页| 永久毛片在线播| 亚洲色成人www在线观看| 91精品国产福利| 国产精品无码久久久久AV| 日韩小视频在线观看| 国产嫩草在线观看| 欧美另类一区| 99久久这里只精品麻豆| 无码精品国产VA在线观看DVD| 国产免费a级片| 国产精品林美惠子在线播放| 亚洲国产欧美目韩成人综合| 青草娱乐极品免费视频| 伊人久久青草青青综合| 欧美色伊人| 国产激情无码一区二区免费| 日本高清有码人妻| 国产噜噜噜| 欧美啪啪一区| 无码人中文字幕| 国产91无码福利在线| 国产无码精品在线| 好久久免费视频高清| 无码综合天天久久综合网| 香蕉久久永久视频| 亚洲av无码久久无遮挡| 国产青榴视频| 国产成人精彩在线视频50| 久久综合五月婷婷| 亚洲人成网站在线播放2019| 国产精品一区二区在线播放| 99在线国产| 麻豆精品在线| 欧美精品亚洲精品日韩专区| 日韩精品资源| 国产一级毛片在线| 特级做a爰片毛片免费69| 亚洲欧美激情小说另类| A级全黄试看30分钟小视频| 日本一本正道综合久久dvd | 香蕉蕉亚亚洲aav综合| 久久久精品国产亚洲AV日韩| 在线看片中文字幕| 国产99热| 91精品情国产情侣高潮对白蜜| 91免费国产高清观看| 免费国产好深啊好涨好硬视频| 日韩精品亚洲人旧成在线| 福利一区在线| 午夜老司机永久免费看片| 成年人国产网站| 国产精品自在拍首页视频8| 国产一区二区三区日韩精品| 色爽网免费视频| 国产精品视频导航| 亚洲男人的天堂网|