999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于證據(jù)相似性度量的沖突性區(qū)間證據(jù)融合方法

2012-07-25 04:11:32馮海山徐曉濱文成林
電子與信息學報 2012年4期
關(guān)鍵詞:融合方法

馮海山 徐曉濱 文成林

(杭州電子科技大學自動化學院 杭州 310018)

1 引言

經(jīng)典的 Dempster-Shafer(D-S)證據(jù)理論以其對不確定性、非精確性信息綜合處理方面的優(yōu)勢,在信息融合、模式識別和故障診斷等領(lǐng)域中得到了廣泛的應(yīng)用[1-13]。基于該理論的信息融合方法,能夠?qū)崿F(xiàn)空間或時間上冗余信息和互補信息的融合,獲得被測對象的一致性描述,有效降低決策中的不確定性。經(jīng)典的D-S證據(jù)理論融合的是單值的基本概率分布函數(shù)(BPA),然而,單值BPA對不確定性信息或模糊信息的度量并不全面,且可能會丟失很多有用的信息,故文獻[6-9]提出了區(qū)間型 BPA(IBPA),即區(qū)間證據(jù)。該種形式的證據(jù)能夠較為全面地度量信息的不確定性,而且符合人的常性思維和主觀概念[11]。所以,區(qū)間證據(jù)理論成為當今不確定性理論研究的熱點,并集中于以下兩個方向:一是區(qū)間證據(jù)融合規(guī)則的構(gòu)造方法,二是區(qū)間證據(jù)的獲取方法[12,13]。在研究方向一中,文獻[10]首次定義了IBPA的乘法和加法運算,但是其運算規(guī)則中引入了主觀性因子,故運算結(jié)果具有較大的不確定性和主觀性。Denoeux[6,7]構(gòu)造二次規(guī)劃模型以融合多個IBPA,給出IBPA有效性和歸一化準則,首次提出 Demspter區(qū)間證據(jù)組合規(guī)則。該規(guī)則中的融合和歸一化是分開進行的,所以得到的融合結(jié)果并不是全局最優(yōu)的,而是次優(yōu)的,且融合結(jié)果的置信區(qū)間較寬,不易用其進行決策。基于Denoeux定義的IBPA有效性和歸一化準則,Wang等人[8,9]提出最優(yōu)的 Demspter區(qū)間證據(jù)組合公式。對于有效且歸一化的 IBPA,該公式遍歷所有滿足約束條件的單值BPA進行融合,然后對單值融合結(jié)果求極值得到融合后的IBPA。該融合過程是一步進行的,從而保證了融合結(jié)果的最優(yōu)性;而研究方向二,主要是討論如何在實際中獲取區(qū)間證據(jù),這也是應(yīng)用融合規(guī)則的前提。在旋轉(zhuǎn)機械故障診斷中,基于傳感器提供的故障特征數(shù)據(jù),Xu等人[11]提出一種獲取區(qū)間型診斷證據(jù)的方法,并利用最優(yōu) Demspter區(qū)間證據(jù)組合公式將多個傳感器提供的證據(jù)融合,并通過融合結(jié)果定位故障。與基于BPA融合的故障定位方法相比,利用區(qū)間型證據(jù)可顯著提高融合診斷系統(tǒng)的確診率,并且文中通過診斷實例進一步驗證其所提出的證據(jù)獲取及融合診斷方法的可靠性和準確性。

但是,在實際應(yīng)用中,由于傳感器測量誤差、環(huán)境噪聲干擾和監(jiān)測數(shù)據(jù)不完整等因素,使得從不同傳感器獲取的單值或區(qū)間證據(jù)之間常常會存在沖突[5]。當今,對于單值BPA沖突處理方法的研究已經(jīng)相對成熟。其中,一類是修改組合規(guī)則的方法,它們以不同的方式重新分配空集的置信度賦值,但其組合規(guī)則只適用于解決某些具體問題,缺乏普適性[2,3];另一類是修改模型的方法,亦即修改證據(jù)的方法。該類方法主要考慮到證據(jù)之間的關(guān)聯(lián)性,用證據(jù)權(quán)重修改原始證據(jù),然后還利用原始的Dempster組合規(guī)則對修改后的證據(jù)進行融合[2-5],此類方法從證據(jù)之間相似性的角度反映了沖突的本質(zhì)。

