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

    群體智能的系統(tǒng)辨識

    2014-03-17 03:55:45陳文雯劉友寬孫建平
    云南電力技術 2014年3期
    關鍵詞:傳遞函數(shù)慣性乘法

    陳文雯,劉友寬,孫建平

    (1.華北電力大學自動化系,河北 保定 071000;2.華北電力大學云南電網(wǎng)公司研究生工作站,昆明 650217; 3.云南電網(wǎng)公司電力研究院,昆明 650217)

    群體智能的系統(tǒng)辨識

    陳文雯1,2,劉友寬3,孫建平1

    (1.華北電力大學自動化系,河北 保定 071000;2.華北電力大學云南電網(wǎng)公司研究生工作站,昆明 650217; 3.云南電網(wǎng)公司電力研究院,昆明 650217)

    概括了系統(tǒng)辨識的方法,重點介紹了最小二乘法、群體智能算法中的粒子群算法和改進的粒子群算法,給出了估計模型的選擇方法,并結合某1000MW火電機組實例,運用兩種方法進行了系統(tǒng)辨識和仿真。仿真結果表明,最小二乘法可以完成對系統(tǒng)的辨識,但存在較大偏差;采用粒子群算法辨識結果良好。

    系統(tǒng)辨識;最小二乘法;粒子群算法;仿真

    1 前言

    系統(tǒng)識別是建模的一種方法。建立數(shù)學模型有兩種方法:解析法和系統(tǒng)辨識。

    系統(tǒng)辨識的過程實質上就是函數(shù)擬合的過程,這里包括傳遞函數(shù)的結構和參數(shù)。因此,所要面臨的是結構優(yōu)化和參數(shù)優(yōu)化的問題。如果已經對系統(tǒng)有了一定的了解,那么可以先給出系統(tǒng)模型描述函數(shù)的結構,然后辨識出函數(shù)中的參數(shù)即可,即把結構 (函數(shù))優(yōu)化問題轉化成參數(shù)優(yōu)化問題。

    文中從傳統(tǒng)系統(tǒng)辨識方法和群體智能系統(tǒng)辨識方法兩方面入手,闡述系統(tǒng)辨識的方法,辨識實際工業(yè)控制系統(tǒng),并對其進行仿真。

    2 最小二乘法

    最小二乘法 (LS)是估計理論的奠基石,其思想是使各次實際觀測值和計算值之間差值的平方乘以度量其精確度的數(shù)值以后的和為最小。最小二乘法容易理解和掌握,利用其所擬定的識別算法在實施上比較簡單。但由于最小二乘估計是非一致的、有偏差的,因而為了克服它的不足,形成了一些以最小二乘法為基礎的辨識方法,如廣義最小二乘法、輔助變量法和增廣矩陣法,以及將一般的最小二乘法與其它方法相結合的方法,譬如最小二乘兩步法和隨機逼近算法[1]。

    3 粒子群算法

    3.1 PSO基本思想

    PSO首先初始化一群隨機粒子 (初始速度、位移及其決定的適應值都隨機化),然后通過迭代搜索最優(yōu)解。在每一次迭代中,粒子通過跟蹤兩個最優(yōu)值來跟新自己,第一個就是粒子本身目前所找到的最優(yōu)解Xbesti,即個體最優(yōu)值。每個粒子都具有記憶功能,Xbesti是它們記住的各自曾經達到的最好位置。另一個最優(yōu)值是整個種群目前找到的最優(yōu)解Xbesti,即全局最優(yōu)值 (假設群體之間存在著某種通信方式,每個粒子都能夠記住目前為止整個群體的最好位置)。

    下面以求某一函數(shù)Q(Xi)的極小值為例,介紹基本粒子群算法的實現(xiàn)方法[2]。

    假設在一個N維的目標搜索空間中 (N相當于Q(Xi)中未知因子個數(shù),也就是優(yōu)化參數(shù)個數(shù)),有m個粒子組成的一個群體 (即m組可能解),其中第 i個粒子的位置表示為向量 Xi= (xi1,xi2,……,xiN);i=1,2,……,m其速度也是一個N維的向量,記為Vi=(vi1,vi2,……, viN)。隨機產生一組 Xi,作為第一代初始種群,將Xi帶入目標函數(shù)Q(Xi)就可以計算出其適應值,根據(jù)適應值的大小衡量Xi的優(yōu)劣。對于最小化問題,目標函數(shù)值越小,對應的適應值越好。設粒子 i迄今為止經歷的最有位置為 Xbesti= (xi1,xi2,……,xiN),相應的適應值記為Qbesti,則粒子的當前最好位置可表示為:

    尋優(yōu)過程中粒子群經歷的最優(yōu)位置記為 Xbestg= (xg1,xg2,…,xgN)其對應適應值即全局最優(yōu)解為Qbestg。則粒子根據(jù)式(2)來更新自己的速度,即:

    式 (2)中,i=1,2,…,m,n=1,2,…,N,t表示第t代。在速度更新時,不應該超出給定的速度范圍,即要求Vi∈[-Vmax,Vmax],單步前進的最大值Vmax根據(jù)粒子的取值區(qū)間長度來確定。

    然后按式 (3)來更新位置向量,即:

    xin(t+1)=xin(t)+vin(t+1)(3)

    式中的變量意義同前。根據(jù)實際問題來確定粒子的取值范圍xin∈[xinmin,xinmax]。

    這樣一代一代地執(zhí)行下去直到達到要求,取得極值。

    對公式中一些符號意義作用的幾點說明:①c1表示認知因子,c2表示社會因子。它們分別代表了向自身極值和全局極值推進的加速極值。實驗結果表明,一般取c1=c2=2比較好,但實際上加速權值是可以變化的,而且如何變化將直接影響尋優(yōu)過程;②r1,r2是0~1之間的隨機變量。

    PSO算法流程:

    1)初始化,包括定義初始種群 (速度—位移模型以及種群大小等),進化代數(shù),還有一些修正改進算法中可能用到的常量。

    2)評價種群。計算初始種群各個粒子的適應度。

    3)求出當前的Qbesti和Qbestg。

    4)進行速度和位置更新。

    5)評價種群。計算新種群中粒子的適應度。

    6)比較Xbesti和Qbestg,若優(yōu)越則替換。

    7)判斷算法結束條件 (包括精度要求和進化代數(shù)要求),滿足則跳出循環(huán),不滿足則跳到 (4)繼續(xù)執(zhí)行。

    以上為基本粒子群算法[3-5]。

    3.2 粒子群算法的改進

    為了改善基本PSO算法的收斂性能,本文提出了帶權重的粒子群算法,即標準粒子群算法。這是對基本粒子群算法最早的一種改進。

    標準 PSO主要是在式 (2)中引入了慣性權重ω,即:

    慣性權重是為了平衡全局搜索和局部搜索而引入的,慣性權重代表了原來速度在下一次迭代中所占的比例,ω較大時,前一速度的影響較大,全局搜索能力比較強;ω較小時,前一速度的影響較小,局部搜索能力比較強。合適的ω值在搜索速度和搜索精度方面起著協(xié)調作用。因此,一般采用慣性權重遞減策略,即在算法的初期取較大的慣性權值ω,使對整個問題空間進行有效的搜索,算法進行后期取較小慣性權值ω,以有利于算法的收斂。慣性權重遞減公式為:

    式中,ωmax和ωmin分別為ω的最大最小值,ω的取值范圍在 [0,1.4]比較合適,但通常取在 [0.8, 1.2];Tmax、t分別是最大的迭代數(shù)和當前的迭代數(shù)。

    另外,Clerc提出的收縮因子法也是一種標準的PSO算法。他是把基本的速度公式,即式 (2)改變?yōu)椋?/p>

    φ=c1+c2,φ>4。通常情況下取 c1=c2=2.05,φ= 4.1,此時γ=0.7298。實驗結果表明,兩種方法差不多,收縮因子更有效率,但是在有些情況下無法得到全局極值點[6-9]。

    4 算法的應用

    4.1 估計模型的選擇

    確定模型的結構是決定模型質量關鍵性的一步[10-12]。

    實際上,一個實際的物理表象,可以用無窮多的數(shù)學模型來描述,物理表象與數(shù)學模型不存在一一對應的關系,我們所能做的就是從各種數(shù)學模型中選擇出一種來近似描述實際的物理表象。這一特性給選擇模型結構帶來困難。我們只好用對各過程領域的先驗知識來假想一個模型結構。專家總結出的經驗模型是:

    式中:K為系統(tǒng)增益;子為純遲延時間常數(shù);T為系統(tǒng)慣性時間常數(shù);β為微分時間常數(shù);當系統(tǒng)為無自平衡時,有自平衡時;為慣性部分的階次。

    由式 (7)系統(tǒng)可以簡化成以下幾種模型結構。

    1)高階對象

    絕大多數(shù)的熱工對象有自平衡能力,并且屬于多階慣性環(huán)節(jié)。一般可以認為它是等容多階對象,定義它為I型對象,即:

    當求出的階次n不是整數(shù)時,用近似的整數(shù)代替。

    2)多容慣性對象

    如果想描述有自平衡對象的細節(jié),則可以用多容慣性對象,定義它為II型對象,即:

    3)具有純遲延的高階慣性對象

    當系統(tǒng)存在純遲延時,可以加入純遲延環(huán)節(jié),定義它為III型對象,即:

    當使用高階對象時,可能會遇到困難,這時可以對其進行降階處理。此外,純遲延對象并不適合于系統(tǒng)分析,這時可以與高階對象互換。如果不要求有特別高的精度,則可以用下面的方法進行升降系統(tǒng)階次與純遲延之間的關系相互轉換來處理。

    如果原傳遞函數(shù)如式 (10)所示,則可把它簡化成:

    兩式中的參數(shù)關系為:

    4)無自平衡能力對象

    對于汽包水位系統(tǒng)等少數(shù)無自平衡能力對象,其傳遞函數(shù)可用式 (13)來描述,并定義它為IV型對象,即:

    5)零穩(wěn)態(tài)對象

    對于具有微分作用的對象,當系統(tǒng)趨于穩(wěn)態(tài)時,輸出趨近于零,把這種對象定義為V型對象,即:

    6)逆向響應系統(tǒng)

    在工程中,存在一種逆向相應系統(tǒng),它的表征是,在階躍擾動作用下,系統(tǒng)的輸出先朝著與最終趨向相反的方向變化,然后再朝著最終趨向變化。汽包鍋爐的蒸汽量階躍擾動引起的汽包水位變化就是逆向響應過程,在熱工里被稱為 “虛假水位”;循環(huán)流化床鍋爐一次風階躍擾動引起的床溫變化也是一個典型的逆向響應過程。逆向響應系統(tǒng)的傳遞函數(shù)如式 (15)和式(16)所示,定義它為VI和VII型對象,即:

    4.2 系統(tǒng)辨識與仿真

    基于上述理論,用粒子群算法和最小二乘法辨識以下兩個系統(tǒng)。

    I.已知某火電機組在負荷下得到的蒸汽量變化對應汽包水位變化的傳遞函數(shù)為:(系統(tǒng)I)(17)

    II.某1000MW超超臨界火電機組,在負荷時,負荷發(fā)生變動后的各變量響應曲線如圖1所示,從該圖中截取的一段數(shù)據(jù)可用于系統(tǒng)辨識。辨識送風量和引風量變化引起爐膛壓力變化的數(shù)學模型。(系統(tǒng)II)

    圖1 1000 MW機組負荷圖

    4.2.1 系統(tǒng)I辨識

    系統(tǒng)I采用最小二乘法辨識的結果為:

    V1=0.150 4,V2=4.945 0,V2=3.615 5,V4=3.685 4,V5=5.504 5。辨識后的傳遞函數(shù)為式(18):

    系統(tǒng)I采用粒子群法辨識的結果為:V1= 0.184 5,V2=8.971 6,V3=3.492 8,V4=2.830 3,V5=5.721 1。辨識后的傳遞函數(shù)為式 (19):

    最小二乘法辨識結果如圖2所示。

    圖2 系統(tǒng)I最小二乘法辨識結果

    把仿真步距設為500,在同一方波信號下觀察比較最小二乘法和粒子群把對系統(tǒng)I的辨識結果。最小二乘辨識見圖2,粒子群算法辨識見圖3。

    圖3 同一方波信號下系統(tǒng)I粒子群辨識結果

    比較圖2、圖3可知,采用最小二乘法辨識的辨識曲線大概能與實測數(shù)據(jù)擬合,但有明顯誤差。采用粒子群發(fā)進行系統(tǒng)辨識時,只要選擇合適的參數(shù),辨識曲線與實測數(shù)據(jù)的擬合效果很好,幾乎完全重合。

    4.2.2 系統(tǒng)II辨識

    讀取的現(xiàn)場數(shù)據(jù)見圖4。系統(tǒng)辨識模型結構為:

    圖4 讀取的現(xiàn)場數(shù)據(jù)

    系統(tǒng)II采用粒子群算法辨識結果見圖5。

    圖5 系統(tǒng)II粒子群算法辨識結果

    辨識結果為:V1=0.577 3,V2=2.126 9,V3=1.859 5。近似為:

    最小二乘法辨識。V1=0.524 3,V2=1.968 7,V3=2.2951,辨識結果見圖6。

    圖6 系統(tǒng)II最小二乘法辨識結果

    5 結束語

    本文介紹了系統(tǒng)辨識的有關概念和方法,系統(tǒng)辨識實質上是結構優(yōu)化和參數(shù)優(yōu)化。重點介紹了粒子群算法和最小二乘法,并根據(jù)火電機組實例,應用這兩種方法進行了系統(tǒng)辨識與仿真。

    粒子群算法是一種全局隨機優(yōu)化算法,具有參數(shù)少、速度快的特點,從仿真結果可以看出粒子群優(yōu)化算法辨識精度較好、速度快,且對輸入輸出觀測量的要求可以根據(jù)實際情況方便選擇。最小二乘法是傳統(tǒng)的辨識方法,收斂速度較快,但存在局部最小問題,易出現(xiàn)產生非奇異矩陣導致收斂失敗的情況。

    [1] 王琳,馬平.系統(tǒng)辨識方法綜述 [J].電力情報,2001, (4):63-66.

    [2] 韓璞,等.智能控制理論及應用 [M].北京:中國電力出版社,2013.

    [3] T.H.Kim,I.Maruta,T.Sugie.Robust PID controller tuning based on the constrained particle swarm optimization[J]. Automatica,2008,44(4):1104-1110.

    [4] 王俊偉,汪定偉.粒子群算法中慣性權重的調整策略[J].計算機工程,2007,33(11):193-195.

    [5] J.Chen,T.C.Huang.Applying neural networks to on-line updated PID controllers for nonlinear process control[J].Journal of Process Control,2004,14(2):211-230.

    [6] 何小賢,朱云龍,王玫.群體智能中的知識涌現(xiàn)與復雜適應性問題綜述研究 [J].信息與控制,2005,34(5): 560-566.

    [7] 李凱.基于微粒群優(yōu)化算法的結構系統(tǒng)識別 [D].同濟大學土木工程學院,2008.

    [8] J.Q.Han.From PID to active disturbance rejection control [J].IEEE Transactions on Industrial Electronics,2009,56(3):900-906.

    [9] V.A.Oliveira,L.V.Cossi,M.C.M.Teixeira,A.M.F.Silva. Synthesis of PID controllers for a class of time delay system [J].Automatica,2009,45(7):1778-1782.

    [10] 劉金琨,沈曉蓉,趙龍.系統(tǒng)辨識理論及MATLAB仿真[M].北京:電子工業(yè)出版社,2013.

    [11] 田謙益.基于群體智能算法的非線性系統(tǒng)參數(shù)辨識[J].佳木斯大學學報 (自然科學版),2010,28(4):502-505.

    [12] 韓璞,朱希彥.自動控制系統(tǒng)數(shù)字仿真 [M].北京:中國電力出版社,1996.

    System Identification Based on Swarm Intelligence

    CHEN Wenwen1,2,LIU Youkuan3,SUN Jianping1
    (1.Department of Automation,North China Electric Power University,Baoding,Hebei 071000; 2.North China Electric Power University Graduate Student Workstations of Yunnan Power Grid Corporation,Kunming 650217; 3.Smart Grid Department of Yunnan Electric Power Research Institute,Kunming 650217)

    Establishment of the system model is necessary when it comes to study the control system.Therefore,system identification plays a crucial role in the study of control system,of which the essences are structural optimization and parameter optimization. The article summarizes the methods of system identification,focusing on the least squares method,the particle swarm optimization (PSO)which included in swarm intelligence algorithm and the improved particle swarm optimization,gives the selection method of estimation model and combine with a 1000MW thermal power instance,use two methods to identify and simulate the control system. Simulation results show that the least squares method of identification can be done on the system,but there is a big deviation and it can get good recognition results by using particle swarm optimization.

    System Identification;Least Squares;PSO;Simulation

    TP15;TP18

    B

    1006-7345(2014)03-0010-05

    2014-05-06

    陳文雯 (1989),女,碩士研究生,華北電力大學云南電網(wǎng)公司研究生工作站,從事控制理論與控制工程、新能源等領域研究 (email)961689821@qq.com。

    劉友寬 (1973),男,云南電網(wǎng)公司電力研究院,長期從事水火電機組控制研究、新能源遠程控制系統(tǒng)研究。

    孫建平 (1961),男,教授,華北電力大學自動化系,長期從事控制理論與控制工程、模式識別領域的教學與研究。

    猜你喜歡
    傳遞函數(shù)慣性乘法
    你真的了解慣性嗎
    算乘法
    沖破『慣性』 看慣性
    我們一起來學習“乘法的初步認識”
    《整式的乘法與因式分解》鞏固練習
    把加法變成乘法
    基于LabVIEW的火焰?zhèn)鬟f函數(shù)測量系統(tǒng)
    測控技術(2018年7期)2018-12-09 08:58:46
    無處不在的慣性
    紅外成像系統(tǒng)的調制傳遞函數(shù)測試
    普遍存在的慣性
    少妇熟女aⅴ在线视频| 久久久久精品国产欧美久久久| 熟妇人妻久久中文字幕3abv| 日韩欧美精品v在线| 精品久久久久久久毛片微露脸| 久久午夜亚洲精品久久| 日韩大尺度精品在线看网址| 国产精品久久久久久久电影 | 亚洲性夜色夜夜综合| 麻豆成人午夜福利视频| 成人欧美大片| 亚洲av成人不卡在线观看播放网| 欧美日韩中文字幕国产精品一区二区三区| bbb黄色大片| 久久久精品大字幕| 1000部很黄的大片| 成人特级av手机在线观看| av女优亚洲男人天堂 | 综合色av麻豆| 成人鲁丝片一二三区免费| 亚洲av成人一区二区三| 免费看a级黄色片| 午夜福利免费观看在线| 国产精品亚洲av一区麻豆| 欧美乱码精品一区二区三区| 2021天堂中文幕一二区在线观| 90打野战视频偷拍视频| 成人鲁丝片一二三区免费| 欧美日本视频| 97人妻精品一区二区三区麻豆| 我要搜黄色片| 校园春色视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆| 香蕉丝袜av| 蜜桃久久精品国产亚洲av| 亚洲欧美日韩东京热| 国产成人系列免费观看| 香蕉丝袜av| 最新在线观看一区二区三区| 成人国产综合亚洲| 精品国产美女av久久久久小说| 亚洲五月天丁香| 亚洲精品乱码久久久v下载方式 | 国产精品久久久久久精品电影| 国产成人精品久久二区二区91| 亚洲欧美精品综合久久99| 国产精品乱码一区二三区的特点| 国产不卡一卡二| 久久久久久久精品吃奶| 欧洲精品卡2卡3卡4卡5卡区| av天堂中文字幕网| 亚洲精品美女久久av网站| 最近在线观看免费完整版| 中文字幕久久专区| 日韩人妻高清精品专区| 国产精品日韩av在线免费观看| 国产精品美女特级片免费视频播放器 | 午夜福利18| 少妇裸体淫交视频免费看高清| 色视频www国产| 无限看片的www在线观看| 少妇人妻一区二区三区视频| 夜夜爽天天搞| 丰满人妻熟妇乱又伦精品不卡| 日韩精品中文字幕看吧| 在线永久观看黄色视频| 免费观看的影片在线观看| 欧美大码av| 熟女电影av网| 深夜精品福利| 日本撒尿小便嘘嘘汇集6| 99久久99久久久精品蜜桃| 亚洲真实伦在线观看| 99国产综合亚洲精品| 成人鲁丝片一二三区免费| 一个人看视频在线观看www免费 | 精品久久久久久久末码| 免费观看精品视频网站| 女人高潮潮喷娇喘18禁视频| 我要搜黄色片| 日韩欧美国产在线观看| 精品电影一区二区在线| 狠狠狠狠99中文字幕| 国产精品久久电影中文字幕| 身体一侧抽搐| 九九热线精品视视频播放| 日本撒尿小便嘘嘘汇集6| 老司机深夜福利视频在线观看| 99re在线观看精品视频| 两性夫妻黄色片| 成人一区二区视频在线观看| 成人三级黄色视频| 亚洲国产欧美一区二区综合| 国产激情欧美一区二区| 在线视频色国产色| 日本一二三区视频观看| 免费一级毛片在线播放高清视频| 丁香欧美五月| 嫩草影院入口| 亚洲 欧美一区二区三区| 此物有八面人人有两片| 母亲3免费完整高清在线观看| av黄色大香蕉| 此物有八面人人有两片| 母亲3免费完整高清在线观看| 中文字幕最新亚洲高清| 国产乱人视频| 夜夜爽天天搞| 久久久久久久久免费视频了| 男女午夜视频在线观看| 免费看十八禁软件| 窝窝影院91人妻| 97碰自拍视频| 久久天躁狠狠躁夜夜2o2o| 国产综合懂色| 99国产极品粉嫩在线观看| 国产伦人伦偷精品视频| 久久热在线av| 精品久久久久久,| 啪啪无遮挡十八禁网站| 国产激情偷乱视频一区二区| 一区福利在线观看| 一进一出抽搐动态| 网址你懂的国产日韩在线| 在线免费观看的www视频| 精品熟女少妇八av免费久了| 又黄又粗又硬又大视频| 97超视频在线观看视频| 欧美中文日本在线观看视频| 国产一区二区三区视频了| 亚洲无线在线观看| 特大巨黑吊av在线直播| 欧美一区二区精品小视频在线| 久久中文字幕人妻熟女| 一本精品99久久精品77| 午夜影院日韩av| 老鸭窝网址在线观看| 亚洲天堂国产精品一区在线| 欧美一区二区精品小视频在线| 欧美性猛交╳xxx乱大交人| 精品日产1卡2卡| 亚洲成人免费电影在线观看| 在线观看美女被高潮喷水网站 | 性色avwww在线观看| 国产精品电影一区二区三区| 99re在线观看精品视频| 国产成人aa在线观看| 欧美绝顶高潮抽搐喷水| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 999久久久精品免费观看国产| 欧美xxxx黑人xx丫x性爽| 精品一区二区三区av网在线观看| 精华霜和精华液先用哪个| 欧美xxxx黑人xx丫x性爽| 99热精品在线国产| 久久亚洲精品不卡| 在线观看日韩欧美| 一二三四在线观看免费中文在| 亚洲精品美女久久av网站| 美女高潮的动态| 国产精品1区2区在线观看.| 国产成年人精品一区二区| 午夜激情福利司机影院| 婷婷精品国产亚洲av| АⅤ资源中文在线天堂| 午夜成年电影在线免费观看| 午夜免费成人在线视频| 欧美成狂野欧美在线观看| 亚洲国产日韩欧美精品在线观看 | 欧美xxxx黑人xx丫x性爽| 成人一区二区视频在线观看| 少妇人妻一区二区三区视频| 给我免费播放毛片高清在线观看| 国内精品一区二区在线观看| 国产一区二区在线av高清观看| 成人一区二区视频在线观看| 午夜a级毛片| 精品欧美国产一区二区三| 国产午夜精品论理片| 久久99热这里只有精品18| 国产淫片久久久久久久久 | 久久草成人影院| 女同久久另类99精品国产91| 免费观看人在逋| 国产 一区 欧美 日韩| 欧美又色又爽又黄视频| 亚洲成av人片免费观看| 白带黄色成豆腐渣| 听说在线观看完整版免费高清| 亚洲精品乱码久久久v下载方式 | 黄频高清免费视频| 麻豆一二三区av精品| 在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久国产a免费观看| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区激情视频| 日日干狠狠操夜夜爽| 怎么达到女性高潮| 免费av毛片视频| 亚洲人与动物交配视频| 久久香蕉精品热| 欧美性猛交╳xxx乱大交人| 99久久成人亚洲精品观看| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 国产成+人综合+亚洲专区| 精品国产美女av久久久久小说| 亚洲人成电影免费在线| 国产亚洲精品久久久久久毛片| 少妇丰满av| 嫩草影院入口| 国产午夜精品论理片| 亚洲精品在线观看二区| 亚洲精华国产精华精| 国产精品野战在线观看| 在线免费观看的www视频| 欧美3d第一页| 国产三级中文精品| 男人和女人高潮做爰伦理| 特级一级黄色大片| 欧美性猛交黑人性爽| 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 俺也久久电影网| 好看av亚洲va欧美ⅴa在| 欧美黄色片欧美黄色片| 精品99又大又爽又粗少妇毛片 | 亚洲国产中文字幕在线视频| 亚洲美女黄片视频| 99在线视频只有这里精品首页| 欧美激情久久久久久爽电影| 国产精品女同一区二区软件 | 嫩草影院精品99| 伦理电影免费视频| 国产精品久久久av美女十八| 国产美女午夜福利| 亚洲一区二区三区色噜噜| 国产成人精品久久二区二区91| 99久久精品一区二区三区| 日本 av在线| bbb黄色大片| 国产亚洲精品久久久久久毛片| 国产伦人伦偷精品视频| 国产亚洲精品av在线| 99精品欧美一区二区三区四区| 中出人妻视频一区二区| 日本与韩国留学比较| 后天国语完整版免费观看| 亚洲欧美日韩东京热| 搡老熟女国产l中国老女人| 国产精品久久久久久精品电影| 成年免费大片在线观看| 国产三级中文精品| 欧美日韩福利视频一区二区| 999久久久精品免费观看国产| 69av精品久久久久久| 亚洲av片天天在线观看| 亚洲中文av在线| 天堂影院成人在线观看| 久久这里只有精品中国| 婷婷丁香在线五月| 亚洲国产精品sss在线观看| 熟女人妻精品中文字幕| 免费在线观看影片大全网站| 十八禁人妻一区二区| 成人高潮视频无遮挡免费网站| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕久久专区| 国产精品一及| 欧美极品一区二区三区四区| 国产精品久久久av美女十八| 美女高潮的动态| 欧美一区二区精品小视频在线| 一边摸一边抽搐一进一小说| 国产精品一区二区精品视频观看| 久久久久久久精品吃奶| 在线十欧美十亚洲十日本专区| 九九热线精品视视频播放| 黑人欧美特级aaaaaa片| 免费电影在线观看免费观看| 成人亚洲精品av一区二区| 99国产极品粉嫩在线观看| 亚洲欧美日韩高清专用| 亚洲熟女毛片儿| 免费搜索国产男女视频| 色老头精品视频在线观看| 亚洲性夜色夜夜综合| 国产精品电影一区二区三区| 国产欧美日韩精品亚洲av| 精品久久久久久久久久久久久| 五月伊人婷婷丁香| 久久草成人影院| 欧美又色又爽又黄视频| 午夜激情欧美在线| 日韩免费av在线播放| 青草久久国产| 999久久久精品免费观看国产| 最近最新免费中文字幕在线| 国产精品98久久久久久宅男小说| 夜夜爽天天搞| 国产乱人伦免费视频| 一夜夜www| 亚洲国产日韩欧美精品在线观看 | 午夜激情欧美在线| 日韩人妻高清精品专区| 国产美女午夜福利| 大型黄色视频在线免费观看| tocl精华| 每晚都被弄得嗷嗷叫到高潮| 男插女下体视频免费在线播放| 网址你懂的国产日韩在线| 国产黄色小视频在线观看| 亚洲av片天天在线观看| 亚洲狠狠婷婷综合久久图片| 黑人操中国人逼视频| 成年女人永久免费观看视频| 可以在线观看毛片的网站| 亚洲18禁久久av| 国内久久婷婷六月综合欲色啪| 国产一区二区在线观看日韩 | 观看免费一级毛片| 九九久久精品国产亚洲av麻豆 | 人人妻人人澡欧美一区二区| 夜夜夜夜夜久久久久| 日韩精品青青久久久久久| 成人性生交大片免费视频hd| 两个人的视频大全免费| 麻豆国产av国片精品| 成年女人看的毛片在线观看| 一本一本综合久久| 12—13女人毛片做爰片一| 成年人黄色毛片网站| 中文亚洲av片在线观看爽| 男女视频在线观看网站免费| 黄色成人免费大全| 桃红色精品国产亚洲av| 亚洲专区国产一区二区| 99riav亚洲国产免费| 免费av不卡在线播放| 99久久无色码亚洲精品果冻| av福利片在线观看| 麻豆av在线久日| 国产精品1区2区在线观看.| 国产精品久久久久久人妻精品电影| 不卡av一区二区三区| 欧美日韩乱码在线| 欧美日韩瑟瑟在线播放| 久久亚洲真实| 欧美日韩综合久久久久久 | 国产精品国产高清国产av| 亚洲精品中文字幕一二三四区| 又黄又粗又硬又大视频| 岛国在线免费视频观看| 网址你懂的国产日韩在线| 色播亚洲综合网| 手机成人av网站| 91久久精品国产一区二区成人 | av在线蜜桃| 久久午夜亚洲精品久久| 久久久久久人人人人人| 在线免费观看不下载黄p国产 | 精品99又大又爽又粗少妇毛片 | 国产激情欧美一区二区| 九九热线精品视视频播放| 一区二区三区国产精品乱码| 俄罗斯特黄特色一大片| 日本熟妇午夜| 亚洲国产欧美一区二区综合| 亚洲熟女毛片儿| 国产精华一区二区三区| 12—13女人毛片做爰片一| 久9热在线精品视频| 色哟哟哟哟哟哟| 久久久国产成人免费| 国产高潮美女av| 欧美日韩福利视频一区二区| 午夜免费成人在线视频| 精品久久蜜臀av无| 草草在线视频免费看| 蜜桃久久精品国产亚洲av| 国产高清激情床上av| 丰满人妻熟妇乱又伦精品不卡| 精品国产美女av久久久久小说| 国产高清videossex| 宅男免费午夜| 欧美最黄视频在线播放免费| a级毛片a级免费在线| 欧美日本视频| 最近最新中文字幕大全免费视频| 亚洲欧美精品综合一区二区三区| 禁无遮挡网站| 男人舔女人的私密视频| 亚洲人成伊人成综合网2020| 国产精品久久久av美女十八| 久久伊人香网站| 欧美午夜高清在线| 免费观看精品视频网站| 中文字幕熟女人妻在线| 国产精品乱码一区二三区的特点| 日本 欧美在线| 在线观看免费视频日本深夜| 亚洲成av人片免费观看| 国产乱人视频| 国产单亲对白刺激| 久久婷婷人人爽人人干人人爱| 亚洲欧洲精品一区二区精品久久久| 久久中文看片网| 日韩大尺度精品在线看网址| 国产精品爽爽va在线观看网站| 久久久久免费精品人妻一区二区| 国产高清有码在线观看视频| 精品不卡国产一区二区三区| 亚洲国产色片| 村上凉子中文字幕在线| 国产成人啪精品午夜网站| 香蕉av资源在线| 夜夜躁狠狠躁天天躁| 床上黄色一级片| 亚洲成av人片在线播放无| 观看美女的网站| 曰老女人黄片| 给我免费播放毛片高清在线观看| 伊人久久大香线蕉亚洲五| 夜夜夜夜夜久久久久| 国产真人三级小视频在线观看| 日本在线视频免费播放| 岛国视频午夜一区免费看| 99精品在免费线老司机午夜| 亚洲九九香蕉| 亚洲欧美精品综合一区二区三区| 好男人在线观看高清免费视频| 亚洲 国产 在线| 欧美丝袜亚洲另类 | 最近最新中文字幕大全电影3| 长腿黑丝高跟| 亚洲一区二区三区不卡视频| 男女视频在线观看网站免费| tocl精华| 欧美+亚洲+日韩+国产| 欧美丝袜亚洲另类 | 国产激情欧美一区二区| 亚洲狠狠婷婷综合久久图片| 国产毛片a区久久久久| 久久天躁狠狠躁夜夜2o2o| 国产一区二区三区视频了| 国产成人av教育| 亚洲国产精品999在线| 午夜福利成人在线免费观看| 1024香蕉在线观看| 99国产精品99久久久久| 两个人视频免费观看高清| 久久国产精品影院| 可以在线观看的亚洲视频| 亚洲美女黄片视频| 精品久久久久久久久久久久久| 亚洲av成人av| 久久精品亚洲精品国产色婷小说| 欧美xxxx黑人xx丫x性爽| av在线蜜桃| 97超级碰碰碰精品色视频在线观看| 午夜福利欧美成人| 久久久精品大字幕| 偷拍熟女少妇极品色| 网址你懂的国产日韩在线| 亚洲成av人片免费观看| 最好的美女福利视频网| 欧美中文综合在线视频| 亚洲精品美女久久久久99蜜臀| 亚洲中文日韩欧美视频| 色播亚洲综合网| 999久久久精品免费观看国产| 亚洲人成伊人成综合网2020| 亚洲片人在线观看| 欧美最黄视频在线播放免费| 性色av乱码一区二区三区2| 亚洲专区国产一区二区| 亚洲成人久久爱视频| 一级a爱片免费观看的视频| 欧美大码av| 成人特级av手机在线观看| 少妇熟女aⅴ在线视频| 男人的好看免费观看在线视频| 嫩草影视91久久| 亚洲av电影不卡..在线观看| 一二三四社区在线视频社区8| 在线观看舔阴道视频| a级毛片在线看网站| 一区二区三区高清视频在线| 99久久综合精品五月天人人| 五月玫瑰六月丁香| 亚洲专区字幕在线| 嫁个100分男人电影在线观看| 午夜福利高清视频| 精品99又大又爽又粗少妇毛片 | 日本熟妇午夜| 成年人黄色毛片网站| 人人妻人人澡欧美一区二区| 精品久久久久久久久久久久久| 国产亚洲欧美在线一区二区| 白带黄色成豆腐渣| 亚洲人成电影免费在线| 男人舔女人下体高潮全视频| 成人一区二区视频在线观看| 成人午夜高清在线视频| 麻豆久久精品国产亚洲av| 色尼玛亚洲综合影院| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| 黄片大片在线免费观看| 欧美日本亚洲视频在线播放| 制服人妻中文乱码| 叶爱在线成人免费视频播放| 午夜两性在线视频| 成年免费大片在线观看| 黄片小视频在线播放| 国产亚洲精品av在线| 免费看十八禁软件| 国产av在哪里看| 午夜久久久久精精品| a级毛片a级免费在线| 黄色日韩在线| 久久这里只有精品19| 九九久久精品国产亚洲av麻豆 | 深夜精品福利| 两性午夜刺激爽爽歪歪视频在线观看| 久久午夜亚洲精品久久| 巨乳人妻的诱惑在线观看| 亚洲欧美日韩高清在线视频| 国产高清有码在线观看视频| 成人av一区二区三区在线看| 欧美日韩亚洲国产一区二区在线观看| 床上黄色一级片| 成人国产综合亚洲| bbb黄色大片| 男女做爰动态图高潮gif福利片| 国产伦精品一区二区三区四那| 日日干狠狠操夜夜爽| 网址你懂的国产日韩在线| bbb黄色大片| 九九久久精品国产亚洲av麻豆 | 十八禁网站免费在线| 天天一区二区日本电影三级| 亚洲国产色片| 国产主播在线观看一区二区| 丁香六月欧美| 一个人观看的视频www高清免费观看 | 亚洲精品粉嫩美女一区| 精品福利观看| 国产成人av教育| 国产精品久久久久久精品电影| 91在线观看av| 国产亚洲精品久久久com| 国产成人系列免费观看| 99精品在免费线老司机午夜| 亚洲国产欧美一区二区综合| av天堂在线播放| 男女那种视频在线观看| 日本a在线网址| 国产精品综合久久久久久久免费| 午夜成年电影在线免费观看| 最近视频中文字幕2019在线8| 在线a可以看的网站| 色视频www国产| 中文字幕人妻丝袜一区二区| 不卡av一区二区三区| 在线观看日韩欧美| 黄色女人牲交| 好看av亚洲va欧美ⅴa在| 午夜精品在线福利| 男人舔女人下体高潮全视频| 久久久久久久久中文| 一本久久中文字幕| 97人妻精品一区二区三区麻豆| 欧美最黄视频在线播放免费| 嫩草影院入口| 欧美日本亚洲视频在线播放| 精品国产乱子伦一区二区三区| 日韩av在线大香蕉| 久久久久亚洲av毛片大全| 国产高清三级在线| 18禁美女被吸乳视频| 12—13女人毛片做爰片一| 欧美国产日韩亚洲一区| 亚洲狠狠婷婷综合久久图片| 精品一区二区三区视频在线观看免费| 精品日产1卡2卡| 综合色av麻豆| 精品国产超薄肉色丝袜足j| 欧美极品一区二区三区四区| 亚洲av电影在线进入| 国产极品精品免费视频能看的| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久大精品| 午夜福利在线观看免费完整高清在 | 91久久精品国产一区二区成人 | 久久99热这里只有精品18| 国产午夜精品论理片| 久久久国产成人免费| 国产精品久久久人人做人人爽| 久久精品亚洲精品国产色婷小说| 国产真实乱freesex| 免费av毛片视频| 亚洲欧美日韩东京热| 热99re8久久精品国产| 深夜精品福利| 一个人看的www免费观看视频| 日本黄色片子视频| 久久香蕉国产精品| 日韩中文字幕欧美一区二区| 欧美日韩亚洲国产一区二区在线观看| 嫩草影院精品99|