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

    基于α-β-γ濾波的迫擊炮定位

    2015-02-17 03:54:06張強(qiáng)
    電訊技術(shù) 2015年6期
    關(guān)鍵詞:龍格濾波

    ?

    基于α-β-γ濾波的迫擊炮定位

    張強(qiáng)**

    (中國西南電子技術(shù)研究所,成都 610036)

    摘要:針對(duì)迫擊炮的偵察定位問題,介紹了一種基于α-β-γ濾波的定位方法。該方法首先根據(jù)外彈道學(xué)確定彈丸飛行的質(zhì)心運(yùn)動(dòng)方程,然后利用常增益的α-β-γ濾波跟蹤求解彈丸的運(yùn)動(dòng)參數(shù),最后通過龍格-庫塔(Runge-Kutta)方法求解質(zhì)心運(yùn)動(dòng)方程,并外推求得炮位位置。仿真結(jié)果表明,該方法在保證較高迫擊炮定位精度的同時(shí)實(shí)時(shí)性能較傳統(tǒng)方法大幅提升,有顯著的工程意義。

    關(guān)鍵詞:迫擊炮定位;質(zhì)心運(yùn)動(dòng)方程;α-β-γ濾波;龍格-庫塔方法

    1引言

    迫擊炮具有構(gòu)造簡(jiǎn)單、操作靈活、造價(jià)低廉、彈道彎曲、最小射程近、射速快、可實(shí)時(shí)伴隨步兵作戰(zhàn)等優(yōu)點(diǎn),對(duì)開闊地及掩體內(nèi)目標(biāo)、各種野戰(zhàn)工事甚至高大障礙物(如山坡)背后目標(biāo)都有著良好的毀傷破壞作用,作為步兵近距離火力支援的有效武器,現(xiàn)今仍被各國軍隊(duì)大量裝備[1-3]。為了提高我方炮兵的作戰(zhàn)反應(yīng)速度和火炮命中概率,各國都在不斷加強(qiáng)火炮偵察定位雷達(dá)的研究,以實(shí)現(xiàn)對(duì)敵方炮火的有效壓制,典型的炮位偵察雷達(dá)有美國AN/TPQ-36 、AN/TPQ-37“火力發(fā)現(xiàn)者”[1,4]和輕型反迫擊炮雷達(dá)(LCMR)。火炮定位算法作為炮位偵察定位雷達(dá)的關(guān)鍵技術(shù)也日益得到重視。文獻(xiàn)[5-6]所提出的外彈道仿真方法未通過濾波跟蹤算法對(duì)彈丸的運(yùn)動(dòng)參數(shù)進(jìn)行解算,因此反推炮位的起點(diǎn)易受隨機(jī)噪聲影響,反推炮位的誤差較大。文獻(xiàn)[7]引入粒子濾波的方法來辨識(shí)彈道參數(shù)和彈丸的運(yùn)動(dòng)參數(shù),然后再外推敵方炮位,該方法定位精度高,但是算法復(fù)雜、實(shí)時(shí)性差,不能滿足工程應(yīng)用的需要。實(shí)踐表明,卡爾曼濾波[8]雖然復(fù)雜程度不及文獻(xiàn)[7]中的粒子濾波方法,但是它的運(yùn)算時(shí)間也比較長,實(shí)時(shí)性仍然不夠好。因此,有必要尋求一種在保證一定濾波精度條件下(即工程上允許的近似)運(yùn)算量小、速度快的濾波器。理論研究表明,在一定的限定條件下,α-β-γ濾波器與卡爾曼濾波器具有相一致的最佳增益[9],但是運(yùn)算量卻大大減小。本文研究α-β-γ 濾波在迫擊炮定位中的應(yīng)用。

    2迫擊炮彈丸運(yùn)動(dòng)數(shù)學(xué)模型

    迫擊炮彈丸的運(yùn)動(dòng)彈道是一種類拋物線,但是由于受到空氣阻力等因素影響,實(shí)際彈道與拋物線又有較大的差異,特別是對(duì)低射角的飛行彈道而言,若用拋物線方程來逼近,產(chǎn)生的定位誤差將不可接受,甚至可能是發(fā)散的,無法外推,因此,本文采用彈丸的質(zhì)心運(yùn)動(dòng)方程來描述其飛行軌跡[6]。

    迫擊炮彈飛行軌跡示意圖如圖1所示,彈丸的飛行軌跡位于XYZ三維笛卡爾坐標(biāo)系內(nèi),炮位位于H點(diǎn),炮位偵察雷達(dá)位于X軸方向。

    圖1 迫擊炮彈飛行軌跡示意圖

    那么以時(shí)間t為自變量,在三維直角坐標(biāo)系中彈丸質(zhì)心運(yùn)動(dòng)方程可描述為

    (1)

    彈丸的質(zhì)心運(yùn)動(dòng)方程是一種理想情況,它把彈丸看作質(zhì)點(diǎn),在空中飛行時(shí)只受到重力和空氣阻力的影響,且不考慮馬格洛斯力、升力和科氏力等,它是一種氣動(dòng)參數(shù)少但能較準(zhǔn)確反映飛行彈丸運(yùn)動(dòng)規(guī)律的彈道方程。彈丸質(zhì)心運(yùn)動(dòng)方程雖然在全彈道上與真實(shí)彈道有一定差距,但是完全能夠滿足炮位偵察的要求。

    3α-β-γ濾波解算目標(biāo)運(yùn)動(dòng)參數(shù)

    為了獲得準(zhǔn)確的彈丸運(yùn)動(dòng)參數(shù),有必要通過濾波跟蹤對(duì)雷達(dá)量測(cè)數(shù)據(jù)進(jìn)行平滑,并獲得后續(xù)炮位反推的起點(diǎn)。典型的濾波跟蹤方法就是卡爾曼濾波,但是卡爾曼濾波的增益計(jì)算過程較為復(fù)雜,需要耗費(fèi)大量時(shí)間,影響算法的實(shí)時(shí)性。α-β-γ濾波作為卡爾曼濾波的特殊形式,由于它是一種常增益濾波器,計(jì)算復(fù)雜度極大地下降,但是針對(duì)勻加速運(yùn)動(dòng)模型的目標(biāo)卻和傳統(tǒng)卡爾曼濾波有相一致的最佳增益[9]。而迫擊炮彈丸由于僅受到萬有引力和空氣阻力作用可以被認(rèn)為處于勻加速運(yùn)動(dòng)狀態(tài),所以引入α-β-γ濾波器來對(duì)雷達(dá)量測(cè)數(shù)據(jù)進(jìn)行平滑并解算目標(biāo)的運(yùn)動(dòng)參數(shù)。對(duì)于笛卡爾坐標(biāo)系下九維向量的α-β-γ平滑濾波可描述如下。

    彈丸運(yùn)動(dòng)的狀態(tài)方程可描述為

    X(k+1)=F(k)X(k)+Γ(k)v(k),k=0,1,2,…,99。

    (2)

    式中,狀態(tài)向量為X(k)=[x vxaxy vyayz vzaz]′;狀態(tài)轉(zhuǎn)移矩陣

    過程噪聲分布矩陣

    Γ(k)=[T2/2 T 1 T2/2 T 1 T2/2 T 1]′;v(k)為過程噪聲,是零均值高斯白噪聲。

    彈丸的量測(cè)方程為

    Z(k)=H(k)X(k)+W(k)。

    (3)

    式中,量測(cè)矩陣

    W(k)為量測(cè)噪聲,是零均值的白噪聲,與過程噪聲序列相互獨(dú)立。

    狀態(tài)一步預(yù)測(cè)為

    (4)

    常增益的表達(dá)式為

    K(k+1)=

    [αβ/Tγ/T2αβ/Tγ/T2αβ/Tγ/T2]。

    (5)

    公式(5)中的常增益是機(jī)動(dòng)指標(biāo)的函數(shù),工程上常用與采樣時(shí)刻k有關(guān)的方法來確定α、β和γ的值,如公式(6)~(8)所示:

    (6)

    (7)

    (8)

    對(duì)α來說,從k=1開始取值;對(duì)β來說,從k=2開始取值;對(duì)γ來說,從k=3開始取值。α、β、γ值與k的關(guān)系如表1所示。

    表1 α、β、γ值與k的關(guān)系

    最后可得α-β-γ濾波的狀態(tài)更新方程為

    (9)

    由于基于α-β-γ濾波的目標(biāo)解算方法在計(jì)算增益時(shí)不需要計(jì)算新息協(xié)方差和更新協(xié)方差向量等過程,計(jì)算復(fù)雜度大大降低。

    4龍格-庫塔法炮位反演

    彈丸的質(zhì)心運(yùn)動(dòng)方程(1)是一個(gè)微分方程,可用龍格-庫塔(Runge-Kutta)方法來實(shí)現(xiàn)微分方程的解[6-7]。

    對(duì)于微分方程如公式(10),那么用四階龍格-庫塔方法的求解如公式(11)所示:

    x′=f(t,x),x(t0)=x0;

    (10)

    xi+1=xi+hφ(ti,xi,h),i=0,1,…。

    (11)

    對(duì)于用龍格-庫塔方法進(jìn)行微分方程數(shù)值求解而言,若x(t+h)-[x(t)+hφ(t,x(t),h)]=O(hp+1),則稱以上迭代公式為p階公式,p的大小反映了截?cái)嗾`差的高低,越高階精度越高。本文所用的四階龍格-庫塔方法其p=4,精度較高。

    龍格-庫塔炮位反推的起點(diǎn)為α-β-γ濾波獲得的精度較高的處于下降段的彈丸軌跡點(diǎn),選取恰當(dāng)?shù)臅r(shí)間間隔h即可獲得較高的炮位反演精度,當(dāng)炮位迭代反推至z=z0時(shí)的彈丸坐標(biāo)(x0,y0)即為迫擊炮位置所在,其中z0為彈丸飛行曲線與地形坐標(biāo)交點(diǎn)的高度值。

    5算法流程

    基于α-β-γ濾波的迫擊炮定位方法的詳細(xì)算法流程圖如圖2所示。

    圖2 定位算法流程圖

    5.1彈道參數(shù)的初始化

    質(zhì)心運(yùn)動(dòng)方程主要彈道參數(shù)根據(jù)《彈箭外彈道學(xué)》[10]確定:

    (1)彈道系數(shù)c=0.8613;

    (2)地面標(biāo)準(zhǔn)虛溫值τon=15℃,z高度處的虛溫值τ=τon-0.006328z;

    (3)空氣密度函數(shù)

    H(z)=(1-2.1904×10-5z)4.399;

    (4)空氣阻力函數(shù)

    式中,v的單位為m/s;

    (5)重力加速度g=9.8 m/s2。

    5.2量測(cè)值的初始化

    對(duì)于偵察定位雷達(dá)每次給出的目標(biāo)測(cè)量值為斜距、方位角和俯仰角,即(R,α,β)。為了方便算法實(shí)現(xiàn),將其轉(zhuǎn)換到雷達(dá)笛卡爾坐標(biāo)系下為Z=(x,y,z),那么系統(tǒng)的初始狀態(tài)向量可利用前三個(gè)時(shí)刻的雷達(dá)測(cè)量值Z(0)、Z(1)和Z(2)確定,即

    5.3濾波平滑

    偵察雷達(dá)需要給出13個(gè)測(cè)量值,其中時(shí)間間隔為T=0.1s,系統(tǒng)根據(jù)第3節(jié)介紹的常增益α-β-γ濾波從k=3時(shí)開始濾波平滑,最終獲得k=12的彈丸運(yùn)動(dòng)參數(shù)(x1vx1ax1y1vy1ay1z1vz1az1),即炮位反推的起點(diǎn)。

    5.4炮位反推

    選取時(shí)間間隔h=0.1s,根據(jù)第4節(jié)介紹的龍格-庫塔方法進(jìn)行炮位反推,當(dāng)炮位迭代反推至z=z0時(shí)的彈丸坐標(biāo)(x0,y0)即為迫擊炮位置所在,其中z0為彈丸飛行曲線與地形坐標(biāo)交點(diǎn)的高度值。

    6仿真與分析

    以美國某型81mm口徑的迫擊炮為例,用Matlab進(jìn)行仿真實(shí)驗(yàn)。設(shè)定彈丸在YOZ平面內(nèi)飛行,地形與水平面相同,那么可從網(wǎng)絡(luò)查閱獲得該型迫擊炮在發(fā)射角為45°和60°時(shí)的相關(guān)運(yùn)動(dòng)參數(shù),如表2所示。

    表2 迫擊炮相關(guān)參數(shù)

    仿真實(shí)驗(yàn)中首先根據(jù)表2中迫擊炮相關(guān)參數(shù)得到彈丸飛行路徑圖,然后取飛行路徑下降段的連續(xù)13點(diǎn)采樣值作為雷達(dá)的量測(cè)值,并據(jù)此進(jìn)行α-β-γ濾波和炮位反推,同時(shí)將本文方法與基于卡爾曼濾波[8]和粒子濾波[7]的方法進(jìn)行比較,反推的炮位坐標(biāo)與真實(shí)坐標(biāo)如表3所示, 可見雖然本文炮位反演方法的絕對(duì)誤差分別為36 m和29 m,高于其他兩種方法,但是仍具有很高的炮位反推精度,對(duì)實(shí)際工程應(yīng)用而言迫擊炮定位精度小于50 m就是非常有意義的。

    表3 定位精度比較

    本文算法與基于卡爾曼濾波和粒子濾波的炮位反演算法計(jì)算時(shí)間比較如表4所示,表中以本文方法在45°發(fā)射角的反推時(shí)間為單位1,其余時(shí)間為單位1的倍數(shù),可見本文方法在計(jì)算時(shí)間上優(yōu)勢(shì)明顯。

    表4 計(jì)算時(shí)間比較

    仿真分析中的彈道參數(shù)設(shè)置是具有通用性的,只是未充分考慮彈丸目標(biāo)的起伏特性和雷達(dá)測(cè)量噪聲對(duì)定位結(jié)果的影響,在實(shí)際工程中定位精度還有降低的可能,但該方法亦不失為一種性能較好的迫擊炮定位方法。

    7結(jié)束語

    本文介紹了一種基于α-β-γ濾波的迫擊炮偵察定位方法。該方法利用α-β-γ濾波跟蹤獲得迫擊炮彈丸飛行的相關(guān)參數(shù),然后用四階龍格-庫塔方法來解算彈丸的軌跡運(yùn)動(dòng)方程,實(shí)現(xiàn)對(duì)炮位的準(zhǔn)確外推。仿真結(jié)果表明,該方法雖然偵察定位精度略微有所下降,但是實(shí)時(shí)性能相對(duì)文獻(xiàn)[7-8]的方法有極大提升,有顯著的工程應(yīng)用價(jià)值。在后續(xù)的工作中有必要進(jìn)一步研究該方法推廣到其他艦炮、坦克炮等火炮的定位效果。

    參考文獻(xiàn):

    [1]朱元亮.國柱陣炮位偵察雷達(dá)關(guān)鍵技術(shù)的研究[D]. 西安:西安電子科技大學(xué),2010:2-15.

    ZHU Yuangliang. Study on the key techniques of artillery location radar with the cylindrical array antenna[D]. Xi′an:Xidian University,2010:2-15.(in Chinese)

    [2]郭蔚.炮位偵校雷達(dá)信號(hào)處理關(guān)鍵技術(shù)研究[D]. 南京:南京理工大學(xué),2013.

    GUO Wei. Key technology of Artillery locating radar signal processing[D]. Nanjing :Nanjing University of Science and Technology,2013.(in Chinese)

    [3]楊翔,王雨時(shí),聞泉.迫擊炮彈空氣動(dòng)力特性攻角系數(shù)數(shù)值研究[J].彈箭與制導(dǎo)學(xué)報(bào),2014,34(2):139-142.

    YANG Xiang,WANG Yushi,WEN Quan. Numerical simulation on aerodynamic characteristics about angle of attack coefficient of mortar projectiles[J]. Journal of Projectiles,Rockets,Missiles and Cuidance,2014,34(2):139-142.(in Chinese)

    [4]王東,張鵬,梁開莉.武器定位雷達(dá)偵察模式下多目標(biāo)能力研究[J].現(xiàn)代雷達(dá),2011,33(11): 6-10.

    WANG Dong,ZHANG Peng,LIANG Kaili. A study on multiple-target locating capability of weapon locating radar in reconnaissance mode[J]. Modern Radar,2011,33(11): 6-10.(in Chinese)

    [5]楊獻(xiàn)輝,趙俊利.基于MATLAB的迫擊炮外彈道模型仿真研究[J].裝備制造技術(shù),2011(2):24-26.

    YANG Xianhui,ZHAO Junli. A simulation research on the external trajectory of mortars based on MATLAB[J]. Equipment Manufacturing Technology,2011(2):24-26.(in Chinese)

    [6]伍建輝,董亮.基于Matlab迫擊炮外彈道仿真[J].火控雷達(dá)技術(shù),2014,43(2):39-42.

    WU Jianhui,DONG Liang. Simulation on external trajectory of mortar based on matlab[J]. Fire Control Radar Technology,2014,43(2):39-42.(in Chinese)

    [7]胡獻(xiàn)君,王航宇,周德超.粒子濾波跟蹤彈道下降段在炮位偵察中的應(yīng)用[J].現(xiàn)代防御技術(shù),2013,41(1):104-19.

    HU Xianjun,WANG Hangyu,ZHOU Dechao. Application of tracking descend trajectory with particle filter on detecting cannon location[J]. Modern Defence Technology,2013,41(1):104-19.(in Chinese)

    [8]唐吉祥,沈云山. 卡爾曼濾波在炮位偵察雷達(dá)中的應(yīng)用[J].電光與控制,2004,11(1):48-51.

    TANG Jixiang,SHEN Yunshan.Application of Kalman filter in cannon locating radar[J].Electronics optics & Control,2004,11(1):48-51.(in Chinese)

    [9]何友,修建娟,張晶煒,等.雷達(dá)數(shù)據(jù)處理及應(yīng)用[M].北京:電子工業(yè)出版社,2012:15-50.

    HE You,XIU Jianjuan,ZHANG Jingwei,et al. Radar Data Processing with Applications[M]. Beijing: Publishing House of Electronics Industry,2012:15-50.(in Chinese)

    [10]韓子鵬.彈箭外彈道學(xué)[M].北京:北京理工大學(xué)出版社,2008.

    HAN Zipeng. Exterior Ballistics of Projectiles and Rockets[M]. Beijing:Beijing Institute of Technology Press,2008.(in Chinese)

    張強(qiáng)(1986—),男,四川鄰水人,2012年于中國科學(xué)院獲博士學(xué)位,現(xiàn)為工程師,主要研究方向?yàn)閳D像處理、目標(biāo)識(shí)別和雷達(dá)系統(tǒng)設(shè)計(jì)。

    ZHANG Qiang was born in Linshui,Sichuan Province,in 1986. He received the Ph.D. degree from Chinese Academy of Sciences in 2012. He is now an engineer. His research interests include image processing,target recognition,radar system design.

    Email: zhangqiang_ioe@126.com

    引用格式:張強(qiáng).基于α-β-γ濾波的迫擊炮定位[J].電訊技術(shù),2015,55(6):634-638.[ZHANG Qiang.Trench Mortar Locating Based on α-β-γ Filter[J].Telecommunication Engineering,2015,55(6):634-638.]

    Trench Mortar Locating Based on α-β-γ Filter

    ZHANG Qiang

    (Southwest China Institute of Electronic Technology,Chengdu 610036,China)

    Abstract:For the locating problem of trench mortar,a locating method based on α-β-γ filter is presented. First,centroid moving equations are determined according to the exterior ballistics. Then the motion parameters are solved by the α-β-γ tracking filter with constant gain. Finally, the trench mortar position is deduced by solving the centroid moving equations with the Runge-Kutta method. The simulation results show that the proposed method has good locating precision with much shorter time compared with conventional method,and it is valuable for engineering.

    Key words:trench mortar locating;centroid moving equation;α-β-γ filter;Runge-Kutta method

    作者簡(jiǎn)介:

    中圖分類號(hào):TN97

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1001-893X(2015)06-0634-05

    通訊作者:*zhangqiangioe@126.comCorresponding author:zhangqiangioe@126.com

    收稿日期:*2014-09-28;修回日期:2015-04-17Received date:2014-09-28;Revised date:2015-04-17

    doi:10.3969/j.issn.1001-893x.2015.06.009

    猜你喜歡
    龍格濾波
    懸索橋吊索尾流致振非定常理論分析
    四階龍格-庫塔方法的程序設(shè)計(jì)與應(yīng)用
    基于Labview的龍格—庫塔法估算Pu內(nèi)照射攝入量
    非結(jié)構(gòu)網(wǎng)格的有限體積法研究
    一種新的InSAR干涉相位濾波方法
    一種GLONASS衛(wèi)星軌道快速計(jì)算方法
    錨段關(guān)節(jié)式電分相過電壓的龍格-庫塔解法及抑制
    基于自適應(yīng)Kalman濾波的改進(jìn)PSO算法
    RTS平滑濾波在事后姿態(tài)確定中的應(yīng)用
    基于線性正則變換的 LMS 自適應(yīng)濾波
    少妇熟女aⅴ在线视频| av福利片在线观看| 欧美bdsm另类| 99九九线精品视频在线观看视频| 人妻少妇偷人精品九色| 免费人成在线观看视频色| 亚洲精品一卡2卡三卡4卡5卡| 人人妻人人澡欧美一区二区| 国内毛片毛片毛片毛片毛片| 免费在线观看成人毛片| 精品日产1卡2卡| 少妇裸体淫交视频免费看高清| 午夜精品在线福利| 国内精品久久久久精免费| 亚洲色图av天堂| 亚洲精品一区av在线观看| 亚洲国产欧洲综合997久久,| 亚洲国产色片| 国产探花在线观看一区二区| 91久久精品电影网| 色综合站精品国产| 在线免费观看不下载黄p国产 | 日日摸夜夜添夜夜添av毛片 | 国产麻豆成人av免费视频| av在线老鸭窝| 国产视频一区二区在线看| 久久亚洲精品不卡| 中国美女看黄片| 亚洲男人的天堂狠狠| 国产伦一二天堂av在线观看| 午夜福利在线在线| 欧美日韩乱码在线| 久久久国产成人精品二区| 午夜日韩欧美国产| 国产精品福利在线免费观看| av在线天堂中文字幕| 狂野欧美白嫩少妇大欣赏| 久久久精品大字幕| 国产爱豆传媒在线观看| 淫秽高清视频在线观看| 国产精品综合久久久久久久免费| 国产男人的电影天堂91| 亚洲av免费高清在线观看| 日韩欧美国产在线观看| 国产精品亚洲一级av第二区| 一进一出好大好爽视频| 99久久精品国产国产毛片| 亚洲第一电影网av| 动漫黄色视频在线观看| 变态另类丝袜制服| 国内精品久久久久久久电影| 免费看光身美女| 亚洲自拍偷在线| 嫩草影院精品99| 国产精品美女特级片免费视频播放器| 看片在线看免费视频| 99久久久亚洲精品蜜臀av| 五月玫瑰六月丁香| 亚洲人与动物交配视频| 国产伦在线观看视频一区| av在线观看视频网站免费| 国产一区二区三区在线臀色熟女| 成年人黄色毛片网站| 九色国产91popny在线| 99视频精品全部免费 在线| 久久久成人免费电影| 国产精品久久视频播放| 国产成人a区在线观看| 国产熟女欧美一区二区| 精品乱码久久久久久99久播| 亚洲欧美日韩高清专用| 人人妻人人看人人澡| 国产高潮美女av| 哪里可以看免费的av片| 久久久久性生活片| 久久久久国产精品人妻aⅴ院| 亚洲成人中文字幕在线播放| 中国美女看黄片| 中文在线观看免费www的网站| 日本成人三级电影网站| 变态另类丝袜制服| 亚洲精品456在线播放app | 免费观看的影片在线观看| 老女人水多毛片| 成人特级黄色片久久久久久久| 国产精品国产三级国产av玫瑰| 久久九九热精品免费| 久久久精品大字幕| 亚洲乱码一区二区免费版| 超碰av人人做人人爽久久| 最新在线观看一区二区三区| 内射极品少妇av片p| 精华霜和精华液先用哪个| 岛国在线免费视频观看| 久久精品国产鲁丝片午夜精品 | 可以在线观看的亚洲视频| 日日摸夜夜添夜夜添小说| 小说图片视频综合网站| 亚洲欧美日韩东京热| 国产一区二区在线观看日韩| 色精品久久人妻99蜜桃| 99热精品在线国产| 人人妻人人澡欧美一区二区| 在线观看午夜福利视频| 日韩大尺度精品在线看网址| 午夜福利18| 99国产精品一区二区蜜桃av| 亚洲人成网站在线播| 在线观看免费视频日本深夜| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕精品亚洲无线码一区| 在线播放无遮挡| 日本色播在线视频| 国产欧美日韩一区二区精品| 韩国av一区二区三区四区| 精品久久久久久,| 亚洲七黄色美女视频| 精品99又大又爽又粗少妇毛片 | 亚洲国产色片| 日韩国内少妇激情av| 色噜噜av男人的天堂激情| 亚洲精品影视一区二区三区av| 精品久久久久久久久久久久久| 成人特级黄色片久久久久久久| 99在线人妻在线中文字幕| 简卡轻食公司| 看免费成人av毛片| 亚洲精品色激情综合| 国产成人av教育| 两性午夜刺激爽爽歪歪视频在线观看| 我的女老师完整版在线观看| 一级黄片播放器| 久久久久久久久久成人| 精品久久久久久久久亚洲 | 欧美一区二区精品小视频在线| 琪琪午夜伦伦电影理论片6080| 99热这里只有是精品在线观看| 一个人看视频在线观看www免费| 亚洲精品成人久久久久久| 国产一区二区亚洲精品在线观看| 少妇裸体淫交视频免费看高清| 国产成人av教育| 99热这里只有是精品50| 97碰自拍视频| 最近中文字幕高清免费大全6 | 无遮挡黄片免费观看| 夜夜爽天天搞| 日日夜夜操网爽| 久久精品人妻少妇| 嫩草影院入口| 尤物成人国产欧美一区二区三区| 国产乱人视频| 少妇的逼好多水| 亚洲精品国产成人久久av| 亚洲真实伦在线观看| 少妇人妻一区二区三区视频| 久久人人爽人人爽人人片va| 欧美黑人巨大hd| 最近在线观看免费完整版| 国产精品久久久久久久久免| 久久人妻av系列| 色5月婷婷丁香| 男女啪啪激烈高潮av片| 男女做爰动态图高潮gif福利片| 国产av不卡久久| 97热精品久久久久久| 国产69精品久久久久777片| 1024手机看黄色片| 亚洲男人的天堂狠狠| 国产黄a三级三级三级人| 国内少妇人妻偷人精品xxx网站| 国产爱豆传媒在线观看| 国产黄a三级三级三级人| 亚洲欧美日韩高清专用| 亚洲人与动物交配视频| 欧美激情国产日韩精品一区| 国产三级中文精品| 男女边吃奶边做爰视频| 别揉我奶头~嗯~啊~动态视频| 久久久久久久久久久丰满 | .国产精品久久| 国产私拍福利视频在线观看| 国产亚洲av嫩草精品影院| 又爽又黄a免费视频| 黄色日韩在线| 他把我摸到了高潮在线观看| 一边摸一边抽搐一进一小说| 老师上课跳d突然被开到最大视频| 欧美+日韩+精品| 婷婷色综合大香蕉| 亚洲性夜色夜夜综合| 久久精品国产99精品国产亚洲性色| 中文字幕高清在线视频| 悠悠久久av| 成人特级av手机在线观看| 日本黄色片子视频| 伦精品一区二区三区| av国产免费在线观看| 制服丝袜大香蕉在线| 欧美一区二区亚洲| 欧美色视频一区免费| 欧美丝袜亚洲另类 | 免费不卡的大黄色大毛片视频在线观看 | 天天躁日日操中文字幕| 女的被弄到高潮叫床怎么办 | 免费电影在线观看免费观看| 色哟哟·www| 嫁个100分男人电影在线观看| 久久精品国产亚洲av天美| 国内精品一区二区在线观看| 五月伊人婷婷丁香| 国产精品精品国产色婷婷| 午夜激情欧美在线| 国产欧美日韩精品一区二区| 悠悠久久av| 女生性感内裤真人,穿戴方法视频| 午夜精品久久久久久毛片777| 窝窝影院91人妻| 国内少妇人妻偷人精品xxx网站| 国产黄a三级三级三级人| 麻豆成人午夜福利视频| 日日摸夜夜添夜夜添av毛片 | 日韩欧美在线乱码| 久久99热这里只有精品18| 熟女人妻精品中文字幕| 国产视频一区二区在线看| 免费大片18禁| 日日摸夜夜添夜夜添小说| 国产主播在线观看一区二区| 午夜爱爱视频在线播放| 日本-黄色视频高清免费观看| 欧美一区二区亚洲| av专区在线播放| 一a级毛片在线观看| 亚洲四区av| 免费看光身美女| 九色成人免费人妻av| 欧美极品一区二区三区四区| 欧美日本亚洲视频在线播放| 精品久久久久久久末码| 久久精品国产亚洲av涩爱 | 色综合亚洲欧美另类图片| xxxwww97欧美| 日本免费一区二区三区高清不卡| 给我免费播放毛片高清在线观看| 成年版毛片免费区| 偷拍熟女少妇极品色| 国产色爽女视频免费观看| 亚洲自拍偷在线| 国产男靠女视频免费网站| 日本-黄色视频高清免费观看| 俄罗斯特黄特色一大片| 小蜜桃在线观看免费完整版高清| 国产亚洲91精品色在线| 久久精品国产亚洲av天美| 午夜免费成人在线视频| 麻豆国产av国片精品| 少妇高潮的动态图| 欧美一区二区亚洲| 亚洲成av人片在线播放无| 99久久中文字幕三级久久日本| 成人二区视频| 国产男靠女视频免费网站| 真人一进一出gif抽搐免费| 国产精品99久久久久久久久| 亚洲av日韩精品久久久久久密| 久久精品国产亚洲av天美| 国产精品福利在线免费观看| 国产中年淑女户外野战色| 熟女电影av网| 亚洲五月天丁香| 内射极品少妇av片p| 色哟哟·www| 国产 一区 欧美 日韩| 亚洲欧美日韩无卡精品| 很黄的视频免费| 身体一侧抽搐| 亚洲va在线va天堂va国产| 亚洲性久久影院| 国产欧美日韩精品亚洲av| 亚洲第一电影网av| 草草在线视频免费看| 免费看日本二区| 动漫黄色视频在线观看| 嫁个100分男人电影在线观看| 国产精品国产高清国产av| 久久久久久久午夜电影| av黄色大香蕉| 精品国产三级普通话版| 又爽又黄无遮挡网站| 日本黄大片高清| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 一个人免费在线观看电影| 少妇的逼水好多| 国产精品国产三级国产av玫瑰| 中文字幕免费在线视频6| 国产精品一区二区性色av| 97热精品久久久久久| 欧美精品啪啪一区二区三区| 搡老妇女老女人老熟妇| av在线亚洲专区| 亚洲内射少妇av| 真实男女啪啪啪动态图| 亚洲一级一片aⅴ在线观看| 精品午夜福利视频在线观看一区| 日韩欧美精品v在线| 亚洲国产色片| 91久久精品电影网| 在线免费十八禁| 免费av不卡在线播放| 校园人妻丝袜中文字幕| 国产伦一二天堂av在线观看| 欧美一区二区亚洲| 欧美人与善性xxx| 国产一区二区亚洲精品在线观看| 99久久精品一区二区三区| 美女黄网站色视频| 男女那种视频在线观看| 国产精品一及| 在线播放国产精品三级| 亚洲美女视频黄频| 色视频www国产| 久久欧美精品欧美久久欧美| 韩国av在线不卡| 中文字幕av成人在线电影| 日韩亚洲欧美综合| 成人特级黄色片久久久久久久| 中亚洲国语对白在线视频| 老司机深夜福利视频在线观看| 91久久精品国产一区二区三区| 99热这里只有是精品50| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片 | 又爽又黄a免费视频| 一区福利在线观看| 色尼玛亚洲综合影院| 精华霜和精华液先用哪个| 女人被狂操c到高潮| 午夜免费成人在线视频| 三级毛片av免费| 国产女主播在线喷水免费视频网站 | 亚洲av第一区精品v没综合| 亚洲精华国产精华液的使用体验 | 91麻豆精品激情在线观看国产| 欧洲精品卡2卡3卡4卡5卡区| av国产免费在线观看| 国产精品伦人一区二区| 亚洲乱码一区二区免费版| 免费在线观看成人毛片| 村上凉子中文字幕在线| 两人在一起打扑克的视频| 精品人妻1区二区| 可以在线观看的亚洲视频| 国产高清视频在线播放一区| 无人区码免费观看不卡| 国产男靠女视频免费网站| 国产日本99.免费观看| 大型黄色视频在线免费观看| 韩国av在线不卡| 波多野结衣高清无吗| 在线免费十八禁| 91久久精品电影网| 美女 人体艺术 gogo| 日本免费一区二区三区高清不卡| 一级黄片播放器| 神马国产精品三级电影在线观看| 简卡轻食公司| 日日啪夜夜撸| 最新在线观看一区二区三区| 欧美日韩精品成人综合77777| 91麻豆av在线| 美女大奶头视频| 亚洲精品456在线播放app | 99热这里只有是精品50| 床上黄色一级片| 日韩一区二区视频免费看| 色哟哟·www| 亚洲av中文字字幕乱码综合| a级毛片a级免费在线| 欧美高清成人免费视频www| 国产精品自产拍在线观看55亚洲| 成人精品一区二区免费| 国内毛片毛片毛片毛片毛片| 亚洲国产精品合色在线| 一区二区三区免费毛片| 熟女电影av网| 亚洲国产色片| 午夜免费成人在线视频| 99久国产av精品| 午夜福利欧美成人| 欧洲精品卡2卡3卡4卡5卡区| 国产精品国产三级国产av玫瑰| 亚洲在线自拍视频| 免费人成在线观看视频色| 十八禁网站免费在线| 啦啦啦韩国在线观看视频| 天堂动漫精品| 免费人成在线观看视频色| 国国产精品蜜臀av免费| 国产亚洲精品综合一区在线观看| 欧美性感艳星| 国产精品一区二区三区四区免费观看 | 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久久毛片| 日本一二三区视频观看| 人妻夜夜爽99麻豆av| 精品人妻熟女av久视频| 日本与韩国留学比较| 亚洲一区高清亚洲精品| eeuss影院久久| 日本三级黄在线观看| www.色视频.com| 亚洲人成网站高清观看| 变态另类丝袜制服| 尾随美女入室| 欧美成人性av电影在线观看| 少妇高潮的动态图| 午夜福利视频1000在线观看| 免费在线观看影片大全网站| xxxwww97欧美| 久久精品久久久久久噜噜老黄 | 亚洲精品国产成人久久av| 亚洲无线观看免费| 岛国在线免费视频观看| 熟女电影av网| 亚洲人与动物交配视频| 91av网一区二区| 亚洲国产精品久久男人天堂| 国产精品国产三级国产av玫瑰| 中文字幕久久专区| 国产一区二区亚洲精品在线观看| 久久久久久久久久成人| 午夜福利在线观看吧| 亚洲人成网站在线播放欧美日韩| 真人一进一出gif抽搐免费| 亚洲 国产 在线| 啦啦啦啦在线视频资源| 一区福利在线观看| 成人美女网站在线观看视频| 91久久精品国产一区二区三区| 搞女人的毛片| 亚洲精品久久国产高清桃花| ponron亚洲| 国内精品久久久久久久电影| 黄片wwwwww| 久久久久久久久久久丰满 | 在线播放国产精品三级| 久久精品夜夜夜夜夜久久蜜豆| 女的被弄到高潮叫床怎么办 | 麻豆一二三区av精品| 久久久精品欧美日韩精品| 国产精品乱码一区二三区的特点| 久久久精品欧美日韩精品| 国产视频内射| 五月伊人婷婷丁香| 国产高清三级在线| 亚洲精品粉嫩美女一区| 欧美日本视频| 亚洲国产精品成人综合色| 熟女电影av网| 无人区码免费观看不卡| 一a级毛片在线观看| 观看免费一级毛片| 综合色av麻豆| 人妻制服诱惑在线中文字幕| 午夜福利视频1000在线观看| 欧美又色又爽又黄视频| 蜜桃久久精品国产亚洲av| 深夜精品福利| 国产精品人妻久久久久久| 免费观看人在逋| 国产高清不卡午夜福利| av女优亚洲男人天堂| 午夜影院日韩av| 日本三级黄在线观看| 国产男人的电影天堂91| 亚洲,欧美,日韩| 一个人观看的视频www高清免费观看| 免费观看在线日韩| 精品无人区乱码1区二区| 91av网一区二区| 亚洲天堂国产精品一区在线| 国产成人福利小说| 亚洲欧美日韩东京热| 日本一二三区视频观看| 国产av在哪里看| 日韩欧美在线二视频| 国产国拍精品亚洲av在线观看| 成年女人看的毛片在线观看| aaaaa片日本免费| 久久精品夜夜夜夜夜久久蜜豆| 黄色欧美视频在线观看| 欧美性感艳星| 久久香蕉精品热| 少妇人妻精品综合一区二区 | 99久久九九国产精品国产免费| 给我免费播放毛片高清在线观看| 2021天堂中文幕一二区在线观| 中文字幕免费在线视频6| 波多野结衣高清作品| netflix在线观看网站| 男女视频在线观看网站免费| 老熟妇仑乱视频hdxx| 干丝袜人妻中文字幕| 国产免费av片在线观看野外av| 日日摸夜夜添夜夜添小说| 国产高清有码在线观看视频| 三级国产精品欧美在线观看| 久久99热6这里只有精品| 丰满的人妻完整版| 看片在线看免费视频| 69av精品久久久久久| 午夜福利在线观看免费完整高清在 | 精品久久久久久久久av| 国产免费男女视频| 男插女下体视频免费在线播放| 亚洲精品在线观看二区| 99久久成人亚洲精品观看| 九九热线精品视视频播放| 亚洲最大成人中文| 久久久国产成人免费| 国产午夜精品久久久久久一区二区三区 | 午夜精品久久久久久毛片777| 成人三级黄色视频| 人妻丰满熟妇av一区二区三区| 欧美性感艳星| 国产精品久久久久久久电影| 国产亚洲欧美98| 亚洲综合色惰| 日日摸夜夜添夜夜添av毛片 | 国产精品久久久久久亚洲av鲁大| 免费黄网站久久成人精品| 有码 亚洲区| 欧美不卡视频在线免费观看| 日韩 亚洲 欧美在线| 麻豆久久精品国产亚洲av| 国产黄色小视频在线观看| 成年女人看的毛片在线观看| 日本黄色片子视频| 欧美黑人欧美精品刺激| 国产探花极品一区二区| 日韩亚洲欧美综合| 精品久久久久久久久亚洲 | 国产精品女同一区二区软件 | 很黄的视频免费| 99热网站在线观看| 搡老熟女国产l中国老女人| videossex国产| 精品人妻1区二区| 国产综合懂色| 深爱激情五月婷婷| 日本黄大片高清| 亚洲av不卡在线观看| 免费高清视频大片| 伦理电影大哥的女人| 51国产日韩欧美| 国产高清激情床上av| av专区在线播放| 欧美成人a在线观看| 久久精品人妻少妇| 精品人妻一区二区三区麻豆 | 在线观看午夜福利视频| 免费看a级黄色片| 国产又黄又爽又无遮挡在线| 午夜亚洲福利在线播放| 国产女主播在线喷水免费视频网站 | 国产综合懂色| 国产伦人伦偷精品视频| 亚洲中文字幕一区二区三区有码在线看| 一级av片app| 国产精品一区二区三区四区久久| 成人国产综合亚洲| 国产精品美女特级片免费视频播放器| 99久久精品热视频| 日韩在线高清观看一区二区三区 | 国产乱人伦免费视频| x7x7x7水蜜桃| 国产私拍福利视频在线观看| 欧美绝顶高潮抽搐喷水| 成人性生交大片免费视频hd| 国产免费av片在线观看野外av| 尾随美女入室| 国产欧美日韩一区二区精品| 国产主播在线观看一区二区| 中国美白少妇内射xxxbb| 亚洲四区av| 最新中文字幕久久久久| 波多野结衣巨乳人妻| 桃红色精品国产亚洲av| 欧美日韩精品成人综合77777| 麻豆一二三区av精品| a在线观看视频网站| 亚洲自拍偷在线| 窝窝影院91人妻| 在线观看午夜福利视频| 日本色播在线视频| 国产在线男女| 日韩高清综合在线| 精品久久久久久,| 99久久精品一区二区三区| 午夜老司机福利剧场| 亚洲精品日韩av片在线观看| 免费一级毛片在线播放高清视频| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区av在线 | 国产色爽女视频免费观看| 亚洲欧美日韩卡通动漫| 91午夜精品亚洲一区二区三区 | 黄色视频,在线免费观看| 久久精品91蜜桃| 内地一区二区视频在线| 欧美一区二区亚洲| 欧美一区二区精品小视频在线| 中文字幕人妻熟人妻熟丝袜美|