在區(qū)間證據(jù)的融合當中,Demspter區(qū)間證據(jù)組合公式的融合機理在于,選取區(qū)間證據(jù)中的單值BPA,用經(jīng)典的 Demspter組合規(guī)則融合,進而統(tǒng)計生成區(qū)間型BPA。雖然經(jīng)典的Demspter融合規(guī)則具有聚焦作用,但在融合高沖突區(qū)間證據(jù)時,融合后的BPA賦值過于分散,并不能很好地聚焦于某些焦元上,導致IBPA的區(qū)間寬度過大,不易用于決策。如何結(jié)合單值證據(jù)沖突處理方法,以減小融合后所得區(qū)間證據(jù)的區(qū)間寬度,使得區(qū)間上下邊界聚焦于同一個焦元上,利于決策,這是區(qū)間證據(jù)沖突融合中一個十分重要的新問題。

針對以上問題,本文基于證據(jù)相似性度量,提出一種沖突性區(qū)間證據(jù)融合的新方法。首先定義擴展型Pignistic概率轉(zhuǎn)換,將區(qū)間證據(jù)轉(zhuǎn)換為區(qū)間型Pignistic概率。利用區(qū)間模糊集的歸一化歐式距離,求取區(qū)間型Pignistic概率之間的相似性,以此確定兩兩證據(jù)間的相似度矩陣,從中獲取區(qū)間證據(jù)的置信度。然后,基于該置信度對原始的區(qū)間證據(jù)進行加權(quán)平均得到新的區(qū)間證據(jù),利用 Demspter區(qū)間證據(jù)組合公式對其進行融合。該方法可以有效地減弱高沖突性區(qū)間證據(jù)在組合規(guī)則中的作用,從而減小融合后所得區(qū)間證據(jù)的寬度,最終可降低決策中的不確定性。最后通過多個典型算例驗證了經(jīng)沖突處理后再對區(qū)間證據(jù)進行融合,要比直接融合能夠產(chǎn)生更為合理和可靠的結(jié)果。

本文章節(jié)安排如下:第2節(jié)介紹了區(qū)間證據(jù)理論的基礎(chǔ),并以實例展示了直接對沖突證據(jù)進行融合所產(chǎn)生的不合理結(jié)果,并分析了造成該結(jié)果的原因;第3節(jié)基于經(jīng)典的Pignistic概率轉(zhuǎn)換,提出擴展型 Pignistic概率轉(zhuǎn)換,用于區(qū)間證據(jù) Pignistic的轉(zhuǎn)換,基于此求取區(qū)間證據(jù)間的歐氏距離,用其度量它們之間的相似度及置信度;第4節(jié)通過多個典型算例驗證本文所提方法,能夠有效地降低沖突區(qū)間證據(jù)融合后的區(qū)間寬度,而且融合后所得區(qū)間證據(jù)的上下界可聚焦于同一個焦元上,更利于決策。

2 區(qū)間證據(jù)理論

令Θ= {θ1,θ2,…,θn}是一個非空的有限子集,其中的元素兩兩互斥,稱其為辨識框架。Θ中包含了所有可能的命題。

定義1區(qū)間基本概率賦值 (Interval Basic Probability Assignment, IBPA)[6-8,11-13]對于Θ的N個子集Ai(i=1,2,…,N),則其IBPA定義為

其中0 ≤ai≤bi≤ 1 。m(Ai)是m(Ai)中的一個元素,若IBPA同時滿足

定義2歸一化準則[8,9]若m是一個有效的IBPA,且m(Ai)=[ai,bi](0 ≤ai≤bi≤ 1 ),如果ai和bi同時滿足

其中i,j=1,2,…,N,則稱m為歸一化的IBPA。

m可能是一個有效的 IBPA,但不一定是歸一化的IBPA,定義2是判斷一個有效IBPA是否歸一化的判據(jù),若未歸一化,則用式(4)進行歸一化處理

式(4)對有效的IBPA進行歸一化處理,來減小區(qū)間的寬度,降低冗余,得到簡潔等效的IBPA。

定義3Demspter區(qū)間證據(jù)組合公式[8,9]若m1和m2是有效且歸一化的 IBPA,分別為[ai,bi](0≤ai≤m1(Ai) ≤bi≤ 1 ,i=1,2,…,N)和[cj,dj](0≤cj≤m2(Aj) ≤dj≤ 1 ,j=1,2,…,N),融合結(jié)果標記為m1⊕m2,其為區(qū)間值

其中(m1⊕m2)-(C)和(m1⊕m2)+(C)分別為如下融合公式的最小值與最大值。

