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

    基于稀疏重構(gòu)的L型陣列MIMO雷達(dá)降維DOA估計(jì)

    2015-06-05 15:33:40浩,崔琛,代林,余

    梁 浩,崔 琛,代 林,余 劍

    (電子工程學(xué)院通信對(duì)抗系,安徽合肥230037)

    基于稀疏重構(gòu)的L型陣列MIMO雷達(dá)降維DOA估計(jì)

    梁 浩,崔 琛,代 林,余 劍

    (電子工程學(xué)院通信對(duì)抗系,安徽合肥230037)

    針對(duì)L型陣列多輸入多輸出(multiple-input multiple-output,MIMO)雷達(dá)二維空間角估計(jì)問題,提出一種基于協(xié)方差矩陣聯(lián)合稀疏重構(gòu)的降維波達(dá)方向(direction of arrival,DOA)估計(jì)算法。該算法根據(jù)L型陣列MIMO雷達(dá)聯(lián)合流型矢量的特點(diǎn),通過降維矩陣的設(shè)計(jì)及回波數(shù)據(jù)的降維變換,最大程度地去除了所有的冗余數(shù)據(jù);通過協(xié)方差矩陣聯(lián)合構(gòu)造稀疏線性模型,將2維角參量空間映射到1維空間,極大降低字典長度和求解復(fù)雜度的同時(shí),不犧牲陣列孔徑,實(shí)現(xiàn)了二維空間角度的有效估計(jì)和參數(shù)的自動(dòng)配對(duì)。理論分析與實(shí)驗(yàn)仿真表明:與RD_MUSIC算法相比,本文降維處理有效提高陣元利用率的同時(shí),最大程度地降低了回波數(shù)據(jù)的維數(shù);與傳統(tǒng)子空間類算法相比,基于協(xié)方差矩陣聯(lián)合構(gòu)造的稀疏線性模型充分利用了陣列孔徑,無需預(yù)先估計(jì)目標(biāo)數(shù)目,參數(shù)估計(jì)性能在低信噪比及小快拍數(shù)據(jù)長度下優(yōu)勢明顯。最后,仿真結(jié)果驗(yàn)證了本文理論分析的正確性和算法的有效性。

    多輸入多輸出雷達(dá);L型陣列;稀疏重構(gòu);波達(dá)方向估計(jì)

    0 引 言

    多輸入多輸出(multiple-input multiple-output,MIMO)雷達(dá)[1-2]能有效提高雷達(dá)參數(shù)識(shí)別性能[3]及目標(biāo)空間分辨力;同時(shí)因其具有數(shù)字化、多通道化、低輻射等特點(diǎn),已成為現(xiàn)代新型雷達(dá)發(fā)展趨勢的綜合體現(xiàn)和典型代表。目標(biāo)角度估計(jì)[25]是雷達(dá)信號(hào)處理中一個(gè)重要的研究方向,目前對(duì)單基地MIMO雷達(dá)的研究大多僅限于一維參數(shù)估計(jì);事實(shí)上,當(dāng)收、發(fā)陣列均采用二維天線構(gòu)型時(shí),意味著目標(biāo)參數(shù)維度的擴(kuò)展,對(duì)目標(biāo)的描述也就更準(zhǔn)確,相應(yīng)的參數(shù)估計(jì)也就更為復(fù)雜;在眾多的平面天線構(gòu)型中,L型陣列結(jié)構(gòu)簡單、陣列冗余度較小,能夠?qū)崿F(xiàn)空域目標(biāo)二維角度定位,且已被證明具有優(yōu)于其他交叉陣列結(jié)構(gòu)的波達(dá)方向(direction of arrival,DOA)估計(jì)性能[67];因此深入研究L型陣列MIMO雷達(dá)的參數(shù)估計(jì)問題具有重要意義和實(shí)用價(jià)值。

    關(guān)于L型陣列MIMO雷達(dá)的研究主要集中在傳統(tǒng)子空間類算法的擴(kuò)展應(yīng)用。文獻(xiàn)[8- 9]基于單基地MIMO雷達(dá)研究了不同平面流型配置下陣列整體的流型敏感性以及陣元的重要程度,為MIMO雷達(dá)系統(tǒng)的天線設(shè)計(jì)提供了依據(jù)。文獻(xiàn)[10]基于DOA矩陣思想,實(shí)現(xiàn)了目標(biāo)二維角度的估計(jì),并進(jìn)一步提出聯(lián)合對(duì)角化DOA矩陣方法解決了角度兼并問題;但該模型采用收發(fā)陣列垂直分置的L型配置,因此算法并不適用于收發(fā)均為L型陣列的共置MIMO雷達(dá)。文獻(xiàn)[11]建立了L型陣列配置的單基地MIMO雷達(dá),并基于MIMO-Capon算法實(shí)現(xiàn)了目標(biāo)方位角和俯仰角的二維估計(jì),但需要二維的譜搜索;文獻(xiàn)[12]針對(duì)文獻(xiàn)[11]算法計(jì)算復(fù)雜度較高的問題,提出一種基于MUSIC算法的L型陣列多輸入多輸出雷達(dá)降維DOA估計(jì)算法。算法通過構(gòu)造降維矩陣對(duì)回波數(shù)據(jù)進(jìn)行降維預(yù)處理后,利用二次優(yōu)化方法將二維DOA估計(jì)分解為兩個(gè)一維DOA估計(jì),一定程度上降低了運(yùn)算復(fù)雜度,但降維矩陣的設(shè)計(jì)并沒有最大程度地降低回波數(shù)據(jù)的維數(shù),沒有去除所有重復(fù)的虛擬陣元,回波數(shù)據(jù)中仍存在冗余;此外在利用二次優(yōu)化進(jìn)行降維求解過程中,對(duì)方向向量中各元素的約束較弱[13],造成估計(jì)精度較差,同時(shí)協(xié)方差矩陣的構(gòu)建和特征分解以及兩次一維譜搜索仍存在較高的運(yùn)算量。

    盡管傳統(tǒng)子空間類算法能夠?qū)崿F(xiàn)二維空間角的有效估計(jì),但普遍存在以下問題:①需要預(yù)先進(jìn)行目標(biāo)數(shù)目的估計(jì),其準(zhǔn)確性會(huì)嚴(yán)重影響后續(xù)子空間的劃分及算法的性能;②需要多次脈沖的數(shù)據(jù)來提高與理論協(xié)方差矩陣的逼近程度,在低信噪比及有限數(shù)據(jù)長度下,對(duì)信號(hào)子空間的估計(jì)會(huì)存在較大偏差,使得估計(jì)精度和穩(wěn)健性能嚴(yán)重惡化。利用遠(yuǎn)場目標(biāo)的空域稀疏性,基于稀疏重構(gòu)理論能夠很好的解決以上問題,文獻(xiàn)[14- 15]利用L型陣列的流型結(jié)構(gòu),通過二維稀疏線性模型分解和重新組合,利用其構(gòu)造冗余字典可以將方位角和俯仰角的組合從2維空間映射到1維空間,極大降低字典的長度和求解的復(fù)雜度;文獻(xiàn)[16]更是將這種思想推廣到了頻域,實(shí)現(xiàn)了寬頻帶內(nèi)目標(biāo)方位與頻率的聯(lián)合估計(jì)。以上算法本質(zhì)上是基于陣列流型,將二維參量解耦合,分維兩次基于傳統(tǒng)一維稀疏重構(gòu)算法實(shí)現(xiàn)二維參量的估計(jì);當(dāng)MIMO雷達(dá)收發(fā)均為L型陣列時(shí),這類算法僅能利用到部分虛擬陣元,造成孔徑的損失,這顯然與MIMO雷達(dá)通過虛擬孔徑提高分辨力的初衷是背馳的。本文針對(duì)L型陣列MIMO雷達(dá)二維空間角估計(jì)問題,提出一種降維預(yù)處理后基于協(xié)方差矩陣聯(lián)合稀疏重構(gòu)(reduceddimensional joint sparse recovery of covariance matrix,RD_ JSRCM)的L型陣列MIMO雷達(dá)DOA估計(jì)算法。

    1 問題建模

    圖1 L型陣列MIMO雷達(dá)結(jié)構(gòu)及角度配置關(guān)系

    遠(yuǎn)場空域內(nèi)存在K個(gè)目標(biāo),第k(k=1,2,…,K)個(gè)目標(biāo)空間角為(?k,φk),滿足cos?=cosθcosφ,cosφ=cosθsinφ;其中(θ,φ)對(duì)應(yīng)目標(biāo)的俯仰角與方位角,(?,φ)對(duì)應(yīng)目標(biāo)分別與X、Y軸的夾角。則第q{q=1,…,Q}次脈沖下接收端的輸出為

    對(duì)yq進(jìn)行列堆棧,即可得到Q次脈沖下L型陣列MIMO雷達(dá)的回波接收數(shù)據(jù):

    式中,η=[β1,…,βQ]對(duì)應(yīng)目標(biāo)的散射系數(shù);N=[n1,…,nQ]

    2 算法描述

    2.1 接收數(shù)據(jù)降維預(yù)處理

    L型陣列MIMO雷達(dá)虛擬擴(kuò)展如圖2所示,顯然Mx+My-1個(gè)實(shí)陣元經(jīng)過MIMO雷達(dá)的虛擬擴(kuò)展后,得到(Mx+My-1)2個(gè)虛擬陣元,然而其中僅有MxMy+Mx+My-2為有效陣元,其余均為冗余陣元。則由信號(hào)模型可得

    圖2 虛擬擴(kuò)展及示意圖

    式中

    對(duì)應(yīng)置換矩陣,構(gòu)造方式與文獻(xiàn)[12]類似;顯然,g(?,φ)與ρ(?,φ)成線性關(guān)系,同時(shí)

    則存在以下線性變換

    式中,Πxy為置換矩陣,與Πx、Πy類似。由式(5)、式(7)和式(8)可得,存在降維矩陣Γ滿足

    則式(4)中回波接收數(shù)據(jù)可表示為

    式(15)證明略。從式(14)、式(15)可得,降維后回波數(shù)據(jù)Z等效長度MxMy+Mx+My-2的加權(quán)平面陣的回波信號(hào),權(quán)值為(ΓHΓ)1/2的對(duì)角元素。由H(?,φ)及圖2可知,本文通過設(shè)計(jì)降維矩陣,將(Mx+My-1)2維的回波信號(hào)降至MxMy+Mx+My-2,去除了所有的重復(fù)陣元,而文獻(xiàn)[12]降維多重信號(hào)分類(reduced-dimensional multiple signal classification,RD_MUSIC)RD_MUSIC算法僅僅是將(Mx+My-1)2維的雷達(dá)回波信號(hào)降至MxMy+Mx+2My-4,降維后仍存在My-2個(gè)重復(fù)陣元。因此與之相比,本文降維預(yù)處理最大程度地降低了回波數(shù)據(jù)的維數(shù),有效去除了原始回波數(shù)據(jù)中所有重復(fù)量,達(dá)到了降維的目的。

    2.2 對(duì)稱置換預(yù)處理

    為后續(xù)處理方便,首先對(duì)降維后的回波數(shù)據(jù)Z進(jìn)行權(quán)值歸一化操作

    在得到預(yù)處理的回波數(shù)據(jù)Zx和Zy后,即可利用子空間類算法(如文獻(xiàn)[12]RD_MUSIC、2D_MUSIC等)獲取對(duì)信號(hào)子空間H(?,φ)的估計(jì),實(shí)現(xiàn)二維空間角的估計(jì)。考慮到子空間類算法在低信噪比及小數(shù)據(jù)快拍長度下估計(jì)性能較差的缺陷,利用遠(yuǎn)場空域目標(biāo)稀疏的先驗(yàn)信息,將DOA估計(jì)轉(zhuǎn)化為從測量矢量中重構(gòu)稀疏信號(hào)的問題,能夠獲得目標(biāo)方位信息的超分辨估計(jì)。直接構(gòu)建二維過完備字典,則意味著后續(xù)高維復(fù)雜度的迭代重構(gòu);由圖2可知,經(jīng)過MIMO雷達(dá)虛擬后,陣列仍呈現(xiàn)為近似L型,通過二維稀疏線性模型分解和重新組合,利用其構(gòu)造冗余字典可以將二維空間映射到一維空間,但問題是這樣只能利用MIMO雷達(dá)的部分虛擬陣元(只利用到X軸上的2Mx-1陣元和Y軸上的2My-1陣元),降低了雷達(dá)整體的孔徑利用率,從而不能很好發(fā)揮MIMO雷達(dá)的高分辨性能;為此本文利用預(yù)處理后的回波數(shù)據(jù)Zx和Zy構(gòu)建相應(yīng)的稀疏線性模型,提出一種基于協(xié)方差矩陣聯(lián)合稀疏重構(gòu)的L型陣列MIMO雷達(dá)降維DOA估計(jì)算法。

    2.3 聯(lián)合稀疏模型建立

    式中,Vx=diag{Wx}。對(duì)Rxx進(jìn)行列向量化操作得

    式中,k=1,…,K;n=1,…,Nx。則式(29)可以重寫為

    同樣可以利用Zy構(gòu)建角參量φ的稀疏線性模型

    2.4 模型求解

    式(33)和式(34)中,由于ex與ey均服從高斯分布,因此可以采用2范數(shù)來約束擬合誤差,則稀疏線性模型的目標(biāo)優(yōu)化函數(shù)可表示為:

    λ對(duì)應(yīng)正則化參數(shù),平衡控制著擬合誤差與μ稀疏性,可以通過L-Curce、廣義的SURE(generation SURE,GSURE)等方法進(jìn)行調(diào)整和選取。直接求解式(35)是困難的,一方面考慮擬合誤差服從高斯分布,可以采用文獻(xiàn)[19]中的平滑零范數(shù),通過一類高斯函數(shù)近似逼近零范數(shù)求解最小化問題,重構(gòu)出相應(yīng)的稀疏信號(hào),但該方法重構(gòu)性能以及精度與高斯函數(shù)和映射函數(shù)的選取相關(guān);另一方面通過凸松弛類算法采用1范數(shù)來代替0范數(shù)進(jìn)行凸松弛。本文考慮凸松弛類算法將非凸優(yōu)化問題轉(zhuǎn)化為以下的凸優(yōu)化問題

    可以分別利用內(nèi)點(diǎn)法求解,考慮式(36)中參數(shù)均為復(fù)數(shù),則對(duì)應(yīng)的二階錐規(guī)劃形式可表示為

    顯然,關(guān)于角參量?和φ的求解是分開進(jìn)行的,得到對(duì)應(yīng)的估計(jì)值和后,定義低維的聯(lián)合過完備字典D()為

    可得重組后的聯(lián)合稀疏線性模型為

    2.5 算法復(fù)雜度分析

    本文降維預(yù)處理后采用文獻(xiàn)[15]分維稀疏重構(gòu)算法簡稱為RD_JEVSD;本文降維預(yù)處理后基于協(xié)方差矩陣聯(lián)合稀疏重構(gòu)的方法稱為RD_JSRCM算法。RD_JEVSD與RD_JSRCM算法在后續(xù)基于凸松弛類算法求解過程中,算法復(fù)雜度相差不大,盡管兩者構(gòu)建的過完備字典的維數(shù)不同,但都要遠(yuǎn)小于完備字典的長度Nx和Ny,因此算法復(fù)雜度主要為O{N+N};而在稀疏線性模型構(gòu)建過程中,RD _JEVSD算法僅利用到(L型陣列MIMO雷達(dá)虛擬擴(kuò)展后的陣列)X、Y軸的等效回波數(shù)據(jù),算法運(yùn)算量集中在協(xié)方差矩陣計(jì)算及特征分解:O{Q(2Mx-1)2+Q(2My-1)2+(2Mx-1)3+(2My-1)3};對(duì)應(yīng)的RD_JSRCM算法運(yùn)算量集中在zx和zy構(gòu)造中協(xié)方差矩陣的計(jì)算:O{Q(2Mx-1)2+Q(2My-1)2+Q(Mx-1)2+Q(My-1)2},因此兩者的運(yùn)算復(fù)雜度相當(dāng),但本文RD_JSRCM算法無需特征分解,不犧牲雷達(dá)孔徑,最大程度地利用了虛擬陣列,因此在估計(jì)性能及穩(wěn)健性方面更優(yōu)。

    3 仿真實(shí)驗(yàn)與分析

    假設(shè)均勻配置下的L型陣列MIMO雷達(dá),陣元間距滿足dx=dy=λ/2,以多相碼為正交發(fā)射信號(hào),分別進(jìn)行以下實(shí)驗(yàn)。

    實(shí)驗(yàn)1 算法的有效性驗(yàn)證。不失一般性,假設(shè)陣元數(shù)滿足Mx=My=8,遠(yuǎn)場空域存在K=3個(gè)目標(biāo),與L型陣列MIMO雷達(dá)X、Y軸的夾角為(50°,60°),(70°,100°),(110°,120°);信噪比為5 d B,快拍數(shù)Q=200;進(jìn)行100次蒙特卡羅實(shí)驗(yàn),驗(yàn)證本文RD_JSRCM算法的有效性,其中過完備字典Dx(?′)和Dy(φ′)均在(0°,180°)內(nèi)以0.1°為間隔均勻采樣構(gòu)建而成,仿真結(jié)果如圖3所示。圖3(a)為100次蒙特卡羅實(shí)驗(yàn)后本文RD_JSRCM算法與文獻(xiàn)RD_JEVSD算法的二維空間角的估計(jì)結(jié)果;圖3(b)為本文RD_JSRCM算法二維空間角經(jīng)過自動(dòng)配對(duì)后的星座圖。

    圖3 目標(biāo)二維空間參量?和φ的估計(jì)結(jié)果

    由圖3可知,本文所提RD_JSRCM算法能夠?qū)崿F(xiàn)對(duì)目標(biāo)二維空間角的有效估計(jì),且能實(shí)現(xiàn)參數(shù)的自動(dòng)配對(duì);同時(shí)由圖3也可以看出,與文獻(xiàn)[15]RD_JEVSD算法相比,本文RD_JSRCM算法空間譜估計(jì)更為準(zhǔn)確和穩(wěn)健,空間角度位置星座圖上的估計(jì)結(jié)果比較集中而沒有出現(xiàn)散布,一定程度上也反映了本文算法的穩(wěn)健性。

    實(shí)驗(yàn)2 算法的估計(jì)性能比較。假設(shè)發(fā)射陣元數(shù)滿足Mx=My=6,遠(yuǎn)場空域目標(biāo)位置及個(gè)數(shù)不變;快拍數(shù)Q=50,信噪比-10~30 dB,比較2D_MUSIC、RD_MUSIC、RD_JEVSD以及本文RD_JSRCM算法的估計(jì)性能,其中2D_MUSIC與文獻(xiàn)[12]RD_MUSIC算法的譜搜索步長均為0.1°,RD_JEVSD及本文RD_JSRCM算法的過完備字典構(gòu)造均以0.1°為間隔采樣,仿真結(jié)果如圖4所示。

    圖4 目標(biāo)二維空間角估計(jì)的RMSE與SNR的關(guān)系

    由圖4可以看出,隨著信噪比(signal-to-noise ratio,SNR)的增大,4種算法對(duì)目標(biāo)二維空間角估計(jì)的均方根誤差(root mean square error,RMSE)都逐漸變小,且信噪比越大,估計(jì)精度越高,這一點(diǎn)很好理解;同時(shí)與子空間類算法(2D_MUSIC和RD_MUSIC算法)相比,基于稀疏重構(gòu)理論的RD_JEVSD和RD_JSRCM算法在低信噪比時(shí),估計(jì)性能較優(yōu),高信噪比時(shí)角度估計(jì)性能趨于一致;與RD_JEVSD算法僅利用到X、Y軸(L型陣列MIMO雷達(dá)虛擬擴(kuò)展后的陣列)的等效回波數(shù)據(jù)相比較,本文提出的RD_JSRCM算法基于協(xié)方差矩陣聯(lián)合構(gòu)建稀疏模型,利用到了所有陣元回波數(shù)據(jù),最大程度地利用了陣列孔徑,具有更高的估計(jì)性能。此外本文RD_JSRCM算法通過降維處理,去除所有重復(fù)項(xiàng)的同時(shí),并沒有改變?cè)肼曁匦?,保留了所有的目?biāo)信息,在利用降維處理后回波數(shù)據(jù)構(gòu)建聯(lián)合稀疏模型時(shí),權(quán)值歸一化操作,一定程度上也降低了噪聲的影響,提高了參數(shù)估計(jì)性能。

    實(shí)驗(yàn)3 算法的性能與參數(shù)之間的關(guān)系。空域目標(biāo)位置、個(gè)數(shù)與實(shí)驗(yàn)2相同,快拍數(shù)Q=50,陣元數(shù)Mx=My=4~12時(shí),信噪比為5 dB,比較不同陣元配置本文RD_JSRCM算法與文獻(xiàn)RD_JEVSD算法的估計(jì)性能,仿真結(jié)果如圖5所示。陣元數(shù)滿足Mx=My=6,Q=50~500變化,比較不同快拍數(shù)據(jù)下本文RD_JSRCM算法與文獻(xiàn)RD_JEVSD算法的估計(jì)性能,仿真結(jié)果如圖6所示。由圖5、圖6可以看出,隨著陣元數(shù)Mx,My的增大,對(duì)應(yīng)陣列雷達(dá)孔徑也逐漸變大,因此兩種算法對(duì)二維空間角的估計(jì)性能也越好;隨著快拍數(shù)Q的增大,統(tǒng)計(jì)意義上稀疏模型的構(gòu)建也就越準(zhǔn)確,重構(gòu)精度也就越高,因此角估計(jì)性能也就更好。同時(shí)由圖5、圖6的仿真結(jié)果可以看出,由于本文RD_ JSRCM算法利用了MIMO雷達(dá)所有虛擬陣元,沒有犧牲陣列孔徑,因此估計(jì)性能要優(yōu)于文獻(xiàn)RD_JEVSD算法。

    圖5 二維空間角估計(jì)的RMSE隨陣元數(shù)的變化

    圖6 二維空間角估計(jì)的RMSE隨快拍數(shù)Q的變化

    4 結(jié) 論

    本文針對(duì)L型陣列MIMO雷達(dá)二維空間角估計(jì)問題,提出一種基于協(xié)方差矩陣聯(lián)合稀疏重構(gòu)的L型陣列MIMO雷達(dá)降維DOA估計(jì)算法。理論分析與仿真表明:①與文獻(xiàn)RD_MUSIC算法相比,本文算法通過降維矩陣的設(shè)計(jì)及回波數(shù)據(jù)的降維變換,最大程度地降低了回波數(shù)據(jù)的維數(shù),進(jìn)一步降低了所需處理的回波的數(shù)據(jù)量,有效提高了陣元的利用率;②基于協(xié)方差矩陣聯(lián)合的稀疏線性模型,將二維角參量空間映射到一維空間,極大降低字典長度和求解復(fù)雜度的同時(shí),充分利用陣列孔徑,實(shí)現(xiàn)了二維參量的聯(lián)合估計(jì)及參數(shù)的自動(dòng)配對(duì);③與傳統(tǒng)子空間類算法相比,本文算法無需預(yù)先估計(jì)目標(biāo)數(shù)目,估計(jì)性能在低信噪比及小快拍數(shù)據(jù)長度下優(yōu)勢明顯。

    [1]Huleihel W,Tabrikian J,Shavit R.Optimal adaptive waveform design for cognitive MIMO radar[J].IEEE Trans.on Signal Processing,2013,61(20):5075- 5089.

    [2]Tang B,Tang J,Zhang Y,et al.Maximum likelihood estimation of DOD and DOA for bistatic MIMO radar[J].Signal Processing,2013,93(5):1349- 1357.

    [3]Arash K.Efficient transmit beamspace design for search-free based DOA estimation in MIMO radar[J].IEEE Trans.on Signal Processing,2014,62(6):1490- 1500.

    [4]Li J F,Zhang X F.A method for joint angle and array gain-phase error estimation in bistatic multiple input multiple output nonlinear arrays[J].IET Signal Processing,2014,8(2):131- 137.

    [5]Jiang H,Zhang Y,Li J,et al.A parafac-based algorithm for multi-dimensional parameter estimation in polarimetric bistatic MIMO radar[J].EURASIP Journal on Advances in Signal Processing,2013,2013(1):1- 14.

    [6]Cheng Q,Hua Y.Further study of the pencil-MUSICalgorithm[J].IEEE Trans.on Aerospace Electronic Systems,1996,32(1):284- 299.

    [7]Hua Y,Sarkar K,Weiner D.An L-shaped array for estimation 2-D directions of wave arrival[J].IEEE Trans.on Antennas and Propagation,1991,39(2):143- 146.

    [8]Chen H W,Li X,Jiang W D,et al.MIMO radar sensitivity analysis of antenna position for direction finding[J].IEEE Trans.on Signal Processing,2012,60(10):5201- 5216.

    [9]Chen H W,Zhou W,Yang J,et al.Manifold sensitivity analysis for MIMO radar[J].IEEE Geoscience and Remote Sensing Letters,2012,9(5):999- 1003.

    [10]Fu W B,Zhao Y B,Su T,et al.DOA matrix method based on MIMO radar with L-shape arrays[J].Journal of Systems Engineering and Electronics,2011,33(11):2398- 2403.(符渭波,趙永波,蘇濤,等.基于L型陣列MIMO雷達(dá)的DOA矩陣方法[J].系統(tǒng)工程與電子技術(shù),2011,33(11):2398- 2403.)

    [11]Xie R,Liu Z,Liu Y F.Multi-target identification and localization in MIMO radar with L-shape arrays[J].Systems Engineering and Electronics,2010,32(1):49- 52.(謝榮,劉錚,劉韻佛.基于L型陣列MIMO雷達(dá)的多目標(biāo)分辨和定位[J].系統(tǒng)工程與電子技術(shù),2010,32(1):49- 52.)

    [12]Wang W,Wang X M,Li X,et al.Reduced-dimensional DOA estimation based on MUSIC algorithm in MIMO radar with L-shaped array[J].Journal of Electronics&Information Technology,2014,36(8):1954- 1959.(王偉,王曉萌,李欣,等.基于MUSIC算法的L型陣列MIMO雷達(dá)降維DOA估計(jì)[J].電子與信息學(xué)報(bào),2014,36(8):1954- 1959.)

    [13]Cai J J,Bao D,Li P,et al.Two-dimensional DOA estimation using reduced-dimensional MUSIC algorithm with strong-constraint optimization[J].Journal of Electronics&Information Technology,2014,36(5):1113- 1118.(蔡晶晶,鮑丹,李鵬,等.強(qiáng)約束優(yōu)化降維MUSIC二維DOA估計(jì)[J].電子與信息學(xué)報(bào),2014,36(5):1113- 1118.)

    [14]Wang L B,Cui C.L-shaped array-based two dimensional direction-of-arrival estimation in a sparse reconstruction framework[J].Journal of Astronautics,2012,33(7):964- 970.(王粒賓,崔琛.基于稀疏重構(gòu)的L型陣列的二維波達(dá)方向估計(jì)[J].宇航學(xué)報(bào),2012,33(7):964- 970.)

    [15]Cui C,W L B.L-shaped array-based joint elevation and azimuth direction finding by spatial sparsity[J].Journal of Circuits and Systems,2013,18(1):297- 303.(崔琛,王粒賓.利用空域稀疏性的L型陣下二維波達(dá)方向估計(jì)[J].電路與系統(tǒng)學(xué)報(bào),2013,18(1):297- 303.)

    [16]Li P F,Zhang M,Zhong Z F.Two dimensional DOA estimation based on sparse representation of space angle[J].Journal of Electronics&Information Technology,2011,33(10):2402- 2406.(李鵬飛,張旻,鐘子發(fā).基于空間角稀疏表示的二維DOA估計(jì)[J].電子與信息學(xué)報(bào),2011,33(10):2402 -2406.)

    [17]Zhang X,Xu D.Low-complexity ESPRIT-based DOA estimation for colocated MIMO radar using reduced-dimension transformation[J].Eletronics Letters,2011,17(47):283- 284.

    [18]Liao B,Chan S C.Direction finding with partly calibrated uniform linear arrays[J].IEEE Trans.on Antennas and Propagation,2012,60(2):922- 929.

    [19]Hyder M M,Mahata K.An improved smoothed L0 approximation algorithm for sparse representation[J].IEEE Trans.on Signal Processing,2010,58(4):2194- 2205.

    Reduced-dimensional DOA estimation based on sparse reconstruction in MIMO radar with L-shaped array

    LIANG Hao,CUI Chen,DAI Lin,YU Jian
    (Department of Communication Countermeasure,Electronic Engineering Institute,Hefei 230037,China)

    Aiming at the problem of two dimensional angles estimation for multiple-input multiple-output(MIMO)radar with L-shaped array,a new reduced-dimensional direction of arrival(DOA)estimation method based on sparse reconstruction is proposed.Giving the steering vector of MIMO radar with L-shaped array,a reduced-dimensional matrix is employed,and data redundancy of high dimensional

    data at the greatest degree can be removed via the reduced-dimensional transformation.Through the joint construction of the twodimensional sparse linear model with covariance matrix,the dimension of the dictionary is reduced to one-dimension from two-dimensional space,and the length of the redundant dictionary and computation complexity is largely reduced.Furthermore,the method,without costing the aperture of array,can realize two dimensional spatial angles estimation with automatic pairing.Compared with reduced-dimensional(RD)MUSIC,the proposed method can reduce the dimension of received data at the greatest degree and enhance sensors efficiency.Compared with the traditional subspace algorithms,the proposed method,which is based on the joint sparse linear model of the covariance matrix,makes the best of all apertures of array and can achieve better estimation performance under lower signal-noise-ratio(SNR)and a few snapshots without pre-estimation for the number of targets.Finally,simulation results verify the correctness of the theoretical analysis and the effect of the proposed algorithm.

    multiple-input multiple-output(MIMO)radar;L-shaped array;sparse reconstruction;direction of arrival(DOA)estimation

    TN 958

    A

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

    梁 浩(198-7- ),男,博士研究生,主要研究方向?yàn)殛嚵行盘?hào)處理以及MIMO雷達(dá)信號(hào)處理。

    E-mail:lhmailhappy@163.com

    崔 琛(196-2- ),男,教授,博士研究生導(dǎo)師,主要研究方向?yàn)槔走_(dá)信號(hào)處理。

    E-mail:kycuichen@163.com

    代 林(198-6- ),男,博士研究生,主要研究方向?yàn)閴嚎s感知雷達(dá)信號(hào)處理。

    E-mail:dailin513@163.com

    余 劍(1980- ),男,講師,博士,主要研究方向?yàn)槔走_(dá)信號(hào)處理、射頻信號(hào)處理。

    E-mail:yujianeei@163.com

    1001-506X(2015)12-2725-08

    2014- 12- 05;

    2015- 04- 15;網(wǎng)絡(luò)優(yōu)先出版日期:2015- 07- 07。

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

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

    1024视频免费在线观看| 久久精品国产a三级三级三级| 日韩欧美在线二视频 | 欧美日韩亚洲高清精品| 一区二区三区激情视频| 国产成人啪精品午夜网站| 午夜久久久在线观看| 亚洲精品av麻豆狂野| netflix在线观看网站| 又黄又粗又硬又大视频| 久久久精品区二区三区| 亚洲片人在线观看| 日本欧美视频一区| 亚洲精品美女久久av网站| 男女下面插进去视频免费观看| 国产色视频综合| 可以免费在线观看a视频的电影网站| 国产日韩一区二区三区精品不卡| av不卡在线播放| 国产欧美日韩一区二区三区在线| 欧美黄色片欧美黄色片| 亚洲av日韩在线播放| 精品一区二区三区四区五区乱码| 国产在线精品亚洲第一网站| 欧美精品亚洲一区二区| 国产在线观看jvid| 91成年电影在线观看| 黄片大片在线免费观看| 成熟少妇高潮喷水视频| 国产亚洲精品一区二区www | netflix在线观看网站| 国产99白浆流出| 国产又爽黄色视频| 国产一区二区三区在线臀色熟女 | 多毛熟女@视频| 无限看片的www在线观看| 999久久久国产精品视频| 在线免费观看的www视频| 99国产精品一区二区蜜桃av | 亚洲精品久久成人aⅴ小说| 俄罗斯特黄特色一大片| 一进一出抽搐动态| 亚洲精品久久午夜乱码| 成人特级黄色片久久久久久久| 国产淫语在线视频| 一夜夜www| 国产麻豆69| 夜夜躁狠狠躁天天躁| 精品少妇一区二区三区视频日本电影| 啦啦啦免费观看视频1| 极品教师在线免费播放| 很黄的视频免费| 国产精品电影一区二区三区 | 夜夜夜夜夜久久久久| 啦啦啦在线免费观看视频4| 亚洲五月色婷婷综合| 精品国产一区二区三区四区第35| 少妇裸体淫交视频免费看高清 | 少妇的丰满在线观看| 热re99久久国产66热| 婷婷成人精品国产| 一本一本久久a久久精品综合妖精| 一边摸一边抽搐一进一小说 | 亚洲精品美女久久av网站| 国产精华一区二区三区| 中文字幕最新亚洲高清| 国产亚洲精品久久久久久毛片 | 色婷婷久久久亚洲欧美| 超碰成人久久| 丰满迷人的少妇在线观看| 1024香蕉在线观看| 高潮久久久久久久久久久不卡| 欧美精品一区二区免费开放| www.999成人在线观看| 精品国产美女av久久久久小说| 成年版毛片免费区| 高潮久久久久久久久久久不卡| 欧美国产精品一级二级三级| 夫妻午夜视频| 久久久国产成人精品二区 | 国产av精品麻豆| 中文字幕人妻丝袜制服| 一级毛片女人18水好多| 成人18禁在线播放| 亚洲成a人片在线一区二区| a在线观看视频网站| 最新的欧美精品一区二区| 十八禁网站免费在线| 日韩中文字幕欧美一区二区| 黑人操中国人逼视频| 美女扒开内裤让男人捅视频| 黄色片一级片一级黄色片| 黄片大片在线免费观看| 久久精品成人免费网站| 一级片免费观看大全| 男女之事视频高清在线观看| 国产有黄有色有爽视频| 国产亚洲欧美在线一区二区| 亚洲熟妇熟女久久| 久久久水蜜桃国产精品网| 成年人免费黄色播放视频| 国产精品久久久人人做人人爽| 亚洲一码二码三码区别大吗| 校园春色视频在线观看| 激情在线观看视频在线高清 | 久久久久久久久免费视频了| 婷婷成人精品国产| 五月开心婷婷网| 国产成人av教育| 久久 成人 亚洲| 999久久久国产精品视频| 亚洲av日韩精品久久久久久密| 国产精品美女特级片免费视频播放器 | 亚洲精品中文字幕在线视频| 亚洲av成人av| 亚洲人成电影免费在线| 少妇猛男粗大的猛烈进出视频| 好男人电影高清在线观看| 久久久久久久午夜电影 | 国产成人精品久久二区二区免费| 精品人妻1区二区| 久久中文字幕人妻熟女| 精品久久蜜臀av无| 美女国产高潮福利片在线看| 国产蜜桃级精品一区二区三区 | 黑人巨大精品欧美一区二区蜜桃| 久久中文字幕一级| 久久久国产一区二区| 精品福利观看| 亚洲精品久久成人aⅴ小说| 亚洲一区二区三区欧美精品| xxx96com| 午夜精品国产一区二区电影| 老司机影院毛片| 国产区一区二久久| 又黄又爽又免费观看的视频| 80岁老熟妇乱子伦牲交| 制服诱惑二区| 成人精品一区二区免费| 国产一区有黄有色的免费视频| 亚洲一码二码三码区别大吗| 国产人伦9x9x在线观看| 身体一侧抽搐| 亚洲av电影在线进入| 久久久久视频综合| 麻豆成人av在线观看| 99热只有精品国产| 99热国产这里只有精品6| netflix在线观看网站| 天天躁狠狠躁夜夜躁狠狠躁| 天天操日日干夜夜撸| av一本久久久久| 亚洲精品在线美女| 午夜福利在线免费观看网站| 久久久久国产一级毛片高清牌| 色老头精品视频在线观看| 亚洲视频免费观看视频| 亚洲人成电影免费在线| 欧美人与性动交α欧美软件| 精品一品国产午夜福利视频| 如日韩欧美国产精品一区二区三区| 国精品久久久久久国模美| www.自偷自拍.com| 亚洲七黄色美女视频| √禁漫天堂资源中文www| 一区二区三区激情视频| 一个人免费在线观看的高清视频| 成人三级做爰电影| 亚洲中文字幕日韩| 久久精品人人爽人人爽视色| 精品久久蜜臀av无| 日韩中文字幕欧美一区二区| 91成年电影在线观看| 婷婷丁香在线五月| 国产99白浆流出| 国产精品av久久久久免费| 免费av中文字幕在线| 久久婷婷成人综合色麻豆| 99精国产麻豆久久婷婷| 国产片内射在线| 三上悠亚av全集在线观看| 狠狠狠狠99中文字幕| 9色porny在线观看| 涩涩av久久男人的天堂| a级毛片在线看网站| 久久青草综合色| 国产野战对白在线观看| 久久久久国内视频| 一区二区三区国产精品乱码| 久久亚洲精品不卡| 午夜老司机福利片| 午夜激情av网站| 大香蕉久久网| 色综合欧美亚洲国产小说| 黄网站色视频无遮挡免费观看| 90打野战视频偷拍视频| www.精华液| 成人av一区二区三区在线看| 狠狠婷婷综合久久久久久88av| 欧美精品av麻豆av| 亚洲美女黄片视频| 国产精品偷伦视频观看了| netflix在线观看网站| 丰满的人妻完整版| x7x7x7水蜜桃| 麻豆乱淫一区二区| 老司机影院毛片| 日韩三级视频一区二区三区| 成人国语在线视频| 精品一区二区三区视频在线观看免费 | 久久中文看片网| 免费不卡黄色视频| 亚洲中文字幕日韩| 亚洲欧美一区二区三区黑人| 怎么达到女性高潮| 国产又爽黄色视频| 在线观看免费午夜福利视频| 老熟妇乱子伦视频在线观看| 色94色欧美一区二区| 午夜成年电影在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 一边摸一边抽搐一进一小说 | 麻豆av在线久日| 精品高清国产在线一区| 18禁国产床啪视频网站| 亚洲精品乱久久久久久| 成人av一区二区三区在线看| 亚洲五月天丁香| 午夜日韩欧美国产| 国产精品免费一区二区三区在线 | 国产欧美日韩一区二区三区在线| 老鸭窝网址在线观看| 久久国产亚洲av麻豆专区| 在线观看免费视频网站a站| 欧美日韩福利视频一区二区| 国产精品99久久99久久久不卡| 久久精品人人爽人人爽视色| 久久热在线av| 免费少妇av软件| 国产欧美日韩一区二区精品| 777久久人妻少妇嫩草av网站| av有码第一页| 久久ye,这里只有精品| 极品少妇高潮喷水抽搐| 欧美成狂野欧美在线观看| 色婷婷久久久亚洲欧美| 人人妻人人澡人人爽人人夜夜| 国产高清videossex| 99国产极品粉嫩在线观看| 国产黄色免费在线视频| 欧美乱色亚洲激情| 欧美日本中文国产一区发布| 美女扒开内裤让男人捅视频| 涩涩av久久男人的天堂| 久9热在线精品视频| 国产精品国产高清国产av | 欧美国产精品va在线观看不卡| 国产欧美亚洲国产| 中文字幕最新亚洲高清| 丁香欧美五月| 国产国语露脸激情在线看| 女同久久另类99精品国产91| 在线看a的网站| 在线观看免费高清a一片| 亚洲视频免费观看视频| 久久久国产一区二区| 亚洲精品在线观看二区| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品合色在线| 国产一区二区三区在线臀色熟女 | 飞空精品影院首页| 1024视频免费在线观看| 免费日韩欧美在线观看| 黄色视频,在线免费观看| xxxhd国产人妻xxx| 亚洲视频免费观看视频| 国产男女内射视频| 久久国产乱子伦精品免费另类| 99香蕉大伊视频| 亚洲欧美日韩另类电影网站| 狠狠狠狠99中文字幕| 丝袜美足系列| 欧美日本中文国产一区发布| 精品第一国产精品| 久久久精品免费免费高清| 99热网站在线观看| 久久久久久人人人人人| 久久午夜亚洲精品久久| 欧美成人午夜精品| 天天操日日干夜夜撸| 91九色精品人成在线观看| 天堂动漫精品| 亚洲精品国产色婷婷电影| 免费少妇av软件| 日本vs欧美在线观看视频| 身体一侧抽搐| 久久精品亚洲av国产电影网| av免费在线观看网站| 性色av乱码一区二区三区2| 在线观看免费高清a一片| 美女高潮喷水抽搐中文字幕| 国产片内射在线| 美女午夜性视频免费| 桃红色精品国产亚洲av| 亚洲中文日韩欧美视频| avwww免费| 国产精品偷伦视频观看了| 久久精品国产综合久久久| 在线免费观看的www视频| 亚洲avbb在线观看| 亚洲国产毛片av蜜桃av| 热99re8久久精品国产| 久久久久久免费高清国产稀缺| 精品福利永久在线观看| 免费在线观看视频国产中文字幕亚洲| 露出奶头的视频| 人人妻人人添人人爽欧美一区卜| 亚洲成人国产一区在线观看| 国产高清国产精品国产三级| 热99久久久久精品小说推荐| 高清av免费在线| 亚洲国产毛片av蜜桃av| 久久久久久亚洲精品国产蜜桃av| 亚洲 国产 在线| 69精品国产乱码久久久| 免费在线观看影片大全网站| 成人黄色视频免费在线看| 在线观看免费视频日本深夜| 麻豆av在线久日| 手机成人av网站| 久久精品国产亚洲av高清一级| netflix在线观看网站| 女性生殖器流出的白浆| 久久 成人 亚洲| 99在线人妻在线中文字幕 | 91av网站免费观看| 久久精品亚洲精品国产色婷小说| 丝袜美腿诱惑在线| 欧美日韩国产mv在线观看视频| 久久精品人人爽人人爽视色| 在线国产一区二区在线| 午夜福利在线免费观看网站| 亚洲午夜理论影院| 成在线人永久免费视频| 色播在线永久视频| 亚洲av电影在线进入| 99国产精品99久久久久| 日韩熟女老妇一区二区性免费视频| 日日夜夜操网爽| 99在线人妻在线中文字幕 | 久久人妻熟女aⅴ| 欧美乱码精品一区二区三区| 亚洲性夜色夜夜综合| 国产1区2区3区精品| 亚洲精品美女久久久久99蜜臀| 制服人妻中文乱码| 国产一卡二卡三卡精品| 欧美精品av麻豆av| 黄片小视频在线播放| 国产aⅴ精品一区二区三区波| 久久久精品免费免费高清| 法律面前人人平等表现在哪些方面| 一级a爱片免费观看的视频| 国产亚洲欧美98| 99国产极品粉嫩在线观看| 老司机午夜十八禁免费视频| 91成人精品电影| 极品人妻少妇av视频| 国产在线一区二区三区精| 国产极品粉嫩免费观看在线| 日韩制服丝袜自拍偷拍| 操美女的视频在线观看| 日韩成人在线观看一区二区三区| 777米奇影视久久| 在线视频色国产色| 欧美一级毛片孕妇| 欧美精品高潮呻吟av久久| 18禁裸乳无遮挡动漫免费视频| 免费黄频网站在线观看国产| 精品久久蜜臀av无| 欧美黑人精品巨大| 国产精品 欧美亚洲| 欧美在线黄色| 国产成人啪精品午夜网站| 成年动漫av网址| 色尼玛亚洲综合影院| 久久亚洲精品不卡| 国产精品一区二区在线不卡| 99热只有精品国产| 国精品久久久久久国模美| av有码第一页| 欧美乱码精品一区二区三区| 亚洲欧美一区二区三区黑人| 97人妻天天添夜夜摸| 午夜福利在线观看吧| 一边摸一边做爽爽视频免费| 男人操女人黄网站| 男人舔女人的私密视频| 十分钟在线观看高清视频www| 在线视频色国产色| 欧美中文综合在线视频| 亚洲,欧美精品.| 精品第一国产精品| 国产三级黄色录像| 91在线观看av| 久久99一区二区三区| 一边摸一边抽搐一进一出视频| 亚洲五月天丁香| 中文亚洲av片在线观看爽 | 91成人精品电影| 久久久精品国产亚洲av高清涩受| 极品少妇高潮喷水抽搐| av网站在线播放免费| 午夜福利乱码中文字幕| 欧美大码av| 日韩欧美三级三区| 看免费av毛片| 国产av精品麻豆| 久久精品亚洲熟妇少妇任你| 一级,二级,三级黄色视频| 亚洲欧美日韩另类电影网站| 一区在线观看完整版| 伊人久久大香线蕉亚洲五| 亚洲成a人片在线一区二区| 国产一区有黄有色的免费视频| 熟女少妇亚洲综合色aaa.| 一进一出好大好爽视频| 国产99久久九九免费精品| 国产亚洲欧美在线一区二区| 国产免费男女视频| 人成视频在线观看免费观看| 超碰成人久久| 国产有黄有色有爽视频| 又黄又爽又免费观看的视频| 日韩大码丰满熟妇| 色老头精品视频在线观看| 极品教师在线免费播放| 亚洲一区中文字幕在线| 日韩大码丰满熟妇| 亚洲全国av大片| 在线av久久热| 久久久精品免费免费高清| 国产精品一区二区免费欧美| 美女高潮到喷水免费观看| 在线观看舔阴道视频| 狠狠狠狠99中文字幕| 久热这里只有精品99| 欧美一级毛片孕妇| 亚洲第一青青草原| 在线国产一区二区在线| 69av精品久久久久久| 又紧又爽又黄一区二区| 欧美人与性动交α欧美精品济南到| 久久精品国产a三级三级三级| 大码成人一级视频| 他把我摸到了高潮在线观看| 日日摸夜夜添夜夜添小说| 久久中文字幕一级| 侵犯人妻中文字幕一二三四区| 两性夫妻黄色片| 老熟女久久久| 亚洲综合色网址| 他把我摸到了高潮在线观看| 看片在线看免费视频| 欧美日韩亚洲综合一区二区三区_| 精品福利观看| 一级作爱视频免费观看| 日韩成人在线观看一区二区三区| 久久ye,这里只有精品| 久久香蕉国产精品| 欧美日韩瑟瑟在线播放| 久久九九热精品免费| 伊人久久大香线蕉亚洲五| 久久久久国产一级毛片高清牌| 一进一出好大好爽视频| 一区二区三区国产精品乱码| 国产成人精品无人区| 麻豆乱淫一区二区| x7x7x7水蜜桃| 人妻 亚洲 视频| 国产99久久九九免费精品| 精品国产乱码久久久久久男人| 欧美精品亚洲一区二区| 欧美日韩福利视频一区二区| 国产精品久久久av美女十八| 欧美精品av麻豆av| 露出奶头的视频| 老熟妇乱子伦视频在线观看| 午夜精品久久久久久毛片777| 涩涩av久久男人的天堂| 中文字幕最新亚洲高清| 露出奶头的视频| 精品福利观看| 国产极品粉嫩免费观看在线| 国产精品香港三级国产av潘金莲| 亚洲欧美日韩高清在线视频| 日韩欧美免费精品| 日本黄色视频三级网站网址 | 两个人免费观看高清视频| 人人妻人人添人人爽欧美一区卜| 午夜亚洲福利在线播放| 久久久精品区二区三区| 久久久久久亚洲精品国产蜜桃av| 18禁国产床啪视频网站| 久久精品国产综合久久久| 亚洲五月天丁香| 18禁观看日本| 熟女少妇亚洲综合色aaa.| 久久久久精品人妻al黑| 免费日韩欧美在线观看| 成人黄色视频免费在线看| tube8黄色片| 两个人看的免费小视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美成人午夜精品| 女人精品久久久久毛片| 精品乱码久久久久久99久播| 黑人操中国人逼视频| 免费看十八禁软件| 国产深夜福利视频在线观看| svipshipincom国产片| 在线天堂中文资源库| 欧美日韩瑟瑟在线播放| 在线观看日韩欧美| 国产真人三级小视频在线观看| 夜夜躁狠狠躁天天躁| 黄频高清免费视频| 国产无遮挡羞羞视频在线观看| 人人妻人人澡人人看| 精品视频人人做人人爽| 看片在线看免费视频| 亚洲av欧美aⅴ国产| 精品电影一区二区在线| 国产精品欧美亚洲77777| 99精品欧美一区二区三区四区| 精品国产乱子伦一区二区三区| 啦啦啦 在线观看视频| 欧美激情久久久久久爽电影 | 80岁老熟妇乱子伦牲交| 午夜日韩欧美国产| a在线观看视频网站| 午夜影院日韩av| 国产aⅴ精品一区二区三区波| 99国产综合亚洲精品| 国产亚洲一区二区精品| 亚洲免费av在线视频| 黄片小视频在线播放| 欧美乱色亚洲激情| 国产欧美日韩一区二区精品| 国产精品香港三级国产av潘金莲| 91精品国产国语对白视频| 最近最新免费中文字幕在线| 久久精品国产a三级三级三级| 国产人伦9x9x在线观看| 91av网站免费观看| 亚洲国产欧美日韩在线播放| 午夜久久久在线观看| 成人亚洲精品一区在线观看| 国产精品亚洲一级av第二区| 亚洲免费av在线视频| 男女床上黄色一级片免费看| 99国产精品一区二区三区| 欧美人与性动交α欧美精品济南到| 精品视频人人做人人爽| 亚洲中文日韩欧美视频| 超色免费av| 丁香欧美五月| 久久 成人 亚洲| 1024视频免费在线观看| 午夜福利视频在线观看免费| 精品国产乱子伦一区二区三区| 人成视频在线观看免费观看| 久久精品国产a三级三级三级| 久久午夜综合久久蜜桃| 日本a在线网址| 久久久久久久国产电影| 国产高清videossex| 香蕉丝袜av| 老汉色∧v一级毛片| 午夜久久久在线观看| 国产男女超爽视频在线观看| 热re99久久精品国产66热6| 人妻久久中文字幕网| 五月开心婷婷网| 女同久久另类99精品国产91| 狂野欧美激情性xxxx| 国产三级黄色录像| 午夜两性在线视频| 亚洲第一欧美日韩一区二区三区| 天堂√8在线中文| 国产精品99久久99久久久不卡| 99热网站在线观看| 老汉色av国产亚洲站长工具| 欧美日韩精品网址| aaaaa片日本免费| 美女高潮喷水抽搐中文字幕| 一边摸一边抽搐一进一出视频| 很黄的视频免费| 99国产精品免费福利视频| 97人妻天天添夜夜摸| 国产97色在线日韩免费| 国产午夜精品久久久久久| 一级a爱视频在线免费观看| 国产精品亚洲av一区麻豆| 捣出白浆h1v1| 一级a爱视频在线免费观看| 日日夜夜操网爽| 夜夜爽天天搞| 男女午夜视频在线观看| av国产精品久久久久影院| 麻豆av在线久日| 在线观看一区二区三区激情| 91精品国产国语对白视频| 国产精华一区二区三区| 亚洲片人在线观看|