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

基于加速度傳感器的中文簽名身份認證

2017-06-27 08:10:42鄭建彬詹恩奇
計算機應用 2017年4期
關鍵詞:用戶系統

劉 威,汪 陽,鄭建彬,詹恩奇

1.武漢理工大學 信息工程學院,武漢 430070; 2.光纖傳感技術與信息處理教育部重點實驗室(武漢理工大學),武漢 430070)(*通信作者電子郵箱1530379001@qq.com)

基于加速度傳感器的中文簽名身份認證

劉 威1,2*,汪 陽1,2,鄭建彬1,2,詹恩奇1,2

1.武漢理工大學 信息工程學院,武漢 430070; 2.光纖傳感技術與信息處理教育部重點實驗室(武漢理工大學),武漢 430070)(*通信作者電子郵箱1530379001@qq.com)

采集用戶在簽名過程中的三軸加速度信息,可用來實現身份認證。中文簽名結構較為復雜,在空中書寫的過程難以被模仿,但同時也會使同一用戶不同次簽名間的差異增大,提高認證難度。傳統的二維簽名或三維手勢認證方法并不能解決這一問題。為了提高中文空中簽名身份認證效果,改進了全局序列對齊(GSA)算法,對匹配后的序列進行插值操作。不同于傳統GSA算法通過最終匹配分數反映序列間相似度,引入兩種距離指標(歐氏距離和絕對值距離)計算序列間的差異。實驗結果表明,基于距離指標的兩種GSA算法均能提高認證精度,與傳統算法相比,二者的系統等誤率(EER)分別降低了37.6%和52.6%。

生物特征;空中簽名;身份認證;序列對齊;加速度傳感器

0 引言

基于人體生物特征的身份認證技術正在成為信息安全領域研究的熱點。生物特征可以劃分為生理生物特征和行為生物特征兩大類[1-2]。基于生理特征的身份認證技術主要利用指紋[3]、人臉[4]、虹膜[5]等生理特征實現身份認證。這些生理特征因其普適性和唯一性得到了普遍應用。基于生物行為特征的身份認證方法,如聲音識別[6]、步態行走識別[7]等在近年也受到了廣泛關注,手寫簽名[8-9]更是長期以來最為人們所接受的行為特征。但具體到應用中,以上生物特征認證方法仍然存在技術成本高、認證效率低、易仿造等問題。本文提出了一種基于加速度傳感器的中文空中簽名認證方法。利用帶有三軸加速度計的電路板,采集用戶在空中書寫簽名時的加速度信息。通過待測簽名與模板簽名的匹配比對,判斷用戶身份真偽。

模式匹配中對序列間差異的度量主要有兩種方案:一種是在動態規劃過程中尋找兩個序列的最大相似性,表現為在匹配過程中構建分數矩陣,尋找最大匹配分數。如最大共同子序列(Longest Common Subsequence, LCS)算法[10]、全局序列對齊(Global Sequence Alignment, GSA)算法[11]。另外一種如動態時間規整(Dynamic Time Warping, DTW)算法[12],通過計算匹配后序列間的累計距離,衡量序列之間最小差異。這些方法在基于生物特征的認證技術中都有廣泛應用:文獻[13]利用Lumia1520手機上自帶的加速度傳感器采集動態手勢數據,在LCS算法基礎上,提出了一種限制匹配窗口的近似判等最長共同子序列算法,實現了身份認證。文獻[14]結合DTW算法與改進的LCS算法,設計了一個基于動態手勢的身份認證系統,在智能手機上實現,并提高了認證精度。文獻[15]將用戶使用加速度傳感器簽名的過程錄像,通過模仿錄像中的動作偽造簽名;這種偽造方式增加了認證難度,他們改進DTW算法,綜合利用三維手勢動作的三軸加速度信息,取得較好的認證效果。文獻[16]限制用戶使用特定手勢,使加速度數據變化主要集中在Y軸,并且在預處理中去除重力加速度干擾,之后運用動態規劃思想,計算加速度序列間的差異。

