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

    基于nested陣列的高分辨DOA估計(jì)*

    2016-10-19 02:10:04劉魯濤
    關(guān)鍵詞:接收數(shù)據(jù)協(xié)方差特征值

    劉魯濤, 韓 麗, 張 明, 張 航

    (哈爾濱工程大學(xué) 信息與通信工程學(xué)院, 哈爾濱 150001)

    ?

    基于nested陣列的高分辨DOA估計(jì)*

    劉魯濤, 韓麗, 張明, 張航

    (哈爾濱工程大學(xué) 信息與通信工程學(xué)院, 哈爾濱 150001)

    針對(duì)nested陣列對(duì)鄰近信號(hào)的分辨力受信噪比和快拍數(shù)等因素限制的問題,提出了基于nested陣列的加權(quán)子空間平滑MUSIC算法.該算法對(duì)協(xié)方差矩陣向量化以提高整個(gè)陣列的自由度,使用空間平滑恢復(fù)新接收數(shù)據(jù)矢量陣的秩,采用校正的噪聲特征值對(duì)噪聲子空間進(jìn)行加權(quán),并對(duì)信號(hào)子空間進(jìn)行空間譜合成,得到新算法的空間譜函數(shù).通過搜索空間譜函數(shù)極大值實(shí)現(xiàn)DOA估計(jì).結(jié)果表明,該算法在低信噪比及小快拍數(shù)條件下,對(duì)間隔較近的信號(hào)具有高分辨力.

    nested陣列; 協(xié)方差矩陣; 向量化; 自由度; 空間平滑; 加權(quán); DOA估計(jì); 高分辨力

    波達(dá)方向(direction-of-arrival,DOA)估計(jì)是陣列信號(hào)處理領(lǐng)域的一個(gè)重要研究方向,被廣泛應(yīng)用于雷達(dá)、聲源定位和通信等領(lǐng)域[1-3].子空間類算法[4-7]已經(jīng)成為DOA估計(jì)中最常用和最經(jīng)典的算法,尤其以多重信號(hào)分類[8-10](multiple signal classification,MUSIC)為典型代表的算法隨著信噪比的提高可以突破瑞利限,但其缺點(diǎn)是自由度不高,估計(jì)目標(biāo)數(shù)最多為N-1(N為陣元數(shù)目).為了提高陣列的自由度,文獻(xiàn)[11]使用最小冗余陣列通過構(gòu)建增強(qiáng)協(xié)方差矩陣來增加陣列的自由度,但在實(shí)際應(yīng)用中,這種陣列的最優(yōu)陣元位置受限于計(jì)算機(jī)建模以及迭代算法的復(fù)雜程度;文獻(xiàn)[12]基于高階累計(jì)量的DOA估計(jì)算法實(shí)現(xiàn)陣列擴(kuò)展,增加虛擬陣元,但是其缺點(diǎn)是需要很大運(yùn)算量來計(jì)算信號(hào)的高階累積量,需要較多的快拍數(shù)才能準(zhǔn)確地估計(jì)出目標(biāo)信號(hào)參數(shù),且不能用于處理高斯信號(hào)的DOA估計(jì);文獻(xiàn)[13]針對(duì)非平穩(wěn)信號(hào)的特性提出了一種Khatrio-Rao子空間估計(jì)算法,實(shí)現(xiàn)了在陣元數(shù)少于信號(hào)源數(shù)目情況下對(duì)準(zhǔn)平穩(wěn)信號(hào)的DOA估計(jì),但其測(cè)角精度不高,且也不適用于平穩(wěn)信號(hào);在此基礎(chǔ)上,文獻(xiàn)[14]利用Khatrio-Rao(KR)積的性質(zhì)提出了一種新的陣列構(gòu)造方式,即nested陣列,使得陣列的自由度大大增加,對(duì)于M個(gè)物理陣元的陣列可以獲得O(M2)個(gè)自由度,并提出了相應(yīng)的基于空間平滑多重信號(hào)分類(spatial smoothing MUSIC,SS MUSIC)的DOA估計(jì)方法.雖然nested陣列可以提高陣列的自由度,解決了在信源數(shù)多于物理陣元數(shù)情況下對(duì)平穩(wěn)信號(hào)的DOA估計(jì),但是其在低信噪比及小快拍數(shù)的條件下,并沒有提高對(duì)鄰近入射信號(hào)的分辨率.

    本文將噪聲子空間與信號(hào)特征值倒數(shù)加權(quán)的信號(hào)子空間相結(jié)合,得到加權(quán)子空間平滑MUSIC(weighted subspace of spatial smoothing,WSSS MUSIC)算法的譜函數(shù),充分利用了信號(hào)的信息,降低了噪聲的影響,提高了nested陣列鄰近信號(hào)的分辨力.

    1 數(shù)據(jù)模型和陣列結(jié)構(gòu)模型

    1.1數(shù)據(jù)模型

    假設(shè)L個(gè)遠(yuǎn)場(chǎng)窄帶信號(hào)入射到M個(gè)陣元的線陣上,信源入射角度分別為θ1,θ2,…,θL,入射信號(hào)彼此相互獨(dú)立.噪聲是高斯白噪聲與信號(hào)相互獨(dú)立,則陣列接收數(shù)據(jù)矢量形式為

    X(t)=AS(t)+N(t)

    (1)

    式中:X(t)為陣列的M×1維陣列接收數(shù)據(jù)矢量;S(t)為空間信號(hào)的L×1維矢量;N(t)為M×1維的噪聲矢量;A=[a(θ1),a(θ2),…,a(θL)]為M×L維的陣列流型,其中

    a(θk)=[exp(-j2πd1sinθk/λ),

    exp(-j2πd2sinθk/λ),…,

    exp(-j2πdMsinθk/λ)]T

    (2)

    式中:di為第i個(gè)陣元相對(duì)陣列參考點(diǎn)的位置;λ為波長(zhǎng),則陣列接收數(shù)據(jù)的協(xié)方差矩陣為

    (3)

    (4)

    Khatri-Rao積的定義:用符號(hào)“⊙”表示Khatri-Rao(KR)積,兩個(gè)矩陣A(N×K)和B(M×K)具有相同的列數(shù),它們的KR積[13]為A⊙B=[a1?b1,a2?b2,…,aK?bK]=

    (5)

    式中:ak,bk分別為矩陣A和B的第k列;vec(·)為將矩陣的每一列元素排成一列,形成一個(gè)新的列向量.

    KR積具有如下性質(zhì):

    若A∈CN×K,B∈CM×K和D=diag(d),且d∈CK×1,則

    vec(ADBH)=(B*⊙A)d

    (6)

    式中,*為復(fù)共軛.

    利用上述性質(zhì),且假設(shè)信號(hào)是相互獨(dú)立的,則式(3)可化簡(jiǎn)為

    z=vec(RXX)=

    (7)

    (8)

    式中,

    (i=1,2,…,M;j=1,2,…,M)

    (9)

    由式(9)可看出,此陣列流型類似于多輸入多輸出(multiple input multiple output,MIMO)雷達(dá)[15],此陣列的陣元位置由集合{di-dj,1≤i,j≤M}(di,dj分別為第i個(gè)和第j個(gè)陣元的位置)中相異的值來決定,其相異值遠(yuǎn)大于物理陣元數(shù)M,因此增大了陣列的自由度.

    1.2二階nested陣列的陣列結(jié)構(gòu)

    在陣元數(shù)固定的前提下,為了使相異值最多(自由度最大),文獻(xiàn)[11]提出了一種基于KR積的陣列構(gòu)造方式,稱其為nested陣列.一個(gè)標(biāo)準(zhǔn)的二階nested陣列是由兩個(gè)相鄰的均勻線陣(uniform linear arrays,ULA)串聯(lián)構(gòu)成的,內(nèi)層的ULA有M1個(gè)陣元,陣元間距為d1,外層的ULA有M2個(gè)陣元,其陣元間距為d2=(M1+1)d1,則第1個(gè)ULA的陣元位置為S1={nd1,n=1,2,…,M1},第2個(gè)ULA的陣元位置為S2={m(M1+1)d1,m=1,2,…,M2},一個(gè)二階nested陣列陣元總數(shù)為M=M1+M2,其結(jié)構(gòu)示意圖如圖1所示.

    圖1 二階nested陣列的結(jié)構(gòu)示意圖Fig.1 Schematic structure of second-order nested array

    一個(gè)陣元總數(shù)為M的nested陣列,由KR積可產(chǎn)生(M2-2)/2+M個(gè)虛擬陣元,此虛擬陣列為陣元間距d1的ULA,其陣元位置集合為

    S={ld1,l=-N,-(N-1),…,N;

    N=M2(M1+1)-1}

    (10)

    這樣便增加了nested陣列的自由度,為實(shí)現(xiàn)信號(hào)源數(shù)目多于物理陣元總數(shù)的DOA估計(jì)提供了可能.

    2 波達(dá)方向估計(jì)

    上文對(duì)陣列接收數(shù)據(jù)協(xié)方差矩陣進(jìn)行了向量化,使得到的新接收數(shù)據(jù)矢量z的秩變?yōu)?,并且由KR積產(chǎn)生的虛擬陣列為ULA,用類似空間平滑算法[16]恢復(fù)矩陣z的秩.首先,從構(gòu)造的新陣列流型(A*⊙A)∈CM2×L中移除相同的行,然后對(duì)其排序,使第i行的元素對(duì)應(yīng)于虛擬陣列的陣元位置為(-M2/4-M/2+i)d1,得到矩陣A1∈C((M2-2)/2+M)×L,這個(gè)操作類似移除矩陣z中相應(yīng)的行,排序得到新的矩陣為

    (11)

    式中,e0∈C((M2-2)/2+M)×1,為除了第(M2/4+M/2)個(gè)元素是1,其余元素均為0的列向量.最后把矩陣z1分成(M2/4+M/2)個(gè)重疊的子陣,其中,每個(gè)子陣有(M2/4+M/2)個(gè)元素,第i個(gè)子陣的陣元位置為

    (12)

    第i個(gè)子陣對(duì)應(yīng)的陣列接收數(shù)據(jù)矢量為

    (13)

    式中,矩陣A1i為(M2/4+M/2)×L維,其包含矩陣A1中從第(M2/4+M/2-i+1)行到第((M2-2)/2+M-i+1)行的元素.

    由上面的分析容易得到

    (14)

    式中,

    Φ=diag(exp(-jφ1),exp(-jφ2),…,exp(-jφL))

    其中,φk=2πd1sinθk/λ.令

    Ri

    (15)

    對(duì)Ri進(jìn)行平均,可得到平滑后的新接收數(shù)據(jù)的協(xié)方差矩陣為

    RSS

    (16)

    平滑后的ULA子陣陣元數(shù)MSS=M2/4+M/2,且平滑后的陣列流型為

    ASS=[aSS(θ1),aSS(θ2),…,aSS(θL)]

    (17)

    式中,

    (18)

    由此可看出,此陣列流型同有MSS個(gè)陣元的ULA陣列流型相同.于是可以構(gòu)造如下SS MUSIC算法的空間譜函數(shù)

    (19)

    式中,EN為最小特征值對(duì)應(yīng)的噪聲子空間.

    3 WSSS MUSIC算法

    (20)

    利用校正后的噪聲特征值對(duì)噪聲子空間進(jìn)行加權(quán)處理,可以提高噪聲特征值信息量的利用率,而且校正值μ不影響信號(hào)特征值,可以保證目標(biāo)信號(hào)的波達(dá)方向估計(jì)穩(wěn)定且準(zhǔn)確.由式(20)可知,校正值μ的選取會(huì)對(duì)算法的估計(jì)性能產(chǎn)生很大的影響,若μ值選取的過小,會(huì)使得對(duì)噪聲的抑制作用不明顯;若μ值選取的過大,會(huì)使信號(hào)特征值與噪聲特征值的差別太小,降低了目標(biāo)信號(hào)波達(dá)方向估計(jì)的準(zhǔn)確率,因此選擇一個(gè)最優(yōu)的校正值μ十分重要.從實(shí)際的大量實(shí)驗(yàn)數(shù)據(jù)可知,噪聲特征值中的最大值與其最小值之比小于2,即

    λL+1/λMSS<2

    (21)

    將式(20)代入式(21)可得

    (22)

    (23)

    (24)

    將式(24)代入式(19)中得到校正后的噪聲子空間的譜函數(shù)為

    (25)

    在非理想環(huán)境下,通過校正值μ可以使噪聲特征值的發(fā)散程度受到抑制,可以在一定程度上抑制噪聲的影響,提高算法的分辨率.對(duì)噪聲子空間進(jìn)行加權(quán)也只是利用了噪聲子空間的信息,卻沒有充分利用信號(hào)子空間的信息.由于特征分解協(xié)方差矩陣得到的信號(hào)子空間受噪聲起伏的影響相對(duì)較小,對(duì)非理想環(huán)境下引起的誤差不敏感,因此,為了充分利用新接收數(shù)據(jù)的信號(hào)子空間所包含的信息,可以利用主特征值的倒數(shù)對(duì)信號(hào)子空間進(jìn)行加權(quán),得到的譜函數(shù)表達(dá)式為

    (26)

    WSSS MUSIC算法不僅降低了噪聲的影響,而且有效地利用了新接收數(shù)據(jù)信號(hào)子空間的信息,使得算法在低信噪比和少快拍數(shù)的非理想情況下,對(duì)鄰近信號(hào)具有較高的分辨率,其空間譜函數(shù)為

    (27)

    算法的具體步驟如下:

    1) 計(jì)算快拍數(shù)據(jù)協(xié)方差矩陣RXX,將其向量化,移除相同的行并排序,得到新的陣列接收數(shù)據(jù)矢量z;

    2) 用空間平滑技術(shù)恢復(fù)矢量陣z的秩,并求得協(xié)方差矩陣RSS;

    3) 對(duì)RSS進(jìn)行特征值分解,將最大特征值對(duì)應(yīng)的特征矢量構(gòu)成的子空間定義為信號(hào)子空間ES=[v1,v2,…,vL],最小特征值對(duì)應(yīng)的特征矢量構(gòu)成的子空間定義為噪聲子空間EN=[vL+1,vL+2,…,vMSS];

    4) 根據(jù)式(27)求出WSSS MUSIC算法的空間譜函數(shù)進(jìn)行譜峰搜索.

    4 仿真實(shí)驗(yàn)

    實(shí)驗(yàn)1考慮兩個(gè)非相干的等強(qiáng)信號(hào)入射二階nested陣列,入射方向分別是10°和12°,信噪比為2 dB,快拍數(shù)為200,本文算法和SS MUSIC算法的空間譜圖比較如圖2所示.

    圖2 兩種算法的空間譜比較Fig.2 Comparison in spatial spectra of two algorithms

    由圖2可以看出,在低信噪比和小快拍數(shù)條件下,本文方法能夠準(zhǔn)確檢測(cè)出兩個(gè)非常鄰近的信號(hào)(10°,12°),而SS MUSIC算法的兩個(gè)鄰近信號(hào)的譜峰在空間譜上混疊成一個(gè)譜峰,算法失效.

    實(shí)驗(yàn)2考慮兩個(gè)非相干等強(qiáng)信號(hào)入射二階nested陣列,入射角分別為θ1=10°,θ2=θ1+Δθ,其中,Δθ∈[1°,8°],信噪比為0 dB,快拍數(shù)為200,本文方法和SS MUSIC算法的分辨概率隨角度間隔變化的曲線如圖3所示.

    圖3 分辨概率隨角度間隔的變化曲線Fig.3 Change curves of resolution probability with angle interval

    由圖3可以看出,在低信噪比及小快拍數(shù)條件下,在角度間隔(小于5°)相同時(shí),本文算法的分辨概率要比SS MUSIC算法高.而在識(shí)別成功率為1的前提下,本文算法有著更好的角度分辨率(角度間隔為4°).

    實(shí)驗(yàn)3兩個(gè)非相干等強(qiáng)信號(hào)源入射方向分別為10°和12°,在使用快拍數(shù)為200的條件下,設(shè)置信噪比從-4 dB按步長(zhǎng)1 dB的間隔增加到10 dB時(shí),進(jìn)行1 000次Monte-Carlo實(shí)驗(yàn).對(duì)本文算法、SS MUSIC算法以及12個(gè)陣元ULA的MUSIC算法成功分辨概率進(jìn)行比較如圖4所示.

    圖4 三種算法成功分辨概率隨信噪比的變化曲線Fig.4 Change curves of success resolution probability with SNR for three algorithms

    由圖4可以看出,在快拍數(shù)一定的條件下,隨著信噪比的增加,三種算法鄰近信號(hào)成功分辨率均隨之增大.但在小信噪比的時(shí)候,本文算法相對(duì)兩個(gè)鄰近信號(hào)分辨的成功概率最高.而SS MUSIC算法由于受信噪比及快拍數(shù)的限制無法使兩個(gè)鄰近信號(hào)(10°和12°)的成功分辨概率達(dá)到1,進(jìn)一步說明了本文算法的分辨性能.

    實(shí)驗(yàn)4信源入射信號(hào)度不變,信噪比為5 dB,設(shè)置快拍數(shù)從1按步長(zhǎng)40增加到601,進(jìn)行1 000次Monte-Carlo實(shí)驗(yàn).本文算法、SS MUSIC算法及12個(gè)陣元ULA的MUSIC算法鄰近信號(hào)成功分辨概率隨快拍數(shù)變化曲線如圖5所示.

    圖5 三種算法成功分辨概率隨快拍數(shù)的變化曲線Fig.5 Change curves of success resolution probability with number of snapshots for three algorithms

    由圖5可以看出,在信噪比一定的條件下,三種算法對(duì)鄰近信號(hào)的成功分辨概率均隨著快拍數(shù)的增加而增大.而在快拍數(shù)較少時(shí),WSSS MUSIC算法對(duì)兩個(gè)鄰近信號(hào)的成功分辨概率明顯優(yōu)于其他兩種算法,在快拍數(shù)較大時(shí),WSSS MUSIC算法的成功分辨概率幾乎和12個(gè)陣元ULA的MUSIC算法同時(shí)達(dá)到1.

    5 結(jié) 論

    本文提出了一種WSSS MUSIC算法,該方法通過對(duì)協(xié)方差矩陣向量化產(chǎn)生新接收數(shù)據(jù),并將噪聲特征值加權(quán)子空間與加權(quán)信號(hào)子空間相結(jié)合,充分利用新接收數(shù)據(jù)協(xié)方差的特征值信息.仿真結(jié)果表明,在低信噪比和小快拍數(shù)條件下,WSSS MUSIC算法對(duì)鄰近信號(hào)的分辨性能優(yōu)于SS MUSIC算法和具有相同物理孔徑的12陣元ULA的MUSIC算法,在實(shí)際應(yīng)用中節(jié)省了陣元數(shù),降低了硬件成本.

    [1]Hassanien A,Vorobyov S.Transmit energy focusing for DOA estimation in MIMO radar with collocated antennas [J].IEEE Transactions on Signal Processing,2011,59(6):2669-2682.

    [2]初萍,司偉建.基于延時(shí)相關(guān)預(yù)處理的MUSIC算法 [J].沈陽工業(yè)大學(xué)學(xué)報(bào),2013,35(1):93-98.

    (CHU Ping,SI Wei-jian.MUSIC algorithm based on delay correlation preprocessing [J].Journal of Shen-yang University of Technology,2013,35(1):93-98.)

    [3]Mousavi P,Shafai L,Veidt B,et al.Feed-reflector design for large adaptive reflector antenna (LAR) [J].IEEE Transactions on Antennas and Propagation,2001,49(8):1142-1154.

    [4]Andersen J B.Array gain and capacity for known random channels with multiple element array at both ends [J].IEEE Journal on Selected Areas in Communications,2000,18(11):2172-2178.

    [5]王進(jìn),盧影,孫開偉,等.基于演化超網(wǎng)絡(luò)的DNA微陣列數(shù)據(jù)分類方法 [J].重慶郵電大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,26(5):679-685.

    (WANG Jin,LU Ying,SUN Kai-wei,et al.DNA microarray data classification based on evolutionary hypernetworks [J].Journal of Chongqing University of Posts and Telecommunications(Natural Science Edition),2014,26(5):679-685.)

    [6]刁鳴,陳超,楊麗麗.四階累積量陣列擴(kuò)展的傳播算子測(cè)向方法 [J].哈爾濱工程大學(xué)學(xué)報(bào),2010,31(5):652-656.

    (DIAO Ming,CHEN Chao,YANG Li-li.Propagator method for direction of arrival estimation based on an array extension of the fourth-order cumulant [J].Journal of Harbin Engineering University,2010,31(5):652-656.)

    [7]Qian C,Huang L,So H C.Improved unitary root-MUSIC for DOA estimation based on pseudo-noise resampling [J].Signal Processing Letters,2014,21(2):140-144.

    [8]Grover R,Pados A,Medley M J.Subspace direction finding with an auxiliary vector basis [J].IEEE Transactions on Signal Processing,2007,55(2):758-763.

    [9]Birot G,Albera L,Chevalier P.Sequential high resolu-tion direction finding from higher order statistics [J].IEEE Transactions on Signal Processing,2010,58(8):4144-4155.

    [10]Feng G Y,Ming J,Liu S,et al.Real-valued MUSIC for efficient direction estimation with arbitrary array geometries [J].IEEE Transaction on Signal Processing,2014,62(6):1548-1560.

    [11]Gu J F,Zhu W P,Swamy M N S.Minimum redundancy linear sparse subarrays for direction of arrival estimation without ambiguity [C]//2011 IEEE International Symposium on Circuits and Systems (ISCAS).Rio de Janeiro,Brazil,2011:390-393.

    [12]Cao S,Ye Z,Hu N,et al.DOA estimation based on fourth order cumulants in the presence of sensor gain-phase errors [J].Signal Processing,2013,93(9):2581-2585.

    [13]Ma W,Hsieh T,Chi C.DOA estimation of quasistationary signals with less sensors than sources and unknown spatial noise covariance:a Khatri-Rao subspace approch [J].IEEE Transaction on Signal Processing,2010,58(4):2168-2180.

    [14]Han K,Nehorai.An improved source number detection and direction estimation with nested arrays and ULAs using jackknifing [J].IEEE Transaction on Signal Processing,2013,61(23):6118-6128.

    [15]Xie R,Liu Z,Zhang Z J.DOA estimation for mono static MIMO radar using polynomial rooting [J].Signal Processing,2010,90(12):3284-3288.

    [16]Cai B,Li Y M,Wang H Y.Forward/backward spatial reconstruction method for directions of arrival estimation of uncorrelated and coherent signals [J].IET Microwaves,Antennas and Propagation,2012,13(6):1498-1505.

    (責(zé)任編輯:景勇英文審校:尹淑英)

    DOA estimation with high resolution based on nested array

    LIU Lu-tao, HAN Li, ZHANG Ming, ZHANG Hang

    (College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China)

    In order to solve the problem that the resolution of nested array for adjacent signals is limited by such factors as signal to noise ratio (SNR) and number of snapshots, a weighted subspace smoothing MUSIC algorithm was proposed based on the nested array. The covariance matrix was vectorized in the algorithm to improve the degrees of freedom of whole array. The spatial smoothing was used to recovery the rank of newly received data vector matrix, and the corrected noise eigenvalues were adopted to weight the noise subspace. In addition, the spatial spectrum of signal subspace was synthesized, and the spatial spectrum function of new algorithm was obtained. Through searching the maximum value of spatial spectrum function, the DOA estimation was realized. The results show that under the condition of low SNR and small number of snapshots, the algorithm has high resolution of signals with relatively close interval.

    nested array; covariance matrix; vectorization; degrees of freedom; spatial smoothing; weighting; DOA estimation; high resolution

    2015-10-15.

    國(guó)家自然科學(xué)基金資助項(xiàng)目(61202410).

    劉魯濤(1977-),男,黑龍江哈爾濱人,副教授,博士,主要從事寬帶信號(hào)處理、檢測(cè)與識(shí)別等方面的研究.

    10.7688/j.issn.1000-1646.2016.05.10

    TN 911.6

    A

    1000-1646(2016)05-0531-06

    *本文已于2016-05-12 13∶56在中國(guó)知網(wǎng)優(yōu)先數(shù)字出版. 網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/21.1189.T.20160512.1356.008.html

    猜你喜歡
    接收數(shù)據(jù)協(xié)方差特征值
    一類帶強(qiáng)制位勢(shì)的p-Laplace特征值問題
    沖激噪聲背景下基于幅度預(yù)處理的測(cè)向新方法*
    單圈圖關(guān)聯(lián)矩陣的特征值
    低復(fù)雜度多輸入多輸出雷達(dá)目標(biāo)角度估計(jì)方法
    單片機(jī)模擬串口數(shù)據(jù)接收程序的實(shí)現(xiàn)及優(yōu)化
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    基于商奇異值分解的一類二次特征值反問題
    關(guān)于兩個(gè)M-矩陣Hadamard積的特征值的新估計(jì)
    縱向數(shù)據(jù)分析中使用滑動(dòng)平均Cholesky分解對(duì)回歸均值和協(xié)方差矩陣進(jìn)行同時(shí)半?yún)?shù)建模
    欧美亚洲日本最大视频资源| 晚上一个人看的免费电影| 国产黄色视频一区二区在线观看| 亚洲精品成人av观看孕妇| 啦啦啦啦在线视频资源| 自线自在国产av| av黄色大香蕉| 亚洲欧美色中文字幕在线| 午夜视频国产福利| 国产在线免费精品| .国产精品久久| 成年美女黄网站色视频大全免费 | 春色校园在线视频观看| 久久久久精品性色| 水蜜桃什么品种好| 久久99热6这里只有精品| 黑人巨大精品欧美一区二区蜜桃 | 亚洲,欧美,日韩| 成人午夜精彩视频在线观看| 久久狼人影院| 亚洲不卡免费看| 香蕉精品网在线| 亚洲美女搞黄在线观看| videosex国产| 精品一区二区免费观看| 成人无遮挡网站| 亚洲图色成人| 免费观看av网站的网址| 国产国拍精品亚洲av在线观看| 99久久精品国产国产毛片| 91成人精品电影| 国产精品一区二区在线观看99| 欧美bdsm另类| 99九九在线精品视频| 自线自在国产av| 久久久久久久久久久免费av| 亚洲美女视频黄频| 欧美 亚洲 国产 日韩一| 精品人妻一区二区三区麻豆| 国产成人精品久久久久久| 七月丁香在线播放| 亚洲精品日本国产第一区| 一个人免费看片子| 欧美精品一区二区大全| 搡老乐熟女国产| 狂野欧美白嫩少妇大欣赏| 在线观看免费视频网站a站| 我的老师免费观看完整版| 国产精品人妻久久久影院| 99久久精品国产国产毛片| 在线观看人妻少妇| 成人毛片a级毛片在线播放| 大片电影免费在线观看免费| 九色亚洲精品在线播放| 欧美激情 高清一区二区三区| 欧美精品人与动牲交sv欧美| 制服人妻中文乱码| 蜜桃国产av成人99| 少妇高潮的动态图| 国产精品秋霞免费鲁丝片| 免费观看无遮挡的男女| 国产欧美日韩一区二区三区在线 | 欧美日韩国产mv在线观看视频| 成年av动漫网址| 少妇高潮的动态图| 亚洲伊人久久精品综合| 国产精品嫩草影院av在线观看| 国产日韩欧美在线精品| 亚洲av国产av综合av卡| 一区二区三区四区激情视频| 日日啪夜夜爽| 午夜福利网站1000一区二区三区| 多毛熟女@视频| 午夜激情福利司机影院| 伊人久久精品亚洲午夜| 少妇被粗大猛烈的视频| 一本久久精品| 免费少妇av软件| 人妻夜夜爽99麻豆av| 精品久久久久久久久av| 久久久精品区二区三区| 免费av不卡在线播放| 久久久久久伊人网av| 丰满少妇做爰视频| 久久精品熟女亚洲av麻豆精品| 亚洲av欧美aⅴ国产| 有码 亚洲区| 欧美xxxx性猛交bbbb| 国产精品人妻久久久久久| 亚洲精品久久久久久婷婷小说| 在线观看国产h片| 热re99久久精品国产66热6| 亚洲精品美女久久av网站| 一级毛片黄色毛片免费观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久久精品电影小说| 黄色欧美视频在线观看| 成人手机av| 国产av国产精品国产| 在线亚洲精品国产二区图片欧美 | 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 亚洲av男天堂| 大香蕉97超碰在线| 午夜免费男女啪啪视频观看| 全区人妻精品视频| 国产熟女欧美一区二区| 大话2 男鬼变身卡| 亚洲精品日韩在线中文字幕| 亚洲av成人精品一二三区| 中国国产av一级| 午夜免费鲁丝| 国内精品宾馆在线| 国产乱来视频区| 亚洲人与动物交配视频| 精品人妻熟女av久视频| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线观看播放| 在线天堂最新版资源| 亚洲色图 男人天堂 中文字幕 | 国产欧美日韩综合在线一区二区| 国产亚洲午夜精品一区二区久久| 大片免费播放器 马上看| 欧美3d第一页| av在线老鸭窝| 九草在线视频观看| 亚洲精品日韩在线中文字幕| 国产精品免费大片| 精品一区二区三区视频在线| 涩涩av久久男人的天堂| 大香蕉久久成人网| 一级毛片我不卡| 日本欧美视频一区| 视频在线观看一区二区三区| 校园人妻丝袜中文字幕| 国产日韩欧美在线精品| 午夜av观看不卡| 丝袜美足系列| 精品人妻在线不人妻| 街头女战士在线观看网站| 少妇猛男粗大的猛烈进出视频| 丝袜喷水一区| 大香蕉久久成人网| 一级毛片我不卡| 在线观看免费视频网站a站| 成人手机av| 亚洲婷婷狠狠爱综合网| 国产av一区二区精品久久| av播播在线观看一区| 日日摸夜夜添夜夜爱| 国产亚洲精品第一综合不卡 | 多毛熟女@视频| 婷婷色av中文字幕| 成人毛片60女人毛片免费| 久久午夜综合久久蜜桃| 久久久国产欧美日韩av| 亚洲色图 男人天堂 中文字幕 | 国产一区二区三区av在线| 婷婷色麻豆天堂久久| 免费播放大片免费观看视频在线观看| 天堂8中文在线网| 人妻少妇偷人精品九色| 国产成人精品一,二区| 成人亚洲欧美一区二区av| 日韩精品有码人妻一区| 国产高清有码在线观看视频| 亚洲av不卡在线观看| 九草在线视频观看| 亚洲天堂av无毛| 在线观看美女被高潮喷水网站| 国产视频内射| 成人午夜精彩视频在线观看| 精品人妻偷拍中文字幕| 有码 亚洲区| a 毛片基地| 亚洲三级黄色毛片| 亚洲欧美精品自产自拍| 久久精品国产自在天天线| 欧美日韩视频高清一区二区三区二| 国产日韩欧美视频二区| 日韩在线高清观看一区二区三区| 国产永久视频网站| a级毛片免费高清观看在线播放| 中文乱码字字幕精品一区二区三区| 国产日韩欧美亚洲二区| 在线播放无遮挡| av在线观看视频网站免费| 少妇的逼好多水| 卡戴珊不雅视频在线播放| 国产69精品久久久久777片| 天天影视国产精品| 久久亚洲国产成人精品v| 国产精品一区二区三区四区免费观看| 国产成人精品久久久久久| 国产一区二区在线观看av| 久久国产亚洲av麻豆专区| 亚洲精品456在线播放app| 在线播放无遮挡| 久久久午夜欧美精品| 欧美国产精品一级二级三级| 精品人妻熟女av久视频| h视频一区二区三区| 欧美精品一区二区大全| 少妇的逼好多水| 国产精品一区www在线观看| 精品视频人人做人人爽| 永久免费av网站大全| 大码成人一级视频| 在线 av 中文字幕| 欧美一级a爱片免费观看看| 午夜福利在线观看免费完整高清在| 一级片'在线观看视频| 国产片特级美女逼逼视频| 亚洲在久久综合| 亚洲精品色激情综合| 国产69精品久久久久777片| 国产精品久久久久久久电影| 久久97久久精品| 国产日韩欧美亚洲二区| 男女啪啪激烈高潮av片| 国产精品久久久久成人av| www.色视频.com| 少妇 在线观看| 亚洲人成77777在线视频| 丰满乱子伦码专区| 中文字幕久久专区| 久久久久久伊人网av| 韩国高清视频一区二区三区| 久久免费观看电影| 在现免费观看毛片| videosex国产| av国产精品久久久久影院| 久久婷婷青草| av又黄又爽大尺度在线免费看| 久久精品久久久久久久性| 免费观看av网站的网址| 18禁在线无遮挡免费观看视频| 国产又色又爽无遮挡免| 久久久国产一区二区| 自线自在国产av| 国产极品天堂在线| 亚洲国产精品成人久久小说| 国产精品久久久久久久电影| 国产综合精华液| 午夜影院在线不卡| 日日撸夜夜添| 欧美日韩亚洲高清精品| 蜜臀久久99精品久久宅男| 我要看黄色一级片免费的| 伊人久久精品亚洲午夜| 美女xxoo啪啪120秒动态图| 天天影视国产精品| 色吧在线观看| 国产精品蜜桃在线观看| 在线观看免费日韩欧美大片 | av网站免费在线观看视频| 精品少妇内射三级| 日韩,欧美,国产一区二区三区| 亚洲欧洲日产国产| 天天影视国产精品| 国产精品国产三级国产av玫瑰| 亚洲精品中文字幕在线视频| 亚洲精品,欧美精品| 久久久久网色| 美女脱内裤让男人舔精品视频| 久久人人爽av亚洲精品天堂| 免费人成在线观看视频色| 精品久久久久久电影网| 久久久久久久大尺度免费视频| 亚洲激情五月婷婷啪啪| 色婷婷久久久亚洲欧美| 日本免费在线观看一区| 99九九在线精品视频| 国产精品人妻久久久久久| 国产 一区精品| 精品久久久久久久久av| 免费高清在线观看视频在线观看| 青青草视频在线视频观看| 中文字幕制服av| 99精国产麻豆久久婷婷| 国产精品欧美亚洲77777| 妹子高潮喷水视频| 久久狼人影院| 国产精品久久久久久久电影| 国产视频首页在线观看| 国产精品久久久久久精品古装| 91精品一卡2卡3卡4卡| 2021少妇久久久久久久久久久| 免费看av在线观看网站| 伦理电影大哥的女人| 国产精品欧美亚洲77777| 99热这里只有是精品在线观看| 国产男女超爽视频在线观看| 有码 亚洲区| 狂野欧美白嫩少妇大欣赏| 一个人看视频在线观看www免费| xxx大片免费视频| 婷婷成人精品国产| 久久人妻熟女aⅴ| 男女国产视频网站| 99热6这里只有精品| 成年女人在线观看亚洲视频| 亚洲综合精品二区| 亚洲精品自拍成人| 丝袜美足系列| 赤兔流量卡办理| 亚洲美女视频黄频| 一级毛片 在线播放| 亚洲精品国产av成人精品| 欧美激情 高清一区二区三区| 观看av在线不卡| 最后的刺客免费高清国语| 欧美变态另类bdsm刘玥| 少妇人妻久久综合中文| a级毛色黄片| 狠狠婷婷综合久久久久久88av| 夫妻午夜视频| 嘟嘟电影网在线观看| 伊人亚洲综合成人网| 观看美女的网站| 久久精品熟女亚洲av麻豆精品| 精品一区二区免费观看| 成人二区视频| 成人毛片a级毛片在线播放| 人人澡人人妻人| 欧美3d第一页| 99视频精品全部免费 在线| 黄片播放在线免费| 国产免费又黄又爽又色| 在线观看免费日韩欧美大片 | 国产精品一二三区在线看| 99视频精品全部免费 在线| 97在线视频观看| 午夜精品国产一区二区电影| 18在线观看网站| 亚洲综合精品二区| 国产精品一区二区三区四区免费观看| 国产男人的电影天堂91| 性色av一级| 成人毛片60女人毛片免费| 一级毛片我不卡| 欧美人与善性xxx| 久久久久久久久大av| 亚洲人与动物交配视频| 免费观看在线日韩| 人人妻人人添人人爽欧美一区卜| 亚洲国产成人一精品久久久| 欧美日韩综合久久久久久| 亚洲精品国产色婷婷电影| 国产无遮挡羞羞视频在线观看| 少妇高潮的动态图| 国产成人精品婷婷| 国产精品久久久久久精品古装| 国产精品久久久久久久电影| 美女xxoo啪啪120秒动态图| 欧美成人午夜免费资源| 久久久亚洲精品成人影院| 在线观看免费高清a一片| 久久婷婷青草| 人人妻人人澡人人看| 午夜激情福利司机影院| 久久久久久久久久人人人人人人| 男女高潮啪啪啪动态图| 国产在线一区二区三区精| 汤姆久久久久久久影院中文字幕| 91精品国产国语对白视频| 国产精品人妻久久久久久| 色吧在线观看| 如何舔出高潮| 一级,二级,三级黄色视频| 亚洲成人手机| 亚洲成人av在线免费| 一区二区三区免费毛片| 精品熟女少妇av免费看| 亚洲国产色片| freevideosex欧美| 啦啦啦啦在线视频资源| 69精品国产乱码久久久| 欧美成人午夜免费资源| 十分钟在线观看高清视频www| 涩涩av久久男人的天堂| 亚洲国产欧美在线一区| 考比视频在线观看| 久久影院123| 亚洲,欧美,日韩| 日本av手机在线免费观看| 亚洲久久久国产精品| 黄色怎么调成土黄色| 国产探花极品一区二区| 卡戴珊不雅视频在线播放| 人人妻人人澡人人爽人人夜夜| 精品国产国语对白av| 九九在线视频观看精品| 美女国产高潮福利片在线看| 欧美精品人与动牲交sv欧美| 国产精品一二三区在线看| 日本av免费视频播放| 婷婷成人精品国产| 亚洲国产av新网站| 国产色爽女视频免费观看| 在线观看三级黄色| 人妻 亚洲 视频| 搡老乐熟女国产| 欧美激情国产日韩精品一区| 国产精品久久久久久av不卡| 免费观看在线日韩| 亚洲人与动物交配视频| 超色免费av| 999精品在线视频| 欧美人与性动交α欧美精品济南到 | 熟妇人妻不卡中文字幕| 极品人妻少妇av视频| 在线观看人妻少妇| √禁漫天堂资源中文www| 国产成人午夜福利电影在线观看| 人人澡人人妻人| 国产欧美亚洲国产| 18禁在线无遮挡免费观看视频| 亚洲经典国产精华液单| 一本一本综合久久| 国产成人精品无人区| 国产黄色视频一区二区在线观看| 国产一区二区三区av在线| 香蕉精品网在线| 国产精品秋霞免费鲁丝片| 日本91视频免费播放| 啦啦啦在线观看免费高清www| 人妻 亚洲 视频| 欧美性感艳星| 亚洲经典国产精华液单| 国产精品熟女久久久久浪| 天天躁夜夜躁狠狠久久av| 午夜激情福利司机影院| 最黄视频免费看| 亚洲国产精品国产精品| 久久精品国产自在天天线| 熟妇人妻不卡中文字幕| 桃花免费在线播放| 日本91视频免费播放| 亚洲av福利一区| 精品亚洲乱码少妇综合久久| 少妇的逼好多水| 欧美激情 高清一区二区三区| av在线播放精品| 成人免费观看视频高清| 大香蕉97超碰在线| 超色免费av| 午夜精品国产一区二区电影| 日韩不卡一区二区三区视频在线| 欧美一级a爱片免费观看看| av一本久久久久| 伦精品一区二区三区| 午夜福利视频精品| 我的老师免费观看完整版| 日韩一区二区三区影片| 亚洲美女视频黄频| 大又大粗又爽又黄少妇毛片口| 一级毛片电影观看| 伊人久久精品亚洲午夜| 最近手机中文字幕大全| 亚洲精品自拍成人| 狠狠精品人妻久久久久久综合| 女人精品久久久久毛片| 80岁老熟妇乱子伦牲交| 午夜免费鲁丝| 亚洲国产日韩一区二区| 最近手机中文字幕大全| 2022亚洲国产成人精品| 欧美成人午夜免费资源| 精品久久久久久久久av| 国产在线视频一区二区| 我要看黄色一级片免费的| 亚洲精品自拍成人| 免费黄网站久久成人精品| 精品午夜福利在线看| 天天影视国产精品| 久久ye,这里只有精品| 日本黄色日本黄色录像| 插阴视频在线观看视频| 亚洲av欧美aⅴ国产| 亚洲精品第二区| 看非洲黑人一级黄片| 老女人水多毛片| 97在线视频观看| 人妻制服诱惑在线中文字幕| 婷婷色综合www| 亚洲精品乱码久久久久久按摩| 晚上一个人看的免费电影| 一本色道久久久久久精品综合| 久久毛片免费看一区二区三区| 中文字幕av电影在线播放| 日本免费在线观看一区| 青春草国产在线视频| 免费播放大片免费观看视频在线观看| 男人添女人高潮全过程视频| 亚洲精品久久午夜乱码| 中文字幕人妻丝袜制服| 国产视频内射| 女人精品久久久久毛片| 亚洲精品一二三| 亚洲婷婷狠狠爱综合网| 男女啪啪激烈高潮av片| a级毛色黄片| 亚洲伊人久久精品综合| 一级爰片在线观看| 人人妻人人添人人爽欧美一区卜| 日本免费在线观看一区| 性色av一级| 一级毛片 在线播放| 18+在线观看网站| 嫩草影院入口| 亚洲精品自拍成人| 伊人久久国产一区二区| 成人影院久久| 国产精品.久久久| 日韩亚洲欧美综合| 国产精品成人在线| av在线老鸭窝| 观看美女的网站| 日本-黄色视频高清免费观看| 国产亚洲一区二区精品| 一本色道久久久久久精品综合| 91精品一卡2卡3卡4卡| 精品国产国语对白av| 视频在线观看一区二区三区| 亚洲av在线观看美女高潮| a级毛片免费高清观看在线播放| 午夜免费观看性视频| 少妇猛男粗大的猛烈进出视频| 成人亚洲精品一区在线观看| 97精品久久久久久久久久精品| 亚洲美女视频黄频| 国产精品国产av在线观看| 免费观看无遮挡的男女| 久久热精品热| 丝瓜视频免费看黄片| 国产片内射在线| 国产男女超爽视频在线观看| 久久精品国产亚洲网站| 成人国产av品久久久| 亚洲内射少妇av| 亚洲精品美女久久av网站| 成年人午夜在线观看视频| 日本色播在线视频| 男女无遮挡免费网站观看| 高清不卡的av网站| 日韩人妻高清精品专区| 免费黄网站久久成人精品| 欧美另类一区| 欧美最新免费一区二区三区| 麻豆乱淫一区二区| 亚洲av中文av极速乱| 精品人妻一区二区三区麻豆| 9色porny在线观看| 久久久国产精品麻豆| 简卡轻食公司| 亚洲怡红院男人天堂| 日韩在线高清观看一区二区三区| 亚洲综合精品二区| 午夜av观看不卡| 青春草视频在线免费观看| 91精品国产国语对白视频| 日韩av在线免费看完整版不卡| 久久精品国产a三级三级三级| 搡老乐熟女国产| av.在线天堂| 晚上一个人看的免费电影| 日韩av免费高清视频| 性高湖久久久久久久久免费观看| 狂野欧美激情性bbbbbb| 亚洲综合色惰| 亚洲无线观看免费| 五月天丁香电影| 美女主播在线视频| 99久久精品一区二区三区| 丝袜美足系列| 国产精品成人在线| 中文字幕久久专区| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| 老司机亚洲免费影院| 久久久久精品久久久久真实原创| 中文字幕人妻熟人妻熟丝袜美| 日本猛色少妇xxxxx猛交久久| 久久久久国产网址| 国产精品一区www在线观看| 国产精品偷伦视频观看了| 亚洲第一av免费看| 久久韩国三级中文字幕| 亚洲一区二区三区欧美精品| 少妇精品久久久久久久| 国产成人免费无遮挡视频| 18在线观看网站| 亚洲第一av免费看| 热re99久久国产66热| 久久99蜜桃精品久久| 免费人成在线观看视频色| av国产久精品久网站免费入址| 国产精品无大码| 99久久精品一区二区三区| 性色avwww在线观看| 观看av在线不卡| 久久久欧美国产精品| 亚洲综合精品二区| av福利片在线| 午夜免费男女啪啪视频观看| 大香蕉97超碰在线| 日韩不卡一区二区三区视频在线| 尾随美女入室| 下体分泌物呈黄色| 99久久中文字幕三级久久日本| 国产日韩欧美亚洲二区| 黄片播放在线免费| 亚洲欧美精品自产自拍| 久久99热6这里只有精品|