日本在线www-日本在线播放一区-日本在线不卡免费视频一区-日本在线不卡视频-成人影院久久久久久影院-成人影院一区二区三区

最新廣告
關(guān)注中國(guó)自動(dòng)化產(chǎn)業(yè)發(fā)展的先行者!
工業(yè)智能邊緣計(jì)算2025年會(huì)
CAIAC 2025
2025工業(yè)安全大會(huì)
OICT公益講堂
當(dāng)前位置:首頁 >> 案例 >> 案例首頁

案例頻道

人工免疫聚類算法在交通時(shí)段自動(dòng)劃分上的應(yīng)用
  • 企業(yè):控制網(wǎng)     領(lǐng)域:儀器儀表     行業(yè):廣播電視    
  • 點(diǎn)擊數(shù):4071     發(fā)布時(shí)間:2008-03-02 20:21:48
  • 分享到:



    劉強(qiáng)(1979-)
男,遼寧工業(yè)大學(xué)研究生,(遼寧工業(yè)大學(xué)信息科學(xué)與工程學(xué)院,遼寧  錦州  121001 ),研究方向?yàn)橹悄芸刂啤?BR>
摘要:
論文基于免疫網(wǎng)絡(luò)理論和克隆選擇原理,建立了一種人工免疫數(shù)據(jù)聚類分析算法,并詳細(xì)闡述了聚類算法在城市交通時(shí)段自動(dòng)劃分中的具體應(yīng)用。實(shí)例分析表明,該算法可以有效減少聚類數(shù)據(jù)的冗余信息,特別適合于解決分級(jí)聚類等傳統(tǒng)方法不適應(yīng)的大數(shù)據(jù)量聚類問題,對(duì)解決城市交通時(shí)段的自動(dòng)劃分等數(shù)據(jù)聚類問題是可行的和有效的。

關(guān)鍵詞:聚類;人工免疫系統(tǒng);模式識(shí)別;交通控制

Abstract: An artificial immune algorithm for data clustering is proposed in this paper. This work is mainly based on the immune network theory and the clone selection principle extracted from vertebrate immune systems. Experiments show that the given clustering algorithm can effectively reduce the redundant information of data, and thus, especially fit to the mass data clustering applications where hierarchical clustering algorithm or other traditional clustering algorithms may be ineffective.

Key words: clustering;artificial immune system;pattern recognition;traffic signal control

1 引言

    目前城市交通信號(hào)控制應(yīng)用較為普遍的是多時(shí)段定時(shí)控制方案。多時(shí)段定時(shí)控制(time of day,TOD)是根據(jù)車流量等交通信息把一天劃分為若干個(gè)時(shí)間段,不同的時(shí)段采用不同的信號(hào)優(yōu)化控制方案。TOD控制的前提是交通時(shí)段的合理劃分。傳統(tǒng)上,交通工程技術(shù)人員根據(jù)收集到的交通流量,繪制流量時(shí)間曲線圖,根據(jù)曲線的特征人工劃分交通時(shí)段。傳統(tǒng)方法具有很大的主觀性,極易產(chǎn)生不合理的時(shí)段劃分,難以滿足現(xiàn)代城市交通模式快速變化的需求。

    以生物免疫系統(tǒng)為基礎(chǔ)的人工免疫系統(tǒng)(artificial immune systems, AIS) [1]已成為目前國(guó)內(nèi)外計(jì)算智能領(lǐng)域一個(gè)新的研究熱點(diǎn)。本文在詳細(xì)研究生物免疫學(xué)和AIS的最新研究成果的基礎(chǔ)上,建立了一種人工免疫數(shù)據(jù)聚類分析算法,該算法可以有效減少聚類數(shù)據(jù)的冗余信息,特別適合于解決分級(jí)聚類等傳統(tǒng)聚類方法不適應(yīng)的大數(shù)據(jù)量聚類問題。論文把建立的人工免疫算法用于城市交通時(shí)段的自動(dòng)劃分[2],克服了人工時(shí)段劃分的不合理性,為TOD交通時(shí)段的劃分及城市交通控制方案的制定提供了一條新的思路[3]。

