CN108957387B - 一种卫星信号二维到达角估计方法及系统 - Google Patents
一种卫星信号二维到达角估计方法及系统 Download PDFInfo
- Publication number
- CN108957387B CN108957387B CN201810488803.4A CN201810488803A CN108957387B CN 108957387 B CN108957387 B CN 108957387B CN 201810488803 A CN201810488803 A CN 201810488803A CN 108957387 B CN108957387 B CN 108957387B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- signal
- arrival angle
- clutter
- obtaining
- 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
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
-
- 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/50—Systems of measurement based on relative movement of target
-
- 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/66—Radar-tracking systems; Analogous systems
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
本发明属于测定方向或测定对预定方向的偏离的系统技术领域,公开了一种卫星信号二维到达角估计方法及系统,利用分数低阶循环相关对非高斯杂波进行初步抑制;利用改进的零阶统计量对非高斯杂波进行二次抑制;通过最优化方法得到信号子空间,并利用多重信号聚类方法得到二维到达角估计的三维图;通过对三维图进行最大值搜索从而得到非高斯杂波场景下二维到达角的估计值。当广义信噪比大于10dB时,方位角和俯仰角估计的均方根误差小于10‑1。由此可见,本发明在低信噪比条件下,对非高斯杂波场景下卫星信号二维到达角具有良好的估计性能;利用本发明可以通过多个卫星辐射源信号对运动目标进行跟踪定位。
Description
技术领域
本发明属于测定方向或测定对预定方向的偏离的系统技术领域,尤其涉及一种卫星信号二维到达角估计方法及系统。
背景技术
目前,业内常用的现有技术是这样的:到达角估计作为阵列信号处理领域的重要研究方向,它是根据阵列天线接收的目标信号入射角来分析目标方位的一种技术,目前广泛应用于定位、通信和预测等各种军事和民用领域。到达角估计的准确程度将直接影响对运动目标的定位精度,所以到达角估计成为运动目标定位的重要环节之一。通常情况下信号处理的技术发展和理论研究都是以高斯分布杂波为背景假设的,因为高斯分布杂波模型比其他杂波模型简单,便于进行理论分析、算法研究和功能实现,这使得基于高斯杂波模型下的到达角估计理论得到飞速发展,相关技术研究也趋于成熟。但是,由于电磁传播环境变得越来越复杂,实际的信号传输过程中往往存在着一类具有脉冲特性的非高斯杂波,它们相对于高斯杂波具有更大幅度的数据突变,这类噪声被称为非高斯杂波,其分布特性符合alpha稳定分布杂波模型。由于非高斯杂波没有二阶及二阶以上的统计量以至于传统在高斯杂波模型下的到达角估计算法在这种环境下将不再适用,所以非高斯杂波模型下的到达角估计逐步成为一个研究热点。为了进一步提高非高斯杂波模型下的到达角估计的准确性和可靠性,研究非高斯杂波场景下基于多卫星辐射源的回波信号二维到达角估计方法具有研究意义和使用价值。现有技术一利用alpha稳定分布杂波下随机变量相关熵的有界性,采用广义的相关熵作为非线性变换,实现了alpha稳定杂波下的DOA估计,虽然克服了构造分数低阶统计量时较依赖噪声特征指数信息的缺点,但是忽略了信号的时间特性。现有技术二提出一种利用ROOT–MUSIC来对DOA进行估计的算法,根据多项式进行求根,从而避免了谱峰搜索,提供更高的分辨能力,可以应用于严重非高斯杂波下的DOA估计,但是该算法只适用于一维DOA估计,不适合二维DOA估计的场景。现有技术三有效地利用了信号协方差矩阵的稀疏模型,提出一种新的算法对信号的DOA进行估计,并提供了正则化参数的解析解,但计算复杂度较高。现有技术四利用信号子阵的协方差矩阵提出了一种alpha稳定分布杂波场景下信号二维DOA的估计方法,但是在噪声的脉冲特性较强时估计性能不好。
综上所述,现有技术存在的问题是:现有的到达角估计方法存在忽略了信号的时间特性,造成DOA估计精度不高;只适用于一维DOA估计,不能用于对二维DOA进行估计的场景;运算量大,导致计算复杂度较高;在噪声特性较强时估计精度不高。
发明内容
针对现有技术存在的问题,本发明提供了一种卫星信号二维到达角估计方法及系统。
本发明是这样实现的,一种卫星信号二维到达角估计方法,所述卫星信号二维到达角估计方法利用分数低阶循环相关对非高斯杂波进行初步抑制;利用改进的零阶统计量对非高斯杂波进行二次抑制;通过最优化方法得到信号子空间,并利用多重信号聚类方法得到二维到达角估计的三维图;通过对三维图进行最大值搜索从而得到非高斯杂波场景下二维到达角的估计值。
进一步,所述卫星信号二维到达角估计方法包括以下步骤:
步骤二,利用改进的零阶统计量对非高斯杂波进行二次抑制;
步骤三,通过最优化方法得到信号子空间;
X=AS+N;
其中X为天线阵列接收到的含杂波信号,S为发射的卫星信号,A为信号S的导向矩阵,N为加性alpha稳定分布杂波;
则接收信号X的分数低阶循环相关矩阵表示为:
式中,p为分数低阶循环相关的阶数,其中(·)<p>=|·|p-1(·)*,ε为循环频率。
进一步,所述步骤二,利用改进的零阶统计量对非高斯杂波进行二次抑制得到基于改进零阶统计量的分数低阶循环相关SIZOS的方法为:
得到卫星信号的基于改进零阶统计量的分数低阶循环相关SIZOS。
进一步,所述步骤三,通过最优化方法得到信号子空间的方法为:第n次迭代的结果为Y(n)和Z(n),其中Y(0)、Z(0)可分别初始化为随机的列满秩和行满秩矩阵,则第n+1次的迭代结果为:
把第n+1次的迭代结果代入代价函数得到:
根据梯度下降法,f(z)的最小值通过以下公式迭代得到:
z(n+1)=z(n)+μ(n)▽f(z(n));
最佳步长由下式得到:
本发明的另一目的在于提供一种实现所述卫星信号二维到达角估计方法的卫星信号二维到达角估计系统,所述卫星信号二维到达角估计系统包括:
二次抑制模块,用于利用改进的零阶统计量对非高斯杂波进行二次抑制;
三维图获取模块,用于通过最优化方法得到信号子空间,并利用多重信号聚类方法得到二维到达角估计的三维图;
估计值计算模块,用于通过对三维图进行最大值搜索从而得到非高斯杂波场景下二维到达角的估计值。
本发明的另一目的在于提供一种实现所述卫星信号二维到达角估计方法的卫星通信系统。
综上所述,本发明的优点及积极效果为:当广义信噪比大于10dB时,方位角和俯仰角估计的均方根误差小于10-1。由此可见,本发明在低信噪比条件下,对非高斯杂波场景下卫星信号二维到达角具有良好的估计性能;利用本发明可以通过多个卫星辐射源信号对运动目标进行跟踪定位。
附图说明
图1是本发明实施例提供的卫星信号二维到达角估计方法流程图。
图2是本发明实施例提供的在非高斯杂波场景下卫星信号二维到达角估计性能示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种非高斯杂波场景下卫星信号二维到达角估计方法,旨在解决非高斯相干杂波场景下卫星信号二维到达角难以估计的问题。
如图1所示,本发明实施例提供的卫星信号二维到达角估计方法包括以下步骤:
S101:利用分数低阶循环相关对非高斯杂波进行初步抑制;利用改进的零阶统计量对非高斯杂波进行二次抑制;
S102:通过最优化方法得到信号子空间,并利用多重信号聚类方法得到二维到达角估计的三维图;
S103:通过对三维图进行最大值搜索从而得到非高斯杂波场景下二维到达角的估计值。
本发明实施例提供的卫星信号二维到达角估计方法具体包括以下步骤:
X=AS+N;
其中X为天线阵列接收到的含杂波信号,S为发射的卫星信号,A为信号S的导向矩阵,N为加性alpha稳定分布杂波。
则接收信号X的分数低阶循环相关矩阵可以表示为:
上式中,p为分数低阶循环相关的阶数,其中(·)<p>=|·|p-1(·)*,ε为循环频率。并且从上式可以看出对信号求分数低阶循环相关时,进行了一次共轭运算,所以信号原来的相位已经被抵消,也就是经过上式后,信号已经不包含原来的相位信息,所以后续的非线性变换不会改变信号原来固有的相位信息。
步骤二,利用改进的零阶统计量对非高斯杂波进行二次抑制得到基于改进零阶统计量的分数低阶循环相关SIZOS的方法为:
得到卫星信号的基于改进零阶统计量的分数低阶循环相关SIZOS。
步骤三,通过最优化方法得到信号子空间的方法为:
如果利用SIZOS计算信号子空间和杂波子空间,根据参数估计理论,当得到最优信号子空间和最优杂波子空间时,包含最优子空间的统计量和原统计量之间存在下列代价函数所述关系:
可以通过矩阵奇异值分解的方式对上式求解,具体表示为:
其中,SVD(·)表示奇异值分解算子,DS和Dn分别表示通过对SIZOS进行奇异值分解后由K个较大奇异值和(MN-K)个较小奇异值构成的对角矩阵,即代表信号子空间的奇异值和杂波子空间的奇异值所构成的对角矩阵。US和VS分别代表对角矩阵DS的左特征矩阵和右特征矩阵。同样地,Un和Vn分别代表对角矩阵Dn的左特征矩阵和右特征矩阵。从上式可知代价函数的最优解为:
根据矩阵的相关理论可知,US张成的线性矩阵即为信号子空间;
为了方便,把最优解式的形式重新写为:
为了求信号子空间,只需要取使上式所示代价函数最小时的Y即可。
把SIZOS代入上式可得到:
则代价函数可以重新写为:
根据最优化准则,通过最小化上式就可以实现Y、Z的估计。本文采用最优化算法实现上式中对未知参量的求解,假设第n次迭代的结果可以写为Y(n)和Z(n),其中Y(0)、Z(0)可分别初始化为随机的列满秩和行满秩矩阵,则第n+1次的迭代结果可以写为:
从上式可以看出,若采用上式的迭代算法,J(Y,Z)算法是收敛的,即:
J(Y(n),Z(n))>J(Y,Z(n+1))>J(Y(n+1),Z(n+1));
把第n+1次的迭代结果代入代价函数便可以得到:
综上可以看出,Z(n+1)和Y(n+1)的求解形式相同,所以将对Z(n+1)进行求解,Y(n+1)的求解方式与式Z(n+1)相同。
从上面分析可知,式Z(n+1)的求解可以转换为以下T个子问题进行求解:
根据梯度下降法,f(z)的最小值可以通过以下公式迭代得到:
z(n+1)=z(n)+μ(n)▽f(z(n));
而最佳步长可以由下式得到:
当得到信号子空间Y后,则可以得到DOA的估计的三维图为:
下面结合仿真实验对本发明的应用效果作详细的描述。
为了测试本方面的检验统计量的性能,参数设置如下:仿真实验采用的信号为GPS、DVB-S和INMARSAT三个卫星信号,它们的载频分别是:fG=1.57GHz、fD=12.38GHz、fI=4.2GHz;TDOA分别是1μs、2μs、3μs;FDOA分别是100Hz、200Hz、300Hz;杂波采用标准alpha稳定分布杂波,其参数为α=1.5,γ=1。接收天线采用5×5的均匀面阵,阵元间距为DVB-S卫星信号频率对应波长的一半。因为三个回波入射的方向一致,所以DOA都设置为(40°,60°),快拍数为600。分数低阶循环相关中参数p=1。每个实验条件下进行1000次Monte-Carlo实验。采用的评估标准是参数估计的均方根误差(RMSE),其具体定义如下:
仿真结果如图2所示,当广义信噪比大于10dB时,方位角和俯仰角估计的均方根误差小于10-1,说明本发明的非高斯杂波场景下卫星信号二维到达角估计方法是有效可行的。由此说明本发明方法在低广义信噪比条件下,对非高斯杂波场景下卫星信号二维到达角具有良好的估计性能。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种卫星信号二维到达角估计方法,其特征在于,所述卫星信号二维到达角估计方法利用分数低阶循环相关对非高斯杂波进行初步抑制;利用改进的零阶统计量对非高斯杂波进行二次抑制;通过最优化方法得到信号子空间,并利用多重信号聚类方法得到二维到达角估计的三维图;通过对三维图进行最大值搜索从而得到非高斯杂波场景下二维到达角的估计值;
所述卫星信号二维到达角估计方法包括以下步骤:
步骤二,利用改进的零阶统计量对非高斯杂波进行二次抑制;
步骤三,通过最优化方法得到信号子空间;通过最优化方法得到信号子空间的方法为:第n次迭代的结果为Y(n)和Z(n),其中Y(0)、Z(0)可分别初始化为随机的列满秩和行满秩矩阵,则第n+1次的迭代结果为:
把第n+1次的迭代结果代入代价函数得到:
根据梯度下降法,f(z)的最小值通过以下公式迭代得到:
最佳步长由下式得到:
6.一种实现权利要求1所述卫星信号二维到达角估计方法的卫星信号二维到达角估计系统,其特征在于,所述卫星信号二维到达角估计系统包括:
二次抑制模块,用于利用改进的零阶统计量对非高斯杂波进行二次抑制;
三维图获取模块,用于通过最优化方法得到信号子空间,并利用多重信号聚类方法得到二维到达角估计的三维图;
估计值计算模块,用于通过对三维图进行最大值搜索从而得到非高斯杂波场景下二维到达角的估计值。
7.一种实现权利要求1~5任意一项所述卫星信号二维到达角估计方法的卫星通信系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810488803.4A CN108957387B (zh) | 2018-05-21 | 2018-05-21 | 一种卫星信号二维到达角估计方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810488803.4A CN108957387B (zh) | 2018-05-21 | 2018-05-21 | 一种卫星信号二维到达角估计方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108957387A CN108957387A (zh) | 2018-12-07 |
CN108957387B true CN108957387B (zh) | 2022-06-21 |
Family
ID=64499308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810488803.4A Active CN108957387B (zh) | 2018-05-21 | 2018-05-21 | 一种卫星信号二维到达角估计方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108957387B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110703145B (zh) * | 2019-09-20 | 2021-09-21 | 广东石油化工学院 | 一种利用多优化理论的变压器振声信号重构方法和系统 |
CN110967665A (zh) * | 2019-10-07 | 2020-04-07 | 西安电子科技大学 | 一种多个外辐射源下运动目标回波的doa估计方法 |
CN111093265B (zh) * | 2019-12-04 | 2020-12-11 | 北京邮电大学 | 一种基于到达角度测距的协作定位方法及装置 |
CN111308530B (zh) * | 2020-02-17 | 2021-12-03 | 中国人民解放军战略支援部队信息工程大学 | 一种基于二维波达方向的短波多站和单星协同直接定位方法 |
CN111551896B (zh) * | 2020-04-25 | 2021-01-26 | 中国人民解放军战略支援部队信息工程大学 | 抑制传感器位置速度先验误差的加权多维标度toa和foa多源协同定位方法 |
CN112731486B (zh) * | 2020-12-24 | 2023-07-04 | 深圳大学 | 信号角度和信号频率的估计方法、装置、设备及存储介质 |
CN118011313A (zh) * | 2024-04-09 | 2024-05-10 | 中国航天科工集团八五一一研究所 | 一种基于高斯混合模型的脉冲测向信息分选方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1046608A (zh) * | 1989-04-20 | 1990-10-31 | 清华大学 | 基于最大熵谱估计的自适应动目标检测器 |
CN105589056A (zh) * | 2015-12-15 | 2016-05-18 | 吉林大学 | 一种多目标远近场混合源定位方法 |
CN107135176A (zh) * | 2017-07-06 | 2017-09-05 | 电子科技大学 | 基于分数低阶循环谱的图域通信信号调制识别方法 |
CN107167785A (zh) * | 2017-07-11 | 2017-09-15 | 吉林大学 | 一种稳健的大阵列mimo雷达目标收发角联合估计方法 |
CN107907854A (zh) * | 2017-10-31 | 2018-04-13 | 大连大学 | 脉冲噪声环境下的doa估计新方法 |
-
2018
- 2018-05-21 CN CN201810488803.4A patent/CN108957387B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1046608A (zh) * | 1989-04-20 | 1990-10-31 | 清华大学 | 基于最大熵谱估计的自适应动目标检测器 |
CN105589056A (zh) * | 2015-12-15 | 2016-05-18 | 吉林大学 | 一种多目标远近场混合源定位方法 |
CN107135176A (zh) * | 2017-07-06 | 2017-09-05 | 电子科技大学 | 基于分数低阶循环谱的图域通信信号调制识别方法 |
CN107167785A (zh) * | 2017-07-11 | 2017-09-15 | 吉林大学 | 一种稳健的大阵列mimo雷达目标收发角联合估计方法 |
CN107907854A (zh) * | 2017-10-31 | 2018-04-13 | 大连大学 | 脉冲噪声环境下的doa估计新方法 |
Non-Patent Citations (1)
Title |
---|
基于分数低阶统计量的DOA和TDOA估计算法研究;刘成材;《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》;20061015(第10期);正文第2.3、3.1.2、4.1-4.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN108957387A (zh) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108957387B (zh) | 一种卫星信号二维到达角估计方法及系统 | |
WO2018094565A1 (zh) | 脉冲噪声下的波束成形方法及装置 | |
Yang et al. | Knowledge‐aided STAP with sparse‐recovery by exploiting spatio‐temporal sparsity | |
CN110320490B (zh) | 一种无直达信号条件下的无线电波达方向估计方法 | |
CN107121665B (zh) | 一种基于稀疏阵的近场相干源的无源定位方法 | |
Lin et al. | Long-time coherent integration for maneuvering target detection based on ITRT-MRFT | |
CN112379327A (zh) | 一种基于秩损估计的二维doa估计与互耦校正方法 | |
CN104749555A (zh) | 一种组合相位差测向和空间谱测向的测向定位系统 | |
CN115856767B (zh) | 一种可重构智能超表面辅助的波到达方向估计方法 | |
CN107991659A (zh) | 基于字典学习的米波雷达低仰角目标测高方法 | |
Chen et al. | Iterative subspace alternating projection method for GNSS multipath DOA estimation | |
CN113567913B (zh) | 基于迭代重加权可降维的二维平面doa估计方法 | |
CN112946577B (zh) | 一种基于宽带压缩感知的超短基线水下声源定位方法 | |
CN107450046A (zh) | 低仰角多径环境下的波达角估计方法 | |
CN110231589B (zh) | 一种大扩散角的多径信号波达方向估计方法 | |
Li et al. | KT and azimuth sub‐region deramp‐based high‐squint SAR imaging algorithm mounted on manoeuvring platforms | |
CN109507634B (zh) | 一种任意传感器阵列下的基于传播算子的盲远场信号波达方向估计方法 | |
Shevgunov et al. | Heuristic global optimization algorithm for estimating a radio source location by a passive radar system | |
CN112881971B (zh) | 电磁方向性互耦效应下相干干扰源测向方法 | |
Sun et al. | Through‐the‐wall radar imaging algorithm for moving target under wall parameter uncertainties | |
Zhang et al. | DOA estimation method for underwater acoustic signals based on two-dimensional power distribution (TPD) for few element array | |
Zhang et al. | Super-resolution Localization of Moving Target Based on Improved Compressive Sensing Algorithm | |
CN112698263A (zh) | 一种基于正交传播算子的单基地互质mimo阵列doa估计算法 | |
CN112579972A (zh) | 方向性电磁耦合效应下空域信息联合估计方法 | |
Jin et al. | Radar target detection via GAMP: A sparse recovery strategy Off the grid |
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 |