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

    一種基于一階差分的野值類型判別及處理方法1

    2015-09-07 08:29:10饒云峰白燕
    時間頻率學(xué)報 2015年4期
    關(guān)鍵詞:野值散點原始數(shù)據(jù)

    饒云峰,白燕,3

    一種基于一階差分的野值類型判別及處理方法1

    饒云峰1,2,白燕1,2,3

    (1.中國科學(xué)院 國家授時中心,西安 710600;2.中國科學(xué)院大學(xué),北京 100049;3.中國科學(xué)院 精密導(dǎo)航定位與定時技術(shù)重點實驗室,西安 710600)

    為了剔除違反規(guī)律的異常測量值,采用一種基于一階差分的野值類型判別及處理方法。仿真結(jié)果表明,該算法能準(zhǔn)確地判別野值類型及其位置,尤其是對于斑點型野值具有較好的效果,同時可以提高數(shù)據(jù)處理精度,縮短數(shù)據(jù)處理時間。

    斑點型野值;散點型野值;一階差分;判別和修正

    0 引言

    測量設(shè)備記錄的測量數(shù)據(jù)往往包含較大的隨機誤差,個別值可能與其他值截然不同,違反了數(shù)據(jù)的整體變化規(guī)律,這些異常測量值就是野值[1]。工程測量中不可避免地存在野值,其對測量精度和數(shù)據(jù)處理的結(jié)果都會帶來嚴重的影響[2],在數(shù)據(jù)處理過程中必須對這這些違反測量規(guī)律的野值進行修正和剔除。野值的表現(xiàn)形式主要有以下2種[3]:

    ①散點型野值:野值以孤立點的形式出現(xiàn)。具體表現(xiàn)是:i時刻的觀測數(shù)據(jù)為野值,而在包含i時刻的某一時段內(nèi)其他觀測數(shù)據(jù)是正常的。

    ②斑點型野值:因相關(guān)性影響,野值成片出現(xiàn),若i時刻的觀測數(shù)據(jù)為野值,其后一段時間內(nèi)的值或絕大部分值也是異常的。

    野值點的成因、野值類型、野值模型等可參見文獻[1]~[3]。野值類型不同,剔除方法也不同。很多野值識別或剔除方法實現(xiàn)起來較復(fù)雜,如文獻[1]~[6]中提到的方法;有些方法只對散點型野值有很好的處理結(jié)果,而對斑點型野值的處理結(jié)果卻并不是很理想,如文獻[7]中提到的統(tǒng)計學(xué)方法和一些非統(tǒng)計學(xué)方法(Tukey53H法和灰自助法)、文獻[8]中提到的相位空間法。文獻[7]~[8]中提到的方法能很好地剔除散點型野值,而對于斑點型野值就很難剔除,尤其是文獻[8]中提到的相位空間法,它能很好地剔除散點型野值,且相比于其他野值剔除算法其最大優(yōu)勢是不需要設(shè)定門限值[7-8],但當(dāng)工程實踐中的測量數(shù)據(jù)出現(xiàn)大量斑點型野值時,即使多次利用相位空間法也無法剔除斑點型野值。為了解決此問題,本文提出一種基于一階差分的野值類型判別及處理方法。該方法可以判別數(shù)據(jù)中是否存在斑點型野值,精確找出斑點型野值在原始數(shù)據(jù)中的位置并加以剔除。

    1 野值類型判別及處理算法

    1.1野值類型判別算法

    當(dāng)測量數(shù)據(jù)中存在野值,則該數(shù)據(jù)序列的一階差分序列會出現(xiàn)明顯的函數(shù)跳躍點(通過與經(jīng)驗閾值比較后確定是否為跳躍點)。若出現(xiàn)單點跳躍點,表明此點前后數(shù)據(jù)中至少有一段數(shù)據(jù)是斑點型野值;若出現(xiàn)連續(xù)2個以上的跳躍點,表明這段連續(xù)跳躍點的第一個和最后一個跳躍點的前后數(shù)據(jù)中至少有一段數(shù)據(jù)是斑點型野值;若僅出現(xiàn)2個連續(xù)的跳躍點,表明原始數(shù)據(jù)中對應(yīng)2跳躍點后者的點為散點。這些跳躍點把原始數(shù)據(jù)分成若干數(shù)據(jù)段和散點型野值點,求每段數(shù)據(jù)的均值,并將得到的均值與對應(yīng)擬合數(shù)據(jù)段的均值進行比較,兩者之差超過經(jīng)驗閾值的這段數(shù)據(jù)即為斑點型野值。這樣就確定了野值類型及野值在原始數(shù)據(jù)中的位置。設(shè)原始數(shù)據(jù)序列為 u(i),1≤i≤n,n為原始測量數(shù)據(jù)的總長度。判斷原始數(shù)據(jù)中野值類型及位置的具體步驟如下:

    ①求 u(i)的一階差分:Δu(i)=u(i+1)-u(i),1≤i≤n;

    ③散點型野值類型判斷:若Mark中的某段相鄰元素有且只有2個元素值是連續(xù)的,則后一個元素對應(yīng)的原始數(shù)據(jù)為散點型野值;

    ④數(shù)據(jù)分段:根據(jù)Mark中相鄰元素值是否連續(xù),可將原始數(shù)據(jù)分段。具體分段方法為:當(dāng)Mark中某個元素的值與其后一個元素的值不連續(xù)時,取該元素加入數(shù)列End;當(dāng)Mark中某段元素僅有2個值連續(xù)時,取前一個元素加入End,后一個元素作為散點型野值點;當(dāng)Mark中某段元素有2個以上的值連續(xù)時,取該段連續(xù)值的第一個元素和最后一個元素加入End;另外將原始數(shù)據(jù)的最后一個數(shù)據(jù)點的序數(shù)也作為最后一個元素加入End;所構(gòu)成的數(shù)列End的長度為m,此時就將原始數(shù)據(jù)分為了m段。對于Mark中僅有2個值連續(xù)的情況,若放入End的元素對應(yīng)為End(j),則其對應(yīng)的第j+1段數(shù)據(jù)的起點為End(j)+2,終點為End(j+1);其余情況時,若第j個數(shù)據(jù)段的終止點為End(j),起始點為End(j -1)+1;另外對于第一個數(shù)據(jù)段,其終止點為End中的第一個元素,起始點為原始數(shù)據(jù)的第一個數(shù)據(jù);

    ⑥對判斷出的野值進行修正。

    通過上述步驟①~⑤就可識別出原始數(shù)據(jù)中的斑點型野值和散點型野值,并精確定位其在原始數(shù)據(jù)中的位置。

    1.2野值修正

    在進行野值修正時,若原始數(shù)據(jù)量大,則可直接剔除野值點,但這樣會導(dǎo)致數(shù)據(jù)的不連續(xù);也可將野值用一定值替代,其他正常值保持不變。一般有以下幾種替代方法[8]:

    ④擬合值;

    ⑤根據(jù)野值兩端的數(shù)值做插值。

    對于靜態(tài)數(shù)據(jù)(無運動趨勢的測量數(shù)據(jù)),上述5種野值替代的方法都可使用,但對于動態(tài)數(shù)據(jù)(數(shù)據(jù)有明顯的運動趨勢),方法①~③反而會使數(shù)據(jù)處理精度變得更差,所以對于動態(tài)數(shù)據(jù)可使用方法④~⑤進行野值修正,也可以把動態(tài)數(shù)據(jù)的運動趨勢去除[9]后,再用上述5種方法中的任意一種替代野值,野值剔除后再把該運動趨勢加回至剔除野值后的數(shù)據(jù)中。當(dāng)野值對整體數(shù)據(jù)影響比較大時,可多次使用本文提出的野值判別及修正方法,直至數(shù)據(jù)中無明顯的粗大值。

    2 數(shù)據(jù)處理及結(jié)果分析

    2.1仿真數(shù)據(jù)處理

    為了進一步說明不同類型的野值判別方法,模擬仿真了一組帶有白噪聲的數(shù)據(jù),其長度n=1000,均值為500,在此基礎(chǔ)上將點80,150,500,760改為散點型野值,點[20,35],[255,285]改為斑點型野值,野值總數(shù)為51,占總數(shù)據(jù)的5.1%,記該組仿真數(shù)據(jù)為data。仿真數(shù)據(jù)data如圖1所示。

    圖1原始數(shù)據(jù)

    對原始數(shù)據(jù)data求一階差分,一階差分的絕對值與閾值1R比較,得到一階差分函數(shù)中的跳躍點。本文設(shè)跳躍點的判決閾值113Rσ=,1σ為一階差分的標(biāo)準(zhǔn)方差。原始數(shù)據(jù)的一階差分如圖2所示。圖2中,實線代表原始數(shù)據(jù)的一階差分,虛線代表跳躍點的判決閾值。

    圖2 原始數(shù)據(jù)的一階差分及其閾值

    圖3為2種類型野值的一階差分局部放大圖,由于篇幅限制,此處只給出其中2段存在野值的數(shù)據(jù)(19,35),(499,500)的情況。一階差分與閾值1R比較,跳躍點的位置為Mark={19,35,79,80,149,150,254,285,499,500,759,760},其中[79,80],[149,150],[499,500],[759,760]是2個連續(xù)跳躍的點(圖3(b)所示),根據(jù)2.1節(jié)中的步驟③可知,原始數(shù)據(jù)data中的點80,150,500,760為散點型野值。一階差分在點19,35,254,285處為單個的跳躍點(圖3(a)所示),由步驟④可得,原始數(shù)據(jù)被分為[1,19],[20,35],[36,79],[81,149],[151,254],[255,285],[286,499],[501,759],[761,1000]9段數(shù)據(jù)和散點80,150,500,760,需進一步對數(shù)據(jù)組判斷,判斷其是否為斑點型的野值。

    圖3 一階差分對應(yīng)野值的局部放大圖

    根據(jù)步驟⑤判斷上述9段數(shù)據(jù)是否為斑點型野值,設(shè)斑點型野值判決閾值為各數(shù)據(jù)段的均值減去對應(yīng)擬合數(shù)據(jù)段的均值的標(biāo)準(zhǔn)方差,即為求標(biāo)準(zhǔn)方差函數(shù)。經(jīng)過比較,只有數(shù)據(jù)組[20,35],[255,285]滿足條件,所以數(shù)據(jù)組[20,35],[255,285]為斑點型野值。

    綜上可得,點80,150,500,760被識別為散點型野值,數(shù)據(jù)組[20,35],[255,285]被識別為斑點型野值,該判別結(jié)果與仿真數(shù)據(jù)中的野值類型及位置一致。

    由于該組仿真數(shù)據(jù)為靜態(tài)數(shù)據(jù),所以這里用原始數(shù)據(jù)的均值代替野值,其他數(shù)據(jù)保持不變。野值修正結(jié)果如圖4所示。

    圖4 原始數(shù)據(jù)與野值剔除結(jié)果比較圖

    圖4中,實線代表原始數(shù)據(jù),帶“*”的曲線為使用一次本文提出算法的數(shù)據(jù)處理結(jié)果。從圖中可以看出,[20,35],[255,285]中的斑點型野值被完全修正,點80,150,500,760處的散點型野值也被完全修正。所以本文提出的算法能準(zhǔn)確判別野值類型,并精確定位野值在原始數(shù)據(jù)中的位置。

    2.2實測數(shù)據(jù)處理結(jié)果

    為了驗證本文提出的算法的可行性和有效性,將本文提出的算法應(yīng)用于處理一組實際無線環(huán)境中的測距數(shù)據(jù),并將經(jīng)過該算法處理后的數(shù)據(jù)與原始數(shù)據(jù)進行比較。同時利用本文提出的算法與相位空間法分別對該組實測數(shù)據(jù)進行處理,并比較2種算法對于測距精度的影響。

    1)基于一階差分的野值類型判別及處理結(jié)果

    圖5為中國科學(xué)院國家授時中心2015年1月的某次無線綜合試驗中一組偽距測量原始數(shù)據(jù)。從圖5可知,該段測量數(shù)據(jù)在點(250,410)之間明顯包含數(shù)段斑點型野值,整段數(shù)據(jù)中也包含一些散點型野值。根據(jù)2.1節(jié)所述步驟對該組數(shù)據(jù)進行處理,取為一階差分的標(biāo)準(zhǔn)方差;σ2為所分數(shù)據(jù)段的均值減去對應(yīng)段擬合數(shù)據(jù)的均值的標(biāo)準(zhǔn)方差,即,野值用6階最小二乘擬合值代替。

    圖6給出利用本文提出的算法對野值進行修正后的結(jié)果,實線代表原始數(shù)據(jù),帶“*”的曲線表示本文提出的算法的數(shù)據(jù)處理結(jié)果。從圖中可以看出,原始數(shù)據(jù)中的點(250,410)內(nèi)的野值被完全修正,處理結(jié)果達到了預(yù)期的目的。

    圖5 實測原始數(shù)據(jù)

    圖6 基于一階差分的野值類型判別及處理方法

    2)相位空間法野值修正結(jié)果

    相位空間法[8]:以原始數(shù)據(jù) u(i)的殘差及原始數(shù)據(jù)的二階導(dǎo)為軸畫一個橢圓,其中,的一階導(dǎo)。若數(shù)據(jù)點在橢圓外,則判為野值,用修正值代替;若數(shù)據(jù)點在橢圓內(nèi),則判為正常的測量數(shù)據(jù),不處理。

    為了與基于一階差分的野值類型判別及處理方法比較,此方法中的野值也用6階最小二乘擬合值代替。

    圖7為循環(huán)使用相位空間法4次修正后的數(shù)據(jù),實線代表原始數(shù)據(jù),帶“o”的曲線代表循環(huán)使用相位空間法4次的野值剔除結(jié)果。對該組數(shù)據(jù),循環(huán)使用相位空間法4次后,所有的數(shù)據(jù)都滿足在橢圓內(nèi)的條件,即該組數(shù)據(jù)最多只能循環(huán)使用相位空間法4次。從圖中可以看出,原始數(shù)據(jù)中的點(250,410)內(nèi)的斑點型野值只有部分被修正,處理結(jié)果并不是很理想,但數(shù)據(jù)中的部分散點型野值得到相應(yīng)的修正。

    圖7 循環(huán)使用相位空間法4次修正后的數(shù)據(jù)

    3)2種方法的比較

    圖8給出了采用2種方法處理后的數(shù)據(jù)與原始數(shù)據(jù)的比對結(jié)果,實線代表原始數(shù)據(jù),帶“o”的曲線代表循環(huán)使用相位空間法4次的野值剔除結(jié)果,帶“*”的曲線代表本文提出的算法的數(shù)據(jù)處理結(jié)果。從圖中可以看出,2種方法都能修正數(shù)據(jù)中的散點型野值,但是對(250,410)區(qū)間內(nèi)的斑點型野值,相位空間法處理結(jié)果并不是很好,而采用本文提出的方法能判斷出該部分的斑點型野值并能很好地修正該部分斑點型野值。

    圖8 2種方法處理結(jié)果與原始數(shù)據(jù)的比較

    另外,不同的數(shù)據(jù)處理方法對測量數(shù)據(jù)的精度會產(chǎn)生不同的影響。對采用上述2種算法進行處理后的數(shù)據(jù)分別求取其均值和標(biāo)準(zhǔn)方差,結(jié)果如表1所示。從表中可看出,采用2種算法分別對該組數(shù)據(jù)進行預(yù)處理后再求其測距精度,相比對沒有經(jīng)過處理的原始數(shù)據(jù)直接求得的測距精度,均有明顯提高;另外,當(dāng)處理次數(shù)相同(1次)的情況下,經(jīng)過本文提出的算法處理后求得的測距精度要比經(jīng)過相位空間法處理后的測距精度高出0.0115;即使循環(huán)使用4次相位空間法,其野值修正后的測距精度也沒有使用本文算法(1次)后得到的測距精度高。因此,本文提出的算法在處理速度上要快于相位空間法,且數(shù)據(jù)處理后的精度也明顯優(yōu)于相位空間法。

    表 1 實測數(shù)據(jù)的處理結(jié)果

    3 結(jié)論

    野值的判別及剔除是測量數(shù)據(jù)處理過程中非常重要的一個環(huán)節(jié),本文提出了一種基于一階差分的野值類型判別及處理方法,該方法能準(zhǔn)確判別數(shù)據(jù)中的野值類型及其在原數(shù)據(jù)中的位置,相比相位空間法,該算法具有更高的數(shù)據(jù)處理精度和數(shù)據(jù)處理速度,可為有關(guān)工程實踐中的數(shù)據(jù)處理提供一定的參考。

    [1]黃家貴,呂紅寧,王安麗.斑點型野值的進一步識別與處理[J].裝備指揮技術(shù)學(xué)院學(xué)報,2002,13(6):56-59.

    [2]胡紹林,孫國基.靶場外測數(shù)據(jù)野值點的統(tǒng)計診斷技術(shù)[J].宇航學(xué)報,1999,20(2):68-73.

    [3]祝轉(zhuǎn)明,秋宏興,李濟生,等.動態(tài)測量數(shù)據(jù)野值的辨識與剔除[J].系統(tǒng)工程與電子技術(shù),2004,26(2):147-149,190.

    [4]胡峰,孫國基.Kalman濾波的抗野值修正[J].自動化學(xué)報,1999,25(5):694-695.

    [5]王佳,王敏,徐曉輝.一種干涉儀數(shù)據(jù)斑點型野值修正方法[J].彈箭與制導(dǎo)學(xué)報,2013,33(2):106-108.

    [6]卓寧.靶場外彈道數(shù)據(jù)處理中野值點剔除方法[J].測試技術(shù)學(xué)報,2008,22(4):313-317.

    [7]楊筱.衛(wèi)星導(dǎo)航系統(tǒng)數(shù)據(jù)與信號質(zhì)量評估技術(shù)研究[D].長沙:國防科學(xué)技術(shù)大學(xué),2009.

    [8]GORING D G.NIKORA V I.Despiking acoustic doppler velocimeter data[J].Journal of Hydraulic Engineering,2002,128(1): 117-126.

    [9]BAHAER.趨勢項提取+萊特準(zhǔn)則的野值剔除方法[EB/OL].(2010-10-10)[2015-01-13].http://blog.sina.com.cn/s/blog_5de f5a660100m3mr.html.

    Amethod based on first-order difference for type-judging and processing of outliers

    RAO Yun-feng1,2,BAI Yan1,2,3
    (1.National Time Service Center,ChineseAcademy of Sciences,Xi′an 710600,China;2.University of ChineseAcademy of Sciences,Beijing 100049,China;3.Key Laboratory of Precision Navigation and Timing Technology,National Time Service Center,ChineseAcademy of Sciences,Xi′an 710600,China)

    For eliminating the abnormal measured values,a method based on the first-order difference is adopted for judging the types of the outliers and processing them.The simulation results show that the type of outliers can be accurately differentiated and the location of outliers in the original data can be precisely pinpointed by using this algorithm,especially for the spotted outliers the processing result is rather good.At the same time,processing precision can be improved and the processing time can be reduced.

    spotted outliers;isolated outliers;first-order difference;differentiation and correction

    TP274

    A

    1674-0637(2015)04-0227-08

    10.13875/j.issn.1674-0637.2015-04-0227-08

    2015-06-25

    國家自然科學(xué)基金資助項目(11203027)

    饒云峰,女,碩士研究生,主要從事Ka頻段星間鏈路信號質(zhì)量性能評估方法的研究。

    猜你喜歡
    野值散點原始數(shù)據(jù)
    項目化學(xué)習(xí):讓知識學(xué)習(xí)從“散點”走向“聚合”
    散點透視“宋詞運用唐詩”
    中華詩詞(2023年2期)2023-07-31 02:18:06
    GOLDEN OPPORTUNITY FOR CHINA-INDONESIA COOPERATION
    基于回歸模型的采集數(shù)據(jù)清洗技術(shù)
    電光與控制(2022年4期)2022-04-07 12:47:06
    受特定變化趨勢限制的傳感器數(shù)據(jù)處理方法研究
    彈道跟蹤數(shù)據(jù)野值剔除方法性能分析
    上海航天(2018年4期)2018-09-07 03:28:46
    沖擊波散點和定點治療肱二頭肌長頭肌腱炎的臨床觀察
    導(dǎo)引頭在線野值剔除與平滑方法研究
    全新Mentor DRS360 平臺借助集中式原始數(shù)據(jù)融合及直接實時傳感技術(shù)實現(xiàn)5 級自動駕駛
    汽車零部件(2017年4期)2017-07-12 17:05:53
    RS422總線數(shù)據(jù)野值剔除方法研究
    国产男人的电影天堂91| 久久精品国产亚洲av涩爱| 亚洲av欧美aⅴ国产| 日本vs欧美在线观看视频| 黄片播放在线免费| 亚洲成人一二三区av| 成年动漫av网址| 美女脱内裤让男人舔精品视频| 少妇被粗大猛烈的视频| 久久精品久久精品一区二区三区| 国产69精品久久久久777片| 免费人成在线观看视频色| 国产熟女欧美一区二区| 成人国产av品久久久| av免费在线看不卡| 又黄又粗又硬又大视频| 在线观看人妻少妇| 成年美女黄网站色视频大全免费| 亚洲av福利一区| 日韩欧美精品免费久久| 在线观看人妻少妇| 边亲边吃奶的免费视频| 啦啦啦在线观看免费高清www| av国产久精品久网站免费入址| 中文字幕亚洲精品专区| 伦理电影大哥的女人| 欧美人与性动交α欧美精品济南到 | 18禁国产床啪视频网站| 中文欧美无线码| 91久久精品国产一区二区三区| 2021少妇久久久久久久久久久| 久久精品国产a三级三级三级| 免费看光身美女| 91国产中文字幕| 国产视频首页在线观看| 麻豆乱淫一区二区| 91精品三级在线观看| 2022亚洲国产成人精品| 一级毛片电影观看| 丰满饥渴人妻一区二区三| 国产午夜精品一二区理论片| 一级片免费观看大全| 啦啦啦在线观看免费高清www| av不卡在线播放| 精品人妻一区二区三区麻豆| 免费女性裸体啪啪无遮挡网站| 一级毛片 在线播放| 精品第一国产精品| 免费观看无遮挡的男女| 亚洲av免费高清在线观看| 国产 精品1| 国产精品嫩草影院av在线观看| 水蜜桃什么品种好| 亚洲精品一区蜜桃| 99精国产麻豆久久婷婷| 大片电影免费在线观看免费| 人妻少妇偷人精品九色| 欧美激情 高清一区二区三区| 这个男人来自地球电影免费观看 | 免费不卡的大黄色大毛片视频在线观看| 欧美日韩国产mv在线观看视频| 香蕉精品网在线| 国产探花极品一区二区| 亚洲国产av新网站| 日韩欧美精品免费久久| 久久精品熟女亚洲av麻豆精品| 国产女主播在线喷水免费视频网站| 精品午夜福利在线看| 在现免费观看毛片| 蜜桃在线观看..| 汤姆久久久久久久影院中文字幕| 精品少妇久久久久久888优播| 如何舔出高潮| 狠狠婷婷综合久久久久久88av| 中国三级夫妇交换| a 毛片基地| 午夜免费观看性视频| 男的添女的下面高潮视频| 欧美国产精品一级二级三级| 亚洲av综合色区一区| 在线精品无人区一区二区三| 精品一区二区三区四区五区乱码 | 大香蕉久久网| 大码成人一级视频| 国产日韩欧美视频二区| 国产精品.久久久| 久久av网站| av在线观看视频网站免费| 成人18禁高潮啪啪吃奶动态图| 麻豆乱淫一区二区| 亚洲国产日韩一区二区| 我的女老师完整版在线观看| 亚洲欧美一区二区三区国产| 在线天堂中文资源库| 免费播放大片免费观看视频在线观看| 欧美 亚洲 国产 日韩一| 日本欧美视频一区| 街头女战士在线观看网站| av国产精品久久久久影院| 精品99又大又爽又粗少妇毛片| 久久女婷五月综合色啪小说| 国产精品三级大全| 黄色 视频免费看| 国产成人精品福利久久| 大片电影免费在线观看免费| 蜜桃国产av成人99| 人体艺术视频欧美日本| 少妇的逼好多水| 最近2019中文字幕mv第一页| 国产一区二区三区综合在线观看 | 成年女人在线观看亚洲视频| 国产精品一国产av| 一本大道久久a久久精品| 美女中出高潮动态图| 精品人妻偷拍中文字幕| 人人澡人人妻人| 纯流量卡能插随身wifi吗| 亚洲人成77777在线视频| 成人国产av品久久久| 91久久精品国产一区二区三区| 激情视频va一区二区三区| 九草在线视频观看| 国产成人精品婷婷| 久久人人爽人人片av| 精品少妇黑人巨大在线播放| 99热网站在线观看| 少妇 在线观看| 中国美白少妇内射xxxbb| 香蕉国产在线看| 卡戴珊不雅视频在线播放| av黄色大香蕉| 日产精品乱码卡一卡2卡三| 在线观看国产h片| 18禁在线无遮挡免费观看视频| 午夜激情av网站| 在线观看一区二区三区激情| 婷婷色麻豆天堂久久| 午夜免费鲁丝| 国产日韩欧美亚洲二区| 亚洲av国产av综合av卡| 嫩草影院入口| 久久毛片免费看一区二区三区| 两性夫妻黄色片 | 久久亚洲国产成人精品v| 99香蕉大伊视频| 欧美日韩国产mv在线观看视频| 精品一区二区三卡| 一区在线观看完整版| av电影中文网址| 日韩电影二区| 2022亚洲国产成人精品| 久久国产亚洲av麻豆专区| 涩涩av久久男人的天堂| 日产精品乱码卡一卡2卡三| 黄色毛片三级朝国网站| 91成人精品电影| 久久免费观看电影| 日本爱情动作片www.在线观看| av线在线观看网站| 国产精品 国内视频| 在线观看美女被高潮喷水网站| 久久人人97超碰香蕉20202| 欧美激情极品国产一区二区三区 | 99热全是精品| 一级毛片黄色毛片免费观看视频| 欧美性感艳星| 国产深夜福利视频在线观看| 日韩人妻精品一区2区三区| 综合色丁香网| 免费黄频网站在线观看国产| 制服诱惑二区| 亚洲四区av| 亚洲,欧美精品.| 成人国产麻豆网| 国产精品嫩草影院av在线观看| 久久av网站| 色婷婷av一区二区三区视频| 97精品久久久久久久久久精品| 国产一区二区三区综合在线观看 | 欧美激情 高清一区二区三区| 亚洲欧洲国产日韩| 亚洲五月色婷婷综合| 国产白丝娇喘喷水9色精品| 国产欧美另类精品又又久久亚洲欧美| 久久久精品区二区三区| 满18在线观看网站| 午夜福利视频精品| 亚洲精华国产精华液的使用体验| 最近最新中文字幕免费大全7| 中文乱码字字幕精品一区二区三区| 亚洲久久久国产精品| 制服诱惑二区| 丰满迷人的少妇在线观看| 少妇熟女欧美另类| 亚洲精品一二三| 夫妻性生交免费视频一级片| 少妇人妻 视频| 日日摸夜夜添夜夜爱| 三级国产精品片| 久久精品熟女亚洲av麻豆精品| 欧美精品人与动牲交sv欧美| 午夜福利视频精品| 国产毛片在线视频| 有码 亚洲区| 国产精品久久久久久久久免| 精品人妻偷拍中文字幕| 春色校园在线视频观看| 在线天堂中文资源库| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久精品古装| 日本免费在线观看一区| 国产av国产精品国产| 男女无遮挡免费网站观看| 中文字幕人妻熟女乱码| 久久精品夜色国产| 大话2 男鬼变身卡| 精品国产国语对白av| 97超碰精品成人国产| 亚洲av成人精品一二三区| 最近的中文字幕免费完整| 波野结衣二区三区在线| 青青草视频在线视频观看| 又黄又爽又刺激的免费视频.| 欧美日韩成人在线一区二区| 在线 av 中文字幕| av一本久久久久| av线在线观看网站| 欧美人与性动交α欧美精品济南到 | 亚洲精品自拍成人| 亚洲av免费高清在线观看| 精品少妇黑人巨大在线播放| 日本猛色少妇xxxxx猛交久久| 天天躁夜夜躁狠狠躁躁| 亚洲精品自拍成人| 美女国产高潮福利片在线看| 伦理电影大哥的女人| 男女国产视频网站| 巨乳人妻的诱惑在线观看| 69精品国产乱码久久久| 97人妻天天添夜夜摸| 久久国产精品大桥未久av| 久久影院123| 国产日韩欧美在线精品| 国产一区二区在线观看av| 国产精品国产三级国产专区5o| 全区人妻精品视频| 99久久综合免费| 宅男免费午夜| √禁漫天堂资源中文www| 久久精品国产a三级三级三级| freevideosex欧美| 美女国产高潮福利片在线看| 免费观看av网站的网址| 午夜91福利影院| 人人妻人人爽人人添夜夜欢视频| 一级毛片电影观看| 久久狼人影院| 亚洲欧洲国产日韩| 中文字幕亚洲精品专区| 亚洲色图 男人天堂 中文字幕 | 日本免费在线观看一区| 最近中文字幕高清免费大全6| 日韩电影二区| 大码成人一级视频| 国产亚洲av片在线观看秒播厂| 亚洲美女黄色视频免费看| 久久久久精品性色| 99视频精品全部免费 在线| 青青草视频在线视频观看| 在线观看免费视频网站a站| 精品国产一区二区久久| 纯流量卡能插随身wifi吗| 最新的欧美精品一区二区| 高清欧美精品videossex| 国产在线免费精品| 综合色丁香网| 人妻一区二区av| 大片电影免费在线观看免费| 亚洲少妇的诱惑av| 日韩欧美一区视频在线观看| 全区人妻精品视频| 啦啦啦中文免费视频观看日本| 午夜免费观看性视频| 日韩大片免费观看网站| 乱码一卡2卡4卡精品| 国产一区二区在线观看日韩| 免费看av在线观看网站| 亚洲精品久久久久久婷婷小说| 国产精品久久久久成人av| 久久久国产精品麻豆| 精品第一国产精品| 亚洲欧美成人精品一区二区| 国产视频首页在线观看| 精品视频人人做人人爽| 国产国拍精品亚洲av在线观看| 丝瓜视频免费看黄片| 国产精品无大码| 水蜜桃什么品种好| 国产永久视频网站| 男人爽女人下面视频在线观看| 人体艺术视频欧美日本| 午夜91福利影院| 黄色一级大片看看| 久久精品久久精品一区二区三区| 久久精品国产亚洲av涩爱| 国产高清国产精品国产三级| 国产又色又爽无遮挡免| 亚洲欧洲精品一区二区精品久久久 | 大码成人一级视频| 一级爰片在线观看| 各种免费的搞黄视频| 亚洲伊人色综图| 一区在线观看完整版| 99热这里只有是精品在线观看| 男的添女的下面高潮视频| 国产精品久久久久成人av| 免费高清在线观看日韩| 丰满少妇做爰视频| 99视频精品全部免费 在线| 在线观看免费高清a一片| 亚洲久久久国产精品| 国产精品一区二区在线不卡| 午夜免费男女啪啪视频观看| 国产亚洲最大av| 免费人成在线观看视频色| 久久久久视频综合| 久久久久久久亚洲中文字幕| 精品久久蜜臀av无| 亚洲欧美色中文字幕在线| av片东京热男人的天堂| 十八禁高潮呻吟视频| 国产国语露脸激情在线看| 飞空精品影院首页| 哪个播放器可以免费观看大片| 亚洲精品成人av观看孕妇| 一级片免费观看大全| 欧美人与性动交α欧美软件 | 久久久久网色| 一边摸一边做爽爽视频免费| 欧美+日韩+精品| 久久这里只有精品19| 黄色毛片三级朝国网站| 免费看不卡的av| 国产成人免费观看mmmm| 国产黄色视频一区二区在线观看| 波野结衣二区三区在线| 天堂8中文在线网| 赤兔流量卡办理| 精品久久蜜臀av无| 国产深夜福利视频在线观看| 亚洲精品视频女| 精品少妇内射三级| 18+在线观看网站| 自拍欧美九色日韩亚洲蝌蚪91| 少妇熟女欧美另类| 免费女性裸体啪啪无遮挡网站| 亚洲欧美日韩卡通动漫| 又黄又爽又刺激的免费视频.| 2022亚洲国产成人精品| 亚洲四区av| 亚洲成人av在线免费| 免费观看性生交大片5| 亚洲精华国产精华液的使用体验| 99热国产这里只有精品6| 亚洲欧美清纯卡通| 亚洲国产日韩一区二区| 中文字幕最新亚洲高清| 新久久久久国产一级毛片| 18在线观看网站| 欧美 亚洲 国产 日韩一| 亚洲人成77777在线视频| 观看美女的网站| 成年女人在线观看亚洲视频| www.色视频.com| 蜜桃国产av成人99| 国产免费福利视频在线观看| 毛片一级片免费看久久久久| 中国美白少妇内射xxxbb| 久久人人爽人人爽人人片va| 妹子高潮喷水视频| 九草在线视频观看| 日本wwww免费看| 男人爽女人下面视频在线观看| 国产日韩欧美在线精品| 国产福利在线免费观看视频| 色哟哟·www| 人妻系列 视频| 午夜福利,免费看| 免费观看a级毛片全部| av国产精品久久久久影院| 狠狠精品人妻久久久久久综合| 亚洲精品美女久久久久99蜜臀 | 汤姆久久久久久久影院中文字幕| 精品第一国产精品| 晚上一个人看的免费电影| 大香蕉久久网| 日韩不卡一区二区三区视频在线| 欧美日韩视频高清一区二区三区二| 国产一级毛片在线| 精品酒店卫生间| 最黄视频免费看| 亚洲高清免费不卡视频| 涩涩av久久男人的天堂| 精品熟女少妇av免费看| 国产激情久久老熟女| 在线免费观看不下载黄p国产| 考比视频在线观看| 日本av免费视频播放| 卡戴珊不雅视频在线播放| 国产亚洲精品第一综合不卡 | 18禁在线无遮挡免费观看视频| videos熟女内射| a级毛片在线看网站| 91在线精品国自产拍蜜月| 久久精品熟女亚洲av麻豆精品| 成人手机av| 国产免费一区二区三区四区乱码| 国产日韩欧美在线精品| 欧美日韩成人在线一区二区| 免费观看无遮挡的男女| 一个人免费看片子| 22中文网久久字幕| 女的被弄到高潮叫床怎么办| 最近最新中文字幕大全免费视频 | av电影中文网址| 日韩,欧美,国产一区二区三区| 中文字幕最新亚洲高清| 亚洲丝袜综合中文字幕| 我的女老师完整版在线观看| 久久久久久人妻| 一区二区三区精品91| 你懂的网址亚洲精品在线观看| 亚洲成人手机| 精品人妻熟女毛片av久久网站| 少妇的丰满在线观看| 久久精品国产a三级三级三级| 国产精品一国产av| 秋霞在线观看毛片| 永久免费av网站大全| 亚洲精品国产av成人精品| 亚洲久久久国产精品| 搡老乐熟女国产| 两个人免费观看高清视频| 18禁裸乳无遮挡动漫免费视频| 久久人人97超碰香蕉20202| 亚洲内射少妇av| 99久久人妻综合| 99热网站在线观看| 老女人水多毛片| 99热网站在线观看| 91成人精品电影| 少妇的丰满在线观看| 国产又爽黄色视频| 国产亚洲欧美精品永久| 午夜老司机福利剧场| 欧美日韩一区二区视频在线观看视频在线| 精品国产国语对白av| 欧美日韩综合久久久久久| 99九九在线精品视频| 熟女电影av网| 中文字幕精品免费在线观看视频 | 国语对白做爰xxxⅹ性视频网站| 欧美日韩综合久久久久久| 午夜福利影视在线免费观看| 国产精品久久久久久久电影| 国产亚洲精品久久久com| 国产免费一级a男人的天堂| 国产精品一区www在线观看| 国产永久视频网站| 亚洲精品第二区| 久久99蜜桃精品久久| 久热久热在线精品观看| 欧美人与性动交α欧美精品济南到 | 精品国产一区二区三区久久久樱花| 99久久综合免费| 99视频精品全部免费 在线| 国产精品久久久久成人av| 国产女主播在线喷水免费视频网站| 黄色视频在线播放观看不卡| 国产黄色免费在线视频| 在线观看www视频免费| 国产精品一国产av| 久久精品久久久久久久性| 欧美精品一区二区免费开放| 亚洲国产欧美日韩在线播放| 亚洲天堂av无毛| 午夜精品国产一区二区电影| 人人妻人人澡人人看| 国产在线一区二区三区精| 一级毛片 在线播放| 天天影视国产精品| 久久久国产一区二区| 国产av一区二区精品久久| 国产精品欧美亚洲77777| 热99久久久久精品小说推荐| 人人澡人人妻人| 国产在线免费精品| 日韩欧美精品免费久久| 亚洲欧美一区二区三区黑人 | 男人舔女人的私密视频| 99热网站在线观看| 久久国产精品男人的天堂亚洲 | 久久久精品免费免费高清| 搡老乐熟女国产| 亚洲精品美女久久av网站| 精品卡一卡二卡四卡免费| 国产熟女欧美一区二区| 成年人午夜在线观看视频| 男女下面插进去视频免费观看 | 制服诱惑二区| 熟妇人妻不卡中文字幕| 99久久中文字幕三级久久日本| 夜夜爽夜夜爽视频| 男女高潮啪啪啪动态图| 熟女av电影| 啦啦啦视频在线资源免费观看| 亚洲成人av在线免费| 男人爽女人下面视频在线观看| 欧美精品人与动牲交sv欧美| 99热这里只有是精品在线观看| 日韩中字成人| 亚洲精品一区蜜桃| 欧美精品国产亚洲| 18+在线观看网站| 啦啦啦中文免费视频观看日本| 国产精品一区二区在线观看99| 热99国产精品久久久久久7| 美女大奶头黄色视频| 久久久精品区二区三区| 日韩精品有码人妻一区| 一个人免费看片子| 欧美少妇被猛烈插入视频| 亚洲精品美女久久久久99蜜臀 | 亚洲第一av免费看| 97在线视频观看| 国产伦理片在线播放av一区| 中文字幕av电影在线播放| 你懂的网址亚洲精品在线观看| 97人妻天天添夜夜摸| 成人黄色视频免费在线看| 久久精品国产a三级三级三级| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精华国产精华液的使用体验| av国产精品久久久久影院| 插逼视频在线观看| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| 99久久人妻综合| 久久久国产一区二区| 欧美变态另类bdsm刘玥| 校园人妻丝袜中文字幕| 亚洲成人av在线免费| 亚洲四区av| 免费高清在线观看视频在线观看| 久久久久久久大尺度免费视频| 一本色道久久久久久精品综合| 少妇的逼水好多| 久久青草综合色| 最近中文字幕2019免费版| 亚洲欧美中文字幕日韩二区| 日本av免费视频播放| av黄色大香蕉| 国国产精品蜜臀av免费| 亚洲精品美女久久久久99蜜臀 | 免费在线观看完整版高清| 国产成人免费无遮挡视频| 欧美日韩一区二区视频在线观看视频在线| 日韩 亚洲 欧美在线| 一级,二级,三级黄色视频| 欧美国产精品va在线观看不卡| 亚洲av国产av综合av卡| 在线亚洲精品国产二区图片欧美| 视频区图区小说| 免费少妇av软件| 黄片播放在线免费| 色网站视频免费| 黄色视频在线播放观看不卡| 欧美人与性动交α欧美精品济南到 | 久久久久久久亚洲中文字幕| 久久久精品94久久精品| 亚洲中文av在线| 亚洲精品日本国产第一区| 亚洲美女黄色视频免费看| 亚洲少妇的诱惑av| 国产精品欧美亚洲77777| 国产探花极品一区二区| 国产在视频线精品| 免费女性裸体啪啪无遮挡网站| 国产 精品1| 99香蕉大伊视频| 天天影视国产精品| 中国国产av一级| 欧美最新免费一区二区三区| 国产亚洲精品久久久com| 91午夜精品亚洲一区二区三区| 国产69精品久久久久777片| 一区二区av电影网| h视频一区二区三区| 久久精品久久精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 国产日韩欧美视频二区| 免费播放大片免费观看视频在线观看| 国产精品久久久av美女十八| 中国国产av一级| 久久国产精品大桥未久av| 国产精品99久久99久久久不卡 | 国产一区二区在线观看av| 在线观看一区二区三区激情| 中国国产av一级| 色网站视频免费| 97精品久久久久久久久久精品| 日韩中字成人| 巨乳人妻的诱惑在线观看| 亚洲情色 制服丝袜| 黑丝袜美女国产一区|