對于滿足定義1和定義2的IBPA,式(6)中的融合規(guī)則遍歷 IBPA區(qū)間中所有滿足約束條件的單值BPA并融合,從而保證了融合結(jié)果的最優(yōu)性。

下面通過例子說明 Demspter區(qū)間證據(jù)組合公式的應(yīng)用,在同一辨識框架Θ= {θ1,θ2,θ3}下有兩組有效且歸一化的IBPA,如下所示:

將兩組區(qū)間證據(jù)用 Demspter區(qū)間證據(jù)組合公式融合后結(jié)果如表1所示。

由表1可知,若區(qū)間證據(jù)同時支持某焦元(如第1組證據(jù)),當利用Demspter區(qū)間證據(jù)組合公式融合時,采用的單值BPA也會同時支持某焦元,這時的融合結(jié)果能夠正確聚焦;若區(qū)間證據(jù)間存在高沖突(如第2組證據(jù)中強烈支持θ1,而強烈支持θ3),利用Demspter組合公式所得的融合結(jié)果中,焦元θ1和θ3的區(qū)間寬度較大且相等,無法決策。

3 基于區(qū)間證據(jù)相似性的沖突證據(jù)度量及融合

在處理單值沖突證據(jù)的修改模型法中,利用權(quán)重(置信度)度量某個證據(jù)和其它證據(jù)之間的沖突程度。也就是說,如果其它證據(jù)支持某個證據(jù)時,則說明該證據(jù)比較可信,其所占權(quán)重較大,對融合結(jié)果的影響也較大;反之,如果某個證據(jù)與其它證據(jù)間的沖突較大時,則該證據(jù)的可信度較低,其所占權(quán)重就較低,對融合結(jié)果的影響也較小。這種方法充分考慮到了證據(jù)之間的相互關(guān)聯(lián)性,減少引起沖突的證據(jù)權(quán)重,提高最終融合結(jié)果的合理性和可靠性。并且應(yīng)用廣泛易于理解,符合客觀情況[2-5]。本文將該思想推廣到?jīng)_突性區(qū)間證據(jù)的處理上。首先,基于單值Pignistic概率轉(zhuǎn)換,提出擴展型Pignistic概率轉(zhuǎn)換,將其應(yīng)用到區(qū)間證據(jù)上,轉(zhuǎn)換后結(jié)果記為IBetP,其中,I表示區(qū)間,BetP代表Pignstic概率轉(zhuǎn)換,此時的IBetP是一個區(qū)間而非單值;其次,計算每個辨識框架單元素的IBetP,并確定區(qū)間IBetP間的相互距離,以間接度量區(qū)間證據(jù)間的沖突程度;最后分析確定區(qū)間證據(jù)的相似度、置信度,更新原始區(qū)間證據(jù),用加權(quán)平均后新區(qū)間證據(jù)替代原始區(qū)間證據(jù),再用Demspter區(qū)間證據(jù)組合公式進行融合,使融合結(jié)果收斂到正確的命題,以便做出決策。

3.1 擴展型Pignistic概率轉(zhuǎn)換

基于期望效用理論,Smets定義了Pignistic概率函數(shù)。其基本思想是,在辨識框架Θ上進行Pignistic轉(zhuǎn)換,將基本概率函數(shù)轉(zhuǎn)換成Pignistic概率函數(shù)。

表1 兩組區(qū)間證據(jù)融合結(jié)果

定義4經(jīng)典Pignistic轉(zhuǎn)換 設(shè)m是在辨識框架Θ下的一個BPA,相應(yīng)的Pignistic概率函數(shù)BetPm:Θ→[0,1]定義為

對于單值的BPA,經(jīng)過以上Pignistic轉(zhuǎn)換后的BetPm仍然是單值的,這種形式的概率轉(zhuǎn)換不能直接用于IBPA上,需對其進行擴展,即區(qū)間證據(jù)的Pignistic轉(zhuǎn)換,表示為IBetPm。

定義5擴展型Pignistic轉(zhuǎn)換 設(shè)m是在辨識框架Θ上的一個有效且歸一化的IBPA,記為[ai,bi](0 ≤ai≤m(Ai) ≤bi≤ 1 ,i=1,2,… ,N), 則 其Pignistic概率函數(shù)IBetPm定義為

例如,在辨識框架Θ= {θ1,θ2,θ3}下的有效且歸一化的IBPA,將其用擴展型 Pignistic轉(zhuǎn)換結(jié)果表示為IBetPm,如表2所示。

