CN110320490B - 一种无直达信号条件下的无线电波达方向估计方法 - Google Patents
一种无直达信号条件下的无线电波达方向估计方法 Download PDFInfo
- Publication number
- CN110320490B CN110320490B CN201910486779.5A CN201910486779A CN110320490B CN 110320490 B CN110320490 B CN 110320490B CN 201910486779 A CN201910486779 A CN 201910486779A CN 110320490 B CN110320490 B CN 110320490B
- Authority
- CN
- China
- Prior art keywords
- signal
- arrival
- generalized
- entropy
- representing
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 9
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 4
- 238000005315 distribution function Methods 0.000 claims description 4
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 2
- 238000003491 array Methods 0.000 claims description 2
- 238000009792 diffusion process Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- LVNGJLRDBYCPGB-LDLOPFEMSA-N (R)-1,2-distearoylphosphatidylethanolamine Chemical compound CCCCCCCCCCCCCCCCCC(=O)OC[C@H](COP([O-])(=O)OCC[NH3+])OC(=O)CCCCCCCCCCCCCCCCC LVNGJLRDBYCPGB-LDLOPFEMSA-N 0.000 abstract description 11
- 238000012545 processing Methods 0.000 abstract description 4
- 230000003595 spectral effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
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/04—Details
- G01S3/10—Means for reducing or compensating for quadrantal, site, or like errors
-
- 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
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
本发明属于信号处理技术领域,涉及到一种无直达信号条件下的无线电波达方向估计方法,尤其能高精度的估计无线电传播过程中存在多径、遮挡和脉冲噪声干扰共存情况下目标的波达方向。本发明在子空间理论的基础上,基于广义自相关熵的DSPE算法,建立分布式信源模型,定义广义自相关熵,将广义自相关熵扩展到DSPE算法中,对线性阵列天线接收信号计算其广义相关熵的协方差矩阵,通过奇异值分解建立二维空间谱,通过谱峰搜索估计出分布式信源的波达方向。本发明的方法具有高鲁棒性,同时也适用于高斯噪声。
Description
技术领域
本发明属于信号处理技术领域,涉及到一种无直达信号条件下的无线电波达方向估计方法,尤其能高精度的估计无线电传播过程中存在多径、遮挡和脉冲噪声干扰共存情况下目标的波达方向。
背景技术
在声呐、无线通信和雷达等阵列信号处理领域,波达方向估计作为被动无线电定位的重要技术,得到了广泛的应用。波达方向估计的基本原理是:对于一个待定位的远场信号源,信号到达接收阵列的不同阵元会产生一个波程差,这个波程差导致各个阵元接收到的信号之间存在一个相位差,利用这个相位差就可以求得信号源的波达方向。
随着科学技术的高速发展,越来越多的电子设备得到应用,导致空间电磁环境越来越复杂;同时随着无线电业务的增长,对阵列信号处理技术提出了更高的要求。传统的波达方向估计方法一般采用点信源模型对目标信号源进行建模,也就是假设目标信号源和观测阵列之间存在着视距传播,在这样的假设情况下,传统波达方向估计算法可以获得很高的估计精度和鲁棒性。但当信号源传播过程中多径数目过多或者存在较多的反射、折射和遮挡情况的时候,传统的点源模型就不能够准确地描述阵列所接收到的数据,导致信号源波达方向估计性能严重下降,甚至无法获得正确结果。分布式信源模型作为有效解决多径和遮挡情况下信源建模被提出并得到了广泛地应用,其准确地描述了空间目标源尺度不能忽略情况下的信源特性,有效地抑制了多径和遮挡情况,可以实现精确地波达方向估计。
以往的研究中,大多将DOA估计中存在的噪声假设为服从高斯(Gaussian)分布的加性噪声,这种假设在大多数情况下是合理的,但研究发现,诸如低频大气噪声、电话线路噪声、水声信号噪声、雷达反向散射回波以及由于海浪和山脉的不连续性引起的阵列输出的信号中都表现出典型的脉冲特性。这些种类的噪声经常严重影响高斯假设下优化的系统的性能;也就是说,如果噪声统计特性偏离高斯分布,则可能发生参数估计性能的严重恶化。研究发现,alpha稳定分布能够更准确地描述脉冲噪声。
针对脉冲噪声,分数低阶矩和相位分数低阶矩的波达方向估计算法可以较好地抑制脉冲噪声,但对系统阶数的设定依赖于alpha稳定分布过程特征指数的先验知识;基于相关熵的信号参数估计和无线定位算法由于对核函数核长的选择依赖于经验,因而在实际应用中均会受到一定的限制。
发明内容
为克服现有方法进行波达方向估计中针对多径、遮挡和脉冲噪声抑制方面的不足,本发明提供了一种无直达信号条件下的无线电波达方向估计方法。该方法采用分布式信源模型对信号进行建模,不仅能够有效地解决多径和遮挡条件下没有直达信号对波达方向估计产生的影响,同时能够有效地抑制脉冲噪声,对高斯噪声也有较好的性能。
本发明的技术方案:
一种无直达信号条件下的无线电波达方向估计方法,在子空间理论的基础上,基于广义自相关熵的DSPE算法,建立分布式信源模型,定义广义自相关熵,将广义自相关熵扩展到DSPE算法中,对线性阵列天线接收信号计算其广义相关熵的协方差矩阵,通过奇异值分解建立二维空间谱,通过谱峰搜索估计出分布式信源的波达方向。
为验证本发明的优势,仿真实验中同DSPE算法进行对比分析。
本发明的有益效果:本发明可以对电磁波传播过程中存在多径和遮挡情况下没有直达波从目标源到达接收传感器并且信道存在脉冲噪声的无线电信号源进行监波达方向估计,且方法具有高鲁棒性,同时也适用于高斯噪声。
附图说明
图1是分布式信源模型图。
图2是本发明与DSPE算法在不同广义信噪比下的对比图。
图3是本发明与DSPE算法在不同特征指数下的对比图。
具体实施方式
以下结合附图和技术方案,进一步说明本发明的具体实施方式。
一种无直达信号条件下的无线电波达方向估计方法,步骤如下:
第一步,建立分布源信号模型。
图1中,等距线性阵列的M个阵元,相邻阵元间距离满足“半波长”条件d≤λ/2,其中波长λ=2πc/ω,c为光的传播速度,ω为信号的中心频率;当L个互不相关的分布式信号入射到等距线阵时,接收到信号的数学模型表示为:
其中,si(θ,ψi)是第i个分布式信源的扩散信号密度,θ∈[-π/2,π/2]是角度方向;ψi是待估计的参数,包含中心波达方向和扩散角;Θ=[-π/2,π/2]为角度的观测域,n(t)=[n1(t),n2(t),…,nM(t)]T是M×1的加性测量噪声,a(θ)=[1,e-j2πd/λsinθ,…,e-j2π(M-1)d/λsinθ]T为入射信号的点源模型下的导向向量,j为复数虚部表示因子。
在分布式信源模型下,分布式信源的扩散信号密度定义为:
si(θ,ψi)=δig(θ,ψi) (2)
其中,g(θ,ψi)表示确定角信号密度,δi是一个随机变量,因此,式(1)转换为式(3)的形式为:
其中,b(ψi)为分布式信源的广义导向向量,见式(4):
b(ψi)=∫θ∈Θa(θ)g(θ,ψi)dθ (4)
在小角度展开的假设下,对于给定的角信号分布函数,广义导向向量的解析表达式为:
其中,a(θ)为点源导向向量,m表示阵列中的第m个阵元,θ′表示积分变量,θi表示点源情况下第i个的波达方向,ρ(θ′;ψi)表示确定性角信号分布函数;
一般情况下,分布式信源的角密度函数符合高斯分布,其概率密度函数表示为:
其中,σi为标准差;
则广义导向向量表示为:
第二步,计算广义自相关熵
定义随机变量X服从alpha稳定分布,则广义自相关熵定义为:
Gσ(X)=E[κσ(|X|-μX)] (9)
其中,E[·]为求数学期望,κσ为核函数,通常选取高斯核或指数核,μX表示为:
其中,N为快拍数,p为(0,2]区间的常数,xi为阵列中第i个传感器接收到的信号序列。
第三步,计算信号的广义相关熵的协方差矩阵
利用等距线性阵列的每个阵元接收到的信号x(t)=[x1(t),x2(t),…,xM(t)]T计算脉冲噪声环境下基于广义相关熵的协方差矩阵:
第四步,根据基于广义相关熵的协方差矩阵计算噪声子空间并估计波达方向,包括中心波达方向和扩散角
(4.1)对第三步得到的基于广义相关熵的协方差矩阵进行奇异值分解,获得噪声子空间;
(4.2)通过步骤(4.1)中的噪声子空间计算二维空间谱函数f(ψ):
其中,Un为噪声子空间,(·)H表示共轭转置,(·)*表示共轭运算,a(θ)为点源模型下的导向向量,g(θ,ψi)表示确定角信号密度。
(4.3)对二维空间谱进行谱峰搜索,估计中心波达方向和扩散角:
在图2中,噪声环境为脉冲噪声,分布源模型下,本发明和DSPE算法进行了对比,随着广义信噪比的从-5dB增加到10dB,两种方法的均方误差都随之下降,但本发明算法性能好于DSPE算法。
在图3中,噪声环境为脉冲噪声,分布源模型下,本发明和DSPE算法进行了对比,随着alpha稳定分布的特征指数从α=1.0增加到α=2.0,两种方法的均方误差都随之下降,但本发明算法性能好于DSPE算法。
Claims (1)
1.一种无直达信号条件下的无线电波达方向估计方法,其特征在于,步骤如下:
第一步,建立分布源信号模型
等距线性阵列的M个阵元,相邻阵元间距离满足半波长条件d≤λ/2,其中,波长λ=2πc/ω,c为光的传播速度,ω为信号的中心频率;当L个互不相关的分布式信号入射到等距线阵时,接收到信号的数学模型表示为:
其中,si(θ,ψi)是第i个分布式信源的扩散信号密度,θ∈[-π/2,π/2]是角度方向;ψi是待估计的参数,包含中心波达方向和扩散角;Θ=[-π/2,π/2]为角度的观测域,n(t)=[n1(t),n2(t),…,nM(t)]T是M×1的加性测量噪声,a(θ)=[1,e-j2πd/λsinθ,…,e-j2π(M-1)d/λsinθ]T为入射信号的点源模型下的导向向量,j为复数虚部表示因子;
在分布式信源模型下,分布式信源的扩散信号密度定义为:
si(θ,ψi)=δig(θ,ψi) (2)
其中,g(θ,ψi)表示确定角信号密度,δi是一个随机变量,因此,式(1)转换为式(3)的形式为:
其中,b(ψi)为分布式信源的广义导向向量,见式(4):
b(ψi)=∫θ∈Θa(θ)g(θ,ψi)dθ (4)
在小角度展开的假设下,对于给定的角信号分布函数,广义导向向量的解析表达式为:
其中,a(θ)为点源导向向量,m表示阵列中的第m个阵元,θ′表示积分变量,θi表示点源情况下第i个的波达方向,ρ(θ′;ψi)表示确定性角信号分布函数;
分布式信源的角密度函数符合高斯分布,其概率密度函数表示为:
其中,σi为标准差;
则广义导向向量表示为:
第二步,计算广义自相关熵
定义随机变量X服从alpha稳定分布,则广义自相关熵定义为:
Gσ(X)=E[κσ(|X|-μX)] (9)
其中,E[·]为求数学期望,κσ为核函数,选取高斯核或指数核;μX表示为:
其中,N为快拍数,p为(0,2]区间的常数,xi为阵列中第i个传感器接收到的信号序列;
第三步,计算信号的广义相关熵的协方差矩阵
利用等距线性阵列的每个阵元接收到的信号x(t)=[x1(t),x2(t),…,xM(t)]T计算脉冲噪声环境下基于广义相关熵的协方差矩阵:
第四步,根据基于广义相关熵的协方差矩阵计算噪声子空间并估计波达方向,包括中心波达方向和扩散角
(4.1)对第三步得到的基于广义相关熵的协方差矩阵进行奇异值分解,获得噪声子空间;
(4.2)通过步骤(4.1)中的噪声子空间计算二维空间谱函数f(ψ):
其中,Un为噪声子空间,(·)H表示共轭转置,(·)*表示共轭运算,a(θ)为点源模型下的导向向量,g(θ,ψi)表示确定角信号密度;
(4.3)对二维空间谱进行谱峰搜索,估计中心波达方向和扩散角:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910486779.5A CN110320490B (zh) | 2019-06-05 | 2019-06-05 | 一种无直达信号条件下的无线电波达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910486779.5A CN110320490B (zh) | 2019-06-05 | 2019-06-05 | 一种无直达信号条件下的无线电波达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110320490A CN110320490A (zh) | 2019-10-11 |
CN110320490B true CN110320490B (zh) | 2020-12-08 |
Family
ID=68120746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910486779.5A Active CN110320490B (zh) | 2019-06-05 | 2019-06-05 | 一种无直达信号条件下的无线电波达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110320490B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111427006A (zh) * | 2020-04-27 | 2020-07-17 | 北京交通大学 | 一种基于相控阵的单通道空间谱测向方法 |
CN112346004A (zh) * | 2020-11-20 | 2021-02-09 | 江苏师范大学 | 一种基于score函数的广义协方差的多信号分类算法 |
CN112834981B (zh) * | 2021-03-15 | 2022-07-15 | 哈尔滨工程大学 | 一种脉冲噪声背景下的零陷阵波达方位估计方法 |
CN113255528A (zh) * | 2021-05-28 | 2021-08-13 | 胡琼 | 一种单多目标分类方法及装置 |
CN113625220A (zh) * | 2021-06-28 | 2021-11-09 | 台州学院 | 一种多径信号波达方向和扩散角快速估计新方法 |
CN113466784B (zh) * | 2021-06-28 | 2024-09-13 | 台州学院 | 一种强脉冲噪声下自适应分布源doa估计方法 |
CN113791379B (zh) * | 2021-09-23 | 2024-02-20 | 南京航空航天大学 | 嵌套阵列非高斯环境下的正交匹配追踪doa估计方法 |
CN114217264B (zh) * | 2021-11-22 | 2024-06-28 | 台州学院 | 一种基于深度学习的无线电信号波达方向估计与定位装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592721B (zh) * | 2009-06-19 | 2011-09-14 | 中国人民解放军理工大学 | 基于特征值重构的相干信号到达角估计方法 |
CN104482925A (zh) * | 2014-12-09 | 2015-04-01 | 中国海洋石油总公司 | 基于分布源模型的多波束测深系统复杂地形的测量方法 |
CN106291449B (zh) * | 2016-08-04 | 2018-10-16 | 大连大学 | 对称稳定分布噪声下波达方向角估计方法 |
CN106980106B (zh) * | 2017-04-21 | 2020-03-17 | 天津大学 | 阵元互耦下的稀疏doa估计方法 |
CN109283485A (zh) * | 2018-10-26 | 2019-01-29 | 大连理工大学 | 一种基于共变的相干分布式非圆信号角度参数联合估计方法 |
CN109738852B (zh) * | 2018-12-15 | 2022-12-16 | 华南理工大学 | 基于低秩矩阵重建的分布式源二维空间谱估计方法 |
-
2019
- 2019-06-05 CN CN201910486779.5A patent/CN110320490B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110320490A (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110320490B (zh) | 一种无直达信号条件下的无线电波达方向估计方法 | |
CN108828522B (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
CN109188344B (zh) | 脉冲噪声环境下基于互循环相关music算法信源个数与来波方向角估计方法 | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN108957387B (zh) | 一种卫星信号二维到达角估计方法及系统 | |
CN108089147B (zh) | 改进的短波单站定位方法 | |
CN103353588B (zh) | 基于天线均匀平面阵的二维波达方向角估计方法 | |
CN109946643B (zh) | 基于music求解的非圆信号波达方向角估计方法 | |
CN112612010A (zh) | 一种基于波瓣分裂预处理的米波雷达低仰角测高方法 | |
CN113189592B (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN110231589B (zh) | 一种大扩散角的多径信号波达方向估计方法 | |
CN110673119A (zh) | 基于压缩感知的非正则化方位估计方法及系统 | |
CN108828586B (zh) | 一种基于波束域的双基地mimo雷达测角优化方法 | |
CN113064147A (zh) | 一种低信噪比下新型匹配场被动定位方法 | |
CN111580042B (zh) | 一种基于相位优化的深度学习测向方法 | |
CN112255629A (zh) | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 | |
CN114239251A (zh) | 一种近端多径条件下阵列测向精度的评估方法 | |
CN114779236A (zh) | 一种基于空间平滑music改进的米波雷达低仰角测高方法 | |
CN106908754A (zh) | L型声矢量传感器阵列esprit解相干参数估计方法 | |
CN111665469B (zh) | 一种基于空间时频分布的水下多径信号参数估计方法 | |
CN113625220A (zh) | 一种多径信号波达方向和扩散角快速估计新方法 | |
CN110389319B (zh) | 一种基于低空多径情况下的mimo雷达doa估计方法 | |
CN111812607A (zh) | 基于波束空间的米波mimo雷达低仰角估计方法 | |
CN113671477B (zh) | 一种基于图信号处理的雷达目标距离估计方法 | |
CN114755628A (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 |