CN103777197A - 单基地mimo雷达中降维传播算子的方位估计方法 - Google Patents
单基地mimo雷达中降维传播算子的方位估计方法 Download PDFInfo
- Publication number
- CN103777197A CN103777197A CN201310717502.1A CN201310717502A CN103777197A CN 103777197 A CN103777197 A CN 103777197A CN 201310717502 A CN201310717502 A CN 201310717502A CN 103777197 A CN103777197 A CN 103777197A
- Authority
- CN
- China
- Prior art keywords
- propagation operator
- dimensionality reduction
- matrix
- algorithm
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/04—Systems determining presence of a target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种单基地多输入多输出(MIMO)雷达中基于降维传播算子的高效方位估计方法,包括如下几个步骤:对接收信号进行降维变换;由降维后接收信号对传播算子进行估计;再根据传播算子构造两个矩阵,最后利用两个矩阵之间的旋转不变性估计波达方向。本算法比传统传播算子算法的复杂度大大降低,同时其角度估计性能要稍好于传统传播算子算法,在单基地MIMO雷达的角度估计中有很好的效果。
Description
技术领域
本发明涉及空间信号到达方向(DOA)估计的技术领域,特别涉及单基地MIMO雷达中降维传播算子的方位估计方法。
背景技术
MIMO雷达利用多天线技术,同时发射独立的信号波束来探测目标,并接收相应的目标反射信号,它比较传统的相控阵雷达有很多潜在的优势。MIMO雷达系统比其他只有一个传输天线的系统有更多自由度,这些额外的自由度能克服衰落效应,提高空间分辨率,增强参数的可识别性并提高目标探测性能。
我们已经研究过MIMO雷达的角度估计算法,但这些算法都有不同的缺陷。Capon算法和多重信号分类(MUSIC)算法需要高复杂度的谱峰搜索;借助旋转不变技术估计信号参数(ESPRIT)算法和求根MUSIC算法需要对互相关矩阵进行特征值分解来获得信号子空间或噪声子空间。平行因子分析(PARAFAC)算法和自适应PARAFAC算法是迭代算法,这通常比子空间算法需要更高的复杂度。由于传播算子算法不需要对互相关矩阵进行特征值分解,也不需要对接收的数据进行奇异值分解,因此它的计算量可以显著减小。
发明内容
本发明目的是针对现有算法存在的缺陷,基于传播算子算法和降维变换提出了一种单基地MIMO雷达下的高效波达方向估计算法。该算法比传统传播算子算法的复杂度更低,同时角度估计性能要稍好于传统传播算子算法。这种算法在不进行谱峰搜索时也十分有效。
本发明一种单基地MIMO雷达中基于降维传播算子的波达方向估计方法包括如下几个步骤:
1)对接收信号进行降维变换:
y(t)=W-1GHx(t)
2)然后由降维后接收信号对传播算子进行估计。对降维后的接收信号进行L次采样并定义Y=[y(1),y(2),…,y(L)],它可由Y=[Y1 T,Y2 T]T表示,其中 进而我们可以得到P的最小二乘估计,
本发明的优点在于:①充分利用了单基地MIMO雷达收发天线均匀线阵特性;②该算法角度估计性能好于传统传播算子算法;③该算法比传统传播算子算法的复杂度更低,无需谱峰搜索。
附图说明
图1是几种算法在M=8,N=8,K=3,且L取不同值时的复杂度对比;
图2是本算法在M=8,N=8,L=100和SNR=10dB时对所有三个目标的角度估计结果;
图3是几种算法在M=8,N=8,L=100时的角度估计性能对比;
图4是本算法在M=8,N=8,且L取不同值时的角度估计性能对比;
图5是本算法在N=8,M取不同值时的角度估计性能比较;
图6是本算法在M=8,N取不同值时的角度估计性能比较;
图7是几种算法在M=8,N=8,L=100时,对两个很近的目标的角度估计性能对比。
具体实施方式
考虑一个包含M个发射天线和N个接收天线的单基地MIMO雷达系统,并且发射和接收阵列都是均匀线阵,天线间隔均为半波长。我们假设有K个不相关的目标,接收端匹配滤波器的输出可以表示为:
其中θk是第k个目标的波达方向;且 fk是多普勒频率,βk是幅度。ar(θk)和at(θk)分别是θk的接收方向矢量和发射方向矢量;是它们的Kronecker算子。n(t)是一个零均值协方差矩阵为σ2IMN的MN×1高斯白噪声向量。我们定义 进行L次采样并定义X=[x(1),x(2),…,x(L)],它可由下式表示
其中AT=[at(θ1),at(θ2),…,at(θK)]∈CM×K是发射方向矩阵;N∈CMN×L是噪声矩阵;S=[s(1),s(2),…,s(L)]T∈CL×K;Φ=diag(exp(-jπsinθ1),exp(-jπsinθ2),…,exp(-jπsinθK))。
其中b(θk)=[1,exp(-jπsinθk),…,exp(-jπ(M+N-2)sinθk)]T,并且矩阵G为:
为接收信号x(t)使用降维变换W-1GH,我们得到:
y(t)=W-1GHx(t)
=W-1W[b(θ1),b(θ2),…,b(θK)]s(t)+W-1GHn(t) (6)
=Bs(t)+W-1GHn(t)
其中B1∈CK×K,B2∈C(M+N-1-K)×K。B1是满秩矩阵。定义B2=PB1,其中矩阵P是传播算子。式(7)中的矩阵B2可以被重新写为
其中p1∈C1×K和pend∈C1×K分别是P的第一行和最后一行。根据式(8-9),我们有
B21=P1B1,B22=P2B1 (10)
算法的几个主要步骤如下:
(1)对接收信号进行降维变换,y(t)=W-1GHx(t);
下面给出仿真结果:
我们提出Monte Carlo仿真来评价本算法的角度估计性能。定义均方根误差(RMSE)为其中是第n次Monte Carlo试验的波达方向θk的估计值。注意:M,N,L和K分别是发射天线,接收天线,快拍和目标的数量。在接下来除了图7的仿真中,我们采用有三个位于θ1=10°,θ2=20°和θ3=30°的非相关目标的MIMO雷达系统。
从图1中我们发现本算法比PM算法的计算量小很多。
图3表明本算法的角度估计性能要稍好于PM算法。
图4表明本算法的角度估计性能随L的增加而提高。
由图5-6可知,本算法的角度估计性能随接收和发射天线数量的增加而提高。由于分集增益,更多的接收和发射天线提高了角度估计性能。
从图7中我们发现在目标相近时本算法比传统传播算子算法在角度估计上要稍好一些。
Claims (3)
1.一种单基地多输入多输出(MIMO)雷达中基于降维传播算子的高效波达方向估计方法,其特征在于包括如下几个步骤:
(1)对接收信号进行降维变换:
y(t)=W-1GHx(t)
其中:x(t)接收端匹配滤波器的输出信号,
(2)然后由降维后接收信号对传播算子进行估计。对降维后的接收信号进行L次采样并定义Y=[y(1),y(2),…,y(L)],它可由Y=[Y1 T,Y2 T]T表示,其中 进而我们可以得到P的最小二乘估计,
2.根据权利要求1所述的一种单基地MIMO雷达中基于降维传播算子的高效波达方向估计方法,其特征在于降维变换前后噪声方差不变,但复杂度大大降低。
3.根据权利要求1所述的一种单基地MIMO雷达中基于降维传播算子的高效波达方向估计方法,其特征在于将降维信号进行分块,再利用最小二乘求传播算子,此方法无需协方差矩阵的求取和其特征值分解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310717502.1A CN103777197A (zh) | 2013-12-24 | 2013-12-24 | 单基地mimo雷达中降维传播算子的方位估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310717502.1A CN103777197A (zh) | 2013-12-24 | 2013-12-24 | 单基地mimo雷达中降维传播算子的方位估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103777197A true CN103777197A (zh) | 2014-05-07 |
Family
ID=50569663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310717502.1A Pending CN103777197A (zh) | 2013-12-24 | 2013-12-24 | 单基地mimo雷达中降维传播算子的方位估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103777197A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104991236A (zh) * | 2015-06-18 | 2015-10-21 | 哈尔滨工程大学 | 一种单基地mimo雷达非圆信号相干源波达方向估计方法 |
CN105490763A (zh) * | 2015-12-04 | 2016-04-13 | 南京航空航天大学 | 一种端到端的宽带移动mimo传播信道模型及建模方法 |
CN106959443A (zh) * | 2017-02-27 | 2017-07-18 | 中国人民解放军空军工程大学 | 一种基于空间差异矩阵重构的低角目标doa估计方法及装置 |
CN107064926A (zh) * | 2017-01-18 | 2017-08-18 | 长江大学 | 空域色噪声背景下的双基地mimo雷达角度估算方法 |
CN107884741A (zh) * | 2017-10-30 | 2018-04-06 | 北京理工大学 | 一种多球阵列多宽带声源快速定向方法 |
CN108344967A (zh) * | 2018-01-20 | 2018-07-31 | 中国人民解放军战略支援部队信息工程大学 | 基于互质面阵的二维波达方向快速估计方法 |
CN110007282A (zh) * | 2019-03-15 | 2019-07-12 | 南京航空航天大学 | 一种连续波体制1-比特雷达目标重构问题降维方法 |
CN111131116A (zh) * | 2019-12-12 | 2020-05-08 | 国网江苏省电力有限公司信息通信分公司 | 一种频率偏移估计方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0535780A1 (en) * | 1991-09-30 | 1993-04-07 | Trw Inc. | A compact, flexible, and integrated millimeter wave radar sensor |
WO1994011751A1 (en) * | 1992-11-17 | 1994-05-26 | Honeywell Inc. | Real time passive threat positioning system |
CN101566683A (zh) * | 2009-03-24 | 2009-10-28 | 西北工业大学 | 基于相位差波束形成的目标方位估计方法 |
CN101799535A (zh) * | 2009-11-27 | 2010-08-11 | 西安电子科技大学 | Mimo雷达目标方向的估计方法 |
CN102169170A (zh) * | 2010-12-29 | 2011-08-31 | 电子科技大学 | 一种相干分布式信号二维波达角的测定方法 |
US20120274499A1 (en) * | 2011-04-29 | 2012-11-01 | Spatial Digital Systems | Radar imaging via spatial spectrum measurement and MIMO waveforms |
CN103091671A (zh) * | 2013-01-15 | 2013-05-08 | 西安电子科技大学 | 基于非同心电磁矢量阵列雷达的两维波达方向估计方法 |
-
2013
- 2013-12-24 CN CN201310717502.1A patent/CN103777197A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0535780A1 (en) * | 1991-09-30 | 1993-04-07 | Trw Inc. | A compact, flexible, and integrated millimeter wave radar sensor |
WO1994011751A1 (en) * | 1992-11-17 | 1994-05-26 | Honeywell Inc. | Real time passive threat positioning system |
CN101566683A (zh) * | 2009-03-24 | 2009-10-28 | 西北工业大学 | 基于相位差波束形成的目标方位估计方法 |
CN101799535A (zh) * | 2009-11-27 | 2010-08-11 | 西安电子科技大学 | Mimo雷达目标方向的估计方法 |
CN102169170A (zh) * | 2010-12-29 | 2011-08-31 | 电子科技大学 | 一种相干分布式信号二维波达角的测定方法 |
US20120274499A1 (en) * | 2011-04-29 | 2012-11-01 | Spatial Digital Systems | Radar imaging via spatial spectrum measurement and MIMO waveforms |
CN103091671A (zh) * | 2013-01-15 | 2013-05-08 | 西安电子科技大学 | 基于非同心电磁矢量阵列雷达的两维波达方向估计方法 |
Non-Patent Citations (2)
Title |
---|
吴海浪等: "十字阵的单基地MIMO雷达中低复杂度的二维DOA估计算法", 《数据采集与处理》 * |
孙中伟等: "双基地MIMO雷达中基于传播算子的DOD和DOA估计算法", 《数据采集与处理》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104991236B (zh) * | 2015-06-18 | 2017-10-31 | 哈尔滨工程大学 | 一种单基地mimo雷达非圆信号相干源波达方向估计方法 |
CN104991236A (zh) * | 2015-06-18 | 2015-10-21 | 哈尔滨工程大学 | 一种单基地mimo雷达非圆信号相干源波达方向估计方法 |
CN105490763A (zh) * | 2015-12-04 | 2016-04-13 | 南京航空航天大学 | 一种端到端的宽带移动mimo传播信道模型及建模方法 |
CN105490763B (zh) * | 2015-12-04 | 2018-06-26 | 南京航空航天大学 | 一种端到端的宽带移动mimo传播信道模型及建模方法 |
CN107064926A (zh) * | 2017-01-18 | 2017-08-18 | 长江大学 | 空域色噪声背景下的双基地mimo雷达角度估算方法 |
CN107064926B (zh) * | 2017-01-18 | 2020-05-19 | 长江大学 | 空域色噪声背景下的双基地mimo雷达角度估算方法 |
CN106959443A (zh) * | 2017-02-27 | 2017-07-18 | 中国人民解放军空军工程大学 | 一种基于空间差异矩阵重构的低角目标doa估计方法及装置 |
CN106959443B (zh) * | 2017-02-27 | 2020-03-27 | 中国人民解放军空军工程大学 | 一种基于空间差异矩阵重构的低角目标doa估计方法及装置 |
CN107884741A (zh) * | 2017-10-30 | 2018-04-06 | 北京理工大学 | 一种多球阵列多宽带声源快速定向方法 |
CN107884741B (zh) * | 2017-10-30 | 2021-01-19 | 北京理工大学 | 一种多球阵列多宽带声源快速定向方法 |
CN108344967A (zh) * | 2018-01-20 | 2018-07-31 | 中国人民解放军战略支援部队信息工程大学 | 基于互质面阵的二维波达方向快速估计方法 |
CN108344967B (zh) * | 2018-01-20 | 2020-05-26 | 中国人民解放军战略支援部队信息工程大学 | 基于互质面阵的二维波达方向快速估计方法 |
CN110007282A (zh) * | 2019-03-15 | 2019-07-12 | 南京航空航天大学 | 一种连续波体制1-比特雷达目标重构问题降维方法 |
CN110007282B (zh) * | 2019-03-15 | 2022-05-13 | 南京航空航天大学 | 一种连续波体制1-比特雷达目标重构问题降维方法 |
CN111131116A (zh) * | 2019-12-12 | 2020-05-08 | 国网江苏省电力有限公司信息通信分公司 | 一种频率偏移估计方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103777197A (zh) | 单基地mimo雷达中降维传播算子的方位估计方法 | |
CN106526530B (zh) | 基于传播算子的2-l型阵列二维doa估计算法 | |
CN104251989B (zh) | 基于压缩空间谱的单基地mimo雷达目标波达方向估计方法 | |
CN107064892B (zh) | 基于张量子空间和旋转不变的mimo雷达角度估计算法 | |
CN107015213B (zh) | 基于music算法的双基地mimo雷达角度估算方法 | |
CN102707264B (zh) | 基于圆形阵列双基地mimo雷达的波达方向估计方法 | |
CN106443610B (zh) | Mimo雷达收发阵列互耦误差自校正方法 | |
Yan et al. | Two-stage ESPRIT for unambiguous angle and range estimation in FDA-MIMO radar | |
CN103983958A (zh) | 基于多测量矢量稀疏表示的mimo雷达连续目标角度估计方法 | |
CN101369014A (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN110133631B (zh) | 一种基于模糊函数的频控阵mimo雷达目标定位方法 | |
CN108562866B (zh) | 基于矩阵填充的双基地mimo雷达角度估算方法 | |
CN101887117B (zh) | 基于三迭代的机载mimo雷达空时降维自适应处理方法 | |
CN108303683B (zh) | 单基地mimo雷达实值esprit非圆信号角度估计方法 | |
CN107064926B (zh) | 空域色噪声背景下的双基地mimo雷达角度估算方法 | |
CN101799535A (zh) | Mimo雷达目标方向的估计方法 | |
CN104833959A (zh) | 基于目标先验信息的mimo雷达波形优化方法 | |
CN105319545A (zh) | 提高stap检测性能的mimo-ofdm雷达波形设计方法 | |
CN104991236A (zh) | 一种单基地mimo雷达非圆信号相干源波达方向估计方法 | |
CN105487054B (zh) | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 | |
CN105182325A (zh) | 基于秩1约束的米波mimo雷达低仰角目标测高方法 | |
CN107132514A (zh) | 改进esprit的双基地mimo雷达角度估算方法 | |
CN103217671B (zh) | 色噪声环境下的多输入多输出雷达收发角度快速估计方法 | |
Li et al. | Parameter estimation based on fractional power spectrum density in bistatic MIMO radar system under impulsive noise environment | |
CN104678362B (zh) | Mimo天波超视距雷达波形优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140507 |