[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN111948599B - 一种角度相关互耦影响下相干信号的高分辨率定位方法 - Google Patents

一种角度相关互耦影响下相干信号的高分辨率定位方法 Download PDF

Info

Publication number
CN111948599B
CN111948599B CN202010819681.XA CN202010819681A CN111948599B CN 111948599 B CN111948599 B CN 111948599B CN 202010819681 A CN202010819681 A CN 202010819681A CN 111948599 B CN111948599 B CN 111948599B
Authority
CN
China
Prior art keywords
array
angle
coefficient
vector
mutual coupling
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.)
Expired - Fee Related
Application number
CN202010819681.XA
Other languages
English (en)
Other versions
CN111948599A (zh
Inventor
杨朝麟
郑植
王文钦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Guangdong Electronic Information Engineering Research Institute of UESTC
Original Assignee
University of Electronic Science and Technology of China
Guangdong Electronic Information Engineering Research Institute of UESTC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China, Guangdong Electronic Information Engineering Research Institute of UESTC filed Critical University of Electronic Science and Technology of China
Priority to CN202010819681.XA priority Critical patent/CN111948599B/zh
Publication of CN111948599A publication Critical patent/CN111948599A/zh
Application granted granted Critical
Publication of CN111948599B publication Critical patent/CN111948599B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • G01S3/02—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 using radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction
    • G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
    • 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/006—Theoretical aspects
    • 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/006—Theoretical aspects
    • 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
    • G01S3/80—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 using ultrasonic, sonic or infrasonic waves
    • G01S3/802—Systems for determining direction or deviation from predetermined direction
    • G01S3/8027—By vectorial composition of signals received by plural, differently-oriented transducers

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种角度相关互耦影响下相干信号的高分辨率定位方法，根据获取的天线阵列，建立阵元角度相关互耦系数模型和相干信号的阵列数据接收模型；然后，根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合；接着，对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合；其次，对中心阵列进行子阵划分，得到线性预测系数；最后根据所述线性预测系数构建估计方程，得到对应的DOA值。本发明不仅去掉了角度相关互耦效应与多径影响，而且在相干源间隔相近时，提供算法的估计精度。

Description

一种角度相关互耦影响下相干信号的高分辨率定位方法
技术领域
本发明涉及阵列信号处理技术领域，尤其涉及一种角度相关互耦影响下相干信号的高分辨率定位方法。
背景技术
波达方向的估计是阵列信号处理中的基本问题，被广泛的应用在雷达、声呐、无线通信以及电子侦查等领域。在过去的几十年中，许多高分辨的DOA估计算法被提出，如MUSIC、ESPRIT和ROOT-MUSIC等。这些算法都需要知道准确的阵列流型矢量。但是在实际的系统中，有两类问题存在使得DOA估计性能下降：一是阵列流型矢量由于阵元间的互耦效应影响而偏离准确值；二是环境中信号的多径传播效应，使协方差矩阵秩缺损。针对互耦效应与信号源相干的问题，许多算法被提出，但都是针对角度不相关互耦情况下的DOA估计。近年来，相继有有学者发现互耦效应与来波角度相关，使上述问题中的信号建模发生改变：针对每一个来波信号，都会有不同的互耦矩阵。为了解决不同信号的互耦矩阵不相同的问题，而在相干源信号下，现有的未知互耦效应影响下相干信号源定位算法存在估计精度降低、估计自由度上减少和估计分辨率减小，导致算法的估计精度降低。
发明内容
本发明的目的在于提供一种角度相关互耦影响下相干信号的高分辨率定位方法，提供算法的估计精度。
为实现上述目的，本发明提供了一种角度相关互耦影响下相干信号的高分辨率定位方法，包括：
根据获取的天线阵列，建立阵元角度相关互耦系数模型和相干信号的阵列数据接收模型；
根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合；
对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合；
对中心阵列进行子阵划分，得到线性预测系数；
根据所述线性预测系数构建估计方程，得到对应的DOA值。
其中，根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合，包括：
根据角度相关互耦矩阵的拓普利兹对称特性，将得到的选择矩阵与所述天线阵列接收的原始数据相乘后，更新得到的去角度相关互耦接收数据中的相干系数，得到对应的阵列流形矢量线性组合。
其中，对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合，包括：
根据得到的所述去角度相关互耦接收数据，计算出在对应快拍数下的协方差矩阵，并对所述协方差矩阵进行特征值分解后，对等式两边同时进行归一化处理，得到对应的归一化线性组合。
其中，对中心阵列进行子阵划分，得到线性预测系数，包括：
对具有多个阵元的中心阵列进行子阵划分，得到多个子阵，并根据所述子阵划分方法，将所述归一化线性组合进行分组，并根据所述归一化线性组合中的元素关系构建前后向线性预测模型。
其中，对中心阵列进行子阵划分，得到线性预测系数，还包括：
将前向线性预测方程和后向线性预测方程进行堆叠，并改写成矩阵形式，并根据得到的线性估计方程计算出线性预测系数。
其中，根据所述线性预测系数构建估计方程，得到对应的DOA值，包括：
将所述归一化线性组合中的任意元素代入所述前向线性预测模型中，并进行简化，得到第一系数，并将第二系数对应的向量代入所述第一系数中，得到所述线性预测系数与划分得到的所述子阵的阵列导向矢量垂直。
其中，根据所述线性预测系数构建估计方程，得到对应的DOA值，还包括：
将与所述阵列导向矢量垂直的所述线性预测系数代入构造的估计方程中，并根据所述估计方程的根计算出对应的DOA值。
本发明的一种角度相关互耦影响下相干信号的高分辨率定位方法，根据获取的天线阵列，建立阵元角度相关互耦系数模型和相干信号的阵列数据接收模型；然后，根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合；接着，对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合；其次，对中心阵列进行子阵划分，得到线性预测系数；最后根据所述线性预测系数构建估计方程，得到对应的DOA值。本发明不仅去掉了角度相关互耦效应与多径影响，而且在相干源间隔相近时，提供算法的估计精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的一种角度相关互耦影响下相干信号的高分辨率定位方法的步骤示意图。
图2是本发明提供的在不同的算法下，角度相关互耦影响下能够估计的相干信号的数目。
图3是本发明提供的在相同快拍数，不同SNR下的不同方法的估计性能比较。
图4是本发明提供的在相同SNR，不同快拍数下的不同方法的估计性能比较。
图5是本发明提供的不同算法的分辨率比较。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
请参阅图1，本发明提供一种角度相关互耦影响下相干信号的高分辨率定位方法，包括：
S101、根据获取的天线阵列，建立阵元角度相关互耦系数模型和相干信号的阵列数据接收模型。
具体的，设置天线阵列，阵列为均匀线阵，阵元数为M，阵元间的间隔为d＝λ/2，λ为信号载波波长。该阵列接收K个窄带相干平稳信号，信号的相干系数为β＝[β1,β2,......βK]。阵元间存在角度相关互耦效应，参考阵元受到第l个阵元第k个信号的互耦系数为Ckl。因此，第k个信号所产生的互耦矩阵为C(θK)＝Toeplitz([cK,O])∈CM×M，具有拓普利兹对称性质。C(θK)中的参数值由行向量ck＝[ck0,ck1,......ckL]决定，L为对第一个阵元存在互耦效应的最大阵元个数。cK的取值满足|ck0|＝1＞|ck1|＞...＞|ckL|。噪声为独立同分布的高斯白噪声，阵列的观测数据矢量为：
Figure BDA0002634010740000041
其中，
Figure BDA0002634010740000042
是收到角度相关互耦效应影响的阵列流形矢量，s1(t)相干信号中的参考信号源，n(t)为噪声项。理想状态下的导向矢量为：
Figure BDA0002634010740000043
S102、根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合。
具体的，根据角度相关互耦矩阵的拓普利兹对称特性，定义选择矩阵：
J＝[0(M-2L)×L,1M-2L,0(M-2L)×L]∈R(M-2L)×M
为保证正确估计信号来波方向，需满足M-2L＞K。将该矩阵与所述天线阵列接收原始数据相乘，得到：
Figure BDA0002634010740000044
其中n1(t)＝Jn1(t)。进而推到可以得到：
Figure BDA0002634010740000045
其中，A1＝[α1(θ1),α1(θ2),...α1(θK)]是消除角度相关互耦效应的孔径缩小2L后中间子阵列的阵列流形矢量，中间子阵的导向矢量为：
Figure BDA0002634010740000046
其中，Λ＝diag(f(θ1),f(θ2),...f(θK))，f(θk)＝ckL+...+ckouL(θk)+...+ckou2L(θk)是消除互耦效应带来的复常数。进一步得到去角度相关互耦效应的接收数据为：
y1(t)＝A1Λβs1(t)+n1(t)＝αs1(t)+n1(t)
得到中心子阵的阵列流形矢量线性组合：
Figure BDA0002634010740000051
其中，
Figure BDA0002634010740000052
是中间子阵接收的相干信号的相干系数。因为去除角度相关互耦效应引入了复常数f(θ1),f(θ2),...f(θK)，所以该系数与原始阵列接收的相干系数不同，看做是中间子阵接收信号的相干系数。
S103、对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合。
具体的，计算解耦后的观测数据y1(t)在快拍数为T的协方差矩阵：
Figure BDA0002634010740000053
其中，
Figure BDA0002634010740000054
与
Figure BDA0002634010740000055
分别为信号功率与噪声功率。与α相对应的归一化值向量为：
Figure BDA0002634010740000056
接着，对协方差矩阵进行特征值分解：
Figure BDA0002634010740000057
其中，∑＝diag(γ1,γ2,...γM-2L)，γi是矩阵
Figure BDA0002634010740000058
的特征值，U＝[ν1,ν2,...νM-2L]是与其相对应特征向量。且满足γ1＞γ2＝...＝γM-2L＝0。
将上面的关于协方差矩阵的等式进行比较，得到：
Figure BDA0002634010740000059
将该等式左右的向量以第一个元素值为基准进行归一化，可得：
Figure BDA00026340107400000510
其中，
Figure BDA00026340107400000511
表示归一化后的特征值向量，[ν1]1为该特征值向量的第一个值。因为等式左右两边的向量进行了归一化，即
Figure BDA00026340107400000512
所以可以得到：
Figure BDA00026340107400000513
进一步，就可以推到得到：
Figure BDA00026340107400000514
即消除角度相关互耦影响后的中心子阵列的阵列流形矢量的线性组合的归一化值与中心阵列的归一化后的特征值向量相等。由此，估计得到了归一化线性组合
Figure BDA00026340107400000515
S104、对中心阵列进行子阵划分，得到线性预测系数。
具体的，将具有M-2L个阵元的中心阵列进行子阵划分，每个子阵含有M1个阵元，共有2L-M1+1。其中，M1＞K+1,保证每个子阵的阵元数目大于信号源数。根据该子阵划分的原则，相应的对其阵列流形矢量的所述归一化线性组合进行划分，即
Figure BDA00026340107400000612
中的元素进行分组。然后每一组，对应一个子阵列，该子阵列中用前M1-1个元素去线性表出第M1个元素，或者用后M1-1个元素线性表出第1个元素，得到前后向线性预测模型。对于第n个子阵列，
Figure BDA0002634010740000061
Figure BDA0002634010740000062
其中，
Figure BDA0002634010740000063
Figure BDA0002634010740000064
Figure BDA0002634010740000065
为线性预测系数。
构建的线性预测方程进行堆叠。一共有N个子阵，前后向线性预测方程总共有2N个，将这2N方程写成矩阵形式：
z＝Φρ
其中，
Figure BDA0002634010740000066
产生线性预测系数ρ的估计方程，对ρ进行估计：
Figure BDA0002634010740000067
S105、根据所述线性预测系数构建估计方程，得到对应的DOA值。
具体的，构建D(z)方程进行DOA估计之前，需要说明ρ与所划分的子阵列的阵列导向矢量垂直，这样才能保证D(z)所构建方程的根为z＝ej2πdsinθ/λ，其中θ为需要估计的角度。下面进行证明：
Figure BDA0002634010740000068
归一化后的中心子阵的阵列流形矢量的所述归一化线性组合为：
Figure BDA0002634010740000069
相应的，
Figure BDA00026340107400000610
的第i个元素为
Figure BDA00026340107400000611
其中，
Figure BDA0002634010740000071
Figure BDA0002634010740000072
将
Figure BDA0002634010740000073
的表达式代入前向线性预测模型中，得
Figure BDA0002634010740000074
将该式子进行化简：
Figure BDA0002634010740000075
其中，
Figure BDA0002634010740000076
且满足M1＞K，rank(B)＝K。
由于第一系数B是行满秩的，所以可以得到
Figure BDA0002634010740000077
将第二系数bn所对应的向量代入上述式子中，可得
Figure BDA0002634010740000078
由此可得：
Figure BDA0002634010740000079
将
Figure BDA00026340107400000710
代入D(z)，得到方程：
Figure BDA00026340107400000711
其中，
Figure BDA00026340107400000712
为ρ中的元素，z＝ej2πdsinθ/λ解上述方程，根据z可求得DOA。
图2中所示为验证该实验是为了验证在不同的算法下，角度相关互耦影响下能够估计的相干信号的数目。设有K＝5个相干源到达，来波方向均匀分布在(-80°,30°)之间。角度相关互耦系数对应于不同的来波方向，分别为c1＝[1,0.22+0.59i,0.41-0.33i]，c2＝[1,0.38+0.48i,0.23-0.16i]，c3＝[1,0.47+0.39i,0.34+0.23i]，c4＝[1,0.12+0.79i,0.21+0.27i]，c5＝[1,0.40+0.19i,0.21+0.36i]，。SNR设置为15dB，快拍数为T＝500。将算法ESPRIT-LIKE、算法SS和所提出的算法进行对比。图2中所示为在10次独立重复下的各种算法对K＝5个相干源的估计情况。方框代表的是
Figure BDA0002634010740000081
θk为来波方向准确值；叉代表
Figure BDA0002634010740000082
θk表示的是各种算法下的来波方向的估计值。可以看到所提出的算法成功在角度相关互耦影响下估计出5个相干信号的来波方向。
图3与图4比较所提出的算法与已有算法的RMSE(均方根误差)进行比较。RMSE的定义为：
Figure BDA0002634010740000083
其中，
Figure BDA0002634010740000084
为θk的第i次独立实验估计值，I＝300为独立实验次数。
这两个实验中，设置了两个相干源，来波方向分别为θ1＝12°与θ2＝19°；相对应的互耦系数分别为c1＝[1,0.60+0.51i,0.30-0.30i]，c2＝[1,0.41+0.30i,0.20+0.13i]。
图3中快拍设置为200，SNR为-10dB到14dB。可以发现随着SNR的提高，除了ESPRIT-LIKE方法，其他方法的估计性能都有进一步提升。在-10dB到14dB间，所提出的方法比其余方法具有更高的估计精度。
图4中，SNR设置为5dB，快拍从10变化到600。随着快拍数的增加，SS方法和DDD方法获得了越来越准确的和稳定的估计性能，但是ESPRIT-like算法估计失败。在整个快拍数变化的区间内，可以看到所提出的方法获得了最为准确的估计精度。
图5是为了比较算法的分辨率性能。该实验中考虑两个相干源，其中一个信号源的来波方向被设置为θ1＝5°，另一个来波方向的角度设置为θ2＝θ1+Δθ，Δθ从1°变化到15°。所对应的角度相关互耦系数为c1＝[1,0.60+0.50i,0.30-0.20i]，c2＝μ[1μ,0.41+0.30i,0.20+0.10i]，其中，μ是[0.8,1.2]件满足均匀分布的随机变量。SNR为5dB，快拍数为200。从图中，我们可以看到ESPRIT-like方法估计准确度不随SNR角度变宽而变准确，其余的方法随着角度间隔增大，RMSE下降。本发明所提出的方法能够在角度间隔较小的时候实现准确的DOA估计，并且RMSE性能比其他方法更好。
本发明的一种角度相关互耦影响下相干信号的高分辨率定位方法，根据获取的天线阵列，建立阵元角度相关互耦系数模型和相干信号的阵列数据接收模型；然后，根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合；接着，对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合；其次，对中心阵列进行子阵划分，得到线性预测系数；最后根据所述线性预测系数构建估计方程，得到对应的DOA值。本发明不仅去掉了角度相关互耦效应与多径影响，而且在相干源间隔相近时，提供算法的估计精度。
以上所揭露的仅为本发明一种较佳实施例而已，当然不能以此来限定本发明之权利范围，本领域普通技术人员可以理解实现上述实施例的全部或部分流程，并依本发明权利要求所作的等同变化，仍属于发明所涵盖的范围。

Claims (1)

1.一种角度相关互耦影响下相干信号的高分辨率定位方法，其特征在于，包括：
根据获取的天线阵列，建立阵元角度相关互耦系数模型和相干信号的阵列数据接收模型；
根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合；
对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合；
对中心阵列进行子阵划分，得到线性预测系数；
根据所述线性预测系数构建估计方程，得到对应的DOA值；
根据获取的天线阵列，建立阵元角度相关互耦系数模型和相干信号的阵列数据接收模型，包括：
设置天线阵列，阵列为均匀线阵，阵元数为M，阵元间的间隔为d＝λ/2，λ为信号载波波长；该阵列接收K个窄带相干平稳信号，信号的相干系数为β＝[β1,β2,......βK]；阵元间存在角度相关互耦效应，参考阵元受到第l个阵元第k个信号的互耦系数为Ckl；因此，第k个信号所产生的互耦矩阵为C(θK)＝Toeplitz([cK,O])∈CM×M，具有拓普利兹对称性质，C(θK)中的参数值由行向量ck＝[ck0,ck1,......ckL]决定，L为对第一个阵元存在互耦效应的最大阵元个数；cK的取值满足|ck0|＝1＞|ck1|＞...＞|ckL|；噪声为独立同分布的高斯白噪声，阵列的观测数据矢量为：
Figure FDA0003691269760000011
其中，
Figure FDA0003691269760000012
是收到角度相关互耦效应影响的阵列流形矢量，s1(t)相干信号中的参考信号源，n(t)为噪声项，理想状态下的导向矢量为：
Figure FDA0003691269760000013
根据角度相关互耦矩阵的拓普利兹对称特性，构建阵列流形矢量线性组合，包括：
根据角度相关互耦矩阵的拓普利兹对称特性，定义选择矩阵：
J＝[0(M-2L)×L,1M-2L,0(M-2L)×L]∈R(M-2L)×M
为保证正确估计信号来波方向，M-2L＞K；将该矩阵与所述天线阵列接收原始数据相乘，得到：
Figure FDA0003691269760000014
其中n1(t)＝Jn1(t)，进而推到可以得到：
Figure FDA0003691269760000021
其中，A1＝[α1(θ1),α1(θ2),...α1(θK)]是消除角度相关互耦效应的孔径缩小2L后中间子阵列的阵列流形矢量，中间子阵的导向矢量为：
Figure FDA0003691269760000022
其中，Λ＝diag(f(θ1),f(θ2),...f(θK))，f(θk)＝ckL+...+ckouL(θk)+...+ckou2L(θk)是消除互耦效应带来的复常数，去角度相关互耦效应的接收数据为：
y1(t)＝A1Λβs1(t)+n1(t)＝αs1(t)+n1(t)
得到中心子阵的阵列流形矢量线性组合：
Figure FDA0003691269760000023
其中，
Figure FDA0003691269760000024
是中间子阵接收的相干信号的相干系数，因为去除角度相关互耦效应引入了复常数f(θ1),f(θ2),...f(θK)，所以该系数与原始阵列接收的相干系数不同，看做是中间子阵接收信号的相干系数；
对多个信号的阵列流形矢量进行归一化处理，得到对应的归一化线性组合，包括：
计算解耦后的观测数据y1(t)在快拍数为T的协方差矩阵：
Figure FDA0003691269760000025
其中，
Figure FDA0003691269760000026
与
Figure FDA0003691269760000027
分别为信号功率与噪声功率，与α相对应的归一化值向量为：
Figure FDA0003691269760000028
接着，对协方差矩阵进行特征值分解：
Figure FDA0003691269760000029
其中，∑＝diag(γ1,γ2,...γM-2L)，γi是矩阵
Figure FDA00036912697600000210
的特征值，U＝[ν1,ν2,...νM-2L]是与其相对应特征向量，且满足γ1＞γ2＝...＝γM-2L＝0；
将上面的关于协方差矩阵的等式进行比较，得到：
Figure FDA00036912697600000211
将该等式左右的向量以第一个元素值为基准进行归一化，可得：
Figure FDA00036912697600000212
其中，
Figure FDA00036912697600000213
表示归一化后的特征值向量，[ν1]1为该特征值向量的第一个值；因为等式左右两边的向量进行了归一化，即
Figure FDA00036912697600000214
所以可以得到：
Figure FDA00036912697600000215
得到：
Figure FDA00036912697600000216
消除角度相关互耦影响后的中心子阵列的阵列流形矢量的线性组合的归一化值与中心阵列的归一化后的特征值向量相等；得到了归一化线性组合
Figure FDA00036912697600000217
对中心阵列进行子阵划分，得到线性预测系数，包括：
将具有M-2L个阵元的中心阵列进行子阵划分，每个子阵含有M1个阵元，共有2L-M1+1；其中，M1＞K+1，保证每个子阵的阵元数目大于信号源数；根据该子阵划分的原则，相应的对其阵列流形矢量的所述归一化线性组合进行划分，即
Figure FDA0003691269760000031
中的元素进行分组；然后每一组，对应一个子阵列，该子阵列中用前M1-1个元素去线性表出第M1个元素，或者用后M1-1个元素线性表出第1个元素，得到前后向线性预测模型，对于第n个子阵列，
Figure FDA0003691269760000032
Figure FDA0003691269760000033
其中，
Figure FDA0003691269760000034
Figure FDA0003691269760000035
Figure FDA0003691269760000036
为线性预测系数；
构建的线性预测方程进行堆叠，一共有N个子阵，前后向线性预测方程总共有2N个，将这2N方程写成矩阵形式：
z＝Φρ
其中，
Figure FDA0003691269760000037
产生线性预测系数ρ的估计方程，对ρ进行估计：
Figure FDA0003691269760000038
根据所述线性预测系数构建估计方程，得到对应的DOA值，包括：
构建D(z)方程进行DOA估计之前，ρ与所划分的子阵列的阵列导向矢量垂直，D(z)所构建方程的根为z＝ej2πdsinθ/λ，其中θ为需要估计的角度：
Figure FDA0003691269760000039
归一化后的中心子阵的阵列流形矢量的所述归一化线性组合为：
Figure FDA00036912697600000310
相应的，
Figure FDA00036912697600000311
的第i个元素为
Figure FDA00036912697600000312
其中，
Figure FDA00036912697600000313
Figure FDA00036912697600000314
将
Figure FDA00036912697600000315
的表达式代入前向线性预测模型中，得
Figure FDA00036912697600000316
将该式子进行化简：
Figure FDA0003691269760000041
其中，
Figure FDA0003691269760000042
且满足M1＞K，rank(B)＝K；
由于第一系数B是行满秩的，得到
Figure FDA0003691269760000043
将第二系数bn所对应的向量代入上述式子中，可得
Figure FDA0003691269760000044
由此可得：
Figure FDA0003691269760000045
将
Figure FDA0003691269760000046
代入D(z)，得到方程：
Figure FDA0003691269760000047
其中，
Figure FDA0003691269760000048
为ρ中的元素，z＝ej2πdsinθ/λ解上述方程，根据z可求得DOA。
CN202010819681.XA 2020-08-14 2020-08-14 一种角度相关互耦影响下相干信号的高分辨率定位方法 Expired - Fee Related CN111948599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010819681.XA CN111948599B (zh) 2020-08-14 2020-08-14 一种角度相关互耦影响下相干信号的高分辨率定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010819681.XA CN111948599B (zh) 2020-08-14 2020-08-14 一种角度相关互耦影响下相干信号的高分辨率定位方法

Publications (2)

Publication Number Publication Date
CN111948599A CN111948599A (zh) 2020-11-17
CN111948599B true CN111948599B (zh) 2022-08-19

Family

ID=73342381

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010819681.XA Expired - Fee Related CN111948599B (zh) 2020-08-14 2020-08-14 一种角度相关互耦影响下相干信号的高分辨率定位方法

Country Status (1)

Country Link
CN (1) CN111948599B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113987733B (zh) * 2020-11-27 2024-05-14 猪草微(深圳)电子有限公司 一种基于线性预测的信源方向鲁棒定位算法
CN114624665B (zh) * 2022-03-24 2023-11-07 电子科技大学 基于动态参数迭代优化的互耦误差doa自校正方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1532860A (en) * 1976-04-15 1978-11-22 Int Standard Electric Corp Holographic tele-communication switching system
EP0037725A1 (en) * 1980-04-04 1981-10-14 McDONNELL DOUGLAS CORPORATION Apparatus for determining the position of a source of electromagnetic radiation relative to a remote object
DE69933807D1 (de) * 1998-03-24 2006-12-14 Nec Corp Peilgerät zum Verarbeiten von Messdaten
WO2010072494A1 (fr) * 2008-12-23 2010-07-01 Thales Procede de determination des angles d'arrivee en azimut et en elevation de sources coherentes
CN103076609A (zh) * 2012-11-14 2013-05-01 武汉德澳科技有限公司 一种固态小型微波电控扫描行车雷达装置
CN105068041A (zh) * 2015-08-28 2015-11-18 哈尔滨工程大学 互耦条件下基于协方差矢量稀疏表示的单基地mimo雷达角度估计方法
CN107290730A (zh) * 2017-05-10 2017-10-24 长江大学 互耦条件下双基地mimo雷达角度估算方法
CN107450047A (zh) * 2017-08-28 2017-12-08 电子科技大学 嵌套阵下基于未知互耦信息的压缩感知doa估计方法
CN108141260A (zh) * 2015-08-18 2018-06-08 诺基亚通信公司 人工相互耦合的天线阵列
CN108680891A (zh) * 2018-01-05 2018-10-19 大连大学 非均匀噪声条件下考虑互耦效应的doa估计方法
CN109239649A (zh) * 2018-04-04 2019-01-18 唐晓杰 一种阵列误差条件下的互质阵列doa估计新方法
CN110109051A (zh) * 2019-04-09 2019-08-09 天津大学 基于频控阵的互耦阵列doa估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050195103A1 (en) * 2004-01-13 2005-09-08 Davis Dennis W. Phased arrays exploiting geometry phase and methods of creating such arrays

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1532860A (en) * 1976-04-15 1978-11-22 Int Standard Electric Corp Holographic tele-communication switching system
EP0037725A1 (en) * 1980-04-04 1981-10-14 McDONNELL DOUGLAS CORPORATION Apparatus for determining the position of a source of electromagnetic radiation relative to a remote object
DE69933807D1 (de) * 1998-03-24 2006-12-14 Nec Corp Peilgerät zum Verarbeiten von Messdaten
WO2010072494A1 (fr) * 2008-12-23 2010-07-01 Thales Procede de determination des angles d'arrivee en azimut et en elevation de sources coherentes
CN103076609A (zh) * 2012-11-14 2013-05-01 武汉德澳科技有限公司 一种固态小型微波电控扫描行车雷达装置
CN108141260A (zh) * 2015-08-18 2018-06-08 诺基亚通信公司 人工相互耦合的天线阵列
CN105068041A (zh) * 2015-08-28 2015-11-18 哈尔滨工程大学 互耦条件下基于协方差矢量稀疏表示的单基地mimo雷达角度估计方法
CN107290730A (zh) * 2017-05-10 2017-10-24 长江大学 互耦条件下双基地mimo雷达角度估算方法
CN107450047A (zh) * 2017-08-28 2017-12-08 电子科技大学 嵌套阵下基于未知互耦信息的压缩感知doa估计方法
CN108680891A (zh) * 2018-01-05 2018-10-19 大连大学 非均匀噪声条件下考虑互耦效应的doa估计方法
CN109239649A (zh) * 2018-04-04 2019-01-18 唐晓杰 一种阵列误差条件下的互质阵列doa估计新方法
CN110109051A (zh) * 2019-04-09 2019-08-09 天津大学 基于频控阵的互耦阵列doa估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"A Low Complexity Algorithm for Direction of Arrival Estimation With Direction-Dependent Mutual Coupling";Q. Ge 等;《IEEE Commun. Lett.》;20200131;90-94 *
"Robust Adaptive Beamforming Against Mutual Coupling Based on Mutual Coupling Coefficients Estimation";Zhi Zheng 等;《IEEE Transactions on Vehicular Technology》;20171231;9124-9133 *
"Robust DOA Estimation Against Mutual Coupling With Nested Array";Zhi Zheng 等;《IEEE Signal Processing Letters》;20200722;2917-2921 *
"互耦误差条件下的阵列信号参数估计研究";谭親林;《中国优秀硕士论文全文数据库》;20101015;全文 *

Also Published As

Publication number Publication date
CN111948599A (zh) 2020-11-17

Similar Documents

Publication Publication Date Title
CN109444810B (zh) 一种非负稀疏贝叶斯学习框架下的互质阵列非网格doa估计方法
CN109738854B (zh) 一种天线阵列来波方向的到达角估计方法
CN109490820B (zh) 一种基于平行嵌套阵的二维doa估计方法
CN110244272B (zh) 基于秩一去噪模型的波达方向估计方法
CN113835063B (zh) 一种无人机阵列幅相误差与信号doa联合估计方法
CN111948599B (zh) 一种角度相关互耦影响下相干信号的高分辨率定位方法
CN110174658B (zh) 基于秩一降维模型和矩阵补全的波达方向估计方法
CN109946643B (zh) 基于music求解的非圆信号波达方向角估计方法
CN110806557A (zh) 一种冲击噪声背景下单基地mimo雷达doa估计方法
CN111983554A (zh) 非均匀l阵下的高精度二维doa估计
CN110286350A (zh) 一种l型稀疏阵doa估计的精确配对方法及装置
CN110286351A (zh) 一种基于l型嵌套阵的二维doa估计方法及装置
CN113189592A (zh) 考虑幅相互耦误差的车载毫米波mimo雷达测角方法
CN110212966A (zh) 一种相干源条件下基于重要性重采样的天线互耦校正方法
CN112731280A (zh) 互质阵列混合噪声环境下的esprit-doa估计方法
CN111368256B (zh) 一种基于均匀圆阵的单快拍测向方法
CN114167347A (zh) 冲击噪声环境下互质阵列的幅相误差校正和测向方法
CN114460531A (zh) 一种均匀线阵music空间谱估计方法
CN106877918B (zh) 互耦条件下稳健自适应波束形成方法
CN111366891B (zh) 一种基于伪协方差矩阵的均匀圆阵单快拍测向方法
CN116953600A (zh) 一种基于离网格稀疏贝叶斯学习的doa估计方法
CN112305494B (zh) 一种传感器位置校准方法及装置
Ma et al. Computation-efficient 2-D DOA estimation algorithm with array motion strategy
CN115421098A (zh) 嵌套面阵下降维求根music的二维doa估计方法
CN114184999B (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: 20220819

CF01 Termination of patent right due to non-payment of annual fee