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

    YALMIP工具箱在運籌學實驗教學中的應(yīng)用

    2017-09-03 06:09:04俞武揚
    實驗室研究與探索 2017年8期
    關(guān)鍵詞:運籌學工具箱背包

    俞武揚

    ( 杭州電子科技大學 管理學院,杭州 310018)

    YALMIP工具箱在運籌學實驗教學中的應(yīng)用

    俞武揚

    ( 杭州電子科技大學 管理學院,杭州 310018)

    Matlab平臺對于運籌學中各種經(jīng)典優(yōu)化模型的實驗教學具有重要作用,然而由于Matlab工具箱中函數(shù)的調(diào)用格式限制條件非常嚴格,從而導(dǎo)致將實驗數(shù)據(jù)表示成Matlab所需矩陣形式時容易造成實驗教學效率低下的情況。借助于YALMIP工具箱結(jié)合Matlab軟件實現(xiàn)了運籌學實驗教學中經(jīng)典模型的求解,對線性規(guī)劃、運輸問題、背包問題、最短路問題等等結(jié)合具體例子說明了YALMIP語言的建模方法。利用YALMIP的統(tǒng)一建模語言,很容易對各種運籌學典型問題的模型加以實現(xiàn),從而降低了利用軟件解決模型問題的難度,同時也對解決運籌學實驗中的其它問題有借鑒作用。

    運籌學; 實驗教學; Matlab; YALMIP工具箱

    0 引 言

    線性規(guī)劃、運輸問題、整數(shù)規(guī)劃以及圖論中的最小樹、最短路、最大流等問題都是運籌學中的經(jīng)典問題[1-5],在運籌學教學特別是運籌學實驗中需要學生掌握這些經(jīng)典問題的軟件求解[6-8]。目前國內(nèi)高校在運籌學實驗中最常見的工具主要有Excel、WinQSB、Matlab和Lingo軟件。其中Excel功能相對比較簡單,適用于小型問題的求解演示[9]。WinQSB主要是采用模塊化的程序進行演示,可擴展性較差且在國內(nèi)應(yīng)用相對較少[10]。Matlab軟件功能非常強大,對于一些工科類學科而言具有不可替代的地位,然而在運籌學教學過程中由于Matlab軟件對于輸入約束條件都要求采用矩陣乘積形式,這對于一般形式的優(yōu)化模型造成了不小的轉(zhuǎn)換困難[11-12]。而Lingo軟件在功能方面比較單一,對于優(yōu)化結(jié)果的圖形化顯示方面不如Matlab強大[13-14]。本文介紹在Matlab軟件平臺上開發(fā)的用于求解最優(yōu)化問題的Yalmip工具箱,它不僅可以利用Matlab強大的計算能力,調(diào)用Matlab中各種工具箱中的函數(shù),而且對其他優(yōu)化求解器(如Cplex和Gurobi等)提供了一種統(tǒng)一的建模調(diào)用方案。

    利用YALMIP工具箱結(jié)合Matlab軟件建立相應(yīng)的模型在運籌學實驗教學過程中是一種新的嘗試。本文闡述了利用Yalmip工具箱求解運籌學經(jīng)典問題的方法,并結(jié)合實例說明了Yalmip工具箱在運籌學實驗教學中的具體應(yīng)用。

    1 YALMIP簡介

    YALMIP是Lofberg開發(fā)的一個免費Matlab工具箱,它提供了關(guān)于凸優(yōu)化與非凸優(yōu)化問題的一種高級建模求解語言。YALMIP不但包含基本的線性規(guī)劃求解算法,如linprog(線性規(guī)劃)、bintprog(二值線性規(guī)劃)、bnb(分支定界算法)等,它還提供了對Cplex、Gurobi、Glpk、Lpsolve等求解工具箱的包裝。通常,不同的求解器對于優(yōu)化問題的描述并不一致,在使用時需要較多的學習成本并且容易出錯。而Yalmip真正實現(xiàn)了建模和算法二者的分離,提供了一種簡單而統(tǒng)一的建模語言和編程接口,以實現(xiàn)不同求解器之間的集成。因此,我們只需要知道在Matlab下如何用Yalmip的方式建模,而不需要針對每一種工具箱學習新的建模語法。

    Yalmip工具箱可以通過下載網(wǎng)址http://users.isy.liu.se/johanl/yalmip/pmwiki.php?n=Tutorials.Installation獲得。下載安裝包并解壓后,通過如下步驟在Matlab軟件中設(shè)置引用路徑:① 點擊File中的Set Path菜單;② 點擊Add with subfolders菜單,找到解壓好的文件夾Yalmip;③ 點擊確定、保存。對于外部的求解器,可采用類似方式設(shè)置Matlab路徑。

    Yalmip的建模語法簡單到只有四個最基本的命令:

    (1) 創(chuàng)建決策變量。

    >>x=sdpvar(m,n,[option]): 創(chuàng)建m*n的連續(xù)型決策變量矩陣,option是對矩陣的一些參數(shù)指定;

    >>x=intvar(m,n,[option]): 創(chuàng)建m*n的整數(shù)型決策變量矩陣;

    >>x=binvar(m,n,[option]): 創(chuàng)建m*n的0-1型決策變量矩陣。

    注:若決策變量x是n*n的方陣,如果沒有對決策變量進行參數(shù)指定,則默認x為對稱矩陣,否則需要加以參數(shù)指定,即以x=sdpvar(n,n,’full’)的形式定義,整數(shù)型與0-1型類似。

    (2) 約束條件設(shè)置。

    >>F=set(constraint,[tag]): 創(chuàng)建一個以constraint指定的約束,可選參數(shù)tag可以給該約束指定一個字符串標記。其中constraint的表示極為自然,例如有x1+x2+x3<=10的約束,直接寫成:

    >>x=sdpvar(3,1);

    >>F=set(x(1)+x(2)+x(3)<=10);

    如果要添加多個約束可以直接用+連接:

    >>F=set(constraint1,[tag1])+set(constraint2,[tag2])。

    (3) 參數(shù)配置。

    >>ops=sdpsettings(option1,value1,option2,value2,…);

    例如:>>ops=sdpsettings(‘solver’,’Cplex’,verbose’,2);

    參數(shù)‘solver’指定程序用Cplex求解器(必須已安裝,否則報錯),如果不指定’solver’參數(shù),Yalmip會根據(jù)決策變量類型自動選擇最適合的求解器;’verbose’指定顯示冗余度(冗余度越大,求解過程信息越詳細)。

    (4) 求解。

    >>result=solvesdp(F,f,ops); 求解一個最小化數(shù)學規(guī)劃問題,該問題的目標函數(shù)由f指定,約束由F指定,ops指定求解參數(shù),最后的結(jié)果則存儲在result結(jié)構(gòu)體中,可用double命令顯示。

    >>x_star=double(x);

    >>f_star=double(f);

    2 線性規(guī)劃問題

    線性規(guī)劃是運籌學中最重要的一個分支,幾乎所有的優(yōu)化求解器都可以求解線性規(guī)劃問題,下面先給出線性規(guī)劃問題的數(shù)學模型:

    minZ=CX

    例1 minZ=12x1+5x2+8x3

    與該線性規(guī)劃問題相對應(yīng)的Yalmip程序為:

    C=[12 5 8];

    A=[2 3 1;4 1 5];

    b=[30; 15];

    x=sdpvar(3,1);

    f=C*x;

    F=set(A*x>=b)+set(x>=0);

    result=solvesdp(F,f);

    x_star=double(x)

    f_star=double(f)

    求得結(jié)果為:最優(yōu)解x_star=[0;9.6429;1.0714],最優(yōu)目標函數(shù)值為f_star=56.7857。

    3 運輸問題

    運輸問題模型是運籌學中的經(jīng)典模型之一,運輸問題的一般描述如下:某種物資有m個產(chǎn)地和n個銷地,產(chǎn)地i到銷地j的單位物資運輸費用為cij,各個產(chǎn)地的產(chǎn)量分別為ai,各個銷地的銷量分別為bj,要求在滿足產(chǎn)銷約束條件下使總運輸費用最少的運輸方案。

    例2 現(xiàn)有4個產(chǎn)地,5個銷地,其中cij、ai、bj見表1,求該運輸問題的解。

    表1 運輸參數(shù)表

    相應(yīng)的Yalmip程序為:

    C=[1 3 5 7 13; 6 4 3 14 8; 13 3 1 7 4; 1 10 12 7 11];

    a_i=[40;50;30;80];

    b_j=[10;20;15;18;25];

    x=sdpvar(4,5,'full');

    f=sum(sum(C.*x));

    F=set(sum(x,2)<=a_i)+set(sum(x,1)>=b_j')+

    set(x>=0); %sum(x,2)表示關(guān)于矩陣x按行求和,

    %sum(x,1)表示關(guān)于矩陣x按列求和。

    result=solvesdp(F,f);

    x_star=double(x)

    f_star=double(f)

    求得結(jié)果為:

    4 背包問題

    背包問題是一類重要的整數(shù)規(guī)劃模型,其一般描述如下:設(shè)背包的載重量和容積分別為W和V,第i種物品的重量、體積和數(shù)量分別為wi、vi和ci,該種物品的效用為ei,在背包載重量和容積限制條件下,使得背包裝載物品總效用最大化。

    例3 有10種物品的質(zhì)量wi、體積Vi、數(shù)量ci以及效用ei見表2所示,背包載重量為80,容積為60,求背包裝載物品總效用最大化的方案。

    表2 物品屬性表

    相應(yīng)的Yalmip程序為:

    wi=[17,19,3,19,13,2,6,11,20,20];

    Vi=[2,10,10,5,9,2,5,10,8,10];

    ci=[5,2,4,3,5,4,3,1,5,3];

    ei=[8,1,11,12,9,10,9,5,8,3];

    W=80;

    V=60;

    x=intvar(10,1,'full');

    f=ei*x;

    F=set(wi*x<=W)+set(Vi*x<=V)+

    set(0<=x’<=ci);

    result=solvesdp(F,-f); %最大化目標函數(shù)

    x_star=double(x’)

    f_star=double(f)

    求得結(jié)果為:

    x_star=[1 0 3 1 0 4 3 0 0 0],f_star=120。

    5 指派問題

    指派問題也是一類重要的整數(shù)規(guī)劃問題,其一般描述如下:現(xiàn)要安排n個人去完成n項不同的工作,第i個人完成第j項工作所需時間為cij,要求每人完成一項工作,每項工作只由一人完成,求使得完成所有工作總時間最少的指派方案。

    例4 現(xiàn)有5人完成5項工作所需時間如表3所示,要求每人只做一項工作,每項工作只由一人完成且所用總時間最少的指派方案。

    表3 時間效率表

    相應(yīng)的Yalmip程序為:

    cij=[12 7 9 7 9;8 9 6 6 6;7 17 12 14 9;15 14 6 6 10; 4 10 7 10 9];

    xij=binvar(5,5,'full');

    f=sum(sum(cij.*xij));

    F=set(sum(xij,1)==1)+set(sum(xij,2)==1);

    result=solvesdp(F,f);

    x_star=double(xij)

    f_star=double(f)

    求得結(jié)果為:

    6 最短路問題

    最短路問題的一般提法如下:設(shè)G=(V,E)為連通圖,圖中各邊[vi,vj]有權(quán)l(xiāng)ij(lij=∞表示vi,vj不相鄰),vs,vt為圖中任意兩點,設(shè)μ是G中從vs到vt的一條路定義路μ的權(quán)為μ中所有邊的權(quán)之和,記為ω(μ)。最短路問題就是要從所有從vs到vt的道路中,求一條從vs到vt的所有道路中總權(quán)最小的道路μ0,即

    例5 設(shè)有一批貨物要從V1運到V7,網(wǎng)絡(luò)圖見圖1,邊上的數(shù)字表示該路距離,求最短距離的運輸路線。

    相應(yīng)的Yalmip程序為:

    C=[0 1 4 100 100 100 100; 1 0 2 4 2 5 100; 4 2 0 100 100 1 100; 100 4 100 0 2 100 100; 100 2 1 2 0 3 2; 100 5 1 100 3 0 6; 100 100 100 100 2 6 0];

    x=binvar(7,7,’full’);

    F=set(sum(x(1,:)-sum(x(:,1))==1)+set(sum(x(7,:))-sum(x(:,7))==-1);

    fori=2:6

    F=F+set(sum(x(i,:))-sum(x(:,i))==0);

    end

    f=sum(sum(C.*x));

    result=solvesdp(F,f);

    x_star=double(x)

    f_star=double(f)

    結(jié)果為:x(1,2)=1,x(2,5)=1,x(5,7)=1;相應(yīng)的最短路線為V1→V2→V5→V7,最短距離為5。

    圖1 非線性狀態(tài)觀測器框圖

    7 最大流問題

    設(shè)有向圖G=(V,E),G的每一條邊(vi,vj)上的非負數(shù)cij稱為邊的容量,僅有一個入次為0的點vs稱為發(fā)點(源),一個出次為0的點vt稱為收點(匯),其余的點為中間點,這樣的網(wǎng)絡(luò)G稱為容量網(wǎng)絡(luò),記作:G=(V,E,C)。當每條邊上的實際流量小于等于其容量時,流動才能順利進行。對于給定的容量網(wǎng)絡(luò),如何求出從發(fā)點到收點的最大可行流量就是最大流問題。

    例6 假設(shè)某山區(qū)有6個村莊。村莊與村莊之間雖有公路,但是通行的能力有很大差異,邊上的數(shù)字小則表示通行能力差,反之則表示通行能力強(見圖2)。試求出從村莊1~6之間的最大流。

    圖2

    相應(yīng)的Yalmip程序為:

    C=[0 4 0 0 6 0;0 0 4 4 1 0;0 0 0 2 2 0;0 0 0 0 0 9;0 0 0 6 0 7;0 0 0 0 0 0];

    x=intvar(6,6,'full');

    F=set(sum(x(:,1))==0)+set(sum(x(6,:))==0)+set(0<=x<=C);

    fori=2:5

    F=F+set(sum(x(i,:))-sum(x(:,i))==0);

    end

    f=sum(x(s,:));

    result=solvesdp(F,-f);

    f_star=double(f)

    x_star=double(x)

    結(jié)果為:

    8 最小費用最大流問題

    例7 假設(shè)從倉庫V4到商店V5要運送8個流量的貨物,邊上左邊數(shù)字表示線段最大通過能力,右邊數(shù)字表示通過單位流量的費用(見圖3),試求出費用最小的運輸方案。

    圖3

    相應(yīng)的Yalmip程序為:

    C=[0 0 2 0 7;5 0 10 0 0;0 0 0 0 4;10 8 0 0 0;0 0 0 0 0];

    D=[0 0 6 0 1;2 0 3 0 0;0 0 0 0 2;4 1 0 0 0;0 0 0 0 0];

    x=intvar(5,5,'full');

    F=set(sum(x(:,4))==0)+set(sum(x(5,:))==0)+set(0<=x<=C)+set(sum(x(4,:))==8);

    fori=1:3

    F=F+set(sum(x(i,:))-sum(x(:,i))==0);

    end

    f=sum(sum(D.*x));

    result=solvesdp(F,f);

    f_star=double(f)

    x_star=double(x)

    結(jié)果為:

    9 結(jié) 語

    在運籌學教學過程中結(jié)合實驗教學要求,利用YALMIP工具箱在Matlab環(huán)境中求解問題是一種新的嘗試。由于YALMIP工具箱提供了一種簡潔而統(tǒng)一的建模語法,比Matlab原有語言的表達方法更為易于掌握,可使學生更為關(guān)注問題建模的本質(zhì),而不需要在求解函數(shù)的區(qū)別上分散注意力,因此從運籌學實驗教學的角度來看,結(jié)合YALMIP與Matlab以用于求解運籌學模型值得進一步的推廣。

    [1] 韓大衛(wèi). 管理運籌學[M].6版,大連:大連理工大學出版社,2010.

    [2] 胡運權(quán),郭耀煌. 運籌學教程[M].4版,北京:清華大學出版社,2012.

    [3] 弗雷德里克.S.希利爾. 運籌學導(dǎo)論[M].10版,北京:清華大學出版社,2015.

    [4] 韓伯棠. 管理運籌學[M].4版,北京:高等教育出版社,2015.

    [5] 卜心怡,俞武揚,余福茂,等.管理運籌學[M].北京:電子工業(yè)出版社,2016.

    [6] 夏良玉.案例教學和上機實驗一體化的運籌學教學模式探討[J].科教導(dǎo)刊,2013(14):113-115.

    [7] 趙清俊,陳桂蘭.運籌學實驗軟件在線性規(guī)劃問題教學中的應(yīng)用[J].重慶文理學院學報,2013,32(3):110-112.

    [8] 鄧勝岳,周立前,方四林,等.數(shù)學建模和數(shù)學實驗在《運籌學》課程教學中的應(yīng)用研究[J].湖南理工學院學報(自然科學版),2015,28(1):91-94.

    [9] 于瑛英.EXCEL運籌學規(guī)劃論教學中的應(yīng)用[J].教育教學論壇,2014(10):278-280.

    [10] 鄧 萍,呂俊娜,閆 芳,等.基于QSB軟件的運籌學實驗教學研究[J].課程教育研究(學法教法研究),2015(34):21-21.

    [11] 楊毓玲.基于Matlab的運籌學實驗教學研究[J].科技經(jīng)濟市場,2012(11):115-116.

    [12] 張 明,王文文.Matlab在經(jīng)管類運籌學教學中的探索與實踐[J].大學教育,2012,7(1):81-83.

    [13] 謝金星,薛 毅. 優(yōu)化建模與LindoLingo軟件[M].北京:清華大學出版社,2005.

    [14] 洪 文,朱云鵑,金 震,等. Lingo在運籌學實驗教學中的應(yīng)用[J]. 實驗室研究與探索,2012,31(4):265-270.

    Application of YALMIP in Experiment Teaching of Operations Research

    YU Wuyang

    (Management School, Hangzhou Dianzi University, Hangzhou 310018, China)

    Matlab platform plays an important role in operations research experiment teaching for all kinds of classical optimization models. Since many functions within Matlab toolboxes have very strict format constraints, which result in the poor efficiency of experiment teaching when we represent experimental data into the matrix form of Matlab. Combining the YALMIP toolbox with Matlab can solve the classic models in operations research, such as linear programming, transportation problem, knapsack problem, the most short circuit problem, and so on. Concrete examples are given to illustrate the modeling method of YALMIP language. Using unified modeling language of YALMIP, it is easy to implement these classical optimization models of operations research. The method reduces the difficulty of solving the problem of models by software, and this toolbox can also be used to solve other problems in operations research experiment.

    ooperations research; experiment teaching; Matlab; YALMIP toolbox

    2016-11-22

    杭州電子科技大學“管理科學與工程”重點研究基地項目(ZD06-201601)

    俞武揚(1974-),男,浙江嵊州人,博士,副教授,研究方向:物流建模與優(yōu)化計算。

    Tel.:0571-86878533;E-mail:yu_wuyang@163.com

    O 221.6

    A

    1006-7167(2017)08-0274-05

    猜你喜歡
    運籌學工具箱背包
    大山里的“背包書記”
    一包裝天下 精嘉Alta銳達Sky51D背包體驗
    鼓鼓的背包
    創(chuàng)意西瓜背包
    童話世界(2017年11期)2017-05-17 05:28:26
    Matlab曲線擬合工具箱在地基沉降預(yù)測模型中的應(yīng)用
    運籌學課程教學改革問題研究
    淺談對運籌學專業(yè)教育的一些看法
    山西青年(2016年17期)2016-02-04 21:00:06
    搜狗分號工具箱 輸入更便捷
    小熊的工具箱
    爸爸的工具箱
    娃娃畫報(2009年13期)2009-11-13 03:54:02
    男女下面插进去视频免费观看| 久99久视频精品免费| 婷婷精品国产亚洲av在线| 校园春色视频在线观看| 午夜福利免费观看在线| 免费无遮挡裸体视频| 香蕉国产在线看| 一区福利在线观看| 啦啦啦韩国在线观看视频| 午夜精品久久久久久毛片777| 在线观看免费视频网站a站| 亚洲国产日韩欧美精品在线观看 | 国产精华一区二区三区| 亚洲午夜精品一区,二区,三区| 久久这里只有精品19| 91成人精品电影| 久久午夜综合久久蜜桃| 非洲黑人性xxxx精品又粗又长| 亚洲一码二码三码区别大吗| 十八禁网站免费在线| 欧美乱码精品一区二区三区| 久久久精品国产亚洲av高清涩受| 国产亚洲av嫩草精品影院| 亚洲成人精品中文字幕电影| 多毛熟女@视频| 亚洲精品美女久久久久99蜜臀| 国产免费男女视频| 美女 人体艺术 gogo| 亚洲片人在线观看| 久久精品亚洲熟妇少妇任你| 天堂√8在线中文| 正在播放国产对白刺激| 99国产综合亚洲精品| 韩国精品一区二区三区| 在线播放国产精品三级| 叶爱在线成人免费视频播放| 在线观看66精品国产| 亚洲男人的天堂狠狠| 国产一区二区激情短视频| 亚洲最大成人中文| 两性夫妻黄色片| 精品不卡国产一区二区三区| 高清毛片免费观看视频网站| 亚洲精品国产一区二区精华液| 黄色女人牲交| 熟妇人妻久久中文字幕3abv| 亚洲欧美激情综合另类| 中文字幕人妻丝袜一区二区| 国产成人精品久久二区二区免费| 美女免费视频网站| 黄网站色视频无遮挡免费观看| 女人精品久久久久毛片| 男人操女人黄网站| 一级,二级,三级黄色视频| 嫩草影院精品99| 久久久国产精品麻豆| 黄色视频不卡| 巨乳人妻的诱惑在线观看| 亚洲欧美日韩无卡精品| 久久久久亚洲av毛片大全| 身体一侧抽搐| 99久久久亚洲精品蜜臀av| 成人国产综合亚洲| 久久狼人影院| 18禁美女被吸乳视频| 可以在线观看的亚洲视频| 久久影院123| 欧美不卡视频在线免费观看 | 亚洲五月婷婷丁香| 亚洲第一电影网av| 国语自产精品视频在线第100页| 嫩草影视91久久| 欧美激情 高清一区二区三区| 视频区欧美日本亚洲| 啪啪无遮挡十八禁网站| 性少妇av在线| 久久青草综合色| 亚洲欧美日韩另类电影网站| 韩国精品一区二区三区| 中文字幕人妻丝袜一区二区| 中国美女看黄片| 欧美成人免费av一区二区三区| 国产1区2区3区精品| 亚洲 欧美 日韩 在线 免费| 999久久久国产精品视频| 国产一卡二卡三卡精品| 男女下面插进去视频免费观看| 性少妇av在线| 成人永久免费在线观看视频| 久久人人97超碰香蕉20202| 亚洲,欧美精品.| 一夜夜www| aaaaa片日本免费| 50天的宝宝边吃奶边哭怎么回事| 搡老熟女国产l中国老女人| 亚洲中文字幕日韩| 精品乱码久久久久久99久播| 久久久久久久久久久久大奶| 午夜久久久在线观看| 国产乱人伦免费视频| 午夜免费观看网址| 午夜精品久久久久久毛片777| 波多野结衣巨乳人妻| 午夜免费成人在线视频| 亚洲精华国产精华精| 伦理电影免费视频| 亚洲专区中文字幕在线| 久久久久久久精品吃奶| 亚洲 欧美一区二区三区| 此物有八面人人有两片| 国产精品久久久久久人妻精品电影| 人妻久久中文字幕网| 欧美不卡视频在线免费观看 | 精品高清国产在线一区| 亚洲成人免费电影在线观看| 宅男免费午夜| 日韩高清综合在线| 国产精品一区二区免费欧美| 一区二区三区激情视频| av欧美777| 18禁美女被吸乳视频| 在线av久久热| 久久精品亚洲熟妇少妇任你| 日本 av在线| 一夜夜www| 亚洲人成电影免费在线| 丁香六月欧美| 亚洲 国产 在线| 美女高潮喷水抽搐中文字幕| 成年人黄色毛片网站| 国产人伦9x9x在线观看| 国产精品亚洲一级av第二区| 久久国产精品男人的天堂亚洲| av在线天堂中文字幕| 国产麻豆成人av免费视频| 国产1区2区3区精品| 每晚都被弄得嗷嗷叫到高潮| 欧美激情极品国产一区二区三区| 69精品国产乱码久久久| 中文字幕久久专区| 大型黄色视频在线免费观看| 欧美日韩亚洲综合一区二区三区_| 一本久久中文字幕| 久久热在线av| 老熟妇仑乱视频hdxx| 亚洲av成人av| 午夜精品在线福利| 人人澡人人妻人| 精品一区二区三区视频在线观看免费| 18禁国产床啪视频网站| 如日韩欧美国产精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲三区欧美一区| 日韩欧美三级三区| 男人操女人黄网站| 99久久99久久久精品蜜桃| 中文亚洲av片在线观看爽| 国产91精品成人一区二区三区| 亚洲精华国产精华精| 成人18禁在线播放| 久久精品91蜜桃| 99热只有精品国产| 在线观看午夜福利视频| 亚洲男人的天堂狠狠| 色婷婷久久久亚洲欧美| avwww免费| а√天堂www在线а√下载| 婷婷丁香在线五月| 成年版毛片免费区| 国产精品,欧美在线| 中国美女看黄片| 午夜久久久久精精品| 一级黄色大片毛片| 一进一出好大好爽视频| 免费看十八禁软件| 日韩大尺度精品在线看网址 | 国产精品精品国产色婷婷| 精品欧美国产一区二区三| 亚洲精品国产色婷婷电影| 麻豆成人av在线观看| 久久精品国产亚洲av香蕉五月| 国产精品 国内视频| 中出人妻视频一区二区| 一本综合久久免费| avwww免费| 视频在线观看一区二区三区| 久久热在线av| 嫩草影视91久久| 俄罗斯特黄特色一大片| 多毛熟女@视频| 成人国产综合亚洲| 黄片小视频在线播放| 18美女黄网站色大片免费观看| 久久人妻av系列| 国产午夜福利久久久久久| 欧美日本中文国产一区发布| 精品国产一区二区久久| 啦啦啦观看免费观看视频高清 | 国产精品二区激情视频| 成人手机av| 中文字幕av电影在线播放| 嫩草影视91久久| 国产av精品麻豆| 亚洲一区高清亚洲精品| av电影中文网址| 男女下面进入的视频免费午夜 | 国产高清激情床上av| 电影成人av| 国产成人啪精品午夜网站| 欧美激情久久久久久爽电影 | 国产高清videossex| 日韩精品中文字幕看吧| 午夜福利一区二区在线看| 中亚洲国语对白在线视频| 亚洲精品国产色婷婷电影| 亚洲中文字幕一区二区三区有码在线看 | 久久精品亚洲熟妇少妇任你| 久久欧美精品欧美久久欧美| 国产成人精品无人区| 18禁观看日本| 91国产中文字幕| 亚洲情色 制服丝袜| 国产精品一区二区精品视频观看| videosex国产| 黄色a级毛片大全视频| 色在线成人网| 欧美日韩福利视频一区二区| 一级黄色大片毛片| 18禁美女被吸乳视频| 午夜免费鲁丝| 男女下面进入的视频免费午夜 | 国产成人精品久久二区二区91| 精品国内亚洲2022精品成人| 男女做爰动态图高潮gif福利片 | 此物有八面人人有两片| a在线观看视频网站| 多毛熟女@视频| 99在线人妻在线中文字幕| 久久天躁狠狠躁夜夜2o2o| 亚洲五月婷婷丁香| 久久久久九九精品影院| av超薄肉色丝袜交足视频| 国产精品野战在线观看| 欧美一级a爱片免费观看看 | 亚洲色图综合在线观看| 国产熟女xx| 亚洲国产高清在线一区二区三 | 久久人妻福利社区极品人妻图片| 精品久久久久久成人av| 欧美 亚洲 国产 日韩一| 美女大奶头视频| cao死你这个sao货| 久久人人爽av亚洲精品天堂| 久久久久久免费高清国产稀缺| 亚洲精品美女久久久久99蜜臀| 午夜老司机福利片| 黄片播放在线免费| 欧美最黄视频在线播放免费| 国产一级毛片七仙女欲春2 | 999久久久国产精品视频| 99在线视频只有这里精品首页| 电影成人av| 狂野欧美激情性xxxx| 一级毛片高清免费大全| 美国免费a级毛片| 亚洲中文字幕一区二区三区有码在线看 | 久久精品国产综合久久久| av欧美777| 99riav亚洲国产免费| 久久久久久久久中文| 久久婷婷成人综合色麻豆| 美女高潮喷水抽搐中文字幕| 真人一进一出gif抽搐免费| 午夜福利视频1000在线观看 | 午夜免费激情av| 国产区一区二久久| 久久人人97超碰香蕉20202| www.www免费av| 搞女人的毛片| 久久中文字幕人妻熟女| 久久 成人 亚洲| 91麻豆精品激情在线观看国产| 国产单亲对白刺激| 女性被躁到高潮视频| 国产真人三级小视频在线观看| 99国产精品一区二区三区| 国产精品99久久99久久久不卡| 久久精品国产亚洲av香蕉五月| 国产一区二区在线av高清观看| 免费人成视频x8x8入口观看| 十八禁人妻一区二区| 久久伊人香网站| 欧美日韩瑟瑟在线播放| 涩涩av久久男人的天堂| 大陆偷拍与自拍| 女人被狂操c到高潮| 国产亚洲精品一区二区www| 99国产极品粉嫩在线观看| 黑丝袜美女国产一区| 亚洲免费av在线视频| 一本综合久久免费| 亚洲国产精品999在线| 欧美成人性av电影在线观看| 狂野欧美激情性xxxx| 又黄又爽又免费观看的视频| 欧美在线黄色| 国产精品日韩av在线免费观看 | 国产一级毛片七仙女欲春2 | 国产av一区在线观看免费| 国产成人系列免费观看| 久久久国产成人精品二区| 好男人电影高清在线观看| 成人手机av| 国产精品免费视频内射| 日本 欧美在线| 欧美性长视频在线观看| 免费看美女性在线毛片视频| 搡老熟女国产l中国老女人| 久久久水蜜桃国产精品网| 成人18禁在线播放| 久久婷婷人人爽人人干人人爱 | 99热只有精品国产| 午夜久久久久精精品| 一区福利在线观看| 国产成人啪精品午夜网站| 一级a爱片免费观看的视频| 免费不卡黄色视频| 涩涩av久久男人的天堂| 最新在线观看一区二区三区| 91成人精品电影| 精品国产一区二区久久| 久久久久亚洲av毛片大全| 欧美色视频一区免费| 日韩中文字幕欧美一区二区| 黄片大片在线免费观看| 亚洲 欧美 日韩 在线 免费| 久久九九热精品免费| 国产亚洲av高清不卡| av网站免费在线观看视频| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 色播亚洲综合网| 亚洲欧美一区二区三区黑人| a级毛片在线看网站| 国产97色在线日韩免费| 国产亚洲精品av在线| 91成人精品电影| 久久久久久国产a免费观看| 国产亚洲av嫩草精品影院| 精品一区二区三区av网在线观看| 久久婷婷成人综合色麻豆| 国产亚洲欧美在线一区二区| 午夜老司机福利片| 超碰成人久久| 老司机午夜十八禁免费视频| 91老司机精品| 欧美日本亚洲视频在线播放| 欧美成人性av电影在线观看| 午夜成年电影在线免费观看| 夜夜看夜夜爽夜夜摸| 91成年电影在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 美女国产高潮福利片在线看| 老汉色∧v一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 国产熟女午夜一区二区三区| 成人永久免费在线观看视频| 啦啦啦免费观看视频1| 免费人成视频x8x8入口观看| 九色国产91popny在线| 亚洲少妇的诱惑av| 成人手机av| 老司机在亚洲福利影院| 午夜福利免费观看在线| √禁漫天堂资源中文www| 国产高清视频在线播放一区| 1024香蕉在线观看| 99久久精品国产亚洲精品| 大型黄色视频在线免费观看| x7x7x7水蜜桃| 亚洲成人精品中文字幕电影| 精品日产1卡2卡| 国产不卡一卡二| 美女扒开内裤让男人捅视频| 精品国产美女av久久久久小说| 高清在线国产一区| 亚洲男人天堂网一区| 精品一区二区三区av网在线观看| 国产午夜精品久久久久久| 亚洲男人的天堂狠狠| 久久人人97超碰香蕉20202| 亚洲精品中文字幕一二三四区| 国产成人一区二区三区免费视频网站| 欧美不卡视频在线免费观看 | 啦啦啦观看免费观看视频高清 | 最新在线观看一区二区三区| 一级a爱视频在线免费观看| 精品乱码久久久久久99久播| 国产成人一区二区三区免费视频网站| 一级a爱视频在线免费观看| 欧美日韩亚洲综合一区二区三区_| 精品国产国语对白av| 99国产综合亚洲精品| 无限看片的www在线观看| 亚洲av第一区精品v没综合| 女性被躁到高潮视频| 亚洲黑人精品在线| 亚洲中文字幕日韩| 欧美日韩一级在线毛片| 精品人妻1区二区| 亚洲一区高清亚洲精品| 国产精品一区二区精品视频观看| 国产三级黄色录像| 韩国av一区二区三区四区| 日韩欧美免费精品| 老司机在亚洲福利影院| 大香蕉久久成人网| 精品电影一区二区在线| tocl精华| 97人妻天天添夜夜摸| 久久久久久久久久久久大奶| 亚洲国产高清在线一区二区三 | 一进一出抽搐动态| av天堂在线播放| 国产成人欧美在线观看| 两个人免费观看高清视频| 国产aⅴ精品一区二区三区波| 日本精品一区二区三区蜜桃| 热re99久久国产66热| 日韩欧美免费精品| 黄色成人免费大全| 人成视频在线观看免费观看| 亚洲国产日韩欧美精品在线观看 | 久久久久久久精品吃奶| 变态另类成人亚洲欧美熟女 | av中文乱码字幕在线| 国产97色在线日韩免费| 欧美黑人精品巨大| 精品第一国产精品| 亚洲无线在线观看| 岛国在线观看网站| 黄片大片在线免费观看| 91字幕亚洲| 女人被狂操c到高潮| 国产又爽黄色视频| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区综合在线观看| 国产熟女xx| 国产成人精品无人区| 亚洲欧美精品综合一区二区三区| 国内久久婷婷六月综合欲色啪| avwww免费| 在线观看免费视频网站a站| 欧美日韩亚洲国产一区二区在线观看| 精品国产美女av久久久久小说| 久久久久久久久免费视频了| 90打野战视频偷拍视频| 欧美激情久久久久久爽电影 | 亚洲 国产 在线| 亚洲一区高清亚洲精品| 在线永久观看黄色视频| 一级毛片女人18水好多| 久久中文字幕一级| 麻豆成人av在线观看| 午夜精品久久久久久毛片777| 女人高潮潮喷娇喘18禁视频| √禁漫天堂资源中文www| 国产成年人精品一区二区| 中文字幕高清在线视频| cao死你这个sao货| 色哟哟哟哟哟哟| 97超级碰碰碰精品色视频在线观看| 精品国产一区二区三区四区第35| 黑人欧美特级aaaaaa片| 亚洲精品美女久久av网站| 黄色a级毛片大全视频| 亚洲,欧美精品.| 亚洲情色 制服丝袜| 黄片小视频在线播放| 国产aⅴ精品一区二区三区波| 国产野战对白在线观看| 日韩 欧美 亚洲 中文字幕| 男人操女人黄网站| 黄色成人免费大全| 香蕉丝袜av| 9191精品国产免费久久| 免费高清视频大片| 久久人人爽av亚洲精品天堂| 1024香蕉在线观看| 国产高清videossex| 国产蜜桃级精品一区二区三区| 精品久久久久久成人av| 国产成人系列免费观看| 狂野欧美激情性xxxx| 美女大奶头视频| 国产单亲对白刺激| 久久国产乱子伦精品免费另类| 国产区一区二久久| 91成人精品电影| 日韩欧美在线二视频| 男人舔女人下体高潮全视频| 色综合欧美亚洲国产小说| 国产真人三级小视频在线观看| 精品欧美一区二区三区在线| 美女午夜性视频免费| 成人亚洲精品av一区二区| 国产不卡一卡二| 亚洲精品国产一区二区精华液| 国产又色又爽无遮挡免费看| 一卡2卡三卡四卡精品乱码亚洲| 国产真人三级小视频在线观看| 久久久久久亚洲精品国产蜜桃av| 国产男靠女视频免费网站| 国产97色在线日韩免费| 国产av一区二区精品久久| 国产精品野战在线观看| 日日干狠狠操夜夜爽| 国产一级毛片七仙女欲春2 | 久久午夜综合久久蜜桃| 国产精品自产拍在线观看55亚洲| 久久国产精品影院| 国内精品久久久久久久电影| 丰满的人妻完整版| 久久久久久久久中文| 日本 欧美在线| 久久人人97超碰香蕉20202| 法律面前人人平等表现在哪些方面| 亚洲熟妇熟女久久| 国产精品野战在线观看| 国产精品精品国产色婷婷| 成年版毛片免费区| 日本五十路高清| 国产高清videossex| 精品久久久久久久毛片微露脸| 美国免费a级毛片| 黄色视频,在线免费观看| 国产精品电影一区二区三区| АⅤ资源中文在线天堂| 丝袜美足系列| 免费看十八禁软件| 欧美日本亚洲视频在线播放| xxx96com| 久久久国产精品麻豆| 热re99久久国产66热| 叶爱在线成人免费视频播放| 在线十欧美十亚洲十日本专区| 青草久久国产| 日本欧美视频一区| 成人欧美大片| 色精品久久人妻99蜜桃| 久久久久亚洲av毛片大全| 欧美中文综合在线视频| 女生性感内裤真人,穿戴方法视频| 亚洲最大成人中文| 在线观看免费午夜福利视频| 婷婷六月久久综合丁香| 9191精品国产免费久久| 又大又爽又粗| 天天躁狠狠躁夜夜躁狠狠躁| 人人澡人人妻人| 免费观看人在逋| 色哟哟哟哟哟哟| 欧美日韩亚洲国产一区二区在线观看| 丰满的人妻完整版| 亚洲 欧美一区二区三区| 久久人妻熟女aⅴ| 老司机靠b影院| 免费久久久久久久精品成人欧美视频| 视频区欧美日本亚洲| 一级a爱视频在线免费观看| 老司机午夜十八禁免费视频| 岛国视频午夜一区免费看| 在线十欧美十亚洲十日本专区| 欧美日韩福利视频一区二区| 啦啦啦免费观看视频1| 国产1区2区3区精品| 两个人免费观看高清视频| 麻豆一二三区av精品| 免费看十八禁软件| 亚洲欧美日韩无卡精品| 啦啦啦韩国在线观看视频| 老鸭窝网址在线观看| 婷婷精品国产亚洲av在线| 国产午夜精品久久久久久| www日本在线高清视频| 可以在线观看的亚洲视频| 丰满的人妻完整版| 又黄又粗又硬又大视频| 国产av又大| 国产主播在线观看一区二区| 久久久久久久久中文| 国产精华一区二区三区| 亚洲国产精品999在线| 日日摸夜夜添夜夜添小说| 久久精品成人免费网站| 日韩成人在线观看一区二区三区| 精品一区二区三区av网在线观看| 日韩欧美免费精品| 亚洲专区国产一区二区| 亚洲精华国产精华精| 久久草成人影院| 久久久国产欧美日韩av| 午夜成年电影在线免费观看| 天堂动漫精品| 欧美乱码精品一区二区三区| av网站免费在线观看视频| 激情视频va一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲伊人色综图| 成人亚洲精品av一区二区| 他把我摸到了高潮在线观看| 高潮久久久久久久久久久不卡| 国产av精品麻豆| 97人妻天天添夜夜摸| 真人一进一出gif抽搐免费|