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

    抗差估計的選權(quán)迭代法分析與比較

    2010-11-13 08:03:48陳西強黃張裕
    測繪工程 2010年4期
    關鍵詞:抗差權(quán)函數(shù)迭代法

    陳西強,黃張裕

    (河海大學 測繪科學與工程系,江蘇 南京 210098)

    抗差估計的選權(quán)迭代法分析與比較

    陳西強,黃張裕

    (河海大學 測繪科學與工程系,江蘇 南京 210098)

    在抗差穩(wěn)健估計中,選權(quán)迭代法應用最廣,較易理解,M估計法計算簡單。基于M估計的基本理論,深入分析Huber、Hampel、IGG 3種常用的選權(quán)迭代法,并利用水準網(wǎng)數(shù)據(jù),分析這3種常用選權(quán)迭代法的適用范圍和相應的抗差效果,表明IGG法比其他兩種的抗粗差效果要好,且易實現(xiàn)。

    抗差穩(wěn)健估計;選權(quán)迭代法;Huber法;Hampel法;IGG法

    在抗差穩(wěn)健估計中,最常用的為M估計法,其估計方法有許多種,其中選權(quán)迭代法應用最廣,較易理解,計算簡單。對于選權(quán)迭代法,理論研究相對比較多,但對其中的參數(shù)選取,在理論上有所研究,實際應用研究上存在不夠深入和不足的問題[1-2]。本文在分析M穩(wěn)健估計基本理論的基礎上,討論了3種常用的選權(quán)迭代法。選用水準網(wǎng)實例,深入分析3種選權(quán)迭代法的初始權(quán)選取和條件參數(shù)選取,并從縱橫兩個方面來進行比較、分析這3種常用選權(quán)迭代法的適用范圍和相應抗差效果,表明IGG法比其他兩種方法的抗粗差效果要好,且易實現(xiàn)。

    1 M穩(wěn)健估計

    設有參數(shù)向量X是未知的非隨機量,為了估計未知向量值,進行n次觀測得到向量L的觀測值l,由極大似然估計有[3]

    式中:f是隨機量L的密度函數(shù)。

    Huber提出用ρ(li,x)代替函數(shù)-lnf(li,^x),使其定義廣義化,于是有

    式中:φ(li,^x)=?ρ(li,x)/?x。

    由式(2)對參數(shù)x進行估計,就是廣義極大似然估計,簡稱M估計。測量平差中,觀測量L的殘差為V,M估計的函數(shù)ρ可取為ρ(vi),即M估計的準則為

    通常殘差V為未知數(shù)的函數(shù),將式(2)對未知數(shù)x求導,取極值點,則有

    考慮到誤差方程vi=ai^x-li,ai為A的第i行向量,則有

    令權(quán)函數(shù)pi(vi)=ρ′(vi)/vi,即得

    用矩陣表示為

    式(7)類似最小二乘估計的法方程形式,用權(quán)函數(shù)矩陣P(V)=diag(p1(vi),p2(vi),…,pn(vi))代替觀測權(quán)陣P。為此,可將穩(wěn)健估計的選權(quán)迭代法歸結(jié)為如下模型:

    誤差方程與權(quán)函數(shù)為

    估計準則為

    將平差模型轉(zhuǎn)換為最小二乘估計的模型,采用類似最小二乘估計。只是權(quán)函數(shù)P(V)是殘差的函數(shù),計算前V未知,只能通過給其賦予一定的初值,采用迭代法估計參數(shù)^x。具體計算過程如下:

    1)列誤差方程,令各觀測權(quán)函數(shù)初值均為1,即令p1(vi)=p2(vi)=…=pn(vi)=1。

    2)解算法方程,得出^x和V的第1次估值^x(1)=(ATA)-1ATl,V(1)=A^x(1)-l。

    3)由V(1)確定各觀測權(quán)函數(shù)pi(vi),再解算法方程,通過迭代計算,直到前后兩次解的差值符合要求為止。

    4)最后結(jié)果為

    隨著函數(shù)ρ的選取不同,可以構(gòu)成權(quán)函數(shù)的多種不同的形式,通常權(quán)函數(shù)是一個在平差過程中隨改正數(shù)變化的量,經(jīng)過多次迭代,使含有粗差的奇異觀測的權(quán)函數(shù)為0(或者接近于0),而相應的殘差值在很大程度上反映了其粗差值[4-5]。

    2 幾種常用的選權(quán)迭代法

    2.1 Huber法

    Huber法的極值估計準則是假設觀測值在不利分布情況下求最優(yōu)估計[6]。取M估計的函數(shù)

    由Huber權(quán)函數(shù)可以看出,Huber估計的權(quán)函數(shù)是單調(diào)的,它在c處的斜率有突然變化。當所有改正數(shù)均在-c和c之間時,Huber估計就是經(jīng)典最小二乘法估計[7]。而當部分改正數(shù)大于c時,其p(v)與改正數(shù)成反比,v愈大,對應的p(v)愈小,與此相應該觀測值對參數(shù)估計的影響也愈小。

    2.2 Hampel法

    對于Hampel法,取M估計的函數(shù)為3段截尾函數(shù)

    求得權(quán)函數(shù)pi(vi)為

    式中:a、b、c為3個系數(shù),可適當選取,它們均為標準差的倍數(shù)。

    2.3 IGG法(周文江法)

    IGG法是由中國科學院測量與地球物理研究所周文江教授最初提出的,取M估計的函數(shù)[8]

    求得權(quán)函數(shù)pi(vi)為

    IGG法屬于有淘汰區(qū)的M估計,權(quán)因子之間變化比較平緩,k值一般取一個很小的數(shù)。

    3 實例分析

    水準路線如圖1所示。

    圖1 水準路線

    如表1所示水準網(wǎng),A、B兩點的已知高程分別為HA=5.000 m,HB=6.000 m,觀測了7條水準線路,表中帶*號的表示有粗差的觀測值,對h2賦以粗差0.050 m,h6賦以粗差-0.030 m。根據(jù)水準路線觀測數(shù)據(jù),按最小二乘原理,在不含有粗差和含有粗差的條件下分別解算,如表2所示。

    表1 水準網(wǎng)觀測數(shù)據(jù)

    表2 基于最小二乘解算結(jié)果

    3.1 3種方法分析

    對于 Huber法分析c選取對估計結(jié)果影響。在迭代過程中,σ值隨每次迭代后算出的觀測改正值來確定。在本次數(shù)據(jù)分析中,c分別取 0.4σ、0.8σ、σ、1.2σ、1.6σ、2.0σ、2.5σ,試驗數(shù)據(jù)結(jié)果如表3所示。

    表3 c選取對估計結(jié)果的影響

    從表3的數(shù)據(jù)可以知道,隨著c取值的增大,總體上中誤差越來越大,迭代次數(shù)越來越多,到達一定的時候,迭代效果失效。而且當c為一倍中誤差的時候,Vh2和Vh6的改正值是最能反映出含有粗差的數(shù)據(jù),但總體效果不是很理想。

    對于Hampel法分析a、b、c的取值對估計結(jié)果影響。a、b、c3個系數(shù)分別取4組值進行比較分析,即取a=0.9σ,b=1.8σ,c=2.7σ;a=σ,b=2σ,c=3σ; a=1.1σ,b=2.2σ,c=3.3σ和a=1.2σ,b=2.4σ,c= 3.6σ,試驗數(shù)據(jù)結(jié)果如表4所示。

    表4 a、b、c取值對估計結(jié)果的影響

    從表4數(shù)據(jù)分析可知,在a、b、c三者的相互間隔增大的同時,方案1、2、3的估計結(jié)果沒有大的變化,只是中誤差和迭代次數(shù)都有所增大。當三者之間的間隔超出了一定的范圍時,必將導致估計失效。試驗分析得出間隔應在0.9~1.1之間,而且當間隔太小時,會導致求解過程中的權(quán)陣單極化。

    對于IGG法分析k取值對估計效果影響。分別取k為0.001、0.000 1、0.000 01,試驗所得數(shù)據(jù)如表5所示。

    表5 k取值對估計效果的影響 mm

    從表5的數(shù)據(jù)分析可知,只要k值在一定的范圍之內(nèi),即對估計的效果沒有影響。權(quán)函數(shù)中的條件是隨著每次迭代得到的中誤差,在整個迭代過程中是變化的,而且越來越低,隨著滿足迭代條件而終止。

    另外,對于3種選權(quán)迭代方法,分別分析初始權(quán)選取對估計結(jié)果的影響。在試驗分析中分別采用了初始權(quán)都為1和初始權(quán)按距離定權(quán)的2種方案進行數(shù)據(jù)分析,3種迭代方法計算結(jié)果完全相同,數(shù)據(jù)結(jié)果如表6所示。

    表6 初始權(quán)選取對估計結(jié)果的影響 mm

    從表6的數(shù)據(jù)分析可以看出,初始權(quán)的選取對3種方法的估計結(jié)果都沒有影響。

    3.2 3種方法比較分析

    現(xiàn)對Huber、Hampel、IGG方法作橫向比較,由這3種方法分別對這個已知含有粗差的水準網(wǎng)觀測數(shù)據(jù)進行平差處理,得到的結(jié)果如表7所示。

    表7 不同選權(quán)迭代法比較

    從表7的數(shù)據(jù)分析可知,在粗差較大的情況下,在迭代次數(shù)和最后的精度及計算結(jié)果綜合因素來看,IGG法比 Hampel法和Huber法效果要好。在抗粗差及剔除粗差方面,IGG和 Hampel都能很好地剔除由于h2帶來的粗差0.050 m和h6已賦的粗差-0.030 m,而Huber法相對來說要差點,而且在Huber法中的c值的選取好壞直接影響到抗粗差的能力。

    4 結(jié) 論

    1)穩(wěn)健估計中的選權(quán)迭代法在初始化權(quán)方面都是相同的,即估計結(jié)果不受初始權(quán)的取值影響,且在選權(quán)迭代過程中,對于選權(quán)的標準或條件系數(shù)都是隨著每次迭代完后重新計算的中誤差相關,說明3種選權(quán)迭代的條件系數(shù)是動態(tài)的。

    2)對于 Huber法,隨著c取值的增大,總體上中誤差越來越大,迭代次數(shù)越來越多,到達一定的時候,迭代失效,而且當c為1倍中誤差時,Vh2和Vh6的改正值是最能反映出含有粗差的數(shù)據(jù),但效果不是很理想,因此,在含粗差較大時,建議 Huber法中c的取值范圍為中誤差的0.8~1.2之間。對于Hampel法,a、b、c三者的取值在一定范圍內(nèi)相互間隔增大時,估計結(jié)果沒有大的變化,只是中誤差和迭代次數(shù)都有所增大,當三者之間的間隔超出了一定的范圍,將導致估計失效,試驗結(jié)果表明間隔應在0.9~1.1之間,而且當間隔太小時,會導致求解過程中的權(quán)陣單極化。在IGG法分析中,k的取值足夠小就能滿足要求,如本文試驗k值取0.001~0.000 01。

    3)在這3種穩(wěn)健估計的選權(quán)迭代法中,IGG法比其他兩種的抗粗差效果要好,且易實現(xiàn);在不清楚數(shù)據(jù)分布的情況下,建議采用IGG法。

    [1]陶本藻.穩(wěn)健估計的應用問題[J].地礦測繪,2000(1): 1-3.

    [2]楊元喜.參數(shù)M估計及驗后精度綜述[J].測繪通報, 1991(6):6-10.

    [3]劉大杰 ,陶本藻.實用測量數(shù)據(jù)處理方法[M].北京:測繪出版社,2000:52-58.

    [4]吳富梅.多種權(quán)函數(shù)固定臨界值與可變臨界值抗差估計的比較[J].測繪工程,2006,15(3):19-22.

    [5]李浩軍,唐詩華,黃杰.抗差估計中幾種選權(quán)迭代法常數(shù)選取的探討[J].測繪科學,2006,31(6):70-72.

    [6]周江文.經(jīng)典誤差理論與抗差估計[A].抗差估計論文集[C].北京:測繪出版社,1992:1-6.

    [7]姚宜斌.基于等價分析法的穩(wěn)健估計權(quán)函數(shù)設計[J].測繪工程,2001,10(2):29-31.

    [8]楊元喜,宋力杰,徐天河.大地測量相關觀測抗差估計理論[J].測繪學報,2002,31(2):95-99.

    Analysisand comparison of the selecting weight iteration method in resistance robust estimation

    CHEN Xi-qiang,HUANG Zhang-yu
    (Department of Survey Science and Engineering,Hohai University,Nanjing 210098,China)

    In the robust estimation,selecting weight iteration method is the most w idely understood,more easily and simp le calculation of M estimation.Based on the basic theory of M estimation method,this paper deep ly analyzes three common selecting w eight iteration methodsw hich are Huber、Hampel、IGGmethods,and analyzes the app licable scope and co rresponding poo r effect on the three commonly selecting w eight iteration methods by using the data of the netwo rk level.It show s that the exco rresponding poo r effect of the IGG is better than other two methods,and easy to achieve results.

    resistance robust estimation;selecting weight iteration method;Hubermethod;Hampelmethod;IGGmethod.

    P207

    A

    1006-7949(2010)04-0008-04

    2009-09-03

    陳西強(1986-),男,碩士研究生.

    [責任編輯:劉文霞]

    猜你喜歡
    抗差權(quán)函數(shù)迭代法
    基于改進權(quán)函數(shù)的探地雷達和無網(wǎng)格模擬檢測混凝土結(jié)構(gòu)空洞缺陷工程中的數(shù)學問題
    迭代法求解一類函數(shù)方程的再研究
    一類廣義的十次Freud-型權(quán)函數(shù)
    異徑電磁流量傳感器權(quán)函數(shù)分布規(guī)律研究*
    改善單頻PPP參數(shù)收斂速度的抗差估計方法
    迭代法求解約束矩陣方程AXB+CYD=E
    預條件SOR迭代法的收斂性及其應用
    地形簡化對DEM不確定性的抗差性研究
    基于抗差最小均方估計的輸電線路參數(shù)辨識
    電力建設(2015年2期)2015-07-12 14:15:59
    兩類ω-超廣義函數(shù)空間的結(jié)構(gòu)表示
    国产精品国产三级国产专区5o| 日韩欧美精品免费久久| 卡戴珊不雅视频在线播放| 春色校园在线视频观看| 国产精品嫩草影院av在线观看| 人妻少妇偷人精品九色| 免费黄色在线免费观看| 黄色 视频免费看| 新久久久久国产一级毛片| 免费黄网站久久成人精品| 欧美成人精品欧美一级黄| 午夜精品国产一区二区电影| 久久人人97超碰香蕉20202| 乱人伦中国视频| 一个人免费看片子| 一级爰片在线观看| 极品少妇高潮喷水抽搐| 精品一区二区三卡| 另类亚洲欧美激情| 亚洲,欧美,日韩| 夫妻性生交免费视频一级片| 日韩中字成人| 大片免费播放器 马上看| 欧美 日韩 精品 国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男人操女人黄网站| 日日爽夜夜爽网站| 18禁观看日本| 国产在线免费精品| 婷婷色综合大香蕉| 亚洲国产精品一区三区| 建设人人有责人人尽责人人享有的| 少妇的逼水好多| 成人毛片a级毛片在线播放| 日韩一区二区三区影片| 精品国产一区二区三区四区第35| 少妇的逼水好多| 欧美xxⅹ黑人| av在线观看视频网站免费| 亚洲精品乱久久久久久| 高清不卡的av网站| 一级a爱视频在线免费观看| 国产欧美日韩一区二区三区在线| 午夜福利一区二区在线看| 90打野战视频偷拍视频| 亚洲成人av在线免费| 亚洲av欧美aⅴ国产| 免费人妻精品一区二区三区视频| 日韩伦理黄色片| 成人国语在线视频| 国产精品三级大全| 在线观看美女被高潮喷水网站| 中文字幕色久视频| 欧美另类一区| 国产免费视频播放在线视频| av不卡在线播放| 热re99久久精品国产66热6| 欧美成人午夜免费资源| 国产xxxxx性猛交| 赤兔流量卡办理| 伦精品一区二区三区| videosex国产| 人体艺术视频欧美日本| 母亲3免费完整高清在线观看 | 极品人妻少妇av视频| 成人国语在线视频| 久久午夜综合久久蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 午夜91福利影院| 精品亚洲乱码少妇综合久久| av.在线天堂| 最黄视频免费看| 永久网站在线| 狠狠精品人妻久久久久久综合| 如何舔出高潮| 晚上一个人看的免费电影| 高清在线视频一区二区三区| 成年动漫av网址| 精品人妻熟女毛片av久久网站| 亚洲成人av在线免费| 91午夜精品亚洲一区二区三区| 欧美日韩一级在线毛片| 亚洲av男天堂| 国产精品一区二区在线不卡| 丰满少妇做爰视频| 女人被躁到高潮嗷嗷叫费观| 美女午夜性视频免费| 国产成人精品在线电影| 亚洲国产av影院在线观看| 亚洲国产欧美网| 精品亚洲成国产av| 久久免费观看电影| 久久 成人 亚洲| 国产在线一区二区三区精| 丝袜美腿诱惑在线| 搡老乐熟女国产| 亚洲精品在线美女| 九草在线视频观看| 日韩中文字幕视频在线看片| 麻豆精品久久久久久蜜桃| 亚洲,欧美,日韩| 日韩伦理黄色片| 午夜老司机福利剧场| 欧美精品高潮呻吟av久久| 久久99一区二区三区| 欧美中文综合在线视频| 国产一级毛片在线| av.在线天堂| 午夜久久久在线观看| 人妻系列 视频| 人人妻人人爽人人添夜夜欢视频| 国产成人免费观看mmmm| 一二三四在线观看免费中文在| 宅男免费午夜| 久久久精品国产亚洲av高清涩受| 久久久久网色| 亚洲第一青青草原| 日日撸夜夜添| 一边摸一边做爽爽视频免费| 成人亚洲精品一区在线观看| 在线观看三级黄色| 亚洲av综合色区一区| 国产成人精品福利久久| 欧美日韩亚洲国产一区二区在线观看 | 国产成人精品一,二区| 精品福利永久在线观看| 最近中文字幕高清免费大全6| 999久久久国产精品视频| 国产av精品麻豆| av网站在线播放免费| 亚洲五月色婷婷综合| 亚洲成av片中文字幕在线观看 | www日本在线高清视频| 九九爱精品视频在线观看| 丝袜喷水一区| 在线观看www视频免费| 亚洲精品中文字幕在线视频| 91成人精品电影| 色94色欧美一区二区| 超色免费av| 成年女人毛片免费观看观看9 | 久久精品久久精品一区二区三区| 人妻少妇偷人精品九色| 999久久久国产精品视频| 两个人免费观看高清视频| 免费观看av网站的网址| 最新中文字幕久久久久| 国产精品一区二区在线不卡| 午夜激情av网站| 亚洲三级黄色毛片| 欧美日韩一区二区视频在线观看视频在线| 91国产中文字幕| 国产成人免费无遮挡视频| 国产日韩欧美亚洲二区| 大香蕉久久成人网| 一本大道久久a久久精品| 国产午夜精品一二区理论片| av片东京热男人的天堂| 人妻系列 视频| 国产免费现黄频在线看| 搡老乐熟女国产| 国产亚洲午夜精品一区二区久久| 精品国产超薄肉色丝袜足j| 久久鲁丝午夜福利片| 爱豆传媒免费全集在线观看| 麻豆精品久久久久久蜜桃| 人人妻人人澡人人爽人人夜夜| 一级片'在线观看视频| 三级国产精品片| 看十八女毛片水多多多| 国产人伦9x9x在线观看 | 亚洲伊人久久精品综合| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人 | 欧美另类一区| 成人影院久久| 欧美成人午夜免费资源| 免费人妻精品一区二区三区视频| 久久午夜福利片| 一区福利在线观看| 亚洲色图 男人天堂 中文字幕| av电影中文网址| 午夜福利一区二区在线看| 久久精品国产自在天天线| 不卡av一区二区三区| 尾随美女入室| 一级片'在线观看视频| 老汉色∧v一级毛片| 人人妻人人澡人人爽人人夜夜| 久久精品久久久久久噜噜老黄| 女人精品久久久久毛片| 精品人妻熟女毛片av久久网站| 寂寞人妻少妇视频99o| 日本wwww免费看| 中国国产av一级| 王馨瑶露胸无遮挡在线观看| 黄色 视频免费看| 18禁国产床啪视频网站| 久久99一区二区三区| 一级毛片 在线播放| 国产精品久久久久久av不卡| www.自偷自拍.com| 伦理电影大哥的女人| 久久久精品国产亚洲av高清涩受| 夫妻性生交免费视频一级片| 中文天堂在线官网| 精品一区二区免费观看| 国产综合精华液| 亚洲欧美日韩另类电影网站| 在线精品无人区一区二区三| 搡女人真爽免费视频火全软件| 亚洲国产最新在线播放| 国产爽快片一区二区三区| 亚洲精品美女久久久久99蜜臀 | 这个男人来自地球电影免费观看 | 伊人久久国产一区二区| 亚洲精品日韩在线中文字幕| 不卡av一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一级片'在线观看视频| 美国免费a级毛片| 老熟女久久久| 国产男女超爽视频在线观看| 日韩不卡一区二区三区视频在线| 嫩草影院入口| 免费黄网站久久成人精品| 亚洲中文av在线| 久久精品国产自在天天线| 国产1区2区3区精品| 亚洲精品国产av蜜桃| 在线天堂中文资源库| 99香蕉大伊视频| 日日摸夜夜添夜夜爱| 老熟女久久久| 亚洲欧美精品自产自拍| www.熟女人妻精品国产| 性色av一级| 精品国产一区二区三区久久久樱花| 久久ye,这里只有精品| 久久久国产欧美日韩av| 69精品国产乱码久久久| 精品国产超薄肉色丝袜足j| 国产免费福利视频在线观看| 9191精品国产免费久久| 精品人妻熟女毛片av久久网站| 日韩av免费高清视频| 久久久久久人人人人人| 丝袜人妻中文字幕| 下体分泌物呈黄色| 亚洲av电影在线进入| 国产一区二区 视频在线| 午夜影院在线不卡| 女人精品久久久久毛片| 国产精品蜜桃在线观看| 免费人妻精品一区二区三区视频| 亚洲欧美精品综合一区二区三区 | 91在线精品国自产拍蜜月| 一区二区日韩欧美中文字幕| 精品午夜福利在线看| 亚洲精品国产av蜜桃| freevideosex欧美| 97人妻天天添夜夜摸| 午夜福利在线免费观看网站| 啦啦啦在线免费观看视频4| 热99久久久久精品小说推荐| 亚洲一级一片aⅴ在线观看| 美女视频免费永久观看网站| av在线观看视频网站免费| 色婷婷久久久亚洲欧美| 国产熟女欧美一区二区| 欧美日韩国产mv在线观看视频| 伦精品一区二区三区| 午夜久久久在线观看| 国产野战对白在线观看| 97精品久久久久久久久久精品| 成年动漫av网址| 成人国产麻豆网| 一边摸一边做爽爽视频免费| 日韩一本色道免费dvd| 亚洲精品国产av成人精品| 大香蕉久久网| 久久久久久久久免费视频了| 欧美激情 高清一区二区三区| 777米奇影视久久| 国产在线一区二区三区精| 欧美日韩视频精品一区| 最近2019中文字幕mv第一页| 你懂的网址亚洲精品在线观看| 国精品久久久久久国模美| 国产不卡av网站在线观看| 久久久久国产一级毛片高清牌| 久久久久久人妻| 欧美精品一区二区免费开放| 欧美精品亚洲一区二区| 成人毛片a级毛片在线播放| 午夜福利影视在线免费观看| 美女脱内裤让男人舔精品视频| 亚洲国产日韩一区二区| 黑人欧美特级aaaaaa片| 久久久欧美国产精品| 午夜激情久久久久久久| 久久国产精品男人的天堂亚洲| 亚洲美女视频黄频| 亚洲精品国产av蜜桃| 欧美精品一区二区免费开放| 久久久久久免费高清国产稀缺| 久久久亚洲精品成人影院| 99久久精品国产国产毛片| 亚洲经典国产精华液单| 美女福利国产在线| 日韩熟女老妇一区二区性免费视频| 亚洲情色 制服丝袜| 久久这里有精品视频免费| 国产一区亚洲一区在线观看| av在线老鸭窝| 久久久久视频综合| 国产一区有黄有色的免费视频| 日韩制服骚丝袜av| 老熟女久久久| 极品人妻少妇av视频| 三上悠亚av全集在线观看| 亚洲精品,欧美精品| 男人爽女人下面视频在线观看| 不卡视频在线观看欧美| 久久久久国产一级毛片高清牌| 国产精品一二三区在线看| 日本免费在线观看一区| 欧美国产精品一级二级三级| 久久久久久免费高清国产稀缺| 黑丝袜美女国产一区| 自拍欧美九色日韩亚洲蝌蚪91| 国产无遮挡羞羞视频在线观看| 大陆偷拍与自拍| √禁漫天堂资源中文www| 999精品在线视频| 日本免费在线观看一区| 制服人妻中文乱码| 国产在线一区二区三区精| 中文天堂在线官网| 校园人妻丝袜中文字幕| 熟妇人妻不卡中文字幕| 在线看a的网站| 在线观看国产h片| 久久精品国产鲁丝片午夜精品| 视频在线观看一区二区三区| 丝袜美足系列| 精品国产超薄肉色丝袜足j| 国产无遮挡羞羞视频在线观看| 精品少妇久久久久久888优播| 欧美成人午夜精品| 亚洲婷婷狠狠爱综合网| 女人高潮潮喷娇喘18禁视频| 成年美女黄网站色视频大全免费| xxx大片免费视频| 亚洲三区欧美一区| 少妇被粗大的猛进出69影院| 91国产中文字幕| 精品卡一卡二卡四卡免费| 在线观看一区二区三区激情| 只有这里有精品99| 丰满少妇做爰视频| 熟女av电影| av免费在线看不卡| 在线免费观看不下载黄p国产| 女人久久www免费人成看片| 2018国产大陆天天弄谢| 晚上一个人看的免费电影| 最近最新中文字幕大全免费视频 | 亚洲三级黄色毛片| 有码 亚洲区| 久久婷婷青草| 国产亚洲最大av| 日韩在线高清观看一区二区三区| 搡老乐熟女国产| 熟女av电影| 久久综合国产亚洲精品| 在线天堂中文资源库| 亚洲 欧美一区二区三区| 久久久久精品人妻al黑| 国产精品久久久久久av不卡| 亚洲人成77777在线视频| 99九九在线精品视频| 欧美最新免费一区二区三区| 咕卡用的链子| 亚洲一码二码三码区别大吗| 女性被躁到高潮视频| 人成视频在线观看免费观看| 日韩 亚洲 欧美在线| 亚洲激情五月婷婷啪啪| 婷婷色av中文字幕| 日韩不卡一区二区三区视频在线| 亚洲图色成人| 大码成人一级视频| 考比视频在线观看| 色视频在线一区二区三区| 在线天堂中文资源库| 久久韩国三级中文字幕| 肉色欧美久久久久久久蜜桃| 99久国产av精品国产电影| 免费观看av网站的网址| 国产精品99久久99久久久不卡 | 国产日韩一区二区三区精品不卡| 深夜精品福利| 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| 久久久久精品性色| 中文字幕色久视频| 亚洲欧美清纯卡通| 国产女主播在线喷水免费视频网站| www日本在线高清视频| 中文天堂在线官网| 青春草国产在线视频| 91aial.com中文字幕在线观看| 久热久热在线精品观看| 日韩制服骚丝袜av| 美女午夜性视频免费| 国产精品 欧美亚洲| 两个人免费观看高清视频| 在线观看国产h片| 91精品三级在线观看| 卡戴珊不雅视频在线播放| 大陆偷拍与自拍| 热99国产精品久久久久久7| 边亲边吃奶的免费视频| 国产精品一二三区在线看| 欧美精品人与动牲交sv欧美| 2021少妇久久久久久久久久久| 亚洲中文av在线| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| 中文欧美无线码| 丝袜脚勾引网站| 国产成人精品久久二区二区91 | 人人妻人人添人人爽欧美一区卜| 高清视频免费观看一区二区| 久久人人爽av亚洲精品天堂| 高清视频免费观看一区二区| 观看美女的网站| 在线精品无人区一区二区三| 亚洲欧洲精品一区二区精品久久久 | 久久ye,这里只有精品| 80岁老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 久久精品aⅴ一区二区三区四区 | 另类精品久久| 亚洲精品日韩在线中文字幕| 午夜激情av网站| 男的添女的下面高潮视频| 黄色视频在线播放观看不卡| 最近最新中文字幕大全免费视频 | 成年av动漫网址| 天天躁夜夜躁狠狠久久av| 在现免费观看毛片| 国产av精品麻豆| 亚洲成人手机| 女人久久www免费人成看片| 天堂8中文在线网| 日韩电影二区| av电影中文网址| 国产伦理片在线播放av一区| 久久精品久久精品一区二区三区| 人体艺术视频欧美日本| 欧美精品亚洲一区二区| 日本色播在线视频| 高清在线视频一区二区三区| 男女边吃奶边做爰视频| 不卡视频在线观看欧美| 又黄又粗又硬又大视频| 电影成人av| 韩国高清视频一区二区三区| 精品久久久久久电影网| 日韩欧美一区视频在线观看| 精品一品国产午夜福利视频| 在线观看国产h片| 五月开心婷婷网| 欧美 亚洲 国产 日韩一| 亚洲精品自拍成人| 国产爽快片一区二区三区| xxx大片免费视频| 久久亚洲国产成人精品v| 日日摸夜夜添夜夜爱| 久久人妻熟女aⅴ| 少妇熟女欧美另类| 视频区图区小说| 国产色婷婷99| 校园人妻丝袜中文字幕| 日本免费在线观看一区| 日韩中文字幕欧美一区二区 | 国产日韩欧美在线精品| 亚洲精品乱久久久久久| 欧美精品高潮呻吟av久久| 日韩大片免费观看网站| 精品亚洲成国产av| 桃花免费在线播放| 97在线人人人人妻| 桃花免费在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 少妇的丰满在线观看| 亚洲国产毛片av蜜桃av| 一区二区av电影网| 青春草视频在线免费观看| 青青草视频在线视频观看| 老司机影院毛片| 国产 精品1| 五月伊人婷婷丁香| 欧美在线黄色| xxxhd国产人妻xxx| 寂寞人妻少妇视频99o| 亚洲精品在线美女| 亚洲成国产人片在线观看| 免费看不卡的av| 毛片一级片免费看久久久久| 天天躁夜夜躁狠狠久久av| 三上悠亚av全集在线观看| 亚洲综合色网址| 伦精品一区二区三区| 一本久久精品| 久久99蜜桃精品久久| 91成人精品电影| 国产日韩一区二区三区精品不卡| 99国产综合亚洲精品| 欧美日韩av久久| 久久免费观看电影| 少妇的逼水好多| 亚洲人成网站在线观看播放| 丰满乱子伦码专区| 免费看不卡的av| 久久青草综合色| 亚洲成人av在线免费| 久久久久网色| 久久99精品国语久久久| 亚洲国产成人一精品久久久| 99久久中文字幕三级久久日本| 日韩精品免费视频一区二区三区| 多毛熟女@视频| av视频免费观看在线观看| 中国三级夫妇交换| 男女午夜视频在线观看| 久久精品久久久久久噜噜老黄| 国产精品久久久久久精品古装| 青草久久国产| 国产精品久久久久久精品古装| 女人高潮潮喷娇喘18禁视频| 一级片'在线观看视频| 成年人免费黄色播放视频| 男人爽女人下面视频在线观看| 欧美人与性动交α欧美精品济南到 | 高清视频免费观看一区二区| 亚洲av免费高清在线观看| 男女边摸边吃奶| 伊人久久国产一区二区| 成人漫画全彩无遮挡| 欧美激情高清一区二区三区 | 国产探花极品一区二区| 亚洲一区中文字幕在线| 成人国产麻豆网| 赤兔流量卡办理| 亚洲欧美一区二区三区久久| 2022亚洲国产成人精品| 波多野结衣av一区二区av| 我的亚洲天堂| 永久免费av网站大全| 熟妇人妻不卡中文字幕| 国产视频首页在线观看| 一级毛片黄色毛片免费观看视频| 亚洲精品中文字幕在线视频| 熟女少妇亚洲综合色aaa.| www.av在线官网国产| 最近2019中文字幕mv第一页| 亚洲综合色惰| 欧美人与善性xxx| 国产精品熟女久久久久浪| 五月伊人婷婷丁香| 久热这里只有精品99| 久久精品国产综合久久久| 国产精品成人在线| 十八禁高潮呻吟视频| 两个人看的免费小视频| 肉色欧美久久久久久久蜜桃| 欧美成人精品欧美一级黄| 国产精品三级大全| 亚洲成人av在线免费| 啦啦啦啦在线视频资源| 高清av免费在线| 久久女婷五月综合色啪小说| 国产亚洲一区二区精品| 国产成人精品婷婷| 色94色欧美一区二区| 黄网站色视频无遮挡免费观看| 大码成人一级视频| videosex国产| 国产乱来视频区| 亚洲av福利一区| 色94色欧美一区二区| 久久久精品区二区三区| 2018国产大陆天天弄谢| 亚洲人成77777在线视频| 久久久精品区二区三区| 亚洲精品国产一区二区精华液| 免费高清在线观看日韩| 成人二区视频| 多毛熟女@视频| 精品国产乱码久久久久久男人| 国产精品蜜桃在线观看| 国产亚洲最大av| 制服丝袜香蕉在线| 欧美成人午夜免费资源| 久久综合国产亚洲精品| 欧美精品一区二区大全| 久久精品国产a三级三级三级| 欧美日韩视频高清一区二区三区二| 成人毛片a级毛片在线播放| 精品人妻一区二区三区麻豆| 精品少妇久久久久久888优播| 91午夜精品亚洲一区二区三区|