目前,基于加速度傳感器的身份認證研究多集中在采集手勢運動或書寫特定數字過程中的三軸加速度信息,運用算法比對分析,認證用戶身份。針對中文空中簽名身份認證的研究較少。而簽名是人們日常生活中經常使用的行為特征,即使是手持電路板在空中書寫簽名,用戶仍然能夠保持自己簽名時的習慣。相對于手勢動作,中文簽名過程更為復雜,也更難以被模仿。因此,基于加速度傳感器的中文空中簽名身份認證研究是可行的。但中文簽名本身結構的復雜性,也使得同一用戶不同次空中簽名間的差異增大,導致當前的認證算法不能較好地區分真偽簽名。本文在傳統GSA算法的基礎上,在匹配后的兩簽名序列間插入合適的值,找到序列間最佳對齊方式,減小真實簽名間的差異。插值后的兩序列長度相等,序列中的點一一對應。采用絕對值距離衡量簽名數據間的差異,設定合適的閾值,提高身份認證精度。

1 認證系統流程

基于加速度傳感器的中文空中簽名認證流程主要包含五部分:簽名數據采集、預處理、特征提取、匹配和決策認證,如圖1所示。

圖1 認證系統框圖

簽名數據采集是通過自制電路板的相關模塊采集用戶簽名時的三軸加速度信息。三軸加速度傳感器采用STMicroelectronics公司生產的超小型高性能的LSM303DLHC模塊。實驗時,設定加速度計的采集頻率為100 Hz。

數據的預處理是對傳感器采集的原始加速度數據進行平滑濾波,去除可能存在的人為抖動噪聲和設備本身帶來的工頻噪聲,將信號進行規整,盡可能地凸顯真實簽名和偽造簽名之間的差異性,得到一個差異性明顯且有效的簽名信號。本文提取簽名中的每一個采樣點數據,分析匹配。

匹配算法采用基于絕對值距離的GSA算法。模板數據庫中包含模板樣本和模板匹配結果。進入系統時,每個用戶首先需手持加速度數據采集裝置,在空中連續書寫自己的中文簽名3次,完成模板樣本采集。對3組簽名數據兩兩相互匹配,得到3個匹配結果,取其均值,作為模板匹配結果。模板匹配結果用Ym指代,代表一個用戶三次簽名的相似性。

待測簽名需與模板數據庫中的三組模板樣本分別匹配,同樣得到3個匹配結果,取其均值,作為測試結果,并用Yt指代。測試結果與模板匹配結果的比值q可用來反映待測簽名與真實簽名模板之間的相似程度:

q=Yt/Ym

(1)

q將與最終的系統閾值Q比較,判斷待測簽名真偽,實現用戶身份認證。

2 基于絕對值距離的全局對齊算法

2.1GSA算法基本原理

動態規劃,是對一個多階段的過程進行子階段化并逐一求解找到最優解的一種數學方法[17]。全局序列對齊(GSA)算法也是基于動態規劃思想,構建分數矩陣,在統一規則下尋找最佳匹配路徑以確定序列間的對齊方式。在實現的過程中引入獎懲模式,即不匹配時執行懲罰判分,匹配時執行獎勵判分。本文在構建分數矩陣的過程中,針對兩種狀況實行都加分的策略,但賦予了不同權重。具體實現過程如下:

假設一個標準參考信號序列R,長度為M,Rm為其第m幀。另有一待測信號序列T,長度為N,Tn為其第n幀,如式(1)所示:

(2)

定義一個分數矩陣S,其大小為N×M。分數矩陣的填充是一個動態規整過程,通過以下遞推關系計算得到:

Si, j=max{Si, j-1+h,Si-1, j-1+ξ,Si-1, j+h}

(3)

其中:h是常量。ξ是一個復雜的決策函數,遵從高斯分布,表示兩序列對應點之間的相似度:

(4)

