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

    魚(yú)類群體運(yùn)動(dòng)的元胞自動(dòng)機(jī)模型中的最小勢(shì)能原理

    2021-07-25 17:50:14陸興遠(yuǎn)袁衛(wèi)鋒
    自動(dòng)化學(xué)報(bào) 2021年6期
    關(guān)鍵詞:自動(dòng)機(jī)元胞魚(yú)群

    陸興遠(yuǎn) 袁衛(wèi)鋒

    群體運(yùn)動(dòng)是自然界中一種常見(jiàn)的生物行為.在過(guò)去的幾十年中,研究人員已經(jīng)對(duì)鳥(niǎo)群、獸群、魚(yú)群和細(xì)菌菌落等不同的社會(huì)有機(jī)體的群體現(xiàn)象進(jìn)行了廣泛的研究[1?5].眾多研究表明,某些規(guī)律對(duì)生物的群體行為起到關(guān)鍵作用.其中,“靠近被排斥,遠(yuǎn)離被吸引”的規(guī)則對(duì)群體的分散、聚合和避免個(gè)體間碰撞以及躲避天敵攻擊有重要作用[6?8];流體力學(xué)和動(dòng)力學(xué)原理也影響著群體的運(yùn)動(dòng)狀態(tài).例如,Dai 等研究發(fā)現(xiàn),魚(yú)類群體在運(yùn)動(dòng)中通過(guò)隊(duì)形來(lái)塑造能減小阻力的流場(chǎng),使得個(gè)體耗能最小[9].

    生物的群體運(yùn)動(dòng)是一種個(gè)體間協(xié)同相互作用的結(jié)果,群體中個(gè)體的信息傳遞對(duì)群體運(yùn)動(dòng)狀態(tài)起到重要作用[10].Gueron 和Couzin 等研究發(fā)現(xiàn),群體信息在群體間是通過(guò)逐層傳遞的,個(gè)體的行為決策會(huì)受到相鄰個(gè)體的影響[11?12].群體中個(gè)體間交互協(xié)同的思想為智能機(jī)器人集群的自動(dòng)控制提供了理論基礎(chǔ).例如,Li 等受到生物系統(tǒng)的啟發(fā),發(fā)明了群體顆粒機(jī)器人.雖然每個(gè)個(gè)體的運(yùn)動(dòng)是單一隨機(jī)的運(yùn)動(dòng),但通過(guò)整個(gè)群體的協(xié)同,機(jī)器人群體可以朝著指定的激勵(lì)源方向移動(dòng)[13],雷斌受生物群體行為啟發(fā).應(yīng)用仿生方法研究群體機(jī)器人的編隊(duì)控制,目標(biāo)搜索,實(shí)現(xiàn)群體機(jī)器人間的協(xié)同控制[14].

    在之前的研究中,研究人員大多在二維的假設(shè)條件下對(duì)群體運(yùn)動(dòng)進(jìn)行研究,鮮有在三維的條件下研究群體運(yùn)動(dòng),這與真實(shí)的群體狀態(tài)不符,并且之前的研究者提出的模型大多比較復(fù)雜.本研究構(gòu)建了一個(gè)基于元胞自動(dòng)機(jī)方法和最小勢(shì)能原理的簡(jiǎn)潔三維模型,并通過(guò)數(shù)值模擬研究勢(shì)能在魚(yú)類群體中運(yùn)動(dòng)決策和個(gè)體間相互協(xié)同的作用以及對(duì)魚(yú)群在有限空間中形成旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài)的群體機(jī)制的影響.

    1 模型

    元胞自動(dòng)機(jī)是一種個(gè)體狀態(tài)由其鄰近個(gè)體狀態(tài)決定的算法[15].該方法已經(jīng)被應(yīng)用于流體力學(xué),人流疏散等復(fù)雜場(chǎng)景的研究中[16?17].Pattanayak 等通過(guò)類元胞自動(dòng)機(jī)方法的自航個(gè)體模型(Self-propelled particles)對(duì)群體運(yùn)動(dòng)的狀態(tài)轉(zhuǎn)換進(jìn)行研究[18].本研究假定魚(yú)群的運(yùn)動(dòng)模式受勢(shì)能的支配,相應(yīng)提出了魚(yú)類個(gè)體運(yùn)動(dòng)的勢(shì)函數(shù)并將之融合到元胞自動(dòng)機(jī)中,建立了一種簡(jiǎn)潔的三維模型以模擬魚(yú)群的運(yùn)動(dòng).在真實(shí)魚(yú)類群體中,個(gè)體會(huì)以一定的速度游動(dòng),且速度不是恒定不變的.本研究假設(shè)模擬的虛擬魚(yú)個(gè)體以隨機(jī)的速度運(yùn)動(dòng).虛擬魚(yú)的速度v定義為:v=v+α·Δv,其中α是[?1,1]的隨機(jī)數(shù).在本研究建立的元胞自動(dòng)機(jī)模型中,以半徑為v的球面上的點(diǎn)作為個(gè)體下一時(shí)間步可選擇的位置,如圖1 所示.個(gè)體將會(huì)選擇勢(shì)能最小的位置作為下一時(shí)間步的位置.本模型中每個(gè)位置的勢(shì)能均與個(gè)體的鄰居相關(guān).

    視覺(jué)是動(dòng)物的共性,其對(duì)生物的群體行為和個(gè)體決策有著重要的影響.視覺(jué)影響個(gè)體感知周圍環(huán)境,是個(gè)體獲取信息的重要方式.研究表明,環(huán)境變化、鄰居數(shù)量和鄰居行為等都會(huì)影響個(gè)體行為決策[19].Ballerini 等研究發(fā)現(xiàn)個(gè)體在群體中與其他個(gè)體交流時(shí)并不是與整個(gè)群體交流,而是與其中的6~7 個(gè)鄰居交流[20].本研究假設(shè)個(gè)體的視覺(jué)為有限且不變的常值.如圖1 所示,Ri為個(gè)體的視覺(jué)距離;以θi為圓錐半角的圓錐體區(qū)域?yàn)閭€(gè)體的視角范圍,圓錐體以外的區(qū)域?yàn)閭€(gè)體的盲區(qū).在本模型中,個(gè)體只會(huì)把在其視覺(jué)范圍之內(nèi)距離其最近的N個(gè)鄰居作為參考鄰居,而不是整個(gè)魚(yú)群.個(gè)體在選擇參考鄰居時(shí),如果視覺(jué)范圍內(nèi)個(gè)體總數(shù)大于最大參考鄰居數(shù)N,就會(huì)以距離其最近的N個(gè)鄰居作為參考個(gè)體.

    圖1 魚(yú)的模型Fig.1 Model of fish

    研究人員發(fā)現(xiàn)個(gè)體在群體中會(huì)因?yàn)橄嗷ラg距離太近而相互排斥;如果個(gè)體間的距離太遠(yuǎn),個(gè)體間又會(huì)相互吸引.個(gè)體間存在相對(duì)穩(wěn)定的距離[6?8,21].本模型中,假設(shè)個(gè)體i和參考鄰居j相互之間存在一個(gè)最佳距離R0,在這個(gè)位置上,個(gè)體i的距離勢(shì)能最小.當(dāng)個(gè)體間距離小于R0時(shí),個(gè)體間距離越小,個(gè)體間相互排斥越強(qiáng)烈;相反,當(dāng)個(gè)體間距離大于R0時(shí),隨著距離的增大,個(gè)體間的相互吸引程度增強(qiáng),關(guān)系如式(2)和圖3 (a)所示.

    式(1)中,N為個(gè)體i的參考鄰居數(shù)量.

    式(2)中,H代表(Heviside)階躍函數(shù);uij=dij ?R0,R0為常數(shù).kd1

    式(3)、(4)中,kh、kθ1、kθ2、均為系數(shù),θih表示個(gè)體i下一時(shí)間步的速度vi(T+1)矢量與水平面的夾角;θij表示個(gè)體i下一時(shí)間步的速度vi(T+1)和參考鄰居j當(dāng)前時(shí)間步的速度vj(T)的矢量夾角;θiL表示個(gè)體i下一時(shí)間步的位置B和參考鄰居j當(dāng)前時(shí)間步的位置C的連線與個(gè)體i下一時(shí)間步的速度vi(T+1)的夾角.關(guān)系如圖2 和圖3 (b)所示.圖2 (a)中為參考鄰居j當(dāng)前時(shí)間步的運(yùn)動(dòng)方向vj(T)的平移,與vi(T+1)相交.圖2 (b)中點(diǎn)劃線為vi(T+1)在水平面的投影.

    圖2 三維空間中虛擬魚(yú)中的角度關(guān)系Fig.2 The angular relationship in the potential energy of the fish in three-dimension

    圖3 勢(shì)能函數(shù)的定義Fig.3 Functions defined for potential

    2 元胞自動(dòng)機(jī)規(guī)則

    本文把時(shí)間離散為若干個(gè)時(shí)間步,并把虛擬魚(yú)個(gè)體的下一時(shí)間步可能移動(dòng)的空間范圍劃分為若干個(gè)位置點(diǎn).虛擬魚(yú)群在每個(gè)時(shí)間步中遵循相同的規(guī)則.本文定義元胞自動(dòng)機(jī)的規(guī)則如下面5 個(gè)步驟.

    1)假設(shè)當(dāng)前時(shí)間步為T.虛擬魚(yú)個(gè)體i確定視覺(jué)范圍內(nèi)所有的參考鄰居和下一時(shí)間步所有可能的移動(dòng)位置,如圖1 所示;

    2)計(jì)算下一時(shí)間步T+1 所有可能移動(dòng)位置的勢(shì)能,每個(gè)位置的勢(shì)能關(guān)系均如式(1)~(4)所示;

    3)比較步驟2)中所有位置的勢(shì)能,確定勢(shì)能最小的位置,并把該位置作為下一時(shí)間步T+1的位置;

    4)使用規(guī)則1)~3)遍歷虛擬魚(yú)群中所有的個(gè)體;

    5)更新所有個(gè)體的位置,并重復(fù)步驟1)~5).

    3 實(shí)驗(yàn)

    本研究對(duì)有限空間中紅斑馬魚(yú)群的運(yùn)動(dòng)狀態(tài)進(jìn)行研究,構(gòu)建了一個(gè)長(zhǎng)寬高分別為120 cm,70 cm,80 cm 的玻璃缸和一個(gè)38 cm,25 cm,20 cm 的箱體,并購(gòu)置了360 條紅斑馬魚(yú),為了讓魚(yú)群適應(yīng)實(shí)驗(yàn)室的環(huán)境,在實(shí)驗(yàn)室里飼養(yǎng)超過(guò)一個(gè)月,利用攝像機(jī)對(duì)紅斑馬魚(yú)群的運(yùn)動(dòng)狀態(tài)進(jìn)行攝像記錄.值得注意的是,試驗(yàn)用的大小水箱被放置在同一個(gè)實(shí)驗(yàn)室中,除尺寸外其它條件相同.環(huán)境的變化對(duì)魚(yú)群是一種外來(lái)刺激,當(dāng)把紅斑馬魚(yú)群從大箱轉(zhuǎn)移到小箱中,或者從小箱轉(zhuǎn)移到大箱中時(shí),魚(yú)群都會(huì)出現(xiàn)轉(zhuǎn)圈的現(xiàn)象.如圖4 所示,當(dāng)把魚(yú)群從大箱體轉(zhuǎn)移到小箱體中時(shí),魚(yú)群會(huì)出現(xiàn)有序的旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài).在實(shí)驗(yàn)中可以觀察到,紅斑馬魚(yú)群從無(wú)序的狀態(tài)逐漸轉(zhuǎn)變?yōu)橛行虻臓顟B(tài),最后整個(gè)群體形成有序的旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài).從整個(gè)過(guò)程來(lái)看,魚(yú)群處于無(wú)序的混亂狀態(tài)下,群體中出現(xiàn)個(gè)別個(gè)體發(fā)起有序的運(yùn)動(dòng),然后加入有序運(yùn)動(dòng)的個(gè)體不斷增多,群體運(yùn)動(dòng)狀態(tài)從無(wú)序的混亂狀態(tài)變?yōu)榫植坑行虻膸钸\(yùn)動(dòng)狀態(tài).當(dāng)大部分個(gè)體參與到有序的群體運(yùn)動(dòng)中時(shí),帶狀的群體運(yùn)動(dòng)狀態(tài)首尾相連,整個(gè)魚(yú)群形成環(huán)狀的旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài).有序群體狀態(tài)開(kāi)始時(shí),帶狀的群體形態(tài)的前部個(gè)體數(shù)量較多,較為密集,而群體尾部個(gè)體數(shù)量較少,較為稀疏.當(dāng)最終形成環(huán)狀旋轉(zhuǎn)群體狀態(tài)時(shí),整個(gè)群體中個(gè)體分布較為均勻.

    圖4 真實(shí)紅斑馬魚(yú)群中的旋轉(zhuǎn)群體狀態(tài)Fig.4 Whirling state in real fish schools of red zebrafish

    4 模擬

    本研究對(duì)真實(shí)的紅斑馬魚(yú)群觀察發(fā)現(xiàn),當(dāng)魚(yú)群從大的箱體轉(zhuǎn)移到小的箱體中時(shí),魚(yú)群會(huì)形成旋渦狀的群體形態(tài).為了研究這種群體運(yùn)動(dòng)形成的機(jī)理,本文提出了一種魚(yú)類個(gè)體運(yùn)動(dòng)的勢(shì)函數(shù),并基于最小勢(shì)能原理和元胞自動(dòng)機(jī)方法建立了一個(gè)三維模型,如式(1)~(3)所示,通過(guò)數(shù)值模擬與紅斑馬魚(yú)群的運(yùn)動(dòng)形態(tài)對(duì)比,以驗(yàn)證模型的合理性.

    本模型假設(shè)參考鄰居j分別通過(guò)θiL、θij、uij影響個(gè)體i的位置、方向和相互之間的距離,并假設(shè)θih與個(gè)體i在深度方向上傾斜游動(dòng)的角度相關(guān).當(dāng)相關(guān)系數(shù)處于適當(dāng)?shù)闹禃r(shí),模擬的虛擬魚(yú)群最終會(huì)演化成為旋渦狀的群體形態(tài).模型中的各系數(shù)綜合地影響虛擬魚(yú)群的運(yùn)動(dòng)狀態(tài),只有當(dāng)各系數(shù)處于合適的取值范圍內(nèi)時(shí),模擬結(jié)果才會(huì)與觀察現(xiàn)象較好地吻合.本文中,對(duì)長(zhǎng)寬高分別為30a,30a,20a的方形箱體中的500 條體長(zhǎng)為a的虛擬魚(yú)進(jìn)行數(shù)值模擬.當(dāng)系數(shù)分別設(shè)置為kd1=4.0×105,kd2=2.2×109,kθ1=1.8×109,kθ2=9.0×109,kh=7.0×1010,R0=1.5a,且虛擬魚(yú)群中所有個(gè)體的位置和運(yùn)動(dòng)方向的初始狀態(tài)都隨機(jī)設(shè)置時(shí)(如圖5,第1 步),虛擬魚(yú)群經(jīng)過(guò)若干時(shí)間步的運(yùn)動(dòng)后,演化為旋轉(zhuǎn)的群體形態(tài),如圖5 所示.從整個(gè)演化過(guò)程看,虛擬魚(yú)群開(kāi)始時(shí)處于隨機(jī)無(wú)序狀態(tài),運(yùn)動(dòng)一段時(shí)間后群體中會(huì)形成局部有序的狀態(tài),然后參與到有序運(yùn)動(dòng)的個(gè)體逐漸增多,原來(lái)帶狀的群體形態(tài)逐漸變長(zhǎng),最后帶狀群體形態(tài)首尾相連形成環(huán)形的旋轉(zhuǎn)運(yùn)動(dòng)形態(tài),最終大部分個(gè)體參與到有序的旋轉(zhuǎn)運(yùn)動(dòng)中,整個(gè)虛擬魚(yú)群形成比較穩(wěn)定的旋渦狀群體運(yùn)動(dòng)狀態(tài).值得注意的是,最終的旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài)并不是所有的個(gè)體一直參與其中,而是存在個(gè)別個(gè)體脫離群體,運(yùn)動(dòng)到群體外部,也存在個(gè)別個(gè)體在脫離群體后,重新加入到旋轉(zhuǎn)的群體運(yùn)動(dòng)中.對(duì)比虛擬魚(yú)群與真實(shí)的紅斑馬魚(yú)群的旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài)及其形成過(guò)程,表明了本模型的模擬現(xiàn)象與真實(shí)魚(yú)群現(xiàn)象相符.

    圖5 數(shù)值模擬過(guò)程中旋轉(zhuǎn)的群體狀態(tài)Fig.5 Whirling state in numerical simulation

    群體運(yùn)動(dòng)是個(gè)復(fù)雜的自然現(xiàn)象.仿真雖然可以總體上得到與觀察相吻合的現(xiàn)象,但是在仿真中將個(gè)體的運(yùn)動(dòng)與實(shí)際觀察進(jìn)行定量比較是非常困難的[6],因此,本文通過(guò)對(duì)魚(yú)類群體運(yùn)動(dòng)的定性分析闡明所提出的模型的合理性.必須說(shuō)明的是,勢(shì)能公式中權(quán)重系數(shù)的取值對(duì)魚(yú)群運(yùn)動(dòng)的模擬結(jié)果具有的影響,只有當(dāng)權(quán)重系數(shù)選取合適時(shí),才能達(dá)到模擬結(jié)果和試驗(yàn)觀察吻合的情況.例如,參數(shù)kd2在范圍“8.0×108~3.0×109”區(qū)間時(shí),數(shù)值模擬能夠呈現(xiàn)魚(yú)群轉(zhuǎn)圈的現(xiàn)象,超出這個(gè)范圍,魚(yú)群運(yùn)動(dòng)很難出現(xiàn)穩(wěn)定的旋轉(zhuǎn)狀態(tài).如圖6 所示,在同一個(gè)時(shí)刻(第730 時(shí)間步),當(dāng)kd2處于較小的值時(shí),群體出現(xiàn)過(guò)于密集的狀態(tài);而當(dāng)kd2處于較大的值時(shí),群體狀態(tài)較為混亂無(wú)序;只有當(dāng)kd2處于合適的范圍內(nèi)時(shí),模擬的效果與觀察到的現(xiàn)象比較吻合.另一方面,kd2參數(shù)取值區(qū)間跨越一個(gè)數(shù)量級(jí),說(shuō)明利用本文模型對(duì)魚(yú)群定性研究時(shí),模擬結(jié)果對(duì)參數(shù)的具體數(shù)值不敏感,這種穩(wěn)定性有利于通過(guò)模擬獲取魚(yú)群運(yùn)動(dòng)的共有特性.在本文模型的基礎(chǔ)上,不同的權(quán)重系數(shù)取值會(huì)導(dǎo)致模擬中出現(xiàn)更多的魚(yú)群運(yùn)動(dòng)形態(tài),從而有可能為群體運(yùn)動(dòng)研究提供新的關(guān)注點(diǎn).從這個(gè)意義上講,結(jié)合最小能量原理而建立的元胞自動(dòng)機(jī)模型可以被看作一個(gè)描述魚(yú)類群體運(yùn)動(dòng)的理論框架,而模型中具體參數(shù)的確定,需要更加廣泛和深入的研究.

    圖6 kd2 取值不同時(shí)虛擬魚(yú)群在第730 時(shí)間步的狀態(tài)Fig.6 The status of fish schools at 730th time step subject to different kd2 values

    5 結(jié)論

    本文對(duì)紅斑馬魚(yú)群在有限空間中的運(yùn)動(dòng)形態(tài)進(jìn)行研究,假設(shè)魚(yú)群中每個(gè)個(gè)體在運(yùn)動(dòng)過(guò)程中總是試圖保持自身的勢(shì)能最小,相應(yīng)提出了魚(yú)類個(gè)體運(yùn)動(dòng)的勢(shì)函數(shù)并將之融合到元胞自動(dòng)機(jī)中以模擬魚(yú)群的運(yùn)動(dòng).本研究發(fā)現(xiàn):1)參考鄰居j分別通過(guò)θiL、θij、uij影響個(gè)體i的位置、方向和相互之間的距離;θih與個(gè)體i在深度方向上傾斜游動(dòng)的角度相關(guān).當(dāng)相關(guān)系數(shù)處于適當(dāng)?shù)闹禃r(shí),虛擬魚(yú)群最終可能演化成為旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài).2)參考鄰居數(shù)量N的值與群體中個(gè)體數(shù)量有關(guān),個(gè)體數(shù)量增多,參考鄰居的個(gè)體數(shù)量也會(huì)增多.當(dāng)參考鄰居數(shù)量N在4~10 個(gè)個(gè)體之間時(shí),魚(yú)群均可能出現(xiàn)旋轉(zhuǎn)運(yùn)動(dòng)狀態(tài).本研究實(shí)驗(yàn)觀察到的現(xiàn)象與數(shù)值模擬現(xiàn)象具有較好的一致性,并且本模型具有較好的魯棒性,二者證明了本模型能夠正確反映有限空間中魚(yú)類群體旋轉(zhuǎn)運(yùn)動(dòng)的形成機(jī)理.本模型有望應(yīng)用于智能水下機(jī)器人或空中無(wú)人機(jī)集群控制中.

    本文將能量最小原理和元胞自動(dòng)機(jī)結(jié)合,提出了一種數(shù)值模型,通過(guò)調(diào)試模型中各參數(shù)的數(shù)值,能夠模擬出魚(yú)群運(yùn)動(dòng)的一些規(guī)律,說(shuō)明能量最小原理有可能是魚(yú)類群體運(yùn)動(dòng)需要遵守的一個(gè)基本原理.本研究的意義在于該模型為研究魚(yú)類或其他生物的群體運(yùn)動(dòng)提供了一個(gè)數(shù)值計(jì)算的理論框架,有助于后續(xù)研究在此框架的基礎(chǔ)上進(jìn)一步揭示群體和其所處環(huán)境之間的物理關(guān)系對(duì)群體運(yùn)動(dòng)特性的影響規(guī)律.本文僅針對(duì)包含同種個(gè)體的魚(yú)群的運(yùn)動(dòng)形態(tài)進(jìn)行了研究,而自然界中常有不同魚(yú)類群體共存的現(xiàn)象,對(duì)不同生物群體間相互作用機(jī)理的研究將更具挑戰(zhàn)性和理論價(jià)值.

    猜你喜歡
    自動(dòng)機(jī)元胞魚(yú)群
    {1,3,5}-{1,4,5}問(wèn)題與鄰居自動(dòng)機(jī)
    一種基于模糊細(xì)胞自動(dòng)機(jī)的新型疏散模型
    基于元胞自動(dòng)機(jī)下的交通事故路段仿真
    廣義標(biāo)準(zhǔn)自動(dòng)機(jī)及其商自動(dòng)機(jī)
    魚(yú)群漩渦
    中外文摘(2017年19期)2017-10-10 08:28:41
    基于改進(jìn)魚(yú)群優(yōu)化支持向量機(jī)的短期風(fēng)電功率預(yù)測(cè)
    基于人工魚(yú)群算法的光伏陣列多峰MPPT控制策略
    基于元胞數(shù)據(jù)的多維數(shù)據(jù)傳遞機(jī)制
    多子群并行人工魚(yú)群算法的改進(jìn)研究
    基于AIS的航道移動(dòng)瓶頸元胞自動(dòng)機(jī)模型
    此物有八面人人有两片| 国产精品久久久av美女十八| 亚洲中文字幕一区二区三区有码在线看 | 露出奶头的视频| 欧美黑人巨大hd| 女人被狂操c到高潮| 看片在线看免费视频| 婷婷精品国产亚洲av| 国产av在哪里看| 国产精品九九99| 免费看十八禁软件| 亚洲在线自拍视频| 免费在线观看影片大全网站| 亚洲精品国产精品久久久不卡| 大香蕉久久成人网| 制服人妻中文乱码| 黑人巨大精品欧美一区二区mp4| 美女大奶头视频| 精品久久久久久久毛片微露脸| 国产熟女午夜一区二区三区| 怎么达到女性高潮| 亚洲国产欧美日韩在线播放| 日本精品一区二区三区蜜桃| 久久香蕉激情| 亚洲 欧美一区二区三区| 成人特级黄色片久久久久久久| 国产高清激情床上av| 久久精品国产亚洲av高清一级| 成年免费大片在线观看| 欧美 亚洲 国产 日韩一| 99国产极品粉嫩在线观看| 女生性感内裤真人,穿戴方法视频| 久久九九热精品免费| 级片在线观看| 亚洲人成网站在线播放欧美日韩| 色综合亚洲欧美另类图片| 人成视频在线观看免费观看| 丰满人妻熟妇乱又伦精品不卡| 女性被躁到高潮视频| 色在线成人网| 成人午夜高清在线视频 | 欧美午夜高清在线| 青草久久国产| 成人18禁高潮啪啪吃奶动态图| 久久香蕉精品热| 久久99热这里只有精品18| 国产成人精品久久二区二区91| 成人国语在线视频| 亚洲色图av天堂| 国产又色又爽无遮挡免费看| 色播亚洲综合网| 色综合欧美亚洲国产小说| 一本大道久久a久久精品| 国产男靠女视频免费网站| 露出奶头的视频| 国内精品久久久久久久电影| av电影中文网址| 国产男靠女视频免费网站| 天天添夜夜摸| 精品久久蜜臀av无| 激情在线观看视频在线高清| 国产黄色小视频在线观看| 精品一区二区三区四区五区乱码| 午夜日韩欧美国产| 亚洲第一av免费看| 黄片大片在线免费观看| 淫妇啪啪啪对白视频| 在线观看www视频免费| 波多野结衣高清作品| 给我免费播放毛片高清在线观看| 亚洲人成网站高清观看| 91字幕亚洲| 他把我摸到了高潮在线观看| videosex国产| 欧美日韩亚洲综合一区二区三区_| 19禁男女啪啪无遮挡网站| 一区二区日韩欧美中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产欧美日韩av| 国产精品香港三级国产av潘金莲| 在线永久观看黄色视频| 国产精品久久久久久精品电影 | 美女扒开内裤让男人捅视频| 成人18禁高潮啪啪吃奶动态图| 亚洲午夜理论影院| 老司机在亚洲福利影院| 精品乱码久久久久久99久播| 精品熟女少妇八av免费久了| 色播在线永久视频| 男女做爰动态图高潮gif福利片| 一边摸一边抽搐一进一小说| 黑人欧美特级aaaaaa片| 欧美日韩黄片免| 国产精品九九99| 在线看三级毛片| 不卡av一区二区三区| 中文字幕人成人乱码亚洲影| 大香蕉久久成人网| 色播在线永久视频| 啦啦啦免费观看视频1| 亚洲一码二码三码区别大吗| 久久国产精品人妻蜜桃| 18禁观看日本| 人成视频在线观看免费观看| 国产99久久九九免费精品| 老熟妇仑乱视频hdxx| 精品久久久久久久毛片微露脸| 久久久久久亚洲精品国产蜜桃av| 久久久久国产精品人妻aⅴ院| 午夜福利欧美成人| 国产av又大| 可以在线观看毛片的网站| 大型黄色视频在线免费观看| 免费看美女性在线毛片视频| 欧美大码av| 亚洲精品国产区一区二| 国产精品自产拍在线观看55亚洲| 欧美中文综合在线视频| 夜夜看夜夜爽夜夜摸| 美女高潮喷水抽搐中文字幕| 97超级碰碰碰精品色视频在线观看| 一本综合久久免费| 国产主播在线观看一区二区| 淫秽高清视频在线观看| 女生性感内裤真人,穿戴方法视频| 国产三级在线视频| 国产在线观看jvid| 听说在线观看完整版免费高清| 最新在线观看一区二区三区| 黄片播放在线免费| 国产乱人伦免费视频| 午夜福利18| 国产av不卡久久| 国产精品久久视频播放| 少妇的丰满在线观看| 看片在线看免费视频| 亚洲中文字幕一区二区三区有码在线看 | 日本 欧美在线| 性欧美人与动物交配| 色综合婷婷激情| 99国产精品99久久久久| 国产精品免费视频内射| 成人国产一区最新在线观看| 欧美av亚洲av综合av国产av| 国产精品一区二区三区四区久久 | 久久草成人影院| 精品久久久久久久久久免费视频| 色播在线永久视频| 少妇熟女aⅴ在线视频| 熟女电影av网| 18禁黄网站禁片免费观看直播| 精品乱码久久久久久99久播| 国产伦人伦偷精品视频| 99久久国产精品久久久| 国产精品永久免费网站| 啦啦啦观看免费观看视频高清| 久久性视频一级片| 老司机福利观看| 亚洲午夜精品一区,二区,三区| 精品国产亚洲在线| www日本黄色视频网| 天天添夜夜摸| 18禁黄网站禁片午夜丰满| 99精品在免费线老司机午夜| 亚洲中文字幕一区二区三区有码在线看 | 精品久久久久久久毛片微露脸| 亚洲成人久久爱视频| 日日摸夜夜添夜夜添小说| 在线观看免费视频日本深夜| 国产精品久久久久久亚洲av鲁大| 精品国产一区二区三区四区第35| 在线免费观看的www视频| 亚洲av熟女| 制服诱惑二区| 欧洲精品卡2卡3卡4卡5卡区| 免费在线观看日本一区| 国产精品一区二区精品视频观看| 此物有八面人人有两片| 亚洲熟妇熟女久久| 熟女少妇亚洲综合色aaa.| 国产高清videossex| 国产成人欧美| 男女视频在线观看网站免费 | xxxwww97欧美| 亚洲熟妇中文字幕五十中出| 亚洲av熟女| 在线av久久热| 免费观看精品视频网站| 黑人操中国人逼视频| 欧美国产精品va在线观看不卡| 老熟妇仑乱视频hdxx| 国产一卡二卡三卡精品| 日本a在线网址| 精品一区二区三区四区五区乱码| 亚洲欧洲精品一区二区精品久久久| 天天一区二区日本电影三级| 男人舔女人的私密视频| 丰满的人妻完整版| 黄频高清免费视频| 亚洲人成伊人成综合网2020| 亚洲国产日韩欧美精品在线观看 | 国产在线精品亚洲第一网站| 最近最新中文字幕大全电影3 | 一区福利在线观看| 一进一出好大好爽视频| 国产精品久久视频播放| 精品午夜福利视频在线观看一区| 白带黄色成豆腐渣| 深夜精品福利| 久久精品成人免费网站| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区三| 日韩三级视频一区二区三区| 天天添夜夜摸| 国产欧美日韩一区二区精品| 国产高清videossex| 国产精品国产高清国产av| 人人妻人人澡欧美一区二区| 岛国视频午夜一区免费看| 中国美女看黄片| 国产v大片淫在线免费观看| 国产亚洲精品av在线| 欧美+亚洲+日韩+国产| 在线天堂中文资源库| 亚洲五月婷婷丁香| 国产激情久久老熟女| 久久精品国产亚洲av香蕉五月| 美女高潮喷水抽搐中文字幕| 色av中文字幕| 999精品在线视频| 亚洲精品中文字幕一二三四区| 成人一区二区视频在线观看| 丰满的人妻完整版| 精品日产1卡2卡| 国产成人欧美| 午夜免费激情av| 超碰成人久久| 女生性感内裤真人,穿戴方法视频| 日本三级黄在线观看| 午夜免费成人在线视频| 操出白浆在线播放| 一区二区三区激情视频| 69av精品久久久久久| 国产激情欧美一区二区| 久热这里只有精品99| 无人区码免费观看不卡| 99久久久亚洲精品蜜臀av| 大型av网站在线播放| 一区二区三区国产精品乱码| 亚洲第一青青草原| 99热这里只有精品一区 | 欧美大码av| 亚洲国产精品999在线| 久久香蕉激情| 国产99白浆流出| 日日干狠狠操夜夜爽| 啦啦啦韩国在线观看视频| 国产av一区二区精品久久| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 在线免费观看的www视频| av免费在线观看网站| 天天躁夜夜躁狠狠躁躁| 亚洲aⅴ乱码一区二区在线播放 | 亚洲aⅴ乱码一区二区在线播放 | 亚洲五月天丁香| 国产精品亚洲av一区麻豆| 亚洲国产高清在线一区二区三 | 久久久久久久久中文| 欧美黄色片欧美黄色片| 色尼玛亚洲综合影院| 最好的美女福利视频网| 久99久视频精品免费| 亚洲熟女毛片儿| 免费女性裸体啪啪无遮挡网站| 免费在线观看日本一区| 88av欧美| 中文在线观看免费www的网站 | 啪啪无遮挡十八禁网站| 久久久久国产精品人妻aⅴ院| 淫秽高清视频在线观看| 精品久久久久久久人妻蜜臀av| 精品乱码久久久久久99久播| 欧美中文综合在线视频| 久久国产精品男人的天堂亚洲| 日本一区二区免费在线视频| 黑人巨大精品欧美一区二区mp4| 亚洲欧美一区二区三区黑人| 99久久99久久久精品蜜桃| 哪里可以看免费的av片| 国产精品一区二区免费欧美| 欧美激情 高清一区二区三区| 黄片大片在线免费观看| 婷婷丁香在线五月| 99国产精品一区二区蜜桃av| 日本一区二区免费在线视频| 九色国产91popny在线| 动漫黄色视频在线观看| 国产精品爽爽va在线观看网站 | 精品国内亚洲2022精品成人| 日韩精品青青久久久久久| 18禁美女被吸乳视频| 国产高清videossex| 国内少妇人妻偷人精品xxx网站 | 身体一侧抽搐| 听说在线观看完整版免费高清| 欧美另类亚洲清纯唯美| 丝袜在线中文字幕| 亚洲av第一区精品v没综合| 制服人妻中文乱码| 日韩免费av在线播放| www.www免费av| 免费无遮挡裸体视频| 久久中文字幕一级| 欧美乱码精品一区二区三区| 国产一区二区三区在线臀色熟女| 99国产精品99久久久久| 91成年电影在线观看| 成人一区二区视频在线观看| 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人| 黄色视频不卡| 真人做人爱边吃奶动态| 90打野战视频偷拍视频| 精品国产国语对白av| 国产精品二区激情视频| 91老司机精品| 亚洲久久久国产精品| 精品久久久久久久久久久久久 | aaaaa片日本免费| 成人av一区二区三区在线看| 99在线视频只有这里精品首页| 国产精品亚洲一级av第二区| 视频在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 日韩av在线大香蕉| 欧美不卡视频在线免费观看 | 午夜久久久在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日本一区二区免费在线视频| 曰老女人黄片| 亚洲国产精品sss在线观看| 手机成人av网站| 人人妻人人澡欧美一区二区| 亚洲熟妇熟女久久| 视频在线观看一区二区三区| 我的亚洲天堂| 97人妻精品一区二区三区麻豆 | 狂野欧美激情性xxxx| 久久久久久免费高清国产稀缺| 日韩视频一区二区在线观看| 国产v大片淫在线免费观看| 亚洲三区欧美一区| 免费观看人在逋| 欧美大码av| 亚洲美女黄片视频| 在线观看免费日韩欧美大片| netflix在线观看网站| 亚洲激情在线av| 国内精品久久久久久久电影| 亚洲 欧美 日韩 在线 免费| 欧美日韩乱码在线| 中文资源天堂在线| 搡老妇女老女人老熟妇| www日本在线高清视频| 亚洲欧美一区二区三区黑人| 白带黄色成豆腐渣| a在线观看视频网站| 侵犯人妻中文字幕一二三四区| 99久久无色码亚洲精品果冻| 亚洲欧美精品综合一区二区三区| a级毛片在线看网站| 国产成人精品久久二区二区免费| 一区二区日韩欧美中文字幕| 成人18禁在线播放| 岛国视频午夜一区免费看| 久久久久久久久免费视频了| 亚洲av成人av| 免费一级毛片在线播放高清视频| videosex国产| 中文字幕精品免费在线观看视频| 免费人成视频x8x8入口观看| 在线观看免费视频日本深夜| 男人操女人黄网站| 999久久久国产精品视频| 99久久国产精品久久久| 色在线成人网| www.精华液| 草草在线视频免费看| 法律面前人人平等表现在哪些方面| av欧美777| 欧美日韩亚洲国产一区二区在线观看| 一区二区日韩欧美中文字幕| 久久午夜综合久久蜜桃| 99在线人妻在线中文字幕| 久久青草综合色| 欧美激情极品国产一区二区三区| 亚洲国产日韩欧美精品在线观看 | 日日摸夜夜添夜夜添小说| 琪琪午夜伦伦电影理论片6080| 黄片播放在线免费| 久久久久久久午夜电影| 中文资源天堂在线| 亚洲av熟女| 亚洲精品一卡2卡三卡4卡5卡| 久久精品成人免费网站| 黄色 视频免费看| 99精品久久久久人妻精品| 99久久久亚洲精品蜜臀av| 国产麻豆成人av免费视频| 亚洲专区中文字幕在线| 成人国语在线视频| 国产高清激情床上av| 欧美大码av| 老汉色av国产亚洲站长工具| 久久国产亚洲av麻豆专区| 亚洲无线在线观看| 香蕉久久夜色| 母亲3免费完整高清在线观看| 一二三四在线观看免费中文在| 久久天躁狠狠躁夜夜2o2o| 变态另类成人亚洲欧美熟女| 给我免费播放毛片高清在线观看| www.自偷自拍.com| 十分钟在线观看高清视频www| 99re在线观看精品视频| 美女扒开内裤让男人捅视频| 岛国视频午夜一区免费看| 亚洲精品av麻豆狂野| 久久香蕉国产精品| 精品久久久久久久久久久久久 | 色av中文字幕| 亚洲专区中文字幕在线| or卡值多少钱| 日本三级黄在线观看| 亚洲欧美日韩无卡精品| 男人操女人黄网站| 不卡av一区二区三区| av在线播放免费不卡| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| 久久99热这里只有精品18| 国产成人av教育| 久久婷婷成人综合色麻豆| 久99久视频精品免费| 国产野战对白在线观看| 伊人久久大香线蕉亚洲五| 在线观看免费日韩欧美大片| 亚洲一区中文字幕在线| 欧美精品啪啪一区二区三区| 长腿黑丝高跟| 少妇的丰满在线观看| 99久久久亚洲精品蜜臀av| 久久性视频一级片| 日韩欧美 国产精品| 亚洲人成伊人成综合网2020| 男人舔女人下体高潮全视频| 伦理电影免费视频| 国产亚洲欧美在线一区二区| 在线十欧美十亚洲十日本专区| 亚洲国产欧美一区二区综合| 每晚都被弄得嗷嗷叫到高潮| 99久久精品国产亚洲精品| 在线观看一区二区三区| 两个人看的免费小视频| 琪琪午夜伦伦电影理论片6080| 熟妇人妻久久中文字幕3abv| av超薄肉色丝袜交足视频| 亚洲精品中文字幕在线视频| 日本撒尿小便嘘嘘汇集6| 国产不卡一卡二| 黑人欧美特级aaaaaa片| 亚洲av日韩精品久久久久久密| 亚洲va日本ⅴa欧美va伊人久久| 久久精品成人免费网站| 亚洲全国av大片| 亚洲一区高清亚洲精品| 好男人电影高清在线观看| av超薄肉色丝袜交足视频| 真人一进一出gif抽搐免费| 女性被躁到高潮视频| 亚洲av日韩精品久久久久久密| 亚洲 欧美一区二区三区| 亚洲精品色激情综合| 国产又黄又爽又无遮挡在线| 色综合亚洲欧美另类图片| 国产激情久久老熟女| 在线观看一区二区三区| 女人被狂操c到高潮| 亚洲av成人不卡在线观看播放网| 久久香蕉激情| 一区二区三区激情视频| 亚洲黑人精品在线| 欧美日韩瑟瑟在线播放| 欧美激情高清一区二区三区| 国产精品亚洲美女久久久| 日韩视频一区二区在线观看| 亚洲av片天天在线观看| 亚洲男人天堂网一区| 欧美中文综合在线视频| 丁香欧美五月| 首页视频小说图片口味搜索| 人人妻人人澡欧美一区二区| 日韩一卡2卡3卡4卡2021年| 男人舔女人的私密视频| 国产极品粉嫩免费观看在线| 男人舔女人的私密视频| 无人区码免费观看不卡| 97碰自拍视频| 午夜激情福利司机影院| 欧美日本亚洲视频在线播放| xxxwww97欧美| 精品不卡国产一区二区三区| 国产精品亚洲一级av第二区| 免费女性裸体啪啪无遮挡网站| 男男h啪啪无遮挡| 啦啦啦 在线观看视频| 亚洲天堂国产精品一区在线| 亚洲片人在线观看| 一卡2卡三卡四卡精品乱码亚洲| 久久精品夜夜夜夜夜久久蜜豆 | 视频区欧美日本亚洲| 久久 成人 亚洲| 免费高清在线观看日韩| 亚洲专区国产一区二区| 久久久精品国产亚洲av高清涩受| 午夜福利一区二区在线看| 99精品在免费线老司机午夜| 亚洲国产精品成人综合色| 久久久久久人人人人人| 91麻豆av在线| 亚洲熟妇中文字幕五十中出| 欧美国产精品va在线观看不卡| 欧美乱妇无乱码| 亚洲欧美日韩无卡精品| 两个人看的免费小视频| 中文字幕另类日韩欧美亚洲嫩草| 国产一级毛片七仙女欲春2 | 18禁国产床啪视频网站| 国产伦在线观看视频一区| 在线观看免费日韩欧美大片| 人人妻人人澡欧美一区二区| 久久热在线av| 欧美日韩乱码在线| 亚洲一区二区三区色噜噜| 久久久国产成人精品二区| 亚洲久久久国产精品| 欧美性长视频在线观看| 麻豆成人午夜福利视频| 国产真人三级小视频在线观看| 亚洲人成网站高清观看| 51午夜福利影视在线观看| 麻豆久久精品国产亚洲av| 国产亚洲欧美精品永久| 女性生殖器流出的白浆| 老司机在亚洲福利影院| 国产99久久九九免费精品| 久久久久久久精品吃奶| 亚洲全国av大片| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区三区在线臀色熟女| 欧美成人午夜精品| www国产在线视频色| 亚洲精品美女久久av网站| 成人国产综合亚洲| 非洲黑人性xxxx精品又粗又长| 久久精品国产清高在天天线| 久久精品亚洲精品国产色婷小说| 大香蕉久久成人网| 亚洲自拍偷在线| 亚洲在线自拍视频| 久久久久久久精品吃奶| 自线自在国产av| 91大片在线观看| 搡老妇女老女人老熟妇| 一进一出抽搐gif免费好疼| 国产精品 欧美亚洲| 黄色视频,在线免费观看| 国产午夜福利久久久久久| 人妻久久中文字幕网| 久久久国产精品麻豆| 欧美成人午夜精品| 欧美一级毛片孕妇| 国产一区二区激情短视频| 国产精品99久久99久久久不卡| 久久久精品国产亚洲av高清涩受| 午夜免费观看网址| 国产成+人综合+亚洲专区| 成年人黄色毛片网站| 久久狼人影院| 成人欧美大片| 99riav亚洲国产免费| 亚洲国产精品999在线| 免费一级毛片在线播放高清视频| 这个男人来自地球电影免费观看| 波多野结衣av一区二区av| 日韩大尺度精品在线看网址| a级毛片在线看网站| 婷婷精品国产亚洲av在线| 中文字幕最新亚洲高清| 日韩欧美免费精品| 国产成人欧美| 国产成人啪精品午夜网站| av在线天堂中文字幕| 国产99白浆流出| 精品欧美一区二区三区在线| 哪里可以看免费的av片| 又紧又爽又黄一区二区| 精品欧美一区二区三区在线| 欧美日韩中文字幕国产精品一区二区三区| 又紧又爽又黄一区二区| 精品久久蜜臀av无| 久久久久久大精品| 999久久久精品免费观看国产| 美女扒开内裤让男人捅视频|