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

    基于Excel的遺傳算法在大學(xué)生數(shù)學(xué)建模中的應(yīng)用

    2020-06-05 02:49:42顏清苗壯徐調(diào)斌蔣昌猛賴鑫生
    關(guān)鍵詞:輪盤表格交叉

    顏清,苗壯,徐調(diào)斌,蔣昌猛,賴鑫生

    (上饒師范學(xué)院 數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院,江西 上饒 334001)

    遺傳算法在工程優(yōu)化和工業(yè)控制等領(lǐng)域得到廣泛應(yīng)用。吳立華等針對某機(jī)械的懸臂結(jié)構(gòu)采用遺傳算法進(jìn)行優(yōu)化設(shè)計(jì),有較快的收斂速度和很高的計(jì)算精度,建立了單目標(biāo)非線性優(yōu)化設(shè)計(jì)數(shù)學(xué)模型,為多目標(biāo)非線性結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供了強(qiáng)有力的方法和工具[1]。李月引入免疫算法來改進(jìn)遺傳算法,在TSP應(yīng)用中解決種群迭代次數(shù)過多的復(fù)雜問題中,經(jīng)過疫苗檢驗(yàn)和退火選擇降低算法的退化概率,避免了種群的退化現(xiàn)象[2]。播種機(jī)的播種質(zhì)量、播種效率和能耗由于受播種地形和地域的影響不能發(fā)揮到最佳狀態(tài),秦華等提出了一種基于BTO模式的多目標(biāo)速度控制遺傳算法優(yōu)化模型,確定了BTO模式下播種機(jī)性能優(yōu)化的初始參數(shù),明顯的改善了播種機(jī)的播種質(zhì)量、播種效率和播種能耗,提高了播種機(jī)的作業(yè)性能[3]。針對基本遺傳算法存在收斂速度慢、早熟等不足,王雷等提出了一種改進(jìn)自適應(yīng)遺傳算法,根據(jù)進(jìn)化過程中個體適應(yīng)度值的大小自動調(diào)節(jié)交叉概率和變異概率,使算法能夠跳出局部最優(yōu)解,移動機(jī)器人路徑規(guī)劃的仿真實(shí)驗(yàn)表明,改進(jìn)的遺傳算法有效可行,明顯提高了機(jī)器人路徑規(guī)劃的質(zhì)量[4]。每年的大學(xué)生數(shù)學(xué)建模競賽中遺傳算法也用于解決數(shù)學(xué)建模中的優(yōu)化問題。劉兆仁等利用遺傳算法對公共自行車如何在各大中型城市高效調(diào)度建立了VRP模型,采用標(biāo)準(zhǔn)算例針對算法進(jìn)行檢驗(yàn),結(jié)果表明遺傳算法能夠有效求解自行車的調(diào)度模型[5]。根據(jù)“立竿見影”和竿影日照圖的原理,于鵬等提出了一種太陽影子定位方法。他們首先結(jié)合太陽高度角、太陽赤緯角,以理論影長和實(shí)際影長的相關(guān)系數(shù)最大和其誤差平方和最小為目標(biāo)函數(shù)建立了求太陽影子定位的多目標(biāo)優(yōu)化模型,采用遺傳算法進(jìn)行求解,求解結(jié)果無論在精度還是在收斂速度上都優(yōu)于傳統(tǒng)的枚舉算法[6]。在每年一屆的大學(xué)生數(shù)學(xué)建模教學(xué)及指導(dǎo)中,如何讓大學(xué)生在一段較短時間內(nèi)較快地掌握遺傳算法是一較為棘手的問題。筆者將遺傳算法嵌入在Excel中,利用Excel表格功能、函數(shù)功能及VBA混合編程對遺傳算法的優(yōu)化過程進(jìn)行仿真,模擬遺傳算法的選擇、交叉和變異三個基本操作,充分展現(xiàn)了遺傳算法解決大學(xué)生數(shù)學(xué)建模中優(yōu)化命題的魅力,獲得了較為滿意的應(yīng)用效果,解決了大學(xué)生數(shù)學(xué)建模中的優(yōu)化問題。

    1 遺傳算法的VBA優(yōu)化仿真概述

    遺傳算法即模擬自然界生物群體一代一代的繁衍進(jìn)化,通過迭代操作,產(chǎn)生新一代較優(yōu)的個體作為問題的解,如同生物群體一代一代不斷繁衍、進(jìn)化,最后收斂到最適合生存環(huán)境的一代生物個體,即為問題的最優(yōu)解。問題解的優(yōu)劣通常用適宜度表達(dá),適宜度即表示為生物群體不斷繁衍中個體對生存環(huán)境適宜性的強(qiáng)弱程度,適宜度大的個體被選擇為繁衍下一代父代的概率就大,其染色體就會在一個甚至多個下一代的基因中出現(xiàn),這種有選擇性的將最優(yōu)染色體進(jìn)入下一代的繁衍行為,表現(xiàn)出達(dá)爾文物擇競天、適者生存、優(yōu)勝劣汰的自然界生存法則,這種最優(yōu)選擇策略在遺傳算法中起到?jīng)Q定性作用。

    所謂的遺傳算法VBA優(yōu)化仿真,即在Excel表格中用一組單元格中的{0,1}二進(jìn)制串表示遺傳算法中的一個“染色體”個體,再由Excel表格功能自動仿真計(jì)算“染色體”個體的環(huán)境“適宜度”,由此再通過競爭在表格中產(chǎn)生更適應(yīng)環(huán)境的新一代仿真“染色體”個體群,一個新一代“染色體”個體即為問題的一個可行解。通過Excel表格功能與VBA混合編程,在Excel表格中生動地模擬出“染色體”的競爭過程:“染色體”個體基因如何被選擇、復(fù)制進(jìn)入下一代,如何通過交叉、變異過程產(chǎn)生更適應(yīng)環(huán)境的新一代“染色體”個體群等,如此通過一代一代進(jìn)化,最終收斂到最適應(yīng)環(huán)境的最新一代的“染色體”個體群,從而獲得問題的最優(yōu)解等過程。

    2 算法過程仿真

    以下通過函數(shù)f(x)=11sin(6x)+7cos(5x)為例,在數(shù)據(jù)處理的常用軟件Excel中,通過Excel的表格功能、函數(shù)功能及與自帶的VBA混合編程,即利用Excel遺傳算法求函數(shù)在x∈[0,2π]區(qū)間中最大值,闡述Excel遺傳算法的優(yōu)化仿真過程。

    2.1 Excel遺傳算法VBA優(yōu)化仿真過程

    2.1.1 遺傳算法的初始群體

    求解的精度決定“染色體”個體的{0,1}二進(jìn)制串長度。設(shè)定求解精度x為四位小數(shù),須將0到2π區(qū)間分為6.283 2×104等分,編碼的二進(jìn)制串至少需要16位,即32 768=215<6.283 2×104<216=65 536,即Excel遺傳算法通過16個連續(xù)單元格中的“0”或“1”的二進(jìn)制串作為“染色體”個體的遺傳編碼。20個“染色體”的群即由20行、每行16個“0”或“1”的連續(xù)單元格即可表示。利用Excel的自動計(jì)算功能,在“染色體”個體的{0,1}二進(jìn)制串之后的單元格,可以分別對“染色體”個體的“十進(jìn)制值”“實(shí)際x值”“函數(shù)f(x)值”進(jìn)行仿真計(jì)算[7]。

    2.1.2 初始群體適宜值的計(jì)算與遺傳選擇

    由于f(x)在x∈[0,2π]區(qū)間上有正有負(fù),以20個“染色體”群中最大的正數(shù)與最小的負(fù)數(shù)之間的區(qū)間為1,“染色體”群中最大的正數(shù)其“區(qū)間值”為1,最小的負(fù)數(shù)其“區(qū)間值”為0。根據(jù)適者生存的生物遺傳法則,以函數(shù)值f(x)的“區(qū)間值”作為適宜值,將“染色體”群中“染色體”個體的“區(qū)間值”累加,個體的“區(qū)間值”占“區(qū)間值”總和的百分比值就是“染色體”遺傳進(jìn)入下一代的概率。根據(jù)個體的“區(qū)間值”的百分比值設(shè)計(jì)成輪盤賭(圖1),每個個體的這個百分比值即為“輪盤比例”,在輪盤賭上根據(jù)“染色體”順序劃定確定的區(qū)域,適宜值高的“染色體”個體在輪盤賭上占有的區(qū)域就大,被選擇的機(jī)會就高,充分體現(xiàn)了適者生存的生物遺傳法則。

    圖1即為初始群的Excel遺傳算法仿真計(jì)算表。圖中輪盤賭的仿真數(shù)值區(qū)域在0-1之間,刻度值與時鐘的0點(diǎn)至12點(diǎn)分度對應(yīng),通過隨機(jī)函數(shù)每次選擇一個0-1之間的隨機(jī)數(shù),確定選擇輪盤賭對應(yīng)區(qū)間的“染色體”,可以看出適宜度高的“染色體”由于其“輪盤比例”大而被多次選擇,一些適宜值低的“染色體”在輪盤賭中對應(yīng)區(qū)間小而難以選擇而被淘汰。適宜度最低的11號“染色體”的“輪盤比例”值為0,直接淘汰。圖1所示,由于函數(shù)f(x)的值為“適宜值”,16號“染色體”的“適宜度”最大,“輪盤比例”達(dá)9.63%,選擇次數(shù)也最多,20次的選擇中,被選擇了3次,即16號“染色體”的基因復(fù)制了3次進(jìn)入下一代“染色體”中,占比達(dá)15%;17號“染色體”的“適宜度”較大,僅次于16號“染色體”,“輪盤比例”為9.27%,也被選擇了3次。第1、3、4、5、7、9、11、18、20號“染色體”,平均“適宜度”占比小于2.61%,20次的選擇中,選擇次數(shù)均為零。注意到,圖1中最后兩列的“選擇值”數(shù)據(jù)列僅對應(yīng)于前列的選擇序號n列,其它各列對應(yīng)行數(shù)據(jù)為各“染色體”個體對應(yīng)數(shù)據(jù)。

    圖1 遺傳算法的選擇與輪盤賭

    圖2 遺傳算法的復(fù)制與仿真計(jì)算

    2.2 遺傳算法VBA優(yōu)化仿真混合編程

    2.2.1 選擇、復(fù)制確定下一代“染色體”基因

    圖1、圖2是依次執(zhí)行“初始群”“選擇”和“復(fù)制”三按鈕中的VBA代碼的仿真計(jì)算結(jié)果。經(jīng)輪盤賭選擇后,通過Excel混合編程進(jìn)行復(fù)制。遺傳復(fù)制遵循“優(yōu)勝劣汰”的自然遺傳法則,可以看出,許多“染色體”在遺傳復(fù)制時未被選擇,淘汰率達(dá)45%。圖1與圖2的仿真計(jì)算顯示,個體的“適宜值”大入選概率就大,“適宜值”小的如第3、5、7、9、11、18號“染色體”,“適宜值”占比小于1.35%均被淘汰。如下為VBA混合編程完成復(fù)制操作的程序代碼。

    L=0'將選擇的“染色體”{0,1}二進(jìn)制串由“Excel遺傳算法”表復(fù)制到“Excel遺傳算法1”表

    For i=1 To 20

    m=Sheets(“Excel遺傳算法”).Cells(i + 1, 60)

    If m=0 Then GoTo 10

    For n=1 To m

    For j=1 To 16

    Sheets(“Excel遺傳算法1”).Cells(L + 2, j + 1) = Sheets(“Excel遺傳算法”).Cells(i+1, j+1)

    Next j

    L=L+1

    Next n

    10 Next i

    被選擇的“染色體”由“Excel遺傳算法1”表轉(zhuǎn)移回“Excel遺傳算法”表后完成選擇、復(fù)制

    For i=1 To 20

    For j=1 To 16

    Sheets(“Excel遺傳算法”).Cells(i+1, j+1)=Sheets(“Excel遺傳算法1”).Cells(i+1, j+1)

    Next j

    Next i

    從VBA程序代碼中看出,輪盤賭選擇確定了新一代“染色體”的基因,根據(jù)“初始染色體群”表(圖1)中的BH列列出對應(yīng)序號的仿真“染色體”個體的選擇次數(shù)進(jìn)行復(fù)制,得到第二代“染色體”群的基因(圖2),其平均“適宜值”比“初始染色體群”要高出許多,整個算法的選擇、復(fù)制仿真過程躍然紙上。

    2.2.2 交叉、變異產(chǎn)生新一代“染色體”個體

    圖3為第七代“染色體”群經(jīng)選擇、復(fù)制后的仿真圖,圖4為經(jīng)配對完成交叉操作后的“染色體”群的仿真圖。對照圖3圖4可以看出,交叉配對產(chǎn)生了更為優(yōu)秀的個體,交叉配對前“染色體”個體“適宜值”最高為17.582 882,交叉配對后“染色體”個體“適宜值”最高為17.595 591,未選取交叉配對的3對即1、11,7、17,8、18等三對個體經(jīng)變異后競爭進(jìn)入下一代“染色體”群。更為可喜的是交叉配對前最優(yōu)秀的五個“染色體”個體,其中4個“染色體”經(jīng)交叉操作“適宜值”得到進(jìn)一步提高,未變的一個“染色體”原因是未交叉直接進(jìn)入,由此可見,交叉配對產(chǎn)生的新個體使優(yōu)秀“染色體”進(jìn)一步優(yōu)化。

    圖3 新一代個體的選擇與復(fù)制

    圖4 交叉、變異產(chǎn)生新一代個體仿真圖

    變異概率確定為0.01,即100個長度的{0,1}二進(jìn)制串中,隨機(jī)抽取一個進(jìn)行變異,群體的{0,1}二進(jìn)制串累計(jì)為320位,共有3.2個二進(jìn)制位變異,即3個或4個二進(jìn)制位變異。變異是一個較小的擾動,只對“染色體”群帶來較小的改變。

    歷經(jīng)選擇、交叉、變異等三個算子作用后,下一代仿真“染色體”群宣告誕生。結(jié)合表格計(jì)算功能的VBA混合編程可以形象地讓學(xué)生了解遺傳算法的基本思想。

    2.3 遺傳算法VBA優(yōu)化仿真的終止條件

    遺傳算法的終止條件通常設(shè)為兩類:(1)最大遺傳代數(shù)的設(shè)定;(2)種群中基因多樣性測度。兩個終止條件通常會結(jié)合起來運(yùn)用,若將最大遺傳代數(shù)的設(shè)定為100代,在不到100代時,計(jì)算種群中基因多樣性測度,若所有基因位相似程度急驟上升,此時收斂于概率為1的優(yōu)良個體,遺傳算法即可終止。

    3 遺傳算法VBA優(yōu)化仿真的應(yīng)用

    3.1 遺傳算法VBA優(yōu)化仿真輔助數(shù)學(xué)建模教學(xué)

    大學(xué)生數(shù)學(xué)建模與優(yōu)化問題聯(lián)系非常密切,遺傳算法解決大學(xué)生數(shù)學(xué)建模優(yōu)化問題是較重要途徑。大學(xué)生學(xué)習(xí)遺傳算法VBA優(yōu)化仿真,很容易掌握遺傳算法的基本思想,熟悉選擇、交叉、變異等基本算子的作用,在輔助大學(xué)生數(shù)學(xué)建模及相關(guān)教學(xué)中得到應(yīng)用。

    3.1.1 算法過程的仿真及其可視化

    在遺傳算法VBA優(yōu)化仿真過程中,以Excel表格數(shù)據(jù)及圖形的形式揭示“染色體”在遺傳過程中的演變規(guī)律,利用Excel混合編程完成遺傳算法基本算子的操作,實(shí)現(xiàn)遺傳算法演化過程的可視化。

    3.1.2 “染色體”個體的仿真

    在Excel表格中,利用一組16個連續(xù)單元格表示的{0,1}二進(jìn)制串,作為仿真“染色體”個體的遺傳算法編碼。如{0,1}二進(jìn)制串“1001000101011001”,十進(jìn)制數(shù)為37 209,若用以表示x為0至6.283 2中的數(shù),即x=0+37 209×6.283 2/(216-1)=3.567 4,f(x)的絕對值為9.798 081 7,都可以通過Excel表格直接計(jì)算得到這些結(jié)果,使算法過程直覺化、可視化,“染色體”個體仿真加深了大學(xué)生算法學(xué)習(xí)的感性認(rèn)識。

    3.1.3 選擇算子中輪盤賭的仿真

    根據(jù)各“染色體”適宜值設(shè)計(jì)輪盤賭,配上插圖就很直觀。圖1中看出整個輪盤賭的數(shù)值區(qū)域在0-1之間,通過隨機(jī)函數(shù)每次選擇一個0-1之間的隨機(jī)數(shù),就可以確定每次是哪個“染色體”被選擇,也可以看出適宜值高的“染色體”往往被多次選擇并遺傳至下一代,有些適宜值低的“染色體”很難被選擇而淘汰。

    3.1.4 交叉操作的可視化

    選擇操作只是完成了產(chǎn)生新一代“染色體”群的第一步,接著就是交叉、變異操作。大學(xué)生可以直接從仿真表格中(圖1)看出,前面最優(yōu)秀的5個“染色體”個體,其中4個“染色體”經(jīng)交叉“適宜值”得到進(jìn)一步提高。大學(xué)生通過對照交叉前后的“染色體”,變化一目了然。由于交叉概率選擇0.7,即只有70%的“染色體”進(jìn)行交叉配對,20個“染色體”個體中有14個“染色體”個體要進(jìn)行交叉配對,最優(yōu)秀的5個“染色體”個體有一個“染色體”未交叉配對直接進(jìn)入新一代“染色體”群。由此可見,交叉配對產(chǎn)生的新個體使優(yōu)秀“染色體”進(jìn)一步得到優(yōu)化。圖5、圖6是1至60代中“染色體”群中最優(yōu)“染色體”與“染色體”平均適宜度的仿真計(jì)算結(jié)果,60代時,最優(yōu)解為17.833 39,與答案17.833 99相差極小。

    圖5 1至60代中的最優(yōu)“染色體” 圖6 1至60代“染色體”平均適宜度

    3.1.5 變異操作方法的多樣性

    變異操作是一個較小的擾動,對“染色體”群帶來較小的改變。通過變異操作,打開了新的搜索空間,容易找出“染色體”個體{0,1}二進(jìn)制串中已經(jīng)變化的幾個位置。歷經(jīng)選擇、交叉、變異等三個算子作用后,下一代“染色體”群宣告誕生,很明顯,大學(xué)生對這一過程心領(lǐng)神會,遺傳算法的VBA優(yōu)化仿真在課堂上以過程仿真的形式講授提升了教學(xué)效果。

    3.2 遺傳算法VBA優(yōu)化仿真與數(shù)學(xué)建模

    課堂上遺傳算法的VBA優(yōu)化仿真形象、生動,學(xué)生容易理解,但要掌握遺傳算法的VBA優(yōu)化仿真在大學(xué)生數(shù)學(xué)建模中解決優(yōu)化命題還需要學(xué)生實(shí)踐操作,最為實(shí)用例子是選擇以往大學(xué)生數(shù)學(xué)建模中的優(yōu)化選題。

    3.2.1 應(yīng)用于適應(yīng)度的變化

    對于函數(shù)f(x)= 11sin(6x)+ 7cos(5x),可以在x∈[0,2π]區(qū)間上對函數(shù)求最大值,也可以對函數(shù)求最小值,這就只要在適應(yīng)度列中作適當(dāng)改變就行。如求最小值,即在適應(yīng)度列的每行加負(fù)號即可。若函數(shù)f(x)改變了,且變量只有一個,也只要改變“適應(yīng)度”列(對應(yīng)“輪盤比例”列)的各行。

    3.2.2 應(yīng)用于求解精度與變量個數(shù)的變化

    遺傳算法的VBA優(yōu)化仿真程序在“染色體”個體與“染色體”群所對應(yīng)的列(圖1)中預(yù)留了32列。如要改變求解精度,可以對“染色體”個體與“染色體”群所對應(yīng)的列進(jìn)行增減。f(x) 函數(shù)關(guān)系改變,“染色體”個體與“染色體”群所對應(yīng)的列有48列,適應(yīng)有3個以內(nèi)變量且求解精度為10-4的f(x)函數(shù)。

    3.2.3 應(yīng)用時克服輪盤賭的缺點(diǎn)

    輪盤賭還與“賭”字有關(guān),如圖1中所示,1、4、20號個體的“適宜度”達(dá)到5.00%以上未被選擇,而“適宜度”占比僅4.19%的15號個體被選擇2次,這也是它的缺點(diǎn)所在。因此,為避免最優(yōu)基因的流失,通常加入最優(yōu)基因保留技術(shù)。

    3.2.4 避免早熟與局部尋優(yōu)問題的思考

    遺傳算法VBA優(yōu)化仿真也可一代一代操作,也可連續(xù)操作直接得到結(jié)果。輔助教學(xué)時,往往采用一代一代操作,學(xué)生能清楚地了解算法的過程細(xì)節(jié)。初始群局限了最優(yōu)解的搜索空間[8],一代一代操作容易發(fā)現(xiàn)遺傳算法的早熟與陷入局部尋優(yōu)等問題,通過對操作的調(diào)整,如選擇較好的初始群就可避免早熟與局部尋優(yōu)。變異操作確保了新一代基因的多樣性,打開了新的搜索空間,可以看出擴(kuò)大搜索空間影響收斂速度,故變異概率常取1%或以下。

    4 結(jié)束語

    在Excel表格中形象地完成遺傳算法的選擇、交叉、變異的過程仿真,將Excel表格自動計(jì)算功能與VBA編程相結(jié)合,算法過程形象地呈現(xiàn)在學(xué)生面前。在Excel表格中操作遺傳算法,Excel算法過程的可視化與實(shí)時動態(tài)仿真效果,使大學(xué)生對算法的原理、過程、算子、適宜度等基本要點(diǎn)了解得非常透徹,使大學(xué)生在數(shù)學(xué)建模中應(yīng)用遺傳算法的VBA優(yōu)化仿真得心應(yīng)手,揮灑自如,為遺傳算法在運(yùn)用時進(jìn)行適當(dāng)改進(jìn)打開了思路[9],為解決數(shù)學(xué)建模中優(yōu)化命題提供了嶄新的解決方案。

    猜你喜歡
    輪盤表格交叉
    《現(xiàn)代臨床醫(yī)學(xué)》來稿表格要求
    某型航空發(fā)動機(jī)鈦合金輪盤模擬疲勞試驗(yàn)件設(shè)計(jì)
    統(tǒng)計(jì)表格的要求
    統(tǒng)計(jì)表格的要求
    統(tǒng)計(jì)表格的要求
    “六法”巧解分式方程
    基于ANSYS的輪盤轉(zhuǎn)子模態(tài)影響因素分析
    連一連
    基于Fast-ICA的Wigner-Ville分布交叉項(xiàng)消除方法
    雙線性時頻分布交叉項(xiàng)提取及損傷識別應(yīng)用
    国产熟女午夜一区二区三区| 久久久精品国产亚洲av高清涩受| 久久久久久九九精品二区国产 | 精品乱码久久久久久99久播| 日韩欧美国产一区二区入口| 动漫黄色视频在线观看| 国产精品一区二区三区四区免费观看 | 亚洲专区中文字幕在线| 欧美激情久久久久久爽电影| 国产精品久久久久久久电影 | 国产爱豆传媒在线观看 | 国产伦在线观看视频一区| 婷婷精品国产亚洲av| 97人妻精品一区二区三区麻豆| 视频区欧美日本亚洲| 国产av又大| 国产视频内射| 免费在线观看视频国产中文字幕亚洲| 亚洲国产看品久久| 亚洲五月天丁香| 欧美成人免费av一区二区三区| 国产午夜精品论理片| 好男人电影高清在线观看| 哪里可以看免费的av片| x7x7x7水蜜桃| 欧美性长视频在线观看| 国产高清激情床上av| 哪里可以看免费的av片| 精品一区二区三区四区五区乱码| 香蕉久久夜色| 五月伊人婷婷丁香| 精品久久久久久久末码| www.精华液| 欧美 亚洲 国产 日韩一| 亚洲成av人片在线播放无| 精品久久久久久久久久久久久| 真人一进一出gif抽搐免费| cao死你这个sao货| 真人一进一出gif抽搐免费| 国产97色在线日韩免费| 久久久久性生活片| 免费观看人在逋| 男女之事视频高清在线观看| or卡值多少钱| 国产三级在线视频| 在线观看免费视频日本深夜| 国产av麻豆久久久久久久| 国产aⅴ精品一区二区三区波| 日日摸夜夜添夜夜添小说| 在线国产一区二区在线| 国产黄片美女视频| 亚洲专区中文字幕在线| 搡老岳熟女国产| 亚洲av电影在线进入| 亚洲成人久久爱视频| 精品国产乱子伦一区二区三区| 国产99白浆流出| 在线免费观看的www视频| 日韩大码丰满熟妇| 脱女人内裤的视频| 国产av在哪里看| 国产一区二区三区在线臀色熟女| 久久九九热精品免费| 国产成人欧美在线观看| 亚洲国产日韩欧美精品在线观看 | 国产成人啪精品午夜网站| 亚洲午夜精品一区,二区,三区| av在线播放免费不卡| 国产成人精品久久二区二区91| 国产日本99.免费观看| 精品一区二区三区四区五区乱码| 亚洲国产精品sss在线观看| 日韩精品免费视频一区二区三区| 亚洲男人天堂网一区| 国产亚洲欧美在线一区二区| 午夜亚洲福利在线播放| 99精品久久久久人妻精品| 女人被狂操c到高潮| 亚洲成a人片在线一区二区| 窝窝影院91人妻| 不卡一级毛片| 美女午夜性视频免费| 男女床上黄色一级片免费看| 色精品久久人妻99蜜桃| 老司机午夜十八禁免费视频| 后天国语完整版免费观看| 狂野欧美白嫩少妇大欣赏| 桃色一区二区三区在线观看| 高清在线国产一区| 午夜日韩欧美国产| 亚洲专区国产一区二区| 国产精品久久久久久亚洲av鲁大| 99久久精品热视频| 成在线人永久免费视频| 中亚洲国语对白在线视频| 久久香蕉激情| 午夜成年电影在线免费观看| 国产成人精品久久二区二区免费| 国产精品99久久99久久久不卡| 大型黄色视频在线免费观看| 婷婷六月久久综合丁香| 精品久久久久久久毛片微露脸| 床上黄色一级片| 亚洲一区二区三区色噜噜| 国产精品久久久久久亚洲av鲁大| 国产高清激情床上av| 久久中文字幕一级| 99国产极品粉嫩在线观看| 久久精品人妻少妇| 久久99热这里只有精品18| 99国产综合亚洲精品| 午夜福利视频1000在线观看| 别揉我奶头~嗯~啊~动态视频| bbb黄色大片| 婷婷丁香在线五月| 怎么达到女性高潮| 国产精品一区二区精品视频观看| 伊人久久大香线蕉亚洲五| 亚洲成人国产一区在线观看| 日本免费a在线| 麻豆成人午夜福利视频| 亚洲avbb在线观看| 免费一级毛片在线播放高清视频| 国产精品98久久久久久宅男小说| 人人妻人人澡欧美一区二区| 九色成人免费人妻av| 亚洲黑人精品在线| 国产日本99.免费观看| 免费观看精品视频网站| 久久久久国产精品人妻aⅴ院| 12—13女人毛片做爰片一| 欧美黑人精品巨大| 国产精品精品国产色婷婷| 久久久国产成人免费| 精品久久久久久成人av| 日韩欧美在线乱码| 在线观看66精品国产| 精品电影一区二区在线| 亚洲精品久久成人aⅴ小说| 国产精品日韩av在线免费观看| 一级片免费观看大全| 免费在线观看成人毛片| 一本精品99久久精品77| 欧美日韩一级在线毛片| 免费在线观看亚洲国产| 国产人伦9x9x在线观看| 亚洲国产日韩欧美精品在线观看 | 亚洲av美国av| 日本免费a在线| 久久香蕉激情| 观看免费一级毛片| 老司机深夜福利视频在线观看| 久久婷婷成人综合色麻豆| av天堂在线播放| 在线观看舔阴道视频| 国产欧美日韩一区二区精品| 日韩欧美免费精品| 桃红色精品国产亚洲av| 久久久国产精品麻豆| 啦啦啦免费观看视频1| 国产av不卡久久| 亚洲五月天丁香| 亚洲免费av在线视频| 日韩欧美国产在线观看| 国内精品一区二区在线观看| 12—13女人毛片做爰片一| av欧美777| or卡值多少钱| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色丝袜av网址大全| 亚洲专区字幕在线| 国产一区二区三区视频了| 国内揄拍国产精品人妻在线| 久久中文字幕一级| 午夜精品在线福利| 日本在线视频免费播放| 日本熟妇午夜| 好男人电影高清在线观看| 在线永久观看黄色视频| 成人午夜高清在线视频| 99久久国产精品久久久| 久久人妻av系列| 日韩欧美国产在线观看| 亚洲精品中文字幕一二三四区| 精品欧美国产一区二区三| 亚洲精品久久成人aⅴ小说| 男女之事视频高清在线观看| 桃色一区二区三区在线观看| 韩国av一区二区三区四区| 好男人电影高清在线观看| 久久久久久久午夜电影| 在线观看日韩欧美| 丁香六月欧美| 欧美成人午夜精品| 亚洲成人久久性| 99国产极品粉嫩在线观看| 女生性感内裤真人,穿戴方法视频| 最近在线观看免费完整版| 欧美大码av| 久久香蕉精品热| 国产91精品成人一区二区三区| 这个男人来自地球电影免费观看| www.精华液| 国产精品免费视频内射| 国产蜜桃级精品一区二区三区| av超薄肉色丝袜交足视频| 成人18禁在线播放| 欧美精品亚洲一区二区| 久久精品aⅴ一区二区三区四区| 欧美zozozo另类| 精品午夜福利视频在线观看一区| 亚洲欧美日韩东京热| 国产精品香港三级国产av潘金莲| 欧美成人性av电影在线观看| 亚洲欧美精品综合久久99| 在线十欧美十亚洲十日本专区| 久久亚洲精品不卡| 亚洲国产精品成人综合色| 色尼玛亚洲综合影院| www日本在线高清视频| 波多野结衣高清无吗| 麻豆国产97在线/欧美 | 免费看十八禁软件| 欧美丝袜亚洲另类 | 两人在一起打扑克的视频| 久久国产精品人妻蜜桃| 国产精品乱码一区二三区的特点| 丁香欧美五月| 精品免费久久久久久久清纯| 国产成人欧美在线观看| 欧美性长视频在线观看| 国产不卡一卡二| 一级片免费观看大全| 高清在线国产一区| 亚洲中文字幕日韩| 国产91精品成人一区二区三区| 黄片大片在线免费观看| 国产黄a三级三级三级人| 午夜免费激情av| 久久热在线av| 亚洲一区二区三区色噜噜| 久久精品影院6| 一本一本综合久久| 国产精品电影一区二区三区| 日本免费一区二区三区高清不卡| АⅤ资源中文在线天堂| 免费在线观看黄色视频的| 日日摸夜夜添夜夜添小说| 嫩草影视91久久| a级毛片a级免费在线| 久久久精品国产亚洲av高清涩受| 91在线观看av| 在线视频色国产色| 欧美zozozo另类| 五月伊人婷婷丁香| 我要搜黄色片| 成熟少妇高潮喷水视频| 国产黄片美女视频| 亚洲精品色激情综合| 亚洲欧美日韩高清专用| 成人一区二区视频在线观看| 中文字幕人成人乱码亚洲影| 日本一二三区视频观看| 亚洲午夜理论影院| 美女大奶头视频| 欧美在线一区亚洲| 欧美中文综合在线视频| 在线a可以看的网站| 欧美黄色片欧美黄色片| 亚洲色图 男人天堂 中文字幕| 一级毛片精品| 午夜福利在线在线| 国产亚洲精品久久久久5区| 99久久综合精品五月天人人| 国内精品一区二区在线观看| 在线永久观看黄色视频| tocl精华| 午夜日韩欧美国产| 啦啦啦免费观看视频1| 午夜福利在线观看吧| 亚洲av成人精品一区久久| 午夜视频精品福利| 国产精品av久久久久免费| 午夜成年电影在线免费观看| 黄色a级毛片大全视频| 麻豆av在线久日| 精品久久久久久久末码| 国产一区二区激情短视频| 两性夫妻黄色片| 久久草成人影院| 1024视频免费在线观看| 国产片内射在线| 色av中文字幕| 淫妇啪啪啪对白视频| 麻豆久久精品国产亚洲av| 俄罗斯特黄特色一大片| 丰满人妻熟妇乱又伦精品不卡| 久久婷婷成人综合色麻豆| 亚洲九九香蕉| 午夜成年电影在线免费观看| 国产午夜福利久久久久久| 色综合站精品国产| 欧美一级毛片孕妇| 精品国产超薄肉色丝袜足j| 亚洲av电影在线进入| 精品日产1卡2卡| 亚洲中文字幕一区二区三区有码在线看 | 亚洲专区国产一区二区| 国产黄a三级三级三级人| 99久久综合精品五月天人人| 亚洲精品在线美女| 欧美日韩亚洲综合一区二区三区_| √禁漫天堂资源中文www| 免费在线观看影片大全网站| 欧美日韩亚洲综合一区二区三区_| 国产精品美女特级片免费视频播放器 | 99在线视频只有这里精品首页| 成人精品一区二区免费| 一进一出抽搐gif免费好疼| 精品久久久久久久毛片微露脸| e午夜精品久久久久久久| 国产免费男女视频| 他把我摸到了高潮在线观看| 精品熟女少妇八av免费久了| 一本大道久久a久久精品| 免费电影在线观看免费观看| 中文字幕精品亚洲无线码一区| 精品福利观看| 天堂动漫精品| 一级毛片女人18水好多| 亚洲精品美女久久久久99蜜臀| 午夜久久久久精精品| 久久精品夜夜夜夜夜久久蜜豆 | 大型黄色视频在线免费观看| 国产成人av教育| 久久久久久人人人人人| 可以在线观看的亚洲视频| 亚洲成人久久性| 后天国语完整版免费观看| 精品国产超薄肉色丝袜足j| 高清毛片免费观看视频网站| 国产精品久久久久久久电影 | 欧美在线一区亚洲| 国产精品98久久久久久宅男小说| 18美女黄网站色大片免费观看| 激情在线观看视频在线高清| 夜夜爽天天搞| 女同久久另类99精品国产91| 青草久久国产| 老汉色∧v一级毛片| 黄片大片在线免费观看| 色哟哟哟哟哟哟| 后天国语完整版免费观看| 国产99白浆流出| 成人手机av| 在线观看www视频免费| 国产成人影院久久av| 国产精品1区2区在线观看.| 男女之事视频高清在线观看| 国产精品国产高清国产av| 一个人免费在线观看电影 | 18禁黄网站禁片免费观看直播| 国产亚洲欧美98| 亚洲精品一区av在线观看| x7x7x7水蜜桃| 欧美av亚洲av综合av国产av| 欧美成人免费av一区二区三区| 一级a爱片免费观看的视频| 久久久久国产一级毛片高清牌| 久久国产乱子伦精品免费另类| 最新在线观看一区二区三区| 18禁黄网站禁片免费观看直播| 欧美一级毛片孕妇| 两个人免费观看高清视频| 欧美中文综合在线视频| 国产私拍福利视频在线观看| 欧美中文综合在线视频| 久久久久亚洲av毛片大全| 日本 欧美在线| 国产亚洲精品综合一区在线观看 | 在线播放国产精品三级| 神马国产精品三级电影在线观看 | 久久国产精品影院| 中文字幕熟女人妻在线| 久久热在线av| 国产三级在线视频| 又爽又黄无遮挡网站| 少妇人妻一区二区三区视频| 熟女少妇亚洲综合色aaa.| av欧美777| www.熟女人妻精品国产| 美女黄网站色视频| 亚洲 欧美一区二区三区| 女生性感内裤真人,穿戴方法视频| 波多野结衣巨乳人妻| 久久久精品欧美日韩精品| 欧美日韩国产亚洲二区| 日韩av在线大香蕉| 一边摸一边做爽爽视频免费| 欧美人与性动交α欧美精品济南到| 999久久久精品免费观看国产| 国产精华一区二区三区| 久久久久久免费高清国产稀缺| 禁无遮挡网站| ponron亚洲| 90打野战视频偷拍视频| 亚洲一区二区三区不卡视频| 人成视频在线观看免费观看| 国产精品亚洲一级av第二区| 99在线视频只有这里精品首页| 日本免费一区二区三区高清不卡| 精品一区二区三区四区五区乱码| 国产成人系列免费观看| 丰满人妻熟妇乱又伦精品不卡| 美女黄网站色视频| 麻豆一二三区av精品| 亚洲av电影在线进入| 国产熟女xx| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 国模一区二区三区四区视频 | 日本熟妇午夜| 搞女人的毛片| 亚洲美女视频黄频| 91国产中文字幕| 日本一本二区三区精品| 久久精品国产亚洲av香蕉五月| 99久久综合精品五月天人人| 伊人久久大香线蕉亚洲五| 又大又爽又粗| 韩国av一区二区三区四区| 色综合亚洲欧美另类图片| 精品熟女少妇八av免费久了| av福利片在线观看| а√天堂www在线а√下载| 波多野结衣巨乳人妻| 亚洲av美国av| 久久久久精品国产欧美久久久| 成人三级黄色视频| 亚洲第一电影网av| 19禁男女啪啪无遮挡网站| 日韩欧美在线乱码| 日韩av在线大香蕉| 欧美日韩黄片免| 热99re8久久精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 久久伊人香网站| 亚洲国产欧洲综合997久久,| 母亲3免费完整高清在线观看| 又紧又爽又黄一区二区| 国产av又大| 亚洲av成人不卡在线观看播放网| 99久久综合精品五月天人人| 黄片小视频在线播放| 天堂av国产一区二区熟女人妻 | 熟女少妇亚洲综合色aaa.| 亚洲成人国产一区在线观看| 真人做人爱边吃奶动态| 亚洲av第一区精品v没综合| 三级男女做爰猛烈吃奶摸视频| 午夜精品在线福利| 99国产极品粉嫩在线观看| 成熟少妇高潮喷水视频| 丰满人妻一区二区三区视频av | 99久久精品国产亚洲精品| 午夜福利欧美成人| 一本精品99久久精品77| 黄色视频不卡| 国产精品乱码一区二三区的特点| 桃红色精品国产亚洲av| 日韩精品青青久久久久久| 亚洲欧美精品综合久久99| 18禁国产床啪视频网站| 黄色毛片三级朝国网站| 制服丝袜大香蕉在线| av福利片在线观看| 久99久视频精品免费| 天天躁狠狠躁夜夜躁狠狠躁| 制服人妻中文乱码| 色综合站精品国产| 欧美激情久久久久久爽电影| tocl精华| 亚洲国产精品sss在线观看| 正在播放国产对白刺激| 亚洲av成人av| 特级一级黄色大片| 一本久久中文字幕| 伊人久久大香线蕉亚洲五| 成人三级做爰电影| 国产成年人精品一区二区| 欧美高清成人免费视频www| 免费人成视频x8x8入口观看| 久久婷婷成人综合色麻豆| 欧美黄色淫秽网站| 国产精品,欧美在线| 麻豆一二三区av精品| 99国产精品一区二区蜜桃av| 动漫黄色视频在线观看| 桃色一区二区三区在线观看| 国产精品一及| 色综合婷婷激情| 久久亚洲真实| 老汉色∧v一级毛片| 国产精品亚洲av一区麻豆| 国产免费av片在线观看野外av| 国产99久久九九免费精品| 久久人妻福利社区极品人妻图片| 在线a可以看的网站| 无人区码免费观看不卡| av欧美777| 日日夜夜操网爽| 国产伦人伦偷精品视频| 妹子高潮喷水视频| 欧美中文日本在线观看视频| 丰满人妻一区二区三区视频av | 成年女人毛片免费观看观看9| 亚洲,欧美精品.| 久久久久免费精品人妻一区二区| 1024视频免费在线观看| 久久久久亚洲av毛片大全| 成人18禁在线播放| 一个人免费在线观看的高清视频| 国产乱人伦免费视频| 操出白浆在线播放| 久久亚洲精品不卡| 一区二区三区国产精品乱码| 香蕉丝袜av| 亚洲欧美日韩东京热| 一级片免费观看大全| 亚洲专区中文字幕在线| 老司机福利观看| 两个人免费观看高清视频| 午夜久久久久精精品| 国产精品日韩av在线免费观看| 亚洲av五月六月丁香网| 午夜精品久久久久久毛片777| 国产一区二区在线观看日韩 | 精品一区二区三区四区五区乱码| 丰满的人妻完整版| 亚洲国产日韩欧美精品在线观看 | 夜夜爽天天搞| 亚洲av成人不卡在线观看播放网| 国产高清videossex| 最近视频中文字幕2019在线8| 麻豆成人av在线观看| 国产视频内射| 亚洲男人天堂网一区| 搞女人的毛片| 一级作爱视频免费观看| 国产精品香港三级国产av潘金莲| 精品久久久久久久人妻蜜臀av| av超薄肉色丝袜交足视频| 波多野结衣高清无吗| 天堂影院成人在线观看| 免费观看人在逋| 老汉色∧v一级毛片| 欧美另类亚洲清纯唯美| 亚洲无线在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲免费av在线视频| 99精品欧美一区二区三区四区| 一本大道久久a久久精品| 精品欧美国产一区二区三| 午夜福利欧美成人| 亚洲av中文字字幕乱码综合| 久久精品成人免费网站| 久久久国产欧美日韩av| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品久久二区二区91| 天天躁夜夜躁狠狠躁躁| 成人18禁在线播放| 丝袜美腿诱惑在线| 三级国产精品欧美在线观看 | 三级男女做爰猛烈吃奶摸视频| 97超级碰碰碰精品色视频在线观看| 亚洲,欧美精品.| 国产精品乱码一区二三区的特点| 日韩有码中文字幕| 久久婷婷人人爽人人干人人爱| 99久久99久久久精品蜜桃| 男人的好看免费观看在线视频 | 精品一区二区三区av网在线观看| 亚洲av熟女| 久久精品人妻少妇| 久久人妻av系列| 欧美成人性av电影在线观看| 伊人久久大香线蕉亚洲五| 一本精品99久久精品77| 亚洲精品中文字幕一二三四区| 999久久久国产精品视频| 久久久久精品国产欧美久久久| 免费观看人在逋| ponron亚洲| 欧美大码av| 岛国在线免费视频观看| 波多野结衣高清无吗| 99热这里只有是精品50| 免费搜索国产男女视频| 亚洲精品一区av在线观看| 啦啦啦观看免费观看视频高清| 日本 欧美在线| 精品久久久久久,| 国产av麻豆久久久久久久| 99riav亚洲国产免费| 天天躁狠狠躁夜夜躁狠狠躁| 国产av麻豆久久久久久久| 99riav亚洲国产免费| 欧美在线黄色| 色综合婷婷激情| 国产av又大| 国产精品久久久人人做人人爽| 欧美绝顶高潮抽搐喷水| 女警被强在线播放| 亚洲av美国av| 国产精品一区二区精品视频观看|