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

    基于正交最小二乘多面函數(shù)構(gòu)建局域海洋地磁場模型研究*

    2012-11-14 13:46:04常宜峰柴洪洲
    大地測量與地球動力學(xué) 2012年3期
    關(guān)鍵詞:抗差局域趨勢

    常宜峰 柴洪洲 王 敏

    (解放軍信息工程大學(xué)測繪學(xué)院,鄭州 450052)

    基于正交最小二乘多面函數(shù)構(gòu)建局域海洋地磁場模型研究*

    常宜峰 柴洪洲 王 敏

    (解放軍信息工程大學(xué)測繪學(xué)院,鄭州 450052)

    針對多面函數(shù)擬合局域海洋地磁場模型中質(zhì)量控制和節(jié)點(diǎn)選擇問題,提出使用抗差估計(jì)與正交最小二乘原理相結(jié)合的方法。在質(zhì)量控制過程中,將局部有顯著起伏的采樣點(diǎn)視為異常觀測值,采用抗差估計(jì)減弱異常值點(diǎn)的不利影響。在節(jié)點(diǎn)選擇中,使用基于正交最小二乘原則判斷節(jié)點(diǎn)的貢獻(xiàn)效力大小,按照此原則自適應(yīng)地選擇節(jié)點(diǎn)。結(jié)果表明,基于正交最小二乘原理的多面函數(shù)法的精度和穩(wěn)定性更好,更適于構(gòu)建局域海洋地磁場模型。

    局域海洋地磁場模型;抗差估計(jì);正交最小二乘原理;多面函數(shù);節(jié)點(diǎn)

    1 引言

    地磁場在了解地球內(nèi)部的物質(zhì)分布和輔助導(dǎo)航中有重要用途,其中,地磁場模型起著非常重要的作用。常用的是泰勒多項(xiàng)式模型,優(yōu)點(diǎn)是計(jì)算簡單,但地磁數(shù)據(jù)存在較大跳變性,難以取得較高的精度。1971年Hardy給出了多面函數(shù)擬合法,1976年用于重力異常和大地水準(zhǔn)面差距計(jì)算,1978年用于地殼形變計(jì)算,2009年趙建虎將此方法引入局域海洋地磁場建模中,并提出分區(qū)建模來提高精度[1]。多面函數(shù)是利用一些規(guī)則的數(shù)學(xué)曲面來逼近實(shí)際曲面,從而獲得較高的精度。本文提出將抗差估計(jì)與正交最小二乘相結(jié)合來擬合局域海洋地磁場模型。

    對于多面函數(shù)模型,核函數(shù)的選擇和節(jié)點(diǎn)的確定至關(guān)重要。對于核函數(shù),已有學(xué)者進(jìn)行過深入研究。正雙曲面函數(shù)構(gòu)建的模型在檢查點(diǎn)上的誤差分布和精度統(tǒng)計(jì)參數(shù)均優(yōu)于倒雙曲面函數(shù),并且正雙曲面核函數(shù)構(gòu)造的模型精度隨平滑因子的變化平緩而相對穩(wěn)定,倒雙曲面函數(shù)構(gòu)造的模型精度變化劇烈且最優(yōu)平滑因子相對較大,而過大的平滑因子可能造成局部模型失真,尤其對于地磁變化相對較大的數(shù)據(jù)的擬合,在此取正雙曲面函數(shù),平滑因子取0.1[1]。對于節(jié)點(diǎn)的選擇,通常是直接取采樣點(diǎn),可是選擇節(jié)點(diǎn)不同產(chǎn)生的結(jié)果也不同,具有一定的隨意性。由于每個(gè)節(jié)點(diǎn)對擬合結(jié)果貢獻(xiàn)是不同的,有學(xué)者采用t檢驗(yàn)或逐步回歸的方法來選擇節(jié)點(diǎn),得到的效果較好[2]。本文利用正交最小二乘原理的基本思想,首先將誤差方程系數(shù)正交化,并計(jì)算觀測向量與正交化向量之間的夾角,節(jié)點(diǎn)對曲面擬合的貢獻(xiàn)就用夾角的大小來表示。按照這種原則方法依次選擇合適的節(jié)點(diǎn),完成自主選擇節(jié)點(diǎn)的目的[3]。此方法的優(yōu)點(diǎn)是既保證了擬合的精度又兼顧了擬合的穩(wěn)定性。

    2 基于抗差估計(jì)的質(zhì)量控制

    在測量過程中受到各種因素的影響,可能出現(xiàn)異常觀測數(shù)據(jù),一般的質(zhì)量控制是采用判斷測量值與對應(yīng)趨勢面擬合值之間差值大小的趨勢面法來進(jìn)行的。趨勢面分析方法是將短尺度的、局部的變化看作隨機(jī)的和非結(jié)構(gòu)的噪聲,因此當(dāng)局部區(qū)域變化比較劇烈時(shí),可將其視作異常觀測。但是單純的趨勢面擬合并不具備抗差的能力,故本文引入抗差估計(jì)理論[4,5],采用抗差趨勢面法進(jìn)行質(zhì)量控制。

    趨勢面函數(shù)如下式所示

    其中,αij(i=0,…,n,j=0,…,n)為待求系數(shù),(x,y)為觀測點(diǎn)坐標(biāo)為對應(yīng)的磁力觀測值。如果有m個(gè)觀測點(diǎn),應(yīng)用最小二乘準(zhǔn)則即可解得多項(xiàng)式系數(shù)。

    設(shè)有誤差方程

    根據(jù)抗差估計(jì)理論,以及等價(jià)權(quán)理論可得

    3 基于正交最小二乘法節(jié)點(diǎn)選擇

    經(jīng)過質(zhì)量控制后,再對剩下的殘差進(jìn)行擬合。根據(jù)多面函數(shù)原理有,

    在已知有限測量數(shù)據(jù)的情況下求解方程,重要的是節(jié)點(diǎn)的選擇,一般情況下直接取均勻分布在研究區(qū)域內(nèi)的觀測點(diǎn),如果選擇所有的采樣點(diǎn)作為節(jié)點(diǎn),可能會出現(xiàn)內(nèi)符合精度虛高,推估和內(nèi)插精度無法保證的情況。特別在采樣點(diǎn)受異常誤差污染時(shí),擬合精度會更低。為解決以上問題,本文在此引入基于正交最小二乘原理的思想來建立一種自主選擇節(jié)點(diǎn)的方法。

    在m個(gè)磁力數(shù)據(jù)中選擇n個(gè)數(shù)據(jù)建模,即

    其中,X為模型系數(shù),A為核函數(shù)矩陣,由最小二乘原理得,

    模型精度評價(jià)指標(biāo)為[6]

    其中,n為節(jié)點(diǎn)個(gè)數(shù),(xi,yi)為節(jié)點(diǎn)坐標(biāo),βi為模型系數(shù),Q(x,y,xi,yi)為二次核函數(shù),一般取正雙曲函數(shù)(如下式),δ2為模型的平滑因子(本文取0.1)。

    先假設(shè)選取所有的采樣點(diǎn)為節(jié)點(diǎn),即m=n,此時(shí)構(gòu)建的誤差方程系數(shù)陣為

    其中,

    從式(11)可以看出,每個(gè)ai對求系數(shù)β的決定效力是有差別的,即每一節(jié)點(diǎn)對函數(shù)擬合的影響是不同的,故提供了一種評價(jià)準(zhǔn)則,便于依照節(jié)點(diǎn)貢獻(xiàn)效力的大小來選擇節(jié)點(diǎn)。

    下面就引入一種表示節(jié)點(diǎn)貢獻(xiàn)大小的原則,將式(5)表示成向量形式,

    假設(shè)ak之間是相互正交的,則左乘βk有,

    兩邊求和并移項(xiàng)得

    定義η為選擇t(t≤m)個(gè)節(jié)點(diǎn)的列向量時(shí)的總貢獻(xiàn)

    η的值在0到1之間。并且當(dāng)t越大,η就越大,表明逼近程度越高;當(dāng)t=m時(shí),η=1,此時(shí)選擇了所有的采樣點(diǎn),擬合誤差為零,逼近精度最高。但由于沒有多余觀測,可靠性較差,為此應(yīng)根據(jù)每個(gè)基矢量的貢獻(xiàn)大小來選擇若干節(jié)點(diǎn),使η足夠大且有多余觀測。

    要想使η最大,即需使βk取最大值,在式(12)兩邊同時(shí)左乘,可得

    顯然,要使βk取得最大值,需向量ak與觀測向量Z之間的夾角θakZ最大,即

    其具體計(jì)算步驟為:

    1)將所有的采樣點(diǎn)視為剩余節(jié)點(diǎn),計(jì)算系數(shù)陣A,m個(gè)列向量可記為a1,a2,…,am;

    2)計(jì)算Z與ak之間的夾角,求得最大值

    式(18)表明aj對Z的貢獻(xiàn)最大,于是可將與aj相對應(yīng)的坐標(biāo)(xj,yj)選為節(jié)點(diǎn)。

    3)按照最小二乘原理解算系數(shù)β

    4 實(shí)驗(yàn)分析

    為驗(yàn)證以上所述方法在局域海洋地磁場模型建立中的有效性,采用某測區(qū)海洋磁力測量的數(shù)據(jù)進(jìn)行計(jì)算和分析。該測區(qū)定位系統(tǒng)采用NGD-60型差分GPS定位系統(tǒng),船載磁力儀采用美國Geometrics公司的G-881型海洋磁力儀,共有1 681個(gè)測量數(shù)據(jù)。經(jīng)過了各項(xiàng)的改正,測區(qū)的等值線如圖1所示。均勻選取其中的441個(gè)數(shù)據(jù)作為已知的測量點(diǎn),將剩余的測量數(shù)據(jù)作為外部的檢核點(diǎn),分別使用多面函數(shù)法、正交最小二乘多面函數(shù)法、基于抗差趨勢面的多面函數(shù)法和基于抗差趨勢面的正交最小二乘多面函數(shù)法,4種方法進(jìn)行計(jì)算(表1,其中最大值和最小值是指模型計(jì)算結(jié)果與實(shí)測值的偏離的最大和最小值,平均值是所有偏離量的算數(shù)平均值,均方誤差是對外部檢核點(diǎn)的綜合評價(jià)指標(biāo))。

    圖1 某測區(qū)磁力異常分布Fig.1 Distribution of magnetic abnormity in the studied area

    表1 不同擬合方法的結(jié)果比較(單位:nT)Tab.1 Comparison among results with different methods (unit:nT)

    從表1可以看出,直接使用多面函數(shù)和正交多面函數(shù)進(jìn)行擬合,均方誤差分別可達(dá)17nT和15nT,平均值都在0.6 nT左右,說明模型的集中程度較好,兩者所得到的精度基本相當(dāng),但正交多面函數(shù)法略好于多面函數(shù)法,并且都有所偏大。加上抗差估計(jì)趨勢面進(jìn)行質(zhì)量控制之后,同樣比較兩種方法的計(jì)算結(jié)果,均方誤差分別為12 nT和7 nT,精度提高了41%,平均值分別為-8.4 nT和0.2 nT,基于抗差趨勢面的多面函數(shù)出現(xiàn)了明顯的系統(tǒng)偏差,分析原因是多面函數(shù)法對節(jié)點(diǎn)沒有選擇性,不是選擇一些特征節(jié)點(diǎn)進(jìn)行建模,而是將所有的節(jié)點(diǎn)作為建模節(jié)點(diǎn),從而引入了一些系統(tǒng)誤差,進(jìn)而產(chǎn)生了系統(tǒng)性偏差。比較多面函數(shù)和基于抗差趨勢面的多面函數(shù),增加了抗差以后,精度由17 nT變?yōu)?2 nT;同樣,對正交多面函數(shù)增加了抗差以后,精度由15 nT變?yōu)? nT,說明基于抗差趨勢面的質(zhì)量控制可以減小一部分不良影響。綜合表明基于抗差趨勢面的正交多面函數(shù)法總體上要優(yōu)于以上其他3種方法。

    圖2和圖3分別是使用基于抗差趨勢面的正交多面函數(shù)法和基于抗差趨勢面的多面函數(shù)法擬合的檢驗(yàn)點(diǎn)的殘差分布圖,從圖中可以看出,圖2的殘差分布相對更為集中,且普遍都靠近零點(diǎn)。

    為進(jìn)一步驗(yàn)證按照正交多面函數(shù)法選擇節(jié)點(diǎn)的多少與所得模型擬合精度之間的關(guān)系,依次計(jì)算了選擇個(gè)檢驗(yàn)點(diǎn)時(shí)的均方誤差,并將用正交多面函數(shù)選擇的節(jié)點(diǎn)個(gè)數(shù)與擬合的均方誤差的關(guān)系繪制成圖4。

    圖2 基于抗差趨勢面的正交多面函數(shù)法擬合殘差Fig.2 Fitting residual with orthogonal multi-function method based on robust trend-surface

    從圖4可以看出,隨著選擇的節(jié)點(diǎn)個(gè)數(shù)的增加,均方誤差逐漸減少。在采用250個(gè)節(jié)點(diǎn)時(shí),精度已經(jīng)達(dá)到了17 nT,基本上相當(dāng)于選擇所有節(jié)點(diǎn)進(jìn)行多面函數(shù)法建模時(shí)的精度。隨著選擇節(jié)點(diǎn)個(gè)數(shù)的繼續(xù)增加,精度基本保持穩(wěn)定并有所提高。同時(shí)應(yīng)該指出,針對地磁變化比較復(fù)雜和比較平緩的兩種區(qū)域,數(shù)據(jù)的代表性對于模型的精度有著重要的影響。因此,在變化比較復(fù)雜的區(qū)域可以適當(dāng)?shù)姆艑挼V沟呐袛鄺l件,盡可能多的自主選擇特征點(diǎn)進(jìn)行建模。

    圖3 基于抗差趨勢面的多面函數(shù)法擬合殘差Fig.3 Fitting residual with multi-function method based on robust trend-surface

    圖4 節(jié)點(diǎn)選擇與擬合誤差關(guān)系Fig.4 Relation between node collecting and error fitting

    圖5 兩種方法擬合效果比較Fig.5 Comparion between fitting results with two methods

    為了使計(jì)算結(jié)果更直觀,用正交多面函數(shù)法和一般多面函數(shù)法分別計(jì)算了整個(gè)區(qū)域的網(wǎng)格點(diǎn)地磁場總強(qiáng)度值,并利用計(jì)算所得數(shù)據(jù)分別繪制正交多面函數(shù)法和普通多面函數(shù)法磁力異常圖,與該區(qū)域?qū)嶋H磁力異常場圖進(jìn)行比較。圖5(a)是原始的磁力異常分布圖,(b)是使用基于抗差趨勢面的正交多面函數(shù)法得到的擬合圖,(c)是使用基于抗差趨勢面的多面函數(shù)法得到的結(jié)果。從圖中可以看出,圖5(b)與(a)等高線的分布較為接近,顯現(xiàn)出了較為細(xì)節(jié)的分布信息,保持了較高的一致性;(c)與(a)的一致性較差。從而顯示圖5(b)的建模方法是穩(wěn)定可靠的。

    5 結(jié)論

    1)進(jìn)一步驗(yàn)證了多面函數(shù)法比較適用于局域海洋地磁場建模,在地磁變化復(fù)雜的海域,先進(jìn)行趨勢面擬合,再對剩余的誤差進(jìn)行多面函數(shù)擬合,可以突出局部變化,擬合效果要優(yōu)于直接使用多面函數(shù)進(jìn)行擬合。

    2)趨勢面擬合主要反映整體變化趨勢,當(dāng)局部磁力異常觀測值變化顯著時(shí),采用抗差趨勢面法可以減小對趨勢面擬合的不良影響。在多面函數(shù)擬合前,采用抗差估計(jì)理論進(jìn)行質(zhì)量控制,可以更有效剔除異常觀測值對趨勢面擬合的影響。

    3)使用正交最小二乘多面函數(shù)法得到的擬合精度優(yōu)于多面函數(shù)法,可以按照規(guī)則自適應(yīng)地選擇特征點(diǎn),提高了人工選擇特征點(diǎn)的效率。對系數(shù)矩陣進(jìn)行正交化,不會影響最終的結(jié)果,因?yàn)檎换蟛]有參與多面函數(shù)系數(shù)的求解,僅是用于節(jié)點(diǎn)的選擇,因此可以作為一種自動選擇節(jié)點(diǎn)的參考方法。

    1 趙建虎,等.基于多面函數(shù)的局域地磁場建模方法研究[J].海洋通報(bào),2009,28(4):89-97.(Zhao Jianhu,et al.Study on construction of local marine geomagnetic field model based on multi-surface function[J].Marine Science Bulletin,2009,28(4):89-97)

    2 劉兆平,楊進(jìn),武煒.地球物理數(shù)據(jù)網(wǎng)格化方法的選?。跩].物探與化探,2010,34(1):93-98.(Liu Zhaoping,Yang Jin and Wu Wei.The choice of gridding methods for geophysical data[J].Geophysicalamp;Eochemical Exploration,2010,34(1):93-98)

    3 張菊清,劉平芝.抗差趨勢面與正交多面函數(shù)結(jié)合擬合DEM數(shù)據(jù)[J].測繪學(xué)報(bào),2008,37(4):526-531.(Zhang Juqing and Liu Pingzhi.Combining fitting based on robust trend surface and orthogonal multiquadrics with application in DEM fitting[J].Acta Geodaetica et Cartographica Sinica,2008,37(4):526-531)

    4 楊元喜.抗差估計(jì)理論及其應(yīng)用[M].北京:八一出版社,1993.(Yang Yuanxi.The theory and application of robust estimation[M].Beijing:Bayi Publishing House,1993)

    5 楊元喜.自適應(yīng)抗差最小二乘估計(jì)[J].測繪學(xué)報(bào),1996,25(3):206-211.(Yang Yuanxi.A daptively robust least squares estimation[J].Acta Geodaetica et Cartographica Sinica,1996,25(3):206-211)

    6 黃學(xué)功,等.地磁圖制備方法及其有效性評估[J].北京航空航天大學(xué)學(xué)報(bào),2009,35(7):891-895.(Huang Xuegong,et al.Geomagnetic mapping and validity estimation[J].Journal of Beijing University of Aeronautics And Stronautics,2009,35(7):891-895)

    STUDY ON CONSTRUCTION OF LOCAL MARINE GEOMAGNETIC FIELD MODEL BASED ON ORTHOGONAL LEAST SQUARE MULTI-SURFACE FUNCTION

    Chang Yifeng,Chai Hongzhou and Wang Min
    (Institute of Surveying and Mapping,Information Engineering University,Zhengzhou 450052)

    Aim at the problem of quality control and node choice in modeling local marine geomagnetic field based on multi-surface function,a fitting method with combining robust trend surface and orthogonal least squares is proposed.In order to control the influence of some outstanding points to surface fitting,a robust fitting of the trend surface by using an equivalent weight is adopted and the adaptive node choosing method is proposed based on the effect of every node on the curve fitting calculated by using the orthogonal least squares as well.The results show that this method is more effective in accuracy and stability than normal method.

    local marine geomagnetic field;robust estimate;orthogonal least square;multi-surface function;node

    1671-5942(2012)03-0086-05

    2011-11-09

    常宜峰,男,1986年生,碩士研究生,主要從事測量數(shù)據(jù)處理理論與方法研究.E-mail:changyifeng1986@163.com

    P207

    A

    猜你喜歡
    抗差局域趨勢
    趨勢
    局域積分散列最近鄰查找算法
    電子測試(2018年18期)2018-11-14 02:30:34
    初秋唇妝趨勢
    Coco薇(2017年9期)2017-09-07 21:23:49
    改善單頻PPP參數(shù)收斂速度的抗差估計(jì)方法
    SPINEXPO?2017春夏流行趨勢
    PET成像的高分辨率快速局域重建算法的建立
    地形簡化對DEM不確定性的抗差性研究
    基于抗差最小均方估計(jì)的輸電線路參數(shù)辨識
    基于局域波法和LSSVM的短期負(fù)荷預(yù)測
    電測與儀表(2015年7期)2015-04-09 11:39:50
    基于非正交變換的局域波束空時(shí)自適應(yīng)處理
    女人十人毛片免费观看3o分钟| 三级毛片av免费| 菩萨蛮人人尽说江南好唐韦庄 | 精品不卡国产一区二区三区| 亚洲国产精品国产精品| 国产成人影院久久av| 国产亚洲91精品色在线| 国产成人91sexporn| 人人妻人人澡人人爽人人夜夜 | 免费搜索国产男女视频| 卡戴珊不雅视频在线播放| 色吧在线观看| 永久网站在线| 伊人久久精品亚洲午夜| 亚洲成人av在线免费| 波野结衣二区三区在线| 免费观看精品视频网站| 久久人妻av系列| 亚洲国产精品久久男人天堂| 熟女人妻精品中文字幕| av免费在线看不卡| 久久欧美精品欧美久久欧美| 亚洲成人久久性| 1000部很黄的大片| 搡老岳熟女国产| av在线观看视频网站免费| 日韩精品青青久久久久久| 18禁黄网站禁片免费观看直播| 久久精品久久久久久噜噜老黄 | 国产精品国产三级国产av玫瑰| 欧美又色又爽又黄视频| 免费看光身美女| 大又大粗又爽又黄少妇毛片口| 欧美色欧美亚洲另类二区| 51国产日韩欧美| 简卡轻食公司| 亚洲欧美清纯卡通| 长腿黑丝高跟| 婷婷精品国产亚洲av在线| 婷婷精品国产亚洲av| 亚洲在线观看片| 一级黄片播放器| 亚洲精品影视一区二区三区av| 久久人人爽人人爽人人片va| 精品乱码久久久久久99久播| 在线观看av片永久免费下载| 亚洲国产日韩欧美精品在线观看| 免费黄网站久久成人精品| av在线老鸭窝| 精品无人区乱码1区二区| 舔av片在线| 亚洲国产精品成人综合色| 亚洲自偷自拍三级| 久久国产乱子免费精品| 国产成人91sexporn| 91午夜精品亚洲一区二区三区| 1024手机看黄色片| a级毛片a级免费在线| 日韩中字成人| 精品人妻偷拍中文字幕| 永久网站在线| 看十八女毛片水多多多| av在线老鸭窝| 久久精品国产99精品国产亚洲性色| 午夜免费激情av| 成人漫画全彩无遮挡| 一进一出好大好爽视频| 大型黄色视频在线免费观看| 婷婷亚洲欧美| 国产探花在线观看一区二区| 日韩一区二区视频免费看| 成年女人毛片免费观看观看9| 久久久国产成人免费| 女人十人毛片免费观看3o分钟| 麻豆国产97在线/欧美| 久久婷婷人人爽人人干人人爱| 亚洲激情五月婷婷啪啪| 国语自产精品视频在线第100页| 97超碰精品成人国产| 国产色爽女视频免费观看| 成人二区视频| 在现免费观看毛片| 观看美女的网站| 天堂√8在线中文| 国产黄a三级三级三级人| 色噜噜av男人的天堂激情| 精品人妻偷拍中文字幕| 在线观看66精品国产| 国产精品不卡视频一区二区| 亚洲国产精品久久男人天堂| 国产日本99.免费观看| 高清毛片免费观看视频网站| 我要搜黄色片| 久久人人爽人人爽人人片va| 亚洲av免费在线观看| 国产精品久久久久久精品电影| 少妇被粗大猛烈的视频| 国产三级在线视频| 在线观看av片永久免费下载| 午夜福利高清视频| 国产真实伦视频高清在线观看| 在线观看av片永久免费下载| 97碰自拍视频| 精华霜和精华液先用哪个| 97超视频在线观看视频| 18禁在线播放成人免费| 亚洲va在线va天堂va国产| 国产视频内射| 国产精品永久免费网站| 久久久a久久爽久久v久久| 亚洲无线观看免费| 女同久久另类99精品国产91| 欧美激情国产日韩精品一区| 色综合色国产| av在线蜜桃| 国产一区二区三区av在线 | 中文亚洲av片在线观看爽| 99国产极品粉嫩在线观看| 亚洲av.av天堂| 噜噜噜噜噜久久久久久91| 亚洲五月天丁香| 国产中年淑女户外野战色| 禁无遮挡网站| 国产美女午夜福利| 欧美极品一区二区三区四区| 亚洲av成人精品一区久久| 校园春色视频在线观看| 看片在线看免费视频| 精品99又大又爽又粗少妇毛片| 国产大屁股一区二区在线视频| 国产在线男女| 欧美日韩一区二区视频在线观看视频在线 | 国产黄片美女视频| 国产精品亚洲美女久久久| 身体一侧抽搐| 日日撸夜夜添| 欧美zozozo另类| 乱系列少妇在线播放| 午夜久久久久精精品| 欧美激情久久久久久爽电影| 亚洲精品亚洲一区二区| www.色视频.com| 超碰av人人做人人爽久久| 亚洲精品影视一区二区三区av| 国产精品嫩草影院av在线观看| 99热6这里只有精品| 久久久久久久久久久丰满| 国产亚洲精品综合一区在线观看| 日本黄色视频三级网站网址| 一区福利在线观看| 国产精品一区二区三区四区久久| 欧美成人精品欧美一级黄| av在线亚洲专区| 22中文网久久字幕| 色吧在线观看| 亚洲四区av| 亚洲美女搞黄在线观看 | 99热网站在线观看| 天天躁日日操中文字幕| 亚洲美女搞黄在线观看 | 午夜免费激情av| 十八禁国产超污无遮挡网站| 亚洲人成网站在线播放欧美日韩| 亚洲av中文av极速乱| 99久久精品热视频| 色播亚洲综合网| 国产91av在线免费观看| 免费av不卡在线播放| eeuss影院久久| 变态另类丝袜制服| 国产亚洲精品久久久com| 久久久久精品国产欧美久久久| 精品不卡国产一区二区三区| 亚洲av免费在线观看| 最后的刺客免费高清国语| 免费看日本二区| 亚洲欧美成人精品一区二区| 欧美日韩一区二区视频在线观看视频在线 | 国产欧美日韩精品亚洲av| 18禁裸乳无遮挡免费网站照片| 老司机午夜福利在线观看视频| 国产亚洲精品综合一区在线观看| 22中文网久久字幕| 国产男人的电影天堂91| 欧美xxxx黑人xx丫x性爽| 欧美极品一区二区三区四区| 国产精品亚洲美女久久久| 亚洲真实伦在线观看| 舔av片在线| 久久精品国产亚洲av天美| 国产女主播在线喷水免费视频网站 | 中文在线观看免费www的网站| 国产高清视频在线播放一区| 亚洲aⅴ乱码一区二区在线播放| 丝袜喷水一区| 久久国产乱子免费精品| 人妻少妇偷人精品九色| 婷婷精品国产亚洲av| 有码 亚洲区| 精品日产1卡2卡| 亚洲专区国产一区二区| av在线天堂中文字幕| 国内久久婷婷六月综合欲色啪| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| 亚洲人成网站高清观看| 欧美日韩精品成人综合77777| 欧美激情在线99| 可以在线观看的亚洲视频| 九色成人免费人妻av| 18禁在线播放成人免费| 最近手机中文字幕大全| 中文在线观看免费www的网站| 免费无遮挡裸体视频| 久久精品久久久久久噜噜老黄 | 中国美白少妇内射xxxbb| 欧美又色又爽又黄视频| 露出奶头的视频| 少妇人妻精品综合一区二区 | 美女高潮的动态| 尾随美女入室| 亚洲自拍偷在线| 亚洲欧美日韩卡通动漫| 亚洲人成网站在线播| 国产精品乱码一区二三区的特点| 在线观看免费视频日本深夜| 午夜精品一区二区三区免费看| 国产亚洲精品综合一区在线观看| 12—13女人毛片做爰片一| 国产精品av视频在线免费观看| 激情 狠狠 欧美| 午夜激情欧美在线| 美女被艹到高潮喷水动态| 免费在线观看影片大全网站| 日本精品一区二区三区蜜桃| 成人性生交大片免费视频hd| 国产激情偷乱视频一区二区| 99久久久亚洲精品蜜臀av| 天堂动漫精品| 麻豆国产av国片精品| 一级毛片aaaaaa免费看小| 大型黄色视频在线免费观看| 最近视频中文字幕2019在线8| 免费在线观看成人毛片| 日韩欧美在线乱码| 午夜影院日韩av| 69人妻影院| 午夜激情福利司机影院| 精品一区二区三区人妻视频| 性插视频无遮挡在线免费观看| 久久久欧美国产精品| 国产私拍福利视频在线观看| 精品午夜福利视频在线观看一区| 欧美+亚洲+日韩+国产| 日韩欧美精品免费久久| 日韩精品中文字幕看吧| 国产精品女同一区二区软件| 精品久久久久久久人妻蜜臀av| 国产成人91sexporn| 日本黄大片高清| АⅤ资源中文在线天堂| 久久人妻av系列| 亚洲成人中文字幕在线播放| 国产精品一二三区在线看| 亚洲第一区二区三区不卡| 天美传媒精品一区二区| 精品久久久久久久人妻蜜臀av| 久久久久久久久大av| 色视频www国产| 亚洲欧美精品自产自拍| 18禁在线无遮挡免费观看视频 | 亚洲国产欧美人成| 亚洲人成网站高清观看| 色综合亚洲欧美另类图片| 国产精品一二三区在线看| av在线蜜桃| 性插视频无遮挡在线免费观看| 国产精品不卡视频一区二区| 亚洲av成人av| 亚洲天堂国产精品一区在线| 男女做爰动态图高潮gif福利片| 秋霞在线观看毛片| 国产精品综合久久久久久久免费| 亚洲四区av| 三级经典国产精品| 一个人免费在线观看电影| 精品国内亚洲2022精品成人| 日韩欧美三级三区| 亚洲国产精品国产精品| 一区福利在线观看| 成熟少妇高潮喷水视频| 美女大奶头视频| 最近的中文字幕免费完整| 国产熟女欧美一区二区| 在线观看午夜福利视频| 人人妻人人澡人人爽人人夜夜 | 一级毛片aaaaaa免费看小| 国产成人91sexporn| 又爽又黄无遮挡网站| 精品福利观看| 欧美日韩精品成人综合77777| 日产精品乱码卡一卡2卡三| 毛片一级片免费看久久久久| 国产欧美日韩精品一区二区| eeuss影院久久| 最近在线观看免费完整版| 成人欧美大片| 91av网一区二区| 午夜福利视频1000在线观看| 国产伦精品一区二区三区四那| 一进一出抽搐gif免费好疼| 最近视频中文字幕2019在线8| 久久精品国产99精品国产亚洲性色| 黄色日韩在线| АⅤ资源中文在线天堂| 99热这里只有是精品在线观看| 99久久九九国产精品国产免费| 精品乱码久久久久久99久播| 可以在线观看的亚洲视频| 伊人久久精品亚洲午夜| 高清日韩中文字幕在线| 最近2019中文字幕mv第一页| 不卡视频在线观看欧美| 99久久中文字幕三级久久日本| 国产麻豆成人av免费视频| 伊人久久精品亚洲午夜| 日韩三级伦理在线观看| 人妻久久中文字幕网| 日本与韩国留学比较| 国产亚洲av嫩草精品影院| 赤兔流量卡办理| 男女下面进入的视频免费午夜| 永久网站在线| 美女免费视频网站| 国模一区二区三区四区视频| 天堂网av新在线| 亚洲成人久久性| 观看美女的网站| 最近最新中文字幕大全电影3| 日韩人妻高清精品专区| 女人十人毛片免费观看3o分钟| 国产大屁股一区二区在线视频| 久久精品人妻少妇| 免费人成视频x8x8入口观看| 日韩大尺度精品在线看网址| 亚洲国产精品sss在线观看| 久久九九热精品免费| 亚洲无线观看免费| 美女黄网站色视频| 少妇熟女欧美另类| 免费av不卡在线播放| 91久久精品国产一区二区三区| 免费观看精品视频网站| 亚洲精品粉嫩美女一区| 国产精品电影一区二区三区| 亚洲精品粉嫩美女一区| 少妇的逼水好多| 有码 亚洲区| 啦啦啦观看免费观看视频高清| 国产一区亚洲一区在线观看| 国产探花在线观看一区二区| 国产一区二区激情短视频| 国内精品久久久久精免费| 亚洲自偷自拍三级| 午夜激情福利司机影院| 国内精品美女久久久久久| 最近手机中文字幕大全| 寂寞人妻少妇视频99o| 性色avwww在线观看| av福利片在线观看| 日本一二三区视频观看| 精品一区二区三区视频在线观看免费| 一本一本综合久久| 两个人的视频大全免费| 午夜爱爱视频在线播放| 此物有八面人人有两片| 国产日本99.免费观看| 婷婷精品国产亚洲av在线| 国产午夜福利久久久久久| 一本久久中文字幕| 色哟哟哟哟哟哟| 亚洲最大成人手机在线| 色哟哟哟哟哟哟| 午夜精品在线福利| 国产精品综合久久久久久久免费| 非洲黑人性xxxx精品又粗又长| 神马国产精品三级电影在线观看| 欧美丝袜亚洲另类| 中文字幕av成人在线电影| 国产麻豆成人av免费视频| 91午夜精品亚洲一区二区三区| 久久精品夜色国产| 亚洲真实伦在线观看| 久久99热6这里只有精品| 色在线成人网| 如何舔出高潮| 久久99热这里只有精品18| 国产av一区在线观看免费| 久久久久久伊人网av| 国产aⅴ精品一区二区三区波| 国产精品国产三级国产av玫瑰| 日本一本二区三区精品| 国产av不卡久久| 免费一级毛片在线播放高清视频| 91久久精品国产一区二区三区| 97超碰精品成人国产| 99九九线精品视频在线观看视频| 亚洲最大成人中文| videossex国产| 你懂的网址亚洲精品在线观看 | 在线免费观看的www视频| 老女人水多毛片| 国产精品国产三级国产av玫瑰| 国产成人a区在线观看| 亚洲国产欧洲综合997久久,| 99久久精品一区二区三区| 国产成人福利小说| 午夜日韩欧美国产| 亚洲七黄色美女视频| 国产白丝娇喘喷水9色精品| 成人无遮挡网站| 中文亚洲av片在线观看爽| 久久久色成人| 97碰自拍视频| 极品教师在线视频| 久久久久久久久大av| 亚洲人成网站高清观看| 丰满乱子伦码专区| 日韩国内少妇激情av| 色播亚洲综合网| 成人亚洲欧美一区二区av| 给我免费播放毛片高清在线观看| 中文字幕av成人在线电影| 99热这里只有是精品在线观看| 91av网一区二区| 亚洲美女搞黄在线观看 | 人人妻人人看人人澡| 国产精品女同一区二区软件| 网址你懂的国产日韩在线| 日韩av在线大香蕉| 国产白丝娇喘喷水9色精品| 久久韩国三级中文字幕| 欧美区成人在线视频| 久久久久国内视频| 国产精品国产高清国产av| 超碰av人人做人人爽久久| 偷拍熟女少妇极品色| 亚洲国产精品国产精品| 一a级毛片在线观看| 麻豆av噜噜一区二区三区| 成年女人看的毛片在线观看| 少妇猛男粗大的猛烈进出视频 | 午夜免费激情av| 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看 | 69av精品久久久久久| 中文字幕av在线有码专区| 日本精品一区二区三区蜜桃| 成熟少妇高潮喷水视频| 亚洲中文字幕一区二区三区有码在线看| 丝袜美腿在线中文| 国产成人一区二区在线| 色哟哟·www| 最近最新中文字幕大全电影3| 欧美bdsm另类| 免费在线观看成人毛片| 久久精品91蜜桃| 99久久无色码亚洲精品果冻| 亚洲性久久影院| 综合色av麻豆| 不卡视频在线观看欧美| 日韩欧美精品免费久久| 一进一出抽搐gif免费好疼| 桃色一区二区三区在线观看| 久久这里只有精品中国| 黄色欧美视频在线观看| 噜噜噜噜噜久久久久久91| 少妇熟女aⅴ在线视频| 国产精品久久视频播放| 国产高潮美女av| 啦啦啦啦在线视频资源| 久久99热这里只有精品18| 啦啦啦观看免费观看视频高清| 中国国产av一级| 国产白丝娇喘喷水9色精品| 久久九九热精品免费| 欧美又色又爽又黄视频| 中文字幕免费在线视频6| 国产精品免费一区二区三区在线| 小蜜桃在线观看免费完整版高清| 丝袜美腿在线中文| 亚洲中文字幕一区二区三区有码在线看| 在线免费观看不下载黄p国产| 日韩中字成人| 人妻丰满熟妇av一区二区三区| 亚洲av免费在线观看| 精品人妻一区二区三区麻豆 | 日日摸夜夜添夜夜添小说| 久久99热这里只有精品18| 性欧美人与动物交配| 久久久午夜欧美精品| 热99re8久久精品国产| 久久久久久九九精品二区国产| 毛片女人毛片| 熟女电影av网| 性色avwww在线观看| 一级黄色大片毛片| 国内精品宾馆在线| 麻豆久久精品国产亚洲av| 国产成人freesex在线 | 97超碰精品成人国产| 中文字幕熟女人妻在线| 免费av观看视频| 乱人视频在线观看| 中文资源天堂在线| 激情 狠狠 欧美| 亚洲无线观看免费| 久久综合国产亚洲精品| 人妻制服诱惑在线中文字幕| 亚洲av免费高清在线观看| 最好的美女福利视频网| 精品人妻一区二区三区麻豆 | 精品久久久久久久久久免费视频| 偷拍熟女少妇极品色| 国产精品亚洲一级av第二区| 非洲黑人性xxxx精品又粗又长| 久久亚洲精品不卡| a级一级毛片免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲一区二区三区色噜噜| 美女高潮的动态| 亚洲一级一片aⅴ在线观看| 国产精品人妻久久久久久| 俄罗斯特黄特色一大片| 99视频精品全部免费 在线| 久久久久性生活片| 久久精品国产清高在天天线| 亚洲国产日韩欧美精品在线观看| 一级毛片我不卡| 此物有八面人人有两片| 国产精品永久免费网站| 91久久精品国产一区二区三区| 国产精品1区2区在线观看.| 国产麻豆成人av免费视频| 国产高清激情床上av| 天堂av国产一区二区熟女人妻| 99久久无色码亚洲精品果冻| 又粗又爽又猛毛片免费看| 国产久久久一区二区三区| 亚洲精品国产成人久久av| 变态另类成人亚洲欧美熟女| 韩国av在线不卡| 国产私拍福利视频在线观看| 亚洲性久久影院| 亚洲最大成人手机在线| 亚洲精品成人久久久久久| 精品久久久噜噜| 亚洲av熟女| 久久热精品热| 国产av一区在线观看免费| 国产精品日韩av在线免费观看| 网址你懂的国产日韩在线| 午夜激情欧美在线| 亚洲综合色惰| 99热精品在线国产| 99九九线精品视频在线观看视频| 少妇高潮的动态图| 午夜激情福利司机影院| 99久久中文字幕三级久久日本| 欧美人与善性xxx| 国产欧美日韩一区二区精品| 日韩中字成人| 麻豆一二三区av精品| 一边摸一边抽搐一进一小说| 久久人人精品亚洲av| 成人综合一区亚洲| 少妇猛男粗大的猛烈进出视频 | 激情 狠狠 欧美| av.在线天堂| 晚上一个人看的免费电影| 少妇的逼好多水| 免费av不卡在线播放| 成年免费大片在线观看| 日本精品一区二区三区蜜桃| 国产精品国产高清国产av| 在线免费观看不下载黄p国产| 国产精品,欧美在线| 国产精品99久久久久久久久| 久久久色成人| 搡女人真爽免费视频火全软件 | 有码 亚洲区| 三级男女做爰猛烈吃奶摸视频| 国产日本99.免费观看| 五月玫瑰六月丁香| 白带黄色成豆腐渣| 国产日本99.免费观看| 国产高清视频在线观看网站| 免费无遮挡裸体视频| 亚洲成a人片在线一区二区| 免费高清视频大片| 成年版毛片免费区| 夜夜夜夜夜久久久久| 男女做爰动态图高潮gif福利片| 黄色欧美视频在线观看| 舔av片在线| 男人的好看免费观看在线视频| 亚洲精华国产精华液的使用体验 | 成年女人永久免费观看视频| av在线亚洲专区| 免费人成视频x8x8入口观看| 欧美三级亚洲精品| 久久精品国产亚洲网站| 精品久久久久久久久亚洲| 国产在视频线在精品| 97在线视频观看|