其中:σ是常量;(i,j)時刻μ、x的值取決于上一時刻兩序列的值,即μ=Ti-1,x=Ri-1。參數h和參數σ的取值由后續的實驗確定。本文最終設定的參數為h=0.25,σ=0.6;并且,分數矩陣構建的初始條件為:S1,1=0,Si,1=i×h,S1, j=j×h。

在分數矩陣內部,由式(3)的遞推關系,可以得到一條從S1,1到SN,M的最佳匹配路徑。以采集到的同一人的前后兩次簽名為例,圖2的曲線表示兩簽名x軸加速度數據經GSA算法匹配后得到的匹配路徑。

圖2 x軸加速度數據匹配路徑

路徑的終點SN,M的匹配分數指代兩序列間的最大相似度。該分數值越大,說明序列間的相似性越大。在本文的認證系統中,每位用戶的三組模板簽名之間若利用GSA算法兩兩相互匹配,將得到三組匹配分數。取三者均值作為該用戶的訓練結果,記作S(m)。每一個待測簽名在認證過程中均需與三組模板簽名分別匹配,得到三組匹配分數。取其均值,作為測試結果,并記作S(t)。測試結果與訓練結果的比值可以反映待測簽名與模板簽名間的相似度。將該比值與系統閾值S(k)比較,若滿足S(t)/S(m) ≥S(k),則判定測試簽名與模板簽名成功匹配,測試簽名為真實簽名;否則,匹配失敗,測試簽名被認證為偽造簽名。

2.2 改進絕對值距離GSA算法

通過傳統的GSA算法,可以確定簽名加速度序列間的最大相似度。若要進一步量化序列間的差異值,則要利用相關距離指標計算得到。具體操作中,需首先回溯出構建分數矩陣過程中的最佳匹配路徑。

在分數矩陣中,通往Si, j的前一個點只能在Si, j-1、Si-1, j-1、Si-1, j中產生,如圖3所示。

圖3 GSA算法分數矩陣

式(3)中,當Si, j的取值為Si, j-1+h時,分數矩陣中將產生垂直方向的移動,此時,Rj-1和Rj均與Ti匹配。水平方向上的移動則意味著Ti-1和Ti均與Rj匹配。傾斜前進時,序列間的點是一一對應的。將三個方向的移動分別編號,垂直方向移動記為“1”,傾斜方向移動記為“2”,水平方向移動記為“3”。在填充分數矩陣過程中,可以通過記錄對應編號反映匹配路徑的對應軌跡。分數矩陣構建完成后,根據記錄的方向編號可以回溯出一條從SN,M到S1,1的匹配路徑。水平方向和豎直方向的移動會使序列間產生“一對多”或者“多對一”的匹配結果,增大序列間的差異。在匹配后的序列間插入適當零值點,使匹配后序列實現點與點之間的一一對應。插值規則如下:

1)從路徑的終點SN,M開始回溯,保留Tn與Rm;

2)當路徑回溯到某點Si, j時,有可能發生圖3所示的三個方向的移動。發生垂直移動時,則在T序列內點Ti和點Ti+1間插入一個零值點;

3)發生傾斜方向移動時,僅保留原序列內的點,不需插值;

4)發生水平移動時,在R序列內點Rj和點Rj+1之間插入一個零值點。

通過上述步驟完成插值操作后,得到兩個新的序列R′和T′,二者間的點一一對應,且長度均為L′,并滿足如下條件:

max(M,N)≤L′≤M+N

(5)

在新的序列中,用合適的值替換插入的零值點可以進一步減小序列間的差異,得到最佳對齊序列,這一替代過程稱為零值補償。零值點由兩點的平均值進行修正,這兩點分別為:此零值點相鄰的前一個點和該零值點之后最近的一個非零值點。零值補償的表達式如式(6):

(6)

零值補償后產生的新序列記為R"和T",其長度仍均為L′。兩序列間的點仍滿足一一對應,且所有插入的零值點均已得到幅值補償。對于新的序列,可以引入距離指標來量化兩者之間的差異。本文選擇的距離指標是兩序列的整體絕對值距離,用δ(T",R")表示,δ(T",R")值越小,代表序列間的差異值越小,序列越相似。