2 免疫系統(tǒng)的基本理論

    生物免疫系統(tǒng)是一個(gè)由細(xì)胞、分子和器官組成的復(fù)雜系統(tǒng),其主要功能是限制異物對(duì)肌體的侵害[4]。免疫系統(tǒng)抵御外部入侵,使其機(jī)體免受病原侵害的應(yīng)答反應(yīng)稱為免疫(immunity)。外部有害病原入侵機(jī)體并激活免疫細(xì)胞,誘導(dǎo)其發(fā)生反應(yīng)的過程稱為免疫應(yīng)答。誘導(dǎo)免疫系統(tǒng)產(chǎn)生免疫應(yīng)答的物質(zhì)稱為抗原(antigen,簡(jiǎn)稱Ag),能與抗原進(jìn)行特異性結(jié)合的免疫細(xì)胞稱為抗體(antibody, 簡(jiǎn)稱Ab)。抗原與抗體的結(jié)合強(qiáng)度用親和度(affinity)度量。AIS主要涉及T細(xì)胞和B細(xì)胞的相關(guān)免疫特性。B細(xì)胞的作用是識(shí)別抗原和分泌抗體,T細(xì)胞能夠促進(jìn)和抑制B細(xì)胞的產(chǎn)生與分化。免疫系統(tǒng)的結(jié)構(gòu)示意如圖1所示。生物免疫學(xué)的研究成果是AIS的理論基礎(chǔ)。本文給出的數(shù)據(jù)聚類算法,主要基于生物免疫隱喻的克隆選擇理論[5]和免疫網(wǎng)絡(luò)理論[6]。當(dāng)淋巴細(xì)胞實(shí)現(xiàn)對(duì)抗原的識(shí)別,即抗體與抗原的親和度(affinity)超過一定閾值后,B細(xì)胞被激活并增殖復(fù)制產(chǎn)生細(xì)胞克隆。由于遺傳和免疫細(xì)胞在增殖中的基因突變,形成了免疫細(xì)胞的多樣性。當(dāng)每一種抗原侵入機(jī)體后都能在機(jī)體內(nèi)選擇出能識(shí)別和消滅相應(yīng)抗原的免疫細(xì)胞克隆,使之激活、分化和增殖,進(jìn)行免疫應(yīng)答以最終清除抗原。克隆選擇的主要特征是免疫細(xì)胞在抗原刺激下產(chǎn)生克隆擴(kuò)增,隨后通過遺傳變異分化為多樣性效應(yīng)細(xì)胞和記憶細(xì)胞。



圖1   免疫系統(tǒng)的基本結(jié)構(gòu)

