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

    基于聚類的差分隱私民航旅客數(shù)據(jù)發(fā)布算法

    2022-03-21 10:32:54丁建立杜天天
    關(guān)鍵詞:信息

    丁建立,杜天天

    (中國民航大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,天津 300300)

    0 引 言

    隱私安全問題通常會(huì)阻礙用戶數(shù)據(jù)分享,甚至阻礙分析從數(shù)據(jù)挖掘的價(jià)值信息[1]。數(shù)據(jù)挖掘在生活中有許多應(yīng)用,聚類是一種廣泛使用的數(shù)據(jù)挖掘方法之一,應(yīng)用數(shù)據(jù)挖掘通常都假定可以自由訪問數(shù)據(jù)集,但這是不現(xiàn)實(shí)的。因此需要做好隱私保護(hù),解決由于數(shù)據(jù)發(fā)布導(dǎo)致的個(gè)人隱私泄露這一問題。隱私保護(hù)的數(shù)據(jù)發(fā)布已有多種方法。其中差分隱私技術(shù)可以保證在兩個(gè)僅相差一條記錄的數(shù)據(jù)集上運(yùn)行算法能夠產(chǎn)生相似的輸出。本文目標(biāo)是實(shí)現(xiàn)發(fā)布民航旅客數(shù)據(jù)的同時(shí)最大程度地保護(hù)隱私。差分隱私提供強(qiáng)大的隱私保證。大數(shù)據(jù)時(shí)代,數(shù)據(jù)無處不在,大量有價(jià)值的信息隱藏在數(shù)據(jù)中等待著研究者挖掘。數(shù)據(jù)分析者可通過對(duì)數(shù)據(jù)整理分析[2]提取有價(jià)值信息,使人們生活更加便利美好。然而在增加人民幸福感的同時(shí),其中的問題也相應(yīng)暴露出,在搜集數(shù)據(jù)的同時(shí)會(huì)攜帶大量隱私信息,比如,在分析不文明旅客時(shí),會(huì)收集其身份信息、聯(lián)系方式等等。如果向公眾直接發(fā)布這些信息,有可能對(duì)用戶構(gòu)成相當(dāng)大的威脅。因此,數(shù)據(jù)發(fā)布隱私保護(hù)[3]的重難點(diǎn)為保證發(fā)布的數(shù)據(jù)可用的同時(shí)不會(huì)泄露隱私信息。

    1 研究現(xiàn)狀

    1.1 聚類的研究現(xiàn)狀

    聚類是數(shù)據(jù)挖掘中的重要工具,在數(shù)據(jù)分析、信息檢索和文本挖掘等領(lǐng)域具有許多應(yīng)用。它旨在將有限的、未標(biāo)記的數(shù)據(jù)集劃分為幾個(gè)自然子集,同一簇內(nèi)的數(shù)據(jù)相近,而來自不同簇的數(shù)據(jù)互不相同。為了實(shí)現(xiàn)該目的,現(xiàn)在已經(jīng)提出了許多聚類算法,例如:基于劃分、基于網(wǎng)絡(luò)以及基于模型的聚類等等。這些聚類算法中的大多數(shù)都需要事先指定集群的數(shù)量或隱式集群的數(shù)量控制參數(shù)。對(duì)于某些應(yīng)用,可以根據(jù)用戶的專業(yè)知識(shí)來估計(jì)集群的數(shù)量。但是,在許多情況下,不能確定數(shù)據(jù)集的簇?cái)?shù)。然而,聚類數(shù)量會(huì)大大影響聚類結(jié)果的好壞。因此,確定數(shù)據(jù)集中的簇?cái)?shù)(通常標(biāo)記為k)是聚類分析中的一個(gè)基本問題?,F(xiàn)有很多估計(jì)k值的研究[4]?;跀?shù)據(jù)類型的差異,這些方法通常可以歸類為用于數(shù)值數(shù)據(jù)、分類數(shù)據(jù)和混合數(shù)據(jù)的聚類算法。在數(shù)值領(lǐng)域,文獻(xiàn)[5]提出了針對(duì)海量數(shù)據(jù)的k-means問題的有效近似方法。將整個(gè)數(shù)據(jù)集遞歸地劃分為少量子集,每個(gè)子集均以其代表(質(zhì)心)和權(quán)重(基數(shù))為特征,然后將k-means算法的加權(quán)版本應(yīng)用于此類局部表示,這樣可以大大減少計(jì)算出的距離數(shù)。對(duì)于分類數(shù)據(jù),文獻(xiàn)[6]提出k-modes聚類的性能對(duì)初始聚類中心的選擇特別敏感。從離群值檢測(cè)的角度考慮了k-modes聚類的初始化。通過使用基于傳統(tǒng)的基于距離的離群值檢測(cè)技術(shù)和基于分區(qū)熵的離群值檢測(cè)技術(shù)來計(jì)算每個(gè)對(duì)象的離群度。在初始化過程中,采用了新的距離度量標(biāo)準(zhǔn)-加權(quán)匹配距離度量標(biāo)準(zhǔn)。

    文獻(xiàn)[7]中有幾種算法可以對(duì)混合數(shù)據(jù)進(jìn)行聚類。但是,所有這些算法都需要預(yù)先直接或間接指定聚類數(shù)。本文提出確定混合數(shù)據(jù)集中簇?cái)?shù)的有效方法。該方法由經(jīng)過改進(jìn)的k-prototype組成。

    1.2 隱私保護(hù)的研究現(xiàn)狀

    隱私數(shù)據(jù)保護(hù)自隱私出現(xiàn)一直被人們所研究,并提出了很多保護(hù)方法,k-anonymity及其改進(jìn)的算法[8]是其中較早且經(jīng)典的研究成果。k-anonymity保護(hù)數(shù)據(jù)中的隱私是通過保證每個(gè)記錄與至少k-1個(gè)其它記錄無法區(qū)分。但是在攻擊者可能獲得背景知識(shí)的前提下,很容易受到攻擊,之后提出的l-diversity和t-closeness,可以預(yù)防多樣性攻擊,但仍會(huì)受到其它攻擊。差分隱私模型通過向數(shù)據(jù)中添加一定量的噪聲,保證發(fā)布的數(shù)據(jù)不會(huì)泄露隱私信息,無需考慮攻擊者的背景知識(shí),因此差分隱私有更加強(qiáng)大的隱私保證。

    近年來,差分隱私保護(hù)技術(shù)已經(jīng)成為了一個(gè)備受關(guān)注的研究熱點(diǎn)。差分隱私[9]被越來越多地用作數(shù)據(jù)分析選擇的隱私保護(hù)技術(shù),使數(shù)據(jù)既能被挖掘有價(jià)值信息的同時(shí)又保護(hù)了用戶個(gè)人隱私。差分隱私數(shù)據(jù)發(fā)布有交互式框架和非交互式框架兩種方式,在差分隱私初期,被應(yīng)用于研究者在數(shù)據(jù)庫上查詢數(shù)據(jù)時(shí),返回的查詢數(shù)據(jù)滿足差分隱私,然而這種形式的差分隱私應(yīng)用對(duì)數(shù)據(jù)有嚴(yán)格的限制,交互式場(chǎng)景下,隨著查詢次數(shù)增多,所得的數(shù)據(jù)偏差會(huì)越來越大,最后所查詢數(shù)據(jù)甚至完全沒有參考價(jià)值,這導(dǎo)致它只能接受有限數(shù)量的查詢。在此弊端的影響下人們研究出了非交互式場(chǎng)景下的差分隱私數(shù)據(jù)發(fā)布,這里的數(shù)據(jù)可對(duì)整個(gè)數(shù)據(jù)集進(jìn)行差分隱私處理,使得可以發(fā)布滿足差分隱私的數(shù)據(jù)集,這些數(shù)據(jù)集保證了基本的數(shù)據(jù)統(tǒng)計(jì)特征,研究者可以利用這些統(tǒng)計(jì)結(jié)果進(jìn)行分析和處理,然而如果想要進(jìn)一步分析,可能使用這些數(shù)據(jù)集不能夠得出精確的分析結(jié)果,這是由于差分隱私需要在數(shù)據(jù)集中加入大量噪聲所造成的,這些噪聲既保護(hù)了數(shù)據(jù),相應(yīng)的又破壞了數(shù)據(jù),使得數(shù)據(jù)的可用性降低,其分析價(jià)值降低。而后研究表明,合理分配隱私預(yù)算,在實(shí)現(xiàn)保護(hù)隱私的前提下,適當(dāng)降低查詢敏感度可以大大提高差分隱私數(shù)據(jù)的可用性。目前對(duì)于混合型數(shù)據(jù)發(fā)布的研究相對(duì)較少,而在實(shí)際生活中大部分?jǐn)?shù)據(jù)集為混合型數(shù)據(jù)集。因此,研究針對(duì)混合型(數(shù)據(jù)類型和分類類型)數(shù)據(jù)集的差分隱私數(shù)據(jù)發(fā)布方法更加具有現(xiàn)實(shí)意義,可實(shí)現(xiàn)人們對(duì)于隱私保護(hù)的需求。文獻(xiàn)[10]研究了兩種方法在差分隱私k-means聚類中的有效性。提出了一種非交互式方法EUGkM,該方法發(fā)布了k-means聚類的差分隱私模型,驗(yàn)證了EUGkM算法的有效性。文獻(xiàn)[11]中介紹了最廣泛使用的聚類k-means容易陷入局部最優(yōu)。并提出傳統(tǒng)的聚類方法直接在隱私數(shù)據(jù)上執(zhí)行,但是無法應(yīng)對(duì)針對(duì)攻擊者任意背景知識(shí)的大規(guī)模數(shù)據(jù)挖掘任務(wù)中的惡意攻擊。這將導(dǎo)致侵犯?jìng)€(gè)人隱私,以及通過系統(tǒng)資源和聚類輸出造成泄漏。

    現(xiàn)階段,在非交互式場(chǎng)景下,針對(duì)差分隱私數(shù)據(jù)發(fā)布基本方法是基于直方圖發(fā)布。這種方法可以很直觀地看出數(shù)據(jù)的分布并可以近似計(jì)算數(shù)據(jù)和、差、方差以及平均數(shù)等統(tǒng)計(jì)信息。然而,當(dāng)數(shù)據(jù)中的屬性數(shù)量增加時(shí),這種方法則會(huì)表現(xiàn)出嚴(yán)重的局限性,現(xiàn)在我們可以得到的數(shù)據(jù)量增大,其數(shù)據(jù)類型也在不斷的增加,我們?nèi)绻胍尤轿坏氖占畔?,信息中的屬性相?yīng)的也會(huì)增加。這使得計(jì)算效率大大降低。除此之外,直方圖發(fā)布方法僅僅提供了分區(qū)數(shù)據(jù)的計(jì)數(shù),比如年齡為20~30的旅客有2000個(gè),卻無法提供具體數(shù)據(jù),這使得我們僅僅可以看到經(jīng)過統(tǒng)計(jì)處理的數(shù)據(jù),數(shù)據(jù)的可用性降低。本文提出滿足差分隱私的數(shù)據(jù)集發(fā)布,這將克服數(shù)據(jù)屬性數(shù)量及無具體數(shù)據(jù)值的限制。為進(jìn)一步提高數(shù)據(jù)集的可用性,本文提出利用聚類算法先對(duì)數(shù)據(jù)進(jìn)行聚類,將一個(gè)大數(shù)據(jù)集聚類為幾個(gè)子集,查詢敏感度由一分化為多,使整個(gè)數(shù)據(jù)集的查詢敏感度降低,進(jìn)而減少噪聲添加量。

    針對(duì)于此,本文提出了一種基于聚類滿足差分隱私的數(shù)據(jù)發(fā)布算法,可以對(duì)民航旅客信息數(shù)據(jù)集進(jìn)行差分隱私保護(hù)。分析民航旅客數(shù)據(jù),從中發(fā)現(xiàn)數(shù)據(jù)包括數(shù)值屬性和分類屬性,k-prototype可對(duì)混合數(shù)據(jù)集聚類,因此使用k-prototype聚類算法,該算法是處理混合數(shù)據(jù)類型的典型聚類算法。根據(jù)民航旅客數(shù)據(jù)集特有的性質(zhì)改進(jìn)了k-prototype算法中分類屬性的距離計(jì)算,能夠更好的對(duì)混合屬性數(shù)據(jù)進(jìn)行聚類,最后通過聚類算法的分組,由單一數(shù)據(jù)集分成多組數(shù)據(jù)集,降低差分隱私的敏感度,從而減少需要添加的噪聲量,因此大大提高數(shù)據(jù)的可用性。

    2 差分隱私定義及噪音機(jī)制

    2.1 差分隱私定義

    差分隱私要求數(shù)據(jù)的輸出應(yīng)該大致相同,即使在輸入數(shù)據(jù)庫中任意添加或刪除任何一個(gè)元組也是如此。

    定義1ε-差分隱私:設(shè)有隨機(jī)算法R滿足ε-差分隱私,以及任意兩個(gè)相鄰數(shù)據(jù)集D1和D2(在一條記錄中有所不同)輸出S∈Range(R)。 有

    (1)

    式中:參數(shù)ε是隱私預(yù)算,由公式可得出ε越小,隨機(jī)算法R隱私保護(hù)程度越高。

    差分隱私提供了正式且可量化的隱私保證,而不受對(duì)手的背景知識(shí)和可用的計(jì)算能力的影響。如果對(duì)于整個(gè)輸出空間,對(duì)于任何一對(duì)相鄰輸入,生成相同輸出的概率彼此之間的較小倍數(shù)內(nèi),則認(rèn)為隨機(jī)算法是差分隱私的。這意味著對(duì)于彼此接近的任何兩個(gè)數(shù)據(jù)集,差分隱私算法在兩個(gè)數(shù)據(jù)集上的表現(xiàn)大致相同。無論對(duì)手是否擁有先驗(yàn)知識(shí),該概念都為用戶提供了足夠的隱私保護(hù)。在本文中,當(dāng)且僅當(dāng)數(shù)據(jù)集D1和D2僅相差一條記錄時(shí),我們才將兩個(gè)數(shù)據(jù)集D1和D2視為鄰居,D1+t表示將元組t與數(shù)據(jù)集D1相加而得的數(shù)據(jù)集,我們用D1?D2來表示。這可以保護(hù)任何單個(gè)元組的隱私。

    定義2 敏感度:對(duì)于輸入數(shù)據(jù)集上的任何查詢f,對(duì)于任何兩個(gè)相鄰的數(shù)據(jù)集D1和D2,f的敏感度為

    (2)

    敏感度較低的查詢所需的噪聲較小。差分隱私保護(hù)技術(shù)即向查詢結(jié)果的數(shù)據(jù)添加噪聲,對(duì)數(shù)據(jù)擾動(dòng),最終達(dá)到隱私保護(hù)的目的。

    對(duì)于滿足差分隱私的數(shù)據(jù)發(fā)布來說,需要多次應(yīng)用差分隱私算法,差分隱私的兩個(gè)組合性質(zhì)可以將隱私預(yù)算合理的分配到整個(gè)差分隱私中。

    性質(zhì)1 序列組合性:設(shè)存在n個(gè)隨機(jī)算法K,數(shù)據(jù)集D, {ki(D)|1≤i≤n} 滿足εi-差分隱私,則K在D上整體滿足∑εi-差分隱私。

    性質(zhì)2 并行組合性:設(shè)存在n個(gè)隨機(jī)算法K,作用在互不相交的一組數(shù)據(jù)集 {Di|1≤i≤n}, {ki(Di)|1≤i≤n} 滿足εi-差分隱私,則K在D上整體滿足max(εi)-差分隱私。

    2.2 噪音機(jī)制

    噪聲量會(huì)影響數(shù)據(jù)安全性和可用性,Laplace機(jī)制與指數(shù)機(jī)制為現(xiàn)有常用的噪音機(jī)制。拉普拉斯機(jī)制對(duì)數(shù)值型數(shù)據(jù)進(jìn)行處理,是實(shí)現(xiàn)差分隱私最常用的機(jī)制。本文使用拉普拉斯機(jī)制和指數(shù)機(jī)制滿足差分隱私。

    定理1Laplace機(jī)制:對(duì)于任意數(shù)據(jù)集D和函數(shù)f,若滿足

    (3)

    則算法K滿足ε-差分隱私保護(hù)。其中,Δf為全局敏感度,lap(Δf/ε) 代表添加的噪聲量,根據(jù)以上公式可知,查詢函數(shù)f的全局敏感度越大,所需添加的噪聲量越大。

    定理2 指數(shù)機(jī)制:給定一個(gè)可用性函數(shù)u(D,r), Δu為u(D,r) 的全局敏感度,O表示輸出域,r表示從O中所選取一個(gè)元素,若算法K滿足

    (4)

    則K滿足ε-差分隱私。

    指數(shù)機(jī)制的關(guān)鍵是可用性函數(shù)u(D,r),u評(píng)估輸出值r,指數(shù)機(jī)制是用指數(shù)級(jí)更高的概率選擇評(píng)分更高的元素,u(D,r) 越大,r被輸出的概率越大。

    3 基于k-prototype聚類的差分隱私民航旅客數(shù)據(jù)發(fā)布算法

    3.1 改進(jìn)k-prototype算法

    k-prototype是聚類中常用的模型,該算法從代表k個(gè)組的k個(gè)隨機(jī)選擇的點(diǎn)開始,計(jì)算每個(gè)記錄到初始點(diǎn)的距離,然后將樣本迭代地聚類到最近的簇,并通過聚類到這些點(diǎn)的樣本的眾值來更新這些點(diǎn)。直至所有的記錄全部分到最近的簇。本文根據(jù)民航旅客數(shù)據(jù)的特點(diǎn),改進(jìn)k-prototype聚類算法。首先確定最小k值和最大k值,依次對(duì)每個(gè)k值進(jìn)行實(shí)驗(yàn),隨機(jī)選擇k個(gè)初始類中心點(diǎn),根據(jù)元組間距離計(jì)算方法得到一條記錄到每個(gè)中心點(diǎn)距離,迭代對(duì)數(shù)據(jù)集進(jìn)行計(jì)算,選取距離最小的簇,得到初步聚類結(jié)果。進(jìn)而選取每個(gè)簇中出現(xiàn)最多的元組作類中心,進(jìn)行距離計(jì)算,再次選取距離最小的簇,迭代至每個(gè)元組跟上一次的簇?cái)?shù)一樣,然后計(jì)算聚類有效值,選擇具有最小有效值評(píng)價(jià)的k值,輸出最佳聚類結(jié)果。

    (5)

    分類型屬性通過漢明距離計(jì)算屬性差異度

    (6)

    針對(duì)分類型數(shù)據(jù)距離計(jì)算結(jié)果對(duì)聚類結(jié)果的影響程度,設(shè)置其在算法中的權(quán)重,以達(dá)到提高聚類精度的目的。綜上所述,最終元組至中心點(diǎn)的距離計(jì)算公式為上述兩個(gè)公式(式(5)、式(6))相加。k-prototype聚類中的損失函數(shù)計(jì)算數(shù)值型和分類性到每簇聚類中心點(diǎn)的距離,因此選擇一個(gè)合適的損失函數(shù)y。假定設(shè)定數(shù)據(jù)集的簇為k個(gè),使用0和1表示第m(1≤m≤k) 個(gè)聚類中是否存在元組i,0則表示存在,反之,1表示不存在。綜上損失函數(shù)可以定義為

    (7)

    (8)

    3.2 算法及差分隱私證明

    將訓(xùn)練集隨機(jī)分為多個(gè)子集,對(duì)每個(gè)子集運(yùn)行聚類算法以獲取輸出,然后使用差分隱私對(duì)每個(gè)元組加噪,具體做法為遍歷每一個(gè)元組并確定其所在簇,形成不同簇之后,總結(jié)每個(gè)簇的分類屬性中的屬性值集合,利用屬性值集合使用指數(shù)機(jī)制對(duì)分類型屬性干擾,數(shù)值型采用拉普拉斯機(jī)制方法加噪,生成滿足差分隱私的數(shù)據(jù)集。此方法將查詢函數(shù)的靈敏度由一整個(gè)數(shù)據(jù)集的所有記錄分散至多個(gè)子集的k個(gè)記錄中,一定程度上減少了數(shù)據(jù)中的噪音量,更接近于原始數(shù)據(jù),加噪后的數(shù)據(jù)可用性大大提高。具體算法見表1。

    算法流程如圖1所示。

    步驟:

    步驟1 選擇民航數(shù)據(jù)中敏感信息屬性,數(shù)據(jù)清洗,并分析屬性之間相關(guān)性。

    步驟2 選擇k最小值以及最大值。

    步驟3 隨機(jī)選擇k個(gè)初始中心點(diǎn),計(jì)算每個(gè)元組至初始中心點(diǎn)距離,元組中的數(shù)值型屬性計(jì)算方法為每個(gè)點(diǎn)到k個(gè)中心點(diǎn)的歐氏距離,其中的分類型屬性計(jì)算方法為每個(gè)點(diǎn)

    表1 基于聚類的差分隱私數(shù)據(jù)發(fā)布算法

    圖1 算法流程

    到k個(gè)中心點(diǎn)的漢明距離,最后將其分至最近的中心點(diǎn)。

    步驟4 計(jì)算損失函數(shù),若結(jié)果不為0,則重新選擇中心點(diǎn),迭代至損失函數(shù)結(jié)果為0。

    步驟5 將聚類結(jié)果進(jìn)行有效性評(píng)估,選擇聚類效果最好的k值。得到聚類結(jié)果。

    步驟6 對(duì)于數(shù)值型屬性,采用Laplace機(jī)制加噪。

    步驟7 對(duì)于分類型屬性,生成每個(gè)簇的屬性值集合,每個(gè)簇分別根據(jù)各自屬性值使用指數(shù)機(jī)制選擇輸出值。

    步驟8 生成滿足差分隱私的數(shù)據(jù)集。

    本文算法分為聚類分組階段和數(shù)據(jù)發(fā)布階段,簇的數(shù)量從k在一定范圍,這導(dǎo)致了一系列連續(xù)的聚類結(jié)果。具體而言,在每個(gè)循環(huán)中,該方法的基本步驟包括:①使用改進(jìn)的k-prototype算法和距離度量,將輸入數(shù)據(jù)集劃分為所需的聚類;②根據(jù)聚類評(píng)估聚類結(jié)果有效性指數(shù);③最后,繪制了給定數(shù)據(jù)的聚類有效性指數(shù)與聚類數(shù)量的關(guān)系圖。根據(jù)該圖,可以目測(cè)為給定的混合數(shù)據(jù)集提供最佳的聚類數(shù)。民航旅客數(shù)據(jù)集都是混合數(shù)據(jù)集,它包含數(shù)值屬性和分類屬性。如果要以統(tǒng)一的方式處理混合數(shù)據(jù),一般是將分類屬性轉(zhuǎn)換為數(shù)值屬性,或者將數(shù)字屬性轉(zhuǎn)換為分類屬性。但是,將分類屬性轉(zhuǎn)換為數(shù)值屬性,很難將合適的數(shù)值分配給分類值。例如,如果color屬性采用集合{red,blue,green}中的值,將該集合轉(zhuǎn)換為{1、2、3}。在這種情況下,計(jì)算任何編碼值之間的距離都是不合適的。將數(shù)值屬性轉(zhuǎn)換為分類屬性,需要使用離散算法將實(shí)值變量的值域劃分為幾個(gè)區(qū)間,并為同一區(qū)間中的所有值分配一個(gè)符號(hào)。但是由于沒有考慮這些值對(duì)離散值的隸屬度,所以通常會(huì)導(dǎo)致信息丟失。因此,本文直接對(duì)混合數(shù)據(jù)進(jìn)行聚類。完成聚類后,進(jìn)行第二個(gè)階段,數(shù)據(jù)中的數(shù)值型采用Laplace機(jī)制以一定概率添加噪聲,通過指數(shù)機(jī)制對(duì)分類型數(shù)據(jù)實(shí)現(xiàn)差分隱私,從一個(gè)簇中獲取屬性值,以簇中每個(gè)屬性值出現(xiàn)的概率選擇值,根據(jù)輸入數(shù)據(jù),差分隱私參數(shù)和概率選擇值,選擇隨出現(xiàn)概率呈指數(shù)增長(zhǎng)。

    定理3 DP-k-prototype算法滿足ε-差分隱私。

    根據(jù)差分隱私組合性質(zhì)分析并證明:

    (9)

    (10)

    4 實(shí)驗(yàn)結(jié)果及分析

    本文選取民航數(shù)據(jù)中的PNR數(shù)據(jù)和離港數(shù)據(jù)。在其中選定了15個(gè)屬性,包括身份、航班、位置、時(shí)間等信息組成實(shí)驗(yàn)所需要的數(shù)據(jù)集(表2),并考慮異構(gòu)屬性類型,由48 842個(gè)記錄組成。對(duì)數(shù)據(jù)集進(jìn)行數(shù)據(jù)清洗,刪除數(shù)據(jù)集中包含空屬性以及有明顯錯(cuò)誤(如不合常規(guī)的年齡以及不存在的日期或者機(jī)場(chǎng)三字代碼等等)的整條記錄,對(duì)所有完整記錄規(guī)范其中的數(shù)據(jù),如上清洗后共有30 160個(gè)數(shù)據(jù)記錄。本節(jié)仿真實(shí)驗(yàn)環(huán)境如下:Intel(R) Core(TM)i5-4590 CPU,4 GB內(nèi)存,Windows8操作系統(tǒng),在pycharm環(huán)境下進(jìn)行仿真實(shí)驗(yàn)。

    表2 數(shù)據(jù)集介紹

    為方便實(shí)驗(yàn),選取其中400條記錄實(shí)現(xiàn)聚類,以及滿足差分隱私的算法。每個(gè)屬性之間的相關(guān)性越低,聚類效果越好,以票號(hào)為主屬性,計(jì)算其與另外14個(gè)屬性的相關(guān)性,如圖2所示。

    圖2 屬性之間的相關(guān)性

    從圖2中可看到只有兩個(gè)屬性之間的相關(guān)度達(dá)到了0.6,其余各個(gè)屬性之間的相關(guān)度大都低于0.1,由此可見,屬性之間的相關(guān)性差,依賴低,可以較好地聚類。

    做好以上數(shù)據(jù)預(yù)處理工作后,進(jìn)行聚類實(shí)驗(yàn)。

    4.1 實(shí)驗(yàn)過程

    實(shí)驗(yàn)時(shí)的參數(shù)設(shè)置,其中聚類實(shí)驗(yàn)設(shè)置k最小值為2,最大值為9,因?yàn)槿绻鹝值太大,則會(huì)造成聚類結(jié)果的實(shí)際意義不大,由于實(shí)驗(yàn)數(shù)據(jù)集總共有400條記錄,所以將聚類簇?cái)?shù)控制在10以內(nèi),實(shí)驗(yàn)結(jié)果如圖3所示。

    圖3 不同k值下的聚類效果(指標(biāo)越低越好)

    圖3顯示了k值與聚類指標(biāo)的關(guān)系,有效性指標(biāo)越小,則聚類效果越好,可以看出本文改進(jìn)的聚類算法比開源包的聚類算法效果要好。觀察改進(jìn)的聚類算法,可以明顯看出k=3時(shí),達(dá)到最低點(diǎn),聚類效果最好。所以選取k=3,將數(shù)據(jù)集聚類,并按照聚類結(jié)果對(duì)數(shù)據(jù)集中的每個(gè)簇分別進(jìn)行加噪。其中數(shù)值型數(shù)值進(jìn)行拉普拉斯機(jī)制加噪,分類型數(shù)值進(jìn)行指數(shù)機(jī)制加噪。

    選擇隱私預(yù)算ε為0.1時(shí)的滿足差分隱私的數(shù)據(jù)集與原數(shù)據(jù)集及進(jìn)行對(duì)比,表3展示了部分記錄中的4個(gè)屬性之間的對(duì)比。

    4.2 實(shí)驗(yàn)評(píng)價(jià)指標(biāo)

    數(shù)據(jù)可用性是通過采用差分隱私技術(shù)所產(chǎn)生的信息損失來測(cè)量。信息損失可以通過誤差平方和(SSE)量化,計(jì)算方法為滿足差分隱私數(shù)據(jù)集的記錄與其相對(duì)應(yīng)原始數(shù)據(jù)集中的記錄之間距離的平方和,誤差平方和(SSE)的計(jì)算公式如下

    (11)

    計(jì)算數(shù)據(jù)度量時(shí),對(duì)于數(shù)據(jù)中的數(shù)值使用標(biāo)準(zhǔn)歐幾里德距離;對(duì)于字符采用漢明距離。數(shù)據(jù)隱私受到保護(hù)的標(biāo)準(zhǔn)通過信息披露來衡量。信息披露的計(jì)算方法在本文中采用計(jì)算滿足差分隱私數(shù)據(jù)集正確匹配的原數(shù)據(jù)記錄的百分比,即記錄關(guān)聯(lián)(RL)

    表3 加噪前后數(shù)據(jù)集之間對(duì)比

    (12)

    式中:n表示原數(shù)據(jù)集的總記錄數(shù),差分隱私干擾后的數(shù)據(jù)集記錄t′的RL概率Pr(t′) 為

    (13)

    將本文提出的算法在民航旅客信息數(shù)據(jù)集上進(jìn)行實(shí)驗(yàn)。確定聚類算法中的k值為3,將隱私參數(shù)ε設(shè)置為 {0.001,0.002,0.003,0.004,0.005,0.006,0.007,0.008,0.009,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,2,3,4,5,6,7,8,9,10}, 根據(jù)以上數(shù)據(jù)進(jìn)行數(shù)據(jù)信息損失與信息泄露對(duì)比實(shí)驗(yàn),信息損失SSE與敏感度的關(guān)系結(jié)果如圖4所示。

    圖4 ε不同時(shí)的信息損失

    當(dāng)ε為0.001到0.01時(shí),SSE較大,即使使用聚類算法分簇后再加噪,信息損失仍較高,數(shù)據(jù)可用性低,當(dāng)ε大于0.01時(shí),SSE較小,并逐漸減小后更趨于保持穩(wěn)定,信息損失較低,可用性高。信息披露RL對(duì)比分析如圖5所示。

    圖5 ε不同時(shí)的信息披露

    當(dāng)ε為0.001至0.1時(shí),RL值的變化幅度較小,且RL值趨于低水平;當(dāng)ε為0.1至1時(shí),RL值的變化幅度較大且處于上升的趨勢(shì);當(dāng)ε為1至10時(shí),RL值的變化幅度最大且近乎呈直線上升。圖中的這些變化是因?yàn)楫?dāng)ε越大時(shí),數(shù)據(jù)中添加的噪聲就越少,當(dāng)數(shù)據(jù)中添加的噪聲不能對(duì)數(shù)據(jù)集進(jìn)行保護(hù)時(shí),將會(huì)出現(xiàn)信息披露的情況。噪聲越小,風(fēng)險(xiǎn)就越高,導(dǎo)致隱私保護(hù)能力就越弱。圖4和圖5結(jié)合可以看出,當(dāng)ε=0.1至10時(shí),其數(shù)據(jù)可用性差距不大,而信息披露風(fēng)險(xiǎn)會(huì)越來越高,相應(yīng)隱私保護(hù)能力越來越低;當(dāng)ε小于0.1時(shí),其數(shù)據(jù)可用性越來越低,而信息披露風(fēng)險(xiǎn)則差距不大。因此當(dāng)ε=0.1獲得最佳效果。

    5 結(jié)束語

    本文提出了一種基于聚類滿足差分隱私的數(shù)據(jù)發(fā)布方法,目的是在數(shù)據(jù)發(fā)布做好隱私保護(hù)的同時(shí)最大限度提高數(shù)據(jù)的可用性,對(duì)于實(shí)驗(yàn)輸出的數(shù)據(jù)集滿足差分隱私保護(hù)模型要求已在本文做出完整的數(shù)學(xué)證明。對(duì)于此方法,我們改進(jìn)了k-prototype聚類算法,在原有k-prototype聚類的基礎(chǔ)上,針對(duì)特定的數(shù)據(jù)集特點(diǎn)(本文采用民航旅客信息數(shù)據(jù)集),總結(jié)出民航數(shù)據(jù)分為數(shù)值屬性(比如年齡、證件號(hào)等)和分類屬性(航班號(hào)、飛機(jī)場(chǎng)三字代碼等),對(duì)此采用不同的計(jì)算方法得出屬性差異度。將差異度較小的記錄分為一組,數(shù)據(jù)集經(jīng)過聚類從單一數(shù)據(jù)集變?yōu)槎嘟M數(shù)據(jù)集,且每組數(shù)據(jù)集中的數(shù)據(jù)都是相似的,對(duì)這些原始數(shù)據(jù)采用差分隱私,針對(duì)數(shù)值型使用Laplace機(jī)制,對(duì)數(shù)值進(jìn)行干擾,分類型使用指數(shù)機(jī)制,利用概率對(duì)分類屬性進(jìn)行干擾,最后對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行分析,該方法可以在信息披露較小的情況下盡量使信息損失量最小,從而在提供隱私保護(hù)的同時(shí)提高數(shù)據(jù)可用性。

    猜你喜歡
    信息
    訂閱信息
    中華手工(2017年2期)2017-06-06 23:00:31
    展會(huì)信息
    信息超市
    展會(huì)信息
    展會(huì)信息
    展會(huì)信息
    展會(huì)信息
    展會(huì)信息
    信息
    健康信息
    祝您健康(1987年3期)1987-12-30 09:52:32
    中文欧美无线码| 亚洲av一区综合| 91久久精品国产一区二区成人| 99久国产av精品国产电影| 国产精品美女特级片免费视频播放器| 国产人妻一区二区三区在| 身体一侧抽搐| 一区二区三区高清视频在线| av播播在线观看一区| 日本午夜av视频| 99热网站在线观看| 亚洲av福利一区| 一级a做视频免费观看| 丝袜美腿在线中文| 国产片特级美女逼逼视频| 韩国高清视频一区二区三区| 啦啦啦啦在线视频资源| 天堂网av新在线| 国产一区有黄有色的免费视频 | 亚洲图色成人| 尤物成人国产欧美一区二区三区| 在线 av 中文字幕| 在现免费观看毛片| 国产精品熟女久久久久浪| 99久国产av精品国产电影| 午夜激情欧美在线| 在线观看人妻少妇| 少妇的逼水好多| 国产一区二区三区综合在线观看 | 我的女老师完整版在线观看| 国产一区二区亚洲精品在线观看| 91精品伊人久久大香线蕉| 99热全是精品| 午夜福利成人在线免费观看| 秋霞伦理黄片| 久久久久久久久久人人人人人人| av网站免费在线观看视频 | 国产精品爽爽va在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 国产精品一区二区三区四区久久| 一边亲一边摸免费视频| 国产高清有码在线观看视频| 嫩草影院入口| 亚洲高清免费不卡视频| 色尼玛亚洲综合影院| www.av在线官网国产| 91aial.com中文字幕在线观看| 欧美性猛交╳xxx乱大交人| 狂野欧美白嫩少妇大欣赏| 可以在线观看毛片的网站| 身体一侧抽搐| 综合色av麻豆| 婷婷色av中文字幕| 高清视频免费观看一区二区 | 成人亚洲精品av一区二区| 搞女人的毛片| 久久精品久久久久久噜噜老黄| 日日摸夜夜添夜夜爱| 国产v大片淫在线免费观看| 日韩,欧美,国产一区二区三区| 黄片无遮挡物在线观看| 我要看日韩黄色一级片| 三级国产精品片| 国产亚洲精品av在线| 亚洲av中文av极速乱| 国产一区二区亚洲精品在线观看| 美女高潮的动态| 日日摸夜夜添夜夜添av毛片| 日本色播在线视频| 久久午夜福利片| 亚洲图色成人| 国产精品一区二区性色av| 亚洲18禁久久av| 久久久久网色| 80岁老熟妇乱子伦牲交| 69人妻影院| 国产精品蜜桃在线观看| av福利片在线观看| 国产国拍精品亚洲av在线观看| 久久久色成人| 天堂av国产一区二区熟女人妻| 国产爱豆传媒在线观看| 综合色丁香网| 亚洲精品乱码久久久v下载方式| 三级经典国产精品| 精品一区二区免费观看| 老司机影院成人| 久久久a久久爽久久v久久| 国产精品一区二区在线观看99 | 亚洲国产最新在线播放| 国产一级毛片在线| 亚洲国产日韩欧美精品在线观看| 国内少妇人妻偷人精品xxx网站| 久久久久久国产a免费观看| 精品人妻熟女av久视频| 97超碰精品成人国产| 久久99精品国语久久久| 日日干狠狠操夜夜爽| 精品久久久久久成人av| 亚洲第一区二区三区不卡| av又黄又爽大尺度在线免费看| 五月天丁香电影| 天堂av国产一区二区熟女人妻| 日日啪夜夜撸| 91精品伊人久久大香线蕉| 极品少妇高潮喷水抽搐| 成人毛片a级毛片在线播放| 夫妻午夜视频| 日本黄色片子视频| 黄色日韩在线| 2018国产大陆天天弄谢| 韩国高清视频一区二区三区| 大香蕉久久网| 国产一区亚洲一区在线观看| 青春草国产在线视频| 十八禁网站网址无遮挡 | 免费少妇av软件| 国产乱人视频| 成人午夜高清在线视频| 丝袜喷水一区| 日日啪夜夜撸| 欧美成人a在线观看| 精品人妻一区二区三区麻豆| 纵有疾风起免费观看全集完整版 | 午夜精品一区二区三区免费看| 永久免费av网站大全| 丰满乱子伦码专区| 男人狂女人下面高潮的视频| 成年女人看的毛片在线观看| 激情 狠狠 欧美| 蜜桃亚洲精品一区二区三区| av.在线天堂| 床上黄色一级片| 天堂俺去俺来也www色官网 | 国产精品女同一区二区软件| 免费看日本二区| 熟妇人妻不卡中文字幕| 你懂的网址亚洲精品在线观看| 一级片'在线观看视频| 午夜激情欧美在线| 老司机影院毛片| 久久久国产一区二区| 国产麻豆成人av免费视频| 日韩在线高清观看一区二区三区| 日韩精品有码人妻一区| 亚洲最大成人手机在线| 欧美日本视频| 亚洲国产av新网站| 国产精品久久久久久精品电影小说 | 久久99热这里只频精品6学生| 蜜桃久久精品国产亚洲av| 夜夜看夜夜爽夜夜摸| 亚洲国产欧美人成| 亚洲不卡免费看| 亚洲第一区二区三区不卡| 国产视频首页在线观看| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 久久久精品94久久精品| 久久久久久久大尺度免费视频| 精品不卡国产一区二区三区| 日韩强制内射视频| 午夜精品国产一区二区电影 | 只有这里有精品99| 国产在视频线精品| 久久久久免费精品人妻一区二区| 草草在线视频免费看| 毛片女人毛片| 中国美白少妇内射xxxbb| 日韩伦理黄色片| 内地一区二区视频在线| 亚洲精品自拍成人| videossex国产| 亚洲av中文字字幕乱码综合| 日韩不卡一区二区三区视频在线| 99热6这里只有精品| 伦精品一区二区三区| 高清日韩中文字幕在线| 韩国av在线不卡| 免费无遮挡裸体视频| 天堂√8在线中文| 免费看av在线观看网站| 国产午夜福利久久久久久| 国产亚洲最大av| 精品99又大又爽又粗少妇毛片| 国产午夜福利久久久久久| 美女xxoo啪啪120秒动态图| 亚洲美女搞黄在线观看| 午夜日本视频在线| 一个人观看的视频www高清免费观看| 国产成人a区在线观看| 2021天堂中文幕一二区在线观| 91在线精品国自产拍蜜月| 国产色爽女视频免费观看| 神马国产精品三级电影在线观看| 一区二区三区乱码不卡18| 小蜜桃在线观看免费完整版高清| 日韩在线高清观看一区二区三区| 久久久成人免费电影| 麻豆成人午夜福利视频| 网址你懂的国产日韩在线| 亚洲内射少妇av| 国产一级毛片在线| 亚洲第一区二区三区不卡| 一级毛片久久久久久久久女| 国产色婷婷99| 又大又黄又爽视频免费| 免费在线观看成人毛片| 成年女人看的毛片在线观看| 一个人观看的视频www高清免费观看| 欧美潮喷喷水| 久久久久久伊人网av| 少妇人妻精品综合一区二区| 欧美最新免费一区二区三区| 深爱激情五月婷婷| 欧美97在线视频| 97人妻精品一区二区三区麻豆| 五月天丁香电影| 免费大片黄手机在线观看| 亚洲av成人av| 亚洲天堂国产精品一区在线| 人人妻人人澡欧美一区二区| 日本一二三区视频观看| 亚洲色图av天堂| 黑人高潮一二区| 欧美丝袜亚洲另类| 国产精品一区二区性色av| 久久精品久久精品一区二区三区| 亚洲精品一二三| 99热网站在线观看| 午夜福利成人在线免费观看| 国产三级在线视频| 最近2019中文字幕mv第一页| 人体艺术视频欧美日本| 欧美一区二区亚洲| 中文字幕制服av| 国产精品一区二区在线观看99 | 亚洲国产日韩欧美精品在线观看| 久久99精品国语久久久| 国产av国产精品国产| 日本黄色片子视频| 97超碰精品成人国产| 神马国产精品三级电影在线观看| 国产黄色小视频在线观看| 又爽又黄无遮挡网站| 尾随美女入室| 我的老师免费观看完整版| 一级毛片电影观看| 人妻系列 视频| 97超碰精品成人国产| 男人舔女人下体高潮全视频| 性插视频无遮挡在线免费观看| 久久精品国产亚洲av天美| 国产白丝娇喘喷水9色精品| 国产 一区精品| 国产视频首页在线观看| 麻豆成人午夜福利视频| 久久这里有精品视频免费| 国产精品无大码| 欧美成人a在线观看| 免费看av在线观看网站| 成年版毛片免费区| 日韩制服骚丝袜av| 成人无遮挡网站| 日本一本二区三区精品| 色综合色国产| www.色视频.com| 免费观看无遮挡的男女| 观看免费一级毛片| av.在线天堂| 少妇裸体淫交视频免费看高清| 亚洲人成网站高清观看| 国产精品无大码| 老司机影院成人| 少妇高潮的动态图| 精品久久久久久久久av| 女的被弄到高潮叫床怎么办| 日本黄大片高清| 久久精品久久精品一区二区三区| 青春草视频在线免费观看| 亚洲精品成人久久久久久| 久久久色成人| 99久久九九国产精品国产免费| 国产亚洲午夜精品一区二区久久 | 日韩欧美国产在线观看| 成人高潮视频无遮挡免费网站| 日韩人妻高清精品专区| 日韩大片免费观看网站| 高清欧美精品videossex| 亚洲人成网站在线观看播放| 男的添女的下面高潮视频| 爱豆传媒免费全集在线观看| 亚洲精品影视一区二区三区av| 免费黄网站久久成人精品| 国产老妇女一区| 国产欧美另类精品又又久久亚洲欧美| 免费黄色在线免费观看| 色5月婷婷丁香| 男女下面进入的视频免费午夜| 一夜夜www| 天堂影院成人在线观看| 天堂中文最新版在线下载 | 中文天堂在线官网| 日韩强制内射视频| 色视频www国产| 夫妻性生交免费视频一级片| 高清日韩中文字幕在线| 成年人午夜在线观看视频 | 免费黄频网站在线观看国产| 插逼视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 老司机影院成人| 日韩欧美一区视频在线观看 | 午夜福利在线在线| 欧美性感艳星| 久久精品久久久久久久性| 岛国毛片在线播放| 欧美日韩精品成人综合77777| 亚洲精品久久午夜乱码| 干丝袜人妻中文字幕| 亚洲国产日韩欧美精品在线观看| 国产一区二区三区av在线| 成人漫画全彩无遮挡| 色视频www国产| 人妻一区二区av| 日韩制服骚丝袜av| 在线 av 中文字幕| 赤兔流量卡办理| 午夜福利在线在线| 国产黄色小视频在线观看| 亚洲精品456在线播放app| 日韩国内少妇激情av| 国产伦精品一区二区三区视频9| 久久6这里有精品| 国产亚洲精品久久久com| 免费少妇av软件| 神马国产精品三级电影在线观看| 国产亚洲午夜精品一区二区久久 | 大香蕉97超碰在线| 美女被艹到高潮喷水动态| 少妇高潮的动态图| 亚洲精品乱码久久久久久按摩| 97超视频在线观看视频| 成人午夜高清在线视频| 成年av动漫网址| 国产男女超爽视频在线观看| 97超视频在线观看视频| 国产人妻一区二区三区在| 小蜜桃在线观看免费完整版高清| 看非洲黑人一级黄片| 丝袜美腿在线中文| 一级爰片在线观看| 一区二区三区高清视频在线| 18禁在线播放成人免费| 啦啦啦啦在线视频资源| 老司机影院毛片| 国产国拍精品亚洲av在线观看| 欧美3d第一页| 欧美日本视频| 国产黄频视频在线观看| 欧美日韩视频高清一区二区三区二| 国产成人精品一,二区| 国产精品人妻久久久影院| 国产精品久久久久久av不卡| 亚洲欧美成人精品一区二区| 国产成人91sexporn| 美女高潮的动态| 国产精品一区二区性色av| 69av精品久久久久久| 亚洲在线自拍视频| av在线观看视频网站免费| 91精品伊人久久大香线蕉| 国产又色又爽无遮挡免| 草草在线视频免费看| 久久精品夜色国产| 波野结衣二区三区在线| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 国产白丝娇喘喷水9色精品| 五月玫瑰六月丁香| 成年免费大片在线观看| 国产精品一区二区在线观看99 | 69av精品久久久久久| 国产激情偷乱视频一区二区| 亚洲国产av新网站| 国产精品一二三区在线看| 欧美区成人在线视频| 夫妻午夜视频| 日本三级黄在线观看| av在线观看视频网站免费| 日韩亚洲欧美综合| 伦精品一区二区三区| 日本与韩国留学比较| 久久精品久久精品一区二区三区| av黄色大香蕉| 亚洲精品自拍成人| 综合色丁香网| 色尼玛亚洲综合影院| 亚洲精品视频女| 丝袜美腿在线中文| 欧美97在线视频| av播播在线观看一区| 亚洲国产最新在线播放| 啦啦啦韩国在线观看视频| 久久久午夜欧美精品| 欧美zozozo另类| 黄色欧美视频在线观看| 哪个播放器可以免费观看大片| 18禁裸乳无遮挡免费网站照片| 久久这里有精品视频免费| 亚洲熟女精品中文字幕| 久久久久九九精品影院| 国产激情偷乱视频一区二区| 国产午夜精品久久久久久一区二区三区| 久久6这里有精品| 精品人妻一区二区三区麻豆| 人妻夜夜爽99麻豆av| 最近视频中文字幕2019在线8| 久热久热在线精品观看| 亚洲成人av在线免费| 亚洲av福利一区| a级毛色黄片| 人体艺术视频欧美日本| 亚洲内射少妇av| 久久人人爽人人片av| 伊人久久精品亚洲午夜| 菩萨蛮人人尽说江南好唐韦庄| 国内精品宾馆在线| kizo精华| 精品国产三级普通话版| av在线播放精品| 国产 一区 欧美 日韩| 在线观看av片永久免费下载| 成人漫画全彩无遮挡| 丝瓜视频免费看黄片| 欧美三级亚洲精品| 激情 狠狠 欧美| 成年av动漫网址| 国产精品一区二区性色av| 欧美一区二区亚洲| 极品少妇高潮喷水抽搐| 在线观看av片永久免费下载| 一个人免费在线观看电影| 极品少妇高潮喷水抽搐| ponron亚洲| 免费观看无遮挡的男女| 国产淫语在线视频| 深夜a级毛片| 免费观看的影片在线观看| 欧美日本视频| 美女国产视频在线观看| 97精品久久久久久久久久精品| 一个人看的www免费观看视频| 天堂√8在线中文| 午夜精品国产一区二区电影 | 晚上一个人看的免费电影| 一级二级三级毛片免费看| 校园人妻丝袜中文字幕| 免费黄频网站在线观看国产| 久久99蜜桃精品久久| 最后的刺客免费高清国语| 亚洲欧美成人精品一区二区| 亚洲综合色惰| 午夜激情欧美在线| 国产精品一区二区在线观看99 | 免费在线观看成人毛片| 国产成人精品久久久久久| 少妇猛男粗大的猛烈进出视频 | 欧美日韩一区二区视频在线观看视频在线 | 少妇高潮的动态图| 久久久久久久午夜电影| 国产精品蜜桃在线观看| 天天躁夜夜躁狠狠久久av| 六月丁香七月| 国产亚洲精品av在线| 2018国产大陆天天弄谢| 亚洲国产精品国产精品| 亚洲精品,欧美精品| 日韩亚洲欧美综合| 国产精品一区二区三区四区免费观看| 亚洲av成人av| 久久精品国产亚洲av天美| 哪个播放器可以免费观看大片| 国产精品一区www在线观看| 国产精品久久久久久久电影| 日日干狠狠操夜夜爽| 国产午夜精品一二区理论片| 亚洲精品乱久久久久久| 欧美成人午夜免费资源| 国产真实伦视频高清在线观看| 午夜精品在线福利| 久久鲁丝午夜福利片| 最近最新中文字幕免费大全7| 91精品国产九色| 亚洲av国产av综合av卡| 春色校园在线视频观看| 人人妻人人澡人人爽人人夜夜 | 亚洲成人久久爱视频| 九九在线视频观看精品| 精品久久久久久久末码| 亚洲色图av天堂| 日本爱情动作片www.在线观看| 精品久久久久久久人妻蜜臀av| 美女国产视频在线观看| 成人鲁丝片一二三区免费| 在线观看av片永久免费下载| 午夜免费激情av| 少妇熟女aⅴ在线视频| 久久精品综合一区二区三区| 夜夜爽夜夜爽视频| av又黄又爽大尺度在线免费看| 国产亚洲av片在线观看秒播厂 | av免费在线看不卡| 九九爱精品视频在线观看| 人人妻人人澡欧美一区二区| 亚洲精品第二区| 亚洲精品aⅴ在线观看| 久久久久九九精品影院| 一个人看的www免费观看视频| 日日干狠狠操夜夜爽| 国产伦在线观看视频一区| 国产一区二区亚洲精品在线观看| 在线观看人妻少妇| 日韩不卡一区二区三区视频在线| 三级毛片av免费| 国产免费一级a男人的天堂| 亚洲婷婷狠狠爱综合网| 女人久久www免费人成看片| 91在线精品国自产拍蜜月| 日韩在线高清观看一区二区三区| 国产国拍精品亚洲av在线观看| 国产日韩欧美在线精品| 亚洲国产精品专区欧美| 久久久久久久国产电影| 婷婷色综合大香蕉| 午夜精品国产一区二区电影 | a级一级毛片免费在线观看| 精品久久久久久久久av| 日韩电影二区| 亚洲内射少妇av| 成人美女网站在线观看视频| 亚洲美女视频黄频| 蜜桃久久精品国产亚洲av| 日韩中字成人| 亚洲一区高清亚洲精品| 黄片无遮挡物在线观看| 亚洲精品成人av观看孕妇| 久久精品人妻少妇| 亚洲欧美一区二区三区国产| 国产女主播在线喷水免费视频网站 | 色综合色国产| 免费看美女性在线毛片视频| 免费观看精品视频网站| 免费电影在线观看免费观看| 亚洲精品一区蜜桃| 欧美日韩精品成人综合77777| 男人狂女人下面高潮的视频| 亚洲自拍偷在线| 久久精品久久久久久久性| 久久久午夜欧美精品| 久久人人爽人人片av| 黄片无遮挡物在线观看| 久99久视频精品免费| 99视频精品全部免费 在线| 视频中文字幕在线观看| 中文字幕av在线有码专区| 草草在线视频免费看| 秋霞在线观看毛片| 国内少妇人妻偷人精品xxx网站| 高清欧美精品videossex| 国产av码专区亚洲av| 一级毛片我不卡| 国产精品一区二区在线观看99 | 亚洲av成人av| 亚洲久久久久久中文字幕| 成人毛片a级毛片在线播放| 亚洲国产最新在线播放| 国产精品一区二区在线观看99 | 免费看av在线观看网站| 久久99精品国语久久久| 51国产日韩欧美| 国产成人精品久久久久久| 插阴视频在线观看视频| 欧美日韩综合久久久久久| 91精品伊人久久大香线蕉| 国产成人精品福利久久| 少妇猛男粗大的猛烈进出视频 | 色5月婷婷丁香| 一级片'在线观看视频| 国产亚洲av片在线观看秒播厂 | 久久6这里有精品| 大香蕉久久网| 少妇猛男粗大的猛烈进出视频 | 成年版毛片免费区| 亚洲国产最新在线播放| 人妻夜夜爽99麻豆av| 国内精品宾馆在线| 午夜福利在线观看免费完整高清在| 三级毛片av免费| 亚洲精品,欧美精品| 特大巨黑吊av在线直播| 久久这里只有精品中国| 乱码一卡2卡4卡精品| 国产国拍精品亚洲av在线观看| 午夜免费男女啪啪视频观看| 中国国产av一级| 黄色一级大片看看| 午夜福利在线观看吧| 欧美xxⅹ黑人| 亚洲四区av| 天堂中文最新版在线下载 | 亚洲精品成人久久久久久| 午夜福利成人在线免费观看| 最新中文字幕久久久久| 精品亚洲乱码少妇综合久久|