(7)

應用到認證系統中,當采用改進絕對值GSA算法進行匹配時,匹配結果用距離值表示,代表兩個序列間的最小差異。此時系統的訓練結果和某個測試簽名的測試結果分別用D(m)和D(t)表示,系統的設定閾值為D(k)。與傳統GSA不同,當滿足D(t)/D(m)≤D(k)時,測試簽名與模板簽名匹配成功并被認證為真實簽名。

3 實驗與分析

3.1 實驗數據庫

通過上述簽名采集系統,采集40名用戶的簽名三軸加速度數據,構成實驗數據庫。每名實驗者手持嵌入有加速度計的電路板,并保持加速度計的Z軸垂直指向自己,在空中書寫自己中文簽名共20次,作為真實簽名樣本。每組另需采集20個偽造簽名,從其他用戶中任選三人偽造完成。本文的偽造方式為熟練偽造:采集用戶真實簽名時,用戶需在紙上手寫自己真實簽名作為之后偽造的依據。其他用戶將根據二維手寫簽名記錄模仿偽造空中手寫簽名。

實驗數據庫中共包含1 600個簽名,其中800個真實簽名,800個偽造簽名。具體地,在本文實驗樣本數據庫中,每個實驗者的樣本被細分為模板樣本(3個)、真實測試樣本(17個)、偽造測試樣本(20個)。

3.2 實驗結果

3.2.1 傳統GSA算法實驗

傳統GSA算法通過序列間的匹配分數評價序列間的相似度。為了能更準確地表現序列間的相似度,本文將分數矩陣得到的最終匹配分數除以序列的最終長度L′。數據庫中的每位用戶均包含37組測試樣本,其中17個真實簽名,20個偽造簽名。本實驗中,利用GSA算法將所有測試樣本分別與模板簽名匹配,記錄不同系統閾值下出現的所有錯誤匹配樣本。最終統計40位用戶共1 480組測試樣本的匹配結果,并計算不同閾值下系統的誤拒率(FalseRejectionRate,FRR)、誤納率(FalseAcceptanceRate,FAR)。其中:誤拒率為被認證系統拒絕的真實簽名數與真實簽名總數的比值;誤納率表示被系統接受的偽造簽名數占偽造簽名總數的比例。統計結果如圖4所示。

隨著閾值的增大,認證系統的誤拒率隨之增大,誤納率則逐漸減小。圖4中兩條曲線相交時得到系統的等錯誤率(EqualErrorRate,EER)為10.65%。實際應用中,誤納率是影響系統安全性的最主要因素。實驗應在保持誤拒率在合理范圍內的同時,盡量降低系統誤納率。當設定系統閾值為0.92時,GSA算法取得最佳認證效果。此時,680個真實簽名樣本中被系統錯誤拒絕的個數為79個;800個偽造簽名樣本中,有73個被錯誤接受。對應FRR為11.6%,FAR為9.13%。

3.2.2 基于距離指標的GSA算法實驗

本次實驗是在GSA算法匹配后并完成插值以及零值補償操作的基礎上,分別運用歐氏距離和絕對值距離,量化序列間的差異,記錄并對比兩種方法的最終認證效果。

兩種方法對應的檢測錯誤權衡圖(DetectionErrorTradeoff,DET)曲線如圖5所示。

圖5 兩種改進GSA算法DET曲線圖

直線y=x與各曲線的交點即是該曲線對應算法的EER值。從圖5可以看出,基于絕對值距離的GSA算法的EER值明顯較小,為5.05%。相同FAR條件下,其FRR值也小于基于歐氏距離的GSA算法的對應值。系統閾值的選取過程同傳統GSA算法一樣,盡量降低誤納率以提高系統安全性。兩種算法在各自閾值下最佳認證效果如表1所示。

