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

    基于Matrix Profile的時間序列變長模體挖掘

    2021-05-27 06:51:34朱曉曉王繼民
    計算機(jī)與現(xiàn)代化 2021年5期
    關(guān)鍵詞:模體歐氏等價

    朱 旭,朱曉曉,王繼民

    (河海大學(xué)計算機(jī)與信息學(xué)院,江蘇 南京 211100)

    0 引 言

    作為模式發(fā)現(xiàn)和相似性搜索的交叉主題,模體挖掘最先由加利福尼亞大學(xué)河濱分校的Lin和Keogh等[1]在2002年首次提出,稱這樣的模式為“模體”。在相關(guān)文獻(xiàn)中,模體被定義為重復(fù)模式、頻繁出現(xiàn)的趨勢或近似和重復(fù)的序列、形狀、片段、子序列等。模體挖掘也可為分類的核心[2]、群集[3-4]、異常檢測[5-6]及規(guī)則發(fā)現(xiàn)[7-8]算法。時間序列模體挖掘技術(shù)在遺傳學(xué)、金融、工業(yè)等諸多領(lǐng)域也得到應(yīng)用。例如,金融領(lǐng)域的證券交易數(shù)據(jù)[9]、工業(yè)領(lǐng)域的用電數(shù)據(jù)[10]等。

    定長模體挖掘僅挖掘由用戶指定長度的模體。目前對于定長模體挖掘研究較多。文獻(xiàn)[1]提出的Brute-Force算法是第一個模體挖掘算法,算法采取枚舉的方式計算所有子序列之間的距離,然后利用最近鄰算法找到模體,它具有完全的覆蓋率和準(zhǔn)確率,被廣泛用作基準(zhǔn)算法以衡量其他算法的準(zhǔn)確性。Chiu等[11]結(jié)合SAX算法與概率思想提出了Random Projection算法。Mueen等[12]提出了MK算法,使用早棄技術(shù)對Brute-Force算法進(jìn)行改進(jìn),并提出通過選擇參考序列得到更緊密的下界距離,從而提高算法效率。Yeh等[13]引入STAMP算法,使用快速相似性搜索算法[14]挖掘時間序列模體。Zhu等[15]對STAMP算法進(jìn)行改進(jìn)提出了STOMP算法,其計算速度更快,通過快速傅里葉變換實現(xiàn)子序列間的點(diǎn)積來計算子序列之間的距離,通過重用前一個相鄰子序列的點(diǎn)積來加速當(dāng)前子序列點(diǎn)積的計算。相較于STAMP的隨機(jī)順序計算,STOMP是一種有序搜索。

    變長模體挖掘無需用戶預(yù)先指定模體長度,相比時間序列定長模體挖掘更難以解決,針對此問題的研究也相對較少。變長模體挖掘一般以定長模體挖掘為基礎(chǔ),通過在可行長度范圍內(nèi)迭代,得到不同長度的模體。Nunthanid等[16]提出了VLMD算法,利用MK算法計算所有可能長度的模體,然后通過模體分組和提取模體分組代表找到自適應(yīng)長度的模體。VLMD算法的時間復(fù)雜度較高,可擴(kuò)展性較差。Mueen等[7]提出了MOEN算法,采用下界距離加速Brute-Force算法的效率,以加速枚舉不同長度模體。Tang等[17]擴(kuò)展了K-Motif算法,基于RP算法,定義搜索空間,并使用分段重疊和分段等價類條件從模體分組中獲取變長模體。

    以上算法雖然在不同方面提高了變長模體的發(fā)現(xiàn)速度,但仍然存在參數(shù)過多、算法可擴(kuò)展性差以及時間復(fù)雜度高等問題。本文提出一種基于Matrix Profile[15]的時間序列變長模體挖掘算法稱為Matrix Profile Variable-Length Motif Discovery (MPVLMD),基于VLMD算法的計算框架,采用STOMP算法挖掘所有可能長度下的定長模體,使用結(jié)合增量距離的下界距離技術(shù)提升計算速度,引入模體重疊條件、長度相似性條件和模體分組等價類劃分,該算法能有效地發(fā)現(xiàn)變長模體?;赨CR數(shù)據(jù)集,針對MPVLMD的運(yùn)行效率和準(zhǔn)確性進(jìn)行實驗,實驗結(jié)果表明,本文提出的算法具有較高的計算效率和準(zhǔn)確性。

    1 相關(guān)定義

    定義1時間序列[18]。時間序列T是在時間上有序的值序列,可記為T={t1,t2,…,tn}。其中,ti表示一個觀察值,|T|=n表示時間序列T的長度。

    定義2時間序列子序列[19]。時間序列子序列是T中一個更小的時間序列,記為Ti,m={ti,ti+1,…,ti+m-1}。其中,i是子序列在T中的開始位置,m是子序列的長度,并且m>1。

    定義3時間序列模體[20]。時間序列模體Mw是指時間序列T中,長度為w且彼此相似度最高的一對子序列??蓪⑵涠x為一個四元組:Mw=(d,L1,L2,w),其中,L1和L2為2個子序列的起始位置,滿足L1

    定義4平凡匹配。給定時間序列T,假設(shè)存在子序列C和F,其起始位置分別為p和q。如果p=q,那么C和F屬于平凡匹配。

    定義5模體重疊[16]。2個長度分別為i和j的模體Mi和Mj,當(dāng)且僅當(dāng)它們滿足條件:Mi.L1≤Mj.L1

    定義6模體分組[16]。模體分組MG是模體的集合,且包含的2個連續(xù)模體互相重疊。

    定義7模體分組代表[16]。模體分組代表MR,指模體分組MG所包含的所有不同長度的模體中d值最小的模體。

    定義8模體分組重疊。2個模體分組MGi和MGj,假設(shè)任意Mw∈MGi,Mw∈MGj當(dāng)且僅當(dāng)它們滿足條件Mw.L1=Mx.L1‖Mw.L1=Mx.L2‖Mw.L2=Mx.L1‖Mw.L1=Mx.L2時,稱模體分組MGi和MGj重疊。

    定義9變長模體集合。給定時間序列T和模體長度范圍[lmin,…,lmax],變長模體集合是所有滿足模體長度約束,且彼此不重疊的模體的集合。

    定義10距離矩陣。距離矩陣是時間序列中所有的子序列之間歐氏距離的矩陣。距離矩陣的特征表示為:

    其中,di,j為時間序列中第i個子序列與第j個子序列之間的歐氏距離。

    定義11Matrix Profile[15]。Matrix Profile是時間序列中所有子序列與其最相似的子序列間距離的向量。即距離矩陣中每一列的最小值。向量表示為:(P1,P2,…,Pn-m+1),其中,n表示時間序列長度,m表示子序列長度,Pi表示第i個子序列與所有子序列間距離的最小值Min(Di)。

    定義12Matrix Profile Index[15]。Matrix Profile Index是時間序列中所有子序列與其最相似的子序列的索引向量。向量索引表示為:(I1,I2,…,In-m+1),其中Ii表示第i個子序列的最近的鄰居。

    2 STOMP算法

    2016年,Yeh等提出了基于Matrix Profile的時間序列模體挖掘算法STAMP[13],在其基礎(chǔ)上,Zhu等進(jìn)一步提出了STOMP算法[15],在距離矩陣基礎(chǔ)上提出了新的數(shù)據(jù)結(jié)構(gòu)Matrix Profile,該向量中的最小距離值對應(yīng)的是彼此相似度最高的2個子序列,即為該時間序列的模體。STOMP算法通過快速傅里葉變換計算距離矩陣,從而可以有效地計算給定時間序列的Matrix Profile和Matrix Profile Index[15]。算法包含提取子序列、子序列全連接和模體發(fā)現(xiàn)這3個步驟。

    2.1 提取子序列

    根據(jù)定義2,使用長度為m的滑動窗口對時間序列T={t1,t2,…,tn}進(jìn)行分段,截取所有長度為m的子序列,T1,m={t1,t2,…,tm},T2,m={t2,t3,…,tm+1},…,Tn-m+1,m={tn-m+1,tn-m+2,…,tn},獲取所有子序列Ti,m={ti,ti+1,…,ti+m-1}(1≤i≤n-m+1)。

    2.2 子序列全連接

    子序列全連接用來計算時間序列的距離矩陣。計算所有子序列的均值μ和標(biāo)準(zhǔn)差σ,通過快速傅里葉變換計算子序列Ti,m與所有子序列之間的點(diǎn)積,使用z歸一化歐氏距離公式(1)計算2個子序列間的距離。

    (1)

    其中,m是子序列的長度,μi和μj分別是Ti,m和Tj,m的均值,σi和σj分別是Ti,m和Tj,m的標(biāo)準(zhǔn)差,QTi,j是Ti,m和Tj,m的點(diǎn)積。計算di,j所需的時間只取決于QTi,j計算所需的時間,QTi-1,j-1可以分解為:

    (2)

    并且QTi,j可以分解為:

    (3)

    由式(2)和式(3)可獲得式(4):

    QTi,j=QTi-1,j-1-ti-1tj-1+ti+m-1tj+m-1

    (4)

    當(dāng)計算第i個子序列與T中所有子序列的點(diǎn)積時,可以利用第i-1個子序列與T中所有子序列的點(diǎn)積結(jié)果加速計算。

    求得第i個子序列與T中所有子序列點(diǎn)積之后,使用z歸一化歐氏距離計算子序列之間的距離。

    2.3 模體發(fā)現(xiàn)

    根據(jù)定義10,使用Matrix Profile存儲距離矩陣中每列元素的最小值。同時根據(jù)定義11,將最小值的行號作為索引存入Matrix Profile Index。提取Matrix Profile中的最小值對應(yīng)的索引所指向的子序列是時間序列T的模體。其中,Matrix Profile的結(jié)構(gòu)如圖1所示。

    圖1 Matrix Profile結(jié)構(gòu)圖

    3 基于Matrix Profile的時間序列變長模體挖掘算法MPVLMD

    本文算法基于VLMD算法的計算框架,首先,使用STOMP算法作為子程序,并引入結(jié)合增量計算的下界距離加速策略,提升模體提取速度;其次,加入模體重疊和長度相似性條件進(jìn)行模體分組,踢除過短和存在平凡匹配的模體;再次,加入模體分組重疊條件對模體分組進(jìn)行等價類劃分,剔除提取模體代表時存在的過長模體;最后,提取每個分組等價類中的模體代表,模體代表集合即為變長模體?;贛atrix Profile的時間序列變長模體挖掘算法的流程如圖2所示。

    圖2 MPVLMD算法流程圖

    3.1 模體提取

    STOMP算法只能提取固定長度的模體,本文算法以STOMP作為子程序,在所有可能長度范圍內(nèi)迭代該程序,結(jié)合增量計算的下界距離加速策略,加速模體提取。

    3.1.1 求解z值數(shù)組

    依據(jù)下界距離的計算流程,首先需要計算對應(yīng)模體長度的z值。在提取所有可能長度模體的時候,模體長度范圍默認(rèn)為從2~n/2(其中n為時間序列長度),則需要計算長度為n/2-1的z值數(shù)組,其具體計算流程如算法1所示。

    算法1求解z值數(shù)組

    forj←m0+1 tomxdo

    fori←1 ton-jdo

    Maxj=Y

    return Max

    Output:z值數(shù)組Max//用于計算下界距離

    3.1.2 模體提取

    1)提取模長為m0的模體。

    使用STOMP算法提取模長為m0的模體,生成按距離升序排列的彼此最相似的子序列對列表List。

    2)結(jié)合增量計算的下界距離加速策略提取模長為m0+1的模體。

    使用緩存技術(shù)[21]保存時間序列值的累積和與累積平方和,便于下界距離和增量距離計算。

    3)提取所有可能長度的模體。

    重復(fù)上述步驟,直到獲取所有可能長度的模體。

    下面闡述下界距離和增量距離的具體求解過程。

    通過式(5)求模體長度為m0+1的下界距離:

    (5)

    其中,z=maxi(ti+m0-μi,m0)/σi,m0,(1≤i≤n-m0,當(dāng)模長為m0時,時間序列中第i+m0個元素歸一化后的最大值),d表示模長為m0的子序列間的z歸一化歐氏距離集合List中的最大值。

    (6)

    基于長度為m的tx、ty、(tx)2、(ty)2和txty的運(yùn)行和,就可以增量計算長度為m+1的距離函數(shù)的值。

    算法2模體提取

    Input:時間序列T,z值數(shù)組Max

    List←STOMP(T,m0)//提取長度為m0的模體,以及n-m0+1對相似子序列集合List

    z←Maxm0+1

    forj←m0+1 tomxdo

    forp←1 to List.Count do

    i←List(P).i;k←List(P).k;

    distance(Ti,j,Tk,j)←List(P).d+

    Sum(Ti,j-1,Tk,j-1,(Ti,j-1)2,(Tk,j-1)2,(Ti,j-1,Tk,j-1));

    d←distance(Ti,j,Tk,j);

    NewList.Add (d,i,k)//將距離d,位置i、k存入NewList中

    best←NewList.Min

    z←Maxj

    L1j←i,L2j←k

    output List.Min for lengthj

    else

    List←STOMP(T,j)

    Output:模體候選集合List

    3.2 模體分組

    模體提取階段產(chǎn)生的模體候選集合List中可能存在較長模體覆蓋較短模體的問題,因此對List中的模體進(jìn)行分組。

    遍歷List集合中的模體,按照定義5,將滿足模體重疊條件的2個模體置入相同模體分組中,反之創(chuàng)建新的模體分組,并將其中未分組的一個模體作為首個元素存儲到其中。

    遍歷所有的模體分組,并對同一個分組中的模體,使用長度相似性條件HMw=?n/2w」,修剪過短模體。如果模體Mw的HMw值與其他模體的HMother值不同,剔除模體Mw。其中,n表示時間序列長度,w表示模體長度。

    算法3模體分組

    Input:模體候選集合List

    for eachMi∈List do

    for eachMj∈List,j>i

    ifMioverlaps withMj

    Mj.groupid=Mi.groupid

    for each groupID groupid do

    Put all motifs that have the groupidiinto groupi

    for eachMw∈groupido

    if HMwnot equals to HMother

    deleteMw

    return group

    Output:模體分組集合group

    3.3 模體分組等價類劃分

    完成模體分組后,模體分組中的模體可能存在以下情況:(si,ti)∈groupi,(sj,tj)∈groupj,并且ti=tj。由于si與ti相似,sj與tj相似,則si與sj有很大概率是相似的。對模體分組集合group進(jìn)行模體分組重疊判斷,遍歷group中的每一個模體分組,將滿足定義8的模體分組置入同一個模體分組等價類中。

    算法4模體分組等價類劃分

    Input:模體分組集合group

    for each groupi∈group do

    for each groupj∈group,j>i

    if ?Mx∈groupioverlaps with?My∈groupj

    //如果Mx和My滿足模體分組重疊條件

    groupj.equalclassid=groupi.equalclassid

    return equalclass

    Output:模體分組等價類集合equalclass

    3.4 變長模體提取

    完成模體分組等價類劃分后,每個模體分組等價類將包含平凡匹配的模體。

    1)提取等價類模體代表。

    每個模體分組等價類中有多個模體分組,每個模體分組中有多個模體。提取模體分組等價類中每個模體分組中z歸一化歐氏距離最小的模體,作為每個模體分組的模體代表,每個模體等價類將包含多條模體代表。將這些模體分組的模體代表按照z歸一化歐氏距離正序排列,選擇中間位置模體代表的z歸一化歐氏距離(如果模體分組個數(shù)為奇數(shù)即為中間位置模體代表的z歸一化歐氏距離,如果是偶數(shù)取中間2個模體代表的z歸一化歐氏距離的均值)作為距離最大值,刪除z歸一化歐氏距離大于該最大距離的模體代表(修剪具有較大z歸一化歐氏距離的過長模體)。

    2)輸出變長模體。

    經(jīng)過模體分組,模體等價類劃分等操作剔除過長、過短和平凡匹配模體后,輸出每個模體分組等價類中長度最長的模體代表的集合,所輸出的所有不同長度模體代表的集合即為時間序列的變長模體。

    算法5變長模體提取

    Input:模體分組等價類集合equalclass

    for each equalclassid∈equalclass do

    for each groupi∈equalclassid do

    for eachMw∈groupido//提取每個分組中每個等價類中的模體代表

    //按照z歸一化歐氏距離正序排列

    //提取中間位置的模體作為模體分組等價類的模體代表

    return VMotifList

    Output:變長模體集合VMotifList

    4 實驗與結(jié)果分析

    4.1 實驗數(shù)據(jù)

    以UCR[22]的部分?jǐn)?shù)據(jù)集作為實驗數(shù)據(jù),數(shù)據(jù)集信息如表1所示。

    表1 數(shù)據(jù)集中所有植入模式的詳細(xì)信息

    UCR數(shù)據(jù)集是由事先確定好模式長度的已知模式組成,將UCR數(shù)據(jù)集中已知模式隨機(jī)植入到隨機(jī)游走數(shù)據(jù)中,創(chuàng)建實驗所用數(shù)據(jù)集Dataset。

    4.2 實驗方法

    本文從模體發(fā)現(xiàn)的時間效率和準(zhǔn)確率這2個方面來分析算法。

    實驗1驗證本文算法的準(zhǔn)確性。將本文提出的方法與MN方法[23]以及原始VLMD算法進(jìn)行比較。采用準(zhǔn)確性檢測方法AoD[24]作為準(zhǔn)確性評價指標(biāo),該方法最初用于計算子序列匹配問題中任意子序列對的重疊百分比。本實驗用其計算算法輸出模體與植入模體間的重疊比,以衡量各算法的準(zhǔn)確性。

    (7)

    其中,

    max(Lx,Ly)+1

    (8)

    min(LX,Ly)+1

    (9)

    其中,Lx、Ly為模體開始位置,Sx、Sy為模體長度。

    實驗2驗證不同算法的時間效率。增加算法中待挖掘的時間序列模體長度,對比VLMD的子程序MK算法和MPVLMD的子程序STOMP算法,獲取不同長度模體所需的時間。

    4.3 實驗結(jié)果與分析

    1)準(zhǔn)確性分析。

    為了驗證MPVLMD算法的準(zhǔn)確性,基于2個不同的數(shù)據(jù)集,分別使用本文算法、MN算法和VLMD算法進(jìn)行多次模體挖掘?qū)嶒灒y(tǒng)計各算法挖掘模體的結(jié)果。并基于本文選用的準(zhǔn)確性衡量方法AoD,計算各算法輸出模體與預(yù)先植入模體的重疊比。所得計算結(jié)果如表2所示,具體展示見圖3。

    表2 各數(shù)據(jù)集中不同算法發(fā)現(xiàn)植入模體的準(zhǔn)確率

    圖3 不同算法發(fā)現(xiàn)植入模體的準(zhǔn)確率

    由圖3可知,MPVLMD算法能夠發(fā)現(xiàn)所有的植入模體,其發(fā)現(xiàn)模體的準(zhǔn)確率要優(yōu)于VLMD算法。同時,針對Dataset中的Beef數(shù)據(jù),VLMD算法會出現(xiàn)不能發(fā)現(xiàn)植入模體的情況,這是因為模體候選集中存在過長、過短和平凡匹配模體影響算法發(fā)現(xiàn)模體的整體準(zhǔn)確性,這也表明本文算法引入的長度相似性和模體分組等價類條件在篩除無效模體上的積極作用。雖然MN算法也能有效地發(fā)現(xiàn)所有的植入模體,但是其發(fā)現(xiàn)模體的準(zhǔn)確率整體來看要略低于MPVLMD算法,表明本文算法不僅能夠有效地發(fā)現(xiàn)不同長度的模體,而且具有較高的準(zhǔn)確率。

    2)時間效率分析。

    為了驗證MPVLMD算法所選用的模體發(fā)現(xiàn)算法STOMP較MK算法的優(yōu)勢,基于Dataset數(shù)據(jù)集,取算法中待挖掘的時間序列模體長度分別為64、128、256、512、1024進(jìn)行實驗,記錄STOMP算法和MK算法獲取不同長度模體所需的平均時間,如圖4所示。

    圖4 模體發(fā)現(xiàn)耗時

    由圖4可知,在模體長度相同的情況下,STOMP算法的效率遠(yuǎn)優(yōu)于MK算法。并且隨著模體長度和數(shù)據(jù)集長度的增加,MK算法運(yùn)行所需時間呈現(xiàn)指數(shù)增長趨勢,其效率快速降低,相反STOMP算法始終能夠保持較高的運(yùn)行速度,其運(yùn)行時間對模體和數(shù)據(jù)集長度變化不敏感。

    綜合上述2個實驗的結(jié)果可知,MPVLMD算法能夠有效地發(fā)現(xiàn)時間序列中不同長度的模體,在準(zhǔn)確率和時間效率方面均優(yōu)于原始的VLMD算法。

    5 結(jié)束語

    本文提出了一種基于Matrix Profile的變長時間序列模體挖掘算法MPVLMD,該算法從所有可能長度的模體中自適應(yīng)返回合適長度的模體。通過模體提取、模體分組、模體分組等價類劃分和變長模體提取這4個步驟,將大量可能的滑動窗口長度縮減到真正有意義的幾個模體長度,能夠有效地發(fā)現(xiàn)時間序列中不同長度的模體。

    基于UCR數(shù)據(jù)集的實驗結(jié)果表明,本文算法具有較高的準(zhǔn)確率和時間性能。

    猜你喜歡
    模體歐氏等價
    植入(l, d)模體發(fā)現(xiàn)若干算法的實現(xiàn)與比較
    n次自然數(shù)冪和的一個等價無窮大
    中文信息(2017年12期)2018-01-27 08:22:58
    基于網(wǎng)絡(luò)模體特征攻擊的網(wǎng)絡(luò)抗毀性研究
    基于模體演化的時序鏈路預(yù)測方法
    收斂的非線性迭代數(shù)列xn+1=g(xn)的等價數(shù)列
    環(huán)Fpm+uFpm+…+uk-1Fpm上常循環(huán)碼的等價性
    基于多維歐氏空間相似度的激光點(diǎn)云分割方法
    一種基于信息容量的模體比較非比對度量算法
    麗江“思奔記”(上)
    探索地理(2013年5期)2014-01-09 06:40:44
    關(guān)于環(huán)Fpm+uFpm上常循環(huán)碼的等價性
    午夜激情久久久久久久| 成人漫画全彩无遮挡| 国产成人免费观看mmmm| 欧美亚洲 丝袜 人妻 在线| 国产av码专区亚洲av| 成年人免费黄色播放视频| 国产不卡av网站在线观看| 日本欧美视频一区| 18禁观看日本| 久久人人爽人人片av| 亚洲欧美精品综合一区二区三区| 亚洲,一卡二卡三卡| 男女国产视频网站| 亚洲欧洲国产日韩| 久久性视频一级片| 母亲3免费完整高清在线观看| 亚洲欧美一区二区三区黑人| 久久久久久久久久久免费av| 亚洲精品美女久久av网站| 亚洲精华国产精华液的使用体验| 一级爰片在线观看| 国产乱人偷精品视频| 三上悠亚av全集在线观看| 欧美日韩av久久| 美女福利国产在线| 丁香六月欧美| 91老司机精品| 国产一区二区 视频在线| www.熟女人妻精品国产| 国产成人免费观看mmmm| 免费在线观看黄色视频的| 亚洲精品视频女| 日韩熟女老妇一区二区性免费视频| 亚洲视频免费观看视频| 青春草视频在线免费观看| 国产亚洲av片在线观看秒播厂| 久久久久久人人人人人| 精品少妇一区二区三区视频日本电影 | 色94色欧美一区二区| 国产在线免费精品| 国产不卡av网站在线观看| 丝袜在线中文字幕| 毛片一级片免费看久久久久| 久久久久精品久久久久真实原创| 91aial.com中文字幕在线观看| 国产成人精品久久久久久| 99国产综合亚洲精品| 国产精品av久久久久免费| 99国产精品免费福利视频| 1024香蕉在线观看| 一区二区日韩欧美中文字幕| 日韩电影二区| 欧美精品高潮呻吟av久久| 国产野战对白在线观看| 精品亚洲乱码少妇综合久久| 在线观看人妻少妇| 精品亚洲乱码少妇综合久久| 人人妻,人人澡人人爽秒播 | av在线播放精品| 天美传媒精品一区二区| 日本av手机在线免费观看| 国产人伦9x9x在线观看| 国产精品免费大片| 夫妻午夜视频| av一本久久久久| 亚洲国产欧美在线一区| 日日撸夜夜添| 亚洲av成人精品一二三区| 在线观看免费午夜福利视频| 国产精品成人在线| 免费高清在线观看视频在线观看| 亚洲精品视频女| 国产成人精品久久二区二区91 | 国产黄色免费在线视频| 婷婷色麻豆天堂久久| 精品福利永久在线观看| 天天躁夜夜躁狠狠久久av| 成年美女黄网站色视频大全免费| 国产男女超爽视频在线观看| 日韩熟女老妇一区二区性免费视频| 99国产精品免费福利视频| 十八禁网站网址无遮挡| 成人亚洲欧美一区二区av| videosex国产| 大陆偷拍与自拍| 丝袜美腿诱惑在线| 天天添夜夜摸| 国产伦人伦偷精品视频| 亚洲国产精品一区二区三区在线| 国产福利在线免费观看视频| 满18在线观看网站| 亚洲国产精品一区二区三区在线| 人成视频在线观看免费观看| 国产成人精品无人区| 男女之事视频高清在线观看 | 国产精品蜜桃在线观看| 亚洲国产欧美日韩在线播放| 国产日韩欧美视频二区| 亚洲精品国产区一区二| 丝袜在线中文字幕| 亚洲国产欧美日韩在线播放| 欧美日韩亚洲高清精品| www.自偷自拍.com| 亚洲伊人色综图| 嫩草影视91久久| 男女午夜视频在线观看| 亚洲熟女毛片儿| 丝袜在线中文字幕| 在线观看www视频免费| 免费黄网站久久成人精品| 婷婷成人精品国产| 三上悠亚av全集在线观看| 亚洲少妇的诱惑av| 精品国产一区二区三区四区第35| 考比视频在线观看| 亚洲国产精品999| 18禁国产床啪视频网站| 国产伦理片在线播放av一区| 十八禁网站网址无遮挡| 亚洲中文av在线| 国产高清不卡午夜福利| 午夜福利,免费看| 人人妻人人爽人人添夜夜欢视频| 久久天堂一区二区三区四区| 久久久欧美国产精品| 午夜免费男女啪啪视频观看| 久久婷婷青草| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产毛片av蜜桃av| 国产野战对白在线观看| 国产深夜福利视频在线观看| 久久人人爽人人片av| 人人澡人人妻人| 狂野欧美激情性xxxx| 最新的欧美精品一区二区| 午夜日本视频在线| 国产黄色免费在线视频| 一级片免费观看大全| 最新的欧美精品一区二区| 欧美日韩一区二区视频在线观看视频在线| 中文字幕人妻丝袜制服| 天天操日日干夜夜撸| 狠狠婷婷综合久久久久久88av| 国产免费福利视频在线观看| 夫妻性生交免费视频一级片| 永久免费av网站大全| 国产精品av久久久久免费| 母亲3免费完整高清在线观看| 18在线观看网站| 午夜免费观看性视频| 欧美成人精品欧美一级黄| 伊人久久国产一区二区| 男女无遮挡免费网站观看| 国产又爽黄色视频| 黄片小视频在线播放| 女性生殖器流出的白浆| 天天添夜夜摸| 99精国产麻豆久久婷婷| 美女大奶头黄色视频| 这个男人来自地球电影免费观看 | 国精品久久久久久国模美| 99久久综合免费| 亚洲精品日本国产第一区| 97精品久久久久久久久久精品| 亚洲精品在线美女| 男女午夜视频在线观看| 国产不卡av网站在线观看| 波野结衣二区三区在线| 伊人亚洲综合成人网| 中文天堂在线官网| 一级片免费观看大全| 黄片无遮挡物在线观看| 精品国产乱码久久久久久小说| 午夜老司机福利片| 国产精品三级大全| 精品第一国产精品| 王馨瑶露胸无遮挡在线观看| 国产精品免费视频内射| 国产欧美日韩一区二区三区在线| 亚洲欧美色中文字幕在线| 人妻 亚洲 视频| 久久av网站| 亚洲中文av在线| 亚洲av日韩在线播放| 一区在线观看完整版| 成人三级做爰电影| 精品第一国产精品| 少妇猛男粗大的猛烈进出视频| 天天躁日日躁夜夜躁夜夜| 久久久久精品性色| 国产无遮挡羞羞视频在线观看| 人人妻,人人澡人人爽秒播 | 美国免费a级毛片| 亚洲精品成人av观看孕妇| 老汉色∧v一级毛片| 国产xxxxx性猛交| 岛国毛片在线播放| 国产伦人伦偷精品视频| 波多野结衣一区麻豆| 女人精品久久久久毛片| 国产精品三级大全| 亚洲婷婷狠狠爱综合网| 麻豆精品久久久久久蜜桃| 亚洲少妇的诱惑av| 亚洲欧美成人精品一区二区| 亚洲第一区二区三区不卡| 亚洲av欧美aⅴ国产| 在线观看免费日韩欧美大片| 国产97色在线日韩免费| 综合色丁香网| 国语对白做爰xxxⅹ性视频网站| 一边摸一边做爽爽视频免费| 熟妇人妻不卡中文字幕| 日本欧美视频一区| 精品酒店卫生间| 男女边摸边吃奶| 日本一区二区免费在线视频| 一级毛片电影观看| 亚洲欧美清纯卡通| 啦啦啦啦在线视频资源| 午夜免费观看性视频| 久久精品人人爽人人爽视色| 国产精品无大码| 大码成人一级视频| 99久久99久久久精品蜜桃| 在线观看www视频免费| 国产又爽黄色视频| 五月开心婷婷网| 只有这里有精品99| 精品久久蜜臀av无| 男男h啪啪无遮挡| 热re99久久精品国产66热6| 巨乳人妻的诱惑在线观看| 黄色一级大片看看| 亚洲精品久久久久久婷婷小说| 色综合欧美亚洲国产小说| 久久韩国三级中文字幕| 亚洲三区欧美一区| 久久国产精品男人的天堂亚洲| 午夜福利免费观看在线| 免费高清在线观看视频在线观看| 亚洲精品国产一区二区精华液| 中文字幕高清在线视频| 国产乱来视频区| 搡老岳熟女国产| 人人妻人人澡人人爽人人夜夜| 国产精品偷伦视频观看了| 久久久精品94久久精品| 不卡av一区二区三区| 丝袜脚勾引网站| 欧美精品高潮呻吟av久久| 五月开心婷婷网| 哪个播放器可以免费观看大片| 国产精品三级大全| 欧美 亚洲 国产 日韩一| 亚洲专区中文字幕在线 | av片东京热男人的天堂| 亚洲情色 制服丝袜| 水蜜桃什么品种好| av在线app专区| 国产欧美亚洲国产| 国产精品偷伦视频观看了| 国产精品欧美亚洲77777| 超碰成人久久| 新久久久久国产一级毛片| 女人久久www免费人成看片| 中文欧美无线码| 一级毛片电影观看| 午夜激情av网站| 欧美av亚洲av综合av国产av | 亚洲精品国产一区二区精华液| 婷婷色综合www| 十八禁高潮呻吟视频| av有码第一页| 久久精品亚洲av国产电影网| av天堂久久9| 黄色视频在线播放观看不卡| tube8黄色片| 人体艺术视频欧美日本| 日韩精品有码人妻一区| 久久久久久久精品精品| 久久久精品区二区三区| 亚洲视频免费观看视频| 母亲3免费完整高清在线观看| 国产男人的电影天堂91| 国语对白做爰xxxⅹ性视频网站| 日韩制服丝袜自拍偷拍| 99国产精品免费福利视频| 亚洲人成电影观看| 免费av中文字幕在线| 老司机在亚洲福利影院| 亚洲国产看品久久| 精品国产露脸久久av麻豆| 久久精品亚洲av国产电影网| 久久影院123| 一边摸一边抽搐一进一出视频| tube8黄色片| 国产伦人伦偷精品视频| 日韩制服骚丝袜av| 欧美久久黑人一区二区| 精品少妇一区二区三区视频日本电影 | 9热在线视频观看99| 国产一区有黄有色的免费视频| 免费观看人在逋| 操美女的视频在线观看| 男男h啪啪无遮挡| 老熟女久久久| 亚洲国产日韩一区二区| 久久精品亚洲av国产电影网| 成人漫画全彩无遮挡| 一边摸一边做爽爽视频免费| 深夜精品福利| 哪个播放器可以免费观看大片| 久久青草综合色| 老汉色∧v一级毛片| 老汉色av国产亚洲站长工具| 亚洲美女黄色视频免费看| 成人漫画全彩无遮挡| 一边摸一边做爽爽视频免费| 蜜桃在线观看..| 十八禁人妻一区二区| 国产人伦9x9x在线观看| 亚洲熟女毛片儿| 亚洲欧美成人精品一区二区| 欧美日韩成人在线一区二区| 亚洲av电影在线观看一区二区三区| 国产精品国产三级国产专区5o| 免费在线观看完整版高清| 国产黄频视频在线观看| 黄片播放在线免费| 日韩中文字幕欧美一区二区 | 国产又爽黄色视频| 国产福利在线免费观看视频| av又黄又爽大尺度在线免费看| 在现免费观看毛片| 国产淫语在线视频| 天天躁夜夜躁狠狠久久av| 国产成人啪精品午夜网站| 最近的中文字幕免费完整| 亚洲精品一二三| 亚洲美女黄色视频免费看| 日本爱情动作片www.在线观看| 日韩中文字幕欧美一区二区 | 免费少妇av软件| 中文天堂在线官网| 欧美另类一区| 亚洲欧美激情在线| 国产伦理片在线播放av一区| 国产一区二区三区av在线| 日日爽夜夜爽网站| 国产精品久久久久久久久免| 美女高潮到喷水免费观看| 两性夫妻黄色片| 亚洲精品久久午夜乱码| 亚洲精品国产av成人精品| 亚洲图色成人| 熟女av电影| 永久免费av网站大全| 国产亚洲欧美精品永久| 最近2019中文字幕mv第一页| 一区二区日韩欧美中文字幕| 免费高清在线观看视频在线观看| 天美传媒精品一区二区| 久久人人97超碰香蕉20202| 美女脱内裤让男人舔精品视频| 国产精品人妻久久久影院| 下体分泌物呈黄色| 久久久欧美国产精品| 国产 一区精品| 亚洲国产欧美网| 国产日韩欧美在线精品| 男女边摸边吃奶| 久久99热这里只频精品6学生| 国产爽快片一区二区三区| 欧美日韩综合久久久久久| xxx大片免费视频| svipshipincom国产片| 老司机靠b影院| 久久女婷五月综合色啪小说| 精品久久蜜臀av无| 久久人人爽人人片av| 亚洲国产日韩一区二区| 成人国产麻豆网| 国产在视频线精品| 一区福利在线观看| 超色免费av| 亚洲色图综合在线观看| 国产日韩欧美亚洲二区| 男女之事视频高清在线观看 | 香蕉丝袜av| 男女国产视频网站| 黄色 视频免费看| 亚洲成色77777| 可以免费在线观看a视频的电影网站 | 国产精品久久久久久久久免| 国产在线免费精品| 可以免费在线观看a视频的电影网站 | 亚洲精品第二区| 亚洲国产日韩一区二区| av福利片在线| 午夜91福利影院| 青春草视频在线免费观看| 视频在线观看一区二区三区| 精品少妇久久久久久888优播| 日韩制服丝袜自拍偷拍| 欧美日韩亚洲综合一区二区三区_| 中文字幕色久视频| 最近的中文字幕免费完整| 在线精品无人区一区二区三| 亚洲一级一片aⅴ在线观看| 秋霞在线观看毛片| av有码第一页| 国产一区亚洲一区在线观看| 日韩中文字幕欧美一区二区 | 国产精品香港三级国产av潘金莲 | 国产免费福利视频在线观看| 成人漫画全彩无遮挡| 伊人久久国产一区二区| 成人毛片60女人毛片免费| av一本久久久久| 欧美日韩精品网址| 亚洲男人天堂网一区| 亚洲美女视频黄频| 一边摸一边抽搐一进一出视频| 不卡视频在线观看欧美| 日韩大片免费观看网站| 亚洲五月色婷婷综合| 久久国产精品大桥未久av| 免费av中文字幕在线| 女人被躁到高潮嗷嗷叫费观| 国产av码专区亚洲av| 精品一区在线观看国产| 人人妻,人人澡人人爽秒播 | 亚洲美女搞黄在线观看| 中文天堂在线官网| 国产免费视频播放在线视频| 久久国产精品男人的天堂亚洲| 制服诱惑二区| 欧美国产精品va在线观看不卡| 宅男免费午夜| 成人国产av品久久久| 丰满乱子伦码专区| 中文字幕亚洲精品专区| 国产福利在线免费观看视频| 亚洲色图综合在线观看| 丝袜喷水一区| 国产精品成人在线| 最近手机中文字幕大全| 欧美激情 高清一区二区三区| 操美女的视频在线观看| 亚洲欧美日韩另类电影网站| 国产高清不卡午夜福利| 久久ye,这里只有精品| 伊人久久国产一区二区| 亚洲av欧美aⅴ国产| 丁香六月欧美| 秋霞伦理黄片| 免费不卡黄色视频| 久久久久久免费高清国产稀缺| 捣出白浆h1v1| 国产麻豆69| 久久人人爽人人片av| 国产熟女欧美一区二区| 久久热在线av| 天天躁狠狠躁夜夜躁狠狠躁| 波多野结衣一区麻豆| 亚洲一区中文字幕在线| 国产在线免费精品| 男女高潮啪啪啪动态图| 韩国精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 最近手机中文字幕大全| 久久精品国产亚洲av涩爱| 国产男女内射视频| 熟妇人妻不卡中文字幕| 18在线观看网站| 夜夜骑夜夜射夜夜干| 捣出白浆h1v1| 国产精品 欧美亚洲| 精品人妻一区二区三区麻豆| 久久精品亚洲av国产电影网| 久久久久网色| 女性生殖器流出的白浆| 久久午夜综合久久蜜桃| 一级爰片在线观看| 看非洲黑人一级黄片| 日韩 亚洲 欧美在线| 国产av码专区亚洲av| 丰满迷人的少妇在线观看| 国产免费又黄又爽又色| 欧美日韩一级在线毛片| 亚洲精品成人av观看孕妇| 最近中文字幕2019免费版| 色视频在线一区二区三区| 久久狼人影院| 七月丁香在线播放| 黄网站色视频无遮挡免费观看| 亚洲欧美一区二区三区黑人| 卡戴珊不雅视频在线播放| 午夜福利,免费看| 国产男女内射视频| 91成人精品电影| 观看av在线不卡| 秋霞在线观看毛片| 亚洲欧美一区二区三区黑人| 国产免费视频播放在线视频| 宅男免费午夜| 色综合欧美亚洲国产小说| 亚洲精品一区蜜桃| 两性夫妻黄色片| 免费黄色在线免费观看| 你懂的网址亚洲精品在线观看| 一区二区三区四区激情视频| 欧美日韩视频精品一区| 国产精品一区二区精品视频观看| 午夜福利在线免费观看网站| 人人妻人人添人人爽欧美一区卜| 成年动漫av网址| 可以免费在线观看a视频的电影网站 | 精品亚洲成国产av| 久久精品久久久久久噜噜老黄| 九色亚洲精品在线播放| 亚洲一码二码三码区别大吗| 午夜福利影视在线免费观看| 五月开心婷婷网| 黑人欧美特级aaaaaa片| 考比视频在线观看| av又黄又爽大尺度在线免费看| 精品亚洲乱码少妇综合久久| 亚洲一卡2卡3卡4卡5卡精品中文| 秋霞在线观看毛片| 欧美在线一区亚洲| 老司机深夜福利视频在线观看 | 久久影院123| 考比视频在线观看| 一边亲一边摸免费视频| 亚洲成人免费av在线播放| 亚洲av成人精品一二三区| 三上悠亚av全集在线观看| 在线观看免费日韩欧美大片| 水蜜桃什么品种好| 午夜福利网站1000一区二区三区| 日韩免费高清中文字幕av| 久久国产精品大桥未久av| 一区二区三区四区激情视频| 熟女av电影| 制服人妻中文乱码| 考比视频在线观看| 又粗又硬又长又爽又黄的视频| 久久久国产欧美日韩av| 国产成人91sexporn| 免费黄网站久久成人精品| 哪个播放器可以免费观看大片| 国产精品国产三级国产专区5o| 亚洲国产欧美一区二区综合| 亚洲欧美精品综合一区二区三区| 黑丝袜美女国产一区| 亚洲国产欧美网| 国产成人精品在线电影| 久久青草综合色| 亚洲欧美成人精品一区二区| 成人影院久久| 亚洲人成网站在线观看播放| 亚洲国产中文字幕在线视频| 久久狼人影院| 777米奇影视久久| 免费少妇av软件| 狂野欧美激情性xxxx| 丝袜喷水一区| 色吧在线观看| 国产亚洲av片在线观看秒播厂| 精品久久久久久电影网| 女的被弄到高潮叫床怎么办| 成人黄色视频免费在线看| 亚洲精品国产区一区二| 国产精品成人在线| 免费黄色在线免费观看| 久久久亚洲精品成人影院| 国产精品一区二区精品视频观看| 高清黄色对白视频在线免费看| 婷婷色综合www| 一级毛片黄色毛片免费观看视频| 国语对白做爰xxxⅹ性视频网站| 夫妻性生交免费视频一级片| 99久久综合免费| 欧美最新免费一区二区三区| 99久久99久久久精品蜜桃| 免费高清在线观看日韩| 亚洲欧美一区二区三区久久| 亚洲人成77777在线视频| 国产成人免费无遮挡视频| 五月开心婷婷网| 午夜日韩欧美国产| 欧美最新免费一区二区三区| 日韩欧美精品免费久久| 香蕉丝袜av| 国产精品熟女久久久久浪| 精品免费久久久久久久清纯 | 最近中文字幕2019免费版| 欧美日韩国产mv在线观看视频| 90打野战视频偷拍视频| 9191精品国产免费久久| 叶爱在线成人免费视频播放| 免费在线观看视频国产中文字幕亚洲 | 大香蕉久久成人网| 性色av一级| 日本午夜av视频| 国产爽快片一区二区三区| 亚洲美女视频黄频| 午夜精品国产一区二区电影| 成人国产麻豆网| 嫩草影视91久久| 亚洲国产av新网站| 久久久国产欧美日韩av| 可以免费在线观看a视频的电影网站 |