表2 利用擴展型Pignistic轉(zhuǎn)換后的結(jié)果

3.2 基于區(qū)間歐氏距離的區(qū)間相似性度量

這里采用區(qū)間歸一化歐氏距離公式度量IBetP間的相似性[14],以間接度量IBetP所對應(yīng)的區(qū)間證據(jù)間的相似性。

設(shè)m1,m2是同一辨識框架Θ中的兩個IBPA,Θ含有n個完備且相互獨立的元素,記為Θ= {θi,i=1,2,…,n},經(jīng)過擴展型Pignistic轉(zhuǎn)換后,分別記為IBet和 I Bet,則其距離為

若兩區(qū)間證據(jù)間的距離用式(10)度量,則其相似度可以表示為

若融合系統(tǒng)含有N個區(qū)間證據(jù),通過式(9),式(10)得到區(qū)間證據(jù)mi與mj之間的距離,并由式(11)確定兩者的相似度,記為Sij,則可構(gòu)造相似度矩陣為

則每個區(qū)間證據(jù)的支持度為

支持度Sup(mi)反映的是被其它區(qū)間證據(jù)所支持的程度,是相似性測度的函數(shù)。如果一個區(qū)間證據(jù)與其它證據(jù)相似程度較高,則認為它們相互支持的程度也較高;反之亦然。

對于證據(jù)體mi的置信度為

若有N條區(qū)間證據(jù),將得到的新區(qū)間證據(jù)用Demspter區(qū)間組合公式融合N-1次,得到的最終融合區(qū)間,相比直接用Demspter區(qū)間組合公式獲得的IBPA,區(qū)間寬度較窄,并且區(qū)間上下界聚焦于同一個焦元的作用很明顯。其主要原因是,由于本文方法考慮了區(qū)間證據(jù)之間的相似程度,各個區(qū)間證據(jù)因相互支持度的不同獲得不同的權(quán)重,如果一個區(qū)間證據(jù)被其它區(qū)間證據(jù)所支持的程度較高,該證據(jù)就越可信,對最后的融合結(jié)果影響程度越大。相反,如果某一證據(jù)與其它證據(jù)是高沖突的,它的權(quán)重就很低,對最終的融合結(jié)果影響程度較小。

4 實例分析

為了說明本文方法的有效性,這里給出3個高沖突區(qū)間證據(jù)融合的典型例子,并設(shè)定辨識框架均為Θ= {θ1,θ2,θ3},且所給出的區(qū)間證據(jù)均是有效且歸一化的。例 1是關(guān)于單元素賦值的區(qū)間證據(jù)。例 2中含有對辨識框架下單元素與全集賦值的區(qū)間證據(jù)。例3中含有對辨識框架下單元素與子集賦值的區(qū)間證據(jù)。

例1設(shè)3個區(qū)間證據(jù)m1,m2和m3為

用式(8)對3個區(qū)間證據(jù)進行擴展型Pignistic轉(zhuǎn)換,并用式(9),式(10)計算得到兩兩區(qū)間證據(jù)間的距離及相似度,以此構(gòu)造成相似度矩陣為

根據(jù)式(11)-式(13)得到各個區(qū)間證據(jù)體的置信度分別為

采用加權(quán)平均后,再用 Demspter區(qū)間證據(jù)組合公式進行融合,并與直接采用 Demspter區(qū)間證據(jù)組合公式融合的結(jié)果進行對比,如表3所示。

表3 兩種融合方法的結(jié)果對比

通過表3可以看出,當區(qū)間證據(jù)存在高沖突時,即m1強烈支持θ1,而m2卻強烈的支持θ3,m3卻不能明顯支持一方。若直接用 Demspter區(qū)間證據(jù)組合公式,融合結(jié)果的區(qū)間寬度很大,難以決策。而本文方法充分考慮了區(qū)間證據(jù)之間的相互關(guān)聯(lián)性,降低沖突性區(qū)間證據(jù)在整個融合系統(tǒng)的權(quán)重,減小其在融合過程中的作用,所得的最終融合結(jié)果的區(qū)間寬度較窄,并且區(qū)間上下界同時聚焦于同一個焦元上的作用明顯,從而提高了決策能力。

例2設(shè)定只對辨識框架下單元素和全集賦值的區(qū)間證據(jù)mi(i=1,2,…,5)如表4所示。

表4 單元素和全集賦值的區(qū)間證據(jù)