表1 基于距離指標的兩種GSA算法認證效果

將表1的統計數據和圖4的認證結果對比可知,相對于通過匹配分數判定序列相似度的傳統GSA算法,改進后的基于距離指標的兩種GSA算法認證效果有明顯提升。說明對傳統GSA算法匹配后的序列進行插值和零值補償操作,可以減小真實空中簽名數據間的差異,更好地區分真實和偽造樣本。而采用不同的距離指標來量化序列間的差異度,最終認證效果也不同。選擇絕對值距離時,認證效果最好,FRR為5.29%,FAR為4.50%。

4 結語

本文針對用戶的身份認證問題,提出了一種基于加速度傳感器的空中中文簽名身份認證方法。通過采集和比對用戶在空中書寫中文簽名過程中三軸加速度數據,來確認用戶身份。根據中文空中簽名認證的實際需求,提出了一種基于絕對值距離的GSA算法。該方法在包含有采集自40位用戶,共1 600組簽名的數據集上測試。實驗結果表明本文方法在系統誤拒率、誤納率上都有改進,提高了系統認證精度。

本文的閾值確定過程較為簡單,下一步將繼續對閾值的選擇方法進行研究。此外,在后續的采集過程中,將用戶在空中書寫真實簽名的過程錄像并作為模仿者偽造簽名時的參考,可以增大認證難度,從而更好地測試認證系統的性能。

)

[1]PHILLIPSPJ,MARTINA,WILSONCL,etal.Anintroductionevaluatingbiometricsystems[J].Computer, 2000, 33(2): 56-63.

[2]YAMPOLSKIYRV,GOVINDARAJUV.Behaviouralbiometrics:asurveyandclassification[J].InternationalJournalofBiometrics, 2008, 1(1): 81-113.

[3] 郝玉潔, 馮銀付, 賴攀.基于指紋識別的VPN身份認證研究[J]. 計算機應用, 2009, 29(2):350-352.(HAOYJ,FENGYF,LAIP.VPNidentitycertificationtechnologybasedonfingerprintidentification[J].JournalofComputerApplications, 2009, 29(2): 350-352.)

[4] 張燕昆, 杜平, 劉重慶.基于主元分析與支持向量機的人臉識別方法[J]. 上海交通大學學報, 2002, 36(6):884-886.(ZHANGYK,DUP,LIUCQ.Afacerecognitionmethodbasedonprincipalcomponentanalysisandsupportvectormachine[J].JournalofShanghaiJiaoTongUniversity, 2002, 36(6): 884-886.)

[5]CHOD,PARKKR,RHEEDW,etal.Pupilandirislocalizationforirisrecognitioninmobilephones[C]//ProceedingsoftheSeventhACISInternationalConferenceonSoftwareEngineering,ArtificialIntelligence,Networking,andParallel/DistributedComputing.Piscataway,NJ:IEEE, 2006: 197-201.

[6] 郭武.復雜信道下的說話人識別[D]. 合肥:中國科學技術大學, 2007.(GUOW.Sessionvariabilitytext-independentspeak-recognition[D].Hefei:UniversityofScienceandTechnologyofChina, 2007.)

[7] 高海燕.人體步態識別研究[D]. 北京:北京交通大學, 2010.(GAOHY.Researchonhumangaitrecognition[D].Beijing:BeijingJiaotongUniversity, 2010.)

[8]FENGH,WAHCC.Onlinesignatureverificationusinganewextremepointswarpingtechnique[J].PatternRecognitionLetters, 2003, 24(16): 2943-2951.

[9]KOSTORZI,DOROZR.On-linesignaturerecognitionbasedonreducedsetofpoints[M]//BURDUKR,KURZYNSKIM,WOZNIAKM.ComputerRecognitionSystems4.Berlin:Springer, 2011: 3-11.

[10]BERGROTHL,HAKONENH,RAITAT.Asurveyoflongestcommonsubsequencealgorithms[C]//SPIRE2000:ProceedingsoftheSeventhInternationalSymposiumonStringProcessingInformationRetrieval.Washington,DC:IEEEComputerSociety, 2000: 39-48.

