CN107703478A - 基于互相关矩阵的扩展孔径二维doa估计方法 - Google Patents
基于互相关矩阵的扩展孔径二维doa估计方法 Download PDFInfo
- Publication number
- CN107703478A CN107703478A CN201711026670.0A CN201711026670A CN107703478A CN 107703478 A CN107703478 A CN 107703478A CN 201711026670 A CN201711026670 A CN 201711026670A CN 107703478 A CN107703478 A CN 107703478A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- matrix
- formula
- estimation
- 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.)
- Granted
Links
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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及采用阵列天线估计接收信号到达方向估计技术,为只需对互相关矩阵子块进行线性运算,具有更低的计算复杂度;且实现俯仰角与方位角的自动配对。为此,本发明采用的技术方案是,基于互相关矩阵的扩展孔径二维DOA估计方法,步骤如下:1)构造互相关矩阵;2)划分互相关矩阵得到阶数为信号个数的子块;3)根据子块构建包含自动配对的方位角和俯仰角信息矩阵;4)得到z轴低精度无模糊的方向余弦估计以及对应的特征向量矩阵γ;5)得到z轴高精度模糊的方向余弦;6)重复步骤1‑5,可以得到对应x,y轴上对应的两种方向余弦;7)得到高精度无模糊的方向余弦8)得到自动配对本发明主要应用于阵列天线估计接收信号到达方向估计。
Description
技术领域
本发明涉及采用阵列天线估计接收信号到达方向的技术领域,尤其涉及采用非均匀2-L型天线阵列的信号到达方向估计方法,具体讲,涉及基于互相关矩阵的扩展孔径二维DOA估计方法。
背景技术
空间信号到达方向(Direction of Arrival,DOA)估计是空间谱估计一个主要研究方向,被广泛应用在雷达、声呐、地震、通信等许多领域。DOA估计的基本问题就是确定各个信号到达阵列参考阵元的方向角,简称波达方向。经典的子空间分解类DOA估计算法有多重信号分类算法(MUSIC,Multiple Signal Classification)和基于旋转不变技术的信号参数估计算法(ESPRIT,Estimation of Signal Parameter via Rotational InvitationTechniques)。其中MUSIC算法是噪声子空间类算法,ESPRIT算法是信号子空间类算法,改进的MUSIC算法包括特征矢量法、求根MUSIC法、加权MUSIC算法等,改进的ESPRIT算法包括最小二乘ESPRIT、总体最小二乘ESPRIT、加权ESPRIT算法等。
传统的MUSIC算法和ESPRIT算法等高分辨率算法,虽然具有良好的估计性能,但是由于需要对接收信号协方差矩阵进行特征值分解,因此具有较大的计算量。传播算子算法使用线性运算使用线性运算代替了奇异值分解和特征值分解运算,显著地降低了计算复杂度。由于传播算子算法具有计算复杂度较低的优点,各国学者们对其进行了广泛研究,并提出大量基于传播算子的DOA估计算法。目前,存在大量基于传播算子的L型阵列、2-L型阵列、双平行线阵、三平行线阵等二维DOA估计算法。但是某些基于双平行线阵的传播算子算法在俯仰角为70°~90°的实际移动通信俯仰角度范围内存在角度估计失效问题,有些基于三平行线阵采用传播算子的二维DOA估计算法并没有充分利用所有的阵元信息,有些基于2-L型阵列的采用传播算子的二维DOA估计算法分别利用阵列的两个L型子阵,单独估计信号的方位角和俯仰角,估计性能较差。与基于传播算子的算法相比,基于互相关矩阵的算法能够消除高斯白噪声的影响,提高角度估计性能。另一方面,扩展孔径可以有效地提高阵列的分辨率和角度估计精度,但会出现模糊的现象。有些算法提出了解模糊的算法,但是由于算法本身是基于ESPRIT算法,因此具有较大的计算复杂度。
发明内容
为克服现有技术的不足,扩展孔径以及基于互相关矩阵的算法均能有效地提高DOA估计性能。因此,本发明旨在提出一种基于互相关矩阵的扩展孔径二维DOA估计算法。与已提出的扩展孔径DOA估计算法相比,本文提出的算法构造更低维度的互相关矩阵,只需对互相关矩阵子块进行线性运算,具有更低的计算复杂度;且实现俯仰角与方位角的自动配对。为此,本发明采用的技术方案是,基于互相关矩阵的扩展孔径二维DOA估计方法,步骤如下:
1)构造互相关矩阵;
2)划分互相关矩阵得到阶数为信号个数的子块;
3)根据子块构建包含自动配对的方位角和俯仰角信息矩阵;
4)得到z轴低精度无模糊的方向余弦估计以及对应的特征向量矩阵γ;
5)得到z轴高精度模糊的方向余弦;
6)重复步骤1-5,可以得到对应x,y轴上对应的两种方向余弦;
7)得到高精度无模糊的方向余弦
8)得到自动配对
互相关矩阵的构造及划分步骤:
选取子阵W构造了互相关矩阵Rzw,Rxw,Rpw,Ryw,Rqw,然后,选取了子阵Z,P构造了Rzp;
其中,Rzw具体构造过程如下
Rzw=E[z(t)wH(t)]
=E[(Azs(t)+nz(t))(Aws(t)+nw(t))H]
=E[Azs(t)(Aws(t))H]+E[Azs(t)nw(t)H]
+E[nz(t)(Aws(t))H]+E[nz(t)nw(t)H] (6)
Rxw,Rpw,Ryw,Rqw,Rzp按照相同的方式构造;
由于式(6)中噪声信息互不相关且与s(t)相互独立,则式(6)简化为
Rzw=AzRsAw H (7)
同理,Rxw,Rpw,Ryw,Rqw,Rzp简化为:
Rxw=AxRsAw H
Rpw=ApRsAw H
Ryw=AyRsAw H
Rqw=AqRsAw H
Rzp=AzRsAp H (8)
式中Rs(t)=E[s(t)s(t)H]表示s(t)的自相关矩阵,由式(8)知,互相关矩阵Rzw,Rxw,Rpw,Ryw,Rqw,Rzp均不包含噪声信息;
由式(7)、式(8)知,互相关矩阵的维数均为M×M,为了实现K个信号的角度估计,矩阵的阶数必须为K。Ax,Ay,Az,Ap,Aq,Aw分别按如下形式分块
式中分别为维数K×K,(M-K)×K的矩阵;
与式(9)相对应,将Rzw按如下形式进行分块
Rxw,Rpw,Ryw,Rqw,Rzp按照相同的方式划分。其中,zw21,xw21,pw21,yw21,qw21均为维数(M-K)×K矩阵,pw12为维数K×(M-K)矩阵,zw11,zp11均为维数K×K矩阵,得到阶数为K的子块用于方向余弦的估计。
方向余弦的估计步骤:
为了从子块中得到方位角和俯仰角信息以及实现俯仰角与方位角的自动配对,根据互相关矩阵部分信息构建如下矩阵:
ψxz=xw21zw11 -1=Ax2Az1 -1
ψpz=pw21zw11 -1=Ap2Az1 -1
ψyz=yw21zw11 -1=Ay2Az1 -1
ψqz=qw21zw11 -1=Aq2Az1 -1
ψzz=zw21zw11 -1=Az2Az1 -1 (11)
与以上定义不同的是ψwz=pw12 Hzp11 -1=Aw2Az1 -1。ψzz,ψwz,ψxz,ψpz,ψyz,ψqz均为维数(M-K)×K矩阵,
根据式(11)可得
式(12)直观地给出了算法实现俯仰角和俯仰角自动配对的依据;
由式(4)和式(11)可知,ψzz,ψwz之间包含无模糊的方向余弦估计信息,则ψzz,ψwz之间的关系用如下公式表示
通过对进行特征值分解,得到与Φz对应的特征值以及特征向量矩阵γ,因为则z轴低精度无模糊的方向余弦估计为:
ψzz,ψwz自身分别包含高精度模糊的方向余弦估计信息,因此,ψzz、ψwz分别取前M-K-1行构成维数为2(M-K-1)×K矩阵P,ψzz、ψwz后M-K-1行构成维数为2(M-K-1)×K矩阵Q,则P,Q用如下公式表示:
P=Λ1Az1 -1
Q=Λ2Az1 -1 (15)
式中,Λ1、Λ2分别表示由Az2、Aw2的前M-K-1行和后M-K-1行构成的维数为2(M-K-1)×K的矩阵;
由式(12)可知,定义B1=Pγ,B2=Qγ,即实现z轴方向上低精度无模糊的方向余弦估计与高精度模糊的方向余弦估计的一一对应,则Φz的估计值可通过B1,B2求得,即因为方向余弦的范围{-1≤υ≤1},通过对进行特征值分解,则z轴高精度模糊的方向余弦估计为:
式中,表示不小于i的最小整数,表示不大于i的最大整数,表示的第k个对角元素;
同理,可得到x轴y轴低精度无模糊的方向余弦估计以及对应高精度模糊的方向余弦估计
实现角度估计的具体步骤:
利用解模糊的方法,z轴高精度无模糊的方向余弦估计为
其中,用下式进行估计
同理,x轴,y轴高精度无模糊的方向余弦估计分别为
根据以上的分析,第k个信号的二维DOA估计表达式如下:
本发明的特点及有益效果是:
通过构造6个维度较小矩阵,能够以较低的计算复杂度获得较好的方位角和俯仰角估计性能;扩展孔径也提高了角度估计性能;能够实现方位角和俯仰角的自动配对;在俯仰角为70°~90°的实际移动通信的俯仰角度范围内不会出现角度模糊。
附图说明:
图1天线阵列结构示意图。
图2方位角估计直方图。
图3俯仰角估计直方图。
图4不同角度组合估计联合均方误差图。
图5方位角估计均方误差随信噪比变化情况。
图6俯仰角估计均方误差随信噪比变化情况。
具体实施方式
针对已有DOA估计算法存在的问题,本发明提出了一种基于非均匀2-L型阵列的扩展二维DOA估计算法,该天线阵列为非均匀2-L型阵列,其中在x轴、y轴和z轴上分别有两个阵元数目为M的均匀线阵,分别用X,Y,Z,P,Q,W表示。各子阵中阵元间距为来波信号波长的一半;各坐标轴上的两个子阵阵元间距为ds(ds=tλ,t为正整数)。
本发明采用的技术方案:基于互相关矩阵的扩展孔径二维DOA估计算法,其特征在于,包括以下步骤:
互相关矩阵的构造及划分
1)构造互相关矩阵;2)划分互相关矩阵得到阶数为信号个数的子块。
方向余弦的估计
3)根据子块构建包含自动配对的方位角和俯仰角信息矩阵;4)得到z轴低精度无模糊的方向余弦估计以及对应的特征向量矩阵γ。
角度估计的实现
5)得到z轴高精度模糊的方向余弦;6)重复步骤1-5,可以得到对应x,y轴上对应的两种方向余弦;7)得到高精度无模糊的方向余弦8)得到自动配对
下面将结合附图和实施例对本发明做进一步的描述:
构造如图1所示的非均匀2-L型天线阵列。假设空间中有K个窄带非相关信号入射到阵列上,其中第k个信号的二维波达方向为和θk分别为来波信号的方位角和俯仰角。
1互相关矩阵的构造及划分
我们选取子阵W构造了互相关矩阵Rzw,Rxw,Rpw,Ryw,Rqw。然后,选取了子阵Z,P构造了Rzp。
其中,Rzw具体构造过程如下
Rzw=E[z(t)wH(t)]
=E[(Azs(t)+nz(t))(Aws(t)+nw(t))H]
=E[Azs(t)(Aws(t))H]+E[Azs(t)nw(t)H]
+E[nz(t)(Aws(t))H]+E[nz(t)nw(t)H] (6)
Rxw,Rpw,Ryw,Rqw,Rzp按照相同的方式构造。
由于式(6)中噪声信息互不相关且与s(t)相互独立,则式(6)可简化为
Rzw=AzRsAw H (7)
同理,Rxw,Rpw,Ryw,Rqw,Rzp可简化为:
Rxw=AxRsAw H
Rpw=ApRsAw H
Ryw=AyRsAw H
Rqw=AqRsAw H
Rzp=AzRsAp H (8)
式中Rs(t)=E[s(t)s(t)H]表示s(t)的自相关矩阵。由式(8)可知,互相关矩阵Rzw,Rxw,Rpw,Ryw,Rqw,Rzp均不包含噪声信息。
由式(7)、式(8)可知,互相关矩阵的维数均为M×M。为了实现K个信号的角度估计,矩阵的阶数必须为K。Ax,Ay,Az,Ap,Aq,Aw分别按如下形式分块
式中分别为维数K×K,(M-K)×K的矩阵。
与式(9)相对应,将Rzw按如下形式进行分块
Rxw,Rpw,Ryw,Rqw,Rzp按照相同的方式划分。其中,zw21,xw21,pw21,yw21,qw21均为维数(M-K)×K矩阵,pw12为维数K×(M-K)矩阵,zw11,zp11均为维数K×K矩阵。得到阶数为K的子块用于方向余弦的估计。
2方向余弦的估计
为了从子块中得到方位角和俯仰角信息以及实现俯仰角与方位角的自动配对,根据互相关矩阵部分信息构建如下矩阵:
ψxz=xw21zw11 -1=Ax2Az1 -1
ψpz=pw21zw11-1=Ap2Az1 -1
ψyz=yw21zw11 -1=Ay2Az1 -1
ψqz=qw21zw11 -1=Aq2Az1 -1
ψzz=zw21zw11 -1=Az2Az1 -1 (11)
与以上定义不同的是ψwz=pw12 Hzp11 -1=Aw2Az1 -1。根据相应的互相关矩阵得到相应的信息矩阵ψzz,ψwz,ψxz,ψpz,ψyz,ψqz(下标表示对应的子阵)均为维数(M-K)×K矩阵。
根据式(11)可得
式(12)直观地给出了算法实现俯仰角和俯仰角自动配对的依据。
由式(4)和式(11)可知,ψzz,ψwz之间包含无模糊的方向余弦估计信息,则ψzz,ψwz之间的关系可用如下公式表示
通过对进行特征值分解,可以得到与Φz对应的特征值以及特征向量矩阵γ,因为则z轴低精度无模糊的方向余弦估计为
ψzz,ψwz自身分别包含高精度模糊的方向余弦估计信息。因此,ψzz、ψwz分别取前M-K-1行构成维数为2(M-K-1)×K矩阵P,ψzz、ψwz后M-K-1行构成维数为2(M-K-1)×K矩阵Q,则P,Q可用如下公式表示:
P=Λ1Az1 -1
Q=Λ2Az1 -1 (15)
式中,Λ1、Λ2分别表示由Az2、Aw2的前M-K-1行和后M-K-1行构成的维数为2(M-K-1)×K的矩阵。
由式(12)可知,定义B1=Pγ,B2=Qγ,即可实现z轴方向上低精度无模糊的方向余弦估计与高精度模糊的方向余弦估计的一一对应。则Φz的估计值可通过B1,B2求得,即因为方向余弦的范围{-1≤υ≤1},通过对进行特征值分解,则z轴高精度模糊的方向余弦估计为:
式中,表示不小于i的最小整数,表示不大于i的最大整数,表示的第k个对角元素。
同理,可得到x轴y轴低精度无模糊的方向余弦估计以及对应高精度模糊的方向余弦估计
3角度估计的实现
根据式(12)可知,本文提出算法的方向余弦一一对应,因此只需对nz,nx,ny进行估计即可。
利用解模糊的方法,z轴高精度无模糊的方向余弦估计为
其中,用下式进行估计
同理,x轴,y轴高精度无模糊的方向余弦估计分别为
根据以上的分析,第k个信号的二维DOA估计表达式如下:
结合上述步骤中的实施方式,对本发明的有效性进行仿真验证如下:
仿真中取M=5,即非均匀2-L型阵列共有28个阵元,阵列间距d=0.5λ,ds=5λ,其中λ为信号波长,对于每次仿真实验取快拍数为200,进行M=500次蒙特卡洛仿真。
仿真实验1:假设有K=2个等功率非相关信号入射到天线阵列,其中SNR=10dB,信号的方位角和俯仰角为图2和图3显示了方位角估计直方图和俯仰角估计直方图。从图中可以看出,本文提出的算法能够准确的分辨这两个来波信号。
仿真实验2:假设有K=1个信号入射到天线阵列,SNR=15dB,其中信号的方位角和俯仰角均在12°~86°之间以2°的步长变化。图4为不同角度组合估计联合均方误差图。
仿真实验3:假设有K=2个等功率非相关信号入射到天线阵列,其中SNR在5dB~30dB之间以5dB的步长变化,信号的方位角和俯仰角为(φ1,θ1)=(45°,65°),(φ2,θ2)=(70°,85°)。图5和图6分别为方位角和俯仰角估计均方误差随信噪比的变化情况。从图中可以看出,随着信噪比的增加,方位角和俯仰角均方误差均减小。
Claims (4)
1.一种基于互相关矩阵的扩展孔径二维DOA估计方法,其特征是,步骤如下:
1)构造互相关矩阵;
2)划分互相关矩阵得到阶数为信号个数的子块;
3)根据子块构建包含自动配对的方位角和俯仰角信息矩阵;
4)得到z轴低精度无模糊的方向余弦估计以及对应的特征向量矩阵γ;
5)得到z轴高精度模糊的方向余弦;
6)重复步骤1-5,可以得到对应x,y轴上对应的两种方向余弦;
7)得到高精度无模糊的方向余弦
8)得到自动配对
2.如权利要求1所述的基于互相关矩阵的扩展孔径二维DOA估计方法,其特征是,互相关矩阵的构造及划分步骤:
选取子阵W构造了互相关矩阵Rzw,Rxw,Rpw,Ryw,Rqw,然后,选取了子阵Z,P构造了Rzp;
其中,Rzw具体构造过程如下
Rzw=E[z(t)wH(T)]
=E[(Azs(t)+nz(T))(Aws(t)+nw(T))H]
=E[Azs(T)(Aws(T))H]+E[Azs(t)nw(T)H]
+E[nz(T)(Aws(T))H]+E[nz(T)nw(T)H] (6)
Rxw,Rpw,Ryw,Rqw,Rzp按照相同的方式构造;
由于式(6)中噪声信息互不相关且与s(t)相互独立,则式(6)简化为
Rzw=AzRsAw H (7)
同理,Rxw,Rpw,Ryw,Rqw,Rzp简化为:
Rxw=AxRsAw H
Rpw=ApRsAw H
Ryw=AyRsAw H
Rqw=AqRsAw H
Rzp=AzRsAp H (8)
式中Rs(T)=E[s(T)s(T)H]表示s(T)的自相关矩阵,由式(8)知,互相关矩阵Rzw,Rxw,Rpw,Ryw,Rqw,Rzp均不包含噪声信息;
由式(7)、式(8)知,互相关矩阵的维数均为M×M,为了实现K个信号的角度估计,矩阵的阶数必须为K。Ax,Ay,Az,Ap,Aq,Aw分别按如下形式分块
<mrow>
<msub>
<mi>A</mi>
<mi>&eta;</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<msub>
<mi>&eta;</mi>
<mn>1</mn>
</msub>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<msub>
<mi>&eta;</mi>
<mn>2</mn>
</msub>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中分别为维数K×K,(M-K)×K的矩阵;
与式(9)相对应,将Rzw按如下形式进行分块
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>z</mi>
<mi>w</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>zw</mi>
<mn>11</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>zw</mi>
<mn>12</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>zw</mi>
<mn>21</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>zw</mi>
<mn>22</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
Rxw,Rpw,Ryw,Rqw,Rzp按照相同的方式划分。其中,zw21,xw21,pw21,yw21,qw21均为维数(M-K)×K矩阵,pw12为维数K×(M-K)矩阵,zw11,zp11均为维数K×K矩阵,得到阶数为K的子块用于方向余弦的估计。
3.如权利要求1所述的基于互相关矩阵的扩展孔径二维DOA估计方法,其特征是,方向余弦的估计步骤:
为了从子块中得到方位角和俯仰角信息以及实现俯仰角与方位角的自动配对,根据互相关矩阵部分信息构建如下矩阵:
ψxz=xw21zw11 -1=Ax2Az1 -1
ψpz=pw21zw11 -1=Ap2Az1 -1
ψyz=yw21zw11 -1=Ay2Az1 -1
ψqz=qw21zw11 -1=Aq2Az1 -1
ψzz=zw21zw11=Az2Az1 (11)
与以上定义不同的是ψwz=pw12 Hzp11 -1=Aw2Az1 -1。ψzz,ψwz,ψxz,ψpz,ψyz,ψqz均为维数(M-K)×K矩阵,
根据式(11)可得
<mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>z</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>w</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<msubsup>
<mi>A</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>w</mi>
<mi>z</mi>
</mrow>
</msub>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>x</mi>
<mi>z</mi>
</mrow>
</msub>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>p</mi>
<mi>z</mi>
</mrow>
</msub>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>y</mi>
<mi>z</mi>
</mrow>
</msub>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>q</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
式(12)直观地给出了算法实现俯仰角和俯仰角自动配对的依据;
由式(4)和式(11)可知,ψzz,ψwz之间包含无模糊的方向余弦估计信息,则ψzz,ψwz之间的关系用如下公式表示
<mrow>
<msub>
<mover>
<mi>&Phi;</mi>
<mo>^</mo>
</mover>
<mi>z</mi>
</msub>
<mo>=</mo>
<msup>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>z</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>+</mo>
</msup>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>w</mi>
<mi>z</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>&Phi;</mi>
<mi>z</mi>
</msub>
<msubsup>
<mi>A</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
通过对进行特征值分解,得到与Φz对应的特征值以及特征向量矩阵γ,因为则z轴低精度无模糊的方向余弦估计为:
<mrow>
<mover>
<msubsup>
<mi>&upsi;</mi>
<mi>k</mi>
<mi>r</mi>
</msubsup>
<mo>^</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<mi>arg</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>&phi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>d</mi>
<mo>/</mo>
<mi>&lambda;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
ψzz,ψwz自身分别包含高精度模糊的方向余弦估计信息,因此,ψzz、ψwz分别取前M-K-1行构成维数为2(M-K-1)×K矩阵P,ψzz、ψwz后M-K-1行构成维数为2(M-K-1)×K矩阵Q,则P,Q用如下公式表示:
P=Λ1Az1 -1
Q=Λ2Az1 -1 (15)
式中,Λ1、Λ2分别表示由Az2、Aw2的前M-K-1行和后M-K-1行构成的维数为2(M-K-1)×K的矩阵;
由式(12)可知,定义B1=Pγ,B2=Qγ,即实现z轴方向上低精度无模糊的方向余弦估计与高精度模糊的方向余弦估计的一一对应,则Φz的估计值可通过B1,B2求得,即因为方向余弦的范围{-1≤υ≤1},通过对进行特征值分解,则z轴高精度模糊的方向余弦估计为:
<mrow>
<mover>
<mi>&upsi;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>z</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>z</mi>
</msub>
<mfrac>
<mi>&lambda;</mi>
<msub>
<mi>d</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>arg</mi>
<mo>{</mo>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mover>
<mi>&Phi;</mi>
<mo>^</mo>
</mover>
<mi>z</mi>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>}</mo>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;d</mi>
<mi>s</mi>
</msub>
<mo>/</mo>
<mi>&lambda;</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>K</mi>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,表示不小于i的最小整数,表示不大于i的最大整数,表示的第k个对角元素;
同理,可得到x轴y轴低精度无模糊的方向余弦估计以及对应高精度模糊的方向余弦估计
4.如权利要求1所述的基于互相关矩阵的扩展孔径二维DOA估计方法,其特征是,实现角度估计的具体步骤:
利用解模糊的方法,z轴高精度无模糊的方向余弦估计为
<mrow>
<msub>
<mover>
<mi>&upsi;</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>n</mi>
<mi>z</mi>
<mi>k</mi>
</msubsup>
<mfrac>
<mi>&lambda;</mi>
<msub>
<mi>d</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,用下式进行估计
<mrow>
<msubsup>
<mi>n</mi>
<mi>z</mi>
<mi>k</mi>
</msubsup>
<mo>=</mo>
<munder>
<mi>argmin</mi>
<msub>
<mi>n</mi>
<mi>z</mi>
</msub>
</munder>
<mo>|</mo>
<msubsup>
<mover>
<mi>&upsi;</mi>
<mo>^</mo>
</mover>
<mi>k</mi>
<mi>r</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>n</mi>
<mi>z</mi>
</msub>
<mi>&lambda;</mi>
</mrow>
<msub>
<mi>d</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>|</mo>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mi>K</mi>
<mo>.</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
同理,x轴,y轴高精度无模糊的方向余弦估计分别为
根据以上的分析,第k个信号的二维DOA估计表达式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711026670.0A CN107703478B (zh) | 2017-10-27 | 2017-10-27 | 基于互相关矩阵的扩展孔径二维doa估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711026670.0A CN107703478B (zh) | 2017-10-27 | 2017-10-27 | 基于互相关矩阵的扩展孔径二维doa估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107703478A true CN107703478A (zh) | 2018-02-16 |
CN107703478B CN107703478B (zh) | 2021-04-27 |
Family
ID=61176284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711026670.0A Expired - Fee Related CN107703478B (zh) | 2017-10-27 | 2017-10-27 | 基于互相关矩阵的扩展孔径二维doa估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107703478B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108535682A (zh) * | 2018-06-15 | 2018-09-14 | 华南理工大学 | 一种基于旋转非均匀双l阵的水下二维doa估计方法及装置 |
CN108872930A (zh) * | 2018-08-28 | 2018-11-23 | 天津大学 | 扩展孔径二维联合对角化doa估计方法 |
CN109490820A (zh) * | 2018-11-13 | 2019-03-19 | 电子科技大学 | 一种基于平行嵌套阵的二维doa估计方法 |
CN110018438A (zh) * | 2019-04-23 | 2019-07-16 | 北京邮电大学 | 一种波达方向估计方法和装置 |
CN110244258A (zh) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | 用于双平行阵二维测向中扩展doa矩阵方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101609150A (zh) * | 2009-07-07 | 2009-12-23 | 哈尔滨工程大学 | 一种提高基阵分辨力和增益的快速波束形成方法 |
CN102327131A (zh) * | 2010-07-12 | 2012-01-25 | 通用电气公司 | 用于控制超声系统中数据的通信的方法和系统 |
CN106842196A (zh) * | 2015-12-03 | 2017-06-13 | 中国航空工业集团公司雷华电子技术研究所 | 一种用于sar子孔径成像的方位模糊抑制方法 |
-
2017
- 2017-10-27 CN CN201711026670.0A patent/CN107703478B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101609150A (zh) * | 2009-07-07 | 2009-12-23 | 哈尔滨工程大学 | 一种提高基阵分辨力和增益的快速波束形成方法 |
CN102327131A (zh) * | 2010-07-12 | 2012-01-25 | 通用电气公司 | 用于控制超声系统中数据的通信的方法和系统 |
CN106842196A (zh) * | 2015-12-03 | 2017-06-13 | 中国航空工业集团公司雷华电子技术研究所 | 一种用于sar子孔径成像的方位模糊抑制方法 |
Non-Patent Citations (5)
Title |
---|
C.S. MARINO 等: "High-resolution DOA estimation from synthetic aperture beamforming", 《 2005 IEEE ANTENNAS AND PROPAGATION SOCIETY INTERNATIONAL SYMPOSIUM》 * |
JIN HE 等: "Extended aperture 2-D Direction Finding With a Two-Parallel-Shape-Array Using Propagator Method", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 * |
YUANMING GUO: "Reduced-rank covariance matrix based algorithm for DOA estimation", 《 PROCEEDINGS OF 2014 3RD ASIA-PACIFIC CONFERENCE ON ANTENNAS AND PROPAGATION》 * |
杨晋生 等: "混合信源背景下的DOA估计分步改进算法", 《信息技术》 * |
顾陈 等: "基于扩展孔径波达方向矩阵法的高精度二维DOA估计", 《电子与信息学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108535682A (zh) * | 2018-06-15 | 2018-09-14 | 华南理工大学 | 一种基于旋转非均匀双l阵的水下二维doa估计方法及装置 |
CN108535682B (zh) * | 2018-06-15 | 2023-09-29 | 华南理工大学 | 一种基于旋转非均匀双l阵的水下二维doa估计方法及装置 |
CN108872930A (zh) * | 2018-08-28 | 2018-11-23 | 天津大学 | 扩展孔径二维联合对角化doa估计方法 |
CN108872930B (zh) * | 2018-08-28 | 2022-09-30 | 天津大学 | 扩展孔径二维联合对角化doa估计方法 |
CN109490820A (zh) * | 2018-11-13 | 2019-03-19 | 电子科技大学 | 一种基于平行嵌套阵的二维doa估计方法 |
CN110018438A (zh) * | 2019-04-23 | 2019-07-16 | 北京邮电大学 | 一种波达方向估计方法和装置 |
CN110018438B (zh) * | 2019-04-23 | 2020-09-25 | 北京邮电大学 | 一种波达方向估计方法和装置 |
CN110244258A (zh) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | 用于双平行阵二维测向中扩展doa矩阵方法 |
CN110244258B (zh) * | 2019-06-12 | 2022-10-04 | 南京航空航天大学 | 用于双平行阵二维测向中扩展doa矩阵方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107703478B (zh) | 2021-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106526530B (zh) | 基于传播算子的2-l型阵列二维doa估计算法 | |
CN107703478A (zh) | 基于互相关矩阵的扩展孔径二维doa估计方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN103605108B (zh) | 声矢量阵高精度远程方位估计方法 | |
CN105137399B (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN106054123A (zh) | 一种稀疏l阵及其二维doa估计方法 | |
CN102175989B (zh) | 一种非相干分布式信号二维波达角的测定方法 | |
CN106019234B (zh) | L型天线阵的低计算复杂度二维波达方向估计方法 | |
CN106772224A (zh) | 一种采用时频分析的l型阵列二维波达方向估计算法 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN106772225A (zh) | 基于压缩感知的波束域doa估计 | |
CN108732549A (zh) | 一种基于协方差矩阵重构的阵元缺损mimo雷达doa估计方法 | |
CN105548957B (zh) | 一种未知有色噪声下多目标远近场混合源定位方法 | |
CN105044684B (zh) | 基于射频隐身的mimo跟踪雷达发射波束的形成方法 | |
CN101251597A (zh) | 一种多输入多输出雷达系统阵列误差自校正的方法 | |
CN105445709B (zh) | 一种稀布阵列近场无源定位幅相误差校正方法 | |
CN103383452A (zh) | 分布式阵列目标到达角估计方法 | |
CN105403871B (zh) | 一种双基地mimo雷达阵列目标角度估计和互耦误差校准方法 | |
CN107390197A (zh) | 基于特征空间的雷达自适应和差波束测角方法 | |
CN107462872A (zh) | 一种抗主瓣干扰算法 | |
CN103983952A (zh) | 一种非圆信号双基地mimo雷达低复杂度收发角度联合估计方法 | |
CN106353738B (zh) | 一种新的doa失配条件下稳健自适应波束形成方法 | |
CN106443570A (zh) | 基于多重信号分类算法矢量相关性的波达方向估计方法 | |
CN110161489A (zh) | 一种基于伪框架的强弱信号测向方法 | |
CN103116162A (zh) | 基于目标空间稀疏性的高分辨声呐定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210427 Termination date: 20211027 |
|
CF01 | Termination of patent right due to non-payment of annual fee |