3 人工免疫聚類分析算法

    AIS借鑒了一些生物免疫系統(tǒng)的功能、原理和模型,使用學(xué)習(xí)、記憶和關(guān)聯(lián)檢索來完成識(shí)別與分類任務(wù),具有強(qiáng)大的和魯棒的信息處理能力,是一種有效的數(shù)據(jù)聚類分析方法[7]。AIS將待分析的數(shù)據(jù)作為抗原,通過抗體對(duì)抗原的學(xué)習(xí)、記憶過程來識(shí)別、分類不同數(shù)據(jù),并從中找出數(shù)據(jù)間的相關(guān)性能。

    假設(shè)標(biāo)稱為[0,1]下的Q個(gè)輸入數(shù)據(jù)集構(gòu)成抗原Ag,聚類分析的目的就是要通過免疫學(xué)習(xí),尋找記憶數(shù)據(jù)集M作為輸入數(shù)據(jù)的聚類。

    為描述算法的方便,定義如下變量:

    Ag:輸入抗原,即待聚類的輸入數(shù)據(jù)集合,或研究問題的論域: 

          (1)

    Ab:抗體,即N個(gè)初始網(wǎng)絡(luò)細(xì)胞:

                (2)

    M:記憶數(shù)據(jù)集,即nt個(gè)網(wǎng)絡(luò)記憶細(xì)胞:

                (3)

    :所有抗體與抗原Agj的距離度量矢量和親和力矢量;

    :M中抗體的相似度矢量;

    Nc:激發(fā)細(xì)胞產(chǎn)生的克隆細(xì)胞數(shù)量;

    Cj:激發(fā)細(xì)胞產(chǎn)生的克隆細(xì)胞群體,親和力成熟后變?yōu)镃。

    基于免疫網(wǎng)絡(luò)理論和克隆選擇原理的人工免疫聚類算法由如下步驟構(gòu)成[5,6,7,8]:

    Step1 隨機(jī)初始化N個(gè)抗體Ab,作為初始網(wǎng)絡(luò)細(xì)胞。

    Step2 設(shè)定循環(huán)控制參數(shù)。

    Step3 輸入待聚類數(shù)據(jù)作為抗原,對(duì)抗原Agj(j=1,2,…,Q)進(jìn)行如下操作:

    Step3.1 計(jì)算距離度量矢量Dj和親和力矢量AFj:

    抗原與抗體之間的距離用歐幾里德距離定義,即

       (4)

    根據(jù)聚類分析的特點(diǎn),抗原與抗體之間的親和力采用式(5)定義:

         (5)

    于是

           (6)

        (7)

    Step3.2按照opR(optimal selecting Rate)的比率從Ab中選擇n個(gè)與Agj具有最高親和力的抗體進(jìn)行克隆擴(kuò)增,產(chǎn)生對(duì)應(yīng)的克隆細(xì)胞集合Cj。根據(jù)克隆選擇原理,抗體與抗原的親和力越大,抗體產(chǎn)生的克隆細(xì)胞數(shù)量越多。對(duì)n個(gè)選擇細(xì)胞按親和力從小到大的順序排序,克隆細(xì)胞的數(shù)量按下式計(jì)算。

        (8)

    式中:Nc為n個(gè)抗體產(chǎn)生的克隆細(xì)胞總數(shù);α為乘數(shù)因子,用于控制克隆群體的規(guī)模;Int(.)為取整函數(shù)。

    Step3.3 應(yīng)用式(9)對(duì)克隆的抗體進(jìn)行隨機(jī)變異操作,實(shí)現(xiàn)親和力的成熟,產(chǎn)生具有更高親和力的抗體細(xì)胞C。

         (9)

    式中:rand(CR,NR)為隨機(jī)函數(shù),表示從CR中隨機(jī)抽取NR個(gè)變量;μ為抗體的變異率。生物免疫系統(tǒng)抗體變異的實(shí)質(zhì)是基因片斷的重組,使得子代抗體與抗原的親和力得以迅速提高,而且與抗原具有較大親和力的抗體具有較小的變異率。抗體的變異率按式(10)確定:

        (10)

    式中:AFv為抗體親和力的標(biāo)稱值,k為比例因子,η為衰減控制系數(shù)。k與η的合理取值可以保證變異個(gè)體在其允許的域值范圍[0,1]之內(nèi)。

    Step3.4 計(jì)算C的抗體細(xì)胞與抗原Agj的親和力矢量AFc j。

    Step3.5 按照rsR(re-selecting Rate)的比率,選擇C中若干個(gè)與抗原Agj具有最高親和力的抗體,作為部分記憶細(xì)胞Mp 。

    Step3.6 去除Mp中相似度sij小于閾值σd的抗體,產(chǎn)生新的記憶集Mk ,實(shí)現(xiàn)免疫系統(tǒng)克隆抑制的效果。抗體的相似度sij用抗體間的歐幾里德距離描述,即:

        (11)

    Step3.7 把部分記憶細(xì)胞Mk合并到記憶細(xì)胞集合M(M [M;Mk])。

    Step4 計(jì)算M中各記憶細(xì)胞的相似度矢量S,去除M中相似度sij低于閾值σs的抗體,實(shí)現(xiàn)不同克隆集的網(wǎng)絡(luò)抑制。sij按式(11)計(jì)算。

    Step5 按照wsR(worst selecting Rate)的比率,隨機(jī)產(chǎn)生若干個(gè)抗體替換原抗體中親和力較低的個(gè)體,體現(xiàn)免疫系統(tǒng)的自組織功能。

    Step6 變量替換,返回Step3,進(jìn)行下一代的網(wǎng)絡(luò)學(xué)習(xí),直到達(dá)到要求的學(xué)習(xí)代數(shù),或滿足設(shè)定的聚類迭代要求(如迭代達(dá)到了預(yù)定的記憶細(xì)胞的個(gè)數(shù)或抗原與記憶細(xì)胞的平均親和力達(dá)到了預(yù)定的誤差范圍等)。

    學(xué)習(xí)過程結(jié)束后,算法的輸出M為聚類數(shù)據(jù)的記憶數(shù)據(jù)集,S為記憶數(shù)據(jù)集各抗體之間的相似度。根據(jù)待聚類數(shù)據(jù)與各記憶細(xì)胞的距離,或抗原與記憶細(xì)胞的親和度,可方便的實(shí)現(xiàn)對(duì)抗原數(shù)據(jù)的聚類。上述聚類算法的核心是保持網(wǎng)絡(luò)抗體代克隆、變異及抑制的操作,最終產(chǎn)生記憶數(shù)據(jù)集。

    調(diào)節(jié)抑制閾值σs可以控制網(wǎng)絡(luò)的學(xué)習(xí)規(guī)模和性能,閾值越大,得到的免疫記憶數(shù)據(jù)越多,分類密度越大,數(shù)據(jù)的相似性越高。算法中其他參數(shù)的選擇原則如下:初始化抗體個(gè)數(shù)N根據(jù)經(jīng)驗(yàn)選定,其大小不影響聚類結(jié)果,但選擇過大影響計(jì)算速度,過小則體現(xiàn)不了抗原特征,一般取N=20。克隆乘數(shù)因子α的選擇取決于初始抗體的個(gè)數(shù)N。最佳抗體選擇率opR一般取初始抗體的10%到20%,但opR的取值應(yīng)保證每代至少有一個(gè)最優(yōu)抗體克隆擴(kuò)增。最差抗體選擇率wsR體現(xiàn)免疫系統(tǒng)的自組織特性,一般取值不應(yīng)超過10%。從抗體多樣性的角度考慮,再次選擇率rsR可以取較大的值,但聚類的目的是減少數(shù)據(jù)冗余信息,rsR不易太大,一般取克隆細(xì)胞的10%作為部分記憶細(xì)胞。

