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

    基于SSA-Sn法的GNSS坐標(biāo)時(shí)間序列粗差探測(cè)

    2022-03-09 12:24:10徐華卿魯鐵定賀小星周世健
    關(guān)鍵詞:估計(jì)量原始數(shù)據(jù)測(cè)站

    徐華卿 魯鐵定 賀小星 周世健 陶 蕊

    1 東華理工大學(xué)測(cè)繪工程學(xué)院,南昌市廣蘭大道418號(hào),330013 2 江西理工大學(xué)土木與測(cè)繪工程學(xué)院,江西省贛州市紅旗大道86號(hào),341000 3 南昌航空大學(xué)校長(zhǎng)辦公室,南昌市豐和南大道696號(hào),330063

    以GNSS為主的衛(wèi)星導(dǎo)航定位基準(zhǔn)站網(wǎng)近年來(lái)發(fā)展迅速,為大地測(cè)量的應(yīng)用提供了高精度的空間基準(zhǔn)基礎(chǔ)設(shè)施[1]。GNSS坐標(biāo)時(shí)間序列蘊(yùn)含豐富的物理信息,包括測(cè)站長(zhǎng)期受區(qū)域內(nèi)構(gòu)造應(yīng)力影響產(chǎn)生的線性變化以及地球物理效應(yīng)與其他外部環(huán)境等導(dǎo)致的非線性變化[2]。由于測(cè)站會(huì)受到電離層延遲、多路徑效應(yīng)、環(huán)境負(fù)載、軌道異常等多種誤差影響[3],GNSS坐標(biāo)時(shí)間序列在呈現(xiàn)非線性變化的同時(shí)也存在異常值[4]。這些異常值在整體或局部數(shù)據(jù)中表現(xiàn)出較大的波動(dòng),無(wú)法準(zhǔn)確反映出時(shí)間序列的變化特征,對(duì)數(shù)據(jù)分析與應(yīng)用產(chǎn)生較大的負(fù)面影響。

    經(jīng)典的粗差探測(cè)法有3σ法[5]和MAD(median absolute deviation)法[6-7]。3σ法的原理是利用最小二乘法求得觀測(cè)序列的殘差,再對(duì)殘差序列進(jìn)行粗差探測(cè)。由于最小二乘法本身含有粗差,因此3σ法的殘差中誤差估計(jì)值不夠準(zhǔn)確,從而影響粗差探測(cè)的可靠性[8-11]。MAD方法的運(yùn)算崩潰污染率可達(dá)50%,其異常值探測(cè)相比于3σ法更加穩(wěn)健,但MAD法需要滿足樣本對(duì)稱分布的條件[12-13]?;谏鲜鲇懻?,諸多學(xué)者在此基礎(chǔ)上引入新的準(zhǔn)則和方法以提高粗差探測(cè)的效率和準(zhǔn)確率。明鋒等[8]提出L1范數(shù)與IQR組合探測(cè)方法,使模型化后的殘差探測(cè)準(zhǔn)確率更高;吳浩等[5]提出一種利用小波改進(jìn)的3σ粗差探測(cè)方法,既能準(zhǔn)確提取變形趨勢(shì),也能提高3σ在坐標(biāo)數(shù)據(jù)中的適用性;Wang等[14]發(fā)現(xiàn)有效提取趨勢(shì)項(xiàng)和周期項(xiàng)至關(guān)重要, SSA在這方面具有一定優(yōu)勢(shì),能準(zhǔn)確提取出殘差項(xiàng);蔡曉軍等[15]研究表明,在提取趨勢(shì)項(xiàng)和周期項(xiàng)方面SSA優(yōu)于小波分析法。

    綜合以上異常值探測(cè)方法的優(yōu)勢(shì)與不足并結(jié)合GNSS坐標(biāo)時(shí)間序列的數(shù)據(jù)特征,建立一種SSA與Sn估計(jì)量[14]相結(jié)合的異常值探測(cè)方法。選用模擬的時(shí)間序列數(shù)據(jù)以及中國(guó)部分IGS站的實(shí)測(cè)數(shù)據(jù)進(jìn)行實(shí)驗(yàn),結(jié)果表明,在粗差探測(cè)方面SSA-Sn法的探測(cè)效果優(yōu)于3σ法和MAD法。

    1 原理及方法

    1.1 SSA數(shù)學(xué)模型

    (1)

    C=VΛVT

    (2)

    將C的特征值進(jìn)行降序排列,得到λk,Λ為λk構(gòu)成的對(duì)角矩陣,正交矩陣V的行向量為vk。主成分矩陣為:

    (3)

    A的第k行ak為第k個(gè)主成分,其中ak,i為:

    (4)

    式中,vk,j是vk的第j個(gè)元素。時(shí)間序列中第k個(gè)重構(gòu)分量為:

    (5)

    (6)

    1.2 交叉驗(yàn)證法確定參數(shù)

    采用交叉驗(yàn)證法確定SSA需要的窗口長(zhǎng)度L與主成分個(gè)數(shù)P兩個(gè)參數(shù)。首先給出兩個(gè)參數(shù)的初始值,并在原始數(shù)據(jù)中隨機(jī)挑選部分?jǐn)?shù)據(jù)用于驗(yàn)證,驗(yàn)證數(shù)據(jù)處用0填充,即可得到訓(xùn)練數(shù)據(jù)。隨后將訓(xùn)練數(shù)據(jù)進(jìn)行奇異譜分析,計(jì)算驗(yàn)證部分模型導(dǎo)出值與參考值的均方根誤差RMSE,RMSE最小時(shí)的窗口長(zhǎng)度L與主成分個(gè)數(shù)P即為最優(yōu)參數(shù)。Wang 等[14]表明,2~3 a的窗口長(zhǎng)度適合提取年周期及半年周期分量,因此本文采用2 a作為窗口長(zhǎng)度的初值。

    1.3 Sn尺度估計(jì)量

    Sn是一個(gè)由Rousseeuw和Croux提出的穩(wěn)健統(tǒng)計(jì)量,Sn估計(jì)量具有與MAD法相同的穩(wěn)健性,并且不依賴于對(duì)稱分布。假設(shè)一組數(shù)據(jù)為{x1,x2,…,xn},則Sn估計(jì)量被定義為:

    Sn=cmedi{medi|xi-xj|}

    (7)

    式中,常數(shù)c是為了滿足Fisher一致性的調(diào)節(jié)因子,其值為1.192 6;外中位數(shù)是一個(gè)秩為[(n+1)/2]階的順序統(tǒng)計(jì)量,而內(nèi)中位數(shù)是一個(gè)秩為[n/2]+1的順序統(tǒng)計(jì)量。

    1.4 SSA-Sn粗差探測(cè)原理

    在對(duì)GNSS坐標(biāo)時(shí)間序列進(jìn)行粗差探測(cè)時(shí),可將時(shí)間序列數(shù)據(jù)視為一組被噪聲和粗差污染的趨勢(shì)信號(hào)。SSA相較于小波分析能夠更加準(zhǔn)確地提取趨勢(shì)項(xiàng)、周期項(xiàng)等主要的時(shí)間序列特征,從而分離出殘差項(xiàng);Sn與MAD法對(duì)異常值探測(cè)的穩(wěn)健性相同且都優(yōu)于3σ法,估計(jì)量崩潰點(diǎn)均為50%,不易受到粗差的污染,但Sn無(wú)需依賴數(shù)據(jù)的對(duì)稱分布,因此本文采用SSA-Sn法進(jìn)行粗差探測(cè)。具體步驟如下:

    1)輸入坐標(biāo)時(shí)間序列原始數(shù)據(jù)集X(t),使用交叉驗(yàn)證法確定最優(yōu)的窗口長(zhǎng)度L與主成分個(gè)數(shù)P;

    3)原始數(shù)據(jù)與重構(gòu)數(shù)據(jù)作差,得到含有噪聲以及粗差的殘差序列:

    是“煙富3號(hào)”的早熟株變。2015年通過貴州省品種審定委員會(huì)審定。其主要的優(yōu)良變異體現(xiàn)在成熟期上,在貴州長(zhǎng)順8月中下旬成熟。果實(shí)近圓形,果形指數(shù)0.87,平均單果重210~250克。果面平滑,稀有果粉,果皮底色淡黃,條紅,著色面積85%以上。果肉淡黃,質(zhì)地硬脆,肉質(zhì)細(xì),汁液多,風(fēng)味酸甜,香氣濃,品質(zhì)上等??扇苄怨绦挝?4.7%~15.9%。

    (8)

    式中,v為殘差項(xiàng);

    4)采用Sn估計(jì)量對(duì)殘差序列v進(jìn)行粗差探測(cè),判別式為:

    |vi-median(v)|>3Sn

    (9)

    式中,median(v)為殘差序列的中值。

    2 實(shí)驗(yàn)結(jié)果與分析

    2.1 模擬GNSS站坐標(biāo)時(shí)間序列

    為檢驗(yàn)粗差探測(cè)方法的有效性,根據(jù)GNSS單站、單分量坐標(biāo)時(shí)間序列函數(shù)模型和隨機(jī)模型構(gòu)造包含趨勢(shì)項(xiàng)、周期項(xiàng)和噪聲項(xiàng)的模擬時(shí)間序列,其函數(shù)表達(dá)式為:

    (10)

    式中,ti為觀測(cè)時(shí)間,a為測(cè)站時(shí)間序列的起始位置,b為測(cè)站運(yùn)動(dòng)的線性速度,c、d、e、f分別為測(cè)站周年運(yùn)動(dòng)和半周年運(yùn)動(dòng)的振幅,vi為GPS坐標(biāo)時(shí)間序列中的噪聲。為更加真實(shí)地模擬GPS時(shí)間序列噪聲,采用式(11)生成白噪聲和有色噪聲:

    (11)

    式中,w(ti)是均值為0、方差為1的白噪聲,f1和f2均為0.2。

    表1 模擬數(shù)據(jù)統(tǒng)計(jì)

    圖1 模擬時(shí)間序列Fig.1 Simulated time series

    2.2 模擬數(shù)據(jù)窗口長(zhǎng)度L和主成分個(gè)數(shù)P的確定

    圖2 模擬數(shù)據(jù)的RMSE隨窗口長(zhǎng)度L與主成分個(gè)數(shù)P的變化情況Fig.2 The RMSE of the simulated data varies with the length of the window and the number of principal components

    由圖2可見,模擬數(shù)據(jù)的RMSE對(duì)主成分個(gè)數(shù)P的變化比對(duì)窗口長(zhǎng)度L的變化更加敏感,在主成分個(gè)數(shù)0~10范圍內(nèi)每條曲線的RMSE都達(dá)到了最小值。經(jīng)交叉驗(yàn)證后,選取最優(yōu)窗口長(zhǎng)度L為830、主成分個(gè)數(shù)P為3。

    2.3 模擬數(shù)據(jù)結(jié)果分析

    圖3 原始序列與重構(gòu)序列Fig.3 Original sequence and reconstructed sequence

    圖4 提取的殘差序列Fig.4 The extracted residual sequence

    利用3σ法、MAD法和SSA-Sn法對(duì)模擬數(shù)據(jù)進(jìn)行粗差探測(cè)(圖5)。由圖可見,由于數(shù)據(jù)標(biāo)準(zhǔn)差σ會(huì)受到粗差的污染,不具有穩(wěn)健性,因此3σ法只能探測(cè)出偏離較大的粗差,對(duì)偏離程度較小粗差的探測(cè)效果有限。相較于3σ法,MAD法具有50%的抗差穩(wěn)健性,能夠探測(cè)出絕大部分粗差。為進(jìn)一步體現(xiàn)粗差探測(cè)的效果,進(jìn)行3種方法的探測(cè)效果統(tǒng)計(jì)(表2)。由表可知,3σ法的探測(cè)數(shù)為43、準(zhǔn)確率為38.7%;MAD法的探測(cè)數(shù)為104、準(zhǔn)確率為93.6%;SSA-Sn法的探測(cè)數(shù)為112,雖然該方法存在1個(gè)誤判,但能夠探測(cè)出所有的粗差點(diǎn),其準(zhǔn)確率相較于3σ法和MAD法分別提升61.3百分點(diǎn)和6.4百分點(diǎn),因此SSA-Sn法對(duì)于偏離程度較大或較小的粗差都具有很好的探測(cè)效果。

    圖5 3種方法的粗差探測(cè)結(jié)果Fig.5 Gross error detection results of three methods

    表2 3種方法粗差探測(cè)對(duì)比

    圖6 原始序列與重構(gòu)序列 and reconstructed sequence

    圖7 提取的殘差序列Fig.7 The extracted residual sequence

    圖8 SSA-Sn法粗差探測(cè)結(jié)果Fig.8 SSA-Sn method gross error detection results

    2.4 實(shí)測(cè)數(shù)據(jù)實(shí)驗(yàn)結(jié)果分析

    為進(jìn)一步驗(yàn)證SSA-Sn法探測(cè)粗差的有效性,使數(shù)據(jù)更加貼合實(shí)際,選用SOPAC GPS提供的原始時(shí)間序列數(shù)據(jù),選取BJFS、SHAO、URUM三個(gè)測(cè)站高程U方向上2009~2015年共計(jì)6 a的實(shí)測(cè)數(shù)據(jù)進(jìn)行處理,圖9為3個(gè)測(cè)站U方向的原始時(shí)間序列數(shù)據(jù)。由圖9(a)和圖9(b)可見,原始數(shù)據(jù)中粗差在時(shí)域和大小上分布范圍較廣。

    圖9 3個(gè)測(cè)站原始數(shù)據(jù)Fig.9 Raw data of three stations

    分別采用3σ法、MAD法和SSA-Sn法對(duì)3個(gè)測(cè)站的數(shù)據(jù)進(jìn)行粗差探測(cè)(圖10)。由BJFS站結(jié)果可知,2009~2010年的5個(gè)數(shù)據(jù)與整體數(shù)據(jù)存在明顯偏離,SSA-Sn法成功探測(cè)出2009~2010年處的粗差,但3σ法和MAD法未能全部探測(cè)出;由URUM站結(jié)果可知,3σ法和MAD法將2011年處的數(shù)據(jù)判定為粗差,但SSA-Sn法并未進(jìn)行判定。結(jié)合圖像分析后得知,URUM站2011年處的數(shù)據(jù)符合整體數(shù)據(jù)的周期趨勢(shì),不存在偏離,因此不屬于粗差。綜上所述,3σ法和MAD法均能準(zhǔn)確探測(cè)出偏離程度較大的粗差,而SSA-Sn法對(duì)于偏離程度較小的粗差數(shù)據(jù)的探測(cè)效果更加優(yōu)越。由于利用剔除粗差個(gè)數(shù)來(lái)評(píng)判粗差探測(cè)方法的優(yōu)劣具有一定的局限性,本文采用SSA重構(gòu)后的時(shí)間序列作為參照,以3種方法剔除粗差后的RMSE來(lái)評(píng)判粗差探測(cè)的效果(表3)。由表可見,經(jīng)SSA-Sn法剔除粗差后3個(gè)測(cè)站數(shù)據(jù)的RMSE分別為4.30 mm、4.03 mm和4.90 mm,均小于MAD法和3σ法處理后3個(gè)測(cè)站數(shù)據(jù)的RMSE。因此3種方法都能探測(cè)出一定數(shù)量的粗差,但SSA-Sn法的粗差探測(cè)效果最優(yōu)。

    圖10 3個(gè)測(cè)站3種方法的粗差探測(cè)結(jié)果Fig.10 Gross error detection results of three methods at three stations

    表3 3種方法的粗差探測(cè)數(shù)量與剔除后的RMSE

    3 結(jié) 語(yǔ)

    1)3σ法和MAD法的粗差探測(cè)率分別為38.7%和93.6%,而SSA-Sn法雖然存在1個(gè)誤判,卻能將模擬的粗差全部探測(cè)出來(lái),并且對(duì)于噪聲模型為白噪聲+閃爍噪聲+隨機(jī)游走噪聲這種貼合實(shí)際情況的模擬數(shù)據(jù)同樣具有適用性。

    2)采用SOPAC GPS提供的BJFS、SHAO、URUM三個(gè)測(cè)站高程U方向的原始數(shù)據(jù)進(jìn)行實(shí)驗(yàn),從探測(cè)粗差個(gè)數(shù)以及剔除粗差后數(shù)據(jù)的RMSE兩方面證明,SSA-Sn法的粗差探測(cè)效果優(yōu)于3σ法和MAD法。

    今后需要結(jié)合GNSS時(shí)間序列和奇異譜分析構(gòu)建系統(tǒng)的奇異譜分析方法,從而達(dá)到更高效的粗差探測(cè)效果。

    猜你喜歡
    估計(jì)量原始數(shù)據(jù)測(cè)站
    GNSS鐘差估計(jì)中的兩種測(cè)站選取策略分析
    GOLDEN OPPORTUNITY FOR CHINA-INDONESIA COOPERATION
    受特定變化趨勢(shì)限制的傳感器數(shù)據(jù)處理方法研究
    全球GPS測(cè)站垂向周年變化統(tǒng)計(jì)改正模型的建立
    測(cè)站分布對(duì)GPS解算ERP的影響分析
    全新Mentor DRS360 平臺(tái)借助集中式原始數(shù)據(jù)融合及直接實(shí)時(shí)傳感技術(shù)實(shí)現(xiàn)5 級(jí)自動(dòng)駕駛
    汽車零部件(2017年4期)2017-07-12 17:05:53
    淺談估計(jì)量的優(yōu)良性標(biāo)準(zhǔn)
    基于配網(wǎng)先驗(yàn)信息的諧波狀態(tài)估計(jì)量測(cè)點(diǎn)最優(yōu)配置
    負(fù)極值指標(biāo)估計(jì)量的漸近性質(zhì)
    基于GPS坐標(biāo)殘差序列的全球測(cè)站非線性變化規(guī)律統(tǒng)計(jì)
    人人妻,人人澡人人爽秒播| 久久亚洲真实| 狂野欧美白嫩少妇大欣赏| 亚洲成人免费电影在线观看| 非洲黑人性xxxx精品又粗又长| 欧美黑人欧美精品刺激| 亚洲精品日韩av片在线观看| 亚洲 国产 在线| 此物有八面人人有两片| 精品人妻偷拍中文字幕| 国产精品电影一区二区三区| 日本 欧美在线| 免费av毛片视频| 中文亚洲av片在线观看爽| 久久久精品大字幕| 最后的刺客免费高清国语| 亚洲精品成人久久久久久| 免费看光身美女| 在现免费观看毛片| 亚洲一级一片aⅴ在线观看| 久久久久国产精品人妻aⅴ院| 国产精品伦人一区二区| 最近视频中文字幕2019在线8| 免费电影在线观看免费观看| 久久人人精品亚洲av| 日韩中文字幕欧美一区二区| 中文资源天堂在线| 中文字幕人妻熟人妻熟丝袜美| 一进一出抽搐动态| 免费人成视频x8x8入口观看| 男女边吃奶边做爰视频| 国产精品无大码| 亚洲电影在线观看av| 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看日本一区| 免费黄网站久久成人精品| 国内精品久久久久久久电影| 国产精品嫩草影院av在线观看 | 少妇被粗大猛烈的视频| 亚洲成人免费电影在线观看| 少妇高潮的动态图| 中国美白少妇内射xxxbb| 成人一区二区视频在线观看| 精品久久久久久久久久久久久| 97超级碰碰碰精品色视频在线观看| 免费看美女性在线毛片视频| 五月玫瑰六月丁香| 亚洲av美国av| 亚洲av日韩精品久久久久久密| 18禁黄网站禁片午夜丰满| 久久国内精品自在自线图片| 午夜免费激情av| 天美传媒精品一区二区| 天堂√8在线中文| 久久久成人免费电影| 午夜免费激情av| 人人妻人人看人人澡| 亚洲久久久久久中文字幕| 欧美xxxx性猛交bbbb| 性色avwww在线观看| 国产精品自产拍在线观看55亚洲| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久亚洲av鲁大| 伦理电影大哥的女人| 成人欧美大片| 国产av一区在线观看免费| 国产69精品久久久久777片| 国产真实乱freesex| 日本 欧美在线| 国产精品爽爽va在线观看网站| 午夜a级毛片| 村上凉子中文字幕在线| 三级国产精品欧美在线观看| 老熟妇乱子伦视频在线观看| 少妇的逼水好多| 日韩欧美在线乱码| 很黄的视频免费| 亚洲国产精品成人综合色| 亚洲av免费在线观看| 日日啪夜夜撸| 成人国产综合亚洲| 观看免费一级毛片| 国产欧美日韩一区二区精品| 国内毛片毛片毛片毛片毛片| 丰满的人妻完整版| 免费不卡的大黄色大毛片视频在线观看 | 中出人妻视频一区二区| 午夜爱爱视频在线播放| 国产亚洲av嫩草精品影院| 日韩欧美精品免费久久| 日本撒尿小便嘘嘘汇集6| 免费看av在线观看网站| 午夜福利视频1000在线观看| 亚洲av熟女| 免费搜索国产男女视频| 欧美在线一区亚洲| 99国产精品一区二区蜜桃av| 搞女人的毛片| 女的被弄到高潮叫床怎么办 | 欧美激情国产日韩精品一区| 麻豆一二三区av精品| av专区在线播放| 国产探花在线观看一区二区| 久久人妻av系列| 国产大屁股一区二区在线视频| 亚洲成人久久爱视频| 国产精品一区www在线观看 | 免费看光身美女| 成人高潮视频无遮挡免费网站| 搡老熟女国产l中国老女人| 欧美不卡视频在线免费观看| 欧美日韩综合久久久久久 | av中文乱码字幕在线| 他把我摸到了高潮在线观看| 丰满乱子伦码专区| 亚洲国产精品成人综合色| 久久九九热精品免费| 热99re8久久精品国产| 亚洲人成伊人成综合网2020| 国产亚洲av嫩草精品影院| 日本 欧美在线| 欧美黑人巨大hd| 中文字幕熟女人妻在线| www日本黄色视频网| 国内精品美女久久久久久| 亚洲在线观看片| 蜜桃久久精品国产亚洲av| 如何舔出高潮| 亚洲欧美激情综合另类| 51国产日韩欧美| 久99久视频精品免费| 国产亚洲av嫩草精品影院| 亚洲精品日韩av片在线观看| 嫩草影视91久久| 精品一区二区免费观看| 成人av一区二区三区在线看| 国产亚洲精品久久久久久毛片| 夜夜看夜夜爽夜夜摸| 亚洲中文字幕日韩| 成人特级黄色片久久久久久久| 亚洲成人免费电影在线观看| 国产激情偷乱视频一区二区| 男女那种视频在线观看| 中文在线观看免费www的网站| 亚洲男人的天堂狠狠| 国产精品久久电影中文字幕| 内地一区二区视频在线| 99久国产av精品| 中文字幕免费在线视频6| 亚洲人成网站在线播| 亚洲国产精品久久男人天堂| 尤物成人国产欧美一区二区三区| 亚洲久久久久久中文字幕| 亚洲国产精品久久男人天堂| 一本精品99久久精品77| 九色国产91popny在线| 国产又黄又爽又无遮挡在线| 国产探花在线观看一区二区| 午夜视频国产福利| 亚洲中文日韩欧美视频| 99视频精品全部免费 在线| 亚洲一区二区三区色噜噜| 99久久久亚洲精品蜜臀av| 亚洲av中文av极速乱 | 嫩草影院入口| a在线观看视频网站| 午夜老司机福利剧场| 久久国内精品自在自线图片| 99精品久久久久人妻精品| 免费看av在线观看网站| 老司机福利观看| 联通29元200g的流量卡| 成人永久免费在线观看视频| 黄色配什么色好看| 变态另类成人亚洲欧美熟女| 国产老妇女一区| 日本三级黄在线观看| 很黄的视频免费| 麻豆av噜噜一区二区三区| 麻豆av噜噜一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕人妻熟人妻熟丝袜美| 亚洲在线自拍视频| av专区在线播放| 久久久久免费精品人妻一区二区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品98久久久久久宅男小说| 夜夜看夜夜爽夜夜摸| 精品人妻偷拍中文字幕| 日韩欧美精品v在线| 精华霜和精华液先用哪个| 91麻豆精品激情在线观看国产| 国产一区二区激情短视频| 最近中文字幕高清免费大全6 | 偷拍熟女少妇极品色| 精品无人区乱码1区二区| 国产精品永久免费网站| 国产亚洲欧美98| 色视频www国产| 韩国av在线不卡| 能在线免费观看的黄片| 好男人在线观看高清免费视频| 99国产精品一区二区蜜桃av| 免费观看人在逋| 最近中文字幕高清免费大全6 | 俄罗斯特黄特色一大片| 亚洲熟妇熟女久久| 国产成年人精品一区二区| 88av欧美| 婷婷精品国产亚洲av| 免费看光身美女| 两性午夜刺激爽爽歪歪视频在线观看| 日本熟妇午夜| 欧美绝顶高潮抽搐喷水| 国产精品乱码一区二三区的特点| 欧美人与善性xxx| 人妻制服诱惑在线中文字幕| 亚洲性久久影院| 日韩国内少妇激情av| 久久久久国产精品人妻aⅴ院| 天天一区二区日本电影三级| 免费看光身美女| 伦理电影大哥的女人| .国产精品久久| 熟女电影av网| 成人国产麻豆网| 琪琪午夜伦伦电影理论片6080| 我要搜黄色片| 中文字幕精品亚洲无线码一区| 老司机福利观看| 天堂av国产一区二区熟女人妻| 久久久久国产精品人妻aⅴ院| 老熟妇仑乱视频hdxx| 热99在线观看视频| 丰满乱子伦码专区| 中亚洲国语对白在线视频| 少妇人妻精品综合一区二区 | 在线免费观看的www视频| 日韩精品有码人妻一区| 琪琪午夜伦伦电影理论片6080| 亚洲精华国产精华精| 国产成人av教育| 亚洲人成网站在线播| 免费观看精品视频网站| 伊人久久精品亚洲午夜| 国内精品美女久久久久久| 一区福利在线观看| 99热这里只有是精品50| 免费av观看视频| 伦理电影大哥的女人| 91av网一区二区| 欧美一区二区亚洲| 91麻豆av在线| 美女被艹到高潮喷水动态| 久久香蕉精品热| avwww免费| 成人特级av手机在线观看| videossex国产| 99久国产av精品| 人妻丰满熟妇av一区二区三区| 欧美日本视频| 亚洲电影在线观看av| 深夜a级毛片| 欧美精品国产亚洲| 亚洲国产精品成人综合色| 国产成人aa在线观看| 直男gayav资源| 韩国av一区二区三区四区| 熟女电影av网| 欧美日本亚洲视频在线播放| 五月玫瑰六月丁香| 亚洲七黄色美女视频| 久久久国产成人免费| 日韩欧美在线二视频| netflix在线观看网站| 成人特级黄色片久久久久久久| 身体一侧抽搐| 最好的美女福利视频网| 欧美黑人欧美精品刺激| 99国产精品一区二区蜜桃av| 亚洲人与动物交配视频| 99久久成人亚洲精品观看| 日日摸夜夜添夜夜添av毛片 | 国产伦精品一区二区三区视频9| 久久久久久久午夜电影| 久久久久九九精品影院| 免费看日本二区| 久久久久国内视频| 久久99热6这里只有精品| 精品人妻一区二区三区麻豆 | 久久久久九九精品影院| 久久草成人影院| 真人做人爱边吃奶动态| 亚洲国产高清在线一区二区三| av国产免费在线观看| 丰满乱子伦码专区| 久久久久九九精品影院| 欧美成人免费av一区二区三区| 看黄色毛片网站| 亚洲第一区二区三区不卡| 婷婷色综合大香蕉| 日韩强制内射视频| 亚洲四区av| 最近在线观看免费完整版| 99在线视频只有这里精品首页| 国内久久婷婷六月综合欲色啪| 搞女人的毛片| 此物有八面人人有两片| 天堂动漫精品| 成人无遮挡网站| 成人亚洲精品av一区二区| 欧美日本视频| 亚洲专区中文字幕在线| 国产精品,欧美在线| 久久6这里有精品| 国产亚洲欧美98| 久久久久久久久中文| 亚洲成人久久性| 色吧在线观看| 亚洲精华国产精华精| 久久精品国产99精品国产亚洲性色| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 国产精品久久久久久精品电影| 免费观看人在逋| 一进一出抽搐动态| 久久久精品大字幕| 麻豆成人av在线观看| 成人国产麻豆网| 亚洲七黄色美女视频| 成年免费大片在线观看| 国产免费av片在线观看野外av| 国产精品99久久久久久久久| 老熟妇仑乱视频hdxx| 亚洲自偷自拍三级| videossex国产| 99热这里只有精品一区| 两性午夜刺激爽爽歪歪视频在线观看| 久久久午夜欧美精品| 国国产精品蜜臀av免费| 久久久国产成人免费| 黄色欧美视频在线观看| 在线天堂最新版资源| 国产精品久久久久久精品电影| 久久久久久久精品吃奶| 亚洲欧美精品综合久久99| 中国美白少妇内射xxxbb| 久久精品人妻少妇| 床上黄色一级片| 俺也久久电影网| 日韩欧美国产在线观看| 国产高清不卡午夜福利| 婷婷六月久久综合丁香| 99国产精品一区二区蜜桃av| 特大巨黑吊av在线直播| 国产亚洲精品av在线| 男女啪啪激烈高潮av片| 国产精品三级大全| 午夜爱爱视频在线播放| 老司机福利观看| 特级一级黄色大片| 国产男靠女视频免费网站| 中文字幕熟女人妻在线| xxxwww97欧美| 午夜日韩欧美国产| 国内少妇人妻偷人精品xxx网站| 亚洲av不卡在线观看| 久久国产精品人妻蜜桃| 亚洲,欧美,日韩| 日韩欧美 国产精品| 久久中文看片网| 中国美女看黄片| 亚洲男人的天堂狠狠| 国产精品乱码一区二三区的特点| 神马国产精品三级电影在线观看| 床上黄色一级片| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 九九在线视频观看精品| 国产高清激情床上av| 女人十人毛片免费观看3o分钟| 极品教师在线视频| 又黄又爽又免费观看的视频| eeuss影院久久| 成人午夜高清在线视频| 亚洲av不卡在线观看| 日韩精品中文字幕看吧| 在线观看一区二区三区| 亚洲黑人精品在线| 别揉我奶头 嗯啊视频| 国产黄片美女视频| 国产主播在线观看一区二区| 欧美黑人巨大hd| 99久久成人亚洲精品观看| 在线观看66精品国产| 99久久中文字幕三级久久日本| 最后的刺客免费高清国语| 观看免费一级毛片| 国内精品美女久久久久久| 久久久久久国产a免费观看| 狂野欧美白嫩少妇大欣赏| 人人妻人人看人人澡| 精品人妻1区二区| 国产黄色小视频在线观看| 欧美另类亚洲清纯唯美| 国产精品98久久久久久宅男小说| 久久久久久久久久黄片| a在线观看视频网站| 成年女人毛片免费观看观看9| 日韩亚洲欧美综合| 欧美日本视频| 午夜免费成人在线视频| 免费搜索国产男女视频| 99久久成人亚洲精品观看| 国产精品亚洲美女久久久| 一级黄色大片毛片| 很黄的视频免费| 国产精品久久久久久久久免| 日韩精品中文字幕看吧| 久久久久免费精品人妻一区二区| 搡老妇女老女人老熟妇| 久久久久性生活片| 最近最新中文字幕大全电影3| 变态另类成人亚洲欧美熟女| 久久久久久国产a免费观看| 国产精品久久久久久久久免| 日韩精品中文字幕看吧| 国产高清视频在线观看网站| 国产精品无大码| 色哟哟哟哟哟哟| 亚洲18禁久久av| 日韩高清综合在线| 久久精品国产亚洲av天美| 久久欧美精品欧美久久欧美| 亚洲五月天丁香| 国产色婷婷99| 欧美+亚洲+日韩+国产| 亚洲av免费在线观看| 午夜免费激情av| 亚洲天堂国产精品一区在线| 日本免费一区二区三区高清不卡| 午夜视频国产福利| 日本 av在线| 在现免费观看毛片| 国产欧美日韩精品一区二区| 久久久久久久午夜电影| 不卡视频在线观看欧美| 最后的刺客免费高清国语| 身体一侧抽搐| 亚洲成a人片在线一区二区| 婷婷色综合大香蕉| 最近最新中文字幕大全电影3| 国产国拍精品亚洲av在线观看| 大又大粗又爽又黄少妇毛片口| a级毛片免费高清观看在线播放| 日韩在线高清观看一区二区三区 | 久久久久久久久中文| 亚洲成人免费电影在线观看| 国产人妻一区二区三区在| 好男人在线观看高清免费视频| 国产精品三级大全| 美女免费视频网站| 免费看a级黄色片| 久久精品久久久久久噜噜老黄 | 一边摸一边抽搐一进一小说| 少妇人妻一区二区三区视频| 国产亚洲精品av在线| 国产精品一区二区免费欧美| 婷婷精品国产亚洲av在线| 一个人看视频在线观看www免费| 夜夜夜夜夜久久久久| av专区在线播放| 精品福利观看| 亚洲午夜理论影院| 久久99热6这里只有精品| 老司机福利观看| 午夜视频国产福利| 国产精品一区二区三区四区免费观看 | 精品人妻熟女av久视频| 亚洲av电影不卡..在线观看| 波多野结衣巨乳人妻| 免费看a级黄色片| or卡值多少钱| aaaaa片日本免费| 黄色配什么色好看| www日本黄色视频网| 久久精品国产亚洲av天美| 悠悠久久av| 欧美xxxx性猛交bbbb| 日本黄色视频三级网站网址| 热99re8久久精品国产| 日本成人三级电影网站| 欧美色视频一区免费| 精品人妻1区二区| 真人一进一出gif抽搐免费| 精品人妻偷拍中文字幕| 精品无人区乱码1区二区| 免费在线观看影片大全网站| 亚洲一区高清亚洲精品| 久久亚洲精品不卡| 成年免费大片在线观看| 婷婷六月久久综合丁香| 中文字幕久久专区| 国产男人的电影天堂91| 中国美白少妇内射xxxbb| 日日撸夜夜添| 尤物成人国产欧美一区二区三区| 变态另类丝袜制服| 美女高潮的动态| 草草在线视频免费看| 黄色欧美视频在线观看| 亚洲自偷自拍三级| 尾随美女入室| 国产精品久久视频播放| 老熟妇乱子伦视频在线观看| 久久久色成人| 2021天堂中文幕一二区在线观| 一进一出抽搐动态| 国产精品无大码| 国产真实乱freesex| 久久国产乱子免费精品| 国产 一区 欧美 日韩| 日韩欧美一区二区三区在线观看| 美女被艹到高潮喷水动态| 一夜夜www| 国产成人一区二区在线| 一边摸一边抽搐一进一小说| 大型黄色视频在线免费观看| 51国产日韩欧美| av国产免费在线观看| 麻豆成人午夜福利视频| 欧美+日韩+精品| 婷婷六月久久综合丁香| 中文字幕av在线有码专区| 国产精品久久久久久av不卡| 国产色爽女视频免费观看| 深爱激情五月婷婷| 三级男女做爰猛烈吃奶摸视频| 乱系列少妇在线播放| www日本黄色视频网| 午夜精品一区二区三区免费看| 日韩欧美免费精品| 久久精品91蜜桃| xxxwww97欧美| 丝袜美腿在线中文| 老司机福利观看| 波多野结衣高清作品| 真实男女啪啪啪动态图| 偷拍熟女少妇极品色| 国产精品98久久久久久宅男小说| 最近最新免费中文字幕在线| 国产爱豆传媒在线观看| 嫁个100分男人电影在线观看| 人人妻,人人澡人人爽秒播| 内地一区二区视频在线| 一级a爱片免费观看的视频| 久久久久久久亚洲中文字幕| 日韩欧美免费精品| 不卡视频在线观看欧美| 精品欧美国产一区二区三| 深夜a级毛片| 欧美色欧美亚洲另类二区| 色综合亚洲欧美另类图片| 91在线观看av| 极品教师在线免费播放| 一区二区三区激情视频| 成年免费大片在线观看| av国产免费在线观看| 窝窝影院91人妻| 免费搜索国产男女视频| 三级男女做爰猛烈吃奶摸视频| 在线看三级毛片| 色吧在线观看| 欧美一级a爱片免费观看看| 校园人妻丝袜中文字幕| 国内精品美女久久久久久| 久久午夜亚洲精品久久| 日韩欧美国产一区二区入口| 毛片女人毛片| 亚洲最大成人中文| 国产欧美日韩一区二区精品| 国产精品一区二区性色av| 国产精品一区www在线观看 | 高清日韩中文字幕在线| 直男gayav资源| 亚洲最大成人av| 日本成人三级电影网站| 亚洲精品久久国产高清桃花| 亚洲精品一区av在线观看| 日本一本二区三区精品| 国产精品一区二区三区四区久久| 伊人久久精品亚洲午夜| av女优亚洲男人天堂| 亚洲电影在线观看av| 国产精品伦人一区二区| 久久久久久久久久黄片| 欧美激情国产日韩精品一区| 欧美在线一区亚洲| 国产乱人视频| 免费看av在线观看网站| avwww免费| 亚洲三级黄色毛片| 国产男人的电影天堂91| 国语自产精品视频在线第100页| 国产成人av教育| 91午夜精品亚洲一区二区三区 | 国产私拍福利视频在线观看| 观看免费一级毛片| 狠狠狠狠99中文字幕| 国产爱豆传媒在线观看| 麻豆国产av国片精品| 神马国产精品三级电影在线观看| 中文在线观看免费www的网站| 色综合站精品国产| 特级一级黄色大片|