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

    基于f-散度的復(fù)雜系統(tǒng)涌現(xiàn)度量方法

    2017-07-05 14:07:55何新華陸皖麟
    關(guān)鍵詞:散度窗體概率分布

    屈 強, 何新華, 陸皖麟

    (裝甲兵工程學(xué)院信息工程系, 北京 100072)

    ?

    基于f-散度的復(fù)雜系統(tǒng)涌現(xiàn)度量方法

    屈 強, 何新華, 陸皖麟

    (裝甲兵工程學(xué)院信息工程系, 北京 100072)

    針對目前復(fù)雜系統(tǒng)涌現(xiàn)研究領(lǐng)域缺少量化手段的問題,提出了一種基于f-散度的復(fù)雜系統(tǒng)涌現(xiàn)度量方法。首先,從有序性角度重新定義了涌現(xiàn),提出了涌現(xiàn)強度的概念;然后,從可度量性、收斂性和靈敏度角度,分析比較了常見的5種f-散度,認(rèn)為Hellinger(Hel)散度最適合涌現(xiàn)度量;最后,給出了Hel散度的近似計算方法和涌現(xiàn)度量流程,并通過螞蟻尋食實例仿真驗證了該涌現(xiàn)度量方法的有效性。

    涌現(xiàn)度量; 涌現(xiàn)性; 復(fù)雜系統(tǒng);f-散度; Hel散度

    系統(tǒng)是由具有相同目標(biāo)的單元相互作用構(gòu)成的整體[1],可分為簡單系統(tǒng)、無組織系統(tǒng)和復(fù)雜系統(tǒng)。復(fù)雜系統(tǒng)具有非線性、自組織性、自適應(yīng)性和涌現(xiàn)性等特性,而涌現(xiàn)性是認(rèn)識復(fù)雜系統(tǒng)的重點。涌現(xiàn)性是系統(tǒng)整體具有,分解到局部就不存在的屬性、特征、行為和現(xiàn)象等[2]。目前,涌現(xiàn)還主要停留在定性研究階段,而定量研究是涌現(xiàn)判斷、預(yù)測和控制的前提。

    GABBAI等[3]最早提出用熵來度量涌現(xiàn),DE WOLF等[4]提出用離散熵來分析無中心系統(tǒng)的自組織涌現(xiàn)行為,然而有序性并不等同于涌現(xiàn)性,因此采用有序性指標(biāo)(熵)來直接度量涌現(xiàn)還有待商榷。MNIF等[5]提出采用不同時間點的熵差來度量自組織系統(tǒng)的涌現(xiàn),HOLZER等[6]提出用熵值的商度量涌現(xiàn),兩者的研究本質(zhì)是用有序性增量作為涌現(xiàn)的度量,是對文獻(xiàn)[3-4]研究的改進(jìn),具有很強的借鑒意義。

    有序性增量既可用熵差、熵值的商來度量,也可用f-散度來度量。如:FISCH等[7]引入f-散度來度量涌現(xiàn),提出采用Hellinger(Hel)散度度量涌現(xiàn),并比較了Hel散度與熵差度量涌現(xiàn)的優(yōu)缺點,但未詳細(xì)說明選用Hel散度的原因。針對上述問題,筆者系統(tǒng)論述了基于f-散度度量復(fù)雜系統(tǒng)涌現(xiàn)的理論依據(jù),通過分析比較證明了Hel散度最適合涌現(xiàn)度量,給出了Hel散度近似計算方法和涌現(xiàn)度量流程,并進(jìn)行了有效性驗證。

    1 基本概念

    1.1f-散度

    f-散度首先由CSISZAR等[8]和ALI等[9]在概率論中提出,用以度量概率分布之間的差異。假設(shè)概率分布P(x)和Q(x)定義在概率空間Ω上,且P(x)對Q(x)連續(xù),則對于滿足f(1)=0的凸函數(shù)f(x),Q(x)到P(x)的f-散度定義為

    (1)

    對于式(1),如果概率空間Ω上存在一個參考概率分布μ,使得P(x)和Q(x)均相對于μ連續(xù),則其對應(yīng)的概率密度函數(shù)p(x)、q(x)可分別寫為dP(x)=p(x)dμ和dQ(x)=q(x)dμ。于是,f-散度定義可改寫為[10]

    (2)

    在實際應(yīng)用中,f-散度定義式常轉(zhuǎn)變?yōu)殡x散形式。如果f(x)為[0,∞)區(qū)間上的凸函數(shù),且概率分布集合P={p1,p2,…,pN}和Q={q1,q2,…,qN}滿足

    (3)

    則Q(x)到P(x)的f-散度可定義為

    (4)

    1.2 涌現(xiàn)和涌現(xiàn)強度

    涌現(xiàn)是指系統(tǒng)單元在環(huán)境刺激下相互作用,從無序(或準(zhǔn)無序)狀態(tài)向有序狀態(tài)轉(zhuǎn)變的過程。無序和有序是系統(tǒng)單元狀態(tài)的統(tǒng)計學(xué)特征,其無序來自于系統(tǒng)單元行為的隨機(jī)性,而有序是涌現(xiàn)的結(jié)果。環(huán)境刺激是涌現(xiàn)的外因,單元間的相互作用是涌現(xiàn)的內(nèi)動力。系統(tǒng)單元在環(huán)境作用下定向改變狀態(tài),并向內(nèi)傳導(dǎo),使得狀態(tài)呈現(xiàn)統(tǒng)計意義上的有序。系統(tǒng)涌現(xiàn)出的有序狀態(tài)定義為涌現(xiàn)特征。

    假設(shè)系統(tǒng)S從t1到t2時刻產(chǎn)生涌現(xiàn),某項系統(tǒng)特征x的取值分布從無序變?yōu)橛行颍瑒t有序性增量E(x)為

    E(x)=Rt2(x)-Rt1(x),

    (5)

    式中:Rt1(x)和Rt2(x)分別表示t1和t2時刻系統(tǒng)特征x取值分布的有序性。

    顯然,有序性增量越大,涌現(xiàn)特征越明顯,越容易被觀察分析。因此,可將系統(tǒng)特征有序性增量定義為涌現(xiàn)強度(用E表示),以表征涌現(xiàn)特征的明顯程度。系統(tǒng)單元之間的不同作用將產(chǎn)生不同的涌現(xiàn)特征,作用范圍不同,涌現(xiàn)強度相差較大,為了化繁為簡,設(shè)定涌現(xiàn)強度門限,重點抓住強涌現(xiàn)特征,忽略弱涌現(xiàn)特征。

    系統(tǒng)特征取值完全無序時,近似服從正態(tài)分布。涌現(xiàn)發(fā)生時,系統(tǒng)特征的概率分布將發(fā)生變化,涌現(xiàn)強度越大,概率分布變化量越大。概率分布變化量可用f-散度來表征,因此涌現(xiàn)強度又可用f-散度來表示,即

    (6)

    式中:Pt1(x)、Qt2(x)分別為t1和t2時刻系統(tǒng)特征x取值的概率分布。

    2 f-散度的種類與選取

    2.1f-散度的種類

    由于凸函數(shù)f(x)有很多函數(shù)形式,根據(jù)式(4)可知f-散度也有多種類型,其中5種常見的f-散度,如表1所示[11-14]。

    表1 常見的f-散度

    2.2f-散度的選取

    根據(jù)可度量性、收斂性和靈敏度來比較選擇適用于度量系統(tǒng)涌現(xiàn)特征的f-散度。

    1)可度量性

    設(shè)X為一個集合,R為實數(shù)集,函數(shù)d(·)滿足映射關(guān)系:X×X→R。若對于集合X中的任意元素x、y、z,函數(shù)d(·)滿足:

    (1) 非負(fù)性:d(x,y)≥ 0,當(dāng)且僅當(dāng)x=y時d(x,y)=0。

    (2) 對稱性:d(x,y)=d(y,x)。

    (3) 三角不等式:d(x,y)≤d(x,z)+d(z,y)。

    則稱函數(shù)d(·)為集合X上的完全可度量函數(shù)。不完全滿足以上3個條件的函數(shù)稱為廣義可度量函數(shù)或廣義距離(Generalized-Metrics)[15-16]。

    5種f-散度中,全變差散度、Jensen-Shannon(JS)散度和Hel散度完全滿足上述3個條件,為完全可度量函數(shù),具有完全可度量性;Kullback-Leibler(KL)散度和χ2散度則不滿足對稱性和三角不等式,屬于廣義距離,不具備完全可度量性。

    2)收斂性

    Hel散度的取值范圍為[0,1],用于后續(xù)數(shù)據(jù)處理時程序收斂速度更快。而全變差散度和JS散度的取值范圍為[0,∞),需要進(jìn)行歸一化處理才能確保程序收斂。

    3)靈敏度

    靈敏度是指函數(shù)變量取值微小變化時函數(shù)值的相對變化量。它是函數(shù)值相對于變量值的靈敏程度,具體表達(dá)式為

    (7)

    將5種f-散度對應(yīng)的f(x)函數(shù)歸一化處理后代入式(7)進(jìn)行比較,可知:全變差散度靈敏度最高,KL散度靈敏度最低,Hel散度、JS散度和χ2散度靈敏度居中。

    綜上分析可知:5種f-散度中,Hel散度具有完全可度量性、較好的收斂性和較高的靈敏度,是涌現(xiàn)度量的最佳選擇。

    3 基于Hel散度的涌現(xiàn)度量

    3.1 Hel散度的近似計算

    對于離散分布的Hel散度,其定義式為

    (8)

    Hel散度計算的關(guān)鍵是進(jìn)行概率分布的估計。概率分布估計的方法分為參數(shù)法和非參數(shù)法,其中:參數(shù)法包括最大似然估計法和貝葉斯估計法等;非參數(shù)法包括Parzen窗法和K近鄰法[17]等。本文采用Parzen窗法進(jìn)行概率分布估計,其流程為:首先,建立系統(tǒng)實例,從初始時刻t1開始,連續(xù)以等間隔d采集v個樣本(v>>hN,其中hN為窗體長度);然后,采用滑動數(shù)據(jù)窗估計p和q。

    Parzen窗法的基本形式為

    (9)

    式中:X={x1,x2,…,xN},為樣本空間,其中N為樣本總數(shù);D為數(shù)據(jù)維數(shù);φ(·) 為窗函數(shù),通常選取正態(tài)窗函數(shù),即

    (10)

    將式(10)代入式(9)可得

    (11)

    根據(jù)式(8)、(11)可近似計算Hel散度。

    Parzen窗估計的滑動加窗方法分為2種:1)采用一窗固定、一窗滑動的方法,將第1個窗體q固定在時間軸的t1時刻,以t2時刻為中心滑動第2個窗體p,如圖1(a)所示;2)第1個窗體q加在t1時刻,第2窗體p加在t2時刻,然后以等間距d同時滑動q和p兩窗,如圖1(b)所示。

    圖1 Parzen窗估計的2種滑動加窗方法

    3.2 涌現(xiàn)度量流程

    假設(shè)系統(tǒng)S由m個單元組成,t1時刻開始觀察,t2時刻出現(xiàn)L種系統(tǒng)宏觀特征,進(jìn)行涌現(xiàn)度量,其流程如下:

    1) 構(gòu)建系統(tǒng)S的特征值集X={x1,x2,…,xL},選取其中一個特征xk(k=1,2,…,L);

    2) 分別對t1、t2時刻特征xk的取值進(jìn)行統(tǒng)計,估計其概率分布p(xk)、q(xk);

    3) 根據(jù)式(8)計算t1至t2時刻特征xk的涌現(xiàn)強度Ek;

    4)設(shè)定判據(jù)門限e,若Ek>e,說明系統(tǒng)涌現(xiàn)出了特征xk,可納入系統(tǒng)S的涌現(xiàn)特征集(系統(tǒng)涌現(xiàn)特征構(gòu)成的全集M,M?X);

    5)重復(fù)上述步驟,直至所有特征度量完畢;

    6)根據(jù)涌現(xiàn)強度大小對涌現(xiàn)特征集中的涌現(xiàn)特征進(jìn)行排序,供后續(xù)分析使用。

    4 仿真實例

    在空間全局信息未知的情況下,蟻群能夠共同找到食物和最短路徑,將食物搬運回巢穴。蟻群尋食過程并不是所有螞蟻個體行為的簡單總和,而是螞蟻相互作用后產(chǎn)生的涌現(xiàn),采用NetLogo仿真軟件對此過程進(jìn)行建模仿真。

    以蟻群分布為系統(tǒng)特征,假設(shè)蟻群由100只螞蟻組成,在初始時刻t1,蟻群無組織地隨機(jī)搜索,螞蟻位置坐標(biāo)在觀察空間呈隨機(jī)分布,如圖2(a)所示;經(jīng)過信息素作用,蟻群出現(xiàn)自組織和協(xié)同,t2時刻涌現(xiàn)出如圖2(b)所示的螞蟻位置斑圖。

    圖2 不同時刻螞蟻位置分布圖

    由于螞蟻位置坐標(biāo)為二維數(shù)據(jù),為計算方便,將其轉(zhuǎn)化換為一維樣本數(shù)據(jù),具體方法為:將觀測空間劃分為9塊區(qū)域,相當(dāng)于分辨率設(shè)置為3×3;假設(shè)落在同一塊區(qū)域的螞蟻位置坐標(biāo)相同,依次對觀察區(qū)域進(jìn)行編號(如圖3所示),則落在第j(j=1,2,…,9)塊區(qū)域中螞蟻的個數(shù)為xj,此時樣本空間為X={x1,x2,…,x9},數(shù)據(jù)維數(shù)D=1,N=9。

    圖3 觀察空間劃分與編號

    采用式(11)中的Parzen窗法估計t1、t2時刻螞蟻在觀察空間上的概率分布Pt1(x)和Qt2(x),選取窗體長度hN=5;然后,將Pt1(x)和Qt2(x)代入式(8)計算得出t2時刻的Hel散度值。同理,在[t1,t2]內(nèi)等間隔劃分15個仿真步長(共計17個步數(shù)),分別計算出相應(yīng)時間點的Hel散度值,如圖4所示。

    圖4 Hel散度變化圖

    由圖4可知:1)開始未放置食物,因而螞蟻的分布坐標(biāo)較為隨機(jī),系統(tǒng)為無序;2)從步數(shù)7開始放置食物,Hel散度值明顯增大,這是因為螞蟻在尋找、搬運食物的過程中釋放信息素,系統(tǒng)因信息素的作用而趨向于有序;3)根據(jù)經(jīng)驗,選取判斷門限e=0.1,則從步數(shù)9時刻開始涌現(xiàn)強度可被分辨;4)從步數(shù)11時刻開始,Hel散度達(dá)到最大并趨于穩(wěn)定,即最大涌現(xiàn)強度為0.16。

    上述分析表明:Hel散度能夠有效度量復(fù)雜系統(tǒng)的涌現(xiàn)特征。當(dāng)然,坐標(biāo)空間上的螞蟻位置斑圖只是蟻群系統(tǒng)的涌現(xiàn)特征之一,若選擇不同的屬性進(jìn)行觀測,則可能還有其他的涌現(xiàn)特征,這里不再逐一度量。

    5 結(jié)論

    從可度量性、收斂性和靈敏度角度分析比較了5種常見的f-散度,認(rèn)為Hel散度最適合涌現(xiàn)度量。給出了基于Hel散度的涌現(xiàn)度量方法流程和近似計算方法,并通過蟻群仿真實驗驗證了有效性。

    基于f-散度的度量方法也存在很多問題,例如:使用Parzen窗估計散度時,需要用戶根據(jù)經(jīng)驗設(shè)定窗體長度;在不同的觀察分辨率下對同一涌現(xiàn)過程進(jìn)行度量,將得到不同的涌現(xiàn)強度。因此,如何選取窗體長度和觀察分辨率等是今后需要重點研究的問題。

    [1] KONSTANTINOS P T,PAUL D C.A comprehensive basis for systems engineering theory[C]∥Proceedings of 2014 IEEE International Systems Conference.Ottawa,ON,Canada:IEEE,2014:139-143.

    [2] 苗東升.論涌現(xiàn)[J].河池學(xué)院學(xué)報,2008,28(1):6-12.

    [3] GABBAI J M E,YIN H,WRIGHT W A,et al.Self organization,emergence and multi-agent systems[C]∥Proceedings of 2005 International Conference on Neural Networks and Brain.Beijing,China:IEEE,2005:1858-1863.

    [4] DE WOLF T,SAMAEY G,HOLVOET T,et al.Decentralised autonomic computing:analysing self organising emergent behaviour using advanced numerical methods[C]∥Proceedings of Second International Conference on Autonomic Computing.Seattle,Wa-shington,D.C.,U.S.:IEEE Computer Society,2005:52-63.

    [5] MNIF M,MULLER-SCHLOER C.Quantitative emergence[J].IEEE mountain workshop on adaptive and learning systems,2011(1):39-52.

    [6] HOLZER R,DE MEER H,BETTSTETTER C.On autonomy and emergence in self organizing systems[J].Lecture notes in computer science,2008,5343(1):157-169.

    [7] FISCH D,JANICKE M,SICK B,et al.Quantitative emergence:a refined approach based on divergence measures[C]∥Proceedings of 2010 4th IEEE International Conference on self adaptive and self organizing Systems.San Francisco,California,America:IEEE,2010:94-103.

    [8] CSISZAR I.Information-type measures of difference of probability distributions and indirect observations[J].Studia scientiarum mathematicarum hungarica,1967,2(3):229-318.

    [9] ALI S M,SILVEY S D.A general class of coefficients of divergence of one distribution from another[J].Journal of the royal statistical society,1966,28(1):131-142.

    [10] 任鳳英,李興斯.基于f-散度的一致風(fēng)險度量[J].大連理工大學(xué)學(xué)報,2013,53(5):772-776.

    [11] JOHN J,JOHNSON IV,ANDREAS T,et al.A theory of emergence and entropy in systems of systems[J].Procedia computer science,2013,20:283-289.

    [12] ZENG J S,KRUGER U,GELUK J,et al.Detecting abnormal situations using the Kullback-Leibler divergence[J].Automatica,2014,50(11):2777-2786.

    [13] VICTOR G C,ROCIO R,ENRIQUE A.Class distribution estimation based on the Hellinger distance[J].Information sciences,2013,218:146-164.

    [14] BRIESCH D M,JAJA C E,MOORE T J,et al.Divergence measures tool:an introduction with brief tutorial[R].Adelphi:Army Research Lab,2014.

    [15] KOPPERMAN R,PAJOOHESH H.Intrinsic generalized metrics[J].Algebra universalis,2012,67(1):1-18.

    [16] GIANNOPOULOS P,VELTKAMP R C.A Pseudo-metric for weighted point sets[J].Lecture notes in computer science,2002,2352(1):89-114.

    [17] 李慶,胡捍英.支持向量預(yù)選取的K邊界鄰近法[J].電路與系統(tǒng)學(xué)報,2013(2):91-96.

    (責(zé)任編輯: 尚菲菲)

    A New Approach to Measure the Emergence of Complex System Based onf-divergence

    QU Qiang, HE Xin-hua, LU Wan-lin

    (Department of Information Engineering, Academy of Armored Force Engineering, Beijing 100072, China)

    In complex system research, there is still no mature means of quantitative description of emergence till now. To solve this problem, an emergence quantitative measurement method based onf-divergence is proposed. Firstly, the concept of emergence is re-defined in the view of system order, and the concept of emergence power is given for the first time. Secondly, analysis and comparison for the five kinds of commonf-divergence from the perspective of measurability, convergence and sensibility, Hel-divergence is considered as the best measure for quantitative emergence. Lastly, the process of emergence measurement based on Hel-divergence and approximate calculation method of Hel-divergence are given. The effectiveness of thef-divergence method is verified by simulation example of ant colony.

    quantitative emergence; emergent property; complex system;f-divergence; Hel-divergence

    1672-1497(2017)03-0106-05

    2016-12-23

    軍隊科研計劃項目

    屈 強(1982-),男,博士研究生。

    N941.4

    A

    10.3969/j.issn.1672-1497.2017.03.020

    猜你喜歡
    散度窗體概率分布
    帶勢加權(quán)散度形式的Grushin型退化橢圓算子的Dirichlet特征值的上下界
    離散型概率分布的ORB圖像特征點誤匹配剔除算法
    具有部分BMO系數(shù)的非散度型拋物方程的Lorentz估計
    試談Access 2007數(shù)據(jù)庫在林業(yè)檔案管理中的應(yīng)用
    檔案天地(2019年5期)2019-06-12 05:12:02
    H型群上一類散度形算子的特征值估計
    關(guān)于概率分布函數(shù)定義的辨析
    科技視界(2016年19期)2017-05-18 10:18:46
    基于概率分布的PPP項目風(fēng)險承擔(dān)支出測算
    H?rmander 向量場上散度型拋物方程弱解的Orlicz估計
    基于LayeredWindow的異形窗體局部刷新
    中文信息(2014年2期)2014-03-06 23:49:14
    巧設(shè)WPS窗體控件讓表格填寫更規(guī)范
    91精品国产国语对白视频| 国产av国产精品国产| 下体分泌物呈黄色| 日本91视频免费播放| 19禁男女啪啪无遮挡网站| 交换朋友夫妻互换小说| 国产男人的电影天堂91| 亚洲精品日本国产第一区| 精品人妻一区二区三区麻豆| 亚洲成色77777| 亚洲精品久久久久久婷婷小说| 爱豆传媒免费全集在线观看| 人妻人人澡人人爽人人| 亚洲av在线观看美女高潮| 欧美黑人欧美精品刺激| 在线观看免费日韩欧美大片| 别揉我奶头~嗯~啊~动态视频 | 黑人猛操日本美女一级片| 亚洲精品自拍成人| 这个男人来自地球电影免费观看 | 亚洲精品国产av蜜桃| 国产精品一区二区精品视频观看| 嫩草影视91久久| 国产又爽黄色视频| 看免费成人av毛片| 久久99一区二区三区| 黄片无遮挡物在线观看| 欧美日韩福利视频一区二区| 午夜av观看不卡| 18在线观看网站| 黄片无遮挡物在线观看| 亚洲av福利一区| 校园人妻丝袜中文字幕| 夜夜骑夜夜射夜夜干| 一区二区三区激情视频| 80岁老熟妇乱子伦牲交| 高清视频免费观看一区二区| 国产成人精品久久久久久| 在线观看免费午夜福利视频| 国产男女内射视频| 王馨瑶露胸无遮挡在线观看| 啦啦啦啦在线视频资源| 久久国产精品男人的天堂亚洲| 国产在线免费精品| 只有这里有精品99| 18禁裸乳无遮挡动漫免费视频| 人人澡人人妻人| 午夜福利影视在线免费观看| 麻豆精品久久久久久蜜桃| 亚洲国产精品国产精品| 亚洲欧美成人精品一区二区| av.在线天堂| 制服诱惑二区| svipshipincom国产片| 一本大道久久a久久精品| 亚洲精品国产区一区二| 搡老岳熟女国产| 日本黄色日本黄色录像| 午夜福利网站1000一区二区三区| 免费人妻精品一区二区三区视频| 日韩制服骚丝袜av| 99热网站在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲色图 男人天堂 中文字幕| 男女边摸边吃奶| 又大又黄又爽视频免费| 亚洲av福利一区| 久久ye,这里只有精品| 午夜久久久在线观看| 成年女人毛片免费观看观看9 | a级毛片黄视频| 无限看片的www在线观看| 777米奇影视久久| 一本色道久久久久久精品综合| 日本av免费视频播放| 99热国产这里只有精品6| 这个男人来自地球电影免费观看 | 王馨瑶露胸无遮挡在线观看| 日韩精品有码人妻一区| 欧美日韩综合久久久久久| 人人妻人人澡人人爽人人夜夜| 亚洲精品美女久久久久99蜜臀 | 国产深夜福利视频在线观看| 国产野战对白在线观看| 国产成人精品久久二区二区91 | 国产精品久久久久成人av| 新久久久久国产一级毛片| 欧美激情高清一区二区三区 | kizo精华| 久久精品久久久久久噜噜老黄| 91精品国产国语对白视频| av线在线观看网站| 精品国产乱码久久久久久男人| xxx大片免费视频| 亚洲综合精品二区| 一级爰片在线观看| 国产成人精品久久久久久| 亚洲精品国产一区二区精华液| 欧美日韩精品网址| 18禁国产床啪视频网站| 69精品国产乱码久久久| 热re99久久国产66热| 精品亚洲乱码少妇综合久久| 男女之事视频高清在线观看 | 两个人看的免费小视频| 嫩草影视91久久| 国产有黄有色有爽视频| 欧美精品av麻豆av| 亚洲美女搞黄在线观看| 纵有疾风起免费观看全集完整版| 欧美另类一区| 国产精品久久久久久久久免| 肉色欧美久久久久久久蜜桃| 国产成人啪精品午夜网站| 一边摸一边做爽爽视频免费| 成人免费观看视频高清| 少妇人妻久久综合中文| 国产 精品1| 国产精品久久久人人做人人爽| 午夜日本视频在线| 啦啦啦视频在线资源免费观看| 黑人巨大精品欧美一区二区蜜桃| 国产福利在线免费观看视频| 午夜免费鲁丝| 丁香六月天网| 国产精品成人在线| 91精品伊人久久大香线蕉| 精品亚洲成国产av| 日韩av在线免费看完整版不卡| 制服人妻中文乱码| 涩涩av久久男人的天堂| 一区二区三区精品91| 国产成人免费无遮挡视频| 97在线人人人人妻| 麻豆av在线久日| 人人澡人人妻人| 日韩制服丝袜自拍偷拍| 国产伦人伦偷精品视频| 国产成人欧美在线观看 | 女性被躁到高潮视频| 国产亚洲一区二区精品| 国产一区有黄有色的免费视频| 久久久久久久久久久免费av| 在线观看免费午夜福利视频| 天美传媒精品一区二区| 久久久久国产精品人妻一区二区| 老熟女久久久| 午夜免费男女啪啪视频观看| 老汉色av国产亚洲站长工具| 免费在线观看完整版高清| 久久久欧美国产精品| 熟女av电影| 性高湖久久久久久久久免费观看| 欧美最新免费一区二区三区| 女人久久www免费人成看片| 岛国毛片在线播放| 亚洲人成77777在线视频| 美女大奶头黄色视频| 精品久久久久久电影网| 国产人伦9x9x在线观看| 伦理电影免费视频| av不卡在线播放| 婷婷色av中文字幕| 亚洲三区欧美一区| 极品少妇高潮喷水抽搐| 国产不卡av网站在线观看| videosex国产| 亚洲精品aⅴ在线观看| 国产一区亚洲一区在线观看| 男男h啪啪无遮挡| 一二三四在线观看免费中文在| 婷婷色av中文字幕| 中文字幕人妻丝袜一区二区 | 好男人视频免费观看在线| 女的被弄到高潮叫床怎么办| 一级爰片在线观看| 99re6热这里在线精品视频| 免费观看人在逋| 免费高清在线观看视频在线观看| 成人国语在线视频| 黄色怎么调成土黄色| 五月天丁香电影| 人人妻,人人澡人人爽秒播 | 亚洲精品自拍成人| 国产一区二区在线观看av| 午夜久久久在线观看| 热99国产精品久久久久久7| 免费久久久久久久精品成人欧美视频| 欧美日韩视频高清一区二区三区二| 色播在线永久视频| 老司机在亚洲福利影院| 国产乱人偷精品视频| 蜜桃国产av成人99| 午夜福利乱码中文字幕| 夫妻性生交免费视频一级片| 香蕉丝袜av| 美女脱内裤让男人舔精品视频| 一级片免费观看大全| av网站在线播放免费| 国产成人免费无遮挡视频| 国产乱来视频区| 侵犯人妻中文字幕一二三四区| 亚洲少妇的诱惑av| 国产成人精品久久久久久| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| av线在线观看网站| 成人免费观看视频高清| 欧美人与性动交α欧美精品济南到| 欧美日韩视频高清一区二区三区二| 精品国产超薄肉色丝袜足j| 在线观看一区二区三区激情| 80岁老熟妇乱子伦牲交| 国产不卡av网站在线观看| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美在线一区| 七月丁香在线播放| 久久国产亚洲av麻豆专区| 我要看黄色一级片免费的| 天天躁夜夜躁狠狠躁躁| 亚洲精品久久久久久婷婷小说| 亚洲欧美一区二区三区久久| 亚洲欧美成人精品一区二区| 国产成人91sexporn| 国产精品熟女久久久久浪| 黑人巨大精品欧美一区二区蜜桃| 大陆偷拍与自拍| 最近最新中文字幕大全免费视频 | 成年动漫av网址| 人人妻人人澡人人爽人人夜夜| 街头女战士在线观看网站| 成人国语在线视频| 黄色怎么调成土黄色| 蜜桃国产av成人99| 国产伦理片在线播放av一区| 欧美亚洲 丝袜 人妻 在线| 国产精品一区二区在线不卡| 超碰成人久久| 亚洲专区中文字幕在线 | 日韩成人av中文字幕在线观看| 亚洲欧美清纯卡通| 丝袜喷水一区| 肉色欧美久久久久久久蜜桃| 午夜老司机福利片| 亚洲精品久久午夜乱码| 在线观看www视频免费| 日韩一区二区视频免费看| 国产成人精品无人区| 看免费av毛片| 日韩电影二区| 久久精品久久久久久噜噜老黄| 久久久久国产精品人妻一区二区| 亚洲一码二码三码区别大吗| 日韩欧美精品免费久久| 亚洲欧美成人精品一区二区| 欧美中文综合在线视频| 高清在线视频一区二区三区| 美女主播在线视频| 久久精品久久久久久久性| 国产免费福利视频在线观看| 美女福利国产在线| 免费在线观看完整版高清| 久久久久国产一级毛片高清牌| 女的被弄到高潮叫床怎么办| 1024视频免费在线观看| 国产精品熟女久久久久浪| av有码第一页| 多毛熟女@视频| av在线app专区| 午夜av观看不卡| 老司机影院成人| 啦啦啦 在线观看视频| www.熟女人妻精品国产| 91成人精品电影| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美精品综合一区二区三区| 国产熟女欧美一区二区| 精品少妇内射三级| 在线天堂中文资源库| 国产人伦9x9x在线观看| 欧美黑人精品巨大| 最近最新中文字幕免费大全7| 亚洲精品中文字幕在线视频| 国产极品天堂在线| 亚洲欧美一区二区三区国产| 精品福利永久在线观看| 无限看片的www在线观看| 国产在线视频一区二区| 国产精品女同一区二区软件| 另类亚洲欧美激情| 国产精品欧美亚洲77777| 日本91视频免费播放| 久久精品久久久久久噜噜老黄| 一二三四在线观看免费中文在| 一级毛片黄色毛片免费观看视频| 嫩草影院入口| 亚洲成色77777| 亚洲精品久久午夜乱码| 亚洲在久久综合| 国产亚洲最大av| 国产日韩欧美在线精品| 可以免费在线观看a视频的电影网站 | 赤兔流量卡办理| 成人国产av品久久久| 亚洲国产精品一区二区三区在线| 我要看黄色一级片免费的| 久久精品人人爽人人爽视色| 久久 成人 亚洲| 亚洲成av片中文字幕在线观看| 亚洲欧洲日产国产| 午夜免费观看性视频| 国产欧美日韩综合在线一区二区| av电影中文网址| 国产成人午夜福利电影在线观看| 亚洲国产av新网站| 老司机在亚洲福利影院| 99热全是精品| 色网站视频免费| 日韩中文字幕视频在线看片| 国产精品嫩草影院av在线观看| 午夜免费男女啪啪视频观看| 人成视频在线观看免费观看| 两个人看的免费小视频| 亚洲av成人不卡在线观看播放网 | 欧美人与善性xxx| 久久久久精品久久久久真实原创| 国产免费一区二区三区四区乱码| 亚洲精品久久成人aⅴ小说| 成人手机av| 美国免费a级毛片| 精品人妻在线不人妻| 别揉我奶头~嗯~啊~动态视频 | 国产免费一区二区三区四区乱码| 亚洲av电影在线进入| 亚洲色图综合在线观看| 男女床上黄色一级片免费看| 国产精品一区二区在线不卡| 久久av网站| 久久久精品94久久精品| 久久精品aⅴ一区二区三区四区| 大码成人一级视频| 日韩 欧美 亚洲 中文字幕| 精品国产乱码久久久久久小说| 曰老女人黄片| 看十八女毛片水多多多| 日本vs欧美在线观看视频| 国产成人a∨麻豆精品| 亚洲成人免费av在线播放| 国产成人精品久久二区二区91 | 久久精品久久久久久久性| 90打野战视频偷拍视频| 日本vs欧美在线观看视频| av不卡在线播放| 日韩欧美三级三区| 国产精品亚洲av一区麻豆| 又紧又爽又黄一区二区| 亚洲视频免费观看视频| 亚洲人成电影观看| 亚洲精品国产色婷婷电影| 黑人巨大精品欧美一区二区蜜桃| 黑丝袜美女国产一区| 成人三级做爰电影| 亚洲va日本ⅴa欧美va伊人久久| 亚洲伊人色综图| 露出奶头的视频| 亚洲精品国产一区二区精华液| 亚洲国产欧美日韩在线播放| 久久中文字幕一级| 中出人妻视频一区二区| 亚洲 欧美 日韩 在线 免费| 99国产综合亚洲精品| 免费在线观看亚洲国产| 9热在线视频观看99| 国产精品,欧美在线| 99久久国产精品久久久| 亚洲欧美激情在线| 亚洲欧美一区二区三区黑人| 亚洲精品国产一区二区精华液| 中文亚洲av片在线观看爽| 久久久国产精品麻豆| 很黄的视频免费| 夜夜看夜夜爽夜夜摸| 19禁男女啪啪无遮挡网站| 中文字幕色久视频| 国产一级毛片七仙女欲春2 | 精品国产一区二区三区四区第35| 久久九九热精品免费| 国内精品久久久久精免费| 午夜福利18| 久久久久久久久久久久大奶| 熟女少妇亚洲综合色aaa.| 欧美日韩黄片免| 欧美日韩中文字幕国产精品一区二区三区 | 成年版毛片免费区| 日韩 欧美 亚洲 中文字幕| 成人手机av| 欧美一级a爱片免费观看看 | 禁无遮挡网站| 国产一区二区在线av高清观看| 又大又爽又粗| 精品高清国产在线一区| 亚洲片人在线观看| 一级,二级,三级黄色视频| 丰满人妻熟妇乱又伦精品不卡| 不卡一级毛片| 波多野结衣av一区二区av| 国产亚洲av嫩草精品影院| 亚洲va日本ⅴa欧美va伊人久久| 亚洲熟妇熟女久久| 欧美中文日本在线观看视频| 女人被躁到高潮嗷嗷叫费观| 夜夜爽天天搞| 麻豆一二三区av精品| 亚洲国产看品久久| 国产av又大| 在线观看舔阴道视频| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人一区二区免费高清观看 | 午夜老司机福利片| 激情视频va一区二区三区| 亚洲黑人精品在线| 一区二区三区国产精品乱码| 99精品欧美一区二区三区四区| 免费久久久久久久精品成人欧美视频| 亚洲伊人色综图| 久久香蕉激情| 九色亚洲精品在线播放| 丝袜美腿诱惑在线| 男男h啪啪无遮挡| 中文字幕人妻丝袜一区二区| 国产成人啪精品午夜网站| 亚洲自拍偷在线| bbb黄色大片| 精品国产超薄肉色丝袜足j| 亚洲五月天丁香| 亚洲自偷自拍图片 自拍| 黑丝袜美女国产一区| 成人欧美大片| 天天添夜夜摸| 色婷婷久久久亚洲欧美| 麻豆av在线久日| 亚洲一区二区三区色噜噜| 精品熟女少妇八av免费久了| 女人高潮潮喷娇喘18禁视频| 校园春色视频在线观看| 日本a在线网址| 又黄又爽又免费观看的视频| 亚洲人成电影观看| 成人国产一区最新在线观看| 欧美成人免费av一区二区三区| 丰满的人妻完整版| 丁香六月欧美| 亚洲七黄色美女视频| 亚洲午夜精品一区,二区,三区| 亚洲五月天丁香| av视频在线观看入口| 午夜激情av网站| 两性夫妻黄色片| 91字幕亚洲| 中文字幕色久视频| 18禁国产床啪视频网站| 可以在线观看毛片的网站| 亚洲九九香蕉| 九色亚洲精品在线播放| 久久久国产成人免费| 最好的美女福利视频网| 国内久久婷婷六月综合欲色啪| 99国产精品免费福利视频| 看片在线看免费视频| 叶爱在线成人免费视频播放| 国产精品免费视频内射| 在线观看免费午夜福利视频| 国产又爽黄色视频| 国产精品一区二区免费欧美| 国产精品自产拍在线观看55亚洲| 国产精品一区二区精品视频观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产精品999在线| 欧美国产日韩亚洲一区| 国产精品九九99| 超碰成人久久| 久久性视频一级片| 女人高潮潮喷娇喘18禁视频| www.熟女人妻精品国产| 成人免费观看视频高清| 亚洲一区中文字幕在线| 99在线视频只有这里精品首页| 欧美日本亚洲视频在线播放| 国产高清videossex| 国产三级黄色录像| 免费在线观看视频国产中文字幕亚洲| 午夜免费激情av| 波多野结衣高清无吗| 欧美另类亚洲清纯唯美| 99国产精品免费福利视频| 精品久久久精品久久久| www.www免费av| 纯流量卡能插随身wifi吗| 久久草成人影院| 亚洲成人国产一区在线观看| 看片在线看免费视频| 国产成人精品久久二区二区91| 妹子高潮喷水视频| 十分钟在线观看高清视频www| 他把我摸到了高潮在线观看| 在线国产一区二区在线| 久久天躁狠狠躁夜夜2o2o| 老熟妇仑乱视频hdxx| 老司机午夜福利在线观看视频| 欧美激情 高清一区二区三区| 精品国产国语对白av| 人成视频在线观看免费观看| 国产精品综合久久久久久久免费 | 97碰自拍视频| 亚洲国产精品久久男人天堂| 亚洲中文av在线| 亚洲午夜精品一区,二区,三区| 男人操女人黄网站| 久久国产乱子伦精品免费另类| 国产aⅴ精品一区二区三区波| 亚洲狠狠婷婷综合久久图片| 精品不卡国产一区二区三区| 极品人妻少妇av视频| 在线永久观看黄色视频| 国产亚洲精品综合一区在线观看 | 狠狠狠狠99中文字幕| 一区二区三区激情视频| 国产区一区二久久| 国产欧美日韩一区二区三| 午夜日韩欧美国产| 这个男人来自地球电影免费观看| 国产av在哪里看| 午夜久久久久精精品| 黑丝袜美女国产一区| 一边摸一边做爽爽视频免费| e午夜精品久久久久久久| 亚洲视频免费观看视频| 女人爽到高潮嗷嗷叫在线视频| 一本久久中文字幕| 精品久久久久久,| 国产在线精品亚洲第一网站| 久久人妻熟女aⅴ| 亚洲精品国产区一区二| 成人三级做爰电影| 国产一区二区在线av高清观看| 曰老女人黄片| 三级毛片av免费| 精品国产美女av久久久久小说| 亚洲av成人一区二区三| 午夜免费激情av| 天天躁狠狠躁夜夜躁狠狠躁| 免费观看精品视频网站| 日韩av在线大香蕉| www.自偷自拍.com| 亚洲国产精品久久男人天堂| 久久久精品欧美日韩精品| 国产精品亚洲一级av第二区| 成人国语在线视频| 久久精品影院6| 老熟妇乱子伦视频在线观看| 人人妻人人澡人人看| 免费av毛片视频| 9191精品国产免费久久| a在线观看视频网站| 久久国产乱子伦精品免费另类| 欧美一级a爱片免费观看看 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品国产区一区二| 欧美另类亚洲清纯唯美| 国产成人精品在线电影| 亚洲国产高清在线一区二区三 | 天堂影院成人在线观看| 香蕉丝袜av| 黄色 视频免费看| 法律面前人人平等表现在哪些方面| 久久久久久大精品| 国产亚洲精品久久久久久毛片| 亚洲成国产人片在线观看| 99在线人妻在线中文字幕| 午夜激情av网站| 欧美色欧美亚洲另类二区 | 国产亚洲精品一区二区www| 大香蕉久久成人网| 国产成人系列免费观看| 可以在线观看毛片的网站| 1024视频免费在线观看| 波多野结衣av一区二区av| 亚洲 欧美一区二区三区| 久久精品91无色码中文字幕| 色在线成人网| 亚洲精品av麻豆狂野| 这个男人来自地球电影免费观看| 亚洲国产精品999在线| 亚洲精品av麻豆狂野| 又大又爽又粗| 免费少妇av软件| 中文字幕高清在线视频| 十八禁人妻一区二区| 久久国产精品男人的天堂亚洲| 成人精品一区二区免费| 国产精华一区二区三区| 男女下面进入的视频免费午夜 | 麻豆成人av在线观看| 欧美中文日本在线观看视频| 女生性感内裤真人,穿戴方法视频| 久久久久国产精品人妻aⅴ院| 国产国语露脸激情在线看| 免费在线观看完整版高清| 激情在线观看视频在线高清| 久99久视频精品免费| 99久久久亚洲精品蜜臀av| 国产高清videossex| а√天堂www在线а√下载| 波多野结衣高清无吗| 午夜免费鲁丝|