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

    目標(biāo)角度估計的多輸入多輸出雷達發(fā)射方向圖綜合

    2015-12-28 00:59:35黃中瑞,鄭志東,張劍云
    電波科學(xué)學(xué)報 2015年4期

    目標(biāo)角度估計的多輸入多輸出雷達發(fā)射方向圖綜合

    黃中瑞1鄭志東2張劍云1

    (1.電子工程學(xué)院,安徽 合肥 230037;2.北方電子設(shè)備研究所,北京 100083)

    摘要針對傳統(tǒng)多輸入多輸出(Multiple-Input Multiple-Output,MIMO)雷達發(fā)射功率分散導(dǎo)致目標(biāo)角度估計性能變差的問題,提出了一種新的MIMO雷達發(fā)射方向圖綜合方法.首先,結(jié)合發(fā)射功率的空域聚焦范圍和發(fā)射導(dǎo)向矢量的旋轉(zhuǎn)不變性,嚴(yán)格控制方向圖逼近誤差小于給定的門限,并且最小化發(fā)射方向圖的峰值旁瓣功率;其次,約束協(xié)方差矩陣的對角元素相等,以實現(xiàn)發(fā)射功率效率最大化,在此基礎(chǔ)之上推導(dǎo)了優(yōu)化模型的二階錐規(guī)劃形式,以便采用原對偶內(nèi)點算法進行有效求解;最后,利用平行因子分析算法對空間目標(biāo)進行了收發(fā)角度估計.仿真結(jié)果表明:優(yōu)化所得加權(quán)矩陣不僅能夠使發(fā)射信號在指定空域進行有效聚焦,而且確保了每個陣元的發(fā)射功率相等,因而在發(fā)射功率和角度估計方法相同時,所提方法相對傳統(tǒng)MIMO雷達具有更好的角度估計精度.

    關(guān)鍵詞MIMO雷達;二階錐規(guī)劃;波束聚焦;方向圖綜合;角度估計

    中圖分類號TN958

    文獻標(biāo)志碼A

    文章編號1005-0388(2015)04-0789-08

    AbstractIn order to solve the problem of worse angel estimation which caused by the power being transmitted dispersedly in traditional multiple-input multiple-output(MIMO) radar, a new method for pattern synthesis of MIMO radar is proposed. Firstly, the optimal model is constructed based on second-order cone programming (SOCP) by combining the focus range of the transmit power with the rotational invariance of the transmit steering vector. It can not only constrain the optimal error less than the given value, but also can minimize the peak side lobe power of the transmit pattern. Further, the diagonal elements of the covariance matrix are constrained, which can maximum the power efficiency. Moreover, the SOCP form of the optimal model is given so as to use the primal-dual interior point method effectively. Lastly, simulations illustrate that the

    收稿日期:2014-09-04

    作者簡介

    Transmit pattern synthesis of MIMO radar for

    the angle estimation

    HUANG Zhongrui1ZHENG Zhidong2ZHANG Jianyun1

    (1.ElectronicEngineeringInstitution,AnhuiHefei230037,China;

    2.InstituteofNorthElectronicEquipment,Beijing100083,China)

    資助項目: 國家自然科學(xué)基金(61201279); 安徽省自然科學(xué)基金(1408085MF128)

    聯(lián)系人: 黃中瑞 E-mail: 18756073857@163.com

    beam weight matrix can make the transmit power focus within the desire space and constrain the transmit array power to be equal, which is propitious to the energy usage. Therefore, the angle estimation performance by using the proposed method in this paper is better than the traditional MIMO radar when the transmit power and angle estimation method are same.

    Key words MIMO radar; second-order cone programming; beam focus; pattern synthesis; angle estimation

    引言

    多輸入多輸出(Multiple-Input Multiple-Output,MIMO)雷達是近年來提出的一種新型體制雷達[1].由于MIMO雷達的每個陣元可以獨立發(fā)射波形因而具有更高的自由度,與相控陣相比,同等陣元配置條件下可以辨識更多的目標(biāo)[2-3].MIMO雷達按其天線陣元的配置方式可分為:統(tǒng)計MIMO雷達[4]和相干MIMO雷達[5-6].統(tǒng)計MIMO雷達在空間采用大間距配置方式,它能夠從不同的角度觀測目標(biāo),利用系統(tǒng)的空間分集增益克服目標(biāo)的爍效應(yīng)[7],提高閃爍目標(biāo)的檢測性能;相干MIMO雷達的陣元配置間距較小,它們相對于遠(yuǎn)場目標(biāo)而言處于同一視角下,因而無法獲得空間分集增益,但是相干MIMO雷達可以利用波形分集形成大的虛擬孔徑,提高雷達的角度估計精度,增加最大可識別的目標(biāo)數(shù)[8]等.

    目前,多數(shù)文獻[9-11]研究了發(fā)射信號為理想正交時MIMO雷達的角度估計問題,這些文獻均假設(shè)MIMO雷達信號為全向等功率發(fā)射,其在全空域滿足均勻分布.但是如果目標(biāo)只集中于某個很小的空域范圍內(nèi),由于發(fā)射功率的分散使得目標(biāo)接收到的發(fā)射功率很小,從而降低了每個接收陣元的信噪比,不利于對目標(biāo)收發(fā)角度的估計.

    針對這個問題,本文對相干(雙基地)MIMO雷達基于接收角度(Direction of Arrival,DOA)和發(fā)射角度(Direction of Departure,DOD)估計的發(fā)射方向圖綜合問題展開研究.首先要考慮的就是MIMO雷達發(fā)射方向圖的問題,而方向圖綜合問題可以歸結(jié)為發(fā)射波形的設(shè)計問題,主要分為兩類:正交信號和部分相關(guān)信號.正交信號[12]的優(yōu)點在于接收端通過匹配濾波器能夠方便提出不同傳播路徑的信息,較好地利用了虛擬空間所帶來的自由度,因而具有很好的參數(shù)估計性能和雜波抑制能力.但是MIMO雷達各陣元發(fā)射的正交波形在空中不能進行相干疊加,因而其發(fā)射信號為全向等功率輻射.在目標(biāo)分布區(qū)域已知的情況下,正交波形的發(fā)射功率分散問題會導(dǎo)致回波信噪比嚴(yán)重下降,既減小了MIMO雷達的作用距離,也在一定程度上限制了參數(shù)估計精度.與此不同的是,對MIMO雷達發(fā)射信號的相關(guān)性進行適當(dāng)設(shè)計,即可實現(xiàn)發(fā)射功率的聚集從而有效提高接收端信號處理的性能[13-14].雖然部分相關(guān)信號能夠克服發(fā)射功率發(fā)散的問題,但是其也存在兩大難點:第一,基于部分相關(guān)信號設(shè)計的發(fā)射方向圖綜合問題,實質(zhì)上是優(yōu)化發(fā)射信號的相關(guān)矩陣,而從優(yōu)化得到的相關(guān)矩陣到具體所需的發(fā)射波形是一類非常復(fù)雜的問題;第二,如何從具有相關(guān)性的回波中提取獨立的目標(biāo)信息仍是一個亟待解決的問題.

    綜上所述,正交信號和部分相關(guān)信號是MIMO雷達發(fā)射信號設(shè)計的兩個對立面,這一矛盾的產(chǎn)生主要是由于大部分文獻方法無法將方向圖設(shè)計與波形設(shè)計進行分離,使得兩者相互制約.文獻[15-16]提出了一種獨立于波形設(shè)計的發(fā)射方向圖綜合方法.其主要思想是:采用一組正交波形集作為基信號,而每個陣元的發(fā)射信號為各基信號的線性組合.這樣對MIMO雷達的方向圖綜合問題可以轉(zhuǎn)化為各個基信號線性組合系數(shù)的優(yōu)化問題,進而實現(xiàn)了方向圖設(shè)計問題與發(fā)射波形設(shè)計問題的有效的分離.但是,文獻[15-16]仍存在一定問題,比如文獻[15]在對發(fā)射方向圖進行優(yōu)化設(shè)計時,只能對基信號線性組合系數(shù)的相關(guān)矩陣進行優(yōu)化,無法直接對線性組合系數(shù)進行優(yōu)化并且要求其相關(guān)矩陣具有滿秩性;文獻[16]無法對MIMO雷達每個陣元的發(fā)射功率進行有效約束,使得不同陣元間的發(fā)射功率差異過大.

    鑒于此,本文提出了一種基于DOA和DOD估計的MIMO雷達發(fā)射方向圖綜合方法.在聚焦空間內(nèi)嚴(yán)格約束優(yōu)化波束與期望波束的誤差小于預(yù)設(shè)門限的條件下,利用二階錐規(guī)劃(Second-Order Cone Programming,SOCP)方法對基波束的加權(quán)矩陣進行設(shè)計,從而最小化方向圖的旁瓣功率,同時為了確保功率的充分利用,對MIMO雷達的每個發(fā)射陣元的最大輻射功率進行了一定的限制.該方法不僅能夠?qū)訖?quán)矩陣進行直接求解,而且能夠使新的發(fā)射方向矢量具有旋轉(zhuǎn)不變性,對后文采用平行因子算法對目標(biāo)進行收發(fā)角度估計作了很好的鋪墊.

    需要說明的是,為了表述方便,文中將MIMO雷達每個發(fā)射陣元發(fā)射理想正交信號的模式記為傳統(tǒng)MIMO雷達.

    1發(fā)射信號模型

    考慮一雙基地MIMO雷達,其發(fā)射陣列和接收陣列分別為均勻線陣,陣元數(shù)分別為M和N,陣元間距分別為dt和dr,λ表示載波波長.并假設(shè)遠(yuǎn)場相同距離分辨單元內(nèi)存在P個目標(biāo),相對于發(fā)射陣和接收陣的角度分別為(φp,θp),其中p=1,…,P,陣列配置如圖1所示.

    圖1 雙基地MIMO雷達收發(fā)陣元配置圖

    令Φ(t)=[φ1(t)φ2(t)…φK(t)]T為K個發(fā)射基信號的復(fù)包絡(luò),基信號之間為歸一化正交,即滿足

    (1)

    式中:Tp為基帶波形的脈沖寬度;(·)T、(·)H和(·)*分別表示向量(矩陣)的轉(zhuǎn)置、共軛轉(zhuǎn)置和共軛.

    為了將發(fā)射信號的功率在指定的空域進行聚焦并且獲得方向圖設(shè)計和波形設(shè)計的分離,需要對基信號進行加權(quán),從而獲得具有一定相關(guān)性的發(fā)射信號.值得注意的是,這里的發(fā)射信號在接收端經(jīng)過匹配濾波可以很好地提取空間自由度.令S(t)=[s1(t)s2(t)…sM(t)]T表示MIMO雷達發(fā)射陣元的實際發(fā)射信號矢量,W為加權(quán)矩陣,因而有

    S(t)=WΦ(t),

    (2)

    W為M×K維的矩陣.

    發(fā)射信號為窄帶信號且假設(shè)不考慮各種衰減因素,則在遠(yuǎn)場φ處的輻射信號為

    y(t,φ)=aH(φ)S(t)=aH(φ)(WΦ(t)).

    (3)

    式中,a(φ)為發(fā)射陣列的方向?qū)蚴噶?且有a(φ)=[1 ej2πdsin φ/λ,…,ej2π(M-1)dsin φ/λ]T.根據(jù)式(3)可以得到發(fā)射信號在空間的功率分布

    P(φ) =κE{y(t,φ)yH(t,φ)}

    =κE{aH(φ)(WΦ(t))(aH(φ)(WΦ(t)))H}

    =κaH(φ)E{WΦ(t)ΦH(t)WH}a(φ)

    =κaH(φ)WWHa(φ).

    (4)

    式中,κ為常數(shù)因子,對后面的優(yōu)化問題沒有影響,因而在后續(xù)問題的考慮中將其略去.

    由式(4)可知,MIMO雷達發(fā)射功率的分配問題實質(zhì)上就轉(zhuǎn)化為如何對加權(quán)矩陣W進行設(shè)計.下面基于SOCP給出了W的設(shè)計方法,以便使得發(fā)射信號功率在指定空域內(nèi)進行聚焦,在不感興趣的空域中能量最?。?/p>

    2基于SOCP的加權(quán)矩陣設(shè)計

    由式(3)可知,在接收端用基信號對接收陣元的接收信號進行匹配濾波后,得到的發(fā)射導(dǎo)向矢量將被W加權(quán),從而使得發(fā)射方向?qū)蚴噶縜(φ)不再具有平移不變特性,將不利于接收端的參數(shù)估計,為了克服這個局限性,需要對新的發(fā)射導(dǎo)向矢量WHa(φ)進行約束使其仍然具有旋轉(zhuǎn)不變性.

    首先對空域進行分隔,令ΞZ表示波束需要聚焦的空域,ΞS表示不感興趣的空域,為了方便后續(xù)優(yōu)化模型的建立和仿真的需要,這里對ΞZ和ΞS進行離散化處理,離散點數(shù)分別為Z和S即:φz∈ΞZ,z=1,…,Z,φs∈ΞS,s=1,…,S.

    ‖WHa(φz)-d(φz)‖≤σ,

    φz∈ΞZ,z=1,…,Z.

    (5)

    式中,σ≥0為預(yù)設(shè)門限值,通過對門限值進行合理設(shè)置可以使新的發(fā)射方向?qū)蚴噶拷茲M足旋轉(zhuǎn)不變性.

    對于不感興趣的空域,需要約束其輻射功率最小,這樣可以使得發(fā)射功率更好地在感興趣空域進行聚焦,因而結(jié)合式(4)可得優(yōu)化目標(biāo)函數(shù)為

    (6)

    結(jié)合式(5)和(6)便可得到加權(quán)矩陣W的優(yōu)化準(zhǔn)則,指定空域內(nèi)優(yōu)化波束與期望波束誤差小于預(yù)設(shè)門限值的條件下,最小化發(fā)射信號在不感興趣空域內(nèi)的功率,即

    s.t.‖WHa(φz)-d(φz)‖≤σ,

    φz∈ΞZ,z=1,…,Z.

    (7)

    上述優(yōu)化模型已經(jīng)完成了優(yōu)化波束對期望波束在形狀上的理想逼近,但遺憾的是優(yōu)化后MIMO雷達每個陣元的實際發(fā)射功率相差甚大,這不但有礙于能量的充分利用,而且對整個MIMO雷達系統(tǒng)性能的發(fā)揮也起著不利作用.

    (8)

    聯(lián)合式(7)和(8)便可得到本文所提算法的優(yōu)化模型:

    s.t. ‖WHa(φz)-d(φz)‖≤σ,

    φz∈ΞZ,z=1,…,Z;

    (9)

    式中τ為預(yù)設(shè)值,作為單個陣元發(fā)射功率的門限,由凸優(yōu)化理論可知,式(9)為典型的SOCP問題,為了更好地采用凸優(yōu)化(Convex,CVX)或SeDuMi工具箱中的內(nèi)點算法進行求解,下文基于凸優(yōu)化理論給出其SOCP形式.

    給出與式(9)等價的優(yōu)化模型為

    s.t. ‖WHa(φs)‖≤ξS

    φs∈ΞS,s=1,…,S,

    ‖WHa(φz)-d(φz)‖≤σ,

    φz∈ΞZ,z=1,…,Z;

    ‖wi‖≤τ,i=1,…,M.

    (10)

    式中,wi表示加權(quán)矩陣W的第i行,對WH按行拉直可以表示為ν=vec(W*)=νr+jνi,νr和νi分別表示向量ν的實部和虛部,其維數(shù)為KM×1.

    minυTy

    φs∈ΞS,s=1,…,S;

    φz∈ΞZ,z=1,…,Z;

    i=1,…,M.

    (11)

    從式(11)可知,優(yōu)化模型(10)共包含了M+S+Z個錐約束,利用內(nèi)點算法可以方便地求得優(yōu)化變量y,進而求得加權(quán)矩陣W.

    3接收端匹配處理和DOD估計

    在上節(jié)中利用SOCP對基信號的加權(quán)矩陣進行了優(yōu)化,本節(jié)在此基礎(chǔ)之上,對接收信號的匹配濾波處理和DOA估計作簡要分析.

    3.1接收端信號匹配處理

    基于波束空間聚焦的MIMO雷達發(fā)射信號是其基信號的線性組合,所以發(fā)射信號之間存在一定的相關(guān)性.如果直接采用傳統(tǒng)MIMO雷達的發(fā)射信號作為接收端匹配信號進行濾波,將無法有效提取空間自由度.為了解決這個問題,必須利用基信號的共軛信號Φ(t)對接收信號進行匹配濾波,由文獻[16]可知接收信號為

    (12)

    Ω=[D(φ)⊙B(θ)]ΗT+.

    (13)

    3.2DOD估計

    當(dāng)離線得到基信號加權(quán)矩陣和接收端的數(shù)據(jù)之后,便可采用多重信號分類(Multiple Signal Classification,MUSIC)[17]、子空間旋轉(zhuǎn)不變技術(shù)(Estimating Signal Parameters via Rotational Invariance Techniques,ESPRIT)[18]和平行因子分析算法[16]對目標(biāo)的收發(fā)角度進行估計.這些算法在現(xiàn)有的文獻中已經(jīng)被研究的比較成熟,這里就不再贅述.

    需要說明的是,盡管在優(yōu)化加權(quán)矩陣時已經(jīng)限定了發(fā)射導(dǎo)向矢量的旋轉(zhuǎn)不變性,但在直接利用已有的文獻方法對收發(fā)射角度進行估計時,接收角雖然能夠準(zhǔn)確地得到,但是發(fā)射角的估計仍會存在一定的問題,這是因為真正的發(fā)射導(dǎo)向矢量已經(jīng)被加權(quán).為此,主要闡述采用查表法解決MIMO雷達空間目標(biāo)發(fā)射角度的估計問題.

    在第2節(jié)得到加權(quán)矩陣之后,對于空間某一確定方向φ,可以寫出實際發(fā)射波束導(dǎo)向矢量的各個元素為

    (14)

    式中,Δk(φ)(k=1,…,K)為逼近誤差.對于式(14)的相鄰兩項作除法,即可得到加權(quán)后導(dǎo)向矢量間的相位差,并且這個相位差與空間角度相對應(yīng),將其進行存儲.在接收端采用同樣方式對加權(quán)導(dǎo)向矢量進行運算,將計算的相位值與存儲相位值進行匹配,其中差值最小相位對應(yīng)的角度即可作為空間目標(biāo)相對MIMO雷達的發(fā)射角估計值.

    4仿真實驗

    為了驗證所提方法的有效性,給出了一些仿真結(jié)果.假設(shè)雙基地MIMO雷達的發(fā)射陣元和接收陣元數(shù)目均為10,陣元間的距離均為半波長,發(fā)射接收陣元各向同性輻射,陣元的總輻射功率E=10.

    4.1MIMO雷達發(fā)射方向圖性能比較

    假設(shè)基信號的個數(shù)K=3,指定的目標(biāo)空域ΞZ=[-10°,10°],期望波束d(φ)=[0 e-j4πsin φ,…,e-j8πsin φ]T,方向圖的旁瓣區(qū)為ΞS=[-90°,-20°]∪[20°,90°],約束期望波束的絕對和誤差為0.36,單個陣元的發(fā)射功率上限為0.072(歸一化值),指定目標(biāo)空域和旁瓣區(qū)的離散化點數(shù)分別為201和202,圖2給出了基于SOCP優(yōu)化所得空間波束聚焦方向圖和傳統(tǒng)MIMO雷達的發(fā)射方向圖.

    圖2 基于SOCP聚焦的發(fā)射方向圖和 傳統(tǒng)MIMO雷達全向發(fā)射方向圖比較

    從圖2可知,基于SOCP聚焦的方向圖在指定空域具有很高的增益,在散射系數(shù)和噪聲一定的條件下,能夠有效地提高接收端信號的信噪比,從而有利于目標(biāo)的檢測和參數(shù)估計.圖3給出了本文優(yōu)化所得MIMO雷達發(fā)射陣元的能量分布圖.

    從圖3可以看出,基于本文算法優(yōu)化所得的加權(quán)矩陣可以使MIMO雷達每個發(fā)射陣元的發(fā)射功率近似相同,因而更有益于發(fā)射功率的利用,優(yōu)于文獻[16]中單個陣元的功率分布(圖4).

    4.2算法的有效性驗證

    假設(shè)空間存在三個不相干的目標(biāo),目標(biāo)相對發(fā)射和接收陣列的角度為(-2°,-5°),(0°,0°),(2°,5°),并利用平行因子分析法對目標(biāo)角度進行估計,信噪比為0dB,蒙特卡洛實驗次數(shù)為10,本文算法、文獻[16]算法和傳統(tǒng)MIMO雷達估計角度的星座圖分別為圖5、圖6和圖7所示.

    圖3 本文方法的發(fā)射陣元功率分布圖

    圖4 文獻[16]單個陣元發(fā)射功率分布圖

    圖5 基于本文算法的目標(biāo)角度估計星座圖

    圖5~7可以看出,本文算法和文獻[16]算法由于實現(xiàn)了發(fā)射方向圖在指定空域的聚焦,因而在相同發(fā)射功率的情況下,多次實驗?zāi)繕?biāo)角度的估計值更加集中并且與真實值之間的誤差更小,優(yōu)于傳統(tǒng)MIMO雷達全向等功率發(fā)射時的估計結(jié)果,其中本文算法的聚集性最好.

    為了進一步驗證所提算法的有效性,分別采用本文算法、文獻[16]算法和傳統(tǒng)MIMO雷達對空中兩個近目標(biāo)的角度估計精度和分辨率作仿真分析.假設(shè)空中存在兩個不相干的目標(biāo),目標(biāo)相對發(fā)射和接收陣列的角度為(-1°,-2°),(1°,2°),驗證兩種情況下目標(biāo)的估計均方誤差和成功率隨信噪比的變化情況,設(shè)置信噪比為[-40,20],其離散化點數(shù)為16,蒙特卡洛實驗次數(shù)為50.結(jié)果如圖8和9所示.

    圖6 基于文獻[16]的目標(biāo)角度估計星座圖

    圖7 基于傳統(tǒng)MIMO雷達的目標(biāo)角度估計星座圖

    圖8 目標(biāo)發(fā)射角度成功分辨概率

    由圖8和9可知,本文和文獻[16]所提方法在信噪比為-15dB左右時,就能百分之百地分辨出兩個目標(biāo),而傳統(tǒng)MIMO雷達只有在信噪比為-5dB時分辨成功率才能達到100%.

    圖9 目標(biāo)接收角度成功分辨概率

    圖10和11給出了空間兩個目標(biāo)收發(fā)角度估計均方誤差隨信噪比的變化曲線.由圖可知,在信噪比相同時,本文算法和文獻[16]方法相對傳統(tǒng)MIMO雷達具有更好的角度估計精度,其中本文算法估計出的角度具有最低的均方誤差,這說明在角度估計精度一定時,本文所需的信噪比最低,這對目標(biāo)參數(shù)估計是十分重要.

    圖10 目標(biāo)1的角度估計均方誤差

    圖11 目標(biāo)2的角度估計均方誤差

    由圖8至圖11可以看出:在信噪比低于-20dB時,文獻[16]的性能略好于本文算法;而在信噪比高于-20dB時,則本文算法的性能要好于文獻[16].其主要原因是,文獻[16]主要采用的是嚴(yán)格約束方向圖的旁瓣電平小于某個門限值,從而能夠獲得更大的主瓣增益,因而在低信噪比條件下,有效地提高了接收信號的能量,獲得更好的參數(shù)估計能力.但是由于文獻[16]只對旁瓣電平作了硬性約束,在提高主瓣增益的同時也放寬了對期望波束的逼近誤差,這也是高信噪比條件下,本文算法具有較好參數(shù)估計性能的原因.更重要的是,本文在兼顧方向圖性能的同時,能夠很好地約束陣元的發(fā)射功率更加趨于均勻化.

    5結(jié)論

    為了解決傳統(tǒng)MIMO雷達發(fā)射功率的發(fā)散問題,本文基于DOA和DOD估計提出了一種MIMO雷達的發(fā)射方向圖綜合方法,新方法能夠有效地將發(fā)射功率向感興趣的空間進行聚焦.該方法將空間波束聚焦范圍和實際發(fā)射導(dǎo)向矢量的旋轉(zhuǎn)不變性相結(jié)合,約束優(yōu)化波束和期望波束的誤差嚴(yán)格小于預(yù)設(shè)門限值.同時為了保證發(fā)射功率的充分利用,對每個陣元的實際輻射功率進行限定的前提下,最小化方向圖的旁瓣電平.仿真實驗表明:本文所提方法能夠保證MIMO雷達的發(fā)射功率在指定空域內(nèi)進行有效聚焦,并且相對文獻[16],陣元的發(fā)射功率更趨均勻化.緊接著利用平行因子分析法對目標(biāo)的角度進行估計,由于本文算法提高了接收端的信噪比,因而在估計算法和發(fā)射功率相同的條件下,對目標(biāo)角度的估計性能優(yōu)于傳統(tǒng)MIMO雷達.

    參考文獻

    [1]FISHERE,HAIMOVICHA,BLUMR,etal.MIMOradar:anideawhosetimehascome[C]//ProcIntRasarConf.Philadelphia,April26-29, 2004: 71-78.

    [2]FISHERE,HAIMOVICHA,BLUMR,etal.Spatialdiversityinradar-modelsanddetectionperformance[J].IEEETransactionsonSignalProcessing, 2006, 54(3): 823-838.

    [3]LUXZ,LIJ,STOCIAP.TargetdetectionandparameterestimationforMIMOradarsystems[J].IEEETransactionsonAerospaceandElectronicSystems, 2008, 43(3):927-939.

    [4]HAIMOVICH,BLUMRS,CIMINIL.MIMOradarwithwidelyseparatedantennas[J].IEEESignalProcessingMagazine, 2008, 25(1): 116-129.

    [5]YANGML,CHENBX,YANGXY.ConjugateESPRITalgorithmforbistaticMIMOradar[J].IEEElectronicsLetters, 2010, 46 (25): 1692-1694.

    [6]STOCIAP,LIJ.MIMOradarwithco-locatedantenna:reviewofsomerecentwork[J].IEEESignalProcessingMagazine, 2007, 24: 106-114.

    [7]HEQ,BLUMRS,GODRICHH,etal.TargetvelocityestimationandantennaplacementforMIMOradarwithwidelyseparatedantennas[J].IEEEJournalofSelectedTopicsinSignalProcessing, 2010, 4(1): 79-100.

    [8]ZHANGX,XUD.AngleestimationinMIMOradarusingreduced-dimensionCapon[J].IEEElectronicsLetters, 2010, 46(12): 860-861.

    [9]劉紅明, 何子述, 李軍. 雙基地MIMO雷達發(fā)射陣目標(biāo)角度估計[J]. 電波科學(xué)學(xué)報, 2010, 25(3): 499-504.

    LIUHongming,HEZishu,LIJun.TargetangleestimationrelativetotransmittingarrayinbistaticMIMOradar[J].ChineseJournalofRadioScience, 2010, 25(3): 499-504.(inChinese)

    [10]李樹鋒, 張履謙, 陳杰, 等. 基于完全互補序列的正交MIMO雷達二維DOA估計[J]. 電波科學(xué)學(xué)報, 2010, 25(4): 617-623.

    LIShufeng,ZHANGLüqian,CHENJie,etal.Two-dimensionalDOAestimationfororthogonalMIMOradarbasedoncompletecomplementarysqrence[J].ChineseJournalofRadioScience, 2010, 25(4): 617-623. (inChinese)

    [11]陳顯舟, 楊源, 韓靜靜. 雙基地多入多出雷達收發(fā)方位角聯(lián)合估計算法[J]. 電波科學(xué)學(xué)報, 2013, 28(1): 176-182.

    CHENXianzhou,YANGYuan,HANJingjing,etal.JointDODandDOAestimationusingpolynomialrootingforbistaticMIMOradar[J].ChineseJournalofRadioScience, 2013, 28(1): 176-182. (inChinese)

    [12]DENGH.Discretefrequency-codingwaveformdesignfornettedradarsystems[J].IEEESignalProcessingLetters, 2004, 11(2): 179-182.

    [13]HASSANIENA,VOROBYOVSA.TransmitenergyfocusingforDOAestimationinMIMOradarwithcolocatedantennas[J].IEEETransactionsonSignalProcessing, 2011, 59(6): 2669-2682.

    [14]AHMEDS,THOMPSONJS,MULGREWB,etal.UnconstrainedsynthesisofcovariancematrixforMIMOradartransmitbeampattern[J].IEEETransactionsonSignalProcessing, 2011, 59(8): 3837-3849.

    [15]BERRYPE,YAUD.Optimalfast-timebeamformingwithlinearlyindependentwaveforms[J].SignalProcessing, 2009, 89(4): 492-495.

    [16]鄭志東, 張劍云, 楊瑛. 基于發(fā)射波束域-平行因子分析的MIMO雷達收發(fā)角度估計[J]. 電子與信息學(xué)報, 2011, 33(12): 2875-2880.

    ZHENGZhidong,ZHANGJianyun,YANGYing.JointDOD-DOAestimationofMIMOradarbasedontransmitbeamspace-PARAFAC[J].JournalofElectronics&InformationTechnology, 2011, 33(12): 2875-2880. (inChinese)

    [17]劉曉莉, 廖桂生. 基于MUSIC和ESPRIT的雙基地MIMO雷達角度估計算法[J]. 電子與信息學(xué)報, 2010, 32(9): 2179-2183.

    LIUXiaoli,LIAOGuisheng.CombinedMUSICwithESPRITalgorithmforangleestimationinbistaticMIMOradarsystem[J].JournalofElectronics&InformationTechnology, 2010, 32(9): 2179-2183. (inChinese)

    [18]CHENDF,CHENBX,QINGD.AngleestimationusingESPRITinMIMOradar[J].IETElectronicsLetters, 2008, 44(12): 770-771.

    黃中瑞(1988-),男,河南人,合肥電子工程學(xué)院博士研究生,主要研究方向是陣列信號處理、MIMO雷達信號處理及凸優(yōu)化理論.

    鄭志東(1985-),男,福建人,北方電子設(shè)備研究所工程師,主要研究方向是陣列信號處理、MIMO雷達信號處理.

    張劍云(1963-),男,浙江人,合肥電子工程學(xué)院教授,博士生導(dǎo)師,中國電子學(xué)會高級會員,雷達分委會委員,IEEE會員,主要研究方向是雷達及目標(biāo)環(huán)境模擬、雷達信號處理、高速信號處理.

    李迎林, 何丙發(fā), 張德斌. 雷達陣面系統(tǒng)小型化研究[J]. 電波科學(xué)學(xué)報,2015,30(4):797-802. doi: 10.13443/j.cjors. 2014091902

    LI Yinglin, HE Bingfa, ZHANG Debin, et al. Study on the integration technology of phased arrays[J]. Chinese Journal of Radio Science,2015,30(4):797-802. (in Chinese). doi: 10.13443/j.cjors. 2014091902

    久久国产精品人妻蜜桃| 一级毛片高清免费大全| 成人18禁在线播放| 国产中年淑女户外野战色| 久久久久久大精品| 一个人免费在线观看的高清视频| 国产高清激情床上av| 3wmmmm亚洲av在线观看| 国产精品久久久人人做人人爽| 丰满的人妻完整版| 成年女人永久免费观看视频| 亚洲精华国产精华精| 三级毛片av免费| 亚洲在线自拍视频| 我的老师免费观看完整版| 国产91精品成人一区二区三区| 免费搜索国产男女视频| 级片在线观看| 老汉色∧v一级毛片| 欧美午夜高清在线| 午夜福利视频1000在线观看| 高潮久久久久久久久久久不卡| 亚洲成人免费电影在线观看| 亚洲激情在线av| 香蕉丝袜av| 一边摸一边抽搐一进一小说| 欧美乱妇无乱码| 真人一进一出gif抽搐免费| 免费观看的影片在线观看| 十八禁人妻一区二区| 中国美女看黄片| 毛片女人毛片| 午夜福利视频1000在线观看| 国产精华一区二区三区| 别揉我奶头~嗯~啊~动态视频| 国产真实乱freesex| 在线观看日韩欧美| 国产97色在线日韩免费| 亚洲人成网站在线播| 日本三级黄在线观看| 国产综合懂色| 淫妇啪啪啪对白视频| 色尼玛亚洲综合影院| 国产日本99.免费观看| 在线观看免费午夜福利视频| 91久久精品国产一区二区成人 | 久久香蕉国产精品| 性色av乱码一区二区三区2| 国产精品嫩草影院av在线观看 | 国内少妇人妻偷人精品xxx网站| 亚洲精品日韩av片在线观看 | 桃红色精品国产亚洲av| 三级毛片av免费| 少妇的丰满在线观看| 一本精品99久久精品77| 亚洲国产日韩欧美精品在线观看 | 欧美zozozo另类| 香蕉丝袜av| 亚洲欧美一区二区三区黑人| 久久人人精品亚洲av| 日韩精品青青久久久久久| 高清日韩中文字幕在线| 男人和女人高潮做爰伦理| 九色国产91popny在线| 在线天堂最新版资源| 国产精品一及| 精品熟女少妇八av免费久了| 色综合婷婷激情| 日韩精品青青久久久久久| 淫秽高清视频在线观看| 成年女人毛片免费观看观看9| 熟妇人妻久久中文字幕3abv| 日韩欧美精品v在线| 香蕉丝袜av| 中出人妻视频一区二区| 国产老妇女一区| 国产高清激情床上av| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品成人综合色| 99在线视频只有这里精品首页| bbb黄色大片| 亚洲一区高清亚洲精品| 日本在线视频免费播放| 国产成人啪精品午夜网站| 日本撒尿小便嘘嘘汇集6| 欧美3d第一页| 亚洲国产精品久久男人天堂| 人妻丰满熟妇av一区二区三区| 757午夜福利合集在线观看| 国产亚洲精品综合一区在线观看| 一级毛片高清免费大全| 一级毛片女人18水好多| 日日摸夜夜添夜夜添小说| 亚洲av中文字字幕乱码综合| h日本视频在线播放| 成人国产一区最新在线观看| 日韩欧美一区二区三区在线观看| 日本精品一区二区三区蜜桃| 欧美性猛交╳xxx乱大交人| 床上黄色一级片| 日本 欧美在线| 亚洲av熟女| 欧美中文日本在线观看视频| av专区在线播放| 午夜福利在线观看免费完整高清在 | 我要搜黄色片| 午夜福利视频1000在线观看| 老司机午夜福利在线观看视频| 精品久久久久久久久久免费视频| 亚洲国产欧美网| 少妇熟女aⅴ在线视频| 国产乱人伦免费视频| 国产老妇女一区| 久久国产乱子伦精品免费另类| 久久这里只有精品中国| 亚洲熟妇中文字幕五十中出| 在线观看美女被高潮喷水网站 | 88av欧美| 搡老妇女老女人老熟妇| 热99re8久久精品国产| 国产高清三级在线| 欧美日本亚洲视频在线播放| 嫁个100分男人电影在线观看| 日韩大尺度精品在线看网址| 美女免费视频网站| 亚洲成av人片免费观看| 18+在线观看网站| 国产av在哪里看| 亚洲av第一区精品v没综合| 内地一区二区视频在线| ponron亚洲| 精品久久久久久久毛片微露脸| 国产精品亚洲av一区麻豆| 白带黄色成豆腐渣| 国产国拍精品亚洲av在线观看 | 999久久久精品免费观看国产| 亚洲国产欧美人成| 免费看日本二区| 亚洲av二区三区四区| 99热只有精品国产| 少妇人妻精品综合一区二区 | 久久香蕉精品热| 有码 亚洲区| 夜夜夜夜夜久久久久| 在线十欧美十亚洲十日本专区| 久久亚洲精品不卡| 国产精品美女特级片免费视频播放器| 手机成人av网站| 美女高潮喷水抽搐中文字幕| e午夜精品久久久久久久| 又黄又爽又免费观看的视频| 免费电影在线观看免费观看| 18禁黄网站禁片午夜丰满| 免费看光身美女| 他把我摸到了高潮在线观看| x7x7x7水蜜桃| 99视频精品全部免费 在线| 国产国拍精品亚洲av在线观看 | 又粗又爽又猛毛片免费看| 国产精品女同一区二区软件 | 久久精品国产99精品国产亚洲性色| 免费无遮挡裸体视频| 亚洲av免费在线观看| 黄色片一级片一级黄色片| 国产探花极品一区二区| www.999成人在线观看| 午夜激情欧美在线| 国产老妇女一区| 丰满乱子伦码专区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 搡老岳熟女国产| 欧美成人一区二区免费高清观看| 在线观看舔阴道视频| 国内揄拍国产精品人妻在线| 午夜激情欧美在线| 亚洲 欧美 日韩 在线 免费| 国产伦人伦偷精品视频| 啦啦啦韩国在线观看视频| 97超级碰碰碰精品色视频在线观看| 亚洲一区二区三区不卡视频| 五月玫瑰六月丁香| 亚洲av电影不卡..在线观看| 欧美+亚洲+日韩+国产| 亚洲国产精品成人综合色| 精品一区二区三区av网在线观看| 色噜噜av男人的天堂激情| 美女cb高潮喷水在线观看| 欧美高清成人免费视频www| 国内揄拍国产精品人妻在线| 在线看三级毛片| 变态另类成人亚洲欧美熟女| 久久久国产精品麻豆| 99热这里只有精品一区| 日韩欧美国产一区二区入口| 久久九九热精品免费| 99久久99久久久精品蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 欧美色视频一区免费| 欧美乱色亚洲激情| 亚洲精品色激情综合| av片东京热男人的天堂| 全区人妻精品视频| 亚洲七黄色美女视频| 久久久久久久亚洲中文字幕 | 精品免费久久久久久久清纯| 国产激情偷乱视频一区二区| 香蕉av资源在线| 国产中年淑女户外野战色| 在线观看66精品国产| 久久久久久大精品| 怎么达到女性高潮| 欧美午夜高清在线| 色av中文字幕| av中文乱码字幕在线| 在线观看午夜福利视频| www日本在线高清视频| 88av欧美| 高清在线国产一区| 国产午夜精品久久久久久一区二区三区 | 亚洲av电影不卡..在线观看| 大型黄色视频在线免费观看| 在线观看午夜福利视频| 久久久久国产精品人妻aⅴ院| 天堂动漫精品| 国产在视频线在精品| 亚洲真实伦在线观看| a在线观看视频网站| 国产伦在线观看视频一区| 88av欧美| 99久久成人亚洲精品观看| 97碰自拍视频| 国产欧美日韩一区二区精品| 九九在线视频观看精品| 香蕉丝袜av| 国内精品久久久久精免费| 午夜福利欧美成人| 亚洲成人久久爱视频| 变态另类成人亚洲欧美熟女| 欧美激情在线99| 中文字幕高清在线视频| 日韩欧美三级三区| 国产精品久久久久久久久免 | 美女免费视频网站| 午夜精品久久久久久毛片777| 国产精品国产高清国产av| 亚洲精品国产精品久久久不卡| 亚洲男人的天堂狠狠| 神马国产精品三级电影在线观看| 最近在线观看免费完整版| 国产精品一区二区三区四区免费观看 | 熟女少妇亚洲综合色aaa.| 欧美一区二区精品小视频在线| 中亚洲国语对白在线视频| 国产熟女xx| 国产高清视频在线观看网站| 制服人妻中文乱码| 日本 av在线| 国内精品久久久久久久电影| 国产一级毛片七仙女欲春2| 成人午夜高清在线视频| 国产精品亚洲美女久久久| 日韩欧美三级三区| 亚洲天堂国产精品一区在线| 久久亚洲精品不卡| 午夜福利成人在线免费观看| 无遮挡黄片免费观看| 亚洲人成电影免费在线| 亚洲一区二区三区色噜噜| 女人十人毛片免费观看3o分钟| 一区二区三区国产精品乱码| 亚洲熟妇熟女久久| 亚洲人成网站高清观看| 国产精品久久久人人做人人爽| 女人被狂操c到高潮| 精品人妻偷拍中文字幕| 久久久色成人| 999久久久精品免费观看国产| 91久久精品国产一区二区成人 | 国产欧美日韩精品亚洲av| 99热这里只有精品一区| 久久精品国产自在天天线| 在线观看免费午夜福利视频| 12—13女人毛片做爰片一| 免费看十八禁软件| 亚洲专区国产一区二区| 一个人观看的视频www高清免费观看| 脱女人内裤的视频| 美女大奶头视频| 免费搜索国产男女视频| 久久精品国产99精品国产亚洲性色| 婷婷亚洲欧美| 欧美在线黄色| 一个人观看的视频www高清免费观看| 国产又黄又爽又无遮挡在线| 小说图片视频综合网站| 变态另类丝袜制服| 老司机午夜十八禁免费视频| 757午夜福利合集在线观看| 美女黄网站色视频| 国产成人av教育| 成人高潮视频无遮挡免费网站| 他把我摸到了高潮在线观看| av女优亚洲男人天堂| 在线观看舔阴道视频| 搡女人真爽免费视频火全软件 | 亚洲人成伊人成综合网2020| 成人国产一区最新在线观看| 熟女人妻精品中文字幕| 老熟妇仑乱视频hdxx| 男人的好看免费观看在线视频| 欧美不卡视频在线免费观看| 丝袜美腿在线中文| 国产精品一区二区三区四区免费观看 | 国产精品一区二区免费欧美| 久久精品国产清高在天天线| 国产aⅴ精品一区二区三区波| 精华霜和精华液先用哪个| 精品熟女少妇八av免费久了| 日本免费一区二区三区高清不卡| 麻豆成人午夜福利视频| 久久久久久人人人人人| 一级作爱视频免费观看| 每晚都被弄得嗷嗷叫到高潮| 精品熟女少妇八av免费久了| ponron亚洲| 久久久国产成人免费| 精品国产超薄肉色丝袜足j| 欧美成人一区二区免费高清观看| 欧美日韩综合久久久久久 | 国产蜜桃级精品一区二区三区| 欧美成人免费av一区二区三区| 日本五十路高清| 中文字幕人妻熟人妻熟丝袜美 | 国产在线精品亚洲第一网站| 色吧在线观看| 亚洲在线自拍视频| 久久精品影院6| 久久久精品大字幕| 亚洲av免费在线观看| 制服人妻中文乱码| 亚洲18禁久久av| 国产在线精品亚洲第一网站| 国产精品一区二区三区四区久久| 亚洲欧美日韩高清在线视频| 国产av不卡久久| 国产成人a区在线观看| 丰满人妻一区二区三区视频av | 日韩 欧美 亚洲 中文字幕| 狂野欧美激情性xxxx| 婷婷精品国产亚洲av在线| 日韩有码中文字幕| 日韩欧美国产一区二区入口| 亚洲成a人片在线一区二区| av天堂在线播放| 精品无人区乱码1区二区| 精品电影一区二区在线| 51国产日韩欧美| 成人国产一区最新在线观看| 亚洲欧美激情综合另类| 亚洲欧美日韩无卡精品| av片东京热男人的天堂| 丁香欧美五月| 欧洲精品卡2卡3卡4卡5卡区| 香蕉av资源在线| 级片在线观看| 18禁国产床啪视频网站| 久久久久久久亚洲中文字幕 | 1024手机看黄色片| ponron亚洲| 可以在线观看的亚洲视频| 一进一出抽搐gif免费好疼| 午夜精品在线福利| 国产爱豆传媒在线观看| xxxwww97欧美| 久久国产乱子伦精品免费另类| 欧美丝袜亚洲另类 | 成人特级黄色片久久久久久久| 十八禁人妻一区二区| 在线a可以看的网站| 嫩草影院入口| 欧美不卡视频在线免费观看| 久久久久国内视频| 国产99白浆流出| tocl精华| 日本免费一区二区三区高清不卡| 国产日本99.免费观看| 嫩草影视91久久| 无人区码免费观看不卡| av中文乱码字幕在线| 成人鲁丝片一二三区免费| 美女大奶头视频| 亚洲av不卡在线观看| 久久精品夜夜夜夜夜久久蜜豆| 日本在线视频免费播放| 亚洲av免费在线观看| 国产色婷婷99| 国产精品一区二区三区四区久久| 日本在线视频免费播放| 国内揄拍国产精品人妻在线| 观看免费一级毛片| 男女做爰动态图高潮gif福利片| 国产美女午夜福利| 欧美+亚洲+日韩+国产| aaaaa片日本免费| 日本与韩国留学比较| www国产在线视频色| 少妇丰满av| 国内久久婷婷六月综合欲色啪| 18禁国产床啪视频网站| 国产激情欧美一区二区| 天天添夜夜摸| 国产精品国产高清国产av| 亚洲精华国产精华精| 中文亚洲av片在线观看爽| 欧美乱妇无乱码| 欧美一区二区精品小视频在线| 热99re8久久精品国产| 日本与韩国留学比较| 中出人妻视频一区二区| 精品国产美女av久久久久小说| 国产高清三级在线| 久久久国产精品麻豆| 国产伦在线观看视频一区| 最近最新中文字幕大全免费视频| 免费av不卡在线播放| 十八禁人妻一区二区| 三级国产精品欧美在线观看| 一级黄片播放器| 亚洲欧美日韩无卡精品| www.熟女人妻精品国产| 波多野结衣高清作品| 国产av麻豆久久久久久久| 久久久精品欧美日韩精品| 老司机深夜福利视频在线观看| 日韩av在线大香蕉| 日韩中文字幕欧美一区二区| 国产美女午夜福利| 亚洲成人中文字幕在线播放| 中文字幕人妻熟人妻熟丝袜美 | 一个人看视频在线观看www免费 | 久久久久久九九精品二区国产| 国产中年淑女户外野战色| 黄片大片在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 久久亚洲真实| 老汉色∧v一级毛片| 中文字幕精品亚洲无线码一区| 亚洲真实伦在线观看| 成熟少妇高潮喷水视频| 久久久久久久久久黄片| 午夜日韩欧美国产| 两个人视频免费观看高清| 精品不卡国产一区二区三区| 亚洲成人久久爱视频| 色在线成人网| 精品人妻一区二区三区麻豆 | 日韩欧美国产在线观看| 国产精品日韩av在线免费观看| 国产成+人综合+亚洲专区| 夜夜爽天天搞| 成人av一区二区三区在线看| 欧美中文日本在线观看视频| 欧美午夜高清在线| 一级毛片高清免费大全| 久久久久久久精品吃奶| 亚洲中文字幕一区二区三区有码在线看| 午夜激情欧美在线| 男女下面进入的视频免费午夜| 国产色爽女视频免费观看| 麻豆国产97在线/欧美| 又黄又粗又硬又大视频| 日韩欧美精品免费久久 | 欧美性猛交黑人性爽| 欧美绝顶高潮抽搐喷水| 久久久久国内视频| 97碰自拍视频| 法律面前人人平等表现在哪些方面| 国产成年人精品一区二区| 国产亚洲精品综合一区在线观看| 男女视频在线观看网站免费| 999久久久精品免费观看国产| 一级黄片播放器| 男人舔女人下体高潮全视频| 啦啦啦观看免费观看视频高清| 最近视频中文字幕2019在线8| 成年版毛片免费区| 亚洲国产色片| eeuss影院久久| 精品一区二区三区av网在线观看| 欧美性感艳星| 又紧又爽又黄一区二区| 久99久视频精品免费| 亚洲国产精品sss在线观看| 中文在线观看免费www的网站| 欧美一区二区国产精品久久精品| 亚洲国产精品合色在线| 麻豆国产av国片精品| 人人妻人人澡欧美一区二区| 日韩欧美国产在线观看| 噜噜噜噜噜久久久久久91| 欧美日韩精品网址| 国产精品美女特级片免费视频播放器| 色噜噜av男人的天堂激情| 高潮久久久久久久久久久不卡| 叶爱在线成人免费视频播放| 九九在线视频观看精品| 国产视频内射| 国产精品,欧美在线| 亚洲中文日韩欧美视频| 中亚洲国语对白在线视频| 国内少妇人妻偷人精品xxx网站| 性欧美人与动物交配| 国产私拍福利视频在线观看| 亚洲成av人片在线播放无| 欧美xxxx黑人xx丫x性爽| 国产乱人视频| 欧美国产日韩亚洲一区| 亚洲av日韩精品久久久久久密| 激情在线观看视频在线高清| 日韩欧美国产在线观看| 乱人视频在线观看| 在线天堂最新版资源| 99久久精品国产亚洲精品| www.色视频.com| 人人妻人人看人人澡| 一个人观看的视频www高清免费观看| 1000部很黄的大片| 成熟少妇高潮喷水视频| 亚洲精品国产精品久久久不卡| 岛国在线观看网站| 99在线视频只有这里精品首页| 一级a爱片免费观看的视频| 夜夜躁狠狠躁天天躁| 一夜夜www| 色哟哟哟哟哟哟| 国产v大片淫在线免费观看| 亚洲熟妇中文字幕五十中出| 成人高潮视频无遮挡免费网站| 国产欧美日韩一区二区精品| 色噜噜av男人的天堂激情| 美女被艹到高潮喷水动态| 日韩av在线大香蕉| 91字幕亚洲| 精品99又大又爽又粗少妇毛片 | 国产高清视频在线观看网站| 亚洲va日本ⅴa欧美va伊人久久| 久久久色成人| 少妇丰满av| 麻豆一二三区av精品| 欧美成人a在线观看| 此物有八面人人有两片| 日韩欧美国产在线观看| 我要搜黄色片| 精品一区二区三区人妻视频| 成人亚洲精品av一区二区| 成年女人永久免费观看视频| av福利片在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品,欧美在线| 又黄又粗又硬又大视频| 国产欧美日韩一区二区精品| 一区福利在线观看| 日韩欧美在线乱码| 美女cb高潮喷水在线观看| 91av网一区二区| 成人一区二区视频在线观看| 日本 欧美在线| 女人高潮潮喷娇喘18禁视频| 国产在线精品亚洲第一网站| 亚洲七黄色美女视频| 午夜福利高清视频| 欧美国产日韩亚洲一区| 深爱激情五月婷婷| 国产99白浆流出| 校园春色视频在线观看| 90打野战视频偷拍视频| 国产精品久久久久久久电影 | 日本精品一区二区三区蜜桃| 久久久久国产精品人妻aⅴ院| 18禁黄网站禁片午夜丰满| 精品国产美女av久久久久小说| 国产精品亚洲美女久久久| 亚洲欧美日韩高清专用| 国产不卡一卡二| 欧美高清成人免费视频www| 禁无遮挡网站| 国产精品香港三级国产av潘金莲| 18+在线观看网站| 欧美乱妇无乱码| 久久精品夜夜夜夜夜久久蜜豆| 18+在线观看网站| 亚洲最大成人手机在线| 啦啦啦免费观看视频1| 真人做人爱边吃奶动态| 免费看十八禁软件| 亚洲国产精品sss在线观看| 国产av麻豆久久久久久久| 国模一区二区三区四区视频| 在线视频色国产色| 看片在线看免费视频| 亚洲精品久久国产高清桃花| 欧美av亚洲av综合av国产av| 免费人成在线观看视频色| 一边摸一边抽搐一进一小说| 美女大奶头视频| 波多野结衣巨乳人妻| 久久性视频一级片| 色噜噜av男人的天堂激情| 亚洲国产日韩欧美精品在线观看 | 在线视频色国产色| 国产伦精品一区二区三区四那| 一边摸一边抽搐一进一小说| 亚洲精品日韩av片在线观看 | 十八禁网站免费在线| 亚洲专区国产一区二区|