4 交通時(shí)段的免疫聚類研究

    交通時(shí)段的人工免疫聚類分析分兩步進(jìn)行:第一步將交通數(shù)據(jù)看作免疫系統(tǒng)的抗原;第二步采用常規(guī)的聚類分析方法,對(duì)交通流量的免疫記憶集做聚類分析,實(shí)現(xiàn)交通時(shí)段的自動(dòng)劃分。

    交通工程實(shí)踐中,一般以15min的交通流量為基礎(chǔ)進(jìn)行交通時(shí)段的劃分,將一天的96個(gè)交通流量數(shù)據(jù)作為免疫系統(tǒng)的抗原。

    4.1 交通數(shù)據(jù)的預(yù)處理

    考慮交通數(shù)據(jù)變化較大(如白天峰值時(shí)段的交通流量與夜晚谷值時(shí)段的交通流量相差十分懸殊),聚類前首先應(yīng)用式(12)對(duì)交通數(shù)據(jù)做標(biāo)準(zhǔn)化變換處理:

        (12)

    式中“*”表示交通流量的原始數(shù)據(jù)。聚類實(shí)踐表明,如果直接應(yīng)用每天的交通數(shù)據(jù)進(jìn)行聚類,結(jié)果存在文獻(xiàn)[2]指出的“時(shí)段沖突”或“聚類重疊”,即同一時(shí)間段在不同的工作日(如每周的星期一)歸屬于不同的類。(3節(jié)中各式的p=1,Ag,Ab,M等矩陣均降為同維列向量)。

    4.2 交通時(shí)段聚類分析實(shí)例

    本文結(jié)合現(xiàn)實(shí)應(yīng)用的需要,利用錦州交警支隊(duì)提供的解放路的交通數(shù)據(jù)庫,進(jìn)行基于人工免疫算法的交通聚類研究,實(shí)現(xiàn)城市交通時(shí)段的自動(dòng)劃分。圖2為根據(jù)歷史交通數(shù)據(jù)的96個(gè)平均流量繪制的解放路與士英街路口星期一的交通流量時(shí)間曲線圖。



