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

    彈載SAR末制導(dǎo)段軌跡控制算法

    2016-09-07 01:09:39徐建城吳華興
    關(guān)鍵詞:水平面制導(dǎo)控制算法

    黃 偉, 徐建城, 吳華興

    (1. 西北工業(yè)大學(xué)電子信息學(xué)院, 陜西 西安 710072;

    2. 空軍工程大學(xué)航空航天工程學(xué)院, 陜西 西安 710038)

    ?

    彈載SAR末制導(dǎo)段軌跡控制算法

    黃偉1, 徐建城1, 吳華興2

    (1. 西北工業(yè)大學(xué)電子信息學(xué)院, 陜西 西安 710072;

    2. 空軍工程大學(xué)航空航天工程學(xué)院, 陜西 西安 710038)

    針對機載合成孔徑雷達(synthetic aperture radar, SAR)制導(dǎo)空地導(dǎo)彈的末制導(dǎo)過程,導(dǎo)彈的運動狀態(tài)如何同時滿足成像分辨率和脫靶量要求的問題,提出一種新的軌跡控制算法。通過分析SAR末制導(dǎo)過程,以及彈目空間幾何與成像分辨率的關(guān)系,將成像分辨率要求轉(zhuǎn)化為對導(dǎo)彈前置角的約束,將脫靶量要求轉(zhuǎn)化為對導(dǎo)彈目標(biāo)線角速度的約束。基于逆最優(yōu)控制,建立含有待定系數(shù)的性能指標(biāo)函數(shù),并根據(jù)給定的系統(tǒng)動態(tài)性能指標(biāo)確定該系數(shù),從而得到具有時變權(quán)值系數(shù)的最優(yōu)控制量,通過權(quán)值系數(shù)的調(diào)節(jié)作用以滿足不同制導(dǎo)模式的指標(biāo)要求。仿真實驗證明了該算法的有效性。

    彈載合成孔徑雷達; 末制導(dǎo); 軌跡控制; 成像分辨率; 脫靶量

    0 引 言

    合成孔徑雷達(synthetic aperture radar, SAR)制導(dǎo)技術(shù)是先進成像制導(dǎo)技術(shù)極為重要的發(fā)展方向[1-3]。通過導(dǎo)引頭的高分辨率成像,可實現(xiàn)對地面或海面目標(biāo)要害部位的精確打擊,由于距離向分辨率主要取決于發(fā)射信號的帶寬,通過導(dǎo)引頭發(fā)射大帶寬信號以獲得高的距離向分辨率,在技術(shù)上已不難實現(xiàn);然而,為了獲得所需的方位向分辨率,制導(dǎo)過程中必須對導(dǎo)彈前置角加以約束;由于SAR導(dǎo)引頭不具備前視成像能力,末制導(dǎo)段需要進行制導(dǎo)模式的切換,切換為前視制導(dǎo)模式后,為保證導(dǎo)彈攻擊地面目標(biāo)的準(zhǔn)確性,必須對目標(biāo)線角速度加以約束。

    由此可知,SAR末制導(dǎo)段的軌跡控制或制導(dǎo)律設(shè)計可描述為多約束條件下的最優(yōu)控制問題。目前,國內(nèi)外對于該問題的研究尚處于起步階段,文獻[4-7]對含有成像約束條件的空地導(dǎo)彈的軌跡優(yōu)化問題進行研究,文獻[8-10]分別應(yīng)用遺傳算法、序列二次規(guī)劃(sequence quadratic program,SQP)算法、Radau偽譜法對SAR制導(dǎo)平臺的導(dǎo)彈飛行軌跡進行了優(yōu)化設(shè)計。然而,以上研究成果主要針對空地SAR制導(dǎo)導(dǎo)彈的中制導(dǎo)段,重點解決導(dǎo)彈的運動狀態(tài)如何滿足成像分辨率的問題,所建立的性能指標(biāo)函數(shù)主要包含成像約束,未涉及末制導(dǎo)段的脫靶量約束。SAR導(dǎo)彈末制導(dǎo)過程中,方位向分辨率越高,需要導(dǎo)彈的前置角越大,而前視制導(dǎo)模式為滿足脫靶量要求,目標(biāo)線角速度應(yīng)趨于零。由此可見,末制導(dǎo)過程中制導(dǎo)模式切換前后對前置角和角速度的控制問題是相互制約的,應(yīng)將末制導(dǎo)過程作為整體考慮。

    基于此,本文針對SAR成像制導(dǎo)空地導(dǎo)彈的末制導(dǎo)過程,將導(dǎo)彈三維運動分解為水平面運動和垂直面運動[11]。在水平面內(nèi)根據(jù)制導(dǎo)模式切換前后的方位向分辨率要求和脫靶量要求,基于逆最優(yōu)控制方法,得到帶有時變權(quán)值系數(shù)的最優(yōu)制導(dǎo)律,通過時變權(quán)值系數(shù)的調(diào)節(jié)作用,使導(dǎo)彈的運動狀態(tài)滿足各制導(dǎo)模式的指標(biāo)要求;在垂直面內(nèi),采用相同的方法,得到滿足脫靶量要求的最優(yōu)制導(dǎo)律,從而建立SAR制導(dǎo)導(dǎo)彈末制導(dǎo)段的最優(yōu)軌跡控制算法,為SAR制導(dǎo)導(dǎo)彈制導(dǎo)系統(tǒng)的總體設(shè)計提供依據(jù)。

    1 空地SAR制導(dǎo)導(dǎo)彈制導(dǎo)過程

    空地SAR制導(dǎo)導(dǎo)彈具備在遠(yuǎn)距離、強地雜波背景條件下探測、識別、跟蹤目標(biāo)及精確打擊目標(biāo)的能力。典型的空地SAR制導(dǎo)導(dǎo)彈制導(dǎo)過程如圖1所示。

    圖1 空地SAR制導(dǎo)過程示意Fig.1 Air-to-ground SAR guidance process

    圖1中,以戰(zhàn)斗機發(fā)射導(dǎo)彈空間位置B點在地面的投影點建立地理坐標(biāo)系OXYZ??盏豐AR制導(dǎo)導(dǎo)彈制導(dǎo)過程可分為如下3個階段:

    (1) 目標(biāo)指示搜索階段(A點至B點):戰(zhàn)斗機任務(wù)系統(tǒng)向?qū)椦b訂目標(biāo)參數(shù),SAR導(dǎo)引頭對斜前方目標(biāo)區(qū)域進行搜索成像。導(dǎo)引頭截獲目標(biāo)后,飛行員手動或自動發(fā)射導(dǎo)彈。

    (2) 聚束成像中制導(dǎo)階段(B點至C點): SAR導(dǎo)引頭工作于聚束成像模式,應(yīng)用具有運動補償?shù)某上袼惴╗12-15]對目標(biāo)區(qū)域成像,將目標(biāo)圖像與預(yù)先存儲的參考圖像相比較,結(jié)合彈載傳感器量測的導(dǎo)彈飛行參數(shù),獲取彈目相對運動狀態(tài)信息。該階段導(dǎo)彈通常位于預(yù)定高度飛行,導(dǎo)彈的運動狀態(tài)應(yīng)滿足SAR成像分辨率要求。

    (3) 聚束成像/前視攻擊末制導(dǎo)階段(C點至D點):隨著彈目距離的不斷減小,當(dāng)彈目距離滿足前視成像分辨率要求時,為實現(xiàn)對目標(biāo)的前視攻擊,導(dǎo)引頭工作模式將由聚束成像切換為前視成像[3]。在此期間,導(dǎo)彈仍保持在預(yù)定高度飛行,直到高低角滿足一定條件時,導(dǎo)彈以俯沖運動狀態(tài)接近目標(biāo)。

    2 聚束式SAR成像分辨率分析

    聚束式SAR導(dǎo)引頭發(fā)射的線性調(diào)頻信號歸一化形式為

    (1)

    式中,Tp為導(dǎo)彈發(fā)射脈沖的時間寬度;τ=t-nTr為脈內(nèi)時間,Tr為脈沖重復(fù)周期,n為脈沖數(shù);fc為載波頻率;γ為信號線性調(diào)頻率;rect()為矩形函數(shù)。導(dǎo)引頭接收的回波信號經(jīng)檢波去除高頻項[16]可得

    (2)

    式中,σn為目標(biāo)散射系數(shù);Ts為合成孔徑時間;td=2R(t)/c表示回波信號時延,R(t)表示彈目距離;λ為雷達波長。rect(t/Ts)e-j4πR(t)/λ項反映了導(dǎo)彈相對目標(biāo)的距離史,決定聚束式SAR方位向分辨率。

    圖1中,設(shè)導(dǎo)彈速度矢量為V,合成孔徑時間內(nèi)導(dǎo)彈質(zhì)心與目標(biāo)的瞬時距離為R,V相對于目標(biāo)線Los的前置角為φ′??紤]到SAR導(dǎo)引頭工作于小合成孔徑角情況,由此可得合成孔徑時間內(nèi)多普勒帶寬為

    (3)

    方位向分辨率[16]為

    (4)

    分析式(4)可知:①聚束式SAR方位向分辨率與合成孔徑時間Ts、雷達波長λ以及導(dǎo)彈相對目標(biāo)的運動參數(shù)有關(guān);②導(dǎo)引頭可通過增加Ts和前置角φ′,有效提高方位向分辨率。然而,為滿足成像制導(dǎo)的實時性要求,需進行若干次成像,Ts通常取固定值并且不宜過大,并且,增大導(dǎo)彈前置角將導(dǎo)致目標(biāo)線角速度增大,隨之脫靶量增大。由此,在導(dǎo)引頭所需分辨率ρT和合成孔徑時間Ts給定的前提下,綜合考慮方位向分辨率及脫靶量要求,導(dǎo)彈最小前置角為

    (5)

    3 SAR末制導(dǎo)段軌跡控制算法

    如圖2所示,導(dǎo)彈的三維運動可由正交平面運動的組合表示[11]。由于SAR末制導(dǎo)階段,導(dǎo)彈在水平面的運動狀態(tài)是確保成像分辨率和減小脫靶量的關(guān)鍵,所以本文重點研究SAR制導(dǎo)導(dǎo)彈在水平面的末制導(dǎo)軌跡控制算法,該算法以導(dǎo)彈相對目標(biāo)的質(zhì)點運動學(xué)方程為基礎(chǔ),根據(jù)末制導(dǎo)段各項指標(biāo)要求,建立具有待定權(quán)值系數(shù)的性能指標(biāo)函數(shù),由系統(tǒng)動態(tài)調(diào)節(jié)時間及相關(guān)約束條件計算該系數(shù),逆向得出具有時變權(quán)值系數(shù)的最優(yōu)制導(dǎo)律;在此基礎(chǔ)上,依據(jù)水平軌跡控制算法的建立方法,得到了垂直平面的末制導(dǎo)最優(yōu)制導(dǎo)律。

    3.1水平面內(nèi)導(dǎo)彈與目標(biāo)相對運動關(guān)系建模

    空地SAR制導(dǎo)導(dǎo)彈主要攻擊地面固定或慢速移動的目標(biāo),由此可忽略目標(biāo)運動的影響,并假設(shè)導(dǎo)彈為質(zhì)點。在圖1的基礎(chǔ)上,可得導(dǎo)彈與地面目標(biāo)的相對運動關(guān)系,如圖2所示。

    圖2 導(dǎo)彈與目標(biāo)相對運動關(guān)系Fig.2 Relative motion between missile and the target

    SAR成像制導(dǎo)過程中,導(dǎo)彈位于預(yù)定高度的水平面內(nèi),V為導(dǎo)彈速度矢量;H為導(dǎo)彈飛行高度;R為彈目距離;Rh為彈目水平距離;φ、ε、θ分別為水平面內(nèi)導(dǎo)彈前置角、目標(biāo)線角及偏航角。速度矢量V在地面上投影為Vh,將Vh投影到水平面目標(biāo)線Losh及其法線上,可得

    (6)

    (7)

    (8)

    將式(6)代入式(8)可得

    (9)

    (10)

    由式(9)和式(10)可得

    (11)

    式(11)即為水平面內(nèi)導(dǎo)彈與目標(biāo)相對運動關(guān)系的數(shù)學(xué)模型。

    3.2水平面軌跡控制算法的建立

    (12)

    式(12)可轉(zhuǎn)化為標(biāo)準(zhǔn)形式

    (13)

    建立性能指標(biāo)函數(shù)基于如下考慮:①為滿足聚束成像段方位向分辨率要求,導(dǎo)彈前置角應(yīng)趨于理想前置角;②為滿足前視攻擊段脫靶量要求,目標(biāo)線角速度ωh應(yīng)趨于零;③應(yīng)考慮制導(dǎo)過程的經(jīng)濟性。由此,二次型性能指標(biāo)函數(shù)為

    (14)

    (15)

    P(t)滿足黎卡提方程及邊值條件:

    (16)

    通常情況下,黎卡提方程的求解是十分困難的,加之制導(dǎo)過程中導(dǎo)彈與目標(biāo)運動狀態(tài)的劇烈變化,結(jié)束控制時刻tk難以準(zhǔn)確估計,這將進一步增加黎卡提方程求解的復(fù)雜性。針對這一問題,有關(guān)學(xué)者在進行大量的裝備驗證之后,在文獻[17]中指出:由于制導(dǎo)過程的每一時刻t都可能對應(yīng)結(jié)束時刻tk,即t→tk,在制導(dǎo)律設(shè)計過程中,尤其針對導(dǎo)彈末制導(dǎo)段可直接應(yīng)用P(t)的穩(wěn)態(tài)解[17]

    (17)

    由式(12)~式(17),取理想目標(biāo)線角速度ωhT=0,可得最優(yōu)控制量為

    (18)

    式中,q11、q22、k1均為性能指標(biāo)函數(shù)中待定的時變系數(shù),假設(shè)比值q11/k1和q22/k1的絕對數(shù)值相當(dāng),分析式(18)可知:①SAR成像制導(dǎo)模式下,彈目水平距離Rh較大,目標(biāo)線角速度較小,q11/k1Vh相對q22/k1Rh較大。此時,控制量jh主要取決于水平前置角誤差φ-φT,由于控制量jh的作用,使實際前置角φ趨于φT,滿足方位向分辨率要求;②前視制導(dǎo)模式下,隨著彈目水平距離Rh不斷減小,q22/k1Rh不斷增大,控制量jh主要取決于目標(biāo)線角速度ωh,在控制量的作用下ωh逐步趨于零,從而有效減小導(dǎo)彈脫靶量;③由于時變權(quán)值系數(shù)的調(diào)節(jié)作用,可實現(xiàn)制導(dǎo)模式切換前后控制量的平滑過渡。

    為實現(xiàn)該算法,應(yīng)根據(jù)條件計算水平面理想前置角φT,以及比值系數(shù)q11/k1和q22/k1。

    3.3水平面理想前置角φT的計算

    理想前置角即為滿足所需成像分辨率的導(dǎo)彈最小前置角。如圖2所示,前置角φ′與水平面前置角φ之間的幾何關(guān)系式為

    (19)

    由式(19)和式(5)可得

    (20)

    式(20)即為Ts和ρT給定的前提下,導(dǎo)彈水平面理想前置角表達式。φT與彈目水平距離Rh、導(dǎo)彈飛行高度H、合成孔徑時間Ts、彈目水平接近速度Vh、導(dǎo)引頭所需分辨率ρT及雷達波長λ有關(guān)。φT隨導(dǎo)彈與目標(biāo)的接近而不斷減小。

    3.4比值系數(shù)q11/k1和q22/k1的計算

    比值系數(shù)q11/k1和q22/k1的計算應(yīng)著重考慮導(dǎo)彈運動的限制條件以及系統(tǒng)動態(tài)調(diào)節(jié)時間兩方面因素。令λ1=q11/k1及λ2=q22/k1,計算λ1和λ2時應(yīng)考慮的限制條件為

    (1) 彈目最小接近速度Vhmin;

    (2) 彈目最小距離Rhmin;

    (3) 目標(biāo)線角速度最大誤差Δωhmax;

    (4) 前置角最大誤差Δφmax;

    (5) 導(dǎo)彈目標(biāo)線法向加速度最大值jhmax。

    上述限制條件滿足關(guān)系式

    (21)

    假設(shè)導(dǎo)彈控制系統(tǒng)無遲滯的響應(yīng)jh,將具有待定系數(shù)的控制量算法式(18)代入式(12),可得系統(tǒng)狀態(tài)方程為

    (22)

    根據(jù)制導(dǎo)系統(tǒng)趨于穩(wěn)定的動態(tài)調(diào)節(jié)時間TLim,結(jié)合式(21)和式(22),可計算得出q11/k1及q22/k1。限于篇幅,直接給出結(jié)果

    (23)

    (24)

    將式(23)和式(24)代入算法式(18),即得到具有時變權(quán)值系數(shù)q11/k1Vh和q22/k1Ra的最優(yōu)控制量,正是由于控制量時變權(quán)值系數(shù)的調(diào)節(jié)作用,使之在不同制導(dǎo)模式下各誤差參數(shù)對控制量的影響不同,從聚束成像段前置角誤差控制量占優(yōu),轉(zhuǎn)化為前視攻擊段目標(biāo)線角速度誤差控制量占優(yōu),滿足不同制導(dǎo)模式的指標(biāo)要求。應(yīng)該注意到,系數(shù)q11、q22、k1同時包含于控制量和性能指標(biāo)函數(shù)中,根據(jù)系統(tǒng)動態(tài)性能指標(biāo)計算得到的q11/k1和q22/k1不會影響控制量的最優(yōu)性。

    3.5三維軌跡控制算法的建立

    SAR成像制導(dǎo)模式下,導(dǎo)彈處于預(yù)定高度飛行,切換為前視制導(dǎo)模式后,當(dāng)高低角滿足β≥β1時,導(dǎo)彈以俯沖運動狀態(tài)接近目標(biāo),式(18)建立了水平面內(nèi)的SAR末制導(dǎo)階段的軌跡控制算法,下面主要研究導(dǎo)彈在垂直平面內(nèi)的軌跡控制算法。應(yīng)用本文提出的逆最優(yōu)控制方法,首先建立垂直平面內(nèi)導(dǎo)彈與目標(biāo)的運動學(xué)方程

    (25)

    (26)

    式中,垂直平面內(nèi)理想目標(biāo)線角速度ωvT=0;q3、k2均為待定的時變系數(shù)。根據(jù)式(15)和式(16),可得最優(yōu)控制量為

    (27)

    式(27)中,jv與系數(shù)的比值q3/k2有關(guān),令k2=1。假設(shè)導(dǎo)彈控制系統(tǒng)無遲滯的響應(yīng)jv,將式(27)代入式(25),可得閉環(huán)系統(tǒng)狀態(tài)方程為

    (28)

    (29)

    由式(29)可得

    (30)

    (31)

    將式(31)給出的q3代入式(27)可得

    (32)

    綜上可得SAR制導(dǎo)導(dǎo)彈末制導(dǎo)過程中水平面和垂直面的最優(yōu)制導(dǎo)律為

    (33)

    式(33)中,權(quán)值系數(shù)q11/k1和q22/k1由式(23)和式(24)給出;前置角φT由式(20)給出;jh和jv可轉(zhuǎn)換為圖2地理坐標(biāo)系中3個坐標(biāo)軸向的指令加速度[11]

    (34)

    將式(34)給出的指令加速度由地理坐標(biāo)系轉(zhuǎn)換到彈體坐標(biāo)系,即可得到由導(dǎo)彈制導(dǎo)系統(tǒng)輸出的3個軸向的指令加速度。

    4 仿真實例

    下面對算法式(33)進行仿真驗證。本文所建立的軌跡控制算法主要解決導(dǎo)彈的運動狀態(tài)如何滿足SAR成像分辨率約束和脫靶量約束的問題,適用于導(dǎo)彈制導(dǎo)系統(tǒng)的總體設(shè)計階段,在仿真過程中,假設(shè)導(dǎo)彈為質(zhì)點,導(dǎo)彈的控制系統(tǒng)是無遲滯響應(yīng)的非慣性系統(tǒng),控制量為算法式(34)給出的指令加速度。需要指出,SAR制導(dǎo)空地導(dǎo)彈在水平面的運動狀態(tài)是確保成像分辨率和減小脫靶量的關(guān)鍵,也是仿真驗證的重點內(nèi)容。設(shè)導(dǎo)彈初始位置坐標(biāo)為(0 km,10 km,3 km),初始速度Vx=200 m/s,Vy=600 m/s,Vz=0 m/s。目標(biāo)中心位置坐標(biāo)為(40 km,40 km,0 km),其他參數(shù)如表1所示。

    表1 仿真參數(shù)

    圖3 三維空間導(dǎo)彈軌跡Fig.3 Missile trajectory in 3D space

    為了更加直觀地反映導(dǎo)彈水平面的運動,圖4給出了導(dǎo)彈末制導(dǎo)過程中水平面的運動軌跡,由圖可知導(dǎo)彈彈道平滑。圖5對應(yīng)水平面導(dǎo)彈的目標(biāo)線法向過載,其過載在允許范圍之內(nèi)。

    圖4 水平面導(dǎo)彈軌跡Fig.4 Missile trajectory in horizontal surface

    圖5 水平面導(dǎo)彈過載Fig.5 Missile overload in horizontal surface

    圖6給出了導(dǎo)彈實際前置角與理想前置角的變化規(guī)律。在SAR成像制導(dǎo)階段,實際前置角趨近于理想前置角,滿足方位向分辨率要求,切換為前視制導(dǎo)模式后,實際前置角逐步趨近于零,從而實現(xiàn)對目標(biāo)的前視打擊。

    圖7給出了算法式(18)所包含的權(quán)值系數(shù)隨制導(dǎo)時間的變化情況:前置角誤差權(quán)值系數(shù)在SAR成像段相對較大,并隨制導(dǎo)時間逐漸減小;目標(biāo)線角速度權(quán)值系數(shù)隨制導(dǎo)時間逐漸增大,表明了權(quán)值系數(shù)對各誤差控制量的調(diào)節(jié)規(guī)律。

    圖6 前置角變化規(guī)律Fig.6 Varying law of lead angle

    圖7 權(quán)值系數(shù)變化規(guī)律Fig.7 Varying law of weight coefficients

    通過分析末制導(dǎo)軌跡控制算法仿真結(jié)果可知:①算法在SAR成像制導(dǎo)模式下滿足導(dǎo)引頭方位向分辨率要求;②算法在前視制導(dǎo)模式下滿足導(dǎo)彈脫靶量要求;③算法綜合考慮了成像分辨率和脫靶量要求,通過時變權(quán)值系數(shù)的調(diào)節(jié)作用,使導(dǎo)彈過載變化平穩(wěn),適用于SAR末制導(dǎo)段的軌跡控制。

    5 結(jié) 論

    本文提出的控制算法有效解決了對于SAR制導(dǎo)導(dǎo)彈而言較為特殊的末制導(dǎo)軌跡控制問題,研究總結(jié)如下:①算法綜合考慮了成像分辨率要求和脫靶量要求,通過時變權(quán)值系數(shù)對控制量的調(diào)節(jié)作用,使該算法在不同制導(dǎo)模式下,滿足不同指標(biāo)要求;②基于逆最優(yōu)控制,算法充分考慮了閉環(huán)制導(dǎo)回路的動態(tài)特性,得到了SAR末制導(dǎo)段三維最優(yōu)制導(dǎo)律;③算法適用于多種包含SAR復(fù)合制導(dǎo)模式導(dǎo)彈的制導(dǎo)系統(tǒng)總體設(shè)計,滿足彈載計算的實時性要求,易于工程實現(xiàn)。

    [1]QinYL,WangJT,WangHQ,etal.Overviewonmissile-bornesyntheticapertureradar[J].Signal Processing, 2009, 25(4):630-635.(秦玉亮,王建濤,王宏強,等.彈載合成孔徑雷達技術(shù)研究綜述[J].信號處理,2009,25(4):630-635.)

    [2]ZhangG,ZhuMB,ZhaoZB,etal.Studyonapplicationmodeandkeytechnologyofmissile-bornesyntheticapertureradar[J].Cruise Missile, 2011(9): 67-73. (張剛, 祝明波, 趙振波,等. 彈載SAR應(yīng)用模式及關(guān)鍵技術(shù)探討[J].飛航導(dǎo)彈, 2011(9): 67-73.)

    [3]ZhuMB,YangLB,YangRL. Missile-borne SAR guidance: fundamental and key techniques[M].Beijing:NationalDefenseIndustryPress, 2014: 33-46. (祝明波, 楊立波, 楊汝良. 彈載合成孔徑雷達制導(dǎo)及其關(guān)鍵技術(shù)[M]. 北京:國防工業(yè)出版社, 2014: 33-46.)

    [4]FarooqA,LimebeerDJN.Trajectoryoptimizationforair-to-surfacemissileswithimagingradar[J].Journal of Guidance, Control, and Dynamics, 2002, 25(5): 876-887.

    [5]FarooqA,LimebeerDJN.Optimaltrajectoryregulationforradarimagingguidance[J].Journal of Guidance,Control,and Dynamics, 2008, 31(4): 1076-1092.

    [6]HodgsonJA,LeeDW.TerminalguidanceusingaDopplerbeamsharpeningradar[C]∥Proc.of the AIAA Guidance, Navigation and Control Conference and Exhibit, 2003: 1-11.

    [7]HodgsonJA.Trajectoryoptimazationusingdifferentialinclusiontominimizeuncertaintyintargetlocationestimation[C]∥Proc.of the AIAA Guidance, Navigation and Control Conference and Exhibit, 2005: 1-17.

    [8]XieHY,FanHQ,ZhaoHZ,etal.TrajectorydesignandoptimizationforaSARseeker[J].Systems Engineering and Electronics, 2010, 32(2): 332-337. (謝華英, 范紅旗, 趙宏鐘, 等.SAR成像導(dǎo)引頭的彈道設(shè)計與優(yōu)化[J].系統(tǒng)工程與電子技術(shù), 2010, 32(2): 332-337.)

    [9]ZhangG,ZhuMB,ChenY,etal. 3DtrajectoryoptimizationoftheSARimagingseekerbasedanSQP[J].Computer Engine-ering & Science,2012,34(4):146-150.(張剛, 祝明波, 陳瑜,等. 基于SQP算法的SAR成像導(dǎo)引頭三維彈道優(yōu)化[J].計算機工程與科學(xué), 2012, 34(4): 146-150.)

    [10]WangJ,ZhangG,YangLB,etal.TrajectoryoptimizationofSARimagingguidancebasedanRadaupseudospectralmethod[J].Modern Defence Technology,2013,41(5):41-45.(王杰,張剛,楊立波,等. 基于Radau偽譜法的SAR末制導(dǎo)成像段彈道優(yōu)化[J].現(xiàn)代防御技術(shù), 2013, 41(5): 41-45.)

    [11]XueLH,FanY,SongC. Modern missile guidance[M].Beijing:NationalDefenseIndustryPress,2013:10-40. (薛麗華,范宇,宋闖.現(xiàn)代導(dǎo)彈制導(dǎo)[M].北京:國防工業(yè)出版社,2013:10-40.)

    [12]XingM,JiangX,WuR,etal.MotioncompensationforUAVSARbasedonrawradardata[J].IEEE Trans. on Geoscience and Remote Sensing, 2009, 47(8): 2870-2883.

    [13]FornaroG.TrajectorydeviationinairborneSAR:analysisandcompensation[J].IEEE Trans.on Aerospace and Electronic Systems, 1999, 35(4): 997-1009.

    [14]FornaroG,FranceschettiG,PernaS.Motioncompensationerrors:effectsontheaccuracyofairborneSARimages[J].IEEE Trans. on Aerospace and Electronic Systems, 2005, 41(4): 1338-1352.

    [15]MacedoKAC,ScheiberR.Precisetopographyandaperture-dependentmotioncompensationforairborneSAR[J].IEEE Geoscience and Remote Sensing Letters,2005,2(2):172-176.

    [16]YangRL. High resolution microware imaging[M].Beijing:NationalDefenseIndustryPress,2013:18-81. (楊汝良. 高分辨率微波成像[M]. 北京:國防工業(yè)出版社, 2013: 18-81.)

    [17]KhannaSK,MailerK. Aeronautical radio control system[M].Moscow:RadioTechnologyPress, 2003: 65-81.

    [18]ShenMH,LuoJX,YuL,etal.Acombinedguidancemethodfortheair-to-groundmissile[J].Systems Engineering and Electronics, 2003, 25(4): 469-471. (沈明輝, 羅繼勛, 于雷, 等. 空地導(dǎo)彈的一種復(fù)合引導(dǎo)法[J].系統(tǒng)工程與電子技術(shù), 2003, 25(4): 469-471.)

    Trajectory control arithmetic of missile-borne SAR terminal guidance

    HUANG Wei1, XU Jian-cheng1, WU Hua-xing2

    (1.ElectronicandInformationCollege,NorthwesternPolytechnicUniversity,Xi’an710072,China;2.AeronauticsandAstronauticsEngineeringCollege,AirForceEngineeringUniversity,Xi’an710038,China)

    For the process of airborne synthetic aperture radar (SAR) missile terminal guidance, a new tra-jectory control arithmetic is proposed to meet the requirements of imaging resolution and the missile miss distance for motion states of the missile. By analyzing the SAR terminal guiding process and the relation between missile-target space geometry and imaging resolution, the need for imaging resolution is converted to constraints for motion states of the missile, and the requirement of miss distance is converted to constraints for lead angle of the missile. Based on inverse optimal control, performance index function containing undetermined coefficients is build and the coefficients is calculated on the basis of dynamic performance in the control system, then optimal control variable is calculated with time-variant weight coefficients, which can be tuned to meet the index requirement for different guidance modes. Simulation results demonstrate the validity of the arithmetic.

    missile-borne synthetic aperture radar (SAR); terminal guidance; trajectory control; imaging resolution; miss distance

    2015-06-24;

    2015-12-18;網(wǎng)絡(luò)優(yōu)先出版日期:2016-03-03。

    國家自然科學(xué)基金(61472441)資助課題

    V 249; TJ 765.3

    A

    10.3969/j.issn.1001-506X.2016.09.21

    黃偉(1980-),男,副教授,博士研究生,主要研究方向為SAR成像制導(dǎo)技術(shù)。

    E-mail:huangwei800519@163.com

    徐建城(1957-),男,教授,博士研究生導(dǎo)師,博士,主要研究方向為無線傳感器網(wǎng)絡(luò)。

    E-mail:xujchg@nwpu.edu.cn

    吳華興(1978-),男,講師,博士研究生,主要研究方向為作戰(zhàn)建模與仿真。

    E-mail:dseyzmz@126.com

    網(wǎng)絡(luò)優(yōu)先出版地址:http://www.cnki.net/kcms/detail/11.2422.TN.20160303.1521.002.html

    猜你喜歡
    水平面制導(dǎo)控制算法
    基于ARM+FPGA的模塊化同步控制算法研究
    坡角多大,圓柱體在水平面滾得最遠(yuǎn)
    基于MPSC和CPN制導(dǎo)方法的協(xié)同制導(dǎo)律
    基于在線軌跡迭代的自適應(yīng)再入制導(dǎo)
    水平面內(nèi)勻速圓周運動的四種模型解讀
    帶有攻擊角約束的無抖振滑模制導(dǎo)律設(shè)計
    一種優(yōu)化的基于ARM Cortex-M3電池組均衡控制算法應(yīng)用
    水平面上圓周運動中臨界問題的分析和解題策略
    復(fù)合制導(dǎo)方式確保精確入軌
    太空探索(2014年1期)2014-07-10 13:41:49
    一種非圓旋轉(zhuǎn)工件支撐裝置控制算法
    80岁老熟妇乱子伦牲交| 最黄视频免费看| 婷婷成人精品国产| 美国免费a级毛片| 免费在线观看黄色视频的| 国产成人精品久久久久久| 亚洲内射少妇av| 精品人妻偷拍中文字幕| 纯流量卡能插随身wifi吗| 美女视频免费永久观看网站| 久久久久网色| 久久久精品94久久精品| 免费在线观看黄色视频的| 亚洲色图 男人天堂 中文字幕 | 欧美+日韩+精品| 宅男免费午夜| 国产成人午夜福利电影在线观看| 又黄又爽又刺激的免费视频.| 最近最新中文字幕大全免费视频 | 免费大片黄手机在线观看| 另类亚洲欧美激情| 久久这里有精品视频免费| 午夜精品国产一区二区电影| 热re99久久精品国产66热6| 少妇熟女欧美另类| 国产精品国产三级国产专区5o| a级毛色黄片| 内地一区二区视频在线| 一级爰片在线观看| 考比视频在线观看| 高清视频免费观看一区二区| 捣出白浆h1v1| 国产在线视频一区二区| 成人亚洲精品一区在线观看| 亚洲欧美一区二区三区黑人 | 国产探花极品一区二区| 亚洲国产欧美日韩在线播放| 免费大片18禁| 一级毛片 在线播放| 精品国产露脸久久av麻豆| 夜夜爽夜夜爽视频| 欧美成人午夜精品| 亚洲,一卡二卡三卡| 成人国语在线视频| 大香蕉久久成人网| 看十八女毛片水多多多| 国产又爽黄色视频| 最近的中文字幕免费完整| www日本在线高清视频| 大香蕉97超碰在线| 亚洲精品456在线播放app| 国产淫语在线视频| 国产有黄有色有爽视频| 中国三级夫妇交换| 啦啦啦中文免费视频观看日本| 9热在线视频观看99| 丝瓜视频免费看黄片| 在线精品无人区一区二区三| 街头女战士在线观看网站| 狂野欧美激情性xxxx在线观看| 又黄又粗又硬又大视频| 乱人伦中国视频| 午夜免费观看性视频| 九九在线视频观看精品| 成人黄色视频免费在线看| 黄片播放在线免费| 免费黄色在线免费观看| 久久精品aⅴ一区二区三区四区 | 亚洲熟女精品中文字幕| 国产综合精华液| 欧美成人午夜精品| 欧美老熟妇乱子伦牲交| 亚洲色图 男人天堂 中文字幕 | 日韩精品免费视频一区二区三区 | 精品久久蜜臀av无| 9191精品国产免费久久| 亚洲国产成人一精品久久久| 国产在线免费精品| 波多野结衣一区麻豆| 免费av中文字幕在线| 国产亚洲午夜精品一区二区久久| 少妇的逼好多水| 看免费av毛片| 欧美 亚洲 国产 日韩一| 另类精品久久| 丁香六月天网| 午夜视频国产福利| 女人久久www免费人成看片| 性高湖久久久久久久久免费观看| 色94色欧美一区二区| 天天操日日干夜夜撸| 婷婷色麻豆天堂久久| 国产熟女午夜一区二区三区| 精品久久久久久电影网| 国产国语露脸激情在线看| 国产一区二区三区av在线| 五月玫瑰六月丁香| 男女边吃奶边做爰视频| 午夜久久久在线观看| 亚洲欧洲日产国产| 9色porny在线观看| 色哟哟·www| 国产精品无大码| 亚洲精品视频女| 久久久久久久久久人人人人人人| 国产免费一级a男人的天堂| 少妇被粗大猛烈的视频| 午夜激情久久久久久久| 国精品久久久久久国模美| 2021少妇久久久久久久久久久| 国产一级毛片在线| 久久狼人影院| 色网站视频免费| 成人手机av| 又粗又硬又长又爽又黄的视频| 在线观看美女被高潮喷水网站| 人妻少妇偷人精品九色| 国产日韩欧美在线精品| 成人免费观看视频高清| 一级毛片黄色毛片免费观看视频| 制服诱惑二区| 极品人妻少妇av视频| 高清av免费在线| 91精品国产国语对白视频| 一二三四在线观看免费中文在 | 久久久久久久久久人人人人人人| 制服诱惑二区| 国产视频首页在线观看| 水蜜桃什么品种好| 秋霞在线观看毛片| 韩国精品一区二区三区 | 精品久久蜜臀av无| 亚洲国产精品一区三区| 国产av国产精品国产| 高清不卡的av网站| 色婷婷久久久亚洲欧美| 国产精品99久久99久久久不卡 | 国产亚洲午夜精品一区二区久久| 青青草视频在线视频观看| 一本大道久久a久久精品| 美女国产高潮福利片在线看| 又黄又爽又刺激的免费视频.| 香蕉丝袜av| 亚洲av成人精品一二三区| 欧美精品国产亚洲| 狂野欧美激情性xxxx在线观看| 亚洲欧美日韩卡通动漫| 9色porny在线观看| 最近中文字幕高清免费大全6| 欧美日韩国产mv在线观看视频| 少妇人妻精品综合一区二区| 在线亚洲精品国产二区图片欧美| 国产不卡av网站在线观看| 日本黄大片高清| 亚洲欧洲国产日韩| 制服人妻中文乱码| 老司机亚洲免费影院| 妹子高潮喷水视频| 在线观看一区二区三区激情| 免费av不卡在线播放| 免费观看a级毛片全部| www.熟女人妻精品国产 | 永久网站在线| 日日撸夜夜添| 亚洲欧美清纯卡通| 乱码一卡2卡4卡精品| www.av在线官网国产| 香蕉国产在线看| 青春草亚洲视频在线观看| 精品一区二区三区视频在线| 国产又爽黄色视频| av.在线天堂| 成人国产麻豆网| 少妇的丰满在线观看| 国产精品欧美亚洲77777| 亚洲av免费高清在线观看| 99久久人妻综合| 国产精品国产三级国产专区5o| 久久久国产欧美日韩av| 国产精品.久久久| 男的添女的下面高潮视频| 亚洲精品乱久久久久久| 国产精品久久久久久av不卡| 蜜桃国产av成人99| 两个人免费观看高清视频| 中文字幕精品免费在线观看视频 | 亚洲高清免费不卡视频| 日韩电影二区| 日产精品乱码卡一卡2卡三| 精品一区在线观看国产| 国产日韩欧美亚洲二区| 亚洲熟女精品中文字幕| 国产精品99久久99久久久不卡 | 男人添女人高潮全过程视频| 五月伊人婷婷丁香| 亚洲精品第二区| 国产成人一区二区在线| 亚洲国产精品专区欧美| 天堂俺去俺来也www色官网| 国产精品一区二区在线观看99| 亚洲av.av天堂| 亚洲av免费高清在线观看| 国产伦理片在线播放av一区| 日韩三级伦理在线观看| 国产爽快片一区二区三区| 高清黄色对白视频在线免费看| 国产综合精华液| 乱码一卡2卡4卡精品| 亚洲精品视频女| 亚洲成人手机| 人妻一区二区av| 久久久久久伊人网av| 亚洲精品第二区| 久久久欧美国产精品| 99精国产麻豆久久婷婷| 国产黄频视频在线观看| 亚洲精品色激情综合| 精品午夜福利在线看| 日本爱情动作片www.在线观看| 久久久久精品性色| 亚洲欧美色中文字幕在线| 成人亚洲欧美一区二区av| 欧美精品一区二区免费开放| 一级a做视频免费观看| 免费黄色在线免费观看| 80岁老熟妇乱子伦牲交| 亚洲精品自拍成人| 亚洲国产精品一区三区| 欧美3d第一页| 丰满饥渴人妻一区二区三| 国产免费又黄又爽又色| 日本黄大片高清| 2018国产大陆天天弄谢| 亚洲国产精品成人久久小说| 日韩视频在线欧美| 男女免费视频国产| xxx大片免费视频| 91精品国产国语对白视频| 超碰97精品在线观看| 亚洲中文av在线| 久久人人爽人人爽人人片va| 乱码一卡2卡4卡精品| 国产精品 国内视频| 大片免费播放器 马上看| 日韩欧美一区视频在线观看| 日韩中字成人| 男人爽女人下面视频在线观看| 在线天堂中文资源库| 国国产精品蜜臀av免费| 欧美精品高潮呻吟av久久| 免费观看a级毛片全部| 久久久久精品性色| 人体艺术视频欧美日本| 91精品国产国语对白视频| 国产成人午夜福利电影在线观看| 国产亚洲欧美精品永久| 十八禁高潮呻吟视频| 男女边摸边吃奶| videossex国产| 欧美激情国产日韩精品一区| 久久ye,这里只有精品| 亚洲图色成人| 久久精品国产自在天天线| 亚洲人与动物交配视频| 男女下面插进去视频免费观看 | 老司机影院毛片| 国产精品欧美亚洲77777| 内地一区二区视频在线| 欧美精品一区二区免费开放| 国产日韩欧美视频二区| 一级,二级,三级黄色视频| 亚洲精品国产av成人精品| 亚洲精品美女久久av网站| 日韩av免费高清视频| 久久精品国产亚洲av天美| 只有这里有精品99| 亚洲五月色婷婷综合| 你懂的网址亚洲精品在线观看| 人人澡人人妻人| h视频一区二区三区| 日本色播在线视频| 国产69精品久久久久777片| 成人毛片60女人毛片免费| 9191精品国产免费久久| videosex国产| 亚洲,一卡二卡三卡| 18禁国产床啪视频网站| 激情五月婷婷亚洲| 丝袜喷水一区| 欧美人与性动交α欧美精品济南到 | 精品一区二区免费观看| 久久久精品区二区三区| 久久99精品国语久久久| a 毛片基地| 欧美日韩综合久久久久久| 激情五月婷婷亚洲| 男女边吃奶边做爰视频| 免费女性裸体啪啪无遮挡网站| 成人手机av| 青春草亚洲视频在线观看| 熟女人妻精品中文字幕| 亚洲av在线观看美女高潮| 免费观看在线日韩| 青春草亚洲视频在线观看| 国产精品熟女久久久久浪| 成年女人在线观看亚洲视频| 亚洲性久久影院| 少妇 在线观看| 九九爱精品视频在线观看| www.av在线官网国产| 国产永久视频网站| 大码成人一级视频| 亚洲欧美中文字幕日韩二区| 伦精品一区二区三区| 男女啪啪激烈高潮av片| 久久久久人妻精品一区果冻| 一区二区av电影网| 成年人免费黄色播放视频| 免费在线观看黄色视频的| 国产国拍精品亚洲av在线观看| 久久精品久久久久久久性| 免费av不卡在线播放| 有码 亚洲区| 丝袜在线中文字幕| 久久鲁丝午夜福利片| 色婷婷av一区二区三区视频| 亚洲国产色片| 久久久久久久久久成人| 我的女老师完整版在线观看| 国产免费一区二区三区四区乱码| 日韩精品免费视频一区二区三区 | 99久久中文字幕三级久久日本| 在现免费观看毛片| 人人妻人人爽人人添夜夜欢视频| 国产一区二区在线观看av| www日本在线高清视频| 国产探花极品一区二区| 久久青草综合色| 日韩成人伦理影院| 美女xxoo啪啪120秒动态图| 日韩欧美一区视频在线观看| 男女国产视频网站| 日韩精品免费视频一区二区三区 | 黑人巨大精品欧美一区二区蜜桃 | 亚洲第一av免费看| 国产男人的电影天堂91| 成人综合一区亚洲| 精品久久国产蜜桃| 一区二区日韩欧美中文字幕 | 日日啪夜夜爽| 精品一区二区免费观看| 久久久国产欧美日韩av| 亚洲成国产人片在线观看| 99久久人妻综合| 免费在线观看黄色视频的| 美女视频免费永久观看网站| 欧美日韩精品成人综合77777| 欧美成人午夜精品| 一区二区三区精品91| 男人爽女人下面视频在线观看| 免费观看av网站的网址| 欧美 亚洲 国产 日韩一| 99视频精品全部免费 在线| 国产免费视频播放在线视频| 九草在线视频观看| 国产1区2区3区精品| 捣出白浆h1v1| 免费观看无遮挡的男女| videosex国产| 午夜91福利影院| 欧美日韩综合久久久久久| 丰满少妇做爰视频| 亚洲第一区二区三区不卡| 2021少妇久久久久久久久久久| 久久人妻熟女aⅴ| av不卡在线播放| 免费av不卡在线播放| 亚洲高清免费不卡视频| 中文字幕人妻丝袜制服| 99久国产av精品国产电影| 日本-黄色视频高清免费观看| 一本大道久久a久久精品| 丝袜在线中文字幕| 婷婷色麻豆天堂久久| 人妻系列 视频| 国产精品国产三级专区第一集| 99热6这里只有精品| 欧美日韩成人在线一区二区| 999精品在线视频| 欧美日韩成人在线一区二区| 欧美 日韩 精品 国产| 如日韩欧美国产精品一区二区三区| 亚洲av综合色区一区| 如日韩欧美国产精品一区二区三区| 婷婷成人精品国产| 亚洲美女视频黄频| 又黄又爽又刺激的免费视频.| xxx大片免费视频| 男女无遮挡免费网站观看| 国产精品免费大片| 久久亚洲国产成人精品v| 丰满迷人的少妇在线观看| 人人妻人人爽人人添夜夜欢视频| 久久 成人 亚洲| av有码第一页| 亚洲欧美清纯卡通| 久久久久久久久久久久大奶| 亚洲精品第二区| 91国产中文字幕| 如何舔出高潮| 亚洲国产av新网站| 夜夜爽夜夜爽视频| 在线观看国产h片| 人人妻人人爽人人添夜夜欢视频| 黄色毛片三级朝国网站| 国产精品久久久久久久电影| a 毛片基地| 最新的欧美精品一区二区| 久久久精品区二区三区| 成人国产麻豆网| 在线观看人妻少妇| 国产精品欧美亚洲77777| 宅男免费午夜| 国产免费视频播放在线视频| 狠狠婷婷综合久久久久久88av| 精品一区二区三区视频在线| 欧美日韩成人在线一区二区| 午夜精品国产一区二区电影| 亚洲精品乱久久久久久| 考比视频在线观看| 天美传媒精品一区二区| 国产精品一区二区在线观看99| 夜夜骑夜夜射夜夜干| 日本-黄色视频高清免费观看| 9热在线视频观看99| 久久99精品国语久久久| 久久人人97超碰香蕉20202| 我的女老师完整版在线观看| 精品人妻在线不人妻| 成人影院久久| av有码第一页| 亚洲四区av| 91午夜精品亚洲一区二区三区| 国产精品国产三级国产专区5o| 免费观看a级毛片全部| 乱人伦中国视频| 亚洲av福利一区| 国产又色又爽无遮挡免| 亚洲,一卡二卡三卡| 99国产综合亚洲精品| 国内精品宾馆在线| 国产精品国产三级专区第一集| 飞空精品影院首页| 一本色道久久久久久精品综合| 久久精品熟女亚洲av麻豆精品| 欧美另类一区| 99re6热这里在线精品视频| 大片免费播放器 马上看| 亚洲欧美一区二区三区黑人 | 国产伦理片在线播放av一区| 精品少妇内射三级| 午夜福利在线观看免费完整高清在| 亚洲欧美中文字幕日韩二区| 国产 一区精品| 日韩制服丝袜自拍偷拍| 一边亲一边摸免费视频| 色5月婷婷丁香| 毛片一级片免费看久久久久| 成人毛片60女人毛片免费| 建设人人有责人人尽责人人享有的| 欧美人与性动交α欧美精品济南到 | 欧美成人午夜精品| 日韩中字成人| 巨乳人妻的诱惑在线观看| 性高湖久久久久久久久免费观看| 看免费av毛片| 免费观看av网站的网址| 一级黄片播放器| 欧美日韩精品成人综合77777| 国产亚洲av片在线观看秒播厂| 多毛熟女@视频| 日韩在线高清观看一区二区三区| www.熟女人妻精品国产 | 爱豆传媒免费全集在线观看| 国产欧美日韩一区二区三区在线| 伊人亚洲综合成人网| 中文字幕精品免费在线观看视频 | 午夜福利影视在线免费观看| 大香蕉久久网| 青春草亚洲视频在线观看| 亚洲,欧美精品.| 在线观看国产h片| 熟女av电影| 国产成人精品婷婷| 成年人午夜在线观看视频| 伦精品一区二区三区| 男女边吃奶边做爰视频| 九九爱精品视频在线观看| 观看av在线不卡| 亚洲成人一二三区av| 最新的欧美精品一区二区| 黄色视频在线播放观看不卡| 不卡视频在线观看欧美| 九草在线视频观看| 午夜久久久在线观看| 国产成人欧美| 午夜av观看不卡| 欧美成人精品欧美一级黄| 精品国产一区二区久久| 热re99久久精品国产66热6| 国产精品无大码| 久久精品久久久久久久性| av福利片在线| 国产日韩欧美视频二区| 国产精品 国内视频| 欧美日韩av久久| 久久人妻熟女aⅴ| 夜夜骑夜夜射夜夜干| 香蕉丝袜av| 春色校园在线视频观看| 一级毛片黄色毛片免费观看视频| 丝袜美足系列| 久久久久久久久久久免费av| 久久精品久久精品一区二区三区| 国产午夜精品一二区理论片| 七月丁香在线播放| 日韩制服丝袜自拍偷拍| 草草在线视频免费看| 久久久久久人人人人人| 69精品国产乱码久久久| 黄色 视频免费看| 日本色播在线视频| 国产成人欧美| 男人操女人黄网站| 亚洲国产欧美日韩在线播放| 这个男人来自地球电影免费观看 | 满18在线观看网站| 日韩免费高清中文字幕av| 亚洲欧美一区二区三区黑人 | 国产69精品久久久久777片| 波多野结衣一区麻豆| 欧美亚洲日本最大视频资源| 国产极品粉嫩免费观看在线| 黑人巨大精品欧美一区二区蜜桃 | 国产精品三级大全| 午夜久久久在线观看| 高清视频免费观看一区二区| 午夜福利视频精品| 亚洲欧美精品自产自拍| 99久国产av精品国产电影| 香蕉国产在线看| 亚洲精品日本国产第一区| 亚洲av中文av极速乱| 久热久热在线精品观看| 亚洲美女黄色视频免费看| 18禁观看日本| 久久这里只有精品19| 美女视频免费永久观看网站| 亚洲国产日韩一区二区| 国产免费一区二区三区四区乱码| 一级毛片黄色毛片免费观看视频| 曰老女人黄片| 欧美精品人与动牲交sv欧美| 亚洲欧美中文字幕日韩二区| 国产精品国产三级专区第一集| 最近2019中文字幕mv第一页| 男女下面插进去视频免费观看 | 国产欧美亚洲国产| 高清不卡的av网站| 国产精品.久久久| 我要看黄色一级片免费的| 国产福利在线免费观看视频| 综合色丁香网| 91精品三级在线观看| av国产久精品久网站免费入址| 一级片免费观看大全| 亚洲av综合色区一区| 日韩中字成人| 国产在线一区二区三区精| 18禁裸乳无遮挡动漫免费视频| 午夜视频国产福利| 欧美精品av麻豆av| 亚洲久久久国产精品| 人人澡人人妻人| 在线天堂最新版资源| 亚洲五月色婷婷综合| 搡女人真爽免费视频火全软件| 美女中出高潮动态图| 久久久精品94久久精品| 人人妻人人爽人人添夜夜欢视频| 美女内射精品一级片tv| 中文字幕另类日韩欧美亚洲嫩草| 国产精品无大码| 亚洲少妇的诱惑av| 亚洲精品国产色婷婷电影| 少妇高潮的动态图| 午夜老司机福利剧场| 男女下面插进去视频免费观看 | 久久99精品国语久久久| 丝袜在线中文字幕| 亚洲美女黄色视频免费看| 乱码一卡2卡4卡精品| 久久久久人妻精品一区果冻| 色网站视频免费| 中文乱码字字幕精品一区二区三区| 国产精品 国内视频| 丝袜在线中文字幕| 精品久久蜜臀av无| 一区二区日韩欧美中文字幕 | 亚洲精品成人av观看孕妇| 亚洲美女视频黄频| 人妻 亚洲 视频| 一区二区三区四区激情视频| 国产av精品麻豆| 欧美xxxx性猛交bbbb|