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

    突出重點目標跟蹤的分布式MIMO雷達陣元選取

    2021-11-11 06:15:22蔣春啟鄭娜娥
    系統(tǒng)工程與電子技術 2021年10期
    關鍵詞:復雜度雷達數(shù)量

    蔣春啟, 鄭娜娥, 左 宗, 王 盛, 陳 翔

    (1. 戰(zhàn)略支援部隊信息工程大學數(shù)據(jù)與目標工程學院, 河南 鄭州 450001;2. 電子信息系統(tǒng)復雜電磁環(huán)境效應國家重點實驗室, 河南 洛陽 471003)

    0 引 言

    多輸入多輸出(multiple input multiple output,MIMO)雷達的概念最早提出于2003年[1],主要分為分布式MIMO雷達[2]和集中式MIMO雷達[1]兩種。其所具有的分集增益,能夠有效抑制干擾、提高分辨率、克服目標的雷達散射截面(radar cross section, RCS)閃爍問題,因此得到了諸多學者和科研機構的廣泛關注,研究方向主要集中在波形設計[3-4],檢測性能分析及檢測器研究[5-8],參數(shù)估計[9-11],陣元、目標檢測或跟蹤[12]、功率和帶寬等系統(tǒng)資源的分配[13-21]等。

    MIMO雷達系統(tǒng)性能優(yōu)勢離不開陣元、功率、波形、頻率等系統(tǒng)資源的合理配置。系統(tǒng)資源的合理優(yōu)化分配可以有效提高MIMO雷達的系統(tǒng)性能。本文針對系統(tǒng)資源更為豐富的分布式MIMO雷達系統(tǒng)的陣元選取問題展開研究。

    通過優(yōu)化配置,可以用更少的系統(tǒng)資源達到同樣的系統(tǒng)性能,降低系統(tǒng)計算復雜度,提升系統(tǒng)綜合效能,因此具有重要的研究價值。在這方面已有諸多的研究成果,文獻[15]采用啟發(fā)式算法將陣元選取建模為背包問題(knapsack problem, KP),降低了計算復雜度。文獻[16]采用了陣元聚類分組的方法,每個目標只由對應的陣元子集跟蹤,分別提出了在復雜度上更有優(yōu)勢的貪婪多起點搜索(greedy multi-start local search, GMLS)算法和在跟蹤精度上更有優(yōu)勢的公平多起點搜索(fair multi-start local search, FMLS)算法。文獻[17]選取最少陣元的同時盡量降低系統(tǒng)的代價損耗。以GMLS算法為基礎,文獻[18-19]以GMLS算法和FMLS算法為參考并進行改進,提出了同時兼顧系統(tǒng)性能和計算復雜度的改進GMLS(modified GMLS, MGMLS)算法。文獻[19-20]在研究多目標跟蹤的陣元選取算法時將目標分為普通目標、可疑目標和危險目標,分別就低跟蹤精度要求場景和高跟蹤精度要求場景進行討論,并分別提出了貪婪與公平多起點搜索(greedy and fair multi-start local search, GFMLS)算法和單天線對多目標(one antenna to all targets, OAAT)算法。GFMLS算法數(shù)據(jù)處理簡單,適合處理低跟蹤精度要求下不同優(yōu)先級目標的陣元選取問題,OAAT算法處理數(shù)據(jù)量大,但可以實現(xiàn)更高要求的目標跟蹤性能。

    但是以上研究都是針對單任務,實際情況中雷達系統(tǒng)常需要面臨多任務的情況,最常見的任務為目標檢測和目標跟蹤。雷達在執(zhí)行多任務時,會帶來新的資源分配問題,需要考慮任務間的資源分配,當資源不充分時,如何優(yōu)先保障優(yōu)先級高的任務。因此,需要對分布式MIMO雷達多任務資源分配技術展開研究。Zhang等在文獻[21]中研究了多任務下分布式MIMO雷達的資源分配問題,從提高系統(tǒng)性能和降低復雜度兩個方面著手,綜合考慮目標跟蹤和目標檢測,在系統(tǒng)資源上同時考慮陣元、功率和帶寬的分配,在陣元選取上采用了改進的粒子群算法,將MIMO雷達的資源分配問題從單任務推進到多任務。但文獻[21]在綜合考慮不同任務的性能指標時,沒有將具有不同量綱的檢測概率和定位誤差這兩個性能指標進行去量綱化處理,將不同量綱的參數(shù)強制相加,得出的目標函數(shù)缺乏實際意義;同時在進行檢測概率的計算時,未能覆蓋整個待檢測區(qū)域,僅考慮了其邊界線,實際應用范圍受限;且該文獻平等對待需跟蹤的每個目標,不能適應需要重點跟蹤特定目標的情況。當在監(jiān)視區(qū)域內(nèi)的多個目標中需要重點關注特定目標且系統(tǒng)資源有限時,需要合理進行資源調(diào)度以保證重點目標的跟蹤性能,確保資源的高效利用。本文針對上述問題,在文獻[17, 19-21]的基礎上提出了基于改進FMLS(modified FMLS, MFMLS)的陣元選取模型及求解算法。在計算檢測概率時,全面考慮整個待檢測區(qū)域,同時針對多任務的情景,對于不同任務不同性能指標的不同量綱的問題,進行去量綱化和歸一化處理,使其具有可比性,在多目標跟蹤中突出重點目標的性能保障,使選用的陣元數(shù)量盡量少,以此建立優(yōu)化模型并求解。

    1 系統(tǒng)模型

    假設目標勻速運動,目標的運動方程為

    (1)

    (2)

    (3)

    式中:κ0為噪聲的強度;I2為2×2的單位陣;?為Kronecker積。

    可得目標的運動方程為

    (4)

    s(t)=[s1(t),s2(t),…,sM(t)]T

    (5)

    (6)

    假設系統(tǒng)能夠?qū)崿F(xiàn)同步,第n個接收陣元接收對第p個目標的低通等效信號[22]為

    (7)

    記為

    (8)

    (9)

    2 重要參數(shù)

    2.1 目標定位誤差的貝葉斯克拉美羅界

    (10)

    (11)

    (12)

    (13)

    (14)

    (15)

    (16)

    2.2 目標的檢測概率

    在考慮分布式MIMO雷達對監(jiān)視區(qū)內(nèi)特定點的檢測時,可以不考慮多普勒效應。以文獻[21,25]為參考,在推導檢測概率時可采用Neyman-Pearson準則。為便于研究,將噪聲考慮為高斯白噪聲,雷達發(fā)射單脈沖探測信號,利用最大似然比檢測器,對于MIMO雷達不同路徑,首先進行匹配濾波處理,其次作平方律檢波,最后相加[25]。由第m部發(fā)射陣元發(fā)射由第n部接收陣元接收,故系統(tǒng)應同時有M×N個輸出。

    將匹配濾波器輸出記為Xr,其為一個MN×1向量,當目標不存在時,用g(r(t)|H0)表示信號的概率密度函數(shù),以g(r(t)|H1) 表示目標存在時信號的概率密度函數(shù)。則似然比檢測為

    (17)

    若檢測門限為δthreshold,當‖Xr‖2≥δthreshold時,表示檢測到目標,當‖Xr‖2<δthreshold時表示沒有檢測到目標。

    (18)

    (19)

    似然比檢測門限δthreshold由虛警概率確定。則相應的虛警概率為

    (20)

    可以得到檢測門限δthreshold為

    (21)

    由此可以推導出目標檢測概率為

    (22)

    3 陣元選取模型的建立與求解

    文獻[17]提出的FMLS算法通過遍歷所有的初始陣元對,進行多輪啟發(fā)式選取,能夠獲得較高的跟蹤精度,但沒有考慮不同目標不同跟蹤精度要求的問題。文獻[19-20]提出的OAAT算法將目標分為普通目標、可疑目標和危險目標,就高跟蹤精度要求場景進行了討論,但仍局限于單任務的情形。本文借鑒FMLS算法和OAAT算法,提出了基于MFMLS的陣元選取模型及求解算法。

    3.1 模型的建立

    在資源有限的情況下進行優(yōu)化模型的建立。本文為最大限度地節(jié)約系統(tǒng)資源,減少陣元使用數(shù)量,以發(fā)射和接收陣元數(shù)量之和為目標函數(shù),并使其最小化。在系統(tǒng)的監(jiān)視區(qū)域內(nèi),區(qū)分對待不同目標,限定發(fā)射陣元、接收陣元的數(shù)量總量和最小檢測概率,同時限制每次可以使用的陣元使用率,以上述條件為約束條件,建立優(yōu)化模型如下:

    (23)

    式中:以最小的陣元集為目標函數(shù),ηp是目標p允許的最大位置估計誤差,Pthreshold為檢測概率最小閾值,發(fā)射陣元的數(shù)量不大于M,接收陣元的數(shù)量不大于N,系統(tǒng)總的陣元使用率不大于系統(tǒng)容許的最大陣元使用率σ。

    3.2 參數(shù)預處理

    3.2.1 目標位置估計精度誤差歸一化

    由于系統(tǒng)同時承擔目標跟蹤和目標檢測雙任務,需要同時考慮檢測概率Pd和跟蹤目標的位置估計精度誤差Emse兩項指標。從理論上講,Pd和Emse的取值范圍分別為[0,1] 和[0,+∞),兩者具有不同的量綱,難以將系統(tǒng)的跟蹤性能和檢測性能在同一框架下進行比較和討論,這也對多任務下的資源分配提出了挑戰(zhàn)。一個思路是對不同任務的性能指標參數(shù)進行去量綱化和歸一化處理,以便使Emse和Pd具有可比性。為此需要將Emse進行歸一化處理,使其取值范圍由[0,+∞)轉(zhuǎn)為[0,1]。為了實現(xiàn)Emse取值范圍的轉(zhuǎn)換,引入?yún)?shù)Emax。Emax為系統(tǒng)進行目標跟蹤時可能達到的最大位置估計精度誤差,或者說當系統(tǒng)對某一目標的位置估計誤差超過Emax時將認定為無效跟蹤,失去對目標進行跟蹤的意義。Emax的值可通過以下方法進行確定,一是進行仿真分析,以最少的資源進行目標的跟蹤,得出系統(tǒng)可能達到的目標定位誤差的最大值,以此作為參考,使Emax不小于該值;二是依據(jù)實際需求而定,當目標定位誤差超過一定值時,跟蹤將失去意義,以此作為Emax的值;三是結(jié)合雷達系統(tǒng)對目標跟蹤的經(jīng)驗積累和仿真分析,人為給出Emax的值。下面給出Emse歸一化的具體步驟。

    步驟 1給定Emax,使得Emse∈[0,Emax];

    步驟 2去量綱處理,令E1=Emse/Emax,使得E1∈[0,1];

    經(jīng)過步驟1和步驟2,雖然已經(jīng)實現(xiàn)了Emse的歸一化處理,使其取值范圍與Pd一致。但兩者仍然存在區(qū)別,Pd的值越大代表檢測性能越好,而Emse的值越小表示跟蹤性能越好。為了進一步消除不一致,在前兩個步驟的基礎上繼續(xù)采取如下步驟。

    步驟 3令E2=E1·(-1)+1,使得E2∈[0,1]。

    至此完成了Emse的歸一化處理。

    3.2.2 不同目標的權重處理

    對于系統(tǒng)需要跟蹤的多個目標,重點目標和普通目標的跟蹤精度要求不同,因此在考慮綜合跟蹤性能時需要對不同的目標區(qū)別對待,為此需要引入權重系數(shù)。

    目標的重要程度與其跟蹤精度要求有關,跟蹤精度要求越高,目標的重要程度越高,而跟蹤精度的數(shù)值卻越小,因此表示為目標p重要程度的權重系數(shù)ξp與跟蹤精度數(shù)值ηp成反比,即

    (24)

    式中:ξ為比例系數(shù),在實際應用中ηp為已知條件。假定跟蹤目標的數(shù)量為Np,所有目標的權重系數(shù)之和為1,則有

    (25)

    通過式(24)和式(25)可求得ξ和ξp。若通過上式確定的ξp還不足以突出重點目標,可以引入人為干擾因子進行調(diào)控。例如,當系統(tǒng)需要跟蹤的目標為2個,一個為重點目標,一個為普通目標,通過式(24)和式(25)最終求得重點目標的權重系數(shù)為ξim,普通目標的權重系數(shù)為ξnor,引入人為干擾因子Vξ以進一步加強對重點目標的跟蹤,得到最終點目標的權重系數(shù)為ξim+Vξ,最終普通目標的權重系數(shù)為ξnor-Vξ。

    3.2.3 系統(tǒng)綜合性能評價

    Perk=δ1(Pd-Pthreshold)+

    (26)

    式中:δ1和δ2分別為系統(tǒng)檢測性能和跟蹤性能的權重系數(shù)。當δ1>δ2時,系統(tǒng)更加側(cè)重于檢測任務;當δ1<δ2時,則更加側(cè)重于跟蹤任務;當δ1=δ2時,同樣對待跟蹤任務和檢測任務。在實際應用中δ1和δ2的值由操作人員根據(jù)實際任務需求人為確定。

    當Perk=0時表示系統(tǒng)性能正好達到性能要求;當Perk>0時表示系統(tǒng)性能優(yōu)于性能要求,且Perk越大性能越好;當Perk<0時表示系統(tǒng)性能劣于性能要求,且Perk越小性能越差。

    通過參數(shù)預處理,優(yōu)化模型可重新表示為

    (27)

    3.3 模型求解

    使用窮舉法進行陣元選取可以得到系統(tǒng)的最優(yōu)解,但其計算量達到2(M+N)σ,當陣元數(shù)量較多時,計算量非常大。此處假定系統(tǒng)跟蹤的目標分為重點目標和普通目標兩種,利用MFMLS算法進行陣元選取,具體步驟如下。

    步驟 2依據(jù)陣元集Amin計算各目標的跟蹤性能和檢測性能,得到每個目標的跟蹤精度誤差和檢測概率。若系統(tǒng)性能達到要求,進入步驟13,若系統(tǒng)性能不能達到要求,轉(zhuǎn)步驟3。

    步驟 3若重點目標的跟蹤精度得到滿足,轉(zhuǎn)步驟8,否則轉(zhuǎn)步驟4。

    步驟 4依據(jù)陣元集Amin,并由式計算系統(tǒng)綜合性能,記為Per。

    步驟 5在陣元集B中選取1個陣元X加入到陣元集Amin中,并更新Amin和B,即Amin=Amin∪{X},B=B{X}。根據(jù)更新后的Amin計算系統(tǒng)的綜合性能,記為Per′,并計算Δper=Per′-Per。隨后在Amin去掉陣元X,再次更新Amin和B,Amin=Amin{X},B=B∪{X}。

    步驟 6重復步驟5,遍歷集合B中的每一個陣元,得到使Δper取得最大值的陣元X′,將陣元X′加入到集Amin中,并更新Amin和B,即Amin=Amin∪{X′},B=B{X′}。

    步驟 7重復步驟4~步驟6,直到重點目標的跟蹤性能得到滿足或者Amin中陣元是數(shù)量達到系統(tǒng)容許的最大陣元使用率。若Amin中陣元數(shù)量達到系統(tǒng)容許的最大陣元使用率,轉(zhuǎn)步驟13,否則轉(zhuǎn)步驟8。

    步驟 8依據(jù)陣元集Amin計算對不同特定目標的跟蹤性能和檢測性能,得到每個目標的跟蹤精度誤差和檢測概率。若系統(tǒng)性能達到要求,轉(zhuǎn)步驟13,若系統(tǒng)性能不能達到要求,轉(zhuǎn)步驟9。

    步驟 9依據(jù)陣元集Amin,忽略掉重點目標并由式(26)計算系統(tǒng)綜合性能,記為Per。

    步驟 10在陣元集B中選取1個陣元X加入到陣元集Amin中,并更新Amin和B,即Amin=Amin∪{X},B=B{X}。根據(jù)更新后的Amin忽略掉重點目標并計算系統(tǒng)的綜合性能,記為Per′,并計算Δper=Per′-Per。隨后在Amin去掉陣元X,再次更新Amin和B,Amin=Amin{X},B=B∪{X}。

    步驟 11重復步驟9,遍歷集合B中的每一個陣元,得到使Δper取得最大值的陣元X′,將陣元X′加入到集Amin中,并更新Amin和B,即Amin=Amin∪{X′},B=B{X′}。

    步驟 12重復步驟8至步驟10,直到系統(tǒng)性能得到滿足或者Amin中陣元是數(shù)量達到系統(tǒng)容許的最大陣元使用率。

    步驟 13重復步驟1~步驟12,遍歷所有的初始陣元對,得到諸多陣元集Amin。在諸多的陣元集Amin中選擇在滿足系統(tǒng)性能且陣元數(shù)量最少的陣元集合,若有多個陣元數(shù)量同樣少的陣元集合,則選擇其中綜合性能最好的陣元集合;若系統(tǒng)性能不能得到滿足,則在諸多的陣元集Amin中選擇能夠滿足重點目標跟蹤性能并盡量提高其余性能的陣元集;若重點目標的跟蹤性能不能得到滿足,則選擇系統(tǒng)綜合性能最好的陣元集。

    3.4 算法復雜度分析

    進行陣元選取時,不同算法的復雜度不同,窮舉法需要進行2(M+N)σ次陣元選取。MFMLS算法最多需要進行((M+N)σ-2)((M+N)σ-1)MN/2次陣元選取,通過對比可以看出,當系統(tǒng)的陣元數(shù)量較多時,本文算法復雜度降低明顯。

    4 仿真實驗與分析

    假設分布式MIMO雷達系統(tǒng)分布在6 km×6 km的二維笛卡爾坐標系內(nèi),發(fā)射陣元數(shù)M=8,接收陣元數(shù)N=8,形成圓形雷達陣,雷達使用使用率σ=0.75。假定待跟蹤的目標數(shù)量為NQ=3,目標初始位置為(-2 800,1 000) m、(-150,800) m和(-150,600) m,以60 m/s的速度沿不同方向運動。第一個目標為重點目標,位置估計精度要求為η=9 m,加入人為干擾因子后其權重系數(shù)ξ1=0.9;另兩個目標為普通目標,位置估計精度要求為η=18 m,加入人為干擾因子后其權重系數(shù)ξ2=ξ3=0.05。檢測概率最小閾值要求為Pd≥0.9,計算檢測概率時,平等對待監(jiān)視區(qū)域內(nèi)的每一個點。以Δt=2 s的時間間隔進行取樣觀測,從觀測數(shù)據(jù)中取12幀數(shù)據(jù)進行對比分析。每個發(fā)射陣元用于不同任務的發(fā)射功率均為pmax=1 500 W,信號有效帶寬、有效時寬、波長分別為1.5 MHz、10 μs、0.3 m,脈沖重復頻率為fr=53 kHz。假設系統(tǒng)中目標對雷達各方向散射截面積均相同,取散射系數(shù)均為0.9。此處平等對待不同任務,故取δ1=δ2=1,同時取Emax=50 m。仿真采用Num=500次蒙特卡羅進行統(tǒng)計平均。

    首先給出陣元在系統(tǒng)中的空間位置及目標的運動軌跡,如圖1所示。

    圖1 雷達與目標的位置及目標的運動軌跡Fig.1 Radar and target position and target’s motion trajectory

    為了驗證本文算法的有效性,此處將MFMLS算法與其他算法進行比較。一是最少陣元數(shù)量窮舉法,用窮舉法尋找滿足系統(tǒng)性能且陣元數(shù)量最少的陣元集,將此處的窮舉法稱為最少陣元數(shù)量窮舉法。二是首先利用MFMLS算法得到陣元數(shù)量,其次用該陣元數(shù)量為參數(shù),利用窮舉法進行陣元選取,此時,將其稱為固定陣元數(shù)量窮舉法;三是首先利用MFMLS算法得到陣元數(shù)量,其次用該陣元數(shù)量為參數(shù),隨機選取發(fā)射和接收陣元,將這種方法稱為隨機算法。圖2~圖5分別給出了4種算法在檢測概率和3個跟蹤目標定位誤差方面的對比。可以看到,隨機算法無法滿足系統(tǒng)的檢測性能要求,也不能滿足重點目標的跟蹤精度要求,不區(qū)分地對待3個跟蹤目標;而MFMLS算法在各方面性能上均和固定陣元數(shù)量窮舉法接近,且部分優(yōu)于最少陣元數(shù)量窮舉法,三者均能滿足系統(tǒng)性能要求,突出對重點目標的跟蹤。

    圖2 檢測概率對比Fig.2 Comparison of detection probability

    圖3 目標1位置估計誤差對比Fig.3 Comparison of position estimation error of target 1

    圖4 目標2位置估計誤差對比Fig.4 Comparison of position estimation error of target 2

    圖5 目標3位置估計誤差對比Fig.5 Comparison of position estimation error of target 3

    從以上分析上可知,MFMLS算法能夠滿足系統(tǒng)各個單方面的性能要求,下面進行系統(tǒng)綜合性能分析。圖6給出了4種算法的綜合性能對比圖,可以看到MFMLS算法的綜合性能優(yōu)于隨機法,與固定陣元數(shù)量窮舉法接近,且由于在陣元數(shù)量上不能達到最少,故在性能上相對于最少數(shù)量窮舉法也表現(xiàn)出了一定的優(yōu)越性。除了隨機法沒能達到系統(tǒng)的性能要求外,其余方法均達到了要求,這時MFMLS算法便在計算復雜度上便體現(xiàn)了優(yōu)勢。圖7是計算復雜度對比,可以看出,相對于最少陣元數(shù)量窮舉法和固定陣元數(shù)量窮舉法,MFMLS算法在系統(tǒng)計算量上均能得到大幅降低。進一步由圖7中數(shù)據(jù)分析可知,MFMLS算法在不同幀相對于最少數(shù)量窮舉法的計算復雜度的降幅均在93%以上,相對于固定陣元數(shù)量窮舉法降幅在35%~65%。在陣元數(shù)量上,最少陣元數(shù)量窮舉法能夠在滿足性能的前提下達到最少陣元,由于涉及到多任務而且不同目標的跟蹤精度要求也不同,難以完全做到陣元數(shù)量最少,但500次蒙特卡羅平均值后,MFMLS算法在陣元數(shù)量上接近最少陣元數(shù)量窮舉法。圖8給出了陣元數(shù)量對比圖。

    圖6 綜合性能對比圖Fig.6 Comparison of comprehensive performance

    圖7 計算復雜度對比Fig.7 Comparison of computational complexity

    圖8 陣元數(shù)量對比Fig.8 Comparison of the number of antennas

    5 結(jié) 論

    本文研究了多目標跟蹤和目標檢測多任務下的分布式MIMO雷達的陣元選取問題,且不同目標具有不同的跟蹤精度誤差。為了突出對重點目標的跟蹤并兼顧系統(tǒng)的整體性能,本文提出了基于MFMLS的陣元選取模型及求解算法。通過仿真實驗,利用MFMLS算法進行陣元選取能夠滿足系統(tǒng)性能要求,在系統(tǒng)性能上本文算法優(yōu)于隨機算法,在計算復雜度上,本文算法相對于最少陣元數(shù)量窮舉算法和固定陣元數(shù)量窮舉算法均有大幅下降,由此驗證了本文算法的有效性。但本文的研究同樣存在局限性,需要展開進一步的研究,如在目標運動規(guī)律上,本文的跟蹤目標為勻速運動,多任務下機動目標的跟蹤仍需研究;在目標數(shù)量上,本文考慮的跟蹤目標數(shù)量較少,大規(guī)模集群目標的跟蹤仍需研究;在系統(tǒng)資源上,本文僅考慮了陣元因素,多任務多資源因素下的資源分配仍需研究。

    猜你喜歡
    復雜度雷達數(shù)量
    有雷達
    大自然探索(2023年7期)2023-08-15 00:48:21
    統(tǒng)一數(shù)量再比較
    一種低復雜度的慣性/GNSS矢量深組合方法
    雷達
    求圖上廣探樹的時間復雜度
    頭發(fā)的數(shù)量
    某雷達導51 頭中心控制軟件圈復雜度分析與改進
    基于空時二維隨機輻射場的彈載雷達前視成像
    現(xiàn)代“千里眼”——雷達
    出口技術復雜度研究回顧與評述
    好男人电影高清在线观看| 黄色视频,在线免费观看| 一区二区三区精品91| 中文字幕色久视频| 黑人操中国人逼视频| 久久久精品区二区三区| 高清视频免费观看一区二区| 国产色视频综合| 欧美另类一区| 欧美日韩精品网址| 久久亚洲精品不卡| 国产福利在线免费观看视频| 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| 热99久久久久精品小说推荐| 亚洲精品日韩在线中文字幕| 国产精品一区二区免费欧美 | 午夜福利,免费看| 国产精品 国内视频| 国产麻豆69| 亚洲 欧美一区二区三区| 又紧又爽又黄一区二区| 亚洲精品一二三| 最近最新免费中文字幕在线| 国产人伦9x9x在线观看| 亚洲专区国产一区二区| 国产伦理片在线播放av一区| 免费观看a级毛片全部| 亚洲国产精品一区三区| 欧美日韩黄片免| 国产伦理片在线播放av一区| 精品熟女少妇八av免费久了| 真人做人爱边吃奶动态| 午夜福利视频在线观看免费| 男女下面插进去视频免费观看| 亚洲精品乱久久久久久| 下体分泌物呈黄色| 18禁国产床啪视频网站| 这个男人来自地球电影免费观看| 欧美少妇被猛烈插入视频| 男人爽女人下面视频在线观看| 免费黄频网站在线观看国产| 精品免费久久久久久久清纯 | 亚洲人成电影免费在线| 亚洲欧美一区二区三区久久| 在线观看一区二区三区激情| 男女下面插进去视频免费观看| 婷婷丁香在线五月| 亚洲精品一卡2卡三卡4卡5卡 | 性少妇av在线| 99国产精品一区二区三区| av在线老鸭窝| 狠狠狠狠99中文字幕| 久久久久视频综合| 日韩 亚洲 欧美在线| 欧美黑人精品巨大| 久久精品亚洲av国产电影网| 十八禁网站网址无遮挡| 超色免费av| 国产日韩一区二区三区精品不卡| 激情视频va一区二区三区| 91麻豆精品激情在线观看国产 | 亚洲全国av大片| 大香蕉久久网| 操美女的视频在线观看| 曰老女人黄片| 美女高潮喷水抽搐中文字幕| 亚洲五月婷婷丁香| 久久国产亚洲av麻豆专区| 纯流量卡能插随身wifi吗| 大香蕉久久网| 91麻豆av在线| 亚洲成人免费电影在线观看| www.999成人在线观看| 久久精品亚洲av国产电影网| 丝袜美足系列| 美女中出高潮动态图| 人成视频在线观看免费观看| 搡老乐熟女国产| 欧美成人午夜精品| 欧美激情极品国产一区二区三区| 99热网站在线观看| a级毛片黄视频| 考比视频在线观看| 国产一区二区三区av在线| 午夜免费鲁丝| 97在线人人人人妻| 久久免费观看电影| 日本av手机在线免费观看| 少妇粗大呻吟视频| www.av在线官网国产| 国产欧美日韩综合在线一区二区| 美国免费a级毛片| 亚洲国产毛片av蜜桃av| 午夜福利在线观看吧| av天堂久久9| 国产亚洲一区二区精品| av有码第一页| 国产男人的电影天堂91| 在线观看www视频免费| 男女床上黄色一级片免费看| 久久久精品国产亚洲av高清涩受| 国产亚洲精品一区二区www | 久久国产精品男人的天堂亚洲| 在线看a的网站| 夫妻午夜视频| 午夜视频精品福利| 中文字幕另类日韩欧美亚洲嫩草| 免费av中文字幕在线| 秋霞在线观看毛片| 亚洲精品美女久久av网站| 亚洲人成电影免费在线| 精品久久久精品久久久| 考比视频在线观看| 视频区图区小说| 久久精品国产综合久久久| tube8黄色片| 男女下面插进去视频免费观看| 丰满人妻熟妇乱又伦精品不卡| 欧美激情高清一区二区三区| 亚洲avbb在线观看| 黄网站色视频无遮挡免费观看| 国产一区二区三区综合在线观看| 亚洲伊人久久精品综合| 手机成人av网站| 亚洲精品一二三| 麻豆国产av国片精品| 50天的宝宝边吃奶边哭怎么回事| 国产日韩欧美亚洲二区| 中国美女看黄片| 50天的宝宝边吃奶边哭怎么回事| 少妇精品久久久久久久| 免费高清在线观看视频在线观看| 国产成人精品久久二区二区免费| 一区二区日韩欧美中文字幕| 国产精品偷伦视频观看了| 50天的宝宝边吃奶边哭怎么回事| 在线观看人妻少妇| 久久精品亚洲av国产电影网| 国产亚洲精品一区二区www | 一区福利在线观看| 精品熟女少妇八av免费久了| 亚洲精品国产av成人精品| 午夜免费观看性视频| 桃红色精品国产亚洲av| 人妻人人澡人人爽人人| kizo精华| 桃花免费在线播放| 国产老妇伦熟女老妇高清| 久久精品aⅴ一区二区三区四区| 十八禁人妻一区二区| 亚洲色图综合在线观看| 日韩视频在线欧美| 99香蕉大伊视频| 国产成+人综合+亚洲专区| 日韩 欧美 亚洲 中文字幕| 黄片大片在线免费观看| 亚洲精品国产av成人精品| 国产精品香港三级国产av潘金莲| 波多野结衣一区麻豆| 免费一级毛片在线播放高清视频 | 国产在线一区二区三区精| 一个人免费在线观看的高清视频 | 曰老女人黄片| 国产又爽黄色视频| 国产一区二区三区av在线| 热re99久久国产66热| 精品亚洲成国产av| 这个男人来自地球电影免费观看| 王馨瑶露胸无遮挡在线观看| 日韩电影二区| 汤姆久久久久久久影院中文字幕| xxxhd国产人妻xxx| 午夜影院在线不卡| 免费观看a级毛片全部| 精品卡一卡二卡四卡免费| 涩涩av久久男人的天堂| 嫁个100分男人电影在线观看| 国产99久久九九免费精品| 亚洲精品一区蜜桃| cao死你这个sao货| 亚洲 欧美一区二区三区| 国产日韩欧美亚洲二区| 中国美女看黄片| 久久青草综合色| 国产日韩欧美亚洲二区| 久久久久国内视频| 一本综合久久免费| 精品少妇黑人巨大在线播放| 亚洲成国产人片在线观看| 亚洲精品第二区| 深夜精品福利| 少妇人妻久久综合中文| 在线十欧美十亚洲十日本专区| 欧美变态另类bdsm刘玥| av福利片在线| 国产精品久久久久久精品电影小说| 日韩精品免费视频一区二区三区| 欧美精品人与动牲交sv欧美| 一区二区三区四区激情视频| 精品一区二区三卡| 国产精品 国内视频| 精品一区二区三区av网在线观看 | 蜜桃在线观看..| 三级毛片av免费| 九色亚洲精品在线播放| 亚洲av日韩精品久久久久久密| 欧美中文综合在线视频| 女人高潮潮喷娇喘18禁视频| 男人操女人黄网站| 日韩 亚洲 欧美在线| videosex国产| 国产日韩一区二区三区精品不卡| 精品国产一区二区三区四区第35| 菩萨蛮人人尽说江南好唐韦庄| 久热爱精品视频在线9| 一个人免费在线观看的高清视频 | 天天躁狠狠躁夜夜躁狠狠躁| 99精品久久久久人妻精品| 美女国产高潮福利片在线看| 成人亚洲精品一区在线观看| 丰满迷人的少妇在线观看| 男女午夜视频在线观看| 最黄视频免费看| 手机成人av网站| 日韩熟女老妇一区二区性免费视频| 午夜成年电影在线免费观看| 一个人免费看片子| 久久精品亚洲av国产电影网| 水蜜桃什么品种好| 51午夜福利影视在线观看| 黄色怎么调成土黄色| 在线亚洲精品国产二区图片欧美| 亚洲一区二区三区欧美精品| 老司机影院成人| 日韩,欧美,国产一区二区三区| 亚洲精品乱久久久久久| 好男人电影高清在线观看| 男女边摸边吃奶| 王馨瑶露胸无遮挡在线观看| 国产欧美亚洲国产| 欧美精品av麻豆av| 天天躁狠狠躁夜夜躁狠狠躁| 视频区图区小说| 日韩一区二区三区影片| 香蕉丝袜av| 成人影院久久| 国产在线观看jvid| 久久久国产成人免费| 水蜜桃什么品种好| 久久国产精品大桥未久av| 久久影院123| 国产xxxxx性猛交| 欧美人与性动交α欧美精品济南到| 欧美日韩av久久| 免费高清在线观看日韩| 国产成人免费无遮挡视频| 欧美黑人精品巨大| 日韩中文字幕视频在线看片| 日本五十路高清| 国产av又大| 美女国产高潮福利片在线看| 亚洲性夜色夜夜综合| 丰满迷人的少妇在线观看| 精品国产一区二区三区四区第35| 亚洲成国产人片在线观看| 老熟妇仑乱视频hdxx| 亚洲精品自拍成人| 亚洲午夜精品一区,二区,三区| 国产一区二区三区在线臀色熟女 | 黄色视频,在线免费观看| 黄色a级毛片大全视频| 精品国产一区二区三区四区第35| 99热网站在线观看| 亚洲全国av大片| 亚洲精品日韩在线中文字幕| 国产精品99久久99久久久不卡| 久久久国产精品麻豆| 久久热在线av| 国产xxxxx性猛交| 国产成人影院久久av| 两性午夜刺激爽爽歪歪视频在线观看 | 三上悠亚av全集在线观看| 欧美精品人与动牲交sv欧美| 我要看黄色一级片免费的| 成年人免费黄色播放视频| 欧美另类一区| 日韩精品免费视频一区二区三区| 各种免费的搞黄视频| 国产高清视频在线播放一区 | 免费在线观看影片大全网站| 大香蕉久久成人网| 亚洲精品一卡2卡三卡4卡5卡 | 淫妇啪啪啪对白视频 | 日韩大片免费观看网站| 久久人人97超碰香蕉20202| 91国产中文字幕| 国产精品久久久久久精品电影小说| 最近中文字幕2019免费版| 国产99久久九九免费精品| 1024香蕉在线观看| 欧美少妇被猛烈插入视频| 亚洲成人免费av在线播放| 别揉我奶头~嗯~啊~动态视频 | 性少妇av在线| 精品人妻1区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产欧美日韩在线播放| 欧美人与性动交α欧美精品济南到| 国产亚洲av高清不卡| 高清黄色对白视频在线免费看| 9191精品国产免费久久| 欧美激情极品国产一区二区三区| 国产av国产精品国产| 亚洲,欧美精品.| 午夜久久久在线观看| 国产在线观看jvid| 91麻豆av在线| 日韩电影二区| 法律面前人人平等表现在哪些方面 | 桃红色精品国产亚洲av| 日韩中文字幕欧美一区二区| 桃红色精品国产亚洲av| 建设人人有责人人尽责人人享有的| 人成视频在线观看免费观看| 黄色视频在线播放观看不卡| 中亚洲国语对白在线视频| 成年人午夜在线观看视频| 精品国产一区二区三区久久久樱花| 少妇猛男粗大的猛烈进出视频| 十八禁网站免费在线| 伦理电影免费视频| 国产不卡av网站在线观看| 欧美另类亚洲清纯唯美| 日韩视频在线欧美| 国产激情久久老熟女| 国产在线观看jvid| 老司机福利观看| 十八禁高潮呻吟视频| bbb黄色大片| 中文字幕制服av| 亚洲色图 男人天堂 中文字幕| 亚洲中文字幕日韩| 如日韩欧美国产精品一区二区三区| 亚洲欧美日韩另类电影网站| 欧美日韩成人在线一区二区| av免费在线观看网站| 欧美av亚洲av综合av国产av| av又黄又爽大尺度在线免费看| 成年av动漫网址| 在线观看免费午夜福利视频| 成年av动漫网址| 99热国产这里只有精品6| 欧美日韩亚洲高清精品| 美女主播在线视频| 好男人电影高清在线观看| 欧美日韩福利视频一区二区| 欧美+亚洲+日韩+国产| 女警被强在线播放| av天堂在线播放| 亚洲伊人久久精品综合| 天天躁狠狠躁夜夜躁狠狠躁| 人妻久久中文字幕网| 欧美精品高潮呻吟av久久| 欧美黄色淫秽网站| 国产精品麻豆人妻色哟哟久久| 99热网站在线观看| 在线观看一区二区三区激情| 免费久久久久久久精品成人欧美视频| 久久精品国产亚洲av香蕉五月 | e午夜精品久久久久久久| 自线自在国产av| 亚洲欧美成人综合另类久久久| 老熟妇仑乱视频hdxx| 欧美日韩福利视频一区二区| 亚洲精品一区蜜桃| 狂野欧美激情性xxxx| 亚洲欧美色中文字幕在线| 精品久久蜜臀av无| 日韩电影二区| 久久ye,这里只有精品| 欧美午夜高清在线| 久久久久久久久免费视频了| 国产一区二区三区综合在线观看| 91精品伊人久久大香线蕉| 91成年电影在线观看| 9191精品国产免费久久| 黑丝袜美女国产一区| 91九色精品人成在线观看| 国产精品av久久久久免费| 女人爽到高潮嗷嗷叫在线视频| 俄罗斯特黄特色一大片| 免费日韩欧美在线观看| 久久久久久久久久久久大奶| 国产精品国产三级国产专区5o| 嫁个100分男人电影在线观看| 亚洲精品在线美女| 日本一区二区免费在线视频| 亚洲 国产 在线| 免费不卡黄色视频| 免费少妇av软件| 国产激情久久老熟女| 考比视频在线观看| 人人妻人人添人人爽欧美一区卜| 岛国毛片在线播放| 国产免费一区二区三区四区乱码| 亚洲欧美一区二区三区黑人| 亚洲欧美精品综合一区二区三区| 亚洲美女黄色视频免费看| 婷婷成人精品国产| 精品国产一区二区久久| 精品亚洲成国产av| 欧美中文综合在线视频| 欧美黄色片欧美黄色片| 少妇猛男粗大的猛烈进出视频| 国产一区二区三区在线臀色熟女 | av超薄肉色丝袜交足视频| 精品少妇一区二区三区视频日本电影| 亚洲天堂av无毛| 国产欧美日韩一区二区三区在线| 国产97色在线日韩免费| 亚洲精品久久成人aⅴ小说| www.熟女人妻精品国产| 亚洲av美国av| 91精品国产国语对白视频| 国产成人精品在线电影| 老司机影院成人| 亚洲精华国产精华精| 亚洲中文字幕日韩| 如日韩欧美国产精品一区二区三区| 另类精品久久| 操出白浆在线播放| 女人高潮潮喷娇喘18禁视频| 激情视频va一区二区三区| 亚洲免费av在线视频| 啦啦啦 在线观看视频| 99热国产这里只有精品6| 久热爱精品视频在线9| 中亚洲国语对白在线视频| www日本在线高清视频| 日韩视频一区二区在线观看| av国产精品久久久久影院| 国产91精品成人一区二区三区 | 日本撒尿小便嘘嘘汇集6| 欧美日韩精品网址| 高清黄色对白视频在线免费看| 国产老妇伦熟女老妇高清| 精品免费久久久久久久清纯 | 在线精品无人区一区二区三| 黄片播放在线免费| 90打野战视频偷拍视频| 美女国产高潮福利片在线看| 中文欧美无线码| 三上悠亚av全集在线观看| 欧美黄色淫秽网站| 他把我摸到了高潮在线观看 | 亚洲精品国产av蜜桃| 老司机午夜十八禁免费视频| 精品一品国产午夜福利视频| 青青草视频在线视频观看| 午夜福利免费观看在线| 在线观看舔阴道视频| 国产成人精品久久二区二区91| 一级黄色大片毛片| 午夜精品久久久久久毛片777| 午夜免费成人在线视频| 狂野欧美激情性xxxx| 热99久久久久精品小说推荐| 国产1区2区3区精品| 天堂中文最新版在线下载| 一级a爱视频在线免费观看| 香蕉国产在线看| av天堂在线播放| 免费久久久久久久精品成人欧美视频| 精品少妇黑人巨大在线播放| 亚洲欧美清纯卡通| 国产黄色免费在线视频| 999久久久国产精品视频| 欧美黄色淫秽网站| 国产一卡二卡三卡精品| 真人做人爱边吃奶动态| 欧美在线一区亚洲| 12—13女人毛片做爰片一| 亚洲av男天堂| 欧美午夜高清在线| 99久久综合免费| 亚洲av日韩在线播放| 亚洲精品第二区| 久久精品久久久久久噜噜老黄| 久久久国产欧美日韩av| 伦理电影免费视频| 久久女婷五月综合色啪小说| a级毛片黄视频| a级毛片在线看网站| 国产高清视频在线播放一区 | 在线观看免费视频网站a站| 色视频在线一区二区三区| 丝袜美腿诱惑在线| 97精品久久久久久久久久精品| 亚洲人成电影免费在线| 日韩 欧美 亚洲 中文字幕| 国产成人av激情在线播放| 老司机影院毛片| 国产精品1区2区在线观看. | 久久精品国产亚洲av香蕉五月 | 午夜成年电影在线免费观看| 国产一区二区激情短视频 | 欧美日韩视频精品一区| 国产精品麻豆人妻色哟哟久久| 久久影院123| 日韩视频一区二区在线观看| av不卡在线播放| 欧美亚洲 丝袜 人妻 在线| 丝瓜视频免费看黄片| 国产亚洲午夜精品一区二区久久| 大香蕉久久网| 水蜜桃什么品种好| 18禁黄网站禁片午夜丰满| 日韩精品免费视频一区二区三区| 窝窝影院91人妻| 国产成人精品久久二区二区91| 美女视频免费永久观看网站| 成人三级做爰电影| 亚洲精品美女久久久久99蜜臀| 国产精品久久久久成人av| 亚洲人成电影观看| 国产日韩欧美亚洲二区| 丰满迷人的少妇在线观看| 久久中文字幕一级| 亚洲精品在线美女| 久久久久国产一级毛片高清牌| 亚洲男人天堂网一区| 亚洲专区字幕在线| 成年人黄色毛片网站| 少妇精品久久久久久久| 久久久久久久精品精品| 99国产精品一区二区蜜桃av | 国产精品国产三级国产专区5o| 色婷婷av一区二区三区视频| 亚洲,欧美精品.| 亚洲专区字幕在线| 后天国语完整版免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 日本av手机在线免费观看| 亚洲国产精品一区二区三区在线| 欧美日韩亚洲高清精品| √禁漫天堂资源中文www| 久久久国产精品麻豆| 久久久久国产精品人妻一区二区| 99久久精品国产亚洲精品| 久久亚洲国产成人精品v| 亚洲欧美清纯卡通| 成人黄色视频免费在线看| 久久国产精品男人的天堂亚洲| av国产精品久久久久影院| 纵有疾风起免费观看全集完整版| 好男人电影高清在线观看| 高清在线国产一区| 丝袜美腿诱惑在线| 999精品在线视频| 一级片'在线观看视频| 黄色视频不卡| 亚洲专区国产一区二区| 夜夜夜夜夜久久久久| 久久中文看片网| 我要看黄色一级片免费的| 一区在线观看完整版| 精品一品国产午夜福利视频| 王馨瑶露胸无遮挡在线观看| 国产精品偷伦视频观看了| 在线观看一区二区三区激情| 中文字幕av电影在线播放| 久久天躁狠狠躁夜夜2o2o| 乱人伦中国视频| 9色porny在线观看| 国产成人精品无人区| 久久人人97超碰香蕉20202| 国产人伦9x9x在线观看| 久久人妻熟女aⅴ| 久久久久久人人人人人| 男女国产视频网站| 免费观看人在逋| 亚洲精品中文字幕在线视频| 日本猛色少妇xxxxx猛交久久| 一个人免费在线观看的高清视频 | 电影成人av| 日韩欧美国产一区二区入口| 中国国产av一级| 中文字幕另类日韩欧美亚洲嫩草| 飞空精品影院首页| √禁漫天堂资源中文www| 俄罗斯特黄特色一大片| 黄色 视频免费看| 法律面前人人平等表现在哪些方面 | 国产一区二区三区av在线| 国产精品香港三级国产av潘金莲| 美女福利国产在线| 亚洲五月婷婷丁香| 搡老岳熟女国产| 精品国产乱子伦一区二区三区 | 亚洲五月色婷婷综合| 亚洲免费av在线视频| 久久久水蜜桃国产精品网| 老司机在亚洲福利影院| 桃花免费在线播放| 亚洲国产日韩一区二区| 九色亚洲精品在线播放| 18禁裸乳无遮挡动漫免费视频| 国产在线免费精品| 少妇 在线观看| 大码成人一级视频| 日韩熟女老妇一区二区性免费视频| 国产精品一二三区在线看| 美女扒开内裤让男人捅视频| 叶爱在线成人免费视频播放|