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

    基于 SNVA 的機動目標狀態(tài)估計

    2016-01-06 05:30:12喬雪,白亞騰
    彈道學報 2015年1期

    基于SNVA的機動目標狀態(tài)估計

    喬雪,白亞騰

    (中國礦業(yè)大學 信息與電氣工程學院,江蘇 徐州,221116)

    摘要:利用位置預測估計值與位置濾波估計值之間的偏差進行加速度方差自適應調(diào)節(jié),提出一種基于狀態(tài)噪聲方差自適應(SNVA)的機動目標狀態(tài)估計方法。采用SNVA對目標加速度噪聲方差進行自適應調(diào)整,實現(xiàn)了對當前統(tǒng)計模型的改進;利用擴展卡爾曼濾波算法對目標狀態(tài)進行估計。仿真結(jié)果表明,基于SNVA的擴展卡爾曼濾波算法對機動目標速度估計的絕對誤差小于0.1 m/s,加速度估計的絕對誤差小于0.1 m/s2,能夠?qū)C動目標的狀態(tài)進行準確的估計。

    關(guān)鍵詞:機動目標;當前統(tǒng)計模型;擴展卡爾曼濾波;狀態(tài)估計

    收稿日期:2014-02-24

    基金項目:教育部博士學科點專項科研

    作者簡介:喬雪(1990- ),女,碩士研究生,研究方向為復雜系統(tǒng)分析與控制。E-mail:1219771581@qq.com。

    中圖分類號:TJ765.4文獻標識碼:A

    ManeuveringTargetStateEstimationBasedonSNVA

    QIAOXue,BAIYa-teng

    (SchoolofInformationandElectricalEngineering,ChinaUniversityofMining&Technology,Xuzhou221116,China)

    Abstract:By using the difference between the location forecast estimation and the corrected location estimation,the acceleration variance was adjusted adaptively.A maneuvering target state estimation method based on state noise variance adaptive(SNVA)was proposed.The SNVA was used to adaptively adjust the system noise variance in target state estimation system,which can improve the current statistical model.The extended Kalman filter algorithm was used to estimate the target state.Simulation results show that the extended Kalman filter algorithm based on SNVA can estimate the target state accurately.The absolute estimation error of velocity is less than 0.1 m/s,and the absolute estimation error of acceleration is less than 0.1 m/s2.The target state can be accurately estimated by the proposed method.

    Keywords:maneueringtarget;currentstatisticalmodel;extendedKalmanfilter;stateestimation

    對于機動目標狀態(tài)估計問題的研究是目前自動控制領(lǐng)域的熱點課題之一。將濾波算法與機動目標模型[1-2]相結(jié)合可以對機動目標的狀態(tài)進行估計。在濾波算法的發(fā)展歷史中有很多學者提出了不同的濾波算法,其中,卡爾曼濾波理論是由Kalman于1960年提出的。對于線性系統(tǒng),運用卡爾曼濾波算法可以得到系統(tǒng)狀態(tài)的最小均方差估計。對發(fā)生機動運動的非線性系統(tǒng)通常采用擴展卡爾曼濾波(extendedKalmanfilter,EKF)算法[4-5]。EKF算法首先將非線性狀態(tài)方程通過泰勒級數(shù)展開進行線性化處理,線性化后所得到的濾波算法通常被用于對機動目標的狀態(tài)進行估計。

    在對機動目標的狀態(tài)進行估計之前,首先應該建立機動目標的模型。周宏仁提出的“當前”統(tǒng)計(currentstatistical,CS)模型本質(zhì)上是一種時間相關(guān)模型,機動目標“當前”加速度的預測值即為均值?!爱斍啊苯y(tǒng)計模型是目前最常用、效果較好的一類實用模型,但是該模型也存在一些局限性,主要表現(xiàn)在以下2個方面:①濾波估計的快速性和準確性均會受限于預先給出的加速度預測極限值;②如果目標實際機動的狀態(tài)與估計狀態(tài)不符,則會造成濾波發(fā)散,對機動目標的狀態(tài)估計結(jié)果造成影響。

    因此,針對上述局限性①的問題,本文在基于CS模型和EKF算法的機動目標狀態(tài)估計基礎(chǔ)上,提出一種基于狀態(tài)噪聲方差自適應(statenoisevarianceadaptive,SNVA)的機動目標狀態(tài)估計方法。該方法的基本思想是:通過濾波對機動目標的位置進行估計,并求得位置濾波估計值與預測值之間的偏差,進而可以計算偏差在采樣時間內(nèi)的變化率,用所得的位置偏差變化率對當前時刻的加速度方差進行修正。采用這種方法是為了避免由于對加速度的極限值設(shè)定不當而造成估計精度降低的問題。與仿真結(jié)果的對比證明了這種方法在對強機動目標的狀態(tài)進行估計時比傳統(tǒng)的EKF濾波算法更加有效。

    1基于“當前”統(tǒng)計模型的機動目標狀態(tài)估計

    假設(shè)機動目標系統(tǒng)的離散狀態(tài)空間模型為

    X(k)=Φ[X(k-1),k-1]+

    Γ[X(k-1),k-1]W(k-1)

    (1)

    Y(k)=H[X(k),k]+V(k)

    (2)

    式中:X(k)是n維狀態(tài)向量;Y(k)是m維觀測向量;W(k-1)和V(k)分別是k-1時刻和k時刻的系統(tǒng)噪聲向量;Γ[X(k-1),k-1]是n×p維系統(tǒng)噪聲輸入矩陣;H[X(k),k]是m×n維觀測矩陣。

    ①系統(tǒng)狀態(tài)一步預測估計方程。

    (3)

    (4)

    Γ(k,k-1)=Γ[X(k-1/k-1),k-1]

    (5)

    (6)

    ②一步預測估計誤差方程。

    P(k|k-1)=Φ(k,k-1)P(k-1/k-1)×

    ΦT(k,k-1)+Γ(k,k-1)Q(k)ΓT(k,k-1)

    (7)

    式中:P(k-1/k-1)是后驗估計誤差協(xié)方差矩陣,Q(k)為系統(tǒng)狀態(tài)噪聲的方差陣,計算公式為

    (8)

    (9)

    式中:amax,a-max分別為預先設(shè)定的加速度上、下限。

    至此,根據(jù)式(1)和式(2)以及擴展卡爾曼濾波方程,就可以得到目標的狀態(tài)估計。圖1給出了擴展卡爾曼濾波算法的結(jié)構(gòu)框圖,由于系統(tǒng)的量測方程在直角坐標下是非線性的,則應將量測方程中的觀測矩陣進行線性化處理:

    (10)

    圖1中,Q(k)為系統(tǒng)狀態(tài)噪聲W(k)的方差,R(k)為觀測噪聲V(k)的方差,K(k)為卡爾曼濾波增益。

    圖1 擴展卡爾曼濾波算法的結(jié)構(gòu)框圖

    在利用擴展卡爾曼濾波進行狀態(tài)估計的過程中,一方面要求預先給定機動目標加速度的極限值amax和a-max,這個參數(shù)不能自適應調(diào)節(jié),因此這種狀態(tài)估計算法的估計能力有限;另一方面,如果缺乏必要的先驗信息,機動目標加速度的極限值很難確定,往往與實際系統(tǒng)模型不匹配,因此造成狀態(tài)估計的精度大大降低。

    2狀態(tài)噪聲方差的自適應性

    由式(6)可知,在CS模型中狀態(tài)噪聲方差取決于機動目標加速度的極限值amax和a-max。當加速度極限值的取值較大時,較大的噪聲方差能夠保證對目標機動的快速響應,但是會造成估計精度的降低;當加速度極限值的取值較小時,狀態(tài)估計精度有所提高,但是系統(tǒng)對機動目標的響應速度較慢。因此,系統(tǒng)不能實時反映目標的真實機動情況。為了擴大基于CS模型的狀態(tài)估計算法在機動目標狀態(tài)估計中的應用范圍,提高濾波精度,將上述狀態(tài)估計算法進行改進,提出基于狀態(tài)噪聲方差自適應(state noise variance adaptive,SNVA)的狀態(tài)估計算法。該算法利用位置觀測值和位置濾波估計值之間偏差的變化率對狀態(tài)噪聲方差進行實時調(diào)整,從而避免加速度上、下限對濾波估計的影響。

    在第k次采樣周期T(T是指相鄰2次實測之間的時間間隔)內(nèi)有:

    (11)

    (12)

    (13)

    (14)

    (15)

    (16)

    式中:μ是一個正常數(shù)。聯(lián)立式(11)~式(16),可得加速度狀態(tài)噪聲方差:

    (17)

    圖2 基于SNVA的機動目標狀態(tài)估計結(jié)構(gòu)圖

    3仿真研究

    3.1 目標機動狀態(tài)估計

    實驗主要考察導彈在鉛垂面攻擊目標的情況,vM,vT分別為導彈和目標的瞬時速率;θ為彈道傾角;θT為目標的高低角;xM,yM和xT,yT分別表示導彈和目標在x和y方向上的位移;比例系數(shù)K的取值在10到15之間。導彈與目標在地面坐標系中的位置關(guān)系如圖3所示,下標d表示地面坐標系。

    圖3 導彈與目標在地面坐標系中的關(guān)系

    目標的運動模型為

    (18)

    (19)

    比例導引時導彈與目標相對運動的方程組為

    (20)

    (21)

    式中:r為導彈到目標的相對距離,q是目標線方位角,假設(shè)地面水平,則θT≈0。設(shè)目標是機動的,且存在未知的擾動量ΔθT,ΔvT;將其代入式(18)和式(19),忽略由于ΔθT,ΔvT的存在而引起的Δr,可得到擾動方程組:

    (22)

    令w1=ΔθT,w2=ΔvT,整理導彈與目標運動方程組,得到系統(tǒng)的狀態(tài)方程和輸出方程:

    (23)

    (24)

    式中:

    其中:g為重力加速度,?為俯仰角,FP為發(fā)動機推力,Fx為迎面阻力,Fy為升力,m為導彈的質(zhì)量,α為攻角,Mz為俯仰力矩,Jz為導彈繞Ozd軸轉(zhuǎn)動的轉(zhuǎn)動慣量,ωz為導彈繞Ozd軸的旋轉(zhuǎn)角速度。

    根據(jù)所設(shè)計的模型采用Monte Carlo仿真方法進行計算機仿真,取采樣次數(shù)N=100,得到機動目標沿x軸和y軸方向的速度估計誤差均值曲線,如圖4和圖5所示,圖中,Δv表示速度誤差,在實驗中采樣周期T=2 s。

    圖4 速度誤差曲線

    由圖4可以看出,在100次的采樣結(jié)果中,x軸方向和y軸方向的速度濾波估計值與真實值之間的誤差均小于0.1 m/s。因此可以說明,經(jīng)過SNVA自適應處理后的機動狀態(tài)估計方法能夠較為準確地反映目標的真實運動狀態(tài)。

    在對目標濾波加速度進行的估計中,通過計算機選取400次的采樣結(jié)果,結(jié)合目標的真實運動加速度求出加速度估計均方根誤差σa,并將結(jié)果擬合成曲線,如圖5所示。σa的具體計算公式如下[10]:

    (25)

    圖5 加速度均方根誤差曲線

    由圖5可以看出,在對目標x方向以及y方向上加速度進行的估計中,濾波后的加速度估計值與目標加速度真實值之間的均方根誤差穩(wěn)定在一個比較小的范圍內(nèi),說明估計值與真實值之間的偏差較小。因此可以說明,基于SNVA的狀態(tài)估計算法對機動目標的加速度能夠進行比較準確的估計。

    綜合上述的仿真結(jié)果可知,在對2個方向上機動目標運動狀態(tài)進行的估計中,濾波估計誤差都在誤差允許的范圍內(nèi),并沒有出現(xiàn)濾波發(fā)散現(xiàn)象。由此可以得出以下結(jié)論:基于SNVA的機動目標狀態(tài)估計能夠得到目標機動狀態(tài)的近似無偏估計,并且估計精度高。

    3.2 對比仿真

    為了進一步驗證基于SNVA的狀態(tài)估計算法在機動目標狀態(tài)估計中的優(yōu)勢,將SNVA算法和擴展卡爾曼濾波算法同時應用于3.1節(jié)的仿真實驗中進行對比。在利用擴展卡爾曼濾波對目標狀態(tài)進行估計時預先設(shè)置機動頻率為0.001 Hz,加速度的極限值amax=1.0 km/s2和a-max=-0.5 km/s2。實驗中得到目標的加速度估計以及加速度估計誤差曲線,如圖6和圖7所示。圖中ax和ay分別表示目標在x軸和y軸方向上的加速度瞬時值,N表示采樣次數(shù),采樣周期T=2 s。

    圖6 加速度估計曲線

    由圖6可以看出,在對目標x軸方向和y軸方向上的加速度進行估計的過程中,每次采樣的基于SNVA的加速度估計值都比擴展卡爾曼濾波的加速度估計值更加接近真實值。因此可以說明,在機動目標的狀態(tài)估計問題中,基于SNVA的狀態(tài)估計算法與擴展卡爾曼濾波算法相比能更加準確快速地對目標狀態(tài)進行估計。

    目標的加速度誤差估計曲線如圖7所示,圖中Δax和Δay分別表示目標在x軸和y軸方向上的加速度誤差估計,采樣周期T=2 s。

    圖7 加速度估計誤差曲線

    由圖7可以看出,在對目標的加速度進行估計的過程中,每個采樣周期內(nèi)基于SNVA算法的加速度估計誤差都明顯小于擴展卡爾曼濾波的估計誤差,誤差均值在0附近波動,幅值不大,并沒有發(fā)生濾波發(fā)散。由此說明,在機動目標的狀態(tài)估計問題中,基于SNVA的狀態(tài)估計算法比擴展卡爾曼濾波更能準確地對目標的狀態(tài)進行估計。

    4結(jié)論

    本文提出的SNVA算法避免了在機動目標狀態(tài)估計問題中預先設(shè)定的加速度極限值對于狀態(tài)估計的影響。該方法的特點主要體現(xiàn)在以下2個方面:

    ①當目標做機動性較強的機動運動時,SNVA可以對目標的狀態(tài)進行近似無偏估計;

    ②避免了預先對加速度極限值的設(shè)定,能夠?qū)崟r調(diào)整加速度極限范圍。

    因此,在實際工程對跟蹤目標的狀態(tài)估計的應用中,SNVA的實用性較強,估計精度高,跟蹤效果很好,并且大大拓寬了對機動目標的機動跟蹤范圍。

    參考文獻

    [1]蘭義華,任浩征,張勇,等.一種基于“當前”模型的改進卡爾曼濾波算法.山東大學學報,2012,42(5):12-17.

    LANYi-hua,RENHao-zheng,ZHANGYong,etal.AnimprovedKalmanfilteralgorithmbasedonthe“current”model.JournalofShandongUniversity,2012,42(5):12-17.(inChinese)

    [2]LIUJian-shu,LIRen-hou,LIUYun-long.Interactingmultiplemodelalgorithmbasedonthecurrentstatisticalmodel.SystemsEngineeringandElectronics,2008,30(7):1 351-1 354.

    [3]WANGJin-hua,CAOJie,WUDi,etal.Anobjecttrackingalgorithmbasedonthe"current"statisticalmodelandthemulti-featurefusion.JournalofSoftware,2012,7(9):2 000-2 008.

    [4]常思江,王中原,牛春峰.基于卡爾曼濾波的彈箭飛行狀態(tài)估計方法.彈道學報,2010,22(3):94-98.

    CHANGSi-jiang,WANGZhong-yuan,NIUChun-feng.EstimationapproachofflightstateforguidedprojectilesbasedonKalmanfilter.JournalofBallistics,2010,22(3):94-98.(inChinese)

    [5]史金光,劉猛,曹成壯,等.彈道修正落點預報方法研究.彈道學報,2014,26(2):29-33.

    SHIJin-guang,LIUMeng,CAOCheng-zhuang,etal.Researchonmethodoffallingpointpredictionfortrajectorycorrectionprojectile.JournalofBallistics,2014,26(2):29-33.(inChinese)

    [6]周宏仁,敬忠良.機動目標跟蹤.北京:國防工業(yè)出版社,1991.

    ZHOUHong-ren,JINGZhong-liang.Trackingofmaneuveringtargets.Beijing:NationalDefenceIndustryPress,1991.(inChinese)

    [7]王以標,徐毓.基于當前統(tǒng)計模型的改進自適應濾波算法.火力與指揮控制,2010,35(5):86-89.

    WANGYi-biao,XUYu.Studyonimprovedadaptivefilteringalgorithmbasedoncurrentstatisticalmodel.FireControl&CommandControl,2010,35(5):86-89.(inChinese)

    [8]徐琰珂,梁曉庚,賈曉洪.基于ANFIS的自適應機動目標狀態(tài)估計算法.系統(tǒng)工程與電子技術(shù),2013,35(2):250-255.

    XUYan-ke,LIANGXiao-geng,JIAXiao-hong.AdaptivemaneuveringtargetstateestimationalgorithmbasedonANFIS.SystemsEngineeringandElectronics,2013,35(2):250-255.(inChinese)

    [9 ]余躍,劉藻珍.擴展卡爾曼濾波器在導彈制導系統(tǒng)中的應用與仿真//系統(tǒng)仿真技術(shù)及其應用學術(shù)交流會,合肥:中國科學技術(shù)大學出版社,2003:397-400.

    YUYue,LIUZao-zhen.SimulationofhomingmissilesystemusingextendedKalmanfilter//AcademicCommunicationSystemSimulationTechnologyanditsApplication.Hefei:PressofUniversityofScienceandTechnologyofChina,2003:397-400.(inChinese)

    [10] 楊宏,李亞安,李國輝.一種改進擴展卡爾曼濾波新方法.計算機工程與應用,2010,46(19):18-20.

    YANGHong,LIYa-an,LIGuo-hui.AnewmethodofimprovedextendedKalmanfilter.ComputerEngineeringandApplications,2010,46(19):18-20.(inChinese)

    国产精品人妻久久久影院| 亚洲国产精品专区欧美| 久久精品人人爽人人爽视色| 久久久国产欧美日韩av| 国产色婷婷99| 亚洲国产精品999| 国产精品一区二区在线观看99| 天堂中文最新版在线下载| 制服丝袜香蕉在线| 亚洲美女黄色视频免费看| 国产欧美日韩综合在线一区二区| 亚洲美女视频黄频| 日韩精品有码人妻一区| 亚洲综合色网址| 久久 成人 亚洲| 亚洲高清免费不卡视频| 99久国产av精品国产电影| 夜夜骑夜夜射夜夜干| 美女福利国产在线| 多毛熟女@视频| 婷婷色综合大香蕉| 免费高清在线观看视频在线观看| 巨乳人妻的诱惑在线观看| 少妇精品久久久久久久| 在线观看一区二区三区激情| 久久久久精品性色| 中文字幕最新亚洲高清| av女优亚洲男人天堂| 日韩制服丝袜自拍偷拍| 久久精品人人爽人人爽视色| 国产成人免费观看mmmm| 人妻人人澡人人爽人人| 丁香六月天网| 日韩电影二区| 欧美xxxx性猛交bbbb| 亚洲综合精品二区| 国产亚洲精品久久久com| 亚洲国产成人一精品久久久| 亚洲中文av在线| 国产 一区精品| 亚洲欧洲精品一区二区精品久久久 | 另类亚洲欧美激情| 日韩视频在线欧美| 高清黄色对白视频在线免费看| 一区二区三区精品91| 寂寞人妻少妇视频99o| 成年人免费黄色播放视频| 成人免费观看视频高清| 成人国产麻豆网| 日日摸夜夜添夜夜爱| 亚洲精品中文字幕在线视频| 久久鲁丝午夜福利片| av播播在线观看一区| 国产日韩一区二区三区精品不卡| 精品少妇久久久久久888优播| 国产一区亚洲一区在线观看| 亚洲精品,欧美精品| 交换朋友夫妻互换小说| 高清av免费在线| 亚洲av免费高清在线观看| 丁香六月天网| 在线看a的网站| 色吧在线观看| 午夜91福利影院| 日韩免费高清中文字幕av| 久久韩国三级中文字幕| 久久 成人 亚洲| 国产精品一区www在线观看| 国产免费又黄又爽又色| 久久免费观看电影| 两个人看的免费小视频| 香蕉精品网在线| 久久久久久久久久成人| 老司机影院毛片| 18禁国产床啪视频网站| 亚洲av在线观看美女高潮| 免费黄网站久久成人精品| 久久久久久久久久成人| 日韩不卡一区二区三区视频在线| 国产成人a∨麻豆精品| 欧美日韩亚洲高清精品| 一级毛片黄色毛片免费观看视频| 黑人高潮一二区| 人人妻人人添人人爽欧美一区卜| 在线观看国产h片| 国产麻豆69| 亚洲成国产人片在线观看| 国产免费一区二区三区四区乱码| 免费黄频网站在线观看国产| 国产成人精品久久久久久| 亚洲av.av天堂| 在线观看免费高清a一片| 性色av一级| 九九爱精品视频在线观看| 在线亚洲精品国产二区图片欧美| 中文字幕免费在线视频6| 国产淫语在线视频| 18在线观看网站| 久久影院123| 建设人人有责人人尽责人人享有的| 18禁国产床啪视频网站| 久久久久精品久久久久真实原创| 大片免费播放器 马上看| 中文字幕精品免费在线观看视频 | 成人二区视频| 国产1区2区3区精品| 黄色 视频免费看| 看免费成人av毛片| 日本vs欧美在线观看视频| 一本久久精品| 赤兔流量卡办理| 成人手机av| 人人澡人人妻人| 精品少妇内射三级| av国产久精品久网站免费入址| 天堂8中文在线网| 久久精品夜色国产| 看免费成人av毛片| 啦啦啦在线观看免费高清www| 水蜜桃什么品种好| 亚洲国产日韩一区二区| 日韩三级伦理在线观看| 热re99久久精品国产66热6| 少妇的逼水好多| 国产免费现黄频在线看| 亚洲成av片中文字幕在线观看 | 日本91视频免费播放| 满18在线观看网站| 欧美性感艳星| 国产乱人偷精品视频| 熟女电影av网| 精品第一国产精品| 欧美亚洲日本最大视频资源| 蜜桃在线观看..| 欧美 日韩 精品 国产| 国产成人aa在线观看| 一区二区日韩欧美中文字幕 | 夜夜骑夜夜射夜夜干| 男女无遮挡免费网站观看| 精品人妻偷拍中文字幕| 9色porny在线观看| 亚洲图色成人| 十八禁网站网址无遮挡| 在线天堂中文资源库| 亚洲,欧美,日韩| 少妇高潮的动态图| 亚洲欧洲日产国产| 三上悠亚av全集在线观看| 国产成人一区二区在线| 国产一区亚洲一区在线观看| 亚洲精品第二区| 99re6热这里在线精品视频| 91国产中文字幕| 精品久久国产蜜桃| 免费观看性生交大片5| 亚洲av.av天堂| 七月丁香在线播放| a级毛片黄视频| kizo精华| 亚洲一码二码三码区别大吗| 亚洲精品一二三| 伊人久久国产一区二区| 免费看av在线观看网站| 久久久久久久久久成人| 国产欧美另类精品又又久久亚洲欧美| 青春草国产在线视频| 久久鲁丝午夜福利片| 色吧在线观看| 美女视频免费永久观看网站| 欧美日韩一区二区视频在线观看视频在线| 日韩av不卡免费在线播放| 日产精品乱码卡一卡2卡三| 日韩制服骚丝袜av| 大片电影免费在线观看免费| 亚洲精品第二区| 国产欧美亚洲国产| 夫妻午夜视频| 国产精品蜜桃在线观看| 午夜福利在线观看免费完整高清在| 夫妻午夜视频| 久久人人97超碰香蕉20202| 熟女av电影| 亚洲成色77777| 天天操日日干夜夜撸| 国产熟女午夜一区二区三区| 亚洲av免费高清在线观看| 熟女av电影| 亚洲国产av影院在线观看| 免费日韩欧美在线观看| av视频免费观看在线观看| 日本午夜av视频| 精品久久久精品久久久| 亚洲av综合色区一区| 亚洲av综合色区一区| 亚洲第一av免费看| 捣出白浆h1v1| 啦啦啦在线观看免费高清www| 日韩一区二区三区影片| 久久人妻熟女aⅴ| 亚洲欧洲精品一区二区精品久久久 | 日韩一本色道免费dvd| 精品福利永久在线观看| 制服诱惑二区| 色婷婷av一区二区三区视频| 亚洲欧洲精品一区二区精品久久久 | 国产一区二区三区av在线| av在线播放精品| 国产一区有黄有色的免费视频| 亚洲精品自拍成人| 久久久久久久精品精品| 精品国产国语对白av| 韩国高清视频一区二区三区| 乱码一卡2卡4卡精品| 欧美最新免费一区二区三区| 午夜久久久在线观看| 久久久久久人人人人人| 2022亚洲国产成人精品| 春色校园在线视频观看| 看免费成人av毛片| 搡老乐熟女国产| 日韩制服丝袜自拍偷拍| 欧美xxⅹ黑人| av免费观看日本| 亚洲一级一片aⅴ在线观看| 大片电影免费在线观看免费| videossex国产| 日本欧美视频一区| 69精品国产乱码久久久| 老司机亚洲免费影院| 不卡视频在线观看欧美| 欧美国产精品一级二级三级| 日韩制服丝袜自拍偷拍| 日本与韩国留学比较| 日日爽夜夜爽网站| a级毛色黄片| 丝袜喷水一区| 久久婷婷青草| 国产成人精品无人区| 中国国产av一级| 国产亚洲av片在线观看秒播厂| 免费人妻精品一区二区三区视频| 亚洲四区av| 看免费成人av毛片| 97在线视频观看| 高清毛片免费看| 亚洲av欧美aⅴ国产| 男女国产视频网站| www日本在线高清视频| 狠狠精品人妻久久久久久综合| 国产不卡av网站在线观看| 欧美亚洲 丝袜 人妻 在线| 精品国产国语对白av| 色5月婷婷丁香| 大码成人一级视频| videos熟女内射| 亚洲精品,欧美精品| 欧美 亚洲 国产 日韩一| 亚洲欧美成人精品一区二区| 免费高清在线观看视频在线观看| 亚洲av福利一区| 卡戴珊不雅视频在线播放| 中文字幕制服av| 日韩中字成人| 亚洲精品成人av观看孕妇| 国产av码专区亚洲av| 午夜免费鲁丝| 国产精品人妻久久久久久| 黄片无遮挡物在线观看| 亚洲 欧美一区二区三区| 黄片播放在线免费| 亚洲高清免费不卡视频| 免费少妇av软件| a级毛片在线看网站| 国产成人a∨麻豆精品| 熟女av电影| 大片电影免费在线观看免费| 老女人水多毛片| 国产一区二区在线观看日韩| 激情视频va一区二区三区| 国产免费现黄频在线看| 国产精品人妻久久久影院| www.熟女人妻精品国产 | 91成人精品电影| 亚洲精品视频女| 99久久精品国产国产毛片| 国产一区二区激情短视频 | 亚洲成人一二三区av| 天天操日日干夜夜撸| 久久久久人妻精品一区果冻| 国产黄色视频一区二区在线观看| 久久精品久久久久久久性| 免费看光身美女| 日本vs欧美在线观看视频| 婷婷成人精品国产| 91aial.com中文字幕在线观看| 国产 精品1| 九九在线视频观看精品| 丝袜美足系列| 岛国毛片在线播放| 午夜免费鲁丝| 亚洲在久久综合| 一本—道久久a久久精品蜜桃钙片| 丰满少妇做爰视频| 在线观看www视频免费| 少妇的逼水好多| 成人手机av| 日韩 亚洲 欧美在线| 国产av精品麻豆| 中文字幕av电影在线播放| 国产成人午夜福利电影在线观看| 免费av中文字幕在线| 亚洲欧美色中文字幕在线| 极品人妻少妇av视频| 亚洲精品中文字幕在线视频| 国产免费福利视频在线观看| 成人影院久久| 99热6这里只有精品| 人妻 亚洲 视频| 国产精品 国内视频| 啦啦啦视频在线资源免费观看| 18+在线观看网站| 精品久久蜜臀av无| 免费少妇av软件| av免费在线看不卡| av线在线观看网站| 一区二区av电影网| 在线亚洲精品国产二区图片欧美| 看非洲黑人一级黄片| 男人爽女人下面视频在线观看| 亚洲精华国产精华液的使用体验| 2022亚洲国产成人精品| 黑人猛操日本美女一级片| 9热在线视频观看99| 国产精品.久久久| 黄片播放在线免费| 国产精品人妻久久久久久| 国产精品三级大全| 欧美精品高潮呻吟av久久| 日本vs欧美在线观看视频| 久久人人97超碰香蕉20202| 69精品国产乱码久久久| 黄网站色视频无遮挡免费观看| 亚洲中文av在线| 各种免费的搞黄视频| 99久国产av精品国产电影| 久久久a久久爽久久v久久| 亚洲av日韩在线播放| 丰满少妇做爰视频| 香蕉丝袜av| 精品少妇内射三级| 久久久国产精品麻豆| 欧美日韩视频高清一区二区三区二| 考比视频在线观看| 在线观看国产h片| 一区二区日韩欧美中文字幕 | 一边摸一边做爽爽视频免费| 国产免费现黄频在线看| 亚洲精品乱久久久久久| 午夜福利影视在线免费观看| 欧美变态另类bdsm刘玥| 伊人亚洲综合成人网| 久久热在线av| 日韩成人伦理影院| 久久久久人妻精品一区果冻| 色5月婷婷丁香| 美女xxoo啪啪120秒动态图| 午夜免费男女啪啪视频观看| 国产亚洲最大av| 男人添女人高潮全过程视频| 日韩伦理黄色片| 在现免费观看毛片| av视频免费观看在线观看| 波多野结衣一区麻豆| 亚洲国产最新在线播放| 夜夜爽夜夜爽视频| 欧美精品一区二区免费开放| 国产无遮挡羞羞视频在线观看| 在线观看免费日韩欧美大片| 超碰97精品在线观看| 亚洲精品,欧美精品| 丰满饥渴人妻一区二区三| 这个男人来自地球电影免费观看 | 免费观看性生交大片5| 又黄又爽又刺激的免费视频.| 久久免费观看电影| 人人妻人人添人人爽欧美一区卜| 极品少妇高潮喷水抽搐| 欧美bdsm另类| 久久久久久人妻| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 十分钟在线观看高清视频www| 亚洲少妇的诱惑av| 亚洲国产欧美在线一区| 一级爰片在线观看| 久久久国产一区二区| 成人综合一区亚洲| 亚洲精品乱久久久久久| 飞空精品影院首页| 香蕉丝袜av| 成人18禁高潮啪啪吃奶动态图| 亚洲少妇的诱惑av| 亚洲一区二区三区欧美精品| 美国免费a级毛片| 久久精品熟女亚洲av麻豆精品| 亚洲av成人精品一二三区| 成人黄色视频免费在线看| 黑人欧美特级aaaaaa片| 久久久久久久国产电影| 国产精品国产三级专区第一集| 亚洲精品av麻豆狂野| 久热这里只有精品99| 日本av免费视频播放| 亚洲精品国产av成人精品| 成年人免费黄色播放视频| 狂野欧美激情性xxxx在线观看| av在线老鸭窝| 色婷婷av一区二区三区视频| 亚洲成人一二三区av| 又粗又硬又长又爽又黄的视频| 狂野欧美激情性bbbbbb| 成年美女黄网站色视频大全免费| 视频中文字幕在线观看| 99久久人妻综合| 亚洲激情五月婷婷啪啪| 国产成人a∨麻豆精品| 深夜精品福利| 一级a做视频免费观看| 亚洲美女搞黄在线观看| 亚洲精品自拍成人| 精品少妇久久久久久888优播| 欧美 日韩 精品 国产| 男的添女的下面高潮视频| 一本久久精品| 九草在线视频观看| 男女高潮啪啪啪动态图| 狠狠精品人妻久久久久久综合| 69精品国产乱码久久久| videos熟女内射| 亚洲综合色网址| 考比视频在线观看| 久久精品久久精品一区二区三区| 人体艺术视频欧美日本| 咕卡用的链子| 久久久久久久久久久免费av| 制服人妻中文乱码| 久久精品国产亚洲av天美| 视频在线观看一区二区三区| 美女福利国产在线| 国产毛片在线视频| 一级片'在线观看视频| 国产一区有黄有色的免费视频| 亚洲内射少妇av| 亚洲人成网站在线观看播放| 亚洲欧美日韩卡通动漫| 精品国产露脸久久av麻豆| 国产精品久久久久久精品电影小说| 汤姆久久久久久久影院中文字幕| 久久久久久久国产电影| av又黄又爽大尺度在线免费看| 秋霞伦理黄片| 91在线精品国自产拍蜜月| 久久亚洲国产成人精品v| 久久久国产一区二区| 国产精品一二三区在线看| 欧美日韩综合久久久久久| 亚洲色图综合在线观看| 亚洲久久久国产精品| 最近中文字幕高清免费大全6| videosex国产| 曰老女人黄片| 香蕉国产在线看| 久久精品国产a三级三级三级| 777米奇影视久久| 亚洲精品一二三| 国产成人免费无遮挡视频| 肉色欧美久久久久久久蜜桃| 少妇熟女欧美另类| 激情视频va一区二区三区| 视频在线观看一区二区三区| 亚洲精品国产av蜜桃| 亚洲精品乱久久久久久| a级片在线免费高清观看视频| 久久午夜综合久久蜜桃| 久久久欧美国产精品| 中文精品一卡2卡3卡4更新| 乱人伦中国视频| 亚洲精品,欧美精品| 97精品久久久久久久久久精品| 亚洲av日韩在线播放| 中文字幕精品免费在线观看视频 | 少妇猛男粗大的猛烈进出视频| 成年美女黄网站色视频大全免费| 黄网站色视频无遮挡免费观看| 搡女人真爽免费视频火全软件| 成人综合一区亚洲| 天美传媒精品一区二区| 99国产综合亚洲精品| 最近的中文字幕免费完整| 欧美亚洲 丝袜 人妻 在线| 国产精品国产三级国产专区5o| 伊人亚洲综合成人网| 亚洲av免费高清在线观看| 国产男女超爽视频在线观看| 看免费av毛片| 亚洲精品一二三| 国产高清国产精品国产三级| 久久精品国产亚洲av涩爱| 人妻人人澡人人爽人人| 免费高清在线观看视频在线观看| 亚洲精品一区蜜桃| 伦理电影免费视频| 波野结衣二区三区在线| 少妇人妻精品综合一区二区| 国产精品秋霞免费鲁丝片| 亚洲精品久久成人aⅴ小说| 视频中文字幕在线观看| 肉色欧美久久久久久久蜜桃| 国内精品宾馆在线| 在线看a的网站| 色婷婷久久久亚洲欧美| 全区人妻精品视频| 99久久综合免费| 国产白丝娇喘喷水9色精品| 最近中文字幕高清免费大全6| 亚洲,欧美,日韩| 免费观看无遮挡的男女| 午夜福利在线观看免费完整高清在| 性色av一级| 一本—道久久a久久精品蜜桃钙片| 少妇的逼水好多| 精品少妇内射三级| 日韩av在线免费看完整版不卡| 自线自在国产av| 欧美国产精品一级二级三级| 90打野战视频偷拍视频| 亚洲五月色婷婷综合| 只有这里有精品99| 蜜桃国产av成人99| 超色免费av| 亚洲精品第二区| 亚洲精品日韩在线中文字幕| 欧美精品亚洲一区二区| 久久免费观看电影| 熟女人妻精品中文字幕| 伊人亚洲综合成人网| 久久精品久久久久久久性| 亚洲一级一片aⅴ在线观看| 国产成人一区二区在线| 99久国产av精品国产电影| 一区在线观看完整版| 男人添女人高潮全过程视频| 国产1区2区3区精品| a级毛片在线看网站| av线在线观看网站| 在线观看国产h片| 你懂的网址亚洲精品在线观看| 啦啦啦啦在线视频资源| 99热国产这里只有精品6| 亚洲国产精品国产精品| 国产麻豆69| 亚洲精品一区蜜桃| 精品一区二区三区四区五区乱码 | 亚洲少妇的诱惑av| 亚洲av电影在线观看一区二区三区| 国产av码专区亚洲av| 啦啦啦啦在线视频资源| 欧美精品国产亚洲| 国产一区二区三区综合在线观看 | 国产一区二区三区av在线| 天天影视国产精品| 免费观看无遮挡的男女| 香蕉丝袜av| 黄色视频在线播放观看不卡| 亚洲精品日韩在线中文字幕| 国产成人aa在线观看| 视频在线观看一区二区三区| 久久人人爽人人片av| 成人黄色视频免费在线看| 欧美 日韩 精品 国产| 精品一区二区免费观看| 成人午夜精彩视频在线观看| 最后的刺客免费高清国语| 香蕉国产在线看| 免费人成在线观看视频色| 精品少妇内射三级| 在线观看免费日韩欧美大片| 夜夜爽夜夜爽视频| 毛片一级片免费看久久久久| av在线观看视频网站免费| 日本欧美视频一区| 91精品伊人久久大香线蕉| 蜜臀久久99精品久久宅男| 精品人妻一区二区三区麻豆| 亚洲内射少妇av| 18禁国产床啪视频网站| 搡老乐熟女国产| 亚洲精品美女久久av网站| 日本av手机在线免费观看| 国产高清不卡午夜福利| 99久国产av精品国产电影| 母亲3免费完整高清在线观看 | 久久久久久久亚洲中文字幕| 黄色毛片三级朝国网站| 最新的欧美精品一区二区| 亚洲第一av免费看| 毛片一级片免费看久久久久| 国产精品偷伦视频观看了| 日韩电影二区| 成人无遮挡网站| 国产熟女午夜一区二区三区| 赤兔流量卡办理| 最后的刺客免费高清国语| 一边摸一边做爽爽视频免费| 18在线观看网站|