圖2   24小時(shí)平均交通流量曲線

    把上述96個(gè)平均交通流量作為抗原應(yīng)用于3節(jié)的人工免疫聚類算法,進(jìn)行免疫識(shí)別和記憶。算法中的參數(shù)設(shè)置如下:初始化抗體個(gè)數(shù)N=20;最佳抗體選擇率opR=20%;克隆乘數(shù)因子α =1;再次選擇率rsR=20%,wsR=10%;最大迭代次數(shù)NA=15;為保證抗體的多樣性,σd取0.001,σs為0.05。表1是交通流數(shù)據(jù)的人工免疫計(jì)算過程,共分15代。

    表1   聚類數(shù)據(jù)的人工免疫計(jì)算過程:

    從數(shù)據(jù)分類結(jié)果可以看出第9次迭代結(jié)果,記憶數(shù)據(jù)集同抗原數(shù)據(jù)之間具有較大的平均親和力,免疫記憶數(shù)據(jù)集的規(guī)模適中,是一種較為理想的數(shù)據(jù)分類結(jié)構(gòu)。

    為直觀起見,將由人工免疫算法得到的記憶數(shù)據(jù)應(yīng)用譜系圖的方法描述如圖3所示,顯然,當(dāng)距離值取0.1時(shí),第9次迭代結(jié)果的14個(gè)記憶數(shù)據(jù)被清晰地分成了5類。



圖3   交通流記憶細(xì)胞集的聚類譜系圖

    采用重心法計(jì)算各交通數(shù)據(jù)與記憶類重心之間的距離,進(jìn)行交通數(shù)據(jù)的二次聚類,并與其交通時(shí)段對(duì)應(yīng),得到解放路與士英街路口星期一的交通時(shí)段劃分情況如圖4所示。即根據(jù)交通流量可以把全天分為如下幾個(gè)交通時(shí)段:類I時(shí)段22:45-5:15;類II時(shí)段12:15-14:00,20:30-22:45;類III時(shí)段5:15-6:30, 8:15-11:00;類Ⅳ時(shí)段11:00-12:15,16:00- 18:00,19:15-20:30;類V時(shí)段6:30-8:15,14:00-16:00,18:00-19:15。



圖4   TOD交通時(shí)段劃分

5 結(jié)束語

    人工免疫數(shù)據(jù)聚類算法可以有效減少聚類數(shù)據(jù)的冗余信息,特別適合于解決分級(jí)聚類等傳統(tǒng)方法不適應(yīng)的大數(shù)據(jù)量聚類問題。作者把建立的人工免疫聚類算法應(yīng)用于城市交通時(shí)段的自動(dòng)劃分,克服了人工時(shí)段劃分的不合理性。合理設(shè)置參數(shù)的取值范圍可以滿足不同聚類問題的需要。

其它作者:

    王艷秋(1955-),女,教授,遼寧工業(yè)大學(xué)信息科學(xué)與工程學(xué)院院長(zhǎng),博士,碩士研究生導(dǎo)師,研究方向:近代交流調(diào)速系統(tǒng);智能控制。
張健(1981- ),男,遼寧工業(yè)大學(xué)研究生,研究方向:智能控制。

參考文獻(xiàn)

    [1]肖人彬, 王磊.人工免疫系統(tǒng):原理、模型、分析及展望[J]. 計(jì)算機(jī)學(xué)報(bào), 2002, 25 (12) : 1282 - 1291. 

    [2]HAUSER T, SCHERER W. Data mining tools for real time traffic signal decision support and maintenance[C]//Proc of IEEE Int Conf on Systems, Man, and Cybernetics. Tucson, USA: IEEE Press, 2001:1471-1477.

    [3]PARK.B, LEE Do-Hoon, YUN Hsoo. Enhancement of Time of Day Based Traffic Signal Control [C]// Proc of IEEE Int Conf on Systems, Man, and Cybernetics. USA: IEEE Press, 2003:3619 - 3624.

    [4]林學(xué)顏, 張玲. 現(xiàn)代細(xì)胞與分子免疫學(xué)[M].北京:科學(xué)出版社, 1999:46-62. 

    [5]CASTROL, ZUBENF. An evolutionary immune network for data clustering[C]// Proc of the Sixth Brazilian Symposium on Neural Networks. Janeiro: IEEE Press, 2000: 84-89.

    [6]DASGU PTAD, FORRESTS. Artificial immune systems in industrial applications[C]//Proc of Second In t Conf on Intelligent Processing and Manufacturing of Materials. Honolulu, USA: IEEE Press, 1999: 257-267.

    [7]Dasnupta D. Artificial Immune Systems and Their Applications[M].New York: Sprinner~Verlan, 1999.

    [8]TIMM IS J, N EAL M. Artificial immune systems for data analysis [J]. B io System s, 2000, 55 (2): 143 - 150.

    [9]D. K. Merchant, G. L. Nemhauser. Optimality conditions for a dynamic traffic assignment model [J]. Transportation Science B, 1978, 12: 200-207.

