• <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)表示
    午夜视频精品福利| 国产精品 欧美亚洲| 亚洲欧美激情综合另类| 亚洲欧美一区二区三区黑人| 国产精品久久视频播放| 日本 欧美在线| 国产高清有码在线观看视频| 六月丁香七月| 国产黄色视频一区二区在线观看 | 免费黄网站久久成人精品| 高清毛片免费看| videos熟女内射| 国产大屁股一区二区在线视频| 久久人妻av系列| 国产成人精品一,二区| 亚洲人成网站在线观看播放| 中文资源天堂在线| 高清视频免费观看一区二区 | 少妇高潮的动态图| 国模一区二区三区四区视频| 欧美一区二区国产精品久久精品| 日本猛色少妇xxxxx猛交久久| 九色成人免费人妻av| 一边亲一边摸免费视频| 亚洲精品国产成人久久av| 嫩草影院新地址| 亚洲国产最新在线播放| 国产真实乱freesex| 久久99热6这里只有精品| 欧美成人免费av一区二区三区| 亚洲天堂国产精品一区在线| 美女被艹到高潮喷水动态| 国产精品不卡视频一区二区| 啦啦啦观看免费观看视频高清| 中文资源天堂在线| 久久精品熟女亚洲av麻豆精品 | .国产精品久久| 美女脱内裤让男人舔精品视频| 欧美三级亚洲精品| 国产精品一区二区三区四区免费观看| 精品酒店卫生间| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久成人| 夫妻性生交免费视频一级片| 日本免费一区二区三区高清不卡| 亚洲国产精品sss在线观看| 99热网站在线观看| 男女边吃奶边做爰视频| 久久精品久久精品一区二区三区| 99久国产av精品国产电影| 国产麻豆成人av免费视频| 国产一级毛片在线| 日本黄色片子视频| 搡女人真爽免费视频火全软件| 五月伊人婷婷丁香| 日本与韩国留学比较| 国产精品av视频在线免费观看| 国产亚洲91精品色在线| 内射极品少妇av片p| 久久鲁丝午夜福利片| 成人毛片a级毛片在线播放| 亚洲av成人精品一二三区| 欧美一区二区国产精品久久精品| 午夜精品在线福利| 成人亚洲精品av一区二区| 亚洲精品一区蜜桃| 成年免费大片在线观看| 干丝袜人妻中文字幕| 亚洲av福利一区| 久久久久久大精品| 亚洲国产色片| 蜜桃久久精品国产亚洲av| 国产精品麻豆人妻色哟哟久久 | 成人综合一区亚洲| 国产成人一区二区在线| 欧美区成人在线视频| 亚洲成人中文字幕在线播放| 日本猛色少妇xxxxx猛交久久| 亚洲aⅴ乱码一区二区在线播放| 国产精品一二三区在线看| 一级爰片在线观看| 国产av码专区亚洲av| 男人舔女人下体高潮全视频| АⅤ资源中文在线天堂| 日本爱情动作片www.在线观看| 久久久久网色| 99久久人妻综合| 成年版毛片免费区| 久久国产乱子免费精品| 国产黄a三级三级三级人| 久久精品久久精品一区二区三区| 日本一二三区视频观看| 亚洲人与动物交配视频| 国产乱来视频区| 中文在线观看免费www的网站| 久久久午夜欧美精品| 久久鲁丝午夜福利片| 国产白丝娇喘喷水9色精品| 成年版毛片免费区| 国产黄片美女视频| 欧美最新免费一区二区三区| av天堂中文字幕网| 国产成人一区二区在线| 亚洲av二区三区四区| 国产免费一级a男人的天堂| 久久99精品国语久久久| 99久久精品热视频| 国产成人精品久久久久久| 午夜福利在线观看免费完整高清在| 日韩av不卡免费在线播放| 天天躁夜夜躁狠狠久久av| 日产精品乱码卡一卡2卡三| 国产不卡一卡二| 国产成年人精品一区二区| 欧美激情在线99| 日韩 亚洲 欧美在线| 久久久久久久久久久丰满| 亚洲不卡免费看| 免费一级毛片在线播放高清视频| 亚洲色图av天堂| 国产极品天堂在线| 国产成年人精品一区二区| 久久韩国三级中文字幕| 建设人人有责人人尽责人人享有的 | 夜夜看夜夜爽夜夜摸| 中文天堂在线官网| 亚洲欧美日韩无卡精品| 国产成人免费观看mmmm| 免费观看a级毛片全部| 日本与韩国留学比较| 亚洲18禁久久av| 亚洲aⅴ乱码一区二区在线播放| 最近的中文字幕免费完整| 国产色爽女视频免费观看| 色综合色国产| 免费观看性生交大片5| 亚洲最大成人手机在线| 最近的中文字幕免费完整| 日韩视频在线欧美| 欧美极品一区二区三区四区| 亚洲av电影不卡..在线观看| 亚洲av中文字字幕乱码综合| 夫妻性生交免费视频一级片| 国产一区二区三区av在线| 一边亲一边摸免费视频| 成人无遮挡网站| 亚洲精品456在线播放app| 综合色av麻豆| 日本欧美国产在线视频| 内地一区二区视频在线| 国内少妇人妻偷人精品xxx网站| 国产亚洲最大av| 如何舔出高潮| 欧美精品一区二区大全| 99国产精品一区二区蜜桃av| 又黄又爽又刺激的免费视频.| 国产精品人妻久久久久久| av在线蜜桃| 九九热线精品视视频播放| 一级av片app| 一级黄色大片毛片| 欧美性感艳星| 国产一区有黄有色的免费视频 | 熟妇人妻久久中文字幕3abv| 国产高清视频在线观看网站| 简卡轻食公司| 搞女人的毛片| 亚洲av成人av| 欧美极品一区二区三区四区| 国产精品综合久久久久久久免费| 免费观看性生交大片5| 老女人水多毛片| 男人舔奶头视频| 色播亚洲综合网| 久久这里只有精品中国| 床上黄色一级片| 国内精品一区二区在线观看| 别揉我奶头 嗯啊视频| 亚洲四区av| 在线天堂最新版资源| 波野结衣二区三区在线| 久久这里有精品视频免费| 天堂网av新在线| 亚洲美女搞黄在线观看| 久久99热6这里只有精品| 丰满少妇做爰视频| 一夜夜www| av视频在线观看入口| 国产亚洲5aaaaa淫片| 亚洲av日韩在线播放| 国产女主播在线喷水免费视频网站 | 美女脱内裤让男人舔精品视频| 欧美97在线视频| 亚洲三级黄色毛片| 国产又黄又爽又无遮挡在线| 午夜久久久久精精品| 久久草成人影院| 99热这里只有是精品在线观看| 韩国高清视频一区二区三区| 久久精品国产亚洲网站| 精品一区二区三区视频在线| 老司机影院成人| 日韩一区二区三区影片| a级毛色黄片| 免费看av在线观看网站| 日本欧美国产在线视频| 简卡轻食公司| 国产真实乱freesex| 少妇猛男粗大的猛烈进出视频 | 欧美3d第一页| 亚洲国产色片| 国产在视频线精品| 日韩 亚洲 欧美在线| 啦啦啦啦在线视频资源| 久久精品夜色国产| 人妻系列 视频| 精品99又大又爽又粗少妇毛片| 特大巨黑吊av在线直播| 老司机福利观看| 高清在线视频一区二区三区 | kizo精华| 国产精品乱码一区二三区的特点| 亚洲最大成人中文| 九九在线视频观看精品| 国产三级中文精品| 村上凉子中文字幕在线| av国产免费在线观看| 久久99热这里只有精品18| 日本wwww免费看| 成人国产麻豆网| 国产成人免费观看mmmm| 欧美一区二区国产精品久久精品| 国产视频首页在线观看| av在线播放精品| 99久国产av精品国产电影| 少妇人妻精品综合一区二区| 青春草国产在线视频| 欧美3d第一页| 天堂影院成人在线观看| 国产激情偷乱视频一区二区| 伊人久久精品亚洲午夜| 国语自产精品视频在线第100页| .国产精品久久| 久久久亚洲精品成人影院| 欧美日韩精品成人综合77777| 高清日韩中文字幕在线| 国产成人福利小说| 日韩强制内射视频| 国产成人a∨麻豆精品| 亚洲精品影视一区二区三区av| 老司机影院成人| 久久精品国产99精品国产亚洲性色| 国产精品蜜桃在线观看| 亚洲不卡免费看| 国产激情偷乱视频一区二区| h日本视频在线播放| 精品一区二区三区人妻视频| 午夜福利在线在线| 美女大奶头视频| 亚洲av.av天堂| 久久精品影院6| 国产高清国产精品国产三级 | 爱豆传媒免费全集在线观看| 床上黄色一级片| 91在线精品国自产拍蜜月| 亚洲va在线va天堂va国产| 国产视频内射| 精品99又大又爽又粗少妇毛片| 欧美一级a爱片免费观看看| 99久久精品一区二区三区| 九草在线视频观看| 中文天堂在线官网| 久久亚洲精品不卡| 噜噜噜噜噜久久久久久91| 搡老妇女老女人老熟妇| 插阴视频在线观看视频| 麻豆精品久久久久久蜜桃| 国产在线一区二区三区精 | 欧美成人一区二区免费高清观看| 能在线免费看毛片的网站| 欧美xxxx性猛交bbbb| 夜夜爽夜夜爽视频| 中文欧美无线码| 久久99精品国语久久久| 亚洲欧洲国产日韩| 欧美一区二区国产精品久久精品| 精品免费久久久久久久清纯| 国产亚洲av嫩草精品影院| 三级国产精品片| 亚洲三级黄色毛片| 成年版毛片免费区| 男的添女的下面高潮视频| 精品一区二区免费观看| 深爱激情五月婷婷| 亚洲综合精品二区| 晚上一个人看的免费电影| 国产 一区 欧美 日韩| 高清视频免费观看一区二区 | 可以在线观看毛片的网站| 内射极品少妇av片p| 午夜福利网站1000一区二区三区| 热99在线观看视频| 亚洲精品国产成人久久av| 一二三四中文在线观看免费高清| 2021少妇久久久久久久久久久| 啦啦啦啦在线视频资源| 看片在线看免费视频| 亚洲乱码一区二区免费版| 欧美性猛交╳xxx乱大交人| 91精品伊人久久大香线蕉| 国产黄a三级三级三级人| 岛国在线免费视频观看| 男女国产视频网站| 亚洲av免费高清在线观看| 欧美另类亚洲清纯唯美| 午夜免费男女啪啪视频观看| 精品熟女少妇av免费看| 热99在线观看视频| 嫩草影院精品99| 99久国产av精品国产电影| 蜜桃亚洲精品一区二区三区| 老师上课跳d突然被开到最大视频| 精品酒店卫生间| 国产真实伦视频高清在线观看| 国产成人福利小说| 亚洲精品国产av成人精品| 精品久久久久久成人av| 一个人看视频在线观看www免费| 中文天堂在线官网| 日韩三级伦理在线观看| 一级毛片我不卡| 亚洲人成网站在线播| 国产精品99久久久久久久久| 日韩av不卡免费在线播放| 亚洲一级一片aⅴ在线观看| 91久久精品国产一区二区成人| 51国产日韩欧美| 亚洲av中文字字幕乱码综合| 亚洲成人av在线免费| 精品久久久久久久久亚洲| 最近视频中文字幕2019在线8| 亚洲欧美中文字幕日韩二区| 日韩欧美国产在线观看| 不卡视频在线观看欧美| 美女大奶头视频| 成人午夜高清在线视频| 免费看日本二区| 亚洲欧美成人精品一区二区| 日韩精品青青久久久久久| 99久久中文字幕三级久久日本| av视频在线观看入口| 高清在线视频一区二区三区 | 国产三级中文精品| 舔av片在线| 女人被狂操c到高潮| 国产欧美日韩精品一区二区| 身体一侧抽搐| 又粗又硬又长又爽又黄的视频| 91精品伊人久久大香线蕉| 中文字幕免费在线视频6| 亚洲欧美日韩东京热| 国产真实伦视频高清在线观看| 国产黄片视频在线免费观看| 建设人人有责人人尽责人人享有的 | 国产一区二区三区av在线| 精品久久国产蜜桃| 亚洲欧美中文字幕日韩二区| 免费观看精品视频网站| 成年女人看的毛片在线观看| 97超视频在线观看视频| 男人的好看免费观看在线视频| 男女那种视频在线观看| 日本黄色视频三级网站网址| 国产男人的电影天堂91| 国产真实乱freesex| 成人欧美大片| 神马国产精品三级电影在线观看| 舔av片在线| 国产精品一二三区在线看| 伦精品一区二区三区| 久久国产乱子免费精品| 久久精品久久精品一区二区三区| 国产精品国产三级国产专区5o | 亚洲av男天堂| 午夜爱爱视频在线播放| 看黄色毛片网站| 国产精品福利在线免费观看| 国产黄色小视频在线观看| 蜜桃久久精品国产亚洲av| .国产精品久久| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品,欧美精品| 免费观看人在逋| 亚洲精品456在线播放app| 欧美区成人在线视频| 国产国拍精品亚洲av在线观看| 51国产日韩欧美| 亚洲国产精品合色在线| 精品久久久久久电影网 | 一本久久精品| 又黄又爽又刺激的免费视频.| 22中文网久久字幕| 欧美性感艳星| 亚洲不卡免费看| av在线蜜桃| 国产亚洲最大av| 嘟嘟电影网在线观看| .国产精品久久| 变态另类丝袜制服| 一级黄片播放器| 中国美白少妇内射xxxbb| 久久精品国产自在天天线| 欧美成人一区二区免费高清观看| 国产成人精品一,二区| 中文字幕av成人在线电影| 在线观看美女被高潮喷水网站| 直男gayav资源| 久久热精品热| 嫩草影院新地址| 最近最新中文字幕大全电影3| 22中文网久久字幕| 天美传媒精品一区二区| 青春草国产在线视频| 精品久久久久久久人妻蜜臀av| 午夜福利在线在线| 我要看日韩黄色一级片| 麻豆av噜噜一区二区三区| 国产探花极品一区二区| 亚洲精品,欧美精品| 国产精品精品国产色婷婷| 三级国产精品片| 不卡视频在线观看欧美| 亚洲av.av天堂| 欧美xxxx性猛交bbbb| 免费看光身美女| 日本与韩国留学比较| 国产精品久久视频播放| 欧美97在线视频| 午夜福利在线在线| 一个人免费在线观看电影| 亚洲一区高清亚洲精品| 黄色一级大片看看| 丰满人妻一区二区三区视频av| 大香蕉97超碰在线| 免费人成在线观看视频色| 91午夜精品亚洲一区二区三区| 69人妻影院| 高清毛片免费看| 亚洲中文字幕日韩| 国产成人精品婷婷| 51国产日韩欧美| 日韩成人伦理影院| 国产三级在线视频| 国产成人a∨麻豆精品| 精品人妻一区二区三区麻豆| 久久久久久九九精品二区国产| 九九热线精品视视频播放| 亚洲人与动物交配视频| 国产探花极品一区二区| 国产成人91sexporn| 少妇熟女欧美另类| 赤兔流量卡办理| 婷婷色综合大香蕉| 成人性生交大片免费视频hd| 亚洲色图av天堂| 久久久精品欧美日韩精品| 乱码一卡2卡4卡精品| 欧美日韩在线观看h| 国产淫片久久久久久久久| 成年版毛片免费区| 国产亚洲最大av| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美精品专区久久| 毛片女人毛片| 欧美激情久久久久久爽电影| 床上黄色一级片| 久久久久久久久久久丰满| 免费看a级黄色片| 亚洲最大成人手机在线| 伦理电影大哥的女人| 中国国产av一级| 国产成人a区在线观看| 成人三级黄色视频| 亚洲欧美一区二区三区国产| 你懂的网址亚洲精品在线观看 | 特大巨黑吊av在线直播| 我的老师免费观看完整版| 99国产精品一区二区蜜桃av| 欧美变态另类bdsm刘玥| 色综合亚洲欧美另类图片| 国产乱来视频区| 欧美日韩综合久久久久久| 舔av片在线| 联通29元200g的流量卡| 日本午夜av视频| 亚洲国产色片| 色网站视频免费| 视频中文字幕在线观看| 久久精品熟女亚洲av麻豆精品 | 一个人看的www免费观看视频| 国产国拍精品亚洲av在线观看| 午夜精品国产一区二区电影 | 亚州av有码| 午夜精品国产一区二区电影 | 亚洲av免费在线观看| 成人亚洲精品av一区二区| 亚洲精品亚洲一区二区| 禁无遮挡网站| 黄色一级大片看看| 欧美性猛交╳xxx乱大交人| 可以在线观看毛片的网站| 日韩制服骚丝袜av| 少妇人妻精品综合一区二区| 国产精品一区二区三区四区免费观看| 欧美性猛交黑人性爽| 亚洲在线自拍视频| 婷婷色综合大香蕉| 老司机影院成人| 国产高清不卡午夜福利| 岛国在线免费视频观看| 国产伦理片在线播放av一区| 变态另类丝袜制服| av在线老鸭窝| 亚洲一级一片aⅴ在线观看| 国产久久久一区二区三区| 91aial.com中文字幕在线观看| 99在线视频只有这里精品首页| 亚洲欧美成人综合另类久久久 | 白带黄色成豆腐渣| 国产精品三级大全| av在线老鸭窝| 日日撸夜夜添| 九九爱精品视频在线观看| 日本与韩国留学比较| 久久久国产成人免费| 亚洲中文字幕日韩| 综合色丁香网| 国产又色又爽无遮挡免| 亚洲av福利一区| 国产乱人偷精品视频| 精华霜和精华液先用哪个| 青春草国产在线视频| 国产精品一区二区性色av| 亚洲国产最新在线播放| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区三区| 日本一本二区三区精品| 国产精品无大码| 人妻系列 视频| 日本免费一区二区三区高清不卡| 最近手机中文字幕大全| 成人亚洲欧美一区二区av| 欧美zozozo另类| 内射极品少妇av片p| 欧美一区二区亚洲| 91精品一卡2卡3卡4卡| 精品久久国产蜜桃| 国产成人a区在线观看| 日韩大片免费观看网站 | 91精品一卡2卡3卡4卡| 成人特级av手机在线观看| 国产探花在线观看一区二区| 久久人人爽人人爽人人片va| 亚洲在线观看片| 国产成年人精品一区二区| 国产 一区 欧美 日韩| 菩萨蛮人人尽说江南好唐韦庄 | 国产探花极品一区二区| 舔av片在线| 亚洲人成网站在线观看播放| 色哟哟·www| 青春草视频在线免费观看| 国产成年人精品一区二区| 亚洲美女搞黄在线观看| av播播在线观看一区| 日本欧美国产在线视频| 国产精品乱码一区二三区的特点| 国内揄拍国产精品人妻在线| 最近手机中文字幕大全| 久久久a久久爽久久v久久| 中文天堂在线官网| 日韩成人av中文字幕在线观看| av在线老鸭窝| 大话2 男鬼变身卡| 女人十人毛片免费观看3o分钟| 日日摸夜夜添夜夜添av毛片| 狂野欧美激情性xxxx在线观看| 久久亚洲国产成人精品v| 色哟哟·www| 亚洲成av人片在线播放无| av在线天堂中文字幕| 不卡视频在线观看欧美| 一二三四中文在线观看免费高清| 熟妇人妻久久中文字幕3abv| av视频在线观看入口| 日韩国内少妇激情av| 美女xxoo啪啪120秒动态图| 国产成人一区二区在线| 国产高清有码在线观看视频| 久久久久国产网址| av在线播放精品| 免费无遮挡裸体视频| 18禁裸乳无遮挡免费网站照片| 日本黄色视频三级网站网址| videossex国产| 亚洲中文字幕一区二区三区有码在线看| 好男人在线观看高清免费视频| 国产精品综合久久久久久久免费| 亚洲成人久久爱视频| 亚洲精品久久久久久婷婷小说 | 精品国内亚洲2022精品成人| 久久久久久久午夜电影| 国产伦在线观看视频一区| 免费观看人在逋| 久久久久国产网址|