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

    基于2DSVD的多變量時(shí)間序列半監(jiān)督分類①

    2019-11-15 07:07:10單中南翁小清馬超紅
    關(guān)鍵詞:協(xié)方差分類器分類

    單中南,翁小清,馬超紅

    (河北經(jīng)貿(mào)大學(xué) 信息技術(shù)學(xué)院,石家莊 050061)

    時(shí)間序列是指按時(shí)間次序有序排列的一組數(shù)據(jù),任何有次序的實(shí)值序列都可當(dāng)作時(shí)間序列來處理[1].時(shí)間序列數(shù)據(jù)廣泛地存在于金融、醫(yī)學(xué)、交通等領(lǐng)域.建立準(zhǔn)確的分類器需要大量的有類別標(biāo)記的樣本數(shù)據(jù),然而在現(xiàn)實(shí)應(yīng)用領(lǐng)域,存在大量沒有類別標(biāo)記的樣本數(shù)據(jù),有標(biāo)記的樣本數(shù)據(jù)很難獲得,或用人工標(biāo)記樣本數(shù)據(jù)成本很高.半監(jiān)督分類(Semi-Supervised Classification,SSC)使用少量有標(biāo)記的樣本數(shù)據(jù)和大量未標(biāo)記的樣本數(shù)據(jù)建立分類器.目前,絕大多數(shù)時(shí)間序列半監(jiān)督分類的研究工作都集中在單變量時(shí)間序列,對多變量時(shí)間序列(Multivariate Time Series,MTS)的半監(jiān)督分類研究還比較少.在對MTS 進(jìn)行半監(jiān)督分類時(shí),主要遇到兩方面的困難[2]:第一,MTS 中含有多個(gè)變量,且變量之間存在復(fù)雜的相關(guān)系;第二,不同MTS 樣本它們的長度不一定相等,這些困難使得標(biāo)準(zhǔn)的分類器很難直接使用.本文針對MTS 特性,采用二維奇異值分解(Two-Dimensional Singular Value Decomposition,2DSVD) 從MTS 樣本中提取特征矩陣,并與其他MTS 半監(jiān)督分類方法進(jìn)行性能對比,討論該方法在MTS 半監(jiān)督分類的優(yōu)勢.本文第1 節(jié)介紹背景和相關(guān)工作;第2 節(jié)提出了基于2DSVD 的MTS半監(jiān)督分類算法;第3 節(jié)通過實(shí)驗(yàn)將本文提出的方法與其它半監(jiān)督分類方法進(jìn)行比較,并采用威爾克森符號秩檢驗(yàn)(Wilcoxon signed ranks test)對實(shí)驗(yàn)結(jié)果進(jìn)行對比,驗(yàn)證算法的有效性;第4 節(jié)給出了本文結(jié)論.

    1 背景和相關(guān)工作

    1.1 基本概念

    定義1.時(shí)間序列.時(shí)間序列是一段時(shí)間內(nèi)的一系列觀測值,用xi(t)[i=1,2,…,n;t=1,2,…,m]表示,其中m是觀測值的個(gè)數(shù),n是變量的個(gè)數(shù)[2].當(dāng)n=1 時(shí),稱為單變量時(shí)間序列,當(dāng)n≥2 時(shí),稱為多變量時(shí)間序列,通常用m×n矩陣存儲(chǔ)MTS.

    定義2.P集合.P為訓(xùn)練數(shù)據(jù)的一個(gè)集合,包括所有正類標(biāo)記的樣本[3].在訓(xùn)練開始時(shí),P只包含少量的正類樣本,或許只包含一個(gè)正類樣本.隨著學(xué)習(xí)的繼續(xù),先前U中一些沒有標(biāo)記的樣本,被標(biāo)記為正類樣本,并移動(dòng)到了P集合,P集合包含樣本的數(shù)量也隨之增加.最終,集合P既包含原來有標(biāo)記的正類樣本,也包括使用分類器從U 中選擇的樣本.

    定義3.U集合.U是未標(biāo)記樣本的集合[3].U中的樣本可以來自正類或者負(fù)類;通常情況下,U中的絕大多數(shù)樣本來自負(fù)類.

    1.2 二維奇異值分解

    Ding 等[4]對標(biāo)準(zhǔn)奇異值分解(即一維奇異值分解,1DSVD)進(jìn)行了擴(kuò)展,提出了基于行-行協(xié)方差矩陣以及列-列協(xié)方差矩陣的二維奇異值分解方法,2DSVD 是基于二維矩陣而不是基于一維向量[2].2DSVD 使用MTS 樣本構(gòu)造行-行以及列-列的協(xié)方差矩陣,然后計(jì)算行-行及列-列協(xié)方差矩陣的特征向量用于MTS 樣本特征矩陣的提取.使用2DSVD 提取出的MTS 樣本的特征矩陣,它們的行數(shù)以及列數(shù)不僅比原始數(shù)據(jù)低,而且還清晰地考慮了原始數(shù)據(jù)的二維特性.

    其中,||·||為L2 范數(shù).

    1.3 相關(guān)工作

    時(shí)間序列的半監(jiān)督分類方法可大致分為3 類[6,7]:基于實(shí)例、基于聚類以及基于模型的半監(jiān)督分類方法.

    Wei 等[8]針對正類中只有少量有標(biāo)記的樣本,使用歐氏距離建立基于最小最近鄰距離的分類器及停止準(zhǔn)則.Ratanamahatana[9]等使用DTW (Dynamic Time Warping)距離來改進(jìn)樣本的選取并提出了新的停止準(zhǔn)則,該準(zhǔn)則基于未標(biāo)記樣本集中候選樣本與正類樣本的歷史距離;Chen[3]等在SSC 算法中,使用一種基于DTW 和ED 相結(jié)合的特殊距離DTW-D,顯著地提高了分類的性能.Begum 等[10,11]提出了一種基于最小描述長度(Minimum Description Length,MDL)的停止準(zhǔn)則,該準(zhǔn)則利用數(shù)據(jù)的內(nèi)在性質(zhì)去發(fā)現(xiàn)停止點(diǎn);然而,時(shí)間序列在時(shí)間軸可能會(huì)存在扭曲(distortion)現(xiàn)象,出現(xiàn)不匹配點(diǎn),Vinh 等[12,13]針對此問題進(jìn)行了改進(jìn),并增加一個(gè)后處理步驟,使分類器更加精確.Vinh 等[14]還提出了一種基于約束的自訓(xùn)練算法,與正類集合最近的實(shí)例t,必須滿足約束DL(t|H)<DL(t),才能添加到正類集合.另外,Vinh 等還定義了安全距離(safe distance),當(dāng)實(shí)例與正類集合之間的距離小于或等于安全距離,則將該實(shí)例放入正類集合中.

    目前絕大多數(shù)研究工作集中在單變量時(shí)間序列半監(jiān)督分類算法性能的提高,以及停止準(zhǔn)則的改進(jìn)方面,對MTS 半監(jiān)督分類的研究很少.在對MTS 進(jìn)行半監(jiān)督分類時(shí),主要存在變量之間的復(fù)雜相關(guān)關(guān)系以及樣本長度不一致等因素,使得標(biāo)準(zhǔn)分類器很難直接使用.Li 等[15,16]提出了兩種基于標(biāo)準(zhǔn)SVD 的特征提取方法(以下簡稱Li’s first、Li’s second 方法)用于MTS 分類,Li’s first 方法是將第1 個(gè)奇異向量u1與經(jīng)過標(biāo)準(zhǔn)化后由奇異值組成的向量σnormalized相連,作為MTS 樣本的特征表示.Li’s second 方法將加權(quán)以后的第1 奇異向量w1u1與加權(quán)后的第2 奇異向量w2u2相連,作為MTS 樣本的特征表示.這兩種方法本質(zhì)上屬于一維奇異值分解,但是MTS 包含變量維和時(shí)間維兩個(gè)維度,本文提出基于2DSVD 的半監(jiān)督分類方法,從行和列兩個(gè)方向?qū)TS 樣本進(jìn)行降維,清晰地考慮了MTS 樣本的二維特性.

    2 基于2DSVD 的MTS 半監(jiān)督分類算法

    2.1 訓(xùn)練分類器

    本文提出的MTS 半監(jiān)督分類算法主要包括4 個(gè)步驟:第一步,使用未標(biāo)記數(shù)據(jù)集U來計(jì)算變換矩陣Ur以及Vs,獲取每個(gè)訓(xùn)練樣本的特征矩陣;第二步,隨機(jī)選取若干個(gè)正類樣本的特征矩陣作為初始標(biāo)記數(shù)據(jù)P;第三步,計(jì)算集合U中每個(gè)樣本到集合P的歐氏距離,將集合U中與集合P最近的樣本,從集合U中刪除,添加至集合P;第四步,重復(fù)第三步,直到滿足停止標(biāo)準(zhǔn)為止.

    基于2DSVD 的MTS 半監(jiān)督分類算法如算法1 所示.在步驟7 中,本文采用Wei 等[8]提出的停止標(biāo)準(zhǔn),即在迭代過程中,當(dāng)正類樣本的最小最近鄰距離在趨于穩(wěn)定后的第一次顯著下降時(shí),即停止.TWOSVDSSC分為兩個(gè)階段,步驟1-步驟5 為降維階段:設(shè)未標(biāo)記數(shù)據(jù)集U中有M個(gè)MTS 樣本,算法的行-行協(xié)方差矩陣F為m×m矩陣,列-列協(xié)方差矩陣G為n×n矩陣[5],由于對n×n矩陣進(jìn)行奇異值分解的時(shí)間復(fù)雜度為O(n3)[2],所以算法中步驟1-步驟4 的時(shí)間復(fù)雜度為O(m3+n3);步驟5 是計(jì)算未標(biāo)記數(shù)據(jù)集U中每一個(gè)MTS 樣本的特征矩陣,時(shí)間復(fù)雜度為O(M*r*s),由于在MTS 樣本中,變量個(gè)數(shù)n以及參數(shù)r和s往往都遠(yuǎn)小于樣本長度m,因此步驟1-步驟5 的時(shí)間復(fù)雜度主要取決于樣本長度;步驟6-步驟8 為訓(xùn)練分類器階段,時(shí)間復(fù)雜度為O(M2).所以算法的復(fù)雜度為O((m3+n3)+(M*r*s)+M2).

    分類器訓(xùn)練好之后,在使用分類器對待測樣本進(jìn)行分類時(shí),如果待測樣本與任何一個(gè)標(biāo)記為正類樣本之間的距離小于閾值r,則該樣本分類為正類,否則為負(fù)類[8],閾值r為正類樣本與其最近鄰之間距離的平均值.

    算法1.基于2DSVD 的MTS 半監(jiān)督分類算法輸入:P 是初始訓(xùn)練集,包含少量已標(biāo)記正類樣本;U 是未標(biāo)記數(shù)據(jù)集;nSeeds 是初始標(biāo)記為正類樣本的個(gè)數(shù).輸出:訓(xùn)練好的分類器.1.計(jì)算U 中行-行協(xié)方差矩陣F;2.使用SVD 計(jì)算F 的特征向量,由F 的前r 個(gè)主要特征向量組成的變換矩陣Ur;3.計(jì)算U 中列-列協(xié)方差矩陣G;4.使用SVD 計(jì)算G 的特征向量,由G 的前s 個(gè)主要特征向量組成的變換矩陣Vs;5.計(jì)算U 中每個(gè)MTS 樣本的特征矩陣Mi;6.隨機(jī)選取nSeeds 個(gè)正類樣本放入集合P;7.計(jì)算集合U 中每個(gè)樣本到集合P 的歐氏距離,將集合U 中與集合P 最近的樣本,從集合U 中刪除,添加至集合P;8.重復(fù)步驟7,直到滿足停止標(biāo)準(zhǔn)為止.

    2.2 評估分類器

    算法1 僅包含來自U中的正類樣本,屬于一類分類器.本文采用測試集對分類器的性能進(jìn)行測試,測試集中包含一些正類樣本和其他類樣本.采用經(jīng)典的精確度(Precision)和召回率(Recall)來衡量分類器的性能.在本文中,精確度的值等于召回率的值,即假的負(fù)類(False negatives) 數(shù)量與假的正類(False positives)數(shù)量相同.精確度的定義如下所示[3],其中K是指測試集中的正類樣本的個(gè)數(shù),Npositive為在前K個(gè)最接近P集合的樣本中,正類樣本的個(gè)數(shù).

    3 實(shí)驗(yàn)

    3.1 數(shù)據(jù)集描述

    本文實(shí)驗(yàn)數(shù)采用的Lp1、Lp2、Lp4、Lp5 數(shù)據(jù)集[17]包含機(jī)器人在故障檢測后的力和扭矩測量值.每個(gè)故障的特征是在故障檢測后每隔一段時(shí)間收集的15 個(gè)力/扭矩樣本,Lp1、Lp2、Lp4、Lp5 數(shù)據(jù)集中每個(gè)樣本包含6 個(gè)變量;BCI 數(shù)據(jù)集[18,19]中MTS 樣本分為兩種類型:一種是被測試者用左手手指按計(jì)算機(jī)鍵盤時(shí)的腦電圖(EEG)情況,有208 個(gè)樣本;另一種是被測試者用右手手指按計(jì)算機(jī)鍵盤時(shí)的腦電圖情況,也有208 個(gè)樣本.數(shù)據(jù)集中每個(gè)樣本包含28 個(gè)變量;Japanese Vowels 數(shù)據(jù)集[20]記錄9 個(gè)男性在發(fā)日語的元音/ae/,這9 個(gè)男性對應(yīng)的樣本個(gè)數(shù)分別為:61,65,118,74,59,54,70,80 以及59,數(shù)據(jù)集中每個(gè)樣本包含12 個(gè)變量;Wafer 數(shù)據(jù)集[21]記錄真空室傳感器監(jiān)控半導(dǎo)體微電子的制造過程,每一個(gè)硅晶片的生產(chǎn)過程可以用含有6 個(gè)變量的MTS 樣本來描述,并被分為正?;虍惓深悾瑪?shù)據(jù)集中包含327 個(gè)MTS 樣本并被分為2 類:其中正常樣本有200 個(gè),異常樣本有127;AUstralian Sign LANguage(以下簡稱AUSLAN)數(shù)據(jù)集[20]由隨機(jī)選取25 種手勢的MTS 樣本(總共675 個(gè)MTS 樣本)組成,每個(gè)樣本包含22 個(gè)變量;Character Trajectories 數(shù)據(jù)集[22]中所有樣本來自同一位作者,通過書寫單個(gè)字符來記錄筆尖(pen tip)軌跡,記錄時(shí)只考慮帶有單一落筆段的字符,每個(gè)樣本包含x 和y 坐標(biāo)以及筆尖力度這3 個(gè)變量;Gas sensors 數(shù)據(jù)集[23,24]包含由MOX 以及溫度和濕度這三種傳感器組成的氣體傳感器,記錄來自3 種不同氣體所產(chǎn)生的觀測值,數(shù)據(jù)集中每個(gè)樣本包含10 個(gè)變量.表1列出了10 個(gè)MTS 數(shù)據(jù)集的主要特征.2DSVD 要求數(shù)據(jù)集中所有MTS 樣本具有相同長度.對于具有不同長度樣本的MTS 數(shù)據(jù)集,本文采用Rodriguez 等[25]提出的方法,將所有MTS 樣本的長度都延長到該數(shù)據(jù)集中最長MTS 樣本的長度.延長方法如下:如將長度為100 的MTS 樣本延長至120,只需將樣本中每5 個(gè)值中的一個(gè)值復(fù)制即可.該方法使得原樣本中的所有值都保留在延長后的樣本中,不會(huì)損失任何數(shù)據(jù)信息.

    表1 數(shù)據(jù)集描述

    3.2 性能比較

    將本文提出的基于2DSVD 的MTS 特征提取方法,與基于擴(kuò)展Frobenius 范數(shù)的距離DEros[26]、中心序列[27]、以及基于一維SVD 的Li’s first,Li’s second 方法[15,16]分類性能進(jìn)行比較.在實(shí)驗(yàn)中,將數(shù)據(jù)集中類別標(biāo)記為1(class label=1)的樣本選為正類樣本數(shù)據(jù),其它類樣本皆為負(fù)類樣本數(shù)據(jù).在算法2.1 中,初始正類樣本的個(gè)數(shù)nSeeds分別取1、3、5 個(gè),實(shí)驗(yàn)重復(fù)100次,表2、3、4 給出了各種方法100 次實(shí)驗(yàn)的平均Precision.

    表2、表3、表4給出了在10 個(gè)數(shù)據(jù)集上使用不同方法進(jìn)行半監(jiān)督分類的Precision.表中列2 和列3 給出了在數(shù)據(jù)集上使用基于擴(kuò)展Frobenius 范數(shù)的距離DEros[26]以及中心序列[27]的方法進(jìn)行分類的Precision;表中列4 和列5 給出了在數(shù)據(jù)集上使用Li’s first 以及Li’s fecond 方法進(jìn)行分類的Precision;列6 給出了使用2DSVD 進(jìn)行分類時(shí)最高的Precision以及相應(yīng)參數(shù)r和s的值,其中,r和s分別表示使用2DSVD 方法得到對應(yīng)特征矩陣的行及列的個(gè)數(shù).

    從表2可以看出,當(dāng)初始正類樣本的個(gè)數(shù)nSeeds為1 時(shí),2DSVD 在10 個(gè)MTS 數(shù)據(jù)集上分類的平均Precision 為 0.76,DEros的平均值為0.39,中心序列的平均值為0.63,Li’s First 以及Li’s Second 的平均值分別為 0.53 和0.52;從表5中可以看到,2DSVD 與其它4 種方法的Wilcoxon 符號秩檢驗(yàn)的概率p值都小于0.05,說明2DSVD 的分類性能顯著地好于其它四種方法.當(dāng)nSeeds 為3 或5 時(shí),也可以得到相同的結(jié)論.從表2、表3、表4中還可以看出,各種方法的平均Precision隨著nSeeds增大而增大,說明增加初始正類樣本個(gè)數(shù),能夠提高算法的分類性能.

    表2 nSeeds=1 時(shí)各種方法的Precision

    表3 nSeeds=3 時(shí)各種方法的Precision

    表4 nSeeds=5 時(shí)各種方法的Precision

    3.3 參數(shù)對半監(jiān)督分類性能的影響

    本文提出的分類算法有兩個(gè)參數(shù):一個(gè)是行-行協(xié)方差矩陣的主要特征向量個(gè)數(shù)r,另一個(gè)是列-列協(xié)方差矩陣的主要特征向量個(gè)數(shù)s.圖1、圖2分別給出了在AUSLAN、Vowel 數(shù)據(jù)集上,將參數(shù)r固定為1,Precision隨參數(shù)s的變化情況.從圖1和圖2可以看出,當(dāng)s=1 時(shí),Precision最??;隨著s 逐漸增加,算法的Precision快速上升,然后趨于平穩(wěn);所以,在算法的執(zhí)行過程中,可以選取較大的s值來提高分類的Precision.

    表5 Wilcoxon 符號秩檢驗(yàn)

    圖1 AUSLAN 數(shù)據(jù)集Precision 隨列-列協(xié)方差矩陣的主要特征向量個(gè)數(shù)s 的變化

    圖2 Vowel 數(shù)據(jù)集Precision 隨列-列協(xié)方差矩陣的主要特征向量個(gè)數(shù)s 的變化

    圖3給出了在AUSLAN 數(shù)據(jù)集上,將參數(shù)s固定為21,Precision隨參數(shù)r的變化情況.圖4給出了在Vowel 數(shù)據(jù)集上,將參數(shù)s固定為12,Precision隨參數(shù)r的變化情況.從圖3和圖4可以看出,當(dāng)參數(shù)r增加時(shí),分類的Precision趨于平穩(wěn);所以,在算法執(zhí)行過程中,可以選取適當(dāng)?shù)膔值即可.

    圖3 AUSLAN 數(shù)據(jù)集Precision 隨行-行協(xié)方差矩陣的主要特征向量個(gè)數(shù)r 的變化

    圖4 Vowel 數(shù)據(jù)集Precision 隨行-行協(xié)方差矩陣的主要特征向量個(gè)數(shù)r 的變化

    在本文實(shí)驗(yàn)中,參數(shù)r和s的選取方法如下[2]:首先選擇一個(gè)較大的s值,使得這s個(gè)列-列協(xié)方差矩陣的主要特征向量能夠描述列-列之間總變異(total column-column variations)的98%或99%,其次,讓r值從1 增加到m,其中m為觀測值個(gè)數(shù),計(jì)算相對于每一個(gè)r值的所有訓(xùn)練樣本的重構(gòu)誤差平方和,最后根據(jù)重構(gòu)誤差平方和的相對變化情況選取適當(dāng)?shù)膮?shù)r.

    4 結(jié)論與展望

    本文提出了一種基于2DSVD 的MTS 半監(jiān)督分類方法,在10 個(gè)MTS 數(shù)據(jù)集上對該方法進(jìn)行驗(yàn)證,實(shí)驗(yàn)結(jié)果表明,本文提出的算法顯著地好于基于一維SVD的Li’s First、Li’s Second 方法[15,16],基于擴(kuò)展Frobenius范數(shù)的距離DEros[26],以及中心序列[27].雖然本文建立的是一類分類器,因此也可以很容易地修改本文提出的算法以適應(yīng)多類問題.本文提出的算法有兩個(gè)參數(shù)r和s,如何自動(dòng)地選擇最優(yōu)的r和s值以及選取更優(yōu)的分類器和停止標(biāo)準(zhǔn)值得今后進(jìn)一步研究.

    猜你喜歡
    協(xié)方差分類器分類
    分類算一算
    分類討論求坐標(biāo)
    數(shù)據(jù)分析中的分類討論
    BP-GA光照分類器在車道線識(shí)別中的應(yīng)用
    電子測試(2018年1期)2018-04-18 11:52:35
    教你一招:數(shù)的分類
    加權(quán)空-譜與最近鄰分類器相結(jié)合的高光譜圖像分類
    結(jié)合模糊(C+P)均值聚類和SP-V-支持向量機(jī)的TSK分類器
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    基于LLE降維和BP_Adaboost分類器的GIS局部放電模式識(shí)別
    爱豆传媒免费全集在线观看| 国产成人午夜福利电影在线观看| 我要看黄色一级片免费的| 又大又黄又爽视频免费| 午夜福利网站1000一区二区三区| 久久久久久久久久成人| 成人综合一区亚洲| 国产精品久久久久久久久免| 少妇人妻精品综合一区二区| 男男h啪啪无遮挡| 亚洲欧美成人综合另类久久久| 国产精品久久久久久久电影| 午夜免费男女啪啪视频观看| 99热这里只有精品一区| 美女国产高潮福利片在线看| 亚洲国产精品一区三区| 中文精品一卡2卡3卡4更新| 交换朋友夫妻互换小说| 伦精品一区二区三区| 丰满饥渴人妻一区二区三| 丰满乱子伦码专区| 亚洲一区二区三区欧美精品| 国产成人freesex在线| 亚洲精品乱码久久久v下载方式| 日本欧美国产在线视频| 18禁在线播放成人免费| 亚洲欧洲国产日韩| 永久免费av网站大全| 高清视频免费观看一区二区| 亚洲精品一区蜜桃| 简卡轻食公司| 中文字幕av电影在线播放| 精品久久蜜臀av无| av在线观看视频网站免费| 欧美日韩综合久久久久久| 国产午夜精品一二区理论片| 老女人水多毛片| 少妇熟女欧美另类| 一区二区三区四区激情视频| 一级片'在线观看视频| 亚洲美女视频黄频| 日产精品乱码卡一卡2卡三| 国产女主播在线喷水免费视频网站| 国产欧美亚洲国产| 女性生殖器流出的白浆| 三级国产精品片| 男男h啪啪无遮挡| 欧美少妇被猛烈插入视频| 亚洲欧洲国产日韩| 亚洲精品乱码久久久久久按摩| videosex国产| 国产精品人妻久久久久久| 国产深夜福利视频在线观看| 久久久久精品性色| 最近中文字幕高清免费大全6| 亚洲精品色激情综合| 亚洲色图 男人天堂 中文字幕 | 国产亚洲av片在线观看秒播厂| 秋霞伦理黄片| 日韩,欧美,国产一区二区三区| 我的老师免费观看完整版| 国产极品粉嫩免费观看在线 | 免费黄色在线免费观看| 久久久久久久久久久丰满| 日韩成人伦理影院| 国产精品.久久久| 精品人妻在线不人妻| 国产精品人妻久久久久久| 美女中出高潮动态图| 一级毛片aaaaaa免费看小| 免费少妇av软件| 亚洲三级黄色毛片| 亚洲欧美成人精品一区二区| 精品亚洲乱码少妇综合久久| 少妇熟女欧美另类| 国产精品久久久久成人av| 一区二区三区四区激情视频| 亚洲欧洲精品一区二区精品久久久 | 国内精品宾馆在线| 国模一区二区三区四区视频| 午夜福利视频精品| 男男h啪啪无遮挡| 国产黄片视频在线免费观看| 18禁在线播放成人免费| 日本av手机在线免费观看| 国产亚洲午夜精品一区二区久久| 十分钟在线观看高清视频www| 亚洲怡红院男人天堂| 精品久久久久久电影网| 久久久午夜欧美精品| 中文字幕人妻丝袜制服| 黄色配什么色好看| 考比视频在线观看| 国产免费又黄又爽又色| 免费黄网站久久成人精品| 99热6这里只有精品| 老司机影院成人| 大片免费播放器 马上看| 黑人欧美特级aaaaaa片| 一级片'在线观看视频| 亚洲av.av天堂| 最黄视频免费看| 午夜日本视频在线| 插逼视频在线观看| 国产成人freesex在线| 国产av一区二区精品久久| av黄色大香蕉| 久久久久久久大尺度免费视频| 亚洲人成网站在线播| 99热6这里只有精品| 午夜激情av网站| 18禁在线无遮挡免费观看视频| 国产在视频线精品| 日本av手机在线免费观看| 国产男女内射视频| 我的女老师完整版在线观看| 国产 一区精品| 亚洲精品自拍成人| 人成视频在线观看免费观看| 午夜福利在线观看免费完整高清在| 在线天堂最新版资源| 国产有黄有色有爽视频| 熟女av电影| 国产不卡av网站在线观看| 国产精品不卡视频一区二区| 午夜免费观看性视频| 国产伦理片在线播放av一区| 久久久久久久亚洲中文字幕| 久热久热在线精品观看| 99热全是精品| 老司机影院成人| 国模一区二区三区四区视频| 在线观看美女被高潮喷水网站| 观看av在线不卡| 国产av精品麻豆| 3wmmmm亚洲av在线观看| 亚洲精品,欧美精品| 亚洲av国产av综合av卡| 91午夜精品亚洲一区二区三区| 成人毛片60女人毛片免费| av.在线天堂| 欧美日韩一区二区视频在线观看视频在线| 午夜福利视频精品| 人人澡人人妻人| 免费观看性生交大片5| 亚洲综合精品二区| 在线免费观看不下载黄p国产| 中文字幕亚洲精品专区| 午夜影院在线不卡| 91成人精品电影| 亚洲av成人精品一区久久| 精品一区二区三卡| 少妇被粗大的猛进出69影院 | 久久人人爽av亚洲精品天堂| 人人妻人人添人人爽欧美一区卜| av国产精品久久久久影院| 一区二区三区四区激情视频| 美女内射精品一级片tv| 久久国内精品自在自线图片| 丝袜脚勾引网站| 美女视频免费永久观看网站| 日日啪夜夜爽| 夜夜爽夜夜爽视频| 久久毛片免费看一区二区三区| 97超碰精品成人国产| 久久av网站| videosex国产| 国产精品 国内视频| 男人操女人黄网站| 九九久久精品国产亚洲av麻豆| 精品国产一区二区久久| 一区二区av电影网| www.色视频.com| 精品久久久久久电影网| 亚洲第一区二区三区不卡| 看十八女毛片水多多多| 看免费成人av毛片| 十八禁高潮呻吟视频| 天堂俺去俺来也www色官网| 一本大道久久a久久精品| 一本一本综合久久| 老司机影院成人| 黄色视频在线播放观看不卡| 51国产日韩欧美| 韩国高清视频一区二区三区| 欧美丝袜亚洲另类| 中文字幕亚洲精品专区| 18禁动态无遮挡网站| 午夜精品国产一区二区电影| 内地一区二区视频在线| av黄色大香蕉| 日韩视频在线欧美| av视频免费观看在线观看| 中国美白少妇内射xxxbb| 女性生殖器流出的白浆| 精品国产露脸久久av麻豆| 亚洲av欧美aⅴ国产| 国产极品天堂在线| 女人精品久久久久毛片| 久久久久久久久大av| 极品人妻少妇av视频| av又黄又爽大尺度在线免费看| 亚洲一级一片aⅴ在线观看| 激情五月婷婷亚洲| 在线观看美女被高潮喷水网站| 精品卡一卡二卡四卡免费| 丝袜脚勾引网站| 嘟嘟电影网在线观看| 亚洲精品乱码久久久v下载方式| 午夜福利,免费看| 欧美激情国产日韩精品一区| 全区人妻精品视频| 成人毛片a级毛片在线播放| 亚洲精品日韩在线中文字幕| 欧美日韩视频精品一区| 最新中文字幕久久久久| 国产欧美另类精品又又久久亚洲欧美| 午夜91福利影院| 熟妇人妻不卡中文字幕| 人人妻人人澡人人爽人人夜夜| 色94色欧美一区二区| 岛国毛片在线播放| 日本免费在线观看一区| 天天躁夜夜躁狠狠久久av| 亚洲色图 男人天堂 中文字幕 | 午夜免费观看性视频| 99re6热这里在线精品视频| 国产在线免费精品| 欧美精品一区二区大全| 9色porny在线观看| 国产69精品久久久久777片| 十八禁网站网址无遮挡| 国产熟女午夜一区二区三区 | 日本黄色日本黄色录像| 国产男女超爽视频在线观看| 纵有疾风起免费观看全集完整版| 亚洲精品久久久久久婷婷小说| 日韩精品有码人妻一区| 啦啦啦在线观看免费高清www| 丰满迷人的少妇在线观看| 精品卡一卡二卡四卡免费| 欧美日韩av久久| 欧美成人精品欧美一级黄| tube8黄色片| 欧美日本中文国产一区发布| 纵有疾风起免费观看全集完整版| 最新中文字幕久久久久| 青青草视频在线视频观看| 欧美精品亚洲一区二区| 久久久久久久亚洲中文字幕| 黄色一级大片看看| 18在线观看网站| 日韩大片免费观看网站| 久久久久久久精品精品| 日韩一区二区三区影片| 18禁观看日本| 青春草视频在线免费观看| a级片在线免费高清观看视频| www.av在线官网国产| 三级国产精品欧美在线观看| 国产精品一国产av| 一本色道久久久久久精品综合| 久久久久人妻精品一区果冻| 一二三四中文在线观看免费高清| 一级毛片电影观看| 国产免费视频播放在线视频| 简卡轻食公司| 国产深夜福利视频在线观看| 久久久欧美国产精品| 成人亚洲精品一区在线观看| 精品人妻偷拍中文字幕| 免费播放大片免费观看视频在线观看| 午夜激情久久久久久久| 91精品国产国语对白视频| av女优亚洲男人天堂| 香蕉精品网在线| 99久久精品一区二区三区| 大香蕉久久网| 亚洲精品亚洲一区二区| 国产一区二区在线观看日韩| 久久午夜综合久久蜜桃| 午夜91福利影院| 中文字幕免费在线视频6| 精品人妻熟女毛片av久久网站| 超色免费av| 国产精品国产三级国产av玫瑰| 18禁在线播放成人免费| 久久久久网色| 亚洲图色成人| 十八禁网站网址无遮挡| 91午夜精品亚洲一区二区三区| 各种免费的搞黄视频| 久久韩国三级中文字幕| 国产69精品久久久久777片| 91精品一卡2卡3卡4卡| 欧美精品一区二区大全| 国产成人精品婷婷| 99久国产av精品国产电影| av在线老鸭窝| 岛国毛片在线播放| 如日韩欧美国产精品一区二区三区 | 欧美精品一区二区大全| 亚洲精品乱久久久久久| 最后的刺客免费高清国语| 亚洲精品一区蜜桃| 日韩熟女老妇一区二区性免费视频| 最新的欧美精品一区二区| 亚洲精品色激情综合| 中文字幕av电影在线播放| 久久精品国产亚洲网站| 在线观看一区二区三区激情| 性高湖久久久久久久久免费观看| 最新的欧美精品一区二区| a级毛片免费高清观看在线播放| 99久久综合免费| 亚洲美女搞黄在线观看| 久久久国产精品麻豆| 草草在线视频免费看| 最后的刺客免费高清国语| 久久人人爽人人片av| 国产淫语在线视频| 99国产综合亚洲精品| 成人国产麻豆网| 性色avwww在线观看| 国产深夜福利视频在线观看| 欧美激情国产日韩精品一区| 视频在线观看一区二区三区| 国产精品人妻久久久影院| 精品卡一卡二卡四卡免费| 香蕉精品网在线| 人人妻人人爽人人添夜夜欢视频| 男女啪啪激烈高潮av片| 成人影院久久| 在线观看www视频免费| 亚洲性久久影院| 国产精品 国内视频| 国产免费福利视频在线观看| 中文欧美无线码| 日韩强制内射视频| 亚洲精品日本国产第一区| 丝袜喷水一区| 色网站视频免费| 亚洲美女视频黄频| 日日爽夜夜爽网站| 欧美3d第一页| 国产伦理片在线播放av一区| 久久久久久久久大av| 在线 av 中文字幕| 嫩草影院入口| 久久久久久久国产电影| 你懂的网址亚洲精品在线观看| 久久影院123| 免费看av在线观看网站| 丁香六月天网| 高清午夜精品一区二区三区| 最新中文字幕久久久久| 亚洲国产精品一区三区| 韩国av在线不卡| 99久久精品一区二区三区| 寂寞人妻少妇视频99o| 99久久中文字幕三级久久日本| 在线观看美女被高潮喷水网站| 青春草国产在线视频| 国产成人freesex在线| 日韩精品有码人妻一区| 亚洲精品av麻豆狂野| 老熟女久久久| 亚洲av日韩在线播放| 精品久久蜜臀av无| 最后的刺客免费高清国语| 久久国产精品男人的天堂亚洲 | 日韩一区二区三区影片| 亚洲性久久影院| 2022亚洲国产成人精品| 久久久精品94久久精品| 久久人人爽人人片av| 99久久精品一区二区三区| 国产女主播在线喷水免费视频网站| 一本色道久久久久久精品综合| 最后的刺客免费高清国语| 老熟女久久久| 国产精品一区二区在线观看99| 免费播放大片免费观看视频在线观看| av国产精品久久久久影院| 嘟嘟电影网在线观看| 一个人看视频在线观看www免费| 大香蕉久久成人网| 欧美性感艳星| 999精品在线视频| 在线精品无人区一区二区三| 欧美日韩在线观看h| 欧美性感艳星| 久久久精品区二区三区| 国产一区有黄有色的免费视频| 男女免费视频国产| 国产精品欧美亚洲77777| 在现免费观看毛片| 午夜影院在线不卡| 久久久久久久久久久久大奶| 青春草亚洲视频在线观看| 自线自在国产av| .国产精品久久| 日韩不卡一区二区三区视频在线| 午夜91福利影院| 男人添女人高潮全过程视频| 日韩人妻高清精品专区| 91久久精品国产一区二区三区| 人妻系列 视频| 不卡视频在线观看欧美| 亚洲欧洲国产日韩| 亚洲综合精品二区| 亚洲第一av免费看| 91午夜精品亚洲一区二区三区| 91久久精品国产一区二区成人| 又黄又爽又刺激的免费视频.| 母亲3免费完整高清在线观看 | 热99国产精品久久久久久7| av视频免费观看在线观看| 哪个播放器可以免费观看大片| 啦啦啦视频在线资源免费观看| av黄色大香蕉| 精品人妻偷拍中文字幕| av国产精品久久久久影院| 欧美精品高潮呻吟av久久| 免费人妻精品一区二区三区视频| 观看美女的网站| 狠狠精品人妻久久久久久综合| 五月开心婷婷网| 97超视频在线观看视频| 亚洲国产精品成人久久小说| 国产成人午夜福利电影在线观看| 男人操女人黄网站| 国产成人a∨麻豆精品| 久久久久久人妻| 欧美日韩综合久久久久久| 久久久国产一区二区| 国产精品一二三区在线看| 九草在线视频观看| 日日撸夜夜添| 最近的中文字幕免费完整| 亚洲国产精品999| 精品国产乱码久久久久久小说| 最后的刺客免费高清国语| 一区二区日韩欧美中文字幕 | videos熟女内射| 日韩av在线免费看完整版不卡| 天堂8中文在线网| 久久精品国产亚洲av天美| 成人亚洲精品一区在线观看| 黄色配什么色好看| 国产色爽女视频免费观看| 免费高清在线观看视频在线观看| 啦啦啦视频在线资源免费观看| 日韩欧美一区视频在线观看| 人人妻人人澡人人看| 成年人午夜在线观看视频| 人人妻人人爽人人添夜夜欢视频| 少妇被粗大猛烈的视频| 亚洲成人手机| 男人添女人高潮全过程视频| 欧美日韩成人在线一区二区| 国产高清有码在线观看视频| 精品亚洲成a人片在线观看| 中文天堂在线官网| 免费观看性生交大片5| 中文字幕人妻熟人妻熟丝袜美| 精品人妻偷拍中文字幕| 高清不卡的av网站| 久久国产精品男人的天堂亚洲 | 熟女电影av网| 成人漫画全彩无遮挡| 国产色爽女视频免费观看| 国产精品偷伦视频观看了| 能在线免费看毛片的网站| 一级毛片aaaaaa免费看小| 国产极品粉嫩免费观看在线 | 欧美成人精品欧美一级黄| 老司机影院毛片| 国产视频首页在线观看| 中文字幕制服av| 美女大奶头黄色视频| 亚洲精品久久午夜乱码| 如何舔出高潮| 寂寞人妻少妇视频99o| 天美传媒精品一区二区| 久久久国产欧美日韩av| 汤姆久久久久久久影院中文字幕| 毛片一级片免费看久久久久| av视频免费观看在线观看| 精品人妻一区二区三区麻豆| 美女大奶头黄色视频| 91精品伊人久久大香线蕉| 一级爰片在线观看| 日本色播在线视频| 国产精品嫩草影院av在线观看| 亚洲精品美女久久av网站| 中文精品一卡2卡3卡4更新| 久久久久久人妻| 亚洲精品自拍成人| 精品久久久精品久久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男女高潮啪啪啪动态图| av线在线观看网站| 最新的欧美精品一区二区| 一本久久精品| 亚洲精品美女久久av网站| 免费久久久久久久精品成人欧美视频 | 国产成人精品久久久久久| 国产日韩一区二区三区精品不卡 | 欧美日本中文国产一区发布| 一级a做视频免费观看| 午夜福利视频在线观看免费| 亚洲精品,欧美精品| 亚洲欧美清纯卡通| 色5月婷婷丁香| 欧美精品国产亚洲| 伊人久久国产一区二区| 免费高清在线观看视频在线观看| 人妻人人澡人人爽人人| 人妻制服诱惑在线中文字幕| 国产片内射在线| tube8黄色片| 亚洲欧洲国产日韩| 久热久热在线精品观看| 国产永久视频网站| 亚洲婷婷狠狠爱综合网| 人人妻人人爽人人添夜夜欢视频| 欧美精品国产亚洲| 丝瓜视频免费看黄片| 国产精品无大码| 简卡轻食公司| 人体艺术视频欧美日本| 色视频在线一区二区三区| 午夜视频国产福利| 中国三级夫妇交换| 特大巨黑吊av在线直播| 中文欧美无线码| 黄色视频在线播放观看不卡| 国产成人精品在线电影| 成人18禁高潮啪啪吃奶动态图 | 一个人看视频在线观看www免费| 亚洲欧洲精品一区二区精品久久久 | 日本色播在线视频| 最近的中文字幕免费完整| 日韩,欧美,国产一区二区三区| 爱豆传媒免费全集在线观看| 内地一区二区视频在线| 国产 一区精品| 青春草国产在线视频| 男的添女的下面高潮视频| 五月天丁香电影| 熟女电影av网| 啦啦啦视频在线资源免费观看| 成人午夜精彩视频在线观看| 久久久久久久久大av| 欧美日韩亚洲高清精品| 国产午夜精品久久久久久一区二区三区| 亚洲精品日本国产第一区| 色网站视频免费| 国产黄片视频在线免费观看| 午夜免费观看性视频| 国产精品人妻久久久影院| 三上悠亚av全集在线观看| 午夜av观看不卡| 少妇熟女欧美另类| av网站免费在线观看视频| 日本欧美国产在线视频| 国产男女超爽视频在线观看| 亚洲精品av麻豆狂野| 色5月婷婷丁香| 制服人妻中文乱码| 成人午夜精彩视频在线观看| 亚州av有码| 天天影视国产精品| 男女边吃奶边做爰视频| 国产淫语在线视频| 一级毛片电影观看| 天堂俺去俺来也www色官网| 日韩人妻高清精品专区| 五月天丁香电影| 久久国产精品大桥未久av| 一二三四中文在线观看免费高清| 十八禁网站网址无遮挡| 久久久久人妻精品一区果冻| 久久久精品免费免费高清| 夫妻午夜视频| 美女中出高潮动态图| 国产一区亚洲一区在线观看| 中文字幕人妻熟人妻熟丝袜美| 少妇人妻久久综合中文| 免费黄网站久久成人精品| 精品久久久久久久久亚洲| 久久久久久人妻| 国产爽快片一区二区三区| 国产熟女午夜一区二区三区 | 满18在线观看网站| 免费观看的影片在线观看| 中文字幕人妻丝袜制服| 高清欧美精品videossex| 最近中文字幕2019免费版| 18禁在线无遮挡免费观看视频| 久久国产精品大桥未久av| 国模一区二区三区四区视频| 亚洲人成77777在线视频| 国语对白做爰xxxⅹ性视频网站| 精品国产一区二区三区久久久樱花| 亚洲国产精品一区三区| 亚洲精品一区蜜桃| 人人妻人人澡人人看| 美女视频免费永久观看网站| 精品国产乱码久久久久久小说| 三上悠亚av全集在线观看| 日韩不卡一区二区三区视频在线| 免费黄网站久久成人精品| 视频区图区小说|