本文方法與直接用 Demspter區(qū)間證據(jù)組合公式融合結(jié)果進行對比,如表5所示,其中列出了區(qū)間證據(jù)依次融合的過程。

通過表4知,區(qū)間證據(jù)m1,m4和m5同時支持θ1,m2支持θ2,m3未明確支持任何一個焦元,可知m2與m1,m4,m5存在著高沖突。若直接用Demspter區(qū)間證據(jù)組合公式融合,則最終的融合結(jié)果的區(qū)間寬度較大,不易決策。而采用本文方法時,由于采用區(qū)間證據(jù)的置信度度量區(qū)間證據(jù)的相互支持度,因此提高了m1,m4和m5的權(quán)重,同時減小m2的權(quán)重,并弱化m3的權(quán)重,使最終融合結(jié)果的區(qū)間寬度很窄,并且區(qū)間的上下界同時聚焦于θ1的作用更為明顯。

表5 兩種融合方法的結(jié)果對比

例3設(shè)定只對辨識框架中單元素和子集賦值的區(qū)間證據(jù)mi(i=1,2,…,5)如表6所示。

表6 單元素和子集賦值的區(qū)間證據(jù)

本文方法與直接用 Demspter區(qū)間證據(jù)組合公式融合結(jié)果進行對比,如表7所示,其中列出了區(qū)間證據(jù)依次融合的過程。

從表6中可知,區(qū)間證據(jù)m1,m3,m4和m5均支持θ1,只有m2支持θ2,依理推知m2是干擾證據(jù),并且最終融合結(jié)果也是應(yīng)該支持θ1,但從融合結(jié)果表7中可以得知,直接用Demspter區(qū)間證據(jù)組合公式得到的最終融合區(qū)間寬度很大,根本無法做出決策。而采用本文方法,通過降低干擾證據(jù)m2的權(quán)重,最終的融合結(jié)果區(qū)間寬度很窄,并且區(qū)間上下界同時聚焦于θ1上的作用更加明顯。

5 結(jié)論

本文提出了一種基于證據(jù)相似性度量的沖突性區(qū)間證據(jù)融合方法,考慮焦元之間及區(qū)間證據(jù)之間的相互關(guān)聯(lián)性,提出擴展型Pignistic概率轉(zhuǎn)換;并基于區(qū)間歸一化歐氏距離公式來間接度量區(qū)間證據(jù)的相似性,分析區(qū)間證據(jù)相互之間的沖突;然后確定各個區(qū)間證據(jù)的置信度,加權(quán)平均原始區(qū)間證據(jù)以得到新的區(qū)間證據(jù),再用 Demspter區(qū)間證據(jù)組合公式進行融合,減少了沖突證據(jù)在融合過程的作用,降低了干擾對最終融合結(jié)果的影響。最終,融合結(jié)果的區(qū)間寬度較窄,并且區(qū)間上下界同時聚焦于同一個焦元上的作用更加明顯,從而提高了沖突性區(qū)間證據(jù)融合的合理性和可靠性。

表7 兩種融合方法的結(jié)果對比

[1]徐曉濱, 文成林, 王迎昌. 基于模糊故障特征信息的隨機集度量信息融合診斷方法[J]. 電子與信息學報, 2009, 31(7):1635-1640.

Xu Xiao-bin, Wen Cheng-lin, and Wang Ying-chang.Information fusion algorithm of fault diagnosis based on random set metrics of fuzzy fault features [J].Journal of Electronics&Information Technology, 2009, 31(7):1635-1640.

[2]Liu Wei-ru. Analyzing the degree of conflict among belief functions [J].Artificial Intelligence, 2006, 170(11): 909-924.

[3]Yager R R. On the fusion of imprecise uncertainty measures using belief structures[J].Information Sciences, 2011, 181(15):3199-3209.

[4]Deng Yong, Shi Wen-kang,et al.. Combining belief function based on distance of evidence [J].Decision Support Systems,2004, 38(3): 489-493.

[5]徐曉濱, 王玉成, 文成林. 基于診斷證據(jù)可靠性評估的信息融合故障診斷方法[J]. 控制理論與應(yīng)用, 2011, 28(4): 504-510.

Xu Xiao-bin, Wang Yu-cheng, and Wen Cheng-lin.Information-fusion method for fault diagnosis based on reliability evaluation of evidence [J].Control Theory and Applications, 2011, 28(4): 504-510.

[6]Denoeux T. Reasoning with imprecise belief structures [J].International Journal of Approximate Reasoning, 1999, 20(1):79-111.

