CN108120953A - 一种基于波达方向估计的无线电定位方法 - Google Patents
一种基于波达方向估计的无线电定位方法 Download PDFInfo
- Publication number
- CN108120953A CN108120953A CN201711442433.2A CN201711442433A CN108120953A CN 108120953 A CN108120953 A CN 108120953A CN 201711442433 A CN201711442433 A CN 201711442433A CN 108120953 A CN108120953 A CN 108120953A
- Authority
- CN
- China
- Prior art keywords
- signal
- signals
- array
- matrix
- calculating
- 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.)
- Withdrawn
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供了一种基于波达方向估计的无线电定位方法。针对存在高噪声和同频干扰的复杂电磁环境下,定义了广义循环互相关熵函数。基于旋转不变信号子空间的理论,对线性阵列天线接收到的两个信号子序列构建广义循环互相关熵矩阵,通过SVD分解的方式,获得信号源的波达方向,从而对信号源进行监测定位。本发明提供了一种基于广义循环相关熵波达方向估计的无线电监测定位方法,该方法不仅能够有效的降低脉冲噪声同频带干扰信号的影响,并且具有较高的计算焦虑,对高斯噪声也有较好的性能。
Description
技术领域
本发明涉及到一种基于波达方向估计的无线电定位方法,尤其能高精度的定位脉冲噪声和同频干扰下的目标。
背景技术
波达方向(Direction of Arrival,DOA)估计作为无线被动定位的重要技术手段,在雷达、声纳以及通信等阵列信号处理领域有着广泛的应用。随着现代移动通信技术的迅速发展,电磁波的传输环境变得越来越复杂,对脉冲噪声与同频带干扰共存条件下DOA估计技术的要求也越来越高。
基于子空间的DOA估计算法不但能够同时对多个信号进行测向,较之传统算法,测向角度也能达到很高的精度,因此近年来受到广泛的关注。Gardner将循环相关矩阵应用到基于多重空间分类和旋转不变信号子空间算法之中提出了循环MUSIC算法和循环ESPRIT算法。上述两种算法都可以获得较高的谱分辨率,但循环MUSIC算法需要进行谱峰搜索,复杂度较高且降低了算法效率。循环ESPRIT算法则不存在这样的问题,因此更适合于工程应用。
针对脉冲噪声与同频带干扰问题,基于分数低阶统计量的信号处理算法可以较好地抑制脉冲噪声,但其对于同频带干扰没有抑制作用。分数低阶循环统计量算法可以有效抑制脉冲噪声与同频干扰,但对***阶数的设定依赖于Alpha稳定分布过程特征指数的先验知识,因而在实际应用中会受到一定的限制。基于相关熵的信号参数估计和无线定位算法,可以不受信号和噪声先验知识的约束,而有效抑制Alpha稳定分布噪声,但这类算法不能有效抑制同频带干扰。
发明内容
为克服现有方法进行无线电监测定位中在脉冲噪声和同频干扰环境下无法精确定位目标的不足,本发明提供了一种基于广义循环相关熵波达方向估计的无线电监测定位方法,该方法不仅能够有效的降低脉冲噪声同频带干扰信号的影响,并且具有较高的计算焦虑,对高斯噪声也有较好的性能。
本发明解决其技术问题所采用的技术方案是:在旋转不变信号子空间基础上,本发明提出了基于循环相关熵的ESPRIT算法(CCE-ESPRIT),定义了循环互相关熵估计矩阵,对线性天线阵列接收到的两组阵列信号计算其循环相关熵估计矩阵,通过特征值分解算法,估计出待监测信号的波达方向。为验证本发明的优势,仿真实验中同分数低阶ESPRIT算法(FLOM-ESPRIT)进行对比分析。
本发明的有益效果是,可以对存在同频干扰和脉冲噪声的复杂电磁环境下的无线电信号源进行监测和定位,并且算法效率高,同时适用于高斯噪声环境。
附图说明
图1是无线电信号接收天线阵列模型图。
图2是本发明与FLOM-ESPRIT算法在不同广义信噪比下的对比(非同频干扰环境)。
图3是本发明与FLOM-ESPRIT算法在不同广义信噪比下的对比(同频干扰环境)。
图4是本发明与FLOM-ESPRIT算法在不同快拍数下的对比(非同频干扰环境)。
图5是本发明与FLOM-ESPRIT算法在不同快拍数下的对比(同频干扰环境)。
具体实施方式
以下结合附图和技术方案,进一步说明本发明的具体实施方式。
一种基于波达方向估计的无线电定位方法,其特征在于,步骤如下:
第一步,建立阵列与信号模型。图1中,等距线性排列的M个阵列天线阵元、第一组M-1个子阵列、第二组M-1个子阵列。子阵列选取基于旋转不变信号子空间算法,通常选择相差一个天线位置,相邻阵元间距离满足“半波长”条件d≤λ/2,其中λ=2πc/ω,c为光的传播速度,ω为信号的中心频率。当L个互不相关的远场信号入射到阵列时,接收到信号的数学模型可表示为:
X(t)=A(θ)S(t)+U(t) (1)
其中,X(t)=[x1(t),x2(t),…,xM(t)]T为M个阵元上的数据观测矩阵,S(t)=[s1(t),s2(t),…,sL(t)]T为L个入射信号,假定s1(t),s2(t),…,sl(t)是循环频率为ε的循环平稳信号,l<L且非相干,其余的L-l个信号sl+1(t),sl+2(t),…,sL(t)具有不同的频率或者不具有循环平稳特性。A(θ)=[a(θ1),a(θ2),…,a(θL)]为入射信号的导向向量。其中θi为第i个信号的入射角。U(t)=[u1(t),u2(t),…,uM(t)]T为噪声向量,则第k个阵元接收到的信号为:
第二步,获取两组同步信号
两组同步信号表达式为:x(n)=[x1(n),x2(n),…,xM-1(n)]T和y(n)=[y1(n),y2(n),…,yM-1(n)]T,其中y(n)=x(n+1),因此还可以给出y(n)的另一种表达形式y(n)=[x2(n),x3(n),…,xM(n)]T,上述两个信号亦可以表示为向量的形式:x(n)=As(n)+u(n),y(n)=Aψs(n)+u(n+1),其中xi(n+1)=xi+1(n),ui(n+1)=ui+1(n),i=1,2,…,M-1。u(n)=[u1(n),u2(n),…,uM-1(n)]为阵列接收到的噪声向量。A为阵列的M×L响应矩阵:
第三步,计算信号的广义循环互相关熵
3.1)利用接收到的两组同步信号x(n)=[x1(n),x2(n),…,xM(n)]T和y(n)=[y1(n),y2(n),…,yM(n)]T计算广义循环互相关熵:
其中α为循环频率,τ为时延;
3.2)根据信号的广义循环互相关熵计算估计矩阵:
第四步,根据估计矩阵计算信号的波达方向:
4.1)计算选择矩阵两组接收信号为:X1=J1X=[x1(1),x1(2),…,x1(N)]和X2=J2X=[x2(1),x2(2),…,x2(N)],其中x1(n)=[x1(n),x2(n),…,xM-1(N)]T,x2(n)=[x2(n),x3(n),…,xM(N)]T。
则:
4.2)利用观测矩阵建立矩阵束{V1,V2},计算V1的特征值分解
其中∑1包含p个主特征值,计算
4.3)计算矩阵束的特征值分解,求解单位圆上的特征值,既可以得到信号的波达方向。
在图2中,噪声环境为脉冲噪声,信号源非同频干扰情况下,本发明和FLOM-ESPRIT算法进行了对比,随着广义信噪比的从0dB增加到10dB,两种算法的均方误差都随之下降,但本发明算法性能好于FLOM-ESPRIT算法。
在图3中,噪声环境为脉冲噪声,信号源存在同频干扰情况下,本发明和FLOM-ESPRIT算法进行了对比,随着广义信噪比的从0dB增加到10dB,两种算法的均方误差都随之下降,但本发明算法性能好于FLOM-ESPRIT算法。
在图4中,噪声环境为脉冲噪声,信号源非同频干扰情况下,本发明和FLOM-ESPRIT算法进行了对比,随着快拍数从100次增加到1000次,两种算法的均方误差都随之下降,但本发明算法性能好于FLOM-ESPRIT算法。
在图5中,噪声环境为脉冲噪声,信号源存在同频干扰情况下,本发明和FLOM-ESPRIT算法进行了对比,随着快拍数从100次增加到1000次,两种算法的均方误差都随之下降,但本发明算法性能好于FLOM-ESPRIT算法。
Claims (1)
1.一种基于波达方向估计的无线电定位方法,其特征在于,步骤如下:
第一步,建立阵列与信号模型
等距线性排列的M个阵列天线阵元、第一组M-1个子阵列、第二组M-1个子阵列;子阵列选取基于旋转不变信号子空间算法,选择相差一个天线位置,相邻阵元间距离满足“半波长”条件d≤λ/2,其中λ=2πc/ω,c为光的传播速度,ω为信号的中心频率;当L个互不相关的远场信号入射到阵列时,接收到信号的数学模型表示为:
X(t)=A(θ)S(t)+U(t) (1)
其中,X(t)=[x1(t),x2(t),…,xM(t)]T为M个阵列天线阵元上的数据观测矩阵,S(t)=[s1(t),s2(t),…,sL(t)]T为L个入射信号,假定s1(t),s2(t),…,sl(t)是循环频率为ε的循环平稳信号,l<L且非相干,其余的L-l个信号sl+1(t),sl+2(t),…,sL(t)具有不同的频率或不具有循环平稳特性;A(θ)=[a(θ1),a(θ2),…,a(θL)]为入射信号的导向向量;其中θi为第i个信号的入射角;U(t)=[u1(t),u2(t),…,uM(t)]T为噪声向量,则第k个阵列天线阵元接收到的信号为:
第二步,获取两组同步信号
两组同步信号表达式为:x(n)=[x1(n),x2(n),…,xM-1(n)]T和y(n)=[y1(n),y2(n),…,yM-1(n)]T,其中y(n)=x(n+1),因此还给出y(n)的另一种表达形式y(n)=[x2(n),x3(n),…,xM(n)]T,上述两个信号亦表示为向量的形式:x(n)=As(n)+u(n),y(n)=Aψs(n)+u(n+1),其中xi(n+1)=xi+1(n),ui(n+1)=ui+1(n),i=1,2,…,M-1;u(n)=[u1(n),u2(n),…,uM-1(n)]为阵列接收到的噪声向量,A为阵列的M×L响应矩阵:
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>j&omega;</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>j&omega;</mi>
<mn>2</mn>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>j&omega;</mi>
<mi>L</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&omega;</mi>
<mi>L</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
第三步,计算信号的广义循环互相关熵
3.1)利用接收到的两组同步信号x(n)=[x1(n),x2(n),…,xM(n)]T和y(n)=[y1(n),y2(n),…,yM(n)]T计算广义循环互相关熵:
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<msub>
<mi>y</mi>
<mi>j</mi>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo><</mo>
<msub>
<mi>&kappa;</mi>
<mi>&sigma;</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<mi>&tau;</mi>
<mn>2</mn>
</mfrac>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>j</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mfrac>
<mi>&tau;</mi>
<mn>2</mn>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<mi>&tau;</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mfrac>
<mi>&tau;</mi>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&alpha;</mi>
<mi>t</mi>
</mrow>
</msup>
<msub>
<mo>></mo>
<mi>t</mi>
</msub>
</mrow>
其中,α为循环频率,τ为时延;
3.2)根据信号的广义循环互相关熵计算估计矩阵:
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
第四步,根据估计矩阵计算信号的波达方向:
4.1)计算选择矩阵两组接收信号为:X1=J1X=[x1(1),x1(2),…,x1(N)]和X2=J2X=[x2(1),x2(2),…,x2(N)],其中x1(n)=[x1(n),x2(n),…,xM-1(N)]T,x2(n)=[x2(n),x3(n),…,xM(N)]T;
则:
<mrow>
<mi>J</mi>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>U</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
4.2)利用观测矩阵建立矩阵束{V1,V2},计算V1的特征值分解
<mrow>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mi>U&Sigma;Q</mi>
<mi>H</mi>
</msup>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>U</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>U</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&Sigma;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mi>O</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>O</mi>
</mtd>
<mtd>
<msub>
<mi>&Sigma;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>Q</mi>
<mn>1</mn>
<mi>H</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>Q</mi>
<mn>2</mn>
<mi>H</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中∑1包含p个主特征值,计算
4.3)计算矩阵束的特征值分解,求解单位圆上的特征值,既得到信号的波达方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711442433.2A CN108120953A (zh) | 2017-12-27 | 2017-12-27 | 一种基于波达方向估计的无线电定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711442433.2A CN108120953A (zh) | 2017-12-27 | 2017-12-27 | 一种基于波达方向估计的无线电定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108120953A true CN108120953A (zh) | 2018-06-05 |
Family
ID=62232075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711442433.2A Withdrawn CN108120953A (zh) | 2017-12-27 | 2017-12-27 | 一种基于波达方向估计的无线电定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108120953A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161491A (zh) * | 2019-06-28 | 2019-08-23 | 电子科技大学 | 一种针对微弱生命体的测距和呼吸频率估计方法 |
CN111366813A (zh) * | 2020-03-17 | 2020-07-03 | 重庆邮电大学 | 一种脉冲噪声环境下的电缆故障定位方法、装置及*** |
CN112731303A (zh) * | 2020-12-21 | 2021-04-30 | 南昌工程学院 | 非高斯噪声下的干涉阵列米波雷达及稳健测高方法与应用 |
CN113255528A (zh) * | 2021-05-28 | 2021-08-13 | 胡琼 | 一种单多目标分类方法及装置 |
CN115825863A (zh) * | 2022-12-16 | 2023-03-21 | 南京航空航天大学 | 一种冲击噪声下的非圆信号快速直接定位方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105302940A (zh) * | 2015-09-16 | 2016-02-03 | 大连理工大学 | 一种基于循环相关熵的载频估计方法 |
CN105933259A (zh) * | 2016-04-21 | 2016-09-07 | 大连理工大学 | 一种基于压缩感知重建的循环相关熵谱的载频估计方法 |
CN107517173A (zh) * | 2017-09-21 | 2017-12-26 | 大连理工大学 | 基于广义循环相关熵的无线电信号载波频率估计方法 |
-
2017
- 2017-12-27 CN CN201711442433.2A patent/CN108120953A/zh not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105302940A (zh) * | 2015-09-16 | 2016-02-03 | 大连理工大学 | 一种基于循环相关熵的载频估计方法 |
CN105933259A (zh) * | 2016-04-21 | 2016-09-07 | 大连理工大学 | 一种基于压缩感知重建的循环相关熵谱的载频估计方法 |
CN107517173A (zh) * | 2017-09-21 | 2017-12-26 | 大连理工大学 | 基于广义循环相关熵的无线电信号载波频率估计方法 |
Non-Patent Citations (4)
Title |
---|
SEN LI 等: "Robust Cyclic MUSIC Algorithm for Finding Directions in Impulsive Noise Environment", 《INTERNATIONAL JOURNAL OF ANTENNAS AND PROPAGATION》 * |
兰天 等: "脉冲噪声环境下循环ESPRIT新方法", 《通信学报》 * |
赵季红 等: "冲击噪声下基于相关熵的二维ESPRIT算法", 《计算机测量与控制》 * |
陈晓静: "稳定分布噪声下循环平稳信号的参数估计", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110161491A (zh) * | 2019-06-28 | 2019-08-23 | 电子科技大学 | 一种针对微弱生命体的测距和呼吸频率估计方法 |
CN110161491B (zh) * | 2019-06-28 | 2021-01-12 | 电子科技大学 | 一种针对微弱生命体的测距和呼吸频率估计方法 |
CN111366813A (zh) * | 2020-03-17 | 2020-07-03 | 重庆邮电大学 | 一种脉冲噪声环境下的电缆故障定位方法、装置及*** |
CN112731303A (zh) * | 2020-12-21 | 2021-04-30 | 南昌工程学院 | 非高斯噪声下的干涉阵列米波雷达及稳健测高方法与应用 |
CN112731303B (zh) * | 2020-12-21 | 2024-04-09 | 南昌工程学院 | 非高斯噪声下的干涉阵列米波雷达及稳健测高方法与应用 |
CN113255528A (zh) * | 2021-05-28 | 2021-08-13 | 胡琼 | 一种单多目标分类方法及装置 |
CN115825863A (zh) * | 2022-12-16 | 2023-03-21 | 南京航空航天大学 | 一种冲击噪声下的非圆信号快速直接定位方法 |
CN115825863B (zh) * | 2022-12-16 | 2023-12-29 | 南京航空航天大学 | 一种冲击噪声下的非圆信号快速直接定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108120953A (zh) | 一种基于波达方向估计的无线电定位方法 | |
CN109188344B (zh) | 脉冲噪声环境下基于互循环相关music算法信源个数与来波方向角估计方法 | |
CN110927661A (zh) | 基于music算法的单基地展开互质阵列mimo雷达doa估计方法 | |
CN108120967B (zh) | 一种平面阵列doa估计方法及设备 | |
Vesa | Direction of arrival estimation using music and root-music algorithm | |
CN109696657B (zh) | 一种基于矢量水听器的相干声源定位方法 | |
CN109946643B (zh) | 基于music求解的非圆信号波达方向角估计方法 | |
CN110244273A (zh) | 一种基于均匀分布式阵列的目标角度估计方法 | |
CN111983554A (zh) | 非均匀l阵下的高精度二维doa估计 | |
Li et al. | Parameter estimation based on fractional power spectrum density in bistatic MIMO radar system under impulsive noise environment | |
CN113189592A (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN112255629A (zh) | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 | |
CN113075649B (zh) | 一种适用于分布式网络化雷达的信号级直接定位方法 | |
Li et al. | DOD and DOA estimation for MIMO radar based on combined MUSIC and sparse Bayesian learning | |
Abdullah et al. | Comparative Study of Super-Performance DOA Algorithms based for RF Source Direction Finding and Tracking | |
CN112363108A (zh) | 信号子空间加权超分辨的波达方向检测方法及*** | |
CN113341371B (zh) | 一种基于l阵和二维esprit算法的doa估计方法 | |
Qin et al. | Cross-correlation matrix Root-MUSIC algorithm for bistatic multiple-input multiple-output radar | |
Ni et al. | Information-theoretic target localization with compressed measurement using FDA radar | |
Yuri et al. | Performance study of beamspace processing DOA estimation by MUSIC and Capon methods | |
Singh et al. | Near field targets localization using bistatic MIMO system with symmetric arrays | |
Qing et al. | Research on high-order cumulant ESPRIT algorithm based on single AVS | |
CN112698263A (zh) | 一种基于正交传播算子的单基地互质mimo阵列doa估计算法 | |
Nechaev et al. | Impact of beamspace processing on accuracy of DOA estimation using MUSIC and Capon methods | |
Yu et al. | DOD-DOA estimation by exploiting signal cyclostationarity for bistatic MIMO radar |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180605 |