• <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)碼的等價性
    videos熟女内射| 母亲3免费完整高清在线观看 | 一级毛片电影观看| 王馨瑶露胸无遮挡在线观看| 在线亚洲精品国产二区图片欧美| 亚洲欧美色中文字幕在线| 国产在线视频一区二区| 亚洲色图 男人天堂 中文字幕| 男女高潮啪啪啪动态图| 一级a爱视频在线免费观看| 久久精品国产鲁丝片午夜精品| 不卡视频在线观看欧美| 一区在线观看完整版| 青春草亚洲视频在线观看| www日本在线高清视频| 老熟女久久久| 国产成人精品福利久久| 99香蕉大伊视频| 国产成人91sexporn| 高清不卡的av网站| 国产欧美日韩一区二区三区在线| 亚洲,欧美,日韩| 老鸭窝网址在线观看| 中文精品一卡2卡3卡4更新| 丰满少妇做爰视频| 国产亚洲av片在线观看秒播厂| 久久精品人人爽人人爽视色| 日韩中字成人| 校园人妻丝袜中文字幕| 午夜精品国产一区二区电影| 久久女婷五月综合色啪小说| 99热网站在线观看| av电影中文网址| 国产精品久久久久久久久免| 亚洲精品美女久久av网站| 国产精品成人在线| 看免费av毛片| 欧美 日韩 精品 国产| 国产毛片在线视频| 99香蕉大伊视频| 亚洲 欧美一区二区三区| 啦啦啦啦在线视频资源| 99国产精品免费福利视频| 欧美日韩成人在线一区二区| 亚洲国产精品一区三区| 韩国av在线不卡| 国产av一区二区精品久久| 人成视频在线观看免费观看| 国产亚洲欧美精品永久| 久久 成人 亚洲| 毛片一级片免费看久久久久| 国产精品成人在线| av福利片在线| 两性夫妻黄色片| 亚洲av男天堂| 国产一区二区 视频在线| 亚洲一区中文字幕在线| 热re99久久精品国产66热6| 巨乳人妻的诱惑在线观看| 国产白丝娇喘喷水9色精品| 精品人妻偷拍中文字幕| 久久精品久久久久久噜噜老黄| 国产精品欧美亚洲77777| 热re99久久国产66热| 精品亚洲成国产av| 99国产精品免费福利视频| 国产成人精品久久二区二区91 | 国产一区二区在线观看av| 一级毛片 在线播放| 丰满迷人的少妇在线观看| 99久国产av精品国产电影| 99久久人妻综合| 国产探花极品一区二区| 黄片无遮挡物在线观看| 男女国产视频网站| 狂野欧美激情性bbbbbb| 亚洲三级黄色毛片| 黑丝袜美女国产一区| av免费观看日本| www.精华液| 一区二区三区乱码不卡18| 韩国精品一区二区三区| 热99久久久久精品小说推荐| 亚洲成人av在线免费| 欧美日韩成人在线一区二区| 国产精品 欧美亚洲| 热re99久久国产66热| 国产乱来视频区| 一区二区日韩欧美中文字幕| 日本免费在线观看一区| 街头女战士在线观看网站| 精品久久蜜臀av无| 国产免费视频播放在线视频| 黄网站色视频无遮挡免费观看| 丁香六月天网| 国产一级毛片在线| 国产色婷婷99| 秋霞伦理黄片| 免费高清在线观看日韩| freevideosex欧美| 色吧在线观看| 亚洲色图综合在线观看| 婷婷色综合www| 久久久久久人人人人人| 尾随美女入室| 日韩av不卡免费在线播放| 美女国产视频在线观看| 777米奇影视久久| 在线观看国产h片| 久久久久网色| 丝袜喷水一区| 美女高潮到喷水免费观看| av免费观看日本| 亚洲,一卡二卡三卡| 中文字幕精品免费在线观看视频| 老汉色∧v一级毛片| 婷婷成人精品国产| 亚洲av综合色区一区| 久久精品国产亚洲av涩爱| 免费高清在线观看视频在线观看| 欧美精品av麻豆av| 亚洲综合精品二区| 韩国精品一区二区三区| av女优亚洲男人天堂| 又粗又硬又长又爽又黄的视频| 一级毛片我不卡| a 毛片基地| tube8黄色片| 国产精品二区激情视频| 国产成人精品久久二区二区91 | 秋霞伦理黄片| 国产熟女欧美一区二区| 天天躁夜夜躁狠狠躁躁| 日韩人妻精品一区2区三区| 日韩伦理黄色片| 欧美另类一区| 好男人视频免费观看在线| 精品第一国产精品| 久久国内精品自在自线图片| 999精品在线视频| 国产午夜精品一二区理论片| 色网站视频免费| 日本色播在线视频| 美女xxoo啪啪120秒动态图| 男人爽女人下面视频在线观看| 亚洲精品一区蜜桃| 黄片小视频在线播放| 大香蕉久久成人网| 日韩不卡一区二区三区视频在线| 啦啦啦视频在线资源免费观看| 婷婷色av中文字幕| 在线精品无人区一区二区三| 免费人妻精品一区二区三区视频| 亚洲av电影在线观看一区二区三区| 搡老乐熟女国产| 人妻少妇偷人精品九色| 国产又爽黄色视频| 亚洲美女视频黄频| 成人影院久久| 日韩中文字幕欧美一区二区 | 成人午夜精彩视频在线观看| 国产成人精品福利久久| 免费高清在线观看视频在线观看| 日韩大片免费观看网站| 午夜日韩欧美国产| 亚洲国产色片| 精品人妻熟女毛片av久久网站| 国产一区二区在线观看av| 又黄又粗又硬又大视频| 18+在线观看网站| 观看av在线不卡| 美女福利国产在线| 少妇的逼水好多| 欧美日韩一级在线毛片| 亚洲精品日本国产第一区| 亚洲第一青青草原| 国产97色在线日韩免费| 丰满饥渴人妻一区二区三| 成人亚洲精品一区在线观看| 色哟哟·www| 成人二区视频| 国产福利在线免费观看视频| 久久99热这里只频精品6学生| 女人高潮潮喷娇喘18禁视频| 亚洲人成电影观看| 午夜福利影视在线免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91在线精品国自产拍蜜月| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 91精品伊人久久大香线蕉| 日韩视频在线欧美| 最新中文字幕久久久久| 97精品久久久久久久久久精品| 久久鲁丝午夜福利片| 免费不卡的大黄色大毛片视频在线观看| 在线免费观看不下载黄p国产| 国产一区二区激情短视频 | 欧美人与善性xxx| 日韩一卡2卡3卡4卡2021年| 国产精品免费视频内射| 大陆偷拍与自拍| 中文字幕最新亚洲高清| 在线观看一区二区三区激情| 久久精品久久久久久噜噜老黄| 男人爽女人下面视频在线观看| 日韩,欧美,国产一区二区三区| 中文字幕色久视频| 中文字幕制服av| 成人毛片a级毛片在线播放| 久久99热这里只频精品6学生| 国产成人精品一,二区| 91国产中文字幕| 我要看黄色一级片免费的| 欧美老熟妇乱子伦牲交| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av综合色区一区| 天堂中文最新版在线下载| 人人妻人人添人人爽欧美一区卜| 青春草视频在线免费观看| 69精品国产乱码久久久| 一二三四在线观看免费中文在| 国产一区二区激情短视频 | 久久久国产精品麻豆| 欧美少妇被猛烈插入视频| 亚洲图色成人| av网站免费在线观看视频| 男男h啪啪无遮挡| 80岁老熟妇乱子伦牲交| 1024香蕉在线观看| 男女边摸边吃奶| 国产免费视频播放在线视频| 亚洲,欧美精品.| 一本色道久久久久久精品综合| 日韩中文字幕欧美一区二区 | 国产免费现黄频在线看| 欧美日韩视频高清一区二区三区二| 在线观看人妻少妇| 一区在线观看完整版| 热re99久久国产66热| 欧美亚洲日本最大视频资源| 精品一区二区三卡| 午夜福利网站1000一区二区三区| 国产亚洲午夜精品一区二区久久| 黄色视频在线播放观看不卡| 一区福利在线观看| 母亲3免费完整高清在线观看 | 国产精品国产三级国产专区5o| 久久久久久久久久人人人人人人| 日韩成人av中文字幕在线观看| 国产xxxxx性猛交| 26uuu在线亚洲综合色| 水蜜桃什么品种好| 国产免费福利视频在线观看| 午夜91福利影院| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲久久久国产精品| 国产成人精品久久二区二区91 | 国产精品久久久久久av不卡| 黑人巨大精品欧美一区二区蜜桃| 中文字幕人妻丝袜制服| 精品午夜福利在线看| 国产 精品1| 欧美黄色片欧美黄色片| 国产精品偷伦视频观看了| 天美传媒精品一区二区| 国产精品国产三级国产专区5o| 人妻一区二区av| 街头女战士在线观看网站| 最近中文字幕高清免费大全6| 午夜精品国产一区二区电影| 亚洲成色77777| 爱豆传媒免费全集在线观看| 久久毛片免费看一区二区三区| 欧美日韩成人在线一区二区| 国产高清不卡午夜福利| 欧美日韩一级在线毛片| 99九九在线精品视频| 男女无遮挡免费网站观看| 日韩人妻精品一区2区三区| √禁漫天堂资源中文www| 亚洲欧美日韩另类电影网站| 亚洲精品一二三| 波多野结衣一区麻豆| 丝袜美足系列| 久久久国产欧美日韩av| 午夜日本视频在线| 我的亚洲天堂| a级片在线免费高清观看视频| 美女国产视频在线观看| 晚上一个人看的免费电影| 天天躁夜夜躁狠狠久久av| 我要看黄色一级片免费的| 90打野战视频偷拍视频| 卡戴珊不雅视频在线播放| 亚洲,欧美,日韩| 男女高潮啪啪啪动态图| 80岁老熟妇乱子伦牲交| 国产精品欧美亚洲77777| 天天影视国产精品| 最近最新中文字幕免费大全7| 欧美日韩亚洲高清精品| 日本色播在线视频| 国产成人欧美| 国产熟女午夜一区二区三区| 久热久热在线精品观看| 制服诱惑二区| 香蕉精品网在线| av在线观看视频网站免费| 亚洲精品日本国产第一区| 在线天堂最新版资源| 午夜福利视频精品| 人人妻人人澡人人爽人人夜夜| 综合色丁香网| 捣出白浆h1v1| 精品少妇久久久久久888优播| 少妇被粗大的猛进出69影院| 久久久精品94久久精品| 热99久久久久精品小说推荐| 久久久精品94久久精品| 久久毛片免费看一区二区三区| 亚洲精品国产一区二区精华液| tube8黄色片| 亚洲av男天堂| 涩涩av久久男人的天堂| 国产欧美日韩综合在线一区二区| 久久久久久人人人人人| 日韩中字成人| 自线自在国产av| 看非洲黑人一级黄片| 国产成人欧美| 99香蕉大伊视频| 久久久久久免费高清国产稀缺| 2022亚洲国产成人精品| 国产日韩欧美视频二区| 免费人妻精品一区二区三区视频| 精品国产超薄肉色丝袜足j| 这个男人来自地球电影免费观看 | 999精品在线视频| 一区二区三区精品91| 九色亚洲精品在线播放| 午夜日本视频在线| 黄网站色视频无遮挡免费观看| 午夜日本视频在线| 黄网站色视频无遮挡免费观看| 视频在线观看一区二区三区| 青春草国产在线视频| 亚洲精品自拍成人| 韩国av在线不卡| 国产精品 国内视频| 久久久久精品性色| 男女边摸边吃奶| 亚洲欧美精品综合一区二区三区 | 99久久精品国产国产毛片| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区久久| 你懂的网址亚洲精品在线观看| 一区二区三区乱码不卡18| 最近中文字幕2019免费版| 国产免费一区二区三区四区乱码| 亚洲成人av在线免费| 99九九在线精品视频| 黑丝袜美女国产一区| 日韩欧美精品免费久久| 国产成人精品久久久久久| 亚洲av免费高清在线观看| 肉色欧美久久久久久久蜜桃| 国产片特级美女逼逼视频| 男女下面插进去视频免费观看| 亚洲一区中文字幕在线| 亚洲精品视频女| 97精品久久久久久久久久精品| 日韩一本色道免费dvd| 在线免费观看不下载黄p国产| 亚洲一区二区三区欧美精品| 久久久久久免费高清国产稀缺| 国产一区二区三区av在线| 天堂俺去俺来也www色官网| 制服丝袜香蕉在线| 一区二区三区四区激情视频| 亚洲欧美精品综合一区二区三区 | 国产高清不卡午夜福利| 91久久精品国产一区二区三区| 久久久久久人妻| 母亲3免费完整高清在线观看 | 亚洲三级黄色毛片| 国产成人91sexporn| 超碰97精品在线观看| 午夜福利网站1000一区二区三区| 在线观看免费高清a一片| 亚洲男人天堂网一区| 欧美另类一区| 天天影视国产精品| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 亚洲av.av天堂| 久久久久久人妻| 黄片无遮挡物在线观看| 在线免费观看不下载黄p国产| 欧美精品人与动牲交sv欧美| 咕卡用的链子| 黄色一级大片看看| 免费观看性生交大片5| 熟女av电影| 搡老乐熟女国产| 少妇被粗大的猛进出69影院| 免费不卡的大黄色大毛片视频在线观看| 久久人人97超碰香蕉20202| 制服诱惑二区| 亚洲精品第二区| 免费观看在线日韩| 日日爽夜夜爽网站| 国产精品国产av在线观看| av免费在线看不卡| 国产精品国产三级国产专区5o| 精品国产乱码久久久久久男人| 久久99热这里只频精品6学生| 亚洲精品一二三| 国产成人91sexporn| 日韩中文字幕欧美一区二区 | 99久国产av精品国产电影| 久久国产精品大桥未久av| 热re99久久国产66热| 在线天堂中文资源库| 中国国产av一级| 男女无遮挡免费网站观看| 七月丁香在线播放| 国产成人aa在线观看| 你懂的网址亚洲精品在线观看| 日韩成人av中文字幕在线观看| 日韩电影二区| 69精品国产乱码久久久| 免费黄色在线免费观看| 成年美女黄网站色视频大全免费| 免费观看无遮挡的男女| 午夜福利视频在线观看免费| 免费看不卡的av| 精品久久蜜臀av无| 亚洲精品一区蜜桃| 精品久久久精品久久久| 两个人免费观看高清视频| 各种免费的搞黄视频| 欧美 日韩 精品 国产| 久久精品亚洲av国产电影网| 一级a爱视频在线免费观看| 欧美精品一区二区大全| 热99国产精品久久久久久7| 老司机影院毛片| 亚洲色图综合在线观看| 国产毛片在线视频| 久久久久人妻精品一区果冻| 国产av一区二区精品久久| 国产男女内射视频| 巨乳人妻的诱惑在线观看| 成年女人毛片免费观看观看9 | 国产xxxxx性猛交| 欧美bdsm另类| 亚洲成人一二三区av| 亚洲人成电影观看| 水蜜桃什么品种好| 成人毛片60女人毛片免费| 热99国产精品久久久久久7| 狠狠婷婷综合久久久久久88av| 91在线精品国自产拍蜜月| 久久久精品免费免费高清| 99久久中文字幕三级久久日本| 亚洲精品一区蜜桃| 亚洲国产av影院在线观看| 亚洲av在线观看美女高潮| 久久精品久久久久久噜噜老黄| 在线亚洲精品国产二区图片欧美| 日本午夜av视频| 欧美成人精品欧美一级黄| 久久99精品国语久久久| 久久久久国产精品人妻一区二区| 国产福利在线免费观看视频| 美国免费a级毛片| 99精国产麻豆久久婷婷| 少妇猛男粗大的猛烈进出视频| 美女高潮到喷水免费观看| 高清不卡的av网站| 国产成人精品福利久久| av在线播放精品| 亚洲久久久国产精品| 欧美日韩视频高清一区二区三区二| 高清视频免费观看一区二区| 久久综合国产亚洲精品| 欧美激情 高清一区二区三区| 制服人妻中文乱码| 尾随美女入室| 成人18禁高潮啪啪吃奶动态图| 国产精品不卡视频一区二区| 在线观看国产h片| 人人妻人人爽人人添夜夜欢视频| 亚洲情色 制服丝袜| a级片在线免费高清观看视频| 青春草亚洲视频在线观看| 亚洲,欧美精品.| 国产成人av激情在线播放| 国产精品二区激情视频| 人人妻人人添人人爽欧美一区卜| 深夜精品福利| 黑人欧美特级aaaaaa片| 午夜福利,免费看| 亚洲精品久久午夜乱码| 久久久久精品人妻al黑| 亚洲欧洲精品一区二区精品久久久 | 日韩人妻精品一区2区三区| 丝袜脚勾引网站| 成人毛片a级毛片在线播放| 国产片特级美女逼逼视频| 精品人妻偷拍中文字幕| 国产老妇伦熟女老妇高清| 高清黄色对白视频在线免费看| 成年人午夜在线观看视频| 欧美日韩一区二区视频在线观看视频在线| 久久精品国产综合久久久| 岛国毛片在线播放| 大码成人一级视频| tube8黄色片| 精品一区二区免费观看| 在线亚洲精品国产二区图片欧美| 26uuu在线亚洲综合色| 2022亚洲国产成人精品| 亚洲欧美清纯卡通| 高清黄色对白视频在线免费看| www日本在线高清视频| 亚洲欧美精品自产自拍| www日本在线高清视频| 国产福利在线免费观看视频| www日本在线高清视频| 男的添女的下面高潮视频| 欧美精品亚洲一区二区| 交换朋友夫妻互换小说| 一本久久精品| 你懂的网址亚洲精品在线观看| av电影中文网址| 国产有黄有色有爽视频| 超碰成人久久| 伊人亚洲综合成人网| 欧美日韩av久久| 麻豆av在线久日| 侵犯人妻中文字幕一二三四区| 丝袜人妻中文字幕| 97精品久久久久久久久久精品| 久久99精品国语久久久| 黄色 视频免费看| 免费播放大片免费观看视频在线观看| 亚洲图色成人| 精品久久蜜臀av无| 一区二区日韩欧美中文字幕| 一本色道久久久久久精品综合| 中文字幕亚洲精品专区| 精品久久久久久电影网| 天天影视国产精品| 国产亚洲午夜精品一区二区久久| 国产精品国产av在线观看| 亚洲国产精品一区三区| 国产探花极品一区二区| 国产精品久久久久久av不卡| 99国产综合亚洲精品| 中文乱码字字幕精品一区二区三区| 男女免费视频国产| 午夜福利影视在线免费观看| 亚洲一码二码三码区别大吗| a 毛片基地| 亚洲精品,欧美精品| 丝袜在线中文字幕| 亚洲美女视频黄频| 午夜影院在线不卡| 国产精品欧美亚洲77777| 久久久欧美国产精品| 在线观看美女被高潮喷水网站| 亚洲国产日韩一区二区| 91精品伊人久久大香线蕉| 欧美激情高清一区二区三区 | 国产成人aa在线观看| 一区二区三区精品91| 亚洲国产成人一精品久久久| 国产成人免费观看mmmm| 日韩电影二区| 国产精品女同一区二区软件| 欧美激情极品国产一区二区三区| av电影中文网址| 欧美精品人与动牲交sv欧美| 青春草国产在线视频| 国产国语露脸激情在线看| 波多野结衣av一区二区av| 欧美国产精品一级二级三级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一二三四在线观看免费中文在| 婷婷色麻豆天堂久久| 国产成人欧美| 色婷婷av一区二区三区视频| 久久狼人影院| 亚洲av电影在线观看一区二区三区| 亚洲国产精品999| 久久韩国三级中文字幕| 黄频高清免费视频| 免费黄色在线免费观看| 免费观看a级毛片全部| 极品少妇高潮喷水抽搐| 欧美中文综合在线视频| 欧美人与善性xxx| 99九九在线精品视频| 亚洲av欧美aⅴ国产| 亚洲精品在线美女| 亚洲av成人精品一二三区| 亚洲国产毛片av蜜桃av| 国产亚洲最大av| 女人久久www免费人成看片| 久久精品国产综合久久久| 男的添女的下面高潮视频| 热re99久久精品国产66热6| 亚洲综合色惰|