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

    高超聲速飛行器自適應(yīng)魯棒跟蹤濾波算法

    2023-06-23 09:49:32梁新茹高長(zhǎng)生荊武興
    宇航學(xué)報(bào) 2023年5期
    關(guān)鍵詞:后驗(yàn)超聲速先驗(yàn)

    梁新茹,高長(zhǎng)生,荊武興

    (哈爾濱工業(yè)大學(xué)航天學(xué)院,哈爾濱 150001)

    0 引 言

    近些年來,高超聲速飛行器因其自身具有機(jī)動(dòng)能力強(qiáng)、機(jī)動(dòng)范圍廣、能夠長(zhǎng)時(shí)飛行等特點(diǎn),引起了越來越廣泛的關(guān)注[1-3]。與傳統(tǒng)飛行器相比,高超聲速飛行器具有飛行空域復(fù)雜、飛行范圍廣等特性[4-5],導(dǎo)致在高超聲速飛行器跟蹤過程中探測(cè)噪聲獲取的統(tǒng)計(jì)信息不準(zhǔn)確,同時(shí)探測(cè)噪聲容易產(chǎn)生異常值,使得實(shí)際測(cè)量噪聲呈現(xiàn)出非平穩(wěn)非高斯特性。

    針對(duì)實(shí)際噪聲呈現(xiàn)非平穩(wěn)非高斯這一問題,文獻(xiàn)[6]指出閃爍噪聲作為非高斯噪聲中的一種特殊分布,通常是由異常值引起的,其分布具有重尾特性[7]。但是在傳統(tǒng)的濾波框架內(nèi),通常假設(shè)噪聲分布為高斯分布[8-9]。針對(duì)探測(cè)過程中出現(xiàn)的異常值問題,先驗(yàn)的高斯噪聲建模將會(huì)損失建模的精度,使得建模成本和建模誤差大大增加。針對(duì)重尾分布建模問題,目前主要采用學(xué)生t (Student’s t, ST)分布及其魯棒算法進(jìn)行描述。文獻(xiàn)[10]指出多元ST分布能夠?yàn)榫哂挟惓V档臄?shù)據(jù)提供更為完善的描述,ST分布可以把高斯分布當(dāng)作一種特殊情況包含在內(nèi)。文獻(xiàn)[11]提供一種基于ST分布近似噪聲分布的后驗(yàn)概率密度函數(shù)(Probability density function, PDF)的濾波算法。文獻(xiàn)[12]設(shè)計(jì)了一種ST分布與粒子濾波相結(jié)合的策略。但是對(duì)于非平穩(wěn)重尾噪聲分布來說,單一的ST分布難以準(zhǔn)確地描述。文獻(xiàn)[6,13]提出一種混合概率可自適應(yīng)變化的高斯-學(xué)生t混合(Gaussian-student’s t mixture, GSTM)分布以適應(yīng)噪聲的非平穩(wěn)變化。但是上述文獻(xiàn)均沒有考慮噪聲先驗(yàn)信息未知的情況,且需要噪聲的先驗(yàn)信息,才能保證算法的成立。

    高超聲速飛行器飛行空域廣,空域環(huán)境復(fù)雜,導(dǎo)致噪聲的先驗(yàn)統(tǒng)計(jì)信息復(fù)雜,噪聲可能是時(shí)變、多峰的,無法給出準(zhǔn)確的先驗(yàn)統(tǒng)計(jì)信息。針對(duì)這一問題,文獻(xiàn)[14-16]提出利用逆Wishart分布對(duì)噪聲的協(xié)方差進(jìn)行建模,以解決噪聲統(tǒng)計(jì)特性不精確的問題。文獻(xiàn)[17]在協(xié)方差建模的基礎(chǔ)上,進(jìn)一步提出將噪聲的均值和協(xié)方差聯(lián)合分布近似為獨(dú)立的高斯-逆Wishart分布。文獻(xiàn)[18]將逆Wishart分布和ST建模結(jié)合,文獻(xiàn)[19]將高斯-逆Wishart分布應(yīng)用在ST建模中,完成對(duì)噪聲的建模。同時(shí),為了進(jìn)一步提高濾波算法對(duì)于時(shí)變且先驗(yàn)信息未知噪聲的魯棒性,需要解決ST分布中自由度采用固定常值這一問題,文獻(xiàn)[20]不再采用固定的自由度參數(shù),而是用伽馬函數(shù)來刻畫自由度參數(shù),提高了算法對(duì)于測(cè)量異常值的魯棒性。

    本文針對(duì)復(fù)雜空域飛行的高超聲速飛行器探測(cè)過程出現(xiàn)的噪聲異常值干擾,以及異常值的隨機(jī)不確定性導(dǎo)致量測(cè)噪聲的非平穩(wěn)非高斯與統(tǒng)計(jì)信息未知一系列問題,在GSTM分布的基礎(chǔ)上提出了一種改進(jìn)魯棒高斯-學(xué)生t混合分布濾波(Robust Gaussian-student’s t mixture distribution filtering, RGSTMF)算法。首先,引入了GSTM分布,使得算法能夠自適應(yīng)地調(diào)整模型概率來適應(yīng)高超聲速飛行器探測(cè)噪聲的非平穩(wěn)特性。其次,結(jié)合高斯-逆Wishart分布,利用伽馬函數(shù)對(duì)GSTM分布進(jìn)行改進(jìn),改進(jìn)后的GSTM分布可以利用高斯-逆Wishart分布代替GSTM中的零均值假設(shè)來描述實(shí)際應(yīng)用中的均值未知且時(shí)變的噪聲。伽馬函數(shù)代替了原有GSTM分布中固定的自由度以解決噪聲非高斯程度不確定的問題。因此,改進(jìn)后的GSTM分布可以實(shí)現(xiàn)在測(cè)量噪聲先驗(yàn)統(tǒng)計(jì)信息未知的情況下,對(duì)噪聲進(jìn)行準(zhǔn)確建模。最后,將改進(jìn)后的GSTM分布與變分貝葉斯算法相結(jié)合,推導(dǎo)出RGSTMF濾波算法中參數(shù)后驗(yàn)更新的顯式表達(dá)式,從而構(gòu)成完整的濾波算法。

    1 高超聲速飛行器助推滑翔段建模

    本文主要對(duì)高超聲速飛行器助推滑翔段的跟蹤問題進(jìn)行研究,以HTV-2飛行器作為典型代表。參考文獻(xiàn)[21-22],在彈道坐標(biāo)系下對(duì)高超聲速飛行器滑翔段建模,得到公式如下:

    (1)

    (2)

    式中:關(guān)于升力系數(shù)CL和阻力系數(shù)CD的計(jì)算,本文引用文獻(xiàn)[23]中的數(shù)據(jù),將CL和CD擬合成馬赫數(shù)Ma和攻角αT的函數(shù),如下:

    CL=-0.177 3+0.046 4αT+0.085 6e0.032 1Ma+

    (3)

    (4)

    本文的大氣密度采用指數(shù)模型

    ρT(h)=ρ0e-h/H

    (5)

    式中:ρ0取值為ρ0=1.46×10-3kg/m3;H=6 970 m。

    地球自轉(zhuǎn)相關(guān)的高階項(xiàng)如下:

    (6)

    式中:ωe表示地球自轉(zhuǎn)角速率的大小。

    2 RGSTMF算法

    2.1 GSTM分布

    在高超聲速飛行器實(shí)際跟蹤過程中,由于目標(biāo)特性和飛行環(huán)境的復(fù)雜性,跟蹤量測(cè)值極易出現(xiàn)異常值,這將導(dǎo)致噪聲分布的非高斯閃爍特性。因此引入GSTM分布對(duì)噪聲進(jìn)行建模,GSTM分布公式如下:

    pN(x)=N(x;μ,Σ)

    (7)

    pST(x)=ST(x;μ,Σ,v)

    (8)

    pGSTM(x)=δpN(x)+(1-δ)pST(x)

    (9)

    式中:N(·)表示高斯分布;μ表示分布的均值;Σ表示分布的協(xié)方差;ST(·)表示ST分布;v表示ST分布的自由度;δ表示混合概率。

    參考文獻(xiàn)[6],混合概率δ滿足

    p(δ)=B(δ;e,1-e)

    (10)

    式中:B(·)表示Beta分布;e表示Beta的先驗(yàn)形狀參數(shù)。

    從式(9)中可以看出當(dāng)混合概率δ=1時(shí),GSTM分布服從高斯分布,當(dāng)混合概率δ=0時(shí),GSTM分布表現(xiàn)為ST分布,因此高斯分布和ST分布是GSTM分布的一種特殊形式。GSTM分布在應(yīng)對(duì)非平穩(wěn)噪聲方面具有一定的優(yōu)勢(shì),具有更加廣泛的應(yīng)用。

    引入隱變量ξ,式(9)可以被改寫為

    pGSTM(x|ξ)=δN(x;μ,Σ)+(1-δ)N(x;μ,Σ/ξ)

    (11)

    式中:隱變量ξ滿足

    (12)

    因?yàn)樵肼暯y(tǒng)計(jì)特性未知且時(shí)變,固定自由度參數(shù)v的使用會(huì)導(dǎo)致噪聲建模的不準(zhǔn)確性,因此,將自由度參數(shù)v建模成為伽馬函數(shù)[18]。

    引用伯努利隨機(jī)變量y,式(11)可以被改寫為

    pGSTM(x|ξ,y)=[N(x;μ,Σ)]y+

    [N(x;μ,Σ/ξ)](1-y)

    (13)

    式中:伯努利隨機(jī)變量y滿足

    p(y|δ)=δy(1-δ)y

    (14)

    因此,式(9)可以寫作

    p(ξ|v)p(δ)p(v)dξdys.t.y∈{0, 1}

    (15)

    式中:關(guān)于隱變量ξ的積分為定積分,關(guān)于伯努利隨機(jī)變量y的積分為不定積分,使得最終的積分結(jié)果為y的函數(shù),用于積分外側(cè)的求和運(yùn)算。

    2.2 問題描述

    構(gòu)建高超聲速飛行器非線性跟蹤系統(tǒng):

    xk=f(xk-1)+ωk-1

    (16)

    zk=h(xk)+υk

    (17)

    式中:xk∈Rn表示k時(shí)刻的狀態(tài)向量;f(·)表示狀態(tài)函數(shù);h(·)表示量測(cè)函數(shù);zk∈Rm表示k時(shí)刻的量測(cè)向量;wk-1表示系統(tǒng)噪聲,滿足高斯分布;υk表示量測(cè)噪聲。

    對(duì)于高超聲速飛行器來說,由于飛行器飛行空域復(fù)雜,目標(biāo)特性復(fù)雜,在跟蹤過程中,容易出現(xiàn)異常值,這會(huì)導(dǎo)致目標(biāo)的量測(cè)噪聲出現(xiàn)非高斯特性。同時(shí)由于異常值的非平穩(wěn)特性,導(dǎo)致目標(biāo)量測(cè)噪聲存在統(tǒng)計(jì)特性未知的情況。因此,在實(shí)際應(yīng)用過程中,量測(cè)噪聲υk的表達(dá)形式和統(tǒng)計(jì)信息沒有先驗(yàn)的假設(shè),量測(cè)噪聲υk可能是非高斯分布,甚至可能是多峰的、傾斜的分布。本文的目的是在沒有先驗(yàn)統(tǒng)計(jì)特性信息的條件下,對(duì)量測(cè)噪聲進(jìn)行建模,在此基礎(chǔ)上推導(dǎo)完成濾波算法。

    2.3 先驗(yàn)分布的選擇

    根據(jù)貝葉斯公式,為了求解濾波的后驗(yàn)狀態(tài)估計(jì),需要對(duì)一步預(yù)測(cè)和似然函數(shù)進(jìn)行建模。

    由于過程噪聲的高斯特性,將一步預(yù)測(cè)建模為高斯分布:

    p(xk|z1:k-1)=N(xk;xk/k-1,Pk/k-1)

    (18)

    式中:xk/k-1表示k時(shí)刻狀態(tài)估計(jì)的一步預(yù)測(cè)均值;Pk/k-1表示k時(shí)刻狀態(tài)估計(jì)的一步預(yù)測(cè)協(xié)方差,同時(shí),下文中出現(xiàn)的下角標(biāo)(·)k/k-1均表示該參數(shù)(·)的一步預(yù)測(cè)值,下角標(biāo)(·)1:k-1表示該參數(shù)(·)的1~k-1時(shí)刻的值。

    由于測(cè)量過程中異常值的影響,量測(cè)噪聲呈現(xiàn)出非高斯特性,因此利用GSTM分布來刻畫似然函數(shù)。同時(shí)考慮到在目標(biāo)測(cè)量時(shí),飛行器的飛行空域和自身輻射特性的復(fù)雜性會(huì)導(dǎo)致量測(cè)信息統(tǒng)計(jì)特性未知,因此將量測(cè)噪聲建模成均值未知且時(shí)變的GSTM分布,根據(jù)量測(cè)噪聲得到似然函數(shù)建模公式如下:

    [N(xk;h(xk)+Ak,Βk/ξk)](1-yk)·

    p(δk)p(vk)dξkdyks.t.yk∈{0,1}

    (19)

    式中:與式(15)相同,隱變量ξk的積分為定積分,關(guān)于伯努利隨機(jī)變量yk的積分為不定積分,使得最終的積分結(jié)果為yk的函數(shù),用于積分外側(cè)的求和運(yùn)算。Ak和Βk分別表示量測(cè)噪聲均值和協(xié)方差,參考文獻(xiàn)[17],將均值和協(xié)方差的聯(lián)合函數(shù)建模成為高斯-逆Wishart分布,公式如下:

    p(Ak|Bk)=N(Ak;ak/k-1,bk/k-1Bk)

    (20)

    p(Bk|z1:k-1)=IW(Bk;uk/k-1,Uk/k-1)

    (21)

    式中:IW(·)表示逆Wishart分布,公式如下

    IW(B;λ,Ψ)=

    (22)

    根據(jù)2.1小節(jié)中的建模思路,將ST分布中的自由度參數(shù)vk表示為伽馬分布。

    p(vk)=G(vk;ck/k-1,dk/k-1)

    (23)

    混合概率δk滿足p(δk)=B(δk;e0,1-e0),其中,e0∈[0,1]表示PDFp(δk)的先驗(yàn)形狀參數(shù)。

    根據(jù)上文的建模,得到相對(duì)應(yīng)的概率圖模型,如圖1所示。

    圖1 概率圖模型Fig.1 Graphical model

    在Kalman濾波框架內(nèi),已知一步預(yù)測(cè)的建模,利用容積濾波器(Cubature Kalman filter, CKF)濾波方法得到時(shí)間更新方程:

    (24)

    (25)

    式中:Qk-1表示過程噪聲協(xié)方差;2m表示粒子數(shù);ζi和wi的計(jì)算公式如下

    Sk-1/k-1=chol(Pk-1/k-1)

    (26)

    (27)

    (28)

    式中:chol(·)表示Cholesky分解,[1]i表示點(diǎn)集[1]的第i組列向量,[1]=[Im,-Im],Im為m維的單位矩陣。

    在飛行器濾波跟蹤算法的實(shí)際應(yīng)用過程中,p(Bk|Bk-1),p(Ak|Ak-1),p(vk|vk-1)無法獲取,本文參考文獻(xiàn)[24],引入遺忘因子ρ來描述系統(tǒng)的動(dòng)態(tài)不確定性,ρ的取值范圍為[0.9,1)。

    uk/k-1=ρ(uk-1/k-1-md-1)+md+1

    (29)

    Uk/k-1=ρUk-1/k-1

    (30)

    ak/k-1=ak-1/k-1,bk/k-1=ρbk-1/k-1

    (31)

    ck/k-1=ρck-1/k-1,dk/k-1=ρdk-1/k-1

    (32)

    式中:md表示矩陣維度。

    2.4 后驗(yàn)分布的變分逼近

    根據(jù)2.3小節(jié)中對(duì)于一步預(yù)測(cè)和似然函數(shù)的建模,為了得到濾波算法的量測(cè)更新,需要對(duì)聯(lián)合后驗(yàn)PDFp(Θk|z1:k)進(jìn)行推導(dǎo),其中Θk包含一步預(yù)測(cè)和似然函數(shù)中的未知參數(shù),定義Θk為變分變量,Θk={xk,Ak,Βk,yk,δk,ξk,vk}。

    2.4.1變分貝葉斯

    由于聯(lián)合后驗(yàn)PDF中各參數(shù)的耦合,對(duì)于后驗(yàn)參數(shù)的估計(jì)難以得到解析解,因此,引入變分貝葉斯(Variational Bayesian, VB)算法[25],利用自由因子分布近似聯(lián)合分布,聯(lián)合后驗(yàn)PDFp(xk,Ak,Βk,yk,δk,ξk,vk|z1:k)可以寫作:

    p(Θk|z1:k)≈q(xk)q(Ak)q(Βk)q(yk)q(δk)·

    q(ξk)q(vk)

    (33)

    式中:q(·)表示p(·)的近似后驗(yàn)概率分布。

    為了得到后驗(yàn)PDF,需要最小化q(xk),q(Ak),q(Βk),q(yk),q(δk),q(ξk),q(vk)和聯(lián)合后驗(yàn)PDFp(Θk|z1:k)之間的Kullback-Leibler散度(Kullback-Leibler divergence, KLD)[26]。

    {q(xk),q(Ak),q(Βk),q(yk),q(δk),q(ξk),q(vk)}=

    argmin KLD(q(xk),q(Ak),q(Βk),q(yk),

    q(δk),q(ξk),q(vk)||p(Θk|z1:k))

    (34)

    式中:KLD用來衡量?jī)蓚€(gè)概率分布之間的差異程度[27],計(jì)算公式為:

    (35)

    根據(jù)式(35)可以看出,當(dāng)兩個(gè)概率分布之間的差異性越小,KLD越小,因此,可以通過最小化KLD,求取概率分布p(x)的近似分布q(x)。

    根據(jù)VB計(jì)算方法,式(34)的最優(yōu)解滿足下面的公式[24]:

    lnq(θ)=EΘ(-θ)[lnp(Θk,z1:k)]+cθ

    (36)

    式中:θ∈Θ并且滿足Θ(-θ)∪θ=Θ;cθ表示獨(dú)立于θ的常值。

    由于Θ中的各個(gè)參數(shù)之間相互耦合,無法直接得到其解析解,因此需要定點(diǎn)迭代求解式(36)。利用q(j)(Θ(-θ))的結(jié)果,近似求解后驗(yàn)PDFq(θ)在j+1次迭代時(shí)的PDFq(j+1)(θ),直到收斂到式(36)的局部最優(yōu)值[24]。下面在2.4.2小節(jié)中推導(dǎo)各個(gè)參數(shù)的迭代更新過程。

    2.4.2變分迭代過程

    已知聯(lián)合分布p(Θk,z1:k)可以被分解為

    p(Θk,z1:k)=p(z1:k-1)N(xk;xk/k-1,Pk/k-1)·

    [N(xk;h(xk)+Ak,Βk)]yk[N(xk;h(xk)+

    Ak,Βk/ξk)](1-yk)N(Ak;ak/k-1,bk/k-1Bk)·

    IW(Bk;uk/k-1,Uk/k-1)δkyk(1-δk)ykB(δk;e0,

    (37)

    1) 令θ={Ak,Bk}

    在第j+1次循環(huán)時(shí),{Ak,Bk}的聯(lián)合后驗(yàn)PDF更新為高斯-逆Wishart分布:

    (38)

    式中:下角標(biāo)(·)k/k均表示該參數(shù)(·)的估計(jì)值。

    (39)

    (40)

    式中:E(j)[x]表示參數(shù)x在j次迭代后的數(shù)學(xué)期望。

    (41)

    (42)

    (43)

    (44)

    式中:

    (45)

    Sk/k-1=chol(Pk/k-1)

    (46)

    (47)

    (48)

    2) 令θ={yk}

    在第j+1次循環(huán)時(shí),yk的后驗(yàn)PDF更新為Bernoulli分布:

    (49)

    (50)

    (51)

    根據(jù)式(49)~(51)和Bernoulli分布的性質(zhì),可以求解得到:

    (52)

    3) 令θ={δk}

    在第j+1次循環(huán)時(shí),δk的后驗(yàn)PDF更新為Beta分布

    (53)

    式中:

    (54)

    (55)

    根據(jù)式(53)~(55)和Beta分布的性質(zhì),可以求解得到:

    (56)

    (57)

    4) 令θ=(ξk)

    (58)

    式中:

    (59)

    (60)

    根據(jù)式(58)~(60)和伽馬分布的性質(zhì),可以求解得到:

    (61)

    (62)

    5) 令θ={vk}

    在第j+1次循環(huán)時(shí),vk后驗(yàn)PDF更新為

    (63)

    式中:

    (64)

    (65)

    6) 令θ={xk}

    在第j+1次循環(huán)時(shí),xk的后驗(yàn)PDF更新為

    (66)

    式中:

    (67)

    (68)

    (69)

    (70)

    式中:

    (71)

    (72)

    3 仿真校驗(yàn)

    3.1 高超聲速飛行器HTV-2的建模仿真

    根據(jù)小節(jié)1中對(duì)高超聲速飛行器HTV-2助推滑翔段的動(dòng)力學(xué)描述,以及考慮飛行器在實(shí)際飛行過程中的飛行走廊約束,完成對(duì)HTV-2的彈道規(guī)劃[28-29]。

    表1 HTV-2高超聲速飛行器彈道仿真參數(shù)Table 1 Trajectory simulation parameters of hypersonic vehicle HTV-2

    仿真結(jié)果圖彈道如圖2所示:

    圖2 彈道形態(tài)Fig.2 Ballistic configuration

    3.2 仿真場(chǎng)景設(shè)置

    采用地基雷達(dá)觀測(cè),觀測(cè)方程如下:

    (73)

    式中:Aak,Eek,Rrk分別為雷達(dá)探測(cè)的高低角方位角和斜距;υk為量測(cè)噪聲。

    (74)

    式中:w.p.表示“with probability”;p1,p2,p3表示異常值概率。R表達(dá)式如下:

    (75)

    式中:σA=σE=0.03°,σR=20 m。

    (76)

    根據(jù)第1小節(jié)對(duì)高超聲速飛行器的動(dòng)力學(xué)建模分析,在高超聲速飛行器的滑翔段,飛行器主要受重力和氣動(dòng)力的控制,因此可以從動(dòng)力學(xué)的角度出發(fā),設(shè)計(jì)合理的模型參數(shù)對(duì)動(dòng)力學(xué)進(jìn)行建模,利用擴(kuò)維后的動(dòng)力學(xué)模型進(jìn)行狀態(tài)估計(jì)[8,28,30-31]。

    本文參考文獻(xiàn)[28],采用彈道系數(shù)模型作為跟蹤模型。定義彈道系數(shù)為

    (77)

    (78)

    為了比較不同算法對(duì)目標(biāo)的狀態(tài)估計(jì)精度,給出狀態(tài)的均方根誤差(RMSE)和平均均方根誤差(ARMSE)的定義如下:

    RMSEk,X=

    (79)

    ARMSEX=

    (80)

    為了對(duì)比算法的估計(jì)性能,本文對(duì)比了5種濾波算法:傳統(tǒng)的CKF算法,量測(cè)噪聲采用高斯建模的自適應(yīng)變分貝葉斯容積卡爾曼濾波(Adaptive variational Bayesian cubature Kalman filtering, AVBCKF-R)算法[32],采用ST建模的魯棒變分貝葉斯學(xué)生t濾波(Robust variational Bayesian Student’s t filt-ering, RVBST)算法[33],未包含測(cè)量統(tǒng)計(jì)信息建模的GSTMF算法[6],以及本文中提出的改進(jìn)RGSTMF算法。同時(shí),為了更好地比較GSTM分布的建模優(yōu)勢(shì),將AVBCKF-R算法和RVBST算法中的噪聲均值按照本文提出的方法進(jìn)行建模。

    3.2.1仿真場(chǎng)景1

    在仿真場(chǎng)景1中,設(shè)置噪聲異常值變化幅度較小且變化平緩,設(shè)置異常值概率分別為:p1=0.01,p2=0.02,p3=0.04。在場(chǎng)景1中設(shè)置異常值所占的概率較小,噪聲重尾分布較輕。仿真結(jié)果見圖3~圖4和表2。

    表2 不同濾波算法的ARMSEsTable 2 The ARMSEs of different filtering algorithms

    圖3 位置變量的RMSEFig.3 The RMSEs of position

    圖4 速度變量的RMSEFig.4 The RMSEs of velocity

    從仿真圖和表格可以看出,由于缺少對(duì)噪聲的建模,CKF算法精度表現(xiàn)最差,本文提出的RGSTMF算法精度最優(yōu)。同時(shí)可以看出,與噪聲高斯建模的AVBCKF-R算法和ST建模的RVBST算法相比,本文提出的RGSTMF算法精度更高,這說明GSTM分布與高斯分布和ST分布相比,能更好地?cái)M合非平穩(wěn)噪聲,這與文獻(xiàn)[5]中的結(jié)論一致。此外,與文獻(xiàn)[5]中的GSTM算法相比,RGSTMF算法包含對(duì)噪聲統(tǒng)計(jì)特性的描述,因此在面對(duì)統(tǒng)計(jì)信息未知的量測(cè)噪聲時(shí),具有更好的魯棒性和更高的計(jì)算精度。

    3.2.2仿真場(chǎng)景2

    在仿真場(chǎng)景2中,設(shè)置噪聲變化幅度較大,設(shè)置異常值概率分別為:p1=0.1,p2=0.2,p3=0.4。在場(chǎng)景2中設(shè)置異常值所占的概率較大,噪聲重尾程度較重。仿真結(jié)果見圖5~圖6和表3。

    表3 不同濾波算法的ARMSETable 3 The ARMSEs of different filtering algorithms

    圖5 位置變量的RMSEFig.5 The RMSEs of position

    圖6 速度變量的RMSEFig.6 The RMSEs of velocity

    從仿真圖可以看出,在噪聲異常值進(jìn)一步增加的情況下,CKF算法已經(jīng)出現(xiàn)了較大的誤差,GSTM分布的誤差也進(jìn)一步加大,而文中所提出的RGSTMF算法仍然能保持較高的精度。

    綜合仿真1和仿真2,可以看出對(duì)比現(xiàn)有的濾波算法,本文所提出的RGSTMF算法濾波精度最高。與噪聲高斯建模的AVBCKF-R算法和ST建模的RVBST算法相比,所提出的RGSTMF算法無論是在面對(duì)異常值概率較低、變化較為緩慢的噪聲分布,還是異常值占比較高、噪聲突變的噪聲分布,都能夠表現(xiàn)出較為良好的跟蹤精度,這說明GSTM分布對(duì)非平穩(wěn)噪聲具有更好的擬合精度。同時(shí),可以看出,采用對(duì)噪聲均值建模的方法可以進(jìn)一步提升RGSTMF算法的魯棒性,使得其與CKF和GSTMF相比,能在噪聲統(tǒng)計(jì)信息未知的情況下保持較好的跟蹤精度。

    3.3 遺忘因子的影響

    本小節(jié)探討影響算法的可調(diào)節(jié)參數(shù)遺忘因子ρ對(duì)算法精度的影響,仿真結(jié)果見圖7和圖8。

    圖7 場(chǎng)景1中不同的遺忘因子對(duì)RMSE的影響Fig.7 RMSE performances of different forgetting factors in Scenario 1

    圖8 場(chǎng)景2中不同的遺忘因子對(duì)RMSE的影響Fig.8 RMSE performances of different forgetting factors in Scenario 2

    從仿真圖可以看出,當(dāng)遺忘因子選取為ρ=0.9,0.92,0.94,0.96,0.98時(shí),算法的精度沒有出現(xiàn)較大的差異,但是當(dāng)算法采用ρ=1時(shí),濾波誤差開始增大,并且隨著噪聲異常值概率的增加,誤差越明顯。這是因?yàn)楫?dāng)ρ=1,算法p(Bk|Bk-1),p(Ak|Ak-1),p(vk|vk-1)中的參數(shù)采用固定的建模參數(shù),難以應(yīng)對(duì)時(shí)變的噪聲。

    4 結(jié) 論

    針對(duì)高超聲速飛行器飛行過程中出現(xiàn)的統(tǒng)計(jì)特性未知的非平穩(wěn)非高斯探測(cè)噪聲,本文提出一種改進(jìn)RGSTMF算法。利用高斯-逆Wishart分布構(gòu)建時(shí)變且未知的均值模型用于代替原始GSTM中的零均值假設(shè),同時(shí)引入伽馬函數(shù)分布來刻畫不確定的非高斯程度,改進(jìn)的RGSTM分布使得濾波算法對(duì)非平穩(wěn)噪聲的建模精度得到進(jìn)一步的提高;將變分貝葉斯算法與RGSTM分布相結(jié)合,解析推導(dǎo)出了后驗(yàn)概率更新的顯式表達(dá)式。仿真結(jié)果驗(yàn)證了本文提出方案的可行性,與其他先進(jìn)算法相比,本文中所提出的RGSTMF算法能夠在先驗(yàn)信息未知的情況下,更好地刻畫噪聲時(shí)變且未知的均值和非高斯程度,在面對(duì)統(tǒng)計(jì)信息未知的非平穩(wěn)非高斯程度不確定的噪聲時(shí),表現(xiàn)出良好的魯棒性和跟蹤精度,并且隨著非高斯程度的進(jìn)一步增加,算法仍能保持良好的跟蹤精度。

    猜你喜歡
    后驗(yàn)超聲速先驗(yàn)
    高超聲速出版工程
    高超聲速飛行器
    基于對(duì)偶理論的橢圓變分不等式的后驗(yàn)誤差分析(英)
    基于無噪圖像塊先驗(yàn)的MRI低秩分解去噪算法研究
    貝葉斯統(tǒng)計(jì)中單參數(shù)后驗(yàn)分布的精確計(jì)算方法
    超聲速旅行
    基于自適應(yīng)塊組割先驗(yàn)的噪聲圖像超分辨率重建
    一種基于最大后驗(yàn)框架的聚類分析多基線干涉SAR高度重建算法
    基于平滑先驗(yàn)法的被動(dòng)聲信號(hào)趨勢(shì)項(xiàng)消除
    先驗(yàn)的廢話與功能的進(jìn)路
    好男人视频免费观看在线| 又粗又硬又长又爽又黄的视频| 九草在线视频观看| 成年女人在线观看亚洲视频 | 亚洲av不卡在线观看| 国产成人aa在线观看| 国产高清国产精品国产三级 | 天天躁夜夜躁狠狠久久av| 精品酒店卫生间| 丝瓜视频免费看黄片| 国内精品宾馆在线| 男女边吃奶边做爰视频| 男女下面进入的视频免费午夜| 亚洲国产日韩一区二区| 久热这里只有精品99| 亚洲人成网站高清观看| av线在线观看网站| 中文字幕免费在线视频6| 亚洲精品456在线播放app| 日本欧美国产在线视频| 美女cb高潮喷水在线观看| 王馨瑶露胸无遮挡在线观看| 国产精品精品国产色婷婷| 亚洲自拍偷在线| 国产精品国产三级专区第一集| 欧美激情久久久久久爽电影| 国产成人freesex在线| www.av在线官网国产| 永久免费av网站大全| 黄片无遮挡物在线观看| 日韩三级伦理在线观看| 免费不卡的大黄色大毛片视频在线观看| 七月丁香在线播放| 青青草视频在线视频观看| 日韩欧美精品免费久久| 久久人人爽人人片av| 国产亚洲一区二区精品| 黄色视频在线播放观看不卡| 少妇裸体淫交视频免费看高清| 精品人妻视频免费看| 亚洲成人av在线免费| 国产精品不卡视频一区二区| 欧美成人午夜免费资源| 欧美成人午夜免费资源| 九九爱精品视频在线观看| 国产欧美日韩一区二区三区在线 | 中国国产av一级| 一级片'在线观看视频| 一区二区三区精品91| 永久免费av网站大全| 国产亚洲5aaaaa淫片| 97人妻精品一区二区三区麻豆| av福利片在线观看| 全区人妻精品视频| 国产亚洲5aaaaa淫片| 五月伊人婷婷丁香| 亚洲av中文av极速乱| 91久久精品电影网| 在线观看av片永久免费下载| 亚洲欧美日韩无卡精品| 国产真实伦视频高清在线观看| 成人免费观看视频高清| 亚洲最大成人中文| 国产成人一区二区在线| 91狼人影院| 特级一级黄色大片| 特级一级黄色大片| 国产一区二区三区综合在线观看 | 蜜桃亚洲精品一区二区三区| 免费在线观看成人毛片| 51国产日韩欧美| 久久精品夜色国产| 秋霞在线观看毛片| 亚洲精品成人av观看孕妇| 麻豆乱淫一区二区| 三级国产精品片| 久久久久性生活片| 又爽又黄无遮挡网站| 狠狠精品人妻久久久久久综合| 久热久热在线精品观看| 亚洲不卡免费看| 精品人妻一区二区三区麻豆| 涩涩av久久男人的天堂| 国产免费又黄又爽又色| a级一级毛片免费在线观看| 美女被艹到高潮喷水动态| 午夜免费男女啪啪视频观看| 久久精品国产亚洲av涩爱| 国产男女超爽视频在线观看| 尾随美女入室| 国产色婷婷99| 亚洲图色成人| 99热这里只有精品一区| 国产一级毛片在线| 国产中年淑女户外野战色| 欧美日韩亚洲高清精品| av黄色大香蕉| 真实男女啪啪啪动态图| av专区在线播放| 午夜视频国产福利| 国产成人精品福利久久| 成人毛片60女人毛片免费| 国产午夜精品久久久久久一区二区三区| 九九在线视频观看精品| 欧美老熟妇乱子伦牲交| 日韩国内少妇激情av| 免费看光身美女| 国产午夜福利久久久久久| 日日啪夜夜撸| 99热网站在线观看| 在线观看美女被高潮喷水网站| 成人无遮挡网站| 久久精品久久久久久久性| 国产乱人偷精品视频| 久久久久久久亚洲中文字幕| 又黄又爽又刺激的免费视频.| 青春草国产在线视频| 丝袜脚勾引网站| 午夜福利视频1000在线观看| 亚洲婷婷狠狠爱综合网| 免费观看无遮挡的男女| 国产成人免费无遮挡视频| 男人和女人高潮做爰伦理| 香蕉精品网在线| 久久精品国产a三级三级三级| 蜜桃久久精品国产亚洲av| 蜜桃亚洲精品一区二区三区| 伊人久久国产一区二区| 尾随美女入室| 黄色配什么色好看| 欧美变态另类bdsm刘玥| 日韩成人伦理影院| 岛国毛片在线播放| 色网站视频免费| 日韩av在线免费看完整版不卡| 亚洲欧美精品专区久久| 国产男女内射视频| 高清欧美精品videossex| 亚洲不卡免费看| 日韩成人伦理影院| 天堂中文最新版在线下载 | 一级av片app| 免费黄色在线免费观看| 国产成人精品婷婷| 国产 一区精品| 日本免费在线观看一区| 亚洲av日韩在线播放| 午夜福利高清视频| 精华霜和精华液先用哪个| 亚洲精品色激情综合| 搡老乐熟女国产| 午夜免费鲁丝| 一个人观看的视频www高清免费观看| 欧美日韩精品成人综合77777| 交换朋友夫妻互换小说| 国产亚洲5aaaaa淫片| 国产一区亚洲一区在线观看| 91久久精品国产一区二区三区| 国产黄a三级三级三级人| 一本色道久久久久久精品综合| 97超视频在线观看视频| 国产成年人精品一区二区| 听说在线观看完整版免费高清| 一区二区三区精品91| 精品国产乱码久久久久久小说| 精品久久久久久电影网| 国产成人福利小说| 精品一区在线观看国产| 亚洲,一卡二卡三卡| 国产精品一二三区在线看| 别揉我奶头 嗯啊视频| 日韩电影二区| 高清视频免费观看一区二区| 国产有黄有色有爽视频| 2022亚洲国产成人精品| 久久久久久久久大av| 成人特级av手机在线观看| 性色avwww在线观看| 欧美3d第一页| 一级毛片aaaaaa免费看小| 看免费成人av毛片| 1000部很黄的大片| 久久精品国产自在天天线| 国产男女内射视频| 亚洲精品久久午夜乱码| 蜜桃亚洲精品一区二区三区| 伦精品一区二区三区| 九九在线视频观看精品| av免费在线看不卡| 少妇猛男粗大的猛烈进出视频 | 成人高潮视频无遮挡免费网站| 亚洲欧美精品自产自拍| 99久久精品热视频| 啦啦啦中文免费视频观看日本| 亚洲成人久久爱视频| 免费人成在线观看视频色| 成人毛片60女人毛片免费| 青春草亚洲视频在线观看| 亚洲av不卡在线观看| 国内精品美女久久久久久| 午夜免费鲁丝| 国产在线男女| 免费播放大片免费观看视频在线观看| 人妻少妇偷人精品九色| 国产黄片美女视频| 激情 狠狠 欧美| 久久精品国产鲁丝片午夜精品| 天堂中文最新版在线下载 | 真实男女啪啪啪动态图| 五月伊人婷婷丁香| 人妻系列 视频| av天堂中文字幕网| 91在线精品国自产拍蜜月| 欧美97在线视频| 婷婷色综合大香蕉| 日韩一区二区三区影片| 我的老师免费观看完整版| 免费大片18禁| 99久久精品一区二区三区| 在线观看一区二区三区| 三级国产精品片| 全区人妻精品视频| 毛片一级片免费看久久久久| 中国美白少妇内射xxxbb| 亚洲一级一片aⅴ在线观看| 成年免费大片在线观看| 亚洲国产高清在线一区二区三| 夫妻午夜视频| 国产成人精品一,二区| 日本午夜av视频| 热99国产精品久久久久久7| 亚洲国产日韩一区二区| 青春草亚洲视频在线观看| 最近最新中文字幕大全电影3| 肉色欧美久久久久久久蜜桃 | 天堂俺去俺来也www色官网| 亚洲婷婷狠狠爱综合网| 超碰97精品在线观看| 搞女人的毛片| 国产高潮美女av| av免费在线看不卡| 啦啦啦中文免费视频观看日本| 国产成人aa在线观看| 成人午夜精彩视频在线观看| 国产高潮美女av| 有码 亚洲区| 国产精品国产三级国产专区5o| 久久女婷五月综合色啪小说 | 听说在线观看完整版免费高清| 人妻系列 视频| 麻豆成人av视频| 男女下面进入的视频免费午夜| 久久韩国三级中文字幕| 免费黄频网站在线观看国产| 国产免费福利视频在线观看| 五月开心婷婷网| www.av在线官网国产| 日韩欧美精品v在线| 99九九线精品视频在线观看视频| 国产精品熟女久久久久浪| 日韩视频在线欧美| 亚洲精品日韩在线中文字幕| 激情 狠狠 欧美| 男人狂女人下面高潮的视频| 大香蕉久久网| 成人国产av品久久久| 亚洲成人一二三区av| 午夜免费观看性视频| 在线观看av片永久免费下载| a级毛片免费高清观看在线播放| 免费观看的影片在线观看| 男人添女人高潮全过程视频| 国产乱人视频| 搡老乐熟女国产| 丝袜喷水一区| 婷婷色综合大香蕉| 男人狂女人下面高潮的视频| 亚洲真实伦在线观看| 国产老妇伦熟女老妇高清| 白带黄色成豆腐渣| 日本猛色少妇xxxxx猛交久久| 最近手机中文字幕大全| 中国美白少妇内射xxxbb| 亚洲电影在线观看av| av黄色大香蕉| 欧美三级亚洲精品| 色吧在线观看| 久久精品夜色国产| 欧美+日韩+精品| 亚洲精品日韩av片在线观看| av国产精品久久久久影院| 国产爽快片一区二区三区| 精品国产一区二区三区久久久樱花 | 波野结衣二区三区在线| 午夜亚洲福利在线播放| 日韩欧美一区视频在线观看 | 午夜福利在线在线| 日韩 亚洲 欧美在线| 亚洲第一区二区三区不卡| 91久久精品国产一区二区成人| 久久久欧美国产精品| 午夜激情久久久久久久| 免费黄频网站在线观看国产| 色吧在线观看| 乱码一卡2卡4卡精品| 亚洲欧美精品专区久久| 国产亚洲一区二区精品| 国产精品偷伦视频观看了| 最近2019中文字幕mv第一页| 国产永久视频网站| 国产精品99久久99久久久不卡 | 国产亚洲一区二区精品| 自拍偷自拍亚洲精品老妇| 简卡轻食公司| 天天躁夜夜躁狠狠久久av| 国产精品熟女久久久久浪| 色网站视频免费| 一个人观看的视频www高清免费观看| 免费看不卡的av| 国产成人精品婷婷| 久久精品熟女亚洲av麻豆精品| 99久久精品热视频| 亚洲av在线观看美女高潮| 亚洲成人中文字幕在线播放| 晚上一个人看的免费电影| 国产精品偷伦视频观看了| 亚洲精品成人av观看孕妇| 国产永久视频网站| 免费黄色在线免费观看| 国产精品国产三级专区第一集| 国产精品熟女久久久久浪| 国产av码专区亚洲av| 一区二区三区免费毛片| 大陆偷拍与自拍| 国产精品久久久久久av不卡| 成年版毛片免费区| 亚洲高清免费不卡视频| 最近中文字幕2019免费版| 久久韩国三级中文字幕| 老司机影院毛片| 欧美成人a在线观看| 大香蕉久久网| av卡一久久| 白带黄色成豆腐渣| 18禁在线无遮挡免费观看视频| 99re6热这里在线精品视频| 在线观看国产h片| 99久久精品一区二区三区| 欧美精品国产亚洲| 联通29元200g的流量卡| 国产男女内射视频| 看黄色毛片网站| av一本久久久久| 小蜜桃在线观看免费完整版高清| 夫妻午夜视频| 亚洲精华国产精华液的使用体验| 国产一区亚洲一区在线观看| 天堂网av新在线| 日韩三级伦理在线观看| 青春草视频在线免费观看| 免费不卡的大黄色大毛片视频在线观看| 蜜桃亚洲精品一区二区三区| 在线观看一区二区三区| 亚洲国产精品国产精品| 久久久久久久久大av| av黄色大香蕉| 国产一区亚洲一区在线观看| 国产欧美日韩精品一区二区| 91久久精品国产一区二区三区| 成人一区二区视频在线观看| 亚洲欧洲日产国产| 国产精品不卡视频一区二区| 黄色怎么调成土黄色| 99re6热这里在线精品视频| 综合色丁香网| 国产精品.久久久| 少妇猛男粗大的猛烈进出视频 | 毛片女人毛片| 身体一侧抽搐| 少妇 在线观看| 亚洲精品日韩在线中文字幕| 免费看av在线观看网站| 26uuu在线亚洲综合色| 女的被弄到高潮叫床怎么办| 日本wwww免费看| 亚洲精品一区蜜桃| 国产一区亚洲一区在线观看| 人妻 亚洲 视频| 亚洲精品国产色婷婷电影| 亚洲欧美成人精品一区二区| 久久99热这里只频精品6学生| 亚洲性久久影院| 日韩成人av中文字幕在线观看| 国产伦理片在线播放av一区| 观看美女的网站| 久久鲁丝午夜福利片| 亚洲色图综合在线观看| 亚洲成人av在线免费| 春色校园在线视频观看| 22中文网久久字幕| 精品人妻熟女av久视频| 黄片无遮挡物在线观看| 欧美高清性xxxxhd video| 亚洲国产精品专区欧美| 国产色婷婷99| 亚洲最大成人av| av在线亚洲专区| 一级毛片久久久久久久久女| 久久久久久久久久久免费av| 国语对白做爰xxxⅹ性视频网站| 最近中文字幕2019免费版| 国产亚洲av嫩草精品影院| 欧美三级亚洲精品| 久久精品夜色国产| 久久国产乱子免费精品| 伊人久久国产一区二区| 久久久国产一区二区| 少妇的逼水好多| 免费大片18禁| 黄色日韩在线| 热99国产精品久久久久久7| 91久久精品电影网| 日本黄色片子视频| 亚洲美女视频黄频| 成人美女网站在线观看视频| 国产精品爽爽va在线观看网站| 18+在线观看网站| 色视频在线一区二区三区| 精品酒店卫生间| 午夜老司机福利剧场| 亚洲欧美精品自产自拍| 久久人人爽人人片av| 亚洲熟女精品中文字幕| 尤物成人国产欧美一区二区三区| 精品一区二区免费观看| 成人美女网站在线观看视频| 噜噜噜噜噜久久久久久91| 麻豆成人av视频| 色播亚洲综合网| 国产精品久久久久久久电影| 国产人妻一区二区三区在| 亚洲真实伦在线观看| 中文字幕亚洲精品专区| 亚洲av免费在线观看| 人妻 亚洲 视频| 一级毛片久久久久久久久女| av天堂中文字幕网| 麻豆成人午夜福利视频| 亚洲成人精品中文字幕电影| 久久精品国产鲁丝片午夜精品| 七月丁香在线播放| 80岁老熟妇乱子伦牲交| 久久精品久久精品一区二区三区| 岛国毛片在线播放| 亚洲av免费高清在线观看| 22中文网久久字幕| 香蕉精品网在线| 久久久久久久久大av| 免费人成在线观看视频色| 亚洲av电影在线观看一区二区三区 | 别揉我奶头 嗯啊视频| 国产伦精品一区二区三区视频9| 直男gayav资源| 777米奇影视久久| 免费看不卡的av| 久久6这里有精品| 日本免费在线观看一区| 成年人午夜在线观看视频| 汤姆久久久久久久影院中文字幕| 蜜桃亚洲精品一区二区三区| 亚洲一级一片aⅴ在线观看| 人妻 亚洲 视频| 嫩草影院新地址| 九九在线视频观看精品| 亚洲成人久久爱视频| 成人黄色视频免费在线看| av在线app专区| 男女边吃奶边做爰视频| 亚洲国产日韩一区二区| 亚洲四区av| 亚洲av国产av综合av卡| 亚洲不卡免费看| 大片电影免费在线观看免费| 不卡视频在线观看欧美| 毛片女人毛片| 黄片无遮挡物在线观看| 最近最新中文字幕大全电影3| 大片电影免费在线观看免费| 内射极品少妇av片p| 97热精品久久久久久| av网站免费在线观看视频| 久久99热6这里只有精品| 亚洲国产av新网站| 亚洲av中文字字幕乱码综合| 久久6这里有精品| 制服丝袜香蕉在线| 亚洲成人一二三区av| 少妇人妻 视频| 国产精品麻豆人妻色哟哟久久| 男女边吃奶边做爰视频| 99久久九九国产精品国产免费| 亚洲国产最新在线播放| 男人舔奶头视频| 成年av动漫网址| 一级片'在线观看视频| 日韩电影二区| 97在线视频观看| 国产成人一区二区在线| 国产在视频线精品| 成人毛片a级毛片在线播放| 久久久久精品性色| 亚洲精品亚洲一区二区| 直男gayav资源| 欧美另类一区| .国产精品久久| 久久精品夜色国产| www.色视频.com| 亚洲人成网站在线观看播放| 色播亚洲综合网| av网站免费在线观看视频| 美女视频免费永久观看网站| 亚洲av福利一区| 久久精品国产亚洲av涩爱| 毛片女人毛片| 国产 一区精品| 成人美女网站在线观看视频| 亚洲精品久久午夜乱码| 人妻夜夜爽99麻豆av| 香蕉精品网在线| 啦啦啦在线观看免费高清www| 可以在线观看毛片的网站| 女人十人毛片免费观看3o分钟| 夜夜看夜夜爽夜夜摸| 成年免费大片在线观看| 爱豆传媒免费全集在线观看| 精品少妇黑人巨大在线播放| 99久久精品国产国产毛片| 亚洲成人精品中文字幕电影| 国产成人精品一,二区| 久久久久久久精品精品| 美女脱内裤让男人舔精品视频| 亚洲人与动物交配视频| 亚洲国产av新网站| 亚洲精华国产精华液的使用体验| 国产成人a∨麻豆精品| 日韩av不卡免费在线播放| 中文字幕免费在线视频6| 一级a做视频免费观看| 色综合色国产| 久久久久国产精品人妻一区二区| 日韩av不卡免费在线播放| 久久精品国产亚洲av天美| 国产精品爽爽va在线观看网站| 成人免费观看视频高清| 国产中年淑女户外野战色| 久久久久久久久大av| 国产毛片在线视频| 国产精品嫩草影院av在线观看| 国产极品天堂在线| 国产日韩欧美亚洲二区| 少妇裸体淫交视频免费看高清| 九九在线视频观看精品| 欧美3d第一页| 又粗又硬又长又爽又黄的视频| 精品少妇黑人巨大在线播放| 国产精品久久久久久精品电影小说 | 亚洲欧美成人综合另类久久久| 国产淫语在线视频| 六月丁香七月| 免费看a级黄色片| 国产一区有黄有色的免费视频| 国产精品爽爽va在线观看网站| 国产黄片美女视频| 秋霞在线观看毛片| 国产成年人精品一区二区| 久久久久国产网址| 久久久久国产精品人妻一区二区| 日本wwww免费看| 亚洲无线观看免费| 日韩亚洲欧美综合| 亚洲av成人精品一二三区| 日韩成人伦理影院| 久久久久网色| 精品久久久久久电影网| 国产成人一区二区在线| 国产淫片久久久久久久久| 亚洲国产av新网站| 欧美+日韩+精品| 国产精品一区www在线观看| 美女国产视频在线观看| 内射极品少妇av片p| 简卡轻食公司| 亚洲av男天堂| 一个人看视频在线观看www免费| 久久久色成人| 搡女人真爽免费视频火全软件| 欧美变态另类bdsm刘玥| 男插女下体视频免费在线播放| 国产午夜精品久久久久久一区二区三区| 国产大屁股一区二区在线视频| 久久ye,这里只有精品| 午夜福利在线观看免费完整高清在| 又爽又黄无遮挡网站| 精品久久久久久久久av| 免费av不卡在线播放| 97精品久久久久久久久久精品| 少妇人妻一区二区三区视频| 极品少妇高潮喷水抽搐| 精品一区在线观看国产| 亚洲色图av天堂| 国产 精品1| 美女主播在线视频| 中文字幕久久专区| 神马国产精品三级电影在线观看| 午夜福利视频精品|