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

    一種基于廣義延拓逼近的動(dòng)態(tài)濾波方法初步研究*

    2020-10-17 12:03:48高為廣
    天文研究與技術(shù) 2020年4期
    關(guān)鍵詞:歷元卡爾曼濾波廣義

    劉 成,李 芳,高為廣,王 威

    (1.北京跟蹤與通信技術(shù)研究所,北京 100094;2.中國(guó)科學(xué)院國(guó)家天文臺(tái),北京 100101)

    1 概 述

    廣義延拓逼近方法自1986年首次提出以來(lái)[1]展示了強(qiáng)大的生命力,應(yīng)用橫跨天文、力學(xué)、電磁學(xué)、機(jī)械、電子、通信和導(dǎo)航等多個(gè)領(lǐng)域,不斷證明了其科學(xué)性和有效性。起初,廣義延拓逼近方法的提出是為了解決大型天文望遠(yuǎn)鏡設(shè)計(jì)中的實(shí)際問(wèn)題,其后很快在一些光學(xué)與射電望遠(yuǎn)鏡、雷達(dá)天線、彈性力學(xué)以及電磁場(chǎng)理論等領(lǐng)域得到實(shí)踐,并被用于構(gòu)造有限單元內(nèi)的插值函數(shù)[2]。1991年,文[3]首次將廣義延拓逼近方法引入太陽(yáng)物理學(xué),證明了延拓邊界元空間的逼近性質(zhì)和邊界元近似解的收斂性,并于1993年出版[4]。2004年,文[5]重新對(duì)方法進(jìn)行了系統(tǒng)歸納和整理。隨后,廣義延拓逼近方法開始受到科研工作者更加廣泛的關(guān)注,并被推廣至電子[6]、電力[7-8]、機(jī)械[9-12]、特種設(shè)備[13]、氣象觀測(cè)[14]等工程領(lǐng)域以及導(dǎo)彈制導(dǎo)等軍事領(lǐng)域[15-16]。

    同時(shí),廣義延拓逼近方法開始在導(dǎo)航與定位領(lǐng)域廣泛應(yīng)用。文[17]將其用于全球定位系統(tǒng)信號(hào)遮擋情況下,對(duì)偽距測(cè)量值進(jìn)行外推,以維持用戶基本定位能力,結(jié)果表明,可在15 s內(nèi)維持優(yōu)于5 m的偽距外推精度。文[18]對(duì)接收機(jī)鐘漂進(jìn)行建模以增強(qiáng)定位,結(jié)果表明,除維持遮擋情況下的定位外,還可提高正常信號(hào)情況下約24%的定位精度。文[19-24]先后開展了廣義延拓逼近方法在導(dǎo)航衛(wèi)星廣播星歷和精密星歷插值中的應(yīng)用,均表明精度優(yōu)于一般插值方法。文[24]的研究表明,在選取合適參數(shù)時(shí),廣義延拓可達(dá)到毫米級(jí)的軌道插值精度,較拉格朗日(Lagrange)插值法高出56.6%~83.5%,較切比雪夫插值法高出20.9%~73.5%,且避免了傳統(tǒng)插值方法的高階振蕩現(xiàn)象。文[25-27]研究了廣義延拓逼近方法在衛(wèi)星精密鐘差中的插值應(yīng)用,獲得的絕大多數(shù)數(shù)據(jù)優(yōu)于0.03 ns、全部數(shù)據(jù)優(yōu)于0.1 ns的插值精度。文[28-29]利用廣義延拓進(jìn)行高程擬合,結(jié)果表明,精度優(yōu)于傳統(tǒng)二次多項(xiàng)式曲面法和多面函數(shù)法。文[30-33]將其用于導(dǎo)航信號(hào)的頻率及載波跟蹤,證明了廣義延拓方法具有更優(yōu)的信號(hào)跟蹤性能。文[34-35]將其用于對(duì)載波測(cè)量值周跳的檢測(cè)和修復(fù),不同的研究結(jié)果均表明,廣義延拓方法能夠準(zhǔn)確探測(cè)和修復(fù)1周以上的周跳。文[36-37]將其用于接收機(jī)碼鑒相器中的碼相位測(cè)量,結(jié)果表明廣義延拓方法具有更高的精度和穩(wěn)定性。文[38-41]探討了廣義延拓方法在衛(wèi)星導(dǎo)航以及室內(nèi)導(dǎo)航算法中的優(yōu)化應(yīng)用,進(jìn)一步拓展了其應(yīng)用場(chǎng)景。

    上述諸多研究和應(yīng)用均獲得了理想效果,證明了廣義延拓方法的科學(xué)性和正確性。然而,無(wú)論是軌道、鐘差、高程還是偽距、載波、碼相關(guān)峰,其特點(diǎn)都是一維、靜態(tài)數(shù)據(jù)下的插值或外推計(jì)算,局限于廣義延拓方法的一般形式。廣義延拓方法能利用插值多項(xiàng)式實(shí)現(xiàn)對(duì)狀態(tài)量的外推預(yù)測(cè),然而,將其與濾波方法相結(jié)合,應(yīng)用于動(dòng)態(tài)數(shù)據(jù)處理領(lǐng)域,需要解決其測(cè)量修正問(wèn)題并構(gòu)建相應(yīng)的誤差模型,目前鮮有相關(guān)研究性工作。本文在上述研究工作的基礎(chǔ)上對(duì)傳統(tǒng)廣義延拓方法進(jìn)行了拓展,提出將其推廣至動(dòng)態(tài)濾波數(shù)據(jù)處理領(lǐng)域,研究和構(gòu)造了基于廣義延拓插值方程的濾波模型與框架,并推導(dǎo)給出了求解方法與流程。所提出的廣義延拓逼近濾波方法既具備多維狀態(tài)量的數(shù)據(jù)處理能力,也能滿足狀態(tài)量動(dòng)態(tài)更新的要求,豐富了廣義延拓方法的內(nèi)涵。分別利用仿真數(shù)據(jù)和北斗衛(wèi)星導(dǎo)航系統(tǒng)(BeiDouSatellite Navigation System, BDS)實(shí)測(cè)數(shù)據(jù)對(duì)方法進(jìn)行了初步驗(yàn)證,結(jié)果表明,在高斯噪聲、非高斯噪聲以及實(shí)際導(dǎo)航應(yīng)用環(huán)境下均獲得了比卡爾曼濾波更好的性能,具有較大的發(fā)展?jié)摿蛻?yīng)用價(jià)值。

    2 廣義延拓逼近數(shù)學(xué)方法

    廣義延拓逼近方法的數(shù)學(xué)本質(zhì)是基于 “編結(jié)” 思想將數(shù)據(jù)集劃分為多個(gè)相互連結(jié)的單元延拓域,并分段進(jìn)行兼顧插值與擬合的數(shù)值處理,如圖1,記某一單元延拓域Δe內(nèi)的廣義延拓二次插值函數(shù)U(t)為

    圖1 延拓域及其對(duì)應(yīng)函數(shù)

    U(t)=a0+a1t+a2t2,t∈[t1,t2],

    (1)

    其中,ti為數(shù)據(jù)對(duì)應(yīng)的時(shí)刻;a0,a1,a2為廣義延拓多項(xiàng)式系數(shù),可由延拓域中的數(shù)學(xué)模型確定:

    (2)

    也就是說(shuō),一方面要求廣義延拓插值多項(xiàng)式在延拓域邊界點(diǎn)上滿足插值條件,保持邊界點(diǎn)的數(shù)據(jù)特征和約束作用,并使得各段數(shù)據(jù)之間的變化具有更好的協(xié)調(diào)性。另一方面,利用分段插值區(qū)域的周圍結(jié)點(diǎn)(包括內(nèi)點(diǎn))實(shí)現(xiàn)分段區(qū)域內(nèi)外部數(shù)據(jù)的殘差最小擬合,充分結(jié)合了插值和擬合方法的優(yōu)點(diǎn)。

    廣義延拓逼近方法能根據(jù)實(shí)際情況,靈活構(gòu)建具體的數(shù)學(xué)模型。例如,當(dāng)對(duì)某個(gè)邊界域節(jié)點(diǎn)的數(shù)值精度具有較高置信度時(shí),可以將其增加作為插值約束點(diǎn),以獲得更準(zhǔn)確的擬合結(jié)果;反之,可以放棄某個(gè)插值節(jié)點(diǎn),將其退化為擬合點(diǎn),降低對(duì)整體結(jié)果的影響,也可以采用三次或更高次的廣義延拓多項(xiàng)式,以獲得更高的數(shù)值處理精度。

    3 基于廣義延拓逼近的動(dòng)態(tài)濾波方法

    3.1 基本原理

    基于廣義延拓逼近的動(dòng)態(tài)濾波方法的基本思想是利用廣義延拓法構(gòu)造一個(gè)具有一定先驗(yàn)數(shù)據(jù)長(zhǎng)度、可靈活設(shè)置和優(yōu)選插值約束點(diǎn)的移動(dòng)窗口,對(duì)系統(tǒng)狀態(tài)量進(jìn)行高精度、具有抗差能力的逼近和預(yù)測(cè),以取代傳統(tǒng)卡爾曼濾波方法中的狀態(tài)量時(shí)間更新(狀態(tài)轉(zhuǎn)移)過(guò)程,如圖2。

    與一般卡爾曼濾波方法相比,廣義延拓逼近濾波方法使用了數(shù)據(jù)窗口,能夠起到類似事后平滑的效果??柭鼮V波的數(shù)學(xué)本質(zhì)是基于馬爾科夫鏈的一步狀態(tài)轉(zhuǎn)移,利用最新歷元n的狀態(tài)量預(yù)測(cè)下一歷元n+ 1的狀態(tài)量。它的優(yōu)點(diǎn)在于大幅減小了計(jì)算量,但同時(shí)也導(dǎo)致過(guò)于依賴當(dāng)前歷元狀態(tài)量的精度;換言之,即使根據(jù)誤差協(xié)方差矩陣P判斷出當(dāng)前歷元n的狀態(tài)量最優(yōu)估計(jì)值具有很大誤差,也只能繼續(xù)利用該狀態(tài)量估計(jì)值進(jìn)行下一步預(yù)測(cè)。相比之下,廣義延拓逼近濾波利用一個(gè)數(shù)據(jù)窗口進(jìn)行預(yù)測(cè),具有更好的平滑性和穩(wěn)定性。

    與一般的開窗濾波方法相比,廣義延拓逼近濾波方法可以在數(shù)據(jù)窗口內(nèi)靈活選擇插值約束點(diǎn),從而實(shí)現(xiàn)更好的抗差能力。例如,當(dāng)當(dāng)前歷元n的定位誤差較大時(shí),可將其退化作為擬合點(diǎn),以弱化其作用,并從窗口信息中另行選擇更加精確的狀態(tài)量估計(jì)值作為插值約束點(diǎn)(如圖2),從而在符合載體實(shí)際運(yùn)動(dòng)軌跡的同時(shí),充分發(fā)揮可靠先驗(yàn)信息的約束作用,規(guī)避粗差影響,提高狀態(tài)轉(zhuǎn)移的精度和可靠性。

    圖2 廣義延拓動(dòng)態(tài)濾波原理示意圖

    3.2 插值約束方程的構(gòu)造與狀態(tài)量的時(shí)間轉(zhuǎn)移

    3.2.1 基本形式

    若已知一段窗口長(zhǎng)度為L(zhǎng)的先驗(yàn)狀態(tài)量最優(yōu)估計(jì)值x^i(i=n-L+ 1,n-L+ 2, …,n),通過(guò)將窗口范圍內(nèi)的某一歷元j(j∈[n-L+ 1,n])的狀態(tài)量最優(yōu)估計(jì)值x^j作為插值點(diǎn),其他歷元的狀態(tài)量最優(yōu)估計(jì)值作為擬合點(diǎn),可建立一個(gè)具有插值約束條件的廣義延拓逼近方程:

    (3)

    (4)

    (5)

    (6)

    3.2.2 擴(kuò)展形式

    上節(jié)給出了廣義延拓逼近方程的一般形式。實(shí)際上,方程也可采用更高階的插值多項(xiàng)式形式:

    (7)

    或者可以進(jìn)一步增加和采用某另一歷元k(k∈[n-L+1,n]且k≠j)的狀態(tài)量最優(yōu)估計(jì)值,從而構(gòu)成具有多個(gè)插值約束條件的方程形式:

    (8)

    以上擴(kuò)展形式可根據(jù)具體需求靈活構(gòu)造。

    3.3 完整動(dòng)態(tài)濾波框架

    利用廣義延拓逼近方法完成狀態(tài)量的插值平滑更新后,結(jié)合測(cè)量更新處理步驟,即可建立完整的優(yōu)化動(dòng)態(tài)濾波框架模型,如圖3。

    圖3 廣義延拓動(dòng)態(tài)濾波模型

    (9)

    進(jìn)行更新,其中,Q為系統(tǒng)轉(zhuǎn)移過(guò)程噪聲;A為系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣。當(dāng)廣義延拓插值方程采用二次多項(xiàng)式時(shí),對(duì)應(yīng)狀態(tài)量是勻速變化模型,有

    (10)

    其中,Δt為兩時(shí)刻之間的時(shí)間間隔。當(dāng)廣義延拓插值方程采用三次多項(xiàng)式時(shí),對(duì)應(yīng)狀態(tài)量是勻加速變化模型,有

    (11)

    在卡爾曼濾波測(cè)量更新過(guò)程中,分別計(jì)算n+ 1歷元的權(quán)增益矩陣Kn+1、狀態(tài)量最優(yōu)估計(jì)值X^n+1及其協(xié)方差矩陣Pn+1:

    (12)

    其中,C為測(cè)量關(guān)系矩陣;R為測(cè)量噪聲協(xié)方差矩陣。

    4 求解方法與步驟

    對(duì)廣義延拓動(dòng)態(tài)濾波逼近方程的求解,可采用線性化展開的形式。具體地,將(3)式中的插值約束方程改寫為

    (13)

    將(13)式代入最小二乘公式可得

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    5 仿真計(jì)算

    5.1 信號(hào)去噪處理仿真

    分別對(duì)正弦信號(hào)在白噪聲和有色噪聲情況下的濾波性能進(jìn)行仿真。仿真時(shí)間T=0~30 s,步長(zhǎng)Δt=0.1 s,重復(fù)計(jì)算1 000次。利用蒙特卡洛模擬對(duì)精度結(jié)果進(jìn)行統(tǒng)計(jì),通過(guò)與廣泛應(yīng)用的卡爾曼濾波器進(jìn)行比較,對(duì)廣義延拓逼近濾波的性能進(jìn)行分析和評(píng)估。

    5.1.1 白噪聲濾波性能

    首先,加入滿足r~N(0, 0.1)分布的高斯白噪聲進(jìn)行數(shù)值仿真,考察和分析濾波方法的基本特性。兩種濾波器的初始狀態(tài)量為X0=[0 0],初始協(xié)方差為P0=diag([0.1 2]),初始過(guò)程噪聲Q0=[0.0 0.000 5; 0.000 5 0.01],初始測(cè)量噪聲為R0=0.01,廣義延拓逼近濾波的窗口長(zhǎng)度L=20,插值約束點(diǎn)為當(dāng)前歷元觀測(cè)點(diǎn)。兩種濾波方法的效果如圖4,誤差概率統(tǒng)計(jì)如圖5。其中,卡爾曼濾波的均方根誤差約0.083,廣義延拓逼近濾波的均方根誤差約0.069,提升約17%。

    圖4 白噪聲下的濾波仿真結(jié)果

    圖5 白噪聲下的濾波誤差概率統(tǒng)計(jì)

    5.1.2 有色噪聲濾波性能

    在上述高斯白噪聲基礎(chǔ)上,加入有色噪聲e(k)=0.2(x(k)+0.9x(k-1)),濾波器誤差參數(shù)和窗口數(shù)據(jù)長(zhǎng)度均保持不變。廣義延拓逼近方法和卡爾曼濾波方法均基于噪聲高斯分布的理想假設(shè),但在實(shí)際應(yīng)用中,理想高斯分布的情況很難真實(shí)存在,因此,考察兩種方法在有色噪聲影響下的性能具有實(shí)際意義。兩種濾波方法的效果如圖6,誤差概率統(tǒng)計(jì)如圖7。其中,卡爾曼濾波的均方根誤差約0.113,廣義延拓逼近濾波的均方根誤差約0.102,提升約10%。

    圖6 疊加有色噪聲后的濾波仿真結(jié)果

    圖7 疊加有色噪聲后的濾波誤差概率統(tǒng)計(jì)

    5.2 衛(wèi)星導(dǎo)航實(shí)測(cè)數(shù)據(jù)處理驗(yàn)證

    2020年1月,采用自主研制的全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System, GNSS)多模導(dǎo)航接收機(jī)在北京地區(qū)測(cè)量得到的北斗衛(wèi)星導(dǎo)航系統(tǒng)觀測(cè)數(shù)據(jù)(包括偽距、多普勒、載波和廣播星歷)對(duì)方法進(jìn)行實(shí)測(cè)驗(yàn)證,并與卡爾曼濾波器以及經(jīng)典希格(Sage)開窗濾波器處理結(jié)果進(jìn)行比較。幾種方法的定位輸出結(jié)果如圖8,其中,廣義延拓逼近濾波的窗口長(zhǎng)度取L=15。

    從圖8可以看出,廣義延拓濾波方法有能力獲得比經(jīng)典卡爾曼濾波和希格開窗濾波更優(yōu)的性能,具有應(yīng)用潛力。以實(shí)時(shí)動(dòng)態(tài)(Real-time Kinemaic, RTK)載波定位解算結(jié)果作為基準(zhǔn),統(tǒng)計(jì)3種方法的東西向平面誤差,結(jié)果見表1。

    表1 濾波定位誤差結(jié)果對(duì)比

    圖8 北斗導(dǎo)航濾波定位結(jié)果對(duì)比

    6 結(jié) 論

    本文在廣義延拓逼近理論的基礎(chǔ)上,提出了一種動(dòng)態(tài)濾波更新方法。研究和構(gòu)造了基于廣義延拓插值方程的濾波模型與框架,并推導(dǎo)給出了其求解方法與流程。方法既具備多維狀態(tài)量的處理能力,也能夠滿足高低階狀態(tài)量動(dòng)態(tài)更新的要求。在一段移動(dòng)平滑窗口內(nèi)構(gòu)造和求解廣義延拓插值多項(xiàng)式,實(shí)現(xiàn)對(duì)系統(tǒng)狀態(tài)量的時(shí)間轉(zhuǎn)移,具有更好的平滑效果,拓展了經(jīng)典廣義延拓?cái)?shù)學(xué)方法的內(nèi)涵。與一般開窗濾波方法相比,廣義延拓逼近濾波方法數(shù)據(jù)窗口內(nèi)插值點(diǎn)的位置及數(shù)量均可靈活選取,從而可通過(guò)鎖定最新觀測(cè)點(diǎn)或最準(zhǔn)確的先驗(yàn)觀測(cè)點(diǎn),實(shí)現(xiàn)更高的精度和抗差能力。分別利用正弦信號(hào)數(shù)值仿真和北斗衛(wèi)星導(dǎo)航實(shí)測(cè)數(shù)據(jù)對(duì)方法進(jìn)行了驗(yàn)證,結(jié)果表明,在高斯噪聲和非高斯噪聲仿真下均獲得了比卡爾曼濾波更優(yōu)的性能,驗(yàn)證了將廣義延拓逼近理論與濾波方法相結(jié)合的可行性。另外需要指出的是,目前初步構(gòu)造廣義延拓濾波模型仍是基于卡爾曼濾波誤差假設(shè)條件的近似,相信隨著后續(xù)理論和模型的完善,方法將具備更優(yōu)的性能與更大的應(yīng)用價(jià)值。

    猜你喜歡
    歷元卡爾曼濾波廣義
    Rn中的廣義逆Bonnesen型不等式
    歷元間載波相位差分的GPS/BDS精密單點(diǎn)測(cè)速算法
    從廣義心腎不交論治慢性心力衰竭
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    Recent advances of TCM treatment of childhood atopic dermatitis
    有限群的廣義交換度
    Clinical observation of Huatan Huoxue Formula in treating coronary heart disease with hyperlipidemia
    Mechanism of sex hormone level in biological clock disorder induced acne and analysis of TCM Pathogenesis
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    基于擴(kuò)展卡爾曼濾波的PMSM無(wú)位置傳感器控制
    搡老妇女老女人老熟妇| 啪啪无遮挡十八禁网站| 啦啦啦 在线观看视频| 亚洲aⅴ乱码一区二区在线播放 | 在线免费观看的www视频| 99国产综合亚洲精品| 亚洲一卡2卡3卡4卡5卡精品中文| 91麻豆av在线| 欧美国产精品va在线观看不卡| 亚洲三区欧美一区| 日韩 欧美 亚洲 中文字幕| 大陆偷拍与自拍| 亚洲五月婷婷丁香| 久久精品国产99精品国产亚洲性色 | 国产精品美女特级片免费视频播放器 | 天堂√8在线中文| 免费在线观看影片大全网站| 亚洲专区国产一区二区| 国产熟女xx| 在线国产一区二区在线| 午夜免费鲁丝| 波多野结衣巨乳人妻| 啦啦啦韩国在线观看视频| 亚洲中文日韩欧美视频| 亚洲自偷自拍图片 自拍| 欧美日韩瑟瑟在线播放| 成人永久免费在线观看视频| 亚洲欧美精品综合一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 久热这里只有精品99| 久久精品国产清高在天天线| 国产亚洲精品第一综合不卡| 国产av精品麻豆| 国产国语露脸激情在线看| 亚洲情色 制服丝袜| 好看av亚洲va欧美ⅴa在| 国产亚洲av嫩草精品影院| 18禁美女被吸乳视频| 国产午夜精品久久久久久| 亚洲第一青青草原| 99热只有精品国产| 变态另类丝袜制服| 叶爱在线成人免费视频播放| 欧美日韩乱码在线| 18美女黄网站色大片免费观看| 欧美成人性av电影在线观看| 黄片大片在线免费观看| 亚洲成av人片免费观看| 中亚洲国语对白在线视频| 18禁裸乳无遮挡免费网站照片 | 黄色片一级片一级黄色片| 久久人人爽av亚洲精品天堂| 亚洲精品在线美女| 国产一区在线观看成人免费| 99国产精品一区二区蜜桃av| 成熟少妇高潮喷水视频| 一本大道久久a久久精品| 神马国产精品三级电影在线观看 | 欧美日本亚洲视频在线播放| 露出奶头的视频| 18禁黄网站禁片午夜丰满| 精品久久久久久久人妻蜜臀av | 国产色视频综合| 最新美女视频免费是黄的| 在线观看66精品国产| 久热这里只有精品99| 午夜福利在线观看吧| 日韩欧美在线二视频| 日韩精品中文字幕看吧| 国产蜜桃级精品一区二区三区| 电影成人av| 黄色丝袜av网址大全| 在线观看66精品国产| 亚洲激情在线av| 久久久久国内视频| 亚洲欧美日韩高清在线视频| 好男人电影高清在线观看| 国产精品九九99| 一进一出抽搐动态| 久久久久久久久中文| 日韩欧美国产一区二区入口| 91av网站免费观看| 校园春色视频在线观看| 夜夜看夜夜爽夜夜摸| 欧美激情 高清一区二区三区| 欧美乱妇无乱码| 50天的宝宝边吃奶边哭怎么回事| 亚洲中文字幕一区二区三区有码在线看 | 国产区一区二久久| 女人爽到高潮嗷嗷叫在线视频| 神马国产精品三级电影在线观看 | 久久久久久久午夜电影| 久久精品国产99精品国产亚洲性色 | 在线视频色国产色| 1024视频免费在线观看| 91老司机精品| 精品人妻1区二区| 欧美日韩亚洲综合一区二区三区_| 最新在线观看一区二区三区| 午夜福利影视在线免费观看| 午夜福利一区二区在线看| 1024视频免费在线观看| 亚洲 欧美一区二区三区| www.自偷自拍.com| 99国产精品一区二区蜜桃av| 久久久久久久久免费视频了| 久久久久亚洲av毛片大全| 国产av精品麻豆| 婷婷丁香在线五月| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人免费无遮挡视频| 国产精品久久久久久亚洲av鲁大| 免费人成视频x8x8入口观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲专区中文字幕在线| 成人亚洲精品一区在线观看| 欧美另类亚洲清纯唯美| 亚洲 欧美 日韩 在线 免费| 天天添夜夜摸| √禁漫天堂资源中文www| 亚洲欧美日韩另类电影网站| 久久久久久国产a免费观看| 精品无人区乱码1区二区| 香蕉久久夜色| 1024香蕉在线观看| 一进一出抽搐动态| 老司机在亚洲福利影院| 十八禁网站免费在线| 久99久视频精品免费| 在线国产一区二区在线| 免费少妇av软件| 母亲3免费完整高清在线观看| 亚洲精品国产色婷婷电影| 亚洲精品在线美女| 淫秽高清视频在线观看| 美女扒开内裤让男人捅视频| 每晚都被弄得嗷嗷叫到高潮| 久久人妻熟女aⅴ| 午夜久久久在线观看| 精品不卡国产一区二区三区| 国产精品爽爽va在线观看网站 | 精品第一国产精品| 亚洲少妇的诱惑av| 电影成人av| 人妻丰满熟妇av一区二区三区| 亚洲伊人色综图| 亚洲七黄色美女视频| 中文字幕最新亚洲高清| 此物有八面人人有两片| 免费久久久久久久精品成人欧美视频| 最新美女视频免费是黄的| 成人永久免费在线观看视频| 韩国精品一区二区三区| 成人av一区二区三区在线看| 日韩免费av在线播放| 日韩欧美一区视频在线观看| 亚洲精品粉嫩美女一区| 最近最新免费中文字幕在线| 极品人妻少妇av视频| 正在播放国产对白刺激| 免费在线观看视频国产中文字幕亚洲| 妹子高潮喷水视频| 宅男免费午夜| 性色av乱码一区二区三区2| 午夜福利在线观看吧| 多毛熟女@视频| 亚洲一区中文字幕在线| 每晚都被弄得嗷嗷叫到高潮| 黄片小视频在线播放| 午夜亚洲福利在线播放| 久久九九热精品免费| 国产精品一区二区免费欧美| 国产亚洲精品久久久久久毛片| 亚洲一区二区三区色噜噜| 在线观看一区二区三区| 叶爱在线成人免费视频播放| 操出白浆在线播放| 国产精品一区二区免费欧美| 999久久久精品免费观看国产| 成人亚洲精品一区在线观看| xxx96com| 亚洲少妇的诱惑av| 国产av又大| 精品久久久精品久久久| 在线免费观看的www视频| 极品人妻少妇av视频| av免费在线观看网站| 99在线人妻在线中文字幕| 亚洲第一欧美日韩一区二区三区| 国产成人欧美| 精品国产乱码久久久久久男人| 国产精品一区二区在线不卡| 欧美av亚洲av综合av国产av| 国产国语露脸激情在线看| 亚洲自偷自拍图片 自拍| 亚洲欧美日韩无卡精品| 制服人妻中文乱码| 淫妇啪啪啪对白视频| 男男h啪啪无遮挡| 丝袜在线中文字幕| 狠狠狠狠99中文字幕| 欧美久久黑人一区二区| 国产精品久久久av美女十八| 在线观看免费午夜福利视频| 一级黄色大片毛片| 一夜夜www| 中文字幕久久专区| 男女做爰动态图高潮gif福利片 | 天天一区二区日本电影三级 | 久久精品国产综合久久久| 国产精品久久久久久精品电影 | 一级毛片高清免费大全| 一进一出好大好爽视频| 久久久久国产一级毛片高清牌| 亚洲三区欧美一区| 满18在线观看网站| 亚洲中文日韩欧美视频| 他把我摸到了高潮在线观看| www国产在线视频色| 亚洲三区欧美一区| 亚洲精品中文字幕一二三四区| 亚洲在线自拍视频| av视频免费观看在线观看| 国产国语露脸激情在线看| 午夜激情av网站| 久久久久久国产a免费观看| 久久伊人香网站| 999久久久国产精品视频| 母亲3免费完整高清在线观看| 男男h啪啪无遮挡| 国产精品久久久久久精品电影 | 欧美性长视频在线观看| 女性生殖器流出的白浆| 91老司机精品| 久久天堂一区二区三区四区| 操美女的视频在线观看| 亚洲专区中文字幕在线| 制服诱惑二区| 在线观看66精品国产| 99国产综合亚洲精品| 91老司机精品| 国产私拍福利视频在线观看| 亚洲国产精品久久男人天堂| 日本欧美视频一区| 精品一区二区三区四区五区乱码| 亚洲一区二区三区色噜噜| cao死你这个sao货| 99国产精品99久久久久| 久久精品影院6| 色尼玛亚洲综合影院| 国产三级在线视频| 男人操女人黄网站| 夜夜看夜夜爽夜夜摸| 男女之事视频高清在线观看| 咕卡用的链子| 午夜福利免费观看在线| 啦啦啦观看免费观看视频高清 | 免费在线观看日本一区| 1024香蕉在线观看| 久久伊人香网站| 午夜福利18| 欧美色视频一区免费| 国产精品一区二区精品视频观看| 日韩欧美免费精品| 中文字幕高清在线视频| 天堂√8在线中文| 亚洲第一电影网av| 精品久久久久久,| 日韩有码中文字幕| 亚洲一区二区三区色噜噜| av超薄肉色丝袜交足视频| 久久精品国产综合久久久| 一二三四社区在线视频社区8| 老司机福利观看| 亚洲精品在线观看二区| 无人区码免费观看不卡| 丝袜在线中文字幕| 国产精华一区二区三区| 精品福利观看| 大型av网站在线播放| 精品久久久精品久久久| 国产精品自产拍在线观看55亚洲| 久久 成人 亚洲| 日韩欧美三级三区| 久久国产精品人妻蜜桃| 国产精品久久久久久精品电影 | 欧美 亚洲 国产 日韩一| 黄色视频不卡| 国内精品久久久久精免费| 午夜福利一区二区在线看| 中出人妻视频一区二区| 亚洲自拍偷在线| 欧美日韩福利视频一区二区| 亚洲一区二区三区色噜噜| 一级,二级,三级黄色视频| 黄色毛片三级朝国网站| 长腿黑丝高跟| 久久青草综合色| 一进一出好大好爽视频| 黄色丝袜av网址大全| 亚洲国产看品久久| 精品久久蜜臀av无| 亚洲电影在线观看av| av有码第一页| 美女国产高潮福利片在线看| 一进一出抽搐gif免费好疼| 国产成年人精品一区二区| 国产成人一区二区三区免费视频网站| 久久精品亚洲熟妇少妇任你| 国产免费av片在线观看野外av| 十八禁人妻一区二区| 免费一级毛片在线播放高清视频 | 青草久久国产| 热re99久久国产66热| 久久狼人影院| 午夜免费激情av| 国产亚洲精品综合一区在线观看 | 麻豆一二三区av精品| 中文字幕人成人乱码亚洲影| 欧美乱妇无乱码| 亚洲一区二区三区不卡视频| 巨乳人妻的诱惑在线观看| 日本免费a在线| 亚洲av片天天在线观看| 免费高清在线观看日韩| 午夜视频精品福利| 亚洲色图 男人天堂 中文字幕| 亚洲专区国产一区二区| 久久久久国内视频| 成人特级黄色片久久久久久久| 亚洲黑人精品在线| 9热在线视频观看99| 欧美成人性av电影在线观看| 亚洲午夜精品一区,二区,三区| 精品免费久久久久久久清纯| 岛国在线观看网站| 两个人视频免费观看高清| 国产熟女xx| 嫩草影院精品99| 国内久久婷婷六月综合欲色啪| 12—13女人毛片做爰片一| 午夜福利视频1000在线观看 | 欧美日韩亚洲综合一区二区三区_| 精品一区二区三区av网在线观看| 亚洲激情在线av| 高清在线国产一区| 欧美激情 高清一区二区三区| 日韩欧美国产在线观看| 91麻豆精品激情在线观看国产| 丰满人妻熟妇乱又伦精品不卡| 国产高清激情床上av| 亚洲欧美精品综合久久99| 男人舔女人的私密视频| 亚洲av日韩精品久久久久久密| 美女午夜性视频免费| 久久人妻熟女aⅴ| 一卡2卡三卡四卡精品乱码亚洲| 777久久人妻少妇嫩草av网站| 国产精品,欧美在线| 非洲黑人性xxxx精品又粗又长| 亚洲精品一卡2卡三卡4卡5卡| av有码第一页| 欧美中文综合在线视频| 成人国产一区最新在线观看| 不卡av一区二区三区| 精品人妻在线不人妻| bbb黄色大片| 一边摸一边做爽爽视频免费| 免费搜索国产男女视频| 亚洲五月色婷婷综合| 亚洲精品久久国产高清桃花| 亚洲人成网站在线播放欧美日韩| 这个男人来自地球电影免费观看| 国产精品免费视频内射| 麻豆av在线久日| 免费无遮挡裸体视频| 亚洲精品中文字幕一二三四区| 国产熟女xx| 俄罗斯特黄特色一大片| 精品高清国产在线一区| 淫秽高清视频在线观看| 色播在线永久视频| 国产一区二区激情短视频| 亚洲欧美日韩高清在线视频| 亚洲熟妇中文字幕五十中出| 久久人人爽av亚洲精品天堂| 国产精品二区激情视频| 亚洲欧美日韩高清在线视频| 亚洲精品美女久久久久99蜜臀| xxx96com| 91精品三级在线观看| 国产人伦9x9x在线观看| 男女下面进入的视频免费午夜 | 亚洲中文av在线| xxx96com| 亚洲久久久国产精品| 黄色丝袜av网址大全| 欧美黑人精品巨大| 色尼玛亚洲综合影院| 国产一卡二卡三卡精品| 精品国内亚洲2022精品成人| 免费高清在线观看日韩| www日本在线高清视频| 最近最新中文字幕大全免费视频| 亚洲五月色婷婷综合| 亚洲五月婷婷丁香| 天天躁夜夜躁狠狠躁躁| 性欧美人与动物交配| 美女大奶头视频| 午夜久久久在线观看| 两个人免费观看高清视频| 中文字幕最新亚洲高清| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区免费| 51午夜福利影视在线观看| 欧美精品啪啪一区二区三区| 亚洲精品久久国产高清桃花| 一区二区三区高清视频在线| 精品国产一区二区久久| 久久国产乱子伦精品免费另类| 中文字幕色久视频| 婷婷六月久久综合丁香| 9色porny在线观看| 国产精品九九99| 在线观看舔阴道视频| 日韩欧美免费精品| 一二三四在线观看免费中文在| 丝袜人妻中文字幕| 亚洲成人久久性| 好男人在线观看高清免费视频 | 日韩视频一区二区在线观看| 国产成+人综合+亚洲专区| 男女下面插进去视频免费观看| 成人国产综合亚洲| 国产成人影院久久av| 国产99白浆流出| 国产成人av教育| 亚洲成人精品中文字幕电影| 亚洲国产中文字幕在线视频| 99国产精品免费福利视频| 侵犯人妻中文字幕一二三四区| 一本久久中文字幕| 成人av一区二区三区在线看| aaaaa片日本免费| 久久国产精品男人的天堂亚洲| 亚洲无线在线观看| 色综合婷婷激情| 嫩草影院精品99| 一区二区三区国产精品乱码| 在线天堂中文资源库| 女人精品久久久久毛片| 日本三级黄在线观看| 久久久久亚洲av毛片大全| 禁无遮挡网站| 亚洲色图av天堂| 精品久久久久久久毛片微露脸| 一级毛片女人18水好多| avwww免费| 99久久综合精品五月天人人| 国产精品久久久久久亚洲av鲁大| 男人操女人黄网站| 亚洲中文字幕日韩| 久久亚洲真实| 久久久久国产精品人妻aⅴ院| 91av网站免费观看| 欧美日本视频| 91成年电影在线观看| 免费在线观看日本一区| 在线免费观看的www视频| 中文字幕最新亚洲高清| 男女做爰动态图高潮gif福利片 | 免费观看精品视频网站| 欧美老熟妇乱子伦牲交| 亚洲男人的天堂狠狠| 成人三级黄色视频| 女人精品久久久久毛片| 亚洲精品一区av在线观看| 精品无人区乱码1区二区| 精品电影一区二区在线| 淫秽高清视频在线观看| 亚洲一区中文字幕在线| 制服诱惑二区| 亚洲自拍偷在线| 高清在线国产一区| 视频区欧美日本亚洲| 亚洲国产精品合色在线| 真人一进一出gif抽搐免费| 亚洲av成人不卡在线观看播放网| 老司机福利观看| 欧美黑人欧美精品刺激| 国产成+人综合+亚洲专区| 免费看十八禁软件| 日韩精品青青久久久久久| 午夜福利成人在线免费观看| 亚洲av成人av| 国产免费男女视频| 一本久久中文字幕| 免费看a级黄色片| 色尼玛亚洲综合影院| 亚洲精品国产精品久久久不卡| 日韩国内少妇激情av| 视频在线观看一区二区三区| 日韩精品青青久久久久久| 精品久久久久久成人av| 亚洲精品美女久久久久99蜜臀| 国产伦一二天堂av在线观看| 久久精品成人免费网站| 日韩欧美国产一区二区入口| 国内久久婷婷六月综合欲色啪| av有码第一页| 国产乱人伦免费视频| 在线观看免费午夜福利视频| 在线观看一区二区三区| 精品国产乱子伦一区二区三区| 国产精品精品国产色婷婷| 侵犯人妻中文字幕一二三四区| 午夜老司机福利片| 国产蜜桃级精品一区二区三区| 久久人人爽av亚洲精品天堂| 少妇裸体淫交视频免费看高清 | 香蕉丝袜av| 后天国语完整版免费观看| 亚洲 欧美 日韩 在线 免费| 免费看十八禁软件| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 麻豆成人av在线观看| 午夜免费激情av| 老汉色∧v一级毛片| 国产av又大| 欧美+亚洲+日韩+国产| 桃色一区二区三区在线观看| 午夜福利高清视频| 久久香蕉激情| 午夜精品国产一区二区电影| 精品欧美国产一区二区三| 日韩中文字幕欧美一区二区| 午夜精品在线福利| 一区二区三区高清视频在线| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 欧美乱妇无乱码| 身体一侧抽搐| tocl精华| 国产高清视频在线播放一区| 久久这里只有精品19| 亚洲精品在线观看二区| 很黄的视频免费| 99国产极品粉嫩在线观看| 日本vs欧美在线观看视频| 欧美一区二区精品小视频在线| 好男人在线观看高清免费视频 | svipshipincom国产片| 国产不卡一卡二| 欧美大码av| 久久中文字幕一级| 亚洲精品国产区一区二| 成人三级做爰电影| 国产一区二区三区综合在线观看| 岛国在线观看网站| 无人区码免费观看不卡| 国产私拍福利视频在线观看| 又黄又粗又硬又大视频| 给我免费播放毛片高清在线观看| 国产xxxxx性猛交| 中文字幕久久专区| 母亲3免费完整高清在线观看| 韩国精品一区二区三区| 久久影院123| 精品免费久久久久久久清纯| 纯流量卡能插随身wifi吗| www.精华液| 成人免费观看视频高清| 欧美成人午夜精品| 成年女人毛片免费观看观看9| 99国产精品免费福利视频| 成人国产综合亚洲| 亚洲人成伊人成综合网2020| 色综合亚洲欧美另类图片| 成人永久免费在线观看视频| 国产熟女午夜一区二区三区| 国产不卡一卡二| 在线国产一区二区在线| 久久性视频一级片| www.自偷自拍.com| 亚洲情色 制服丝袜| 大型av网站在线播放| 成人国产一区最新在线观看| 久久影院123| 欧美成人性av电影在线观看| 成人国产一区最新在线观看| 久久热在线av| 日韩av在线大香蕉| 亚洲精品在线观看二区| 50天的宝宝边吃奶边哭怎么回事| 国产精品久久久久久人妻精品电影| 黑人巨大精品欧美一区二区蜜桃| 久久热在线av| 亚洲国产精品成人综合色| 婷婷精品国产亚洲av在线| 50天的宝宝边吃奶边哭怎么回事| 一个人观看的视频www高清免费观看 | 中文字幕人妻丝袜一区二区| 欧美日韩瑟瑟在线播放| 很黄的视频免费| www.999成人在线观看| 亚洲午夜精品一区,二区,三区| 日本vs欧美在线观看视频| 一二三四社区在线视频社区8| 亚洲精品一卡2卡三卡4卡5卡| 国产精品亚洲av一区麻豆| 久久精品影院6| 国产亚洲av高清不卡| 精品日产1卡2卡|