[11] HUANG X. On global sequence alignment [J]. Computer Applications in the Biosciences: Bioinformatics, 1994, 10(3): 227-235.

[12] 李正欣, 張鳳鳴, 李克武.基于DTW的多元時間序列模式匹配方法[J]. 模式識別與人工智能, 2011, 24(3):425-430.(LI Z X, ZHANG F M, LI K W. DTW based pattern matching method for multivariate time series [J]. Pattern Recognition and Artificial Intelligence, 2011, 24(3): 425-430.)

[13] 高煥之, 郭云楠, 劉晴, 等.基于加速度傳感器的WLALCS身份認證及實現[J]. 湖南大學學報(自然科學版), 2014, 41(8):111-116.(GAO H Z, GUO Y N, LIU Q, et al. Accelerometer based authentication method in WLALCS [J]. Journal of Hunan University (Natural Sciences), 2014, 41(8): 111-116.)

[14] 周治平, 苗敏敏.結合DTW與改進STLCS的動態手勢認證研究[J]. 電子測量與儀器學報, 2015, 29(7):1064-1073.(ZHOU Z P, MIAO M M. Research of dynamic gesture authentication based on DTW and improved STLCS [J]. Journal of Electronic Measurement and Instrumentation, 2015, 29(7): 1064-1073.)

[15] CASANOVA J G, AVILA C S, DE SANTOS SIERRA A, et al. A real-time in-air signature biometric technique using a mobile device embedding an accelerometer[M]// Networked Digital Technologies. Berlin: Springer, 2010: 497-503.

[16] OKUMURA F, KUBOTA A, HATORI Y, et al. A study on biometric authentication based on arm sweep action with acceleration sensor[M]// ISPACS 2006: Proceedings of the International Symposium on Intelligent Signal Processing and Communications. Piscataway, NJ: IEEE, 2006: 219-222.

[17] 李端, 錢富才, 李立, 等.動態規劃問題研究[J]. 系統工程理論與實踐, 2007, 27(8):56-64.(LI D, QIAN F C, LI L, et al. Research on dynamic programming [J]. Systems Engineering — Theory & Practice, 2007, 27(8): 56-64.)

This work is partially supported by the National Natural Science Foundation of China (61303028).

LIU Wei, born in 1993, M. S. candidate. His research interests include pattern recognition, image processing.

WANG Yang, born in 1977, Ph. D., associate professor. His research interests include robot control, embedded system.

ZHENG Jianbin, born in 1966, Ph. D., professor. His research interests include pattern recognition, embedded system.

ZHAN Enqi, born in 1972, Ph. D., associate professor. His research interests include signal processing, pattern recognition.

Chinese signature authentication based on accelerometer

LIU Wei1,2*, WANG Yang1,2, ZHENG Jianbin1,2, ZHAN Enqi1,2

(1. College of Information Engineering, Wuhan University of Technology, Wuhan Hubei 430070, China;2. Key Laboratory of Fiber Optic Sensing Technology and Information Processing of Ministry of Education(Wuhan University of Technology), Wuhan Hubei 430070, China)

Acceleration data in 3 axes during a signature process can be collected to authenticate users. Because of complex structures of Chinese signature, the process of signing in the air is hard to be forged, but it also increases differences between signatures performed by the same user which brings more difficulties in authentication. Classical verification methods applied to 2-D signature or hand gesture cannot solve this problem. In order to improve the performance of in-air Chinese signature verification, the classical Global Sequence Alignment (GSA) algorithm was improved, and the interpolation was applied to matching sequences. Different from classical GSA algorithm which uses matching score to measure similarity between sequences, two distance indexes, Euclidean distance and absolute value distance, were introduced to calculate the differences between sequences after interpolation. Experimental results show that both of the two improved GSA algorithms can improve the accuracy of authentication, the Equal Error Rate (EER) of them are decreased by 37.6% and 52.6% respectively compared with the classical method.