熱點(diǎn)新聞

推薦產(chǎn)品

x
  • 在線反饋
1.我有以下需求:



2.詳細(xì)的需求:
姓名:
單位:
電話:
郵件:
主站蜘蛛池模板: 久久精品久噜噜噜久久-久久精品久久久久-久久精品久久久-久久精品久久精品久久精品-亚洲影视一区-亚洲影视精品 | 欧美精品一区二区精品久久-欧美精品一二区-2019酷路泽大改款-欧美精品亚洲一区二区在线播放-99日韩-99人中文字幕亚洲区 | 国内自拍中文字幕,久久久一本精品99久久精品66,精品400部自拍视频在线播放,国产麻豆精品在线,日韩欧美高清视频,久久久免费精品视频 | 欧美日韩 国产区 在线观看-欧美日操-欧美日本综合一区二区三区-欧美日本中文字幕-欧美日本中文-欧美日本在线一区二区三区 | 亚洲幼女网-亚洲呦女专区-亚洲永久在线观看-亚洲永久免费视频-精品久久香蕉国产线看观看麻豆-精品久久日日躁夜夜躁AV | 国产a一级-国产a网-国产a视频精品免费观看-国产a视频-国产a免费观看-国产a毛片 | 波多野结衣中文字幕教师-波多野结衣中文字幕2022免费-波多野结衣中文在线观看-波多野结衣中文在线播放-波多野结衣中文一区-波多野结衣中文丝袜字幕 | 热久久国产-热99这里只有精品-热99这里有精品综合久久-热99在线观看-国产精品99久久久-国产精品99久久99久久久看片 | 天堂网在线最新版www-天堂网在线资源-天堂网在线看-天堂网在线观看视频-日本精品无码特级毛片-日本精品无码一区二区三区久久久 | 黄色在线免费播放-黄色在线免费-黄色在线看网站-黄色在线观看网站-色综合小说天天综合网-色综合小说久久综合图片 | 久久久久久99精品-久久久久久99-久久久久久91香蕉国产-久久久久久91精品色婷婷-中国一级片免费看-中国一级毛片国产高清 | 莜田优在线观看-尤物最新网址-尤物综合-尤物自拍-久久久久久网址-久久久久久网站 | 国产精品亚洲精品日韩己满十八小-国产精品亚洲精品观看不卡-国产精品亚洲精品不卡-国产精品亚洲精品-四虎影视精品-四虎影视国产精品永久在线 | 亚洲欧美激情另类,国产成人一区二区三区免费观看,一区二区三区在线视频观看,亚洲一区二区三区精品视频,国产乱了真实在线观看,国产播放器一区 | 精品在线视频播放-精品在线免费播放-精品在线观看一区-精品在线观看国产-亚洲 男人 天堂-亚州视频一区二区 | 久久免费视频1-久久免费看视频-久久免费精品一区二区-久久免费精品视频-中文字幕亚洲专区-中文字幕亚洲一区二区三区 | 欧美日韩 国产区 在线观看-欧美日操-欧美日本综合一区二区三区-欧美日本中文字幕-欧美日本中文-欧美日本在线一区二区三区 | 久久综合影视-久久综合伊人77777麻豆-久久综合一区二区三区-久久综合一-亚欧三级-亚欧美综合 | 91香蕉导航-91香蕉成人免费高清网站-91香蕉成人-91午夜视频-91午夜精品亚洲一区二区三区-91网址在线观看 | 久久精品片-久久精品欧美一区二区-久久精品女人毛片国产-久久精品嫩草影院免费看-在线日韩国产-在线日韩不卡 | 亚洲欧美激情另类,国产成人一区二区三区免费观看,一区二区三区在线视频观看,亚洲一区二区三区精品视频,国产乱了真实在线观看,国产播放器一区 | 国产日韩精品欧美一区-国产日韩高清一区二区三区-国产日韩不卡免费精品视频-国产日产欧美精品一区二区三区-午夜国产精品免费观看-午夜国产精品理论片久久影院 | 久久精品一区二区三区不卡牛牛,国产片精品电影www,久久久久久97,www奇米,好姑娘中文在线播放,美女被羞羞视频网站在线 | 热血影院-热热99-热久在线-热久久最新网址-成人毛片一区二区三区-成人毛片视频免费网站观看 | 里番acg纲手的熟蜜姬训练场-里番※琉璃全彩acg奈亚子-李亚男三级-李丽莎与土豪50分钟在线观看-中文字幕中文字幕永久免费-中文字幕在线永久 | 佐藤遥希在线播放一二区-佐山爱巨大肥臀在线-佐山爱痴汉theav-佐良娜被爆漫画羞羞漫画-麻豆视频传媒二区-麻豆视频传媒 | 国产一级一级一级成人毛片-国产一级一级片-国产一级网站-国产一级特黄在线播放-午夜影院一区二区三区-午夜影院小视频 | 欧美精品综合-欧美精品自拍-欧美精品专区免费观看-欧美精品在线一区二区三区-欧美精品在线一区-欧美精品在线视频观看 | 日韩成人A片一区二区三区-日韩成人黄色片-日韩成人在线视频-日韩吃奶摸下AA片免费观看-日韩大片免费看-日韩大片在线观看 | 日本乱中文字幕系列-日本伦理黄色大片在线观看网站-日本伦理片网站-日本伦理网站-韩国三级视频网站-韩国三级视频在线观看 | 无人 视频 在线观看免费高清-无吗视频-无码中文字幕日韩专区-无码中文字幕乱在线观看-国产区1-国产情侣一区二区 | 九九99国产精品视频-九草在线观看-九草在线-精品综合久久久久久蜜月-亚洲厕所凸轮-亚洲播播 | 欧美精品久久久亚洲,欧美国产第一页,亚欧成人一区二区,久久精品一区二区影院,国产成人精品综合久久久,亚洲国产精品一区 | 国产精品成av人在线观看片-国产精品成久久久久三级-国产精品成久久久久三级四虎-国产精品成久久久久三级无码-国产精品成年片在线观看-国产精品成人 | 丰满的日本护士-丰满的日本护士xxx-丰满的少妇一区二区三区免费观看-丰满的熟妇岳中文字幕-丰满多毛少妇做爰视频-丰满多毛少妇做爰视频爽爽和R | 日本漫画母亲口工子全彩-日本漫画大全无翼乌-日本妈妈在线观看中文字幕-日本妈妈xxxx-操他射他影院-操老太太的逼 | 亚洲尹人-亚洲伊人久久综合影院2021-亚洲伊人久久在-亚洲伊人久久网-免费网站在线观看人数在哪里看的-免费网站在线观看国产v片 | 久久久精品国产sm最大网站-久久久精品国产免大香伊-久久久精品国产免费A片胖妇女-久久久精品久久久久久96-久久久精品久久久久久久久久久-久久久精品久久日韩一区综合 | 日本漫画母亲口工子全彩-日本漫画大全无翼乌-日本妈妈在线观看中文字幕-日本妈妈xxxx-操他射他影院-操老太太的逼 | 久草免费在线观看-久草免费在线色站-久草免费在线视频-久草免费资源-天堂网中文在线-天堂网中文字幕 | 欧美成人tv在线观看免费-欧美成人tv-欧美成人h版影院在线播放-欧美成人h版影片在线观看-欧美成人h版白雪公主-欧美成人h版 |