[7]Denoeux T. Modelling vague belief using fuzzy-valued belief structures [J].Fuzzy Sets and Systems, 2000, 116(2): 167-199.

[8]Wang Ying-ming, Yang Jian-bo,et al.. On the combination and normalization of interval-valued belief structures [J].Information Sciences, 2007, 177(5): 1230-1247.

[9]Wang Ying-ming, Yang Jian-bo,et al.. The evidential reasoning approach for multiple attribute decision analysis using interval belief degrees [J].European Journal of Operational Research, 2006, 175(1): 35-66.

[10]Lee E S and Zhu Qing. An interval Dempster-Shafer approach [J].Computers and Mathematics with Applications,1992, 24(7): 89-95.

[11]Xu Xiao-bin, Feng Hai-shan,et al.. An information fusion method of fault diagnosis based on interval basic probability assignment [J].Chinese Journal of Electronics, 2011, 20(2):255-260.

[12]Su Zhi-gang, Wang Pei-hong,et al.. Maximal confidence intervals of the interval-valued belief structure and applications[J].InformationSciences, 2011, 181(9):1700-1721.

[13]Fu Chao and Yang Shan-lin. Analyzing the applicability of Dempster’s rule to the combination of interval-valued belief structures [J].Expert Systems with Applications, 2011, 38(4):4291-4301.

[14]Zeng Wen-yi and Guo Ping. Normalized distance, similarity measure, inclusion measure and entropy of interval-valued fuzzy sets and their relationship [J].Information Sciences,2008, 178(5): 1334-1342.

猜你喜歡
融合方法
一次函數(shù)“四融合”
村企黨建聯(lián)建融合共贏
融合菜
從創(chuàng)新出發(fā),與高考數(shù)列相遇、融合
寬窄融合便攜箱IPFS500
《融合》
學習方法
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
賺錢方法
主站蜘蛛池模板: 无码免费视频| 欧美激情成人网| 素人激情视频福利| 国产香蕉一区二区在线网站| 午夜限制老子影院888| 91丨九色丨首页在线播放| 日韩无码白| 日韩123欧美字幕| 97久久免费视频| 国产精品深爱在线| 日韩第一页在线| 在线日本国产成人免费的| 国产乱子伦一区二区=| 无码中文字幕精品推荐| 国产欧美在线观看精品一区污| 亚洲综合精品香蕉久久网| 国产日本一线在线观看免费| 欧美午夜一区| 亚洲丝袜中文字幕| 午夜国产大片免费观看| 毛片基地美国正在播放亚洲 | 免费无码网站| 久久www视频| 成人在线天堂| 人妻中文久热无码丝袜| 国产成人禁片在线观看| 免费久久一级欧美特大黄| 亚洲日韩国产精品无码专区| 91美女视频在线| 国产福利一区视频| 国内a级毛片| 国产在线视频福利资源站| 日韩区欧美国产区在线观看| 色香蕉影院| 一级毛片基地| 久久久波多野结衣av一区二区| 亚洲欧美在线精品一区二区| 99精品免费欧美成人小视频| 精品成人免费自拍视频| 欧美成人二区| 欧美啪啪网| 国产偷国产偷在线高清| 一级毛片在线播放免费观看| 在线播放精品一区二区啪视频| 日韩美女福利视频| 91久久精品国产| 国产成人高精品免费视频| 国产噜噜噜视频在线观看| 亚洲天堂网在线播放| 日韩专区第一页| 九九九精品成人免费视频7| 日韩在线欧美在线| 国产成人a在线观看视频| 毛片三级在线观看| 五月婷婷丁香综合| 丁香亚洲综合五月天婷婷| 一本视频精品中文字幕| 色妞www精品视频一级下载| 久久久精品国产SM调教网站| 无码高潮喷水在线观看| 久久黄色毛片| 久草热视频在线| 亚洲男人的天堂久久香蕉网| 日韩第九页| 国产精品开放后亚洲| 五月天综合婷婷| 欧美激情首页| 欧美成人综合在线| 成人国产精品2021| 日本三级黄在线观看| 亚洲天堂网2014| 久久久久青草大香线综合精品 | 怡春院欧美一区二区三区免费| 黄色网在线| 久热中文字幕在线观看| 亚洲男女天堂| 日韩一区精品视频一区二区| 国产精品亚欧美一区二区三区 | 伊人久久久久久久久久| 久草网视频在线| 午夜小视频在线| 免费久久一级欧美特大黄|