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

    混合型非線性規(guī)劃的MATLAB實(shí)現(xiàn)方法

    2011-04-12 08:02:10沈文勝熊方方金升平余后強(qiáng)
    關(guān)鍵詞:定界剪枝約束條件

    沈文勝 熊方方 金升平 余后強(qiáng)

    (武漢理工大學(xué)理學(xué)院 武漢 430070)

    0 引 言

    最優(yōu)化問(wèn)題廣泛應(yīng)用于工程應(yīng)用的各部門中.在實(shí)際應(yīng)用中經(jīng)常需要求解混合規(guī)劃問(wèn)題,對(duì)于整數(shù)線性規(guī)劃問(wèn)題,文獻(xiàn)[1]給出了其 MATLAB(本文均指 MATLAB7.0)的實(shí)現(xiàn),而對(duì)于求解一般的非線性混合整數(shù)規(guī)劃的問(wèn)題時(shí),在MATLAB中可以調(diào)用一個(gè)基于分枝定界法編寫(xiě)的現(xiàn)成函數(shù)bnb20()[2].但是對(duì)于有既含有整數(shù)又含有離散型的混合非線性規(guī)劃問(wèn)題,例如在機(jī)械工程的齒輪設(shè)計(jì)中涉及到的齒輪的模數(shù)以及齒數(shù)時(shí),模數(shù)就是統(tǒng)一的標(biāo)準(zhǔn)數(shù)值(如系列(GB1357-87))而齒數(shù)則必須是整數(shù)的問(wèn)題,現(xiàn)在的一些數(shù)學(xué)軟件也沒(méi)有可以解決此類問(wèn)題的方法,本文給出一個(gè)在MATLAB中求解此類問(wèn)題的實(shí)現(xiàn)方法.

    考慮一般的混合型非線性規(guī)劃問(wèn)題如式(1)所示.

    式中:x =(x1,x2,…,xn)T為自變量;UB,LB 為向量x的上下界;Aeq,beq分別為線性約束條件中等式約束條件的系數(shù)矩陣和常數(shù)項(xiàng);A,b分別為線性約束條件中不等式約束條件的系數(shù)矩陣和常數(shù)項(xiàng);Ceq(x)和C(x)分別為非線性約束條件中等式約束條件和不等式約束條件;xi1,xi2,…,xir為離散型變量,取值范圍分別為Di1,Di2,…,Dir,xir+1,…,xik為整數(shù)變量;其余變量均為連續(xù)變量.

    分枝定界算法是一種在問(wèn)題的解空間樹(shù)(二叉樹(shù))上搜索問(wèn)題的解的方法.本文利用分枝定界算法對(duì)問(wèn)題的解空間樹(shù)進(jìn)行搜索,所采用的分枝定界法的思想是:將原始問(wèn)題分解,產(chǎn)生一組子問(wèn)題,分枝(branching)是將一組解分為幾組子解,定界(bounding)是建立這些子組解的目標(biāo)函數(shù)的邊界,如果某一子組的解在這些邊界之外,就將這一子組舍棄(剪枝(prune)).

    1 MATLAB實(shí)現(xiàn)中的技術(shù)處理

    搜索算法按搜索的方式分有兩類:一類是深度優(yōu)先搜索(depth first search,DFS),其耗費(fèi)的時(shí)間取決于所采用的存儲(chǔ)結(jié)構(gòu).當(dāng)用二維數(shù)組表示鄰接矩陣作圖的存儲(chǔ)結(jié)構(gòu)時(shí),查找每個(gè)頂點(diǎn)的鄰接點(diǎn)所需時(shí)間為O(n2),其中n為圖中頂點(diǎn)數(shù).而當(dāng)以鄰接表作圖的存儲(chǔ)結(jié)構(gòu)時(shí),查找鄰接點(diǎn)所需的時(shí)間為O(e),其中e為無(wú)向圖中邊的數(shù)或有向圖中弧的數(shù).由此,當(dāng)以鄰接表作存儲(chǔ)結(jié)構(gòu)時(shí),深度優(yōu)先搜索遍歷圖的時(shí)間復(fù)雜度O(n+e);另一類是廣度優(yōu)先搜索(breadth first searchm,BFS)是樹(shù)的按層次遍歷的推廣,遍歷圖的過(guò)程實(shí)質(zhì)上是通過(guò)邊或弧找鄰接點(diǎn)過(guò)程,因此廣度優(yōu)先搜索遍歷圖的時(shí)間復(fù)雜性和深度優(yōu)先搜索遍歷相同,兩者不同之處僅僅在于對(duì)頂點(diǎn)訪問(wèn)的順序不同[3].

    在MATLAB實(shí)現(xiàn)過(guò)程中不需要保留樹(shù)的結(jié)構(gòu),所以本文利用了兩種常見(jiàn)的數(shù)據(jù)結(jié)構(gòu)來(lái)解決變量存儲(chǔ)問(wèn)題:對(duì)于深度優(yōu)先用棧(stack)數(shù)據(jù)結(jié)構(gòu)來(lái)做其存儲(chǔ)結(jié)構(gòu),這是利用棧是按照后進(jìn)先出的原則存儲(chǔ)數(shù)據(jù)的性質(zhì);對(duì)于廣度優(yōu)先法則用隊(duì)列(queue)來(lái)表示,這是利用它只允許在表的前端(front)進(jìn)行刪除操作,而在表的后端(rear)進(jìn)行插入操作的特性.再結(jié)合MATLAB的特點(diǎn),在實(shí)現(xiàn)過(guò)程中來(lái)節(jié)約時(shí)間和空間開(kāi)銷.文獻(xiàn)[4]中例子的搜索樹(shù)如圖1所示,以此樹(shù)為例介紹DFS和BFS兩種遍歷過(guò)程中變量的存儲(chǔ)情況.首先介紹DFS這種情況見(jiàn)表1.其初始狀態(tài)為P0,但是P0中的要左分枝即第1步;P1滿足約束條件是葉節(jié)點(diǎn)需要剪枝,在表1中就是刪除P1,剪枝后P0右分枝得到第2步;P2中的解不滿足約束條件,再左分枝得到P3表1第3步;P3中的解不滿足約束條件,繼續(xù)左分得到P4表1第4步;P4滿足約束條件,P4小于P1,所以P4是葉節(jié)點(diǎn)需要剪枝,剪枝后返回對(duì)P3右分枝得到P5(表1第5步),其他的以此類推.

    圖1 例題的搜索樹(shù)

    表1 DFS遍歷過(guò)程

    BFS具體的實(shí)際情況見(jiàn)表2.初始狀態(tài)是相同,但是分之后存儲(chǔ)的方式不同,每次分枝后其原來(lái)節(jié)點(diǎn)就被刪除,被兩個(gè)分枝節(jié)點(diǎn)所代替,即P0分枝后得到表2第1步,P1葉節(jié)點(diǎn)需要?jiǎng)h除,得到第2步;P2分枝得到第3步;P3分枝后得到第4步.其他的以此類推.

    表2 BFS遍歷過(guò)程

    由此可見(jiàn)經(jīng)過(guò)這種處理,極大的減少了計(jì)算機(jī)的存儲(chǔ)空間.尤其對(duì)于運(yùn)用在大規(guī)模計(jì)算中,其優(yōu)點(diǎn)是顯然的.

    對(duì)于線性問(wèn)題,文獻(xiàn)[4]中認(rèn)為DFS方法優(yōu)于BFS方法,并且給出了3點(diǎn)理由.基于類似的原因,對(duì)于非線性問(wèn)題,建議使用DFS方法,第二部分的設(shè)計(jì)方法就是以DFS為基礎(chǔ)編制的,也可以用BFS來(lái)編寫(xiě),原則上只是訪問(wèn)的順序不同.

    2 基于MATLAB的設(shè)計(jì)方法

    用一個(gè)3×k的數(shù)組xstatus來(lái)設(shè)定離散或者整數(shù)變量的狀態(tài),其中xstatus(1;:)是自變量中離散或整數(shù)變量的下標(biāo);xstatus(2;:)中的每個(gè)元素取值為1或2,1表示xstatus(1;:)中對(duì)應(yīng)列的自變量為整數(shù)型變量,2表示xstatus(1;:)中對(duì)應(yīng)列的變量為離散變量;xstatus(3;:)是xstatus(1;:)中對(duì)應(yīng)列的離散變量的取值范圍的序號(hào),若變量為整數(shù)則為零.于是xstatus就可以表示為

    構(gòu)造節(jié)點(diǎn)的數(shù)據(jù)結(jié)構(gòu)如下.

    其中:‘Left_or_Right’表示左右節(jié)點(diǎn),‘branch_variable_index’表示分枝變量的下標(biāo),‘xlowe’和‘xupper’分別表示分枝變量的下界和上界,‘solution’表示解向量,‘obj_value’表示目標(biāo)函數(shù)的值,‘LB’和‘UB’分別表示該節(jié)點(diǎn)中變量的下界和上界.分枝的過(guò)程就是不斷地修改‘LB’和‘UB’的過(guò)程.

    對(duì)于此類問(wèn)題所采取的策略是:(1)設(shè)定當(dāng)前最好解 MIDP(1).設(shè)定目標(biāo)函數(shù) MIDP(1).obj_value=inf;(2)利用函數(shù)FMINCON()計(jì)算當(dāng)前節(jié)點(diǎn)松弛后的解.使用MATLAB自帶函數(shù)FMINCON()進(jìn)行求解,若EXITFLAG的返回值小于等于零則求解出現(xiàn)錯(cuò)誤,該節(jié)點(diǎn)需要剪枝.若當(dāng)前的目標(biāo)函數(shù)值大于 MIDP(1).obj_value,則該節(jié)點(diǎn)需剪枝;(3)判斷是否符合離散和整數(shù)(checkIntDisc).若是整數(shù)變量則執(zhí)行判斷該變量是否在某個(gè)整數(shù)的δ(δ是計(jì)算時(shí)允許的誤差限,根據(jù)實(shí)際問(wèn)題選取一較小的值)鄰域內(nèi),若在該鄰域內(nèi)則說(shuō)明該值是一可行解,返回到(3)繼續(xù)判斷下一個(gè)變量;若不在該鄰域內(nèi),則對(duì)該變量做向上舍入(ceil)和向下舍入(floor),且將舍入后的值記為該變量的下一次分枝的區(qū)間的上界和下界,返回該值.結(jié)束此次判斷.若是離散變量則將該變量的所有取值范圍與此時(shí)的值進(jìn)行判斷,判斷此變量是否在其取值范圍中的某個(gè)值的δ鄰域內(nèi),若在該鄰域內(nèi)則說(shuō)明該值是一可行解,返回到(3)繼續(xù)判斷下一個(gè)變量.否則得到該變量的離散取值范圍中的上下值作為下一次分枝的上界和下界;(4)對(duì)(3)的結(jié)果進(jìn)行判斷.若符合整數(shù)和離散的要求,則更新 MIDP(1),此時(shí)該節(jié)點(diǎn)需剪枝;否則執(zhí)行(5);(5)對(duì)于節(jié)點(diǎn)(node)進(jìn)行左分枝(left)或者右分枝(right).設(shè)定分枝開(kāi)關(guān),當(dāng) MIDP(index).Left_or_right=1時(shí)向左進(jìn)行分枝,當(dāng)MIDP(index).Left_or_right=2時(shí)進(jìn)行右分枝.具體方法如下:設(shè)定MIDP(index).Left_or_right=1時(shí)更新自變量的取值區(qū)間,將上一節(jié)點(diǎn)取值區(qū)間的下界記為這一節(jié)點(diǎn)的取值區(qū)間的上界;否則 MIDP(index).Left_or_right=2向右進(jìn)行分枝,將上一節(jié)點(diǎn)取值區(qū)間的上界記為這一節(jié)點(diǎn)的取值區(qū)間的下界,返回(2).

    3 應(yīng)用舉例

    3.1 一個(gè)實(shí)例

    給定一個(gè)求解某齒輪的問(wèn)題,自變量為x=[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14].式中:x1,x2,x3為模數(shù)且取值只能為D1中的數(shù);x4,x5,x6,x7,x8,x9為齒輪的齒數(shù),故取值只能為整數(shù);x10,x11為螺旋角的弧度;x12,x13,x14為嚙合齒齒寬.并且x的上下界和D1分別為選取本問(wèn)題的初始值為

    在MATLAB中運(yùn)用DFS編制的程序來(lái)計(jì)算,計(jì)算過(guò)程如下:在MATLAB中輸入A=[];b=[];Aeq=[];beq=[];xstatus=[1,2,3,4,5,6,7,8,9;2,2,2,1,1,1,1,1,1;1,1,1,0,0,0,0,0,0];以及LB,UB,D1,X0后,運(yùn)行本方法所編軟件,經(jīng)過(guò)約4 min的計(jì)算,得出解向量和目標(biāo)函數(shù)值分別為:[2.75,2.5,6,28,85,11,74,12,53,0.254 8,0.254 8,31.052 5,31.157 8,74.778 7]和1.157 8×107.

    3.2 應(yīng)用于文獻(xiàn)[5]問(wèn)題

    對(duì)于優(yōu)化問(wèn)題,工程人員常用的方法是把先忽略自變量的離散或整數(shù)約束要求,當(dāng)作連續(xù)型的非線性規(guī)劃問(wèn)題進(jìn)行求解,然后將得到的結(jié)果圓整[5],即將自變量選取和它最近的離散值或整數(shù)值.這在數(shù)學(xué)中是沒(méi)有理論依據(jù)的,進(jìn)行圓整所得到的結(jié)果很有可能不滿足約束條件.

    對(duì)于文獻(xiàn)[5]中的目標(biāo)函數(shù)及其約束條件,運(yùn)用本文方法和所開(kāi)發(fā)的程序進(jìn)行求解,設(shè)置參數(shù)如下:D1=[2,2.5,3,4,5];D2=[3,4,5,6];xstatus=[1,2,3,4;2,2,1,1;1,2,0,0].可得到如下解向量和目標(biāo)函數(shù)的值:解向量為[2,4,19,17,5.816 9,0.990 3];目標(biāo)函數(shù)值為351.054 7.

    文獻(xiàn)[5]的解向量和目標(biāo)函數(shù)的值分別是[2.005 0,4.004 2,14.001 7,16.376 6,5.8,0.99]和346.259,其圓整后結(jié)果解向量為[2,4,15,17,6,0.990 3],目標(biāo)函數(shù)值是350.

    經(jīng)過(guò)實(shí)際驗(yàn)算,還發(fā)現(xiàn)文獻(xiàn)[5-7]圓整后的解向量不滿足第一和第二個(gè)約束條件,是不可行解,說(shuō)明其圓整的方法對(duì)所述問(wèn)題是行不通的.而本文的方法勿需"圓整",簡(jiǎn)單易行,且求出的解為最優(yōu)解.

    4 結(jié)束語(yǔ)

    對(duì)于變量中既含離散約束、又含整數(shù)約束的混合非線性規(guī)劃問(wèn)題,本文依據(jù)分枝定界原理給出了MATLAB中的實(shí)現(xiàn)方法,在其實(shí)現(xiàn)過(guò)程中,在樹(shù)的遍歷上運(yùn)用了兩種存儲(chǔ)技巧,減少了運(yùn)行時(shí)間和存儲(chǔ)空間,開(kāi)發(fā)了相應(yīng)的軟件,實(shí)例表明方法的正確性和軟件的可行性.

    [1]潘 君.整數(shù)規(guī)劃的分枝定界法及其MATLAB實(shí)現(xiàn)[J].科技信息,2008,(7):167-168.

    [2]薛定宇,陳陽(yáng)泉.高等應(yīng)用數(shù)學(xué)問(wèn)題的 MATALAB求解[M].北京:清華大學(xué)出版社,2004.

    [3]嚴(yán)蔚敏,吳偉民.數(shù)據(jù)結(jié)構(gòu)[M].北京:清華大學(xué)出版社,1997.

    [4]Vanderbei R J.Linear programming:foundations and extensions[M].Boston:Kluwer Academic Publisher,2001.

    [5]崔樹(shù)平,趙 亮,崔 涵.基于 MATLAB最小體積齒輪減速器的優(yōu)化設(shè)計(jì)[J].機(jī)械管理開(kāi)發(fā),2008,23(6):54-55.

    [6]金全意,孟 航.基于 MATLAB的圓柱齒輪減速器優(yōu)化設(shè)計(jì)[J].遼寧工程技術(shù)大學(xué)學(xué)報(bào):自然科學(xué)版,2010,29(z1):55-59.

    [7]張慧鵬.基于MATLAB的二級(jí)圓柱齒輪減速器優(yōu)化設(shè)計(jì) [J].Machinery Design & Manufacture,2010(4):101-106.

    猜你喜歡
    定界剪枝約束條件
    基于一種改進(jìn)AZSVPWM的滿調(diào)制度死區(qū)約束條件分析
    人到晚年宜“剪枝”
    RTK技術(shù)在土地勘測(cè)定界中的應(yīng)用研究
    基于YOLOv4-Tiny模型剪枝算法
    一類DC規(guī)劃問(wèn)題的分支定界算法
    A literature review of research exploring the experiences of overseas nurses in the United Kingdom (2002–2017)
    基于外定界橢球集員估計(jì)的純方位目標(biāo)跟蹤
    剪枝
    線性規(guī)劃的八大妙用
    一種面向不平衡數(shù)據(jù)分類的組合剪枝方法
    亚洲精华国产精华精| av黄色大香蕉| 欧美日韩国产亚洲二区| 91久久精品电影网| 丰满的人妻完整版| 91久久精品国产一区二区成人 | 国产欧美日韩精品一区二区| av在线蜜桃| 国产免费av片在线观看野外av| 欧美性猛交╳xxx乱大交人| 国模一区二区三区四区视频| 真人一进一出gif抽搐免费| 亚洲人成电影免费在线| 国产精品美女特级片免费视频播放器| 亚洲人与动物交配视频| 欧美日韩一级在线毛片| 免费观看的影片在线观看| 嫁个100分男人电影在线观看| 岛国在线免费视频观看| 母亲3免费完整高清在线观看| 色精品久久人妻99蜜桃| 欧美成人a在线观看| 午夜免费激情av| 麻豆国产av国片精品| 久久久久久久午夜电影| 亚洲国产精品sss在线观看| 有码 亚洲区| 男人的好看免费观看在线视频| 亚洲男人的天堂狠狠| 欧美乱码精品一区二区三区| 婷婷精品国产亚洲av在线| 日韩欧美三级三区| 最近在线观看免费完整版| 日本一本二区三区精品| 日本成人三级电影网站| 国产精品香港三级国产av潘金莲| 成人午夜高清在线视频| 欧美成人免费av一区二区三区| 91在线观看av| 美女大奶头视频| 亚洲性夜色夜夜综合| 熟女人妻精品中文字幕| 国产熟女xx| 男人舔奶头视频| 琪琪午夜伦伦电影理论片6080| 日韩亚洲欧美综合| 国产激情偷乱视频一区二区| 欧美又色又爽又黄视频| 国产伦人伦偷精品视频| 亚洲第一欧美日韩一区二区三区| 婷婷亚洲欧美| 国产精品一及| 久久性视频一级片| 国产探花在线观看一区二区| 精品国产超薄肉色丝袜足j| 久久久久久久久大av| 亚洲黑人精品在线| 级片在线观看| 精品一区二区三区人妻视频| 免费电影在线观看免费观看| 久久性视频一级片| 国产伦精品一区二区三区四那| 不卡一级毛片| 久久久国产精品麻豆| 国产精品美女特级片免费视频播放器| 国产亚洲精品一区二区www| 日韩欧美国产一区二区入口| 一夜夜www| 性欧美人与动物交配| 国产午夜精品论理片| 亚洲精品在线观看二区| 18禁在线播放成人免费| 法律面前人人平等表现在哪些方面| 国产毛片a区久久久久| 国产一区在线观看成人免费| 午夜老司机福利剧场| 精品人妻一区二区三区麻豆 | 少妇的逼好多水| 亚洲成av人片免费观看| 免费一级毛片在线播放高清视频| 国产精品久久久久久亚洲av鲁大| 亚洲精品美女久久久久99蜜臀| 国内毛片毛片毛片毛片毛片| 极品教师在线免费播放| 宅男免费午夜| 19禁男女啪啪无遮挡网站| 欧美性猛交╳xxx乱大交人| 亚洲五月婷婷丁香| 国产日本99.免费观看| 99国产精品一区二区三区| 国产亚洲精品久久久久久毛片| av中文乱码字幕在线| 色综合亚洲欧美另类图片| 噜噜噜噜噜久久久久久91| 久久精品91无色码中文字幕| 美女高潮的动态| 国产探花在线观看一区二区| 国产毛片a区久久久久| 亚洲精品在线观看二区| 夜夜爽天天搞| 级片在线观看| 国内精品美女久久久久久| 88av欧美| 非洲黑人性xxxx精品又粗又长| 亚洲国产日韩欧美精品在线观看 | 成熟少妇高潮喷水视频| 欧美日韩综合久久久久久 | h日本视频在线播放| 国产精品久久久久久精品电影| 日韩高清综合在线| 国产精品电影一区二区三区| 亚洲美女黄片视频| 欧美日韩乱码在线| 精品99又大又爽又粗少妇毛片 | 亚洲国产日韩欧美精品在线观看 | 国产 一区 欧美 日韩| 国产精品久久久久久久电影 | 国产亚洲精品av在线| 色在线成人网| 色播亚洲综合网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产一级毛片七仙女欲春2| 天堂动漫精品| 成人鲁丝片一二三区免费| www日本在线高清视频| 国内揄拍国产精品人妻在线| 国产精品嫩草影院av在线观看 | av女优亚洲男人天堂| 小蜜桃在线观看免费完整版高清| 五月伊人婷婷丁香| 最近最新免费中文字幕在线| 日日夜夜操网爽| 又粗又爽又猛毛片免费看| 叶爱在线成人免费视频播放| 99热这里只有精品一区| 一区二区三区激情视频| 俄罗斯特黄特色一大片| 丰满人妻熟妇乱又伦精品不卡| 欧美三级亚洲精品| a级毛片a级免费在线| e午夜精品久久久久久久| 欧美精品啪啪一区二区三区| 天堂动漫精品| 欧美一级毛片孕妇| 麻豆成人av在线观看| 亚洲精品乱码久久久v下载方式 | 亚洲精品在线观看二区| 亚洲自拍偷在线| 亚洲精品美女久久久久99蜜臀| 99久久无色码亚洲精品果冻| 午夜激情福利司机影院| 天堂网av新在线| 国产亚洲精品一区二区www| 1000部很黄的大片| 久久亚洲精品不卡| 观看美女的网站| 精品欧美国产一区二区三| 色尼玛亚洲综合影院| 久久精品国产亚洲av涩爱 | 日日摸夜夜添夜夜添小说| 女同久久另类99精品国产91| 亚洲精品成人久久久久久| 国产一区二区亚洲精品在线观看| 毛片女人毛片| www.熟女人妻精品国产| 国产爱豆传媒在线观看| 精品免费久久久久久久清纯| 国内精品久久久久精免费| 午夜免费成人在线视频| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美国产在线观看| 人人妻人人澡欧美一区二区| 最后的刺客免费高清国语| 啦啦啦韩国在线观看视频| 最好的美女福利视频网| tocl精华| 色老头精品视频在线观看| 精品乱码久久久久久99久播| 国产精品亚洲av一区麻豆| av福利片在线观看| 嫩草影院精品99| 精品久久久久久成人av| 精品电影一区二区在线| 日本 欧美在线| 国内精品久久久久精免费| 亚洲精品一区av在线观看| 99国产精品一区二区蜜桃av| 国产欧美日韩精品一区二区| 国产精品久久久人人做人人爽| 日韩亚洲欧美综合| 网址你懂的国产日韩在线| 99在线人妻在线中文字幕| 国产97色在线日韩免费| 不卡一级毛片| 天天躁日日操中文字幕| 久久精品国产亚洲av涩爱 | www.www免费av| 午夜视频国产福利| 国内揄拍国产精品人妻在线| 国内精品久久久久久久电影| 日韩中文字幕欧美一区二区| 精品免费久久久久久久清纯| 脱女人内裤的视频| 婷婷精品国产亚洲av在线| 亚洲精华国产精华精| 免费一级毛片在线播放高清视频| 少妇丰满av| 久99久视频精品免费| www日本在线高清视频| 一二三四社区在线视频社区8| 国产成人欧美在线观看| 久久久久久九九精品二区国产| 成人午夜高清在线视频| 午夜日韩欧美国产| 男女做爰动态图高潮gif福利片| 国产真实乱freesex| 日本熟妇午夜| 欧美一区二区国产精品久久精品| 亚洲精品乱码久久久v下载方式 | 黄色丝袜av网址大全| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看| 少妇熟女aⅴ在线视频| 色吧在线观看| 欧美极品一区二区三区四区| 国产三级中文精品| 欧美乱妇无乱码| 国产成+人综合+亚洲专区| 国产精品,欧美在线| xxxwww97欧美| 人妻夜夜爽99麻豆av| 性色av乱码一区二区三区2| 国产色爽女视频免费观看| 国产av在哪里看| 最后的刺客免费高清国语| 韩国av一区二区三区四区| 高清日韩中文字幕在线| 色视频www国产| 精品国产亚洲在线| 久久精品91无色码中文字幕| 99久久99久久久精品蜜桃| 国产私拍福利视频在线观看| 校园春色视频在线观看| 亚洲精品一区av在线观看| 免费观看的影片在线观看| 国产免费男女视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人成网站高清观看| 国产高清激情床上av| 无人区码免费观看不卡| 999久久久精品免费观看国产| 日日干狠狠操夜夜爽| 99精品在免费线老司机午夜| 亚洲美女视频黄频| 欧美中文日本在线观看视频| 日韩欧美国产一区二区入口| 99精品久久久久人妻精品| 国产乱人伦免费视频| 天天添夜夜摸| 日日夜夜操网爽| 国产蜜桃级精品一区二区三区| 欧美色欧美亚洲另类二区| 日本成人三级电影网站| 99久久无色码亚洲精品果冻| 亚洲va日本ⅴa欧美va伊人久久| 老汉色av国产亚洲站长工具| 国产高清三级在线| 一进一出好大好爽视频| 国产久久久一区二区三区| 亚洲美女黄片视频| 一级作爱视频免费观看| 欧美日韩黄片免| 91久久精品国产一区二区成人 | 欧美黑人欧美精品刺激| 日本在线视频免费播放| 日韩 欧美 亚洲 中文字幕| 三级毛片av免费| 欧美日韩综合久久久久久 | 午夜影院日韩av| 久9热在线精品视频| 国语自产精品视频在线第100页| 一个人免费在线观看的高清视频| av专区在线播放| 国产蜜桃级精品一区二区三区| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av在线| 真人一进一出gif抽搐免费| 欧美不卡视频在线免费观看| 51午夜福利影视在线观看| 午夜视频国产福利| 国产精品三级大全| 人妻丰满熟妇av一区二区三区| 久久婷婷人人爽人人干人人爱| 国产视频内射| 此物有八面人人有两片| 亚洲人与动物交配视频| 午夜福利高清视频| 精品熟女少妇八av免费久了| 99久久久亚洲精品蜜臀av| 精品久久久久久久末码| 手机成人av网站| 99国产综合亚洲精品| 国产精品电影一区二区三区| 长腿黑丝高跟| 看免费av毛片| 亚洲精品美女久久久久99蜜臀| www.熟女人妻精品国产| 国产免费av片在线观看野外av| 国产精品久久久久久久久免 | 两性午夜刺激爽爽歪歪视频在线观看| 国内精品久久久久久久电影| 亚洲成av人片免费观看| 亚洲欧美日韩高清在线视频| 亚洲一区二区三区不卡视频| 免费无遮挡裸体视频| 色在线成人网| 搡老妇女老女人老熟妇| 欧洲精品卡2卡3卡4卡5卡区| 九九在线视频观看精品| 九色成人免费人妻av| 日本黄色片子视频| 久久国产精品影院| 国产一区二区在线av高清观看| 亚洲成a人片在线一区二区| 激情在线观看视频在线高清| 99精品欧美一区二区三区四区| bbb黄色大片| 熟妇人妻久久中文字幕3abv| 国产主播在线观看一区二区| 亚洲欧美日韩高清专用| 亚洲国产欧美网| 欧美日本亚洲视频在线播放| 给我免费播放毛片高清在线观看| 久久精品91蜜桃| 在线观看午夜福利视频| 在线a可以看的网站| 日日干狠狠操夜夜爽| 国产伦一二天堂av在线观看| 黄色视频,在线免费观看| 一区二区三区国产精品乱码| 国产精品爽爽va在线观看网站| 成人一区二区视频在线观看| 国产国拍精品亚洲av在线观看 | 999久久久精品免费观看国产| www.www免费av| 国语自产精品视频在线第100页| 日日干狠狠操夜夜爽| 欧美+日韩+精品| 极品教师在线免费播放| 久久精品亚洲精品国产色婷小说| 搡老妇女老女人老熟妇| 欧美激情久久久久久爽电影| 丝袜美腿在线中文| 婷婷精品国产亚洲av| а√天堂www在线а√下载| 无限看片的www在线观看| 俄罗斯特黄特色一大片| 老司机深夜福利视频在线观看| 亚洲av电影在线进入| 看片在线看免费视频| 长腿黑丝高跟| 天堂av国产一区二区熟女人妻| 久久婷婷人人爽人人干人人爱| 中文字幕精品亚洲无线码一区| 麻豆国产97在线/欧美| 国产老妇女一区| 两个人的视频大全免费| 日韩欧美三级三区| 免费在线观看亚洲国产| 亚洲中文日韩欧美视频| 搞女人的毛片| 3wmmmm亚洲av在线观看| 色综合站精品国产| 久久性视频一级片| 可以在线观看毛片的网站| 亚洲精品亚洲一区二区| 国产欧美日韩精品亚洲av| 999久久久精品免费观看国产| 噜噜噜噜噜久久久久久91| 在线播放国产精品三级| 男人舔奶头视频| 国模一区二区三区四区视频| 亚洲午夜理论影院| 国产精品久久久人人做人人爽| 午夜激情福利司机影院| 在线a可以看的网站| av在线天堂中文字幕| 最近最新中文字幕大全电影3| 最好的美女福利视频网| 欧美成狂野欧美在线观看| 色综合婷婷激情| 在线免费观看的www视频| 国产 一区 欧美 日韩| 美女黄网站色视频| 啦啦啦观看免费观看视频高清| 99精品欧美一区二区三区四区| 亚洲在线自拍视频| 变态另类丝袜制服| 中文字幕高清在线视频| 舔av片在线| 亚洲一区二区三区不卡视频| 岛国在线免费视频观看| 在线观看66精品国产| 一区二区三区免费毛片| 国产99白浆流出| 亚洲国产色片| 神马国产精品三级电影在线观看| 美女高潮的动态| 国产成人系列免费观看| 欧美日韩亚洲国产一区二区在线观看| www.色视频.com| 精品人妻一区二区三区麻豆 | 久久久久九九精品影院| 亚洲成人久久性| 亚洲人与动物交配视频| 色在线成人网| 内地一区二区视频在线| 18禁国产床啪视频网站| 国产精品电影一区二区三区| 国产精品精品国产色婷婷| 国产亚洲精品一区二区www| 性色av乱码一区二区三区2| 国产精品国产高清国产av| 夜夜爽天天搞| 亚洲精品在线美女| 美女 人体艺术 gogo| 欧美又色又爽又黄视频| 99久国产av精品| 国产一区二区在线观看日韩 | 日本与韩国留学比较| 国产97色在线日韩免费| 日本 av在线| 欧美大码av| 亚洲av成人av| 91九色精品人成在线观看| 日本免费a在线| 成人国产综合亚洲| 午夜免费激情av| 12—13女人毛片做爰片一| 热99re8久久精品国产| 久久人人精品亚洲av| www国产在线视频色| 国内揄拍国产精品人妻在线| 男女那种视频在线观看| 床上黄色一级片| 国产亚洲精品av在线| 亚洲欧美一区二区三区黑人| 亚洲av五月六月丁香网| 欧美乱妇无乱码| 好男人电影高清在线观看| 美女高潮的动态| av欧美777| 两个人看的免费小视频| 综合色av麻豆| 国产精品一区二区三区四区免费观看 | 成人鲁丝片一二三区免费| 日韩欧美精品v在线| 亚洲av第一区精品v没综合| 欧美激情久久久久久爽电影| 十八禁人妻一区二区| 日韩欧美 国产精品| 国产日本99.免费观看| 精华霜和精华液先用哪个| 欧美国产日韩亚洲一区| 日韩欧美一区二区三区在线观看| 成人鲁丝片一二三区免费| 国产乱人伦免费视频| 精品一区二区三区视频在线观看免费| 精品久久久久久久毛片微露脸| 国产在线精品亚洲第一网站| 男女做爰动态图高潮gif福利片| 国产成人av教育| 欧美三级亚洲精品| 99国产综合亚洲精品| 脱女人内裤的视频| 91av网一区二区| 可以在线观看的亚洲视频| 欧美日韩综合久久久久久 | 日本三级黄在线观看| 啦啦啦免费观看视频1| 嘟嘟电影网在线观看| 2021天堂中文幕一二区在线观| 国产精品久久久久久精品电影小说 | 淫秽高清视频在线观看| 青春草视频在线免费观看| 在线免费观看的www视频| 亚洲久久久久久中文字幕| 婷婷色综合大香蕉| 蜜桃亚洲精品一区二区三区| 免费高清在线观看视频在线观看| 亚洲av免费高清在线观看| 99视频精品全部免费 在线| 六月丁香七月| 欧美不卡视频在线免费观看| 国精品久久久久久国模美| 国产探花极品一区二区| 欧美性感艳星| 18+在线观看网站| 国产亚洲最大av| 免费av观看视频| 色综合站精品国产| 可以在线观看毛片的网站| 日本免费a在线| 熟女电影av网| 国产乱来视频区| 日本wwww免费看| 成人一区二区视频在线观看| 久久99热这里只频精品6学生| 最近中文字幕2019免费版| 国产成人aa在线观看| 久久国产乱子免费精品| 国产一区有黄有色的免费视频 | 亚洲在线观看片| 我要看日韩黄色一级片| 三级经典国产精品| 熟女人妻精品中文字幕| 最新中文字幕久久久久| 成年女人在线观看亚洲视频 | 又爽又黄无遮挡网站| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品 | 国产午夜精品一二区理论片| 在线观看免费高清a一片| 大话2 男鬼变身卡| 国产v大片淫在线免费观看| 精品久久久久久电影网| 老师上课跳d突然被开到最大视频| 日日撸夜夜添| 国产又色又爽无遮挡免| 男插女下体视频免费在线播放| 国产精品1区2区在线观看.| 午夜爱爱视频在线播放| 亚洲av中文av极速乱| 九九爱精品视频在线观看| 床上黄色一级片| 黄色配什么色好看| 一本一本综合久久| 免费观看a级毛片全部| 又黄又爽又刺激的免费视频.| 狂野欧美激情性xxxx在线观看| h日本视频在线播放| 国产美女午夜福利| 亚洲怡红院男人天堂| 波多野结衣巨乳人妻| 日日摸夜夜添夜夜爱| 日韩欧美 国产精品| 亚洲av不卡在线观看| 国产成年人精品一区二区| 伦精品一区二区三区| 大片免费播放器 马上看| 亚洲激情五月婷婷啪啪| 直男gayav资源| 日韩强制内射视频| 精品国产露脸久久av麻豆 | 午夜福利在线观看免费完整高清在| 99热这里只有是精品在线观看| av免费在线看不卡| 欧美+日韩+精品| 国产精品熟女久久久久浪| 97人妻精品一区二区三区麻豆| 精品久久久噜噜| 成人午夜高清在线视频| 午夜福利高清视频| 99久久人妻综合| 国产高清国产精品国产三级 | 亚洲性久久影院| 高清欧美精品videossex| 深夜a级毛片| 国产探花在线观看一区二区| 免费看美女性在线毛片视频| 少妇高潮的动态图| 亚洲va在线va天堂va国产| 亚洲丝袜综合中文字幕| 亚洲av日韩在线播放| 国产精品一区www在线观看| 亚洲精品久久久久久婷婷小说| 久久精品夜色国产| 国产av在哪里看| 亚洲一区高清亚洲精品| 亚洲乱码一区二区免费版| 啦啦啦啦在线视频资源| 欧美bdsm另类| 天堂影院成人在线观看| av又黄又爽大尺度在线免费看| 美女被艹到高潮喷水动态| 搡老乐熟女国产| 日韩国内少妇激情av| 日日摸夜夜添夜夜爱| 你懂的网址亚洲精品在线观看| 国产片特级美女逼逼视频| 国产av在哪里看| 精品久久久久久成人av| 丰满人妻一区二区三区视频av| 成人欧美大片| 我要看日韩黄色一级片| 大香蕉久久网| 亚洲四区av| 国产精品国产三级专区第一集| 国产精品伦人一区二区| 婷婷六月久久综合丁香| 日韩视频在线欧美| 听说在线观看完整版免费高清| 国产毛片a区久久久久| 春色校园在线视频观看| 天天躁夜夜躁狠狠久久av| 日本三级黄在线观看| 久久综合国产亚洲精品| 边亲边吃奶的免费视频| 别揉我奶头 嗯啊视频| 18禁在线播放成人免费| 边亲边吃奶的免费视频| 高清在线视频一区二区三区| 日韩成人av中文字幕在线观看| av网站免费在线观看视频 | 全区人妻精品视频|