• <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
    亚洲情色 制服丝袜| 久久这里只有精品19| 欧美日韩中文字幕国产精品一区二区三区 | 欧美日韩精品网址| 国产男靠女视频免费网站| 国产深夜福利视频在线观看| 国产91精品成人一区二区三区| 国产免费现黄频在线看| 国产色视频综合| 99久久精品国产亚洲精品| 老司机午夜福利在线观看视频| 91大片在线观看| 国产亚洲精品一区二区www | 久久精品亚洲av国产电影网| 欧美乱妇无乱码| aaaaa片日本免费| 好男人电影高清在线观看| 757午夜福利合集在线观看| 精品一区二区三区av网在线观看| 国产人伦9x9x在线观看| 高清欧美精品videossex| 老司机靠b影院| 日本欧美视频一区| 久久人妻福利社区极品人妻图片| 欧美人与性动交α欧美精品济南到| 一级片'在线观看视频| 夜夜躁狠狠躁天天躁| 国产免费男女视频| 欧美 日韩 精品 国产| 国产精品香港三级国产av潘金莲| 精品国产一区二区三区久久久樱花| 制服诱惑二区| 国产高清国产精品国产三级| 亚洲国产精品合色在线| 午夜免费成人在线视频| 欧美黑人欧美精品刺激| 久久99一区二区三区| 一区在线观看完整版| 99久久国产精品久久久| 国产成人av激情在线播放| 亚洲欧美日韩高清在线视频| 一级a爱视频在线免费观看| 日韩欧美国产一区二区入口| 日本五十路高清| 欧美激情 高清一区二区三区| 一个人免费在线观看的高清视频| 午夜亚洲福利在线播放| 在线观看免费午夜福利视频| 欧美久久黑人一区二区| 成年版毛片免费区| 亚洲美女黄片视频| 成年人免费黄色播放视频| 亚洲aⅴ乱码一区二区在线播放 | 9191精品国产免费久久| 身体一侧抽搐| 成年动漫av网址| 一个人免费在线观看的高清视频| 制服诱惑二区| 最近最新中文字幕大全电影3 | 一进一出抽搐gif免费好疼 | 免费不卡黄色视频| 成年女人毛片免费观看观看9 | 男女免费视频国产| 涩涩av久久男人的天堂| 国产精品1区2区在线观看. | 成年人黄色毛片网站| 一进一出抽搐gif免费好疼 | 两性午夜刺激爽爽歪歪视频在线观看 | 免费高清在线观看日韩| 久久国产精品男人的天堂亚洲| 乱人伦中国视频| 身体一侧抽搐| 啦啦啦视频在线资源免费观看| 国产成人免费观看mmmm| 大香蕉久久成人网| 日韩欧美免费精品| 国产欧美日韩精品亚洲av| www.自偷自拍.com| 国产一区二区激情短视频| 美女午夜性视频免费| 亚洲欧美激情在线| 久久午夜综合久久蜜桃| 国产熟女午夜一区二区三区| 国产伦人伦偷精品视频| 交换朋友夫妻互换小说| 中国美女看黄片| 国产色视频综合| 久久性视频一级片| 日韩制服丝袜自拍偷拍| 黄色 视频免费看| 美女国产高潮福利片在线看| 亚洲 欧美一区二区三区| 高清在线国产一区| 国产麻豆69| 别揉我奶头~嗯~啊~动态视频| 麻豆国产av国片精品| 桃红色精品国产亚洲av| av天堂久久9| 最近最新中文字幕大全免费视频| 岛国毛片在线播放| 国产精品秋霞免费鲁丝片| 亚洲国产欧美网| 久久午夜综合久久蜜桃| 国产精品 欧美亚洲| 久久久久久亚洲精品国产蜜桃av| 狠狠狠狠99中文字幕| 亚洲精品国产精品久久久不卡| 亚洲国产精品sss在线观看 | 国产欧美日韩精品亚洲av| 日本欧美视频一区| 久久精品熟女亚洲av麻豆精品| 变态另类成人亚洲欧美熟女 | 人妻丰满熟妇av一区二区三区 | 日韩欧美一区二区三区在线观看 | 亚洲综合色网址| 两个人看的免费小视频| 国产精品乱码一区二三区的特点 | 久久这里只有精品19| 国产国语露脸激情在线看| 夜夜夜夜夜久久久久| 免费在线观看黄色视频的| 三级毛片av免费| 男女高潮啪啪啪动态图| 成人影院久久| 日日夜夜操网爽| 久久久久久久久免费视频了| 窝窝影院91人妻| 好看av亚洲va欧美ⅴa在| 欧美不卡视频在线免费观看 | 777米奇影视久久| 成人三级做爰电影| 色尼玛亚洲综合影院| 国产精品乱码一区二三区的特点 | 久久精品91无色码中文字幕| 99国产精品免费福利视频| 精品福利永久在线观看| 波多野结衣一区麻豆| 色尼玛亚洲综合影院| 欧美乱色亚洲激情| 丝袜在线中文字幕| 亚洲成人国产一区在线观看| 亚洲av美国av| 51午夜福利影视在线观看| 黄色视频,在线免费观看| 国产主播在线观看一区二区| 久久久水蜜桃国产精品网| 中文字幕制服av| 搡老岳熟女国产| 亚洲性夜色夜夜综合| 欧美人与性动交α欧美精品济南到| 99热网站在线观看| 可以免费在线观看a视频的电影网站| 狂野欧美激情性xxxx| 大香蕉久久成人网| 精品无人区乱码1区二区| 午夜久久久在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 中文字幕人妻丝袜制服| 国产精品免费一区二区三区在线 | 别揉我奶头~嗯~啊~动态视频| 国产精品久久电影中文字幕 | 欧美日韩av久久| 岛国在线观看网站| 12—13女人毛片做爰片一| 人人妻人人爽人人添夜夜欢视频| 国产一区二区三区综合在线观看| 视频区欧美日本亚洲| 丝袜人妻中文字幕| 黄色片一级片一级黄色片| xxx96com| 亚洲一区中文字幕在线| 亚洲精品在线观看二区| 天堂动漫精品| tocl精华| 成人影院久久| 亚洲一区二区三区不卡视频| 欧美日韩亚洲高清精品| 亚洲第一av免费看| 日本黄色日本黄色录像| 美女高潮到喷水免费观看| 人人妻人人澡人人爽人人夜夜| 99精国产麻豆久久婷婷| 这个男人来自地球电影免费观看| 成年人黄色毛片网站| 18禁裸乳无遮挡动漫免费视频| 午夜久久久在线观看| 亚洲专区字幕在线| 动漫黄色视频在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲第一青青草原| 免费女性裸体啪啪无遮挡网站| 人人妻人人澡人人看| www.熟女人妻精品国产| 99热国产这里只有精品6| 国产亚洲欧美在线一区二区| 亚洲熟妇中文字幕五十中出 | 一本一本久久a久久精品综合妖精| 国产成人欧美| 久久国产亚洲av麻豆专区| 免费在线观看影片大全网站| 免费黄频网站在线观看国产| 亚洲精品粉嫩美女一区| 在线观看免费视频网站a站| 午夜福利免费观看在线| 国产成人av教育| 中出人妻视频一区二区| 亚洲自偷自拍图片 自拍| 色在线成人网| 99精国产麻豆久久婷婷| 精品国产一区二区久久| 看免费av毛片| av电影中文网址| 久久久精品区二区三区| 精品人妻在线不人妻| 精品久久久久久久毛片微露脸| 国产成人影院久久av| 久久这里只有精品19| 一区二区三区国产精品乱码| 在线观看免费视频日本深夜| 80岁老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 一级片'在线观看视频| 欧美乱码精品一区二区三区| 身体一侧抽搐| 大香蕉久久网| 成人黄色视频免费在线看| 91大片在线观看| 亚洲中文字幕日韩| 国产不卡av网站在线观看| 成人手机av| 国产成人av激情在线播放| 天天躁夜夜躁狠狠躁躁| 超色免费av| av电影中文网址| 亚洲成av片中文字幕在线观看| 亚洲精品乱久久久久久| 成人国产一区最新在线观看| 国产精品成人在线| 国产精品一区二区精品视频观看| 亚洲精品中文字幕一二三四区| 亚洲情色 制服丝袜| 好男人电影高清在线观看| 亚洲一区二区三区欧美精品| 午夜91福利影院| 久久中文看片网| 精品欧美一区二区三区在线| 亚洲av成人不卡在线观看播放网| 亚洲人成电影观看| 国产精品久久久人人做人人爽| 69av精品久久久久久| 黄色a级毛片大全视频| 精品视频人人做人人爽| 欧美丝袜亚洲另类 | xxxhd国产人妻xxx| 午夜福利一区二区在线看| 精品无人区乱码1区二区| 天堂俺去俺来也www色官网| 国产伦人伦偷精品视频| 亚洲成人手机| 日本精品一区二区三区蜜桃| 欧美大码av| 99国产极品粉嫩在线观看| 嫁个100分男人电影在线观看| 久久香蕉激情| 不卡一级毛片| 18禁裸乳无遮挡免费网站照片 | 国产一卡二卡三卡精品| 亚洲精品国产一区二区精华液| 亚洲专区中文字幕在线| tocl精华| 变态另类成人亚洲欧美熟女 | 日韩熟女老妇一区二区性免费视频| www.熟女人妻精品国产| 精品人妻在线不人妻| 国产成人免费观看mmmm| 日本五十路高清| 69精品国产乱码久久久| 老汉色∧v一级毛片| 亚洲精华国产精华精| 香蕉久久夜色| 国产精品av久久久久免费| 看免费av毛片| 亚洲av成人av| 亚洲欧美日韩高清在线视频| 18禁国产床啪视频网站| 欧美中文综合在线视频| 美女午夜性视频免费| 色婷婷av一区二区三区视频| 大型黄色视频在线免费观看| 免费观看a级毛片全部| 亚洲专区国产一区二区| 亚洲五月天丁香| av电影中文网址| 成人免费观看视频高清| 91精品国产国语对白视频| 日本精品一区二区三区蜜桃| 国产精品一区二区免费欧美| 免费在线观看完整版高清| 很黄的视频免费| 激情视频va一区二区三区| 日韩人妻精品一区2区三区| 日韩免费高清中文字幕av| av有码第一页| 黄色毛片三级朝国网站| www.999成人在线观看| 免费av中文字幕在线| 亚洲国产毛片av蜜桃av| 一进一出抽搐动态| ponron亚洲| 久久性视频一级片| 超碰97精品在线观看| 丝袜人妻中文字幕| 中文亚洲av片在线观看爽 | 婷婷成人精品国产| 久久性视频一级片| 热99国产精品久久久久久7| 国产av精品麻豆| 99精品在免费线老司机午夜| 欧美激情极品国产一区二区三区| 一级毛片精品| 国产一区二区三区综合在线观看| 国产激情欧美一区二区| 久久精品国产亚洲av香蕉五月 | 丰满的人妻完整版| 国产精品 欧美亚洲| 亚洲国产看品久久| 成人三级做爰电影| 一区福利在线观看| 熟女少妇亚洲综合色aaa.| 亚洲成人免费av在线播放| 午夜老司机福利片| 自线自在国产av| 男人操女人黄网站| 天堂俺去俺来也www色官网| 视频在线观看一区二区三区| 777久久人妻少妇嫩草av网站| 国产不卡一卡二| 在线观看日韩欧美| 日韩制服丝袜自拍偷拍| 最新的欧美精品一区二区| 亚洲精品av麻豆狂野| 国产91精品成人一区二区三区| 亚洲国产中文字幕在线视频| 一区二区日韩欧美中文字幕| 亚洲国产中文字幕在线视频| 亚洲精品久久成人aⅴ小说| 成熟少妇高潮喷水视频| av在线播放免费不卡| 在线观看免费午夜福利视频| 超碰97精品在线观看| 亚洲精华国产精华精| 欧美日韩av久久| 国产有黄有色有爽视频| netflix在线观看网站| 久久久水蜜桃国产精品网| 亚洲欧洲精品一区二区精品久久久| 99re在线观看精品视频| 国产精品欧美亚洲77777| 国产高清视频在线播放一区| 青草久久国产| 老熟妇仑乱视频hdxx| 国产欧美日韩一区二区精品| 丝袜在线中文字幕| 国产国语露脸激情在线看| 精品电影一区二区在线| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久久久毛片 | 亚洲国产欧美日韩在线播放| 国产av又大| 中文亚洲av片在线观看爽 | 99国产极品粉嫩在线观看| 很黄的视频免费| av中文乱码字幕在线| 两人在一起打扑克的视频| 国产深夜福利视频在线观看| 十分钟在线观看高清视频www| 免费黄频网站在线观看国产| 久久久国产一区二区| 国产一区在线观看成人免费| 丰满人妻熟妇乱又伦精品不卡| tube8黄色片| 首页视频小说图片口味搜索| 巨乳人妻的诱惑在线观看| 成熟少妇高潮喷水视频| 在线永久观看黄色视频| 777米奇影视久久| 午夜福利欧美成人| 成人影院久久| 捣出白浆h1v1| 国产99久久九九免费精品| 亚洲精品美女久久av网站| 国产精品 欧美亚洲| 欧美午夜高清在线| 91老司机精品| 大香蕉久久网| 国产在视频线精品| 亚洲黑人精品在线| 天堂俺去俺来也www色官网| 后天国语完整版免费观看| 午夜两性在线视频| 这个男人来自地球电影免费观看| 精品国产亚洲在线| 黄色片一级片一级黄色片| 午夜免费成人在线视频| 亚洲免费av在线视频| 国产高清视频在线播放一区| 两人在一起打扑克的视频| 国产成人系列免费观看| 男女高潮啪啪啪动态图| 精品人妻熟女毛片av久久网站| 看免费av毛片| 久久中文字幕人妻熟女| 欧美国产精品va在线观看不卡| 久久99一区二区三区| 九色亚洲精品在线播放| 嫁个100分男人电影在线观看| 身体一侧抽搐| 日韩一卡2卡3卡4卡2021年| 涩涩av久久男人的天堂| 国产精品久久电影中文字幕 | 国产亚洲欧美精品永久| 亚洲成人国产一区在线观看| 久久 成人 亚洲| 精品国产一区二区三区久久久樱花| 高清欧美精品videossex| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品美女久久久久99蜜臀| 免费观看精品视频网站| 99精国产麻豆久久婷婷| 9热在线视频观看99| 高清毛片免费观看视频网站 | 巨乳人妻的诱惑在线观看| 多毛熟女@视频| 精品国产美女av久久久久小说| av一本久久久久| 欧美黑人欧美精品刺激| 亚洲av成人一区二区三| 叶爱在线成人免费视频播放| 一边摸一边抽搐一进一出视频| 久久久久国内视频| 欧美 日韩 精品 国产| 日本撒尿小便嘘嘘汇集6| av视频免费观看在线观看| 国产人伦9x9x在线观看| 男女床上黄色一级片免费看| 久久精品国产综合久久久| 精品国产一区二区三区久久久樱花| 女性被躁到高潮视频| a级毛片在线看网站| 欧美日韩精品网址| 亚洲一码二码三码区别大吗| 一进一出好大好爽视频| 免费av中文字幕在线| 91精品三级在线观看| 91精品国产国语对白视频| 中文字幕人妻丝袜制服| 国产亚洲av高清不卡| av福利片在线| 亚洲成人免费av在线播放| 午夜成年电影在线免费观看| 99热只有精品国产| 黑人操中国人逼视频| 亚洲aⅴ乱码一区二区在线播放 | 成年人免费黄色播放视频| 免费一级毛片在线播放高清视频 | 亚洲中文av在线| 咕卡用的链子| 亚洲国产欧美一区二区综合| 少妇 在线观看| 别揉我奶头~嗯~啊~动态视频| 精品亚洲成国产av| 成人影院久久| 90打野战视频偷拍视频| 老熟妇乱子伦视频在线观看| 麻豆国产av国片精品| 91成人精品电影| 亚洲国产中文字幕在线视频| 久久久水蜜桃国产精品网| 亚洲专区字幕在线| 黄片大片在线免费观看| 一进一出抽搐gif免费好疼 | a级毛片黄视频| 久久九九热精品免费| 精品人妻1区二区| 夫妻午夜视频| 欧美日本中文国产一区发布| 天天躁日日躁夜夜躁夜夜| 精品福利永久在线观看| 不卡一级毛片| 国产亚洲欧美98| 另类亚洲欧美激情| 黄色怎么调成土黄色| 国产有黄有色有爽视频| 亚洲三区欧美一区| 精品一区二区三区四区五区乱码| 一级a爱片免费观看的视频| 欧美大码av| 精品久久久久久久毛片微露脸| 成年版毛片免费区| 免费观看a级毛片全部| 999久久久国产精品视频| 成熟少妇高潮喷水视频| 一二三四社区在线视频社区8| 久久国产亚洲av麻豆专区| 国产成人精品无人区| 日本黄色视频三级网站网址 | 国产欧美日韩一区二区三| 少妇猛男粗大的猛烈进出视频| 成人国语在线视频| 色在线成人网| 欧美乱妇无乱码| 最新在线观看一区二区三区| 黄色 视频免费看| 久久久精品国产亚洲av高清涩受| 91字幕亚洲| 视频区图区小说| 十八禁高潮呻吟视频| 日本撒尿小便嘘嘘汇集6| 午夜福利视频在线观看免费| 久久久久视频综合| 日韩精品免费视频一区二区三区| 欧美精品av麻豆av| 丰满的人妻完整版| 久久国产精品人妻蜜桃| 热99久久久久精品小说推荐| 亚洲精品国产一区二区精华液| 精品人妻在线不人妻| 久久久久精品国产欧美久久久| 亚洲少妇的诱惑av| 国产97色在线日韩免费| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| 村上凉子中文字幕在线| 热99re8久久精品国产| 日韩免费av在线播放| 国产av又大| 日韩免费av在线播放| 女人精品久久久久毛片| 欧美av亚洲av综合av国产av| 老司机午夜福利在线观看视频| 国产精品香港三级国产av潘金莲| 欧美中文综合在线视频| 91麻豆精品激情在线观看国产 | 欧美另类亚洲清纯唯美| 久久久精品免费免费高清| 久久人妻av系列| 日本黄色日本黄色录像| 精品高清国产在线一区| 香蕉丝袜av| 老司机亚洲免费影院| 亚洲国产欧美一区二区综合| 亚洲欧美激情在线| 欧美色视频一区免费| 好男人电影高清在线观看| 久久久久久亚洲精品国产蜜桃av| 欧美在线黄色| 亚洲欧洲精品一区二区精品久久久| 免费在线观看完整版高清| 脱女人内裤的视频| 在线天堂中文资源库| 免费日韩欧美在线观看| 国产精品美女特级片免费视频播放器 | 咕卡用的链子| 欧美乱妇无乱码| 日韩 欧美 亚洲 中文字幕| 亚洲 欧美一区二区三区| 夜夜爽天天搞| 免费看a级黄色片| 50天的宝宝边吃奶边哭怎么回事| 99国产精品一区二区蜜桃av | 婷婷丁香在线五月| 亚洲精品美女久久av网站| 欧美日韩精品网址| 99久久精品国产亚洲精品| 欧美黑人欧美精品刺激| 曰老女人黄片| 久久这里只有精品19| 精品乱码久久久久久99久播| www日本在线高清视频| 日本一区二区免费在线视频| 久久精品国产清高在天天线| 9191精品国产免费久久| 午夜精品在线福利| 欧美黑人精品巨大| av电影中文网址| 一级作爱视频免费观看| 国产精品国产av在线观看| 又黄又粗又硬又大视频| 亚洲熟女精品中文字幕| 午夜福利乱码中文字幕| 中文字幕高清在线视频| 免费看十八禁软件| 在线观看一区二区三区激情| 黄色a级毛片大全视频| 久久国产精品影院| 久久人妻熟女aⅴ| 又紧又爽又黄一区二区| 一区在线观看完整版| 欧美av亚洲av综合av国产av| 一级毛片精品| 最近最新免费中文字幕在线| 色在线成人网| av国产精品久久久久影院| 中文字幕另类日韩欧美亚洲嫩草| 精品国产超薄肉色丝袜足j| 99久久综合精品五月天人人| 一二三四社区在线视频社区8| 国产精品1区2区在线观看. | 真人做人爱边吃奶动态| 高清毛片免费观看视频网站 | 一本一本久久a久久精品综合妖精| 亚洲一区中文字幕在线| 色婷婷av一区二区三区视频| 欧美日本中文国产一区发布| 在线观看66精品国产|