biometric; in-air signature; identity verification; sequence alignment; accelerometer

2016- 08- 15;

2016- 12- 26。 基金項目:國家自然科學基金資助項目(61303028)。

劉威(1993—),男,安徽阜陽人,碩士研究生,主要研究方向:模式識別、圖像處理; 汪陽(1977—),男,湖北武漢人,副教授,博士,主要研究方向:機器人控制、嵌入式系統; 鄭建彬(1966—),男,湖北黃岡人,教授,博士,主要研究方向:模式識別、嵌入式系統; 詹恩奇(1972—),男,河南新野人,副教授,博士,主要研究方向:信號處理、模式識別。

1001- 9081(2017)04- 1004- 04

10.11772/j.issn.1001- 9081.2017.04.1004

TP391

A

猜你喜歡
用戶系統
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
基于PowerPC+FPGA顯示系統
半沸制皂系統(下)
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
關注用戶
商用汽車(2016年11期)2016-12-19 01:20:16
關注用戶
商用汽車(2016年6期)2016-06-29 09:18:54
關注用戶
商用汽車(2016年4期)2016-05-09 01:23:12
Camera360:拍出5億用戶
創業家(2015年10期)2015-02-27 07:55:08
主站蜘蛛池模板: 91无码人妻精品一区二区蜜桃| 97国产在线播放| AV网站中文| 国产亚洲精品自在线| 午夜不卡福利| 午夜丁香婷婷| 不卡午夜视频| 亚洲国产第一区二区香蕉| 青青国产在线| 色综合天天视频在线观看| 国产在线91在线电影| 九色国产在线| 精品国产成人三级在线观看| 亚洲不卡无码av中文字幕| 综合色在线| 欧美在线综合视频| 亚洲天堂区| 特级毛片8级毛片免费观看| 亚洲精品日产AⅤ| 国产精品七七在线播放| 麻豆精品在线| 3p叠罗汉国产精品久久| 免费人成网站在线高清| 国产区精品高清在线观看| 免费aa毛片| 伊人婷婷色香五月综合缴缴情| 老汉色老汉首页a亚洲| 欧美日韩国产综合视频在线观看| 玖玖免费视频在线观看| 久久综合国产乱子免费| 青青久视频| 午夜小视频在线| 国产天天色| 亚洲IV视频免费在线光看| 狠狠色综合久久狠狠色综合| 青草午夜精品视频在线观看| 欧美日韩在线观看一区二区三区| 91小视频在线观看| 亚洲综合极品香蕉久久网| 无遮挡国产高潮视频免费观看 | 韩日无码在线不卡| 天天做天天爱夜夜爽毛片毛片| 伊人91在线| 婷婷中文在线| 国产大片黄在线观看| 日本AⅤ精品一区二区三区日| 狠狠做深爱婷婷久久一区| 中文字幕无线码一区| 999在线免费视频| 人妻夜夜爽天天爽| 91人妻日韩人妻无码专区精品| 国产sm重味一区二区三区| 国产成人亚洲综合A∨在线播放| 色哟哟国产精品| 亚洲综合色吧| 亚洲系列无码专区偷窥无码| 四虎永久在线| 97se亚洲综合| 免费 国产 无码久久久| 99re视频在线| 精品久久777| 人妻中文久热无码丝袜| 国产麻豆福利av在线播放| 白浆免费视频国产精品视频| 在线国产三级| 成人亚洲天堂| 亚洲狼网站狼狼鲁亚洲下载| 成人精品在线观看| 第一页亚洲| 欧美性天天| 亚洲精品视频免费| 国产精品福利尤物youwu| 日韩毛片在线视频| 国产精品主播| 国产女人在线观看| 亚洲永久视频| 日本不卡在线| 国产日韩丝袜一二三区| 成人毛片在线播放| 婷婷色狠狠干| 国产福利免费在线观看| 亚洲中久无码永久在线观看软件|