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

    一種全局優(yōu)化的兩階段算法

    2012-12-09 07:41:10張玉芬王永軍
    關(guān)鍵詞:全局梯度粒子

    張玉芬,王永軍

    (1.河北大學(xué)數(shù)學(xué)與計(jì)算機(jī)學(xué)院,河北保定 071002;

    2.東方地球物理公司研究院數(shù)據(jù)處理中心,河北涿州 072751)

    一種全局優(yōu)化的兩階段算法

    張玉芬1,王永軍2

    (1.河北大學(xué)數(shù)學(xué)與計(jì)算機(jī)學(xué)院,河北保定 071002;

    2.東方地球物理公司研究院數(shù)據(jù)處理中心,河北涿州 072751)

    為了提高算法的有效性,利用梯度算法和粒子群算法獨(dú)立的運(yùn)行機(jī)制,采用驅(qū)趕技術(shù)和重新初始化部分群體的技術(shù),提出了一種基于梯度下降法和粒子群算法的兩階段優(yōu)化算法,并對新算法進(jìn)行了理論分析和數(shù)值仿真.數(shù)值結(jié)果顯示新算法比單純梯度算法有更好的全局優(yōu)化能力,比單純粒子群算法有更快的收斂速度和更高的精度.新算法求解質(zhì)量更高,運(yùn)行更穩(wěn)定.

    全局優(yōu)化;兩階段算法;梯度算法;粒子群算法

    全局優(yōu)化問題廣泛應(yīng)用于金融模型、網(wǎng)絡(luò)交通、圖像處理、集成電路設(shè)計(jì)、分子生物學(xué)、數(shù)據(jù)庫、物流配送系統(tǒng)設(shè)計(jì)中.因?yàn)檫@些優(yōu)化問題存在多個不同于全局最優(yōu)解的局部最優(yōu)解,所以傳統(tǒng)的非線性規(guī)劃中的局部優(yōu)化方法不能有效地處理這類問題.而且,隨著科學(xué)的發(fā)展,工程中遇到的優(yōu)化問題的規(guī)模越來越大、復(fù)雜性越來越高,這使得全局優(yōu)化問題的求解變得更加困難.全局優(yōu)化問題的難以解決阻滯了許多學(xué)科的發(fā)展,因而全局優(yōu)化方法的研究顯得尤為重要.

    梯度類算法是一種確定性的局部優(yōu)化方法,它能夠快速高效地找到臨近初始點(diǎn)的局部最優(yōu)值.常見的算法有最速下降法、牛頓法、擬牛頓法、共軛梯度方法等.它們已經(jīng)廣泛地應(yīng)用于許多領(lǐng)域[1-5].此類算法的主要缺點(diǎn)是對初始點(diǎn)的嚴(yán)重依賴性、目標(biāo)函數(shù)和約束區(qū)域的拓?fù)浣Y(jié)構(gòu)的嚴(yán)格要求性.

    粒子群算法是一種群體智能優(yōu)化方法[1],在粒子群算法中,群體中的一個粒子根據(jù)它本身找到的最優(yōu)位置以及整個群體當(dāng)前為止找到的最優(yōu)位置,來不斷地調(diào)整它本身的位置.算法的核心思想是用當(dāng)前最好的已知的位置來引導(dǎo)群體朝搜索空間的最優(yōu)位置移動.粒子群算法設(shè)計(jì)簡單,參數(shù)少,容易實(shí)現(xiàn),且表現(xiàn)出了較好的局部跳躍性,已經(jīng)被應(yīng)用到許多領(lǐng)域.但算法有時搜索到的只是局部解,找到高精度的解需要較長的運(yùn)行時間,解決大規(guī)模優(yōu)化問題至今還沒有涉及.

    文獻(xiàn)[2]提出了一種把梯度信息混合到粒子的速度更新公式中加快算法收斂的方法.這種方法沒有使用現(xiàn)存的梯度算法的有效性和群體搜索中找到的最優(yōu)信息.過于強(qiáng)調(diào)啟發(fā),而改變了傳統(tǒng)粒子群算法的核心思想.文獻(xiàn)[6]提出了一種基于梯度方法和動態(tài)隧道法的混合技術(shù)(用GRDT表示)來解決全局優(yōu)化問題.在GRDT方法中,梯度用來加快局部搜索,動態(tài)方程用來尋找好的初始點(diǎn).但是,該方法需要多次積分動態(tài)方程,對高維函數(shù)耗時較多,填充函數(shù)法和隧道法中的參數(shù)較難調(diào)節(jié).

    1 兩階段算法

    1.1 算法原理

    不同于一般意義下的進(jìn)化算法,即局部搜索限制在進(jìn)化算法的框架內(nèi),本文設(shè)計(jì)了一種新的混和算法.在新算法中,粒子群按照其獨(dú)有的進(jìn)化機(jī)制演化,而梯度算法使用目標(biāo)函數(shù)的梯度信息進(jìn)行獨(dú)立的搜索.粒子群算法起著全局搜索的作用,而梯度搜索起著局部探測的作用.從而把2種算法混合為一個兩階段方法,發(fā)揮出它們各自的優(yōu)點(diǎn)以解決復(fù)雜函數(shù)的優(yōu)化問題是本文的核心思想所在.具體而言,第1步,從可行域中隨機(jī)產(chǎn)生的一點(diǎn)出發(fā),使用梯度算法得到目標(biāo)函數(shù)的一個局部最優(yōu)點(diǎn).其次,以得到的最優(yōu)點(diǎn)為均值,使用高斯分布函數(shù)初始化粒子群.然后,在最大的迭代次數(shù)內(nèi),運(yùn)用粒子群算法找到一個更好的點(diǎn).以找到的點(diǎn)作為初始點(diǎn),回到第1步,重新開始梯度搜索.

    另外,為了利用算法已經(jīng)得到的最優(yōu)值信息(所有已經(jīng)找到的局部最優(yōu)點(diǎn))防止群體早熟現(xiàn)象,把一種驅(qū)趕群體技巧[7]以及部分初始化群體方法加入到新算法中.驅(qū)趕技術(shù)利用當(dāng)前點(diǎn)與已經(jīng)找到的局部最優(yōu)點(diǎn)之間的距離,建立起一種屏障,使得新產(chǎn)生的點(diǎn)不要落到已經(jīng)找到的局部最優(yōu)點(diǎn)太近的范圍內(nèi),這可以使粒子向著更有希望的區(qū)域飛行,同時避免重復(fù)搜索.文獻(xiàn)[7]已經(jīng)證明了此方法可以有效地加快群智能算法的收斂.部分初始化群體方法是在群體搜索不能取得全局進(jìn)展時,使用均勻分布函數(shù)產(chǎn)生的部分群體代替當(dāng)前群體中同等規(guī)模的較差個體的一種做法.

    1.2 算法步驟

    假設(shè)群體規(guī)模是NP,n是問題的規(guī)模,M是一個隨機(jī)生成的數(shù)值,M∈{1,2,…,NP}.本文使用的驅(qū)趕技術(shù)見算法1,更多的信息可以參考文獻(xiàn)[5].本文使用的部分初始化群體技術(shù)執(zhí)行如下:根據(jù)當(dāng)前群體的適應(yīng)值進(jìn)行排序,在可行域中按均勻分布隨機(jī)產(chǎn)生M個樣本點(diǎn),執(zhí)行驅(qū)趕技術(shù),然后代替當(dāng)前群體中M個較差的個體,其余NP-M個個體保留在群體中,進(jìn)而形成同等規(guī)模的新群體.新群體繼續(xù)執(zhí)行粒子群算法搜索.新算法記做GRPSO(gradient based particle swarm optimization method),見算法2.

    算法1 驅(qū)趕技術(shù)[5]

    步驟0.輸入算法找到的極小點(diǎn)集S*={X*j(j=1,2,…,m)},這里的m是找到的極小點(diǎn)的個數(shù).給定小生境的半徑rij與驅(qū)趕因子mij.

    步驟1.對群體中的每一個粒子Xi,i=1,2,…,NP.

    步驟1a.計(jì)算距離d ij=|X i-X*j|j=1,2,m.

    步驟1b.對集合S*中的每一個元素X*j,j=1,2,…,m.

    如果d ij<rij,那么令Zij=(X i-X*j)/d ij.然后置X i=Xi+Pij Z ij.

    步驟2.輸出群體Xi,i=1,2,…,NP.

    算法2 提出的新算法(GRPSO)

    步驟0.隨機(jī)生成可行域中的初始點(diǎn)X0,并置局部搜索計(jì)數(shù)器k=0.給定粒子群內(nèi)部的最大迭代次數(shù)Smax和計(jì)數(shù)器Npso=1、調(diào)用粒子群算法的最大次數(shù)SSmax.設(shè)最優(yōu)點(diǎn)集合S*=Φ.

    步驟1.(局部搜索)從Xk出發(fā),應(yīng)用梯度下降法找到可行域中的1個極小點(diǎn)X*k,同時保留到集合S*中.

    步驟2.(使用粒子群算法的全局搜索)

    步驟2a.根據(jù)高斯分布,以X*k為均值,在可行域中初始化群體,并找出最好個體gbest.若f(X*k)<f(gbest),令gbest=X*k,k=1轉(zhuǎn)步驟1.否則,令gbest=X*k,計(jì)數(shù)器SS=1,轉(zhuǎn)入步驟2b.

    步驟2b.當(dāng)最大迭代次數(shù)Smax沒有達(dá)到時

    1)更新pbest,如果當(dāng)前個體適應(yīng)值好于f(pbest),i=1,2,…,NP.

    2)更新gbest,如果當(dāng)前群體中的最優(yōu)個體適應(yīng)值好于f(gbest).

    3)如果f(gbest)<f(X*k),令X*k=gbest,結(jié)束循環(huán).

    4)根據(jù)速度更新公式(1)和位置更新公式(2)更新群體的速度和位置.

    5)把當(dāng)前群體和極小點(diǎn)集S*作為輸入集,執(zhí)行算法1的驅(qū)趕技術(shù)(Algorithm 1).然后,令Npso=Npso+1.

    步驟2c.如果gbest=X*k,SS=SS+1.

    如果SS>SSmax,停止計(jì)算并輸出gbest=X*k.否則,重新初始化M∈{1,2,…,NP}個粒子,并執(zhí)行算法1的驅(qū)趕技術(shù),保證此M個粒子遠(yuǎn)離找到的局部最優(yōu)點(diǎn).令Npso=+1,轉(zhuǎn)入步驟2b.否則,令X*k+1=gbest,k=k+1,轉(zhuǎn)入步驟1.

    2 實(shí)驗(yàn)結(jié)果與分析

    2.1 實(shí)驗(yàn)結(jié)果

    用文獻(xiàn)[3]中已經(jīng)被GRDT測試過的10個函數(shù)來檢驗(yàn)新算法的優(yōu)化性能.首先,把函數(shù)值的計(jì)算次數(shù)設(shè)定為2 100,比較最終解的質(zhì)量,數(shù)值結(jié)果見表1.用新算法GRPSO對10個標(biāo)準(zhǔn)測試函數(shù)進(jìn)行優(yōu)化,把運(yùn)算的結(jié)果與其他的算法,如GRDT(梯度算法和動態(tài)隧道法的混合方法[3])、HGPSO(混合粒子群算法[6])和專門設(shè)計(jì)處理高維函數(shù)的GRADSA(混合梯度和模擬退火技術(shù)的算法[4])做比較.由于解的質(zhì)量、收斂速度(計(jì)算機(jī)的CPU時間或目標(biāo)函數(shù)值的計(jì)算次數(shù))、算法的尋優(yōu)成功率是衡量算法優(yōu)劣的主要標(biāo)準(zhǔn),對10個復(fù)雜函數(shù)[3]的20次獨(dú)立運(yùn)行后這幾個指標(biāo)的平均結(jié)果總結(jié)如表1.然后,在給定的最大函數(shù)值計(jì)算次數(shù)下,比較算法找到給定精度的解的成功率,數(shù)值結(jié)果見表2和3.

    表1 新算法、GRDT以及其他的文獻(xiàn)找到的函數(shù)最優(yōu)值Tab.1 Optimal value of functions by the new algorithm,GRDT method and other algorithms in bibliography

    表2 新算法和其他文獻(xiàn)提到算法的函數(shù)值平均計(jì)算次數(shù)Tab.2 Average numbers of functions by the new algorithm and other algorithms in bibliographies

    表3 新算法和HGPSO方法的成功率Tab.3 Rate of success in computation by the new algorithm and HGPSO method

    2.2 結(jié)果分析

    從表1可以看出,GRPSO可以找到所有測試函數(shù)的最優(yōu)值,其結(jié)果接近理論上的分析值.與GRDT和其他文獻(xiàn)報道的結(jié)果相比,新算法找到的結(jié)果更好一些,比如函數(shù)H3和SH是2個典型的例子.

    從表2和3可以看到,新算法HGPSO運(yùn)行效果更好,盡管它們都是基于梯度搜索的混PSO算法.同時,數(shù)值結(jié)果意味著在同樣給定的計(jì)算時間內(nèi),新算法可以得到比GRPSO更高精度的解.說明兩階段的獨(dú)立搜索更有助于各個算法發(fā)揮其基本功能.

    從表2可以看出,根據(jù)函數(shù)值的計(jì)算次數(shù),新算法優(yōu)化性能高于GRDT和HGPSO方法,尤其對于GP和H3函數(shù)的優(yōu)化效果更明顯.與其他典型方法如PRS(純隨機(jī)搜索)、MS(多起點(diǎn)方法)、SA1(基于微分方程的模擬退火方法)、SA2(基本模擬退火方法)、TS(禁忌搜索)、TA(樹退火)和TT(信賴方法TRUST)相比,新算法的計(jì)算量大大降低.

    總之,梯度算法幫助新算法加快局部收斂,并找到高精度的解.加入了驅(qū)趕技術(shù)和部分初始化群體的粒子群算法增加了群體多樣性,防止了群體早熟,進(jìn)而使新算法能有效地跳出局部最優(yōu).所有這些顯示新算法對低維的復(fù)雜函數(shù)優(yōu)化效果好.

    3 結(jié)論

    提出了一種基于梯度算法和粒子群算法的兩階段混和算法.從有效的局部搜索(使用梯度算法)開始,然后借助進(jìn)化算法(采用粒子群算法)的全局跳躍性質(zhì)探索更有希望區(qū)域,為局部搜索找到更好的出發(fā)點(diǎn).尋找高精度的解由局部搜索完成,因而省卻了演化算法的進(jìn)化時間,粒子群算法的全局搜索信息得以利用.不同于一般框架下的演化算法,該算法有效地利用了梯度算法和粒子群算法獨(dú)立的運(yùn)行機(jī)制.另外,新算法中采用了驅(qū)趕技術(shù)和重新初始化部分群體的技術(shù).數(shù)值結(jié)果顯示新算法比單純梯度算法有更好的全局優(yōu)化能力,比單純粒子群算法有更快的收斂速度、解的精度.

    [1]KENNEDY J,EBERHART R C.Particle swarm optimization[Z].Proceedings of IEEE International Conference on Neural Networks,Perth Australia,1995.

    [2]ROY C P,SINGH Y P,CHANSARKAR R A.Hybridization of gradient descent algorithms with dynamic tunneling methods for global optimization[J].IEEE T Syst Man Cy,2000,30(3):384-390.

    [3]NOEL M M,JANNETT T C.Simulation of a new hybrid particle swarm optimization algorithm[J].System Symposium,2004,32(4):150-153.

    [4]SUGANTHAN P N.Particle swarm optimiser with neighbourhood operator[Z].Proceedings of the IEEE Congress on Evolutionary Computation,Piscataway,NJ,1999.

    [5]PARSOPOULOS K E,VRAHATIS M N.On the computation of all global optimizers through particles warm optimization[J].IEEE Trans on Evolutionary Computation,2004,8(3):211-224.

    [6]劉星寶,蔡自興,王勇.用于全局優(yōu)化問題的混合免疫進(jìn)化算法[J].西安電子科技大學(xué)學(xué)報,2010,37(5):971-980.

    LIU Xingbao,CAI Zixing,WANG Yong.Hybrid immune evolutionary algorithm for global optimization problems[J].Journal of Xidian University,2010,37(5):971-980.

    [7]路克中,王汝傳,章家順.最優(yōu)化問題全局尋優(yōu)的PSO-BFGS混合算法[J].計(jì)算機(jī)應(yīng)用研究,2007,24(5):17-19.

    LU Kezhong,WANG Ruchuan,ZHANG Jiashun.PSO-BFGS algorithm of global optimum for optimization problems[J].Application Research of Computers,2007,24(5):17-19.

    A two-phase algorithm for global optimization

    ZHANG Yu-fen1,WANG Yong-jun2
    (1.College of Mathematics and Computer Science,Hebei University,Baoding 071002,China;2.Date Processing Center,Geophysical Prospecting Research Institute,Zhuozhou 072751,China)

    To enhance effectiveness of algorithm,on the basis of analyzing the independent operating mechanism of both gradient algorithm and particle swarm algorithm,a two-phase optimization algorithm based on gradient descent and particle swarm algorithm is presented;it adopts the driving technique and the re-initialization technique of part of population.Then,the theoretical analysis and numerical simulation about the new algorithm are made.The numerical simulation shows this new algorithm has better global optimization ability than the gradient algorithm,and it has faster convergences speed and lighter solution accuracy than particle swarm algorithm.This new algorithm produces a lighter quality solution and has more stable operation.

    global optimization;two-phase algorithm;gradient algorithm;particle swarm op timization

    O241.3

    A

    1000-1565(2012)02-0207-05

    2011-09-11

    河北省軟科學(xué)研究計(jì)劃項(xiàng)目(11457250);河北省自然科學(xué)基金資助項(xiàng)目(F2009000236)

    張玉芬(1964-),女,河北肅寧人,河北大學(xué)副教授,主要從事不確定信息處理方向研究.E-mail:happy6002_cn@sina.com

    孟素蘭)

    猜你喜歡
    全局梯度粒子
    Cahn-Hilliard-Brinkman系統(tǒng)的全局吸引子
    量子Navier-Stokes方程弱解的全局存在性
    一個改進(jìn)的WYL型三項(xiàng)共軛梯度法
    一種自適應(yīng)Dai-Liao共軛梯度法
    一類扭積形式的梯度近Ricci孤立子
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    落子山東,意在全局
    金橋(2018年4期)2018-09-26 02:24:54
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    新思路:牽一發(fā)動全局
    河南科技(2014年3期)2014-02-27 14:05:45
    www.熟女人妻精品国产| 一级a爱视频在线免费观看| 一进一出抽搐gif免费好疼 | 男女床上黄色一级片免费看| 两性夫妻黄色片| 一级片'在线观看视频| 91九色精品人成在线观看| 动漫黄色视频在线观看| 热99久久久久精品小说推荐| 51午夜福利影视在线观看| 很黄的视频免费| 国产激情欧美一区二区| 国产精品综合久久久久久久免费 | 欧美日韩黄片免| 国产蜜桃级精品一区二区三区 | 国产亚洲一区二区精品| 欧美午夜高清在线| 桃红色精品国产亚洲av| 夜夜躁狠狠躁天天躁| 亚洲熟女毛片儿| 国产精品久久久人人做人人爽| 亚洲九九香蕉| 久久草成人影院| 国产成+人综合+亚洲专区| 乱人伦中国视频| 岛国毛片在线播放| 日本wwww免费看| 亚洲 欧美一区二区三区| 国产精品1区2区在线观看. | 亚洲精品国产一区二区精华液| 精品久久久精品久久久| 亚洲第一av免费看| 欧美精品人与动牲交sv欧美| 国产精品一区二区在线不卡| 精品国产超薄肉色丝袜足j| 国产一卡二卡三卡精品| 日韩大码丰满熟妇| 精品国产乱子伦一区二区三区| 国产亚洲一区二区精品| 五月开心婷婷网| 99精品久久久久人妻精品| 国内久久婷婷六月综合欲色啪| 色婷婷久久久亚洲欧美| 国产不卡一卡二| 男女免费视频国产| www.熟女人妻精品国产| 亚洲欧美日韩另类电影网站| 国产av一区二区精品久久| 欧美在线一区亚洲| 搡老乐熟女国产| 女同久久另类99精品国产91| 亚洲熟妇中文字幕五十中出 | 精品少妇久久久久久888优播| 中文亚洲av片在线观看爽 | 午夜日韩欧美国产| 大型黄色视频在线免费观看| 成在线人永久免费视频| 操出白浆在线播放| 好看av亚洲va欧美ⅴa在| 亚洲少妇的诱惑av| 天天操日日干夜夜撸| 人妻久久中文字幕网| 美女国产高潮福利片在线看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩人妻精品一区2区三区| 国产精品久久电影中文字幕 | 叶爱在线成人免费视频播放| 一级片'在线观看视频| 激情在线观看视频在线高清 | 99精品在免费线老司机午夜| 极品少妇高潮喷水抽搐| 一夜夜www| 国产成人欧美在线观看 | 欧美亚洲 丝袜 人妻 在线| 搡老乐熟女国产| 国产99白浆流出| 国产一区二区三区综合在线观看| av福利片在线| 日韩免费高清中文字幕av| 国产一区有黄有色的免费视频| 欧美日韩乱码在线| 成人三级做爰电影| 亚洲三区欧美一区| 国产三级黄色录像| 久久久久视频综合| 麻豆乱淫一区二区| 国产精品一区二区在线观看99| 日韩一卡2卡3卡4卡2021年| netflix在线观看网站| 午夜免费成人在线视频| 亚洲七黄色美女视频| 国产亚洲欧美精品永久| 麻豆成人av在线观看| 日本精品一区二区三区蜜桃| 日韩有码中文字幕| 免费日韩欧美在线观看| 人人妻人人爽人人添夜夜欢视频| 日本wwww免费看| 午夜福利影视在线免费观看| 久久天堂一区二区三区四区| 麻豆国产av国片精品| 日日夜夜操网爽| 国产深夜福利视频在线观看| 一夜夜www| √禁漫天堂资源中文www| 国产亚洲av高清不卡| 中文字幕人妻丝袜一区二区| 99久久精品国产亚洲精品| 国产在视频线精品| 国产精品免费一区二区三区在线 | 亚洲欧美精品综合一区二区三区| 19禁男女啪啪无遮挡网站| 国产精品秋霞免费鲁丝片| 91在线观看av| 日韩欧美一区二区三区在线观看 | 变态另类成人亚洲欧美熟女 | 午夜福利一区二区在线看| 国产成人影院久久av| 夜夜躁狠狠躁天天躁| 久久人人97超碰香蕉20202| 亚洲人成77777在线视频| 女人被狂操c到高潮| 美女国产高潮福利片在线看| 在线天堂中文资源库| 色婷婷av一区二区三区视频| 在线观看免费视频日本深夜| 露出奶头的视频| 精品久久蜜臀av无| 久久久久久免费高清国产稀缺| 日本黄色日本黄色录像| 亚洲综合色网址| 成年版毛片免费区| 一二三四社区在线视频社区8| 日本一区二区免费在线视频| 又黄又粗又硬又大视频| 人妻一区二区av| 女警被强在线播放| 97人妻天天添夜夜摸| 高清视频免费观看一区二区| 老司机午夜福利在线观看视频| av免费在线观看网站| 亚洲男人天堂网一区| 免费在线观看黄色视频的| 好看av亚洲va欧美ⅴa在| 欧美黑人精品巨大| 美女 人体艺术 gogo| 高清欧美精品videossex| 久久久精品免费免费高清| 一区在线观看完整版| 色综合婷婷激情| 国产xxxxx性猛交| 又大又爽又粗| 99精品在免费线老司机午夜| 91成人精品电影| 欧美在线一区亚洲| 精品高清国产在线一区| 午夜免费观看网址| 女警被强在线播放| 老司机在亚洲福利影院| 午夜福利免费观看在线| 中文字幕最新亚洲高清| 欧美av亚洲av综合av国产av| 亚洲熟女精品中文字幕| 咕卡用的链子| 精品视频人人做人人爽| 后天国语完整版免费观看| 国产亚洲精品一区二区www | 老司机深夜福利视频在线观看| 国产欧美日韩一区二区三| 久9热在线精品视频| 韩国av一区二区三区四区| 久久久国产欧美日韩av| 欧美人与性动交α欧美精品济南到| netflix在线观看网站| 久久精品aⅴ一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 色播在线永久视频| 欧美老熟妇乱子伦牲交| 久久久精品国产亚洲av高清涩受| 天堂俺去俺来也www色官网| 日本黄色日本黄色录像| 国产成人欧美在线观看 | 夜夜夜夜夜久久久久| 午夜激情av网站| 19禁男女啪啪无遮挡网站| 久久久久国产一级毛片高清牌| 麻豆乱淫一区二区| 国产伦人伦偷精品视频| 欧美大码av| 亚洲 国产 在线| 男女下面插进去视频免费观看| 久久久久久久久免费视频了| 国产又爽黄色视频| 国产精品秋霞免费鲁丝片| 国产乱人伦免费视频| 99在线人妻在线中文字幕 | 啦啦啦免费观看视频1| 久久性视频一级片| 激情在线观看视频在线高清 | 欧美精品啪啪一区二区三区| 国产免费现黄频在线看| 国产精品一区二区免费欧美| 午夜福利欧美成人| aaaaa片日本免费| av免费在线观看网站| 日韩熟女老妇一区二区性免费视频| 亚洲综合色网址| 12—13女人毛片做爰片一| 天堂中文最新版在线下载| 黑人欧美特级aaaaaa片| 国产高清激情床上av| 免费在线观看黄色视频的| 亚洲色图av天堂| 欧美日韩中文字幕国产精品一区二区三区 | 又紧又爽又黄一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 久久中文看片网| 日韩一卡2卡3卡4卡2021年| 国产精品亚洲一级av第二区| 亚洲七黄色美女视频| 看免费av毛片| 午夜免费观看网址| 校园春色视频在线观看| 亚洲精品在线观看二区| 中文字幕人妻熟女乱码| 国产视频一区二区在线看| 99国产综合亚洲精品| 9色porny在线观看| av视频免费观看在线观看| 久久这里只有精品19| 一二三四在线观看免费中文在| 国产免费现黄频在线看| 淫妇啪啪啪对白视频| 精品久久久久久,| 制服诱惑二区| 久久香蕉精品热| 免费女性裸体啪啪无遮挡网站| 人人澡人人妻人| 巨乳人妻的诱惑在线观看| 一级作爱视频免费观看| 国产精品一区二区精品视频观看| 90打野战视频偷拍视频| 精品国内亚洲2022精品成人 | 欧美黄色淫秽网站| 日韩免费av在线播放| 久久亚洲真实| 精品久久久精品久久久| 久久久国产精品麻豆| 成熟少妇高潮喷水视频| 老司机亚洲免费影院| 91成人精品电影| 中文字幕另类日韩欧美亚洲嫩草| 在线av久久热| 韩国av一区二区三区四区| 欧美 日韩 精品 国产| 老汉色av国产亚洲站长工具| 在线观看日韩欧美| 久久国产亚洲av麻豆专区| 久久人人97超碰香蕉20202| 黄色视频,在线免费观看| 亚洲欧美激情综合另类| 又黄又粗又硬又大视频| 99re在线观看精品视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美成狂野欧美在线观看| 一区二区三区国产精品乱码| 涩涩av久久男人的天堂| 18禁美女被吸乳视频| 法律面前人人平等表现在哪些方面| 国产精品 欧美亚洲| 中文字幕色久视频| 精品卡一卡二卡四卡免费| 国产无遮挡羞羞视频在线观看| av有码第一页| 国产精品亚洲av一区麻豆| 亚洲一区二区三区不卡视频| 国产欧美日韩一区二区三区在线| 美女午夜性视频免费| 国产精品自产拍在线观看55亚洲 | 午夜激情av网站| 亚洲成av片中文字幕在线观看| 一级作爱视频免费观看| 精品久久久精品久久久| 日韩 欧美 亚洲 中文字幕| 国产黄色免费在线视频| 50天的宝宝边吃奶边哭怎么回事| 好男人电影高清在线观看| 极品教师在线免费播放| 欧洲精品卡2卡3卡4卡5卡区| 久久香蕉精品热| 国产三级黄色录像| 久久草成人影院| 亚洲成人免费电影在线观看| 亚洲精品粉嫩美女一区| www.熟女人妻精品国产| 国产在线一区二区三区精| 少妇裸体淫交视频免费看高清 | 中文字幕av电影在线播放| 宅男免费午夜| 最新美女视频免费是黄的| 久久精品熟女亚洲av麻豆精品| 日韩一卡2卡3卡4卡2021年| 女人被躁到高潮嗷嗷叫费观| 一级作爱视频免费观看| 欧美 日韩 精品 国产| 男人的好看免费观看在线视频 | 久久精品国产亚洲av香蕉五月 | 欧美黄色片欧美黄色片| 欧美大码av| 久久亚洲精品不卡| 国产精品 欧美亚洲| 欧美av亚洲av综合av国产av| 老司机午夜十八禁免费视频| 亚洲 国产 在线| 一区在线观看完整版| 黄网站色视频无遮挡免费观看| 欧美精品人与动牲交sv欧美| 国产精品99久久99久久久不卡| 亚洲欧美一区二区三区黑人| 免费在线观看视频国产中文字幕亚洲| 久久久国产精品麻豆| 国产成人免费观看mmmm| www日本在线高清视频| 免费少妇av软件| 国产又色又爽无遮挡免费看| 在线观看免费日韩欧美大片| 美国免费a级毛片| 少妇被粗大的猛进出69影院| 麻豆国产av国片精品| 女性被躁到高潮视频| 天堂√8在线中文| 欧美日韩亚洲国产一区二区在线观看 | 国产成人精品久久二区二区免费| 国产免费av片在线观看野外av| 国产97色在线日韩免费| 国产精品98久久久久久宅男小说| 欧美老熟妇乱子伦牲交| tocl精华| 精品卡一卡二卡四卡免费| 热99久久久久精品小说推荐| 日本一区二区免费在线视频| 免费在线观看视频国产中文字幕亚洲| 男女床上黄色一级片免费看| 亚洲免费av在线视频| 少妇粗大呻吟视频| 性色av乱码一区二区三区2| 精品久久久久久久久久免费视频 | 亚洲欧美日韩另类电影网站| 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色 | 丰满的人妻完整版| 国产精华一区二区三区| 一本综合久久免费| 一a级毛片在线观看| bbb黄色大片| 午夜两性在线视频| 国产亚洲精品久久久久5区| 人妻丰满熟妇av一区二区三区 | 岛国毛片在线播放| 成人国语在线视频| 动漫黄色视频在线观看| 成人免费观看视频高清| 少妇的丰满在线观看| 欧美成人午夜精品| 天天添夜夜摸| av天堂久久9| 人人妻,人人澡人人爽秒播| 国产成人精品久久二区二区91| 黄色a级毛片大全视频| 老熟妇仑乱视频hdxx| videosex国产| 欧美日韩亚洲国产一区二区在线观看 | 国产精品乱码一区二三区的特点 | 王馨瑶露胸无遮挡在线观看| 国产亚洲精品第一综合不卡| 一区二区三区精品91| 欧美日韩国产mv在线观看视频| 黄色女人牲交| 久9热在线精品视频| 国产成人精品久久二区二区免费| 两个人免费观看高清视频| 国产国语露脸激情在线看| 亚洲一区中文字幕在线| 亚洲精品美女久久久久99蜜臀| 精品久久久久久电影网| 麻豆成人av在线观看| 极品教师在线免费播放| 制服人妻中文乱码| 国产一区有黄有色的免费视频| 国产精品偷伦视频观看了| 亚洲欧美一区二区三区久久| 午夜免费观看网址| 国产极品粉嫩免费观看在线| 国产乱人伦免费视频| 久久香蕉国产精品| 午夜成年电影在线免费观看| 亚洲av成人不卡在线观看播放网| 中国美女看黄片| 国产成+人综合+亚洲专区| 99re6热这里在线精品视频| aaaaa片日本免费| 亚洲精品国产区一区二| 男女下面插进去视频免费观看| 国产99久久九九免费精品| 一区福利在线观看| 无人区码免费观看不卡| 国产乱人伦免费视频| 欧美日韩亚洲综合一区二区三区_| 精品第一国产精品| av在线播放免费不卡| 黄色片一级片一级黄色片| 国产精品免费视频内射| 国产精品一区二区在线观看99| xxx96com| avwww免费| aaaaa片日本免费| 亚洲精品成人av观看孕妇| 黄频高清免费视频| 9191精品国产免费久久| 国产一区二区激情短视频| xxx96com| 18禁国产床啪视频网站| 国产成人啪精品午夜网站| 自线自在国产av| 大香蕉久久成人网| 国产乱人伦免费视频| 免费高清在线观看日韩| 精品卡一卡二卡四卡免费| 久久精品aⅴ一区二区三区四区| 国产精品久久电影中文字幕 | 精品一区二区三区视频在线观看免费 | 12—13女人毛片做爰片一| 免费人成视频x8x8入口观看| 欧美色视频一区免费| 脱女人内裤的视频| 99国产精品免费福利视频| 少妇的丰满在线观看| 又紧又爽又黄一区二区| 人成视频在线观看免费观看| 欧美乱码精品一区二区三区| 久久亚洲精品不卡| 日本wwww免费看| www.自偷自拍.com| av福利片在线| 日韩免费av在线播放| 欧美色视频一区免费| 亚洲精品在线美女| 国产精品久久久久久人妻精品电影| 国产精品电影一区二区三区 | 国产精品.久久久| 制服人妻中文乱码| 一区二区日韩欧美中文字幕| 中文字幕色久视频| 日本五十路高清| 亚洲人成电影免费在线| 国产成人欧美| 波多野结衣一区麻豆| 久久久精品区二区三区| 青草久久国产| 一区二区三区激情视频| 中国美女看黄片| 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲 | 国产一区二区三区综合在线观看| 伦理电影免费视频| 国产av一区二区精品久久| 国产97色在线日韩免费| 国产无遮挡羞羞视频在线观看| 18在线观看网站| 久久国产精品影院| 精品免费久久久久久久清纯 | 免费高清在线观看日韩| 亚洲精品久久午夜乱码| 国产91精品成人一区二区三区| 午夜久久久在线观看| 村上凉子中文字幕在线| 亚洲欧美激情综合另类| 久久久久久人人人人人| 久久久国产成人免费| 欧美精品亚洲一区二区| 久久亚洲真实| 又黄又爽又免费观看的视频| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 国产精品一区二区在线观看99| 看黄色毛片网站| 亚洲熟妇熟女久久| av国产精品久久久久影院| 18禁国产床啪视频网站| 正在播放国产对白刺激| 国产精品国产av在线观看| 色尼玛亚洲综合影院| 一区二区日韩欧美中文字幕| 久久久久久亚洲精品国产蜜桃av| 99国产极品粉嫩在线观看| 日韩一卡2卡3卡4卡2021年| 久久婷婷成人综合色麻豆| 中文欧美无线码| 一边摸一边做爽爽视频免费| 久久精品成人免费网站| 亚洲中文字幕日韩| 超色免费av| 超碰97精品在线观看| 制服诱惑二区| 精品一区二区三区视频在线观看免费 | 天堂俺去俺来也www色官网| 在线观看午夜福利视频| 国产精品久久久久成人av| 亚洲七黄色美女视频| 天天躁日日躁夜夜躁夜夜| 老司机影院毛片| 欧美日韩国产mv在线观看视频| 久久人妻福利社区极品人妻图片| 亚洲av成人一区二区三| 露出奶头的视频| 99久久99久久久精品蜜桃| 国产男女超爽视频在线观看| 美女福利国产在线| videosex国产| 99国产精品99久久久久| 在线天堂中文资源库| 免费一级毛片在线播放高清视频 | 黑人巨大精品欧美一区二区mp4| 久久久国产精品麻豆| 99国产精品99久久久久| 日韩视频一区二区在线观看| 国产在视频线精品| 俄罗斯特黄特色一大片| 老司机亚洲免费影院| 啪啪无遮挡十八禁网站| 午夜激情av网站| 免费观看人在逋| 黄色毛片三级朝国网站| 国产欧美亚洲国产| 国产精品一区二区免费欧美| 黄色女人牲交| 高清在线国产一区| 午夜福利一区二区在线看| 欧美中文综合在线视频| 村上凉子中文字幕在线| 操美女的视频在线观看| 新久久久久国产一级毛片| 精品欧美一区二区三区在线| 成年动漫av网址| 国产成人啪精品午夜网站| 欧美午夜高清在线| 国产在线精品亚洲第一网站| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品乱码一区二三区的特点 | 亚洲欧美激情在线| 男人的好看免费观看在线视频 | 亚洲一区高清亚洲精品| 97人妻天天添夜夜摸| 男人的好看免费观看在线视频 | 亚洲人成77777在线视频| 成人精品一区二区免费| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美一区二区三区久久| 女人久久www免费人成看片| 久久人妻av系列| 他把我摸到了高潮在线观看| 国产精品久久电影中文字幕 | 亚洲黑人精品在线| 久久精品亚洲精品国产色婷小说| 亚洲成人国产一区在线观看| 高潮久久久久久久久久久不卡| 成年人免费黄色播放视频| 人妻久久中文字幕网| 日韩成人在线观看一区二区三区| 操美女的视频在线观看| 久久国产乱子伦精品免费另类| 亚洲人成伊人成综合网2020| 国产日韩欧美亚洲二区| 精品久久久久久电影网| 欧美精品高潮呻吟av久久| 亚洲免费av在线视频| 免费看a级黄色片| 可以免费在线观看a视频的电影网站| 欧美在线黄色| 正在播放国产对白刺激| 人人妻人人爽人人添夜夜欢视频| 丝袜美足系列| 免费观看人在逋| 水蜜桃什么品种好| 亚洲av电影在线进入| 亚洲午夜精品一区,二区,三区| 国产黄色免费在线视频| 国产精品亚洲一级av第二区| 久久九九热精品免费| 国产精品久久久久久精品古装| 夜夜夜夜夜久久久久| 天堂√8在线中文| 亚洲专区国产一区二区| 国产成人欧美在线观看 | 欧美性长视频在线观看| bbb黄色大片| 在线免费观看的www视频| 久久久久国产精品人妻aⅴ院 | 国产淫语在线视频| 黑人操中国人逼视频| 欧美日韩视频精品一区| 精品人妻在线不人妻| av电影中文网址| 午夜精品久久久久久毛片777| 少妇被粗大的猛进出69影院| 美国免费a级毛片| 欧美日韩一级在线毛片| 1024香蕉在线观看| 免费在线观看视频国产中文字幕亚洲| 国产欧美日韩综合在线一区二区| 亚洲精品一二三| 757午夜福利合集在线观看| 精品久久久久久电影网| 亚洲av成人av| 黑人操中国人逼视频|