CN113608164A - 一种幅相误差自校正的稳健单站直接定位方法 - Google Patents
一种幅相误差自校正的稳健单站直接定位方法 Download PDFInfo
- Publication number
- CN113608164A CN113608164A CN202110595981.9A CN202110595981A CN113608164A CN 113608164 A CN113608164 A CN 113608164A CN 202110595981 A CN202110595981 A CN 202110595981A CN 113608164 A CN113608164 A CN 113608164A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- phase error
- station
- signal
- observation
- 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
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
- 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
- G01S5/0249—Determining position using measurements made by a non-stationary device other than the device whose position is being determined
-
- 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
- G01S5/0205—Details
- G01S5/021—Calibration, monitoring or correction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
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)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种幅相误差自校正的稳健单站直接定位方法,建立幅相误差下运动单站接收的信号模型,构造接收数据的协方差矩阵,对采样协方差矩阵特征分解,根据已知的幅相误差估计值估计发射器位置,根据已知的发射器位置估计值估计幅相误差矢量,多次迭代,将最终估计的辐射源位置作为最终定位结果。本发明基于特征分解子空间正交性质迭代完成目标位置和幅相误差参数的在线联合估计,推导了目标位置和幅相误差参数的闭式解形式,提高定位性能的同时降低计算复杂度,通过推导目标位置和幅相误差参数的闭式解实现了目标位置和幅相误差矢量的在线联合估计,大大提高了定位精度,提升了算法的稳健性。
Description
技术领域
本发明涉及定向定位领域,尤其是一种稳健单站直接定位方法。
背景技术
无线电定位技术在民用和军事应用中发挥着重要作用,传统的无源定位技术采用“中间参数估计+目标位置解算”两步定位模式,这种定位方法在中间参数估计过程中会损失一些目标位置信息,具有一定的局限性,在低信噪比下定位性能难以进一步提高。直接定位技术无需中间参数估计,接收的信号采用单步定位方法直接进行目标位置的估计,在低信噪比下相比两步定位方法具有估计精度高、分辨能力强、无需数据关联等优点。可移动单站直接定位系统具有机动性好、隐蔽性高、无需向中心站上传数据等优点,成为直接定位技术的研究热点。
目前基于直接定位方法的研究大多基于阵列流型精确已知的假设,因而性能优良。然而在实际应用中,阵列误差不可避免,限制了定位方法取得更好的定位性能。当观测站存在由于接收通道内放大器增益不一致引起的通道幅相误差时,阵列模型不准确会导致定位性能急剧下降。
发明内容
为了克服现有技术的不足,本发明提供一种幅相误差自校正的稳健单站直接定位方法,该方法基于特征分解子空间正交性质迭代完成目标位置和幅相误差参数的在线联合估计,推导了目标位置和幅相误差参数的闭式解形式,提高定位性能的同时降低计算复杂度。
本发明解决其技术问题所采用的技术方案包括以下步骤:
a)建立幅相误差下运动单站接收的信号模型;
b)构造接收数据的协方差矩阵:
c)对采样协方差矩阵特征分解:
e)根据已知的幅相误差估计值估计发射器位置:子空间数据融合的直接定位代价函数如下式:
f)根据已知的发射器位置估计值估计幅相误差矢量:
式中
其中αl(pq)=diag(al(pq))∈CM×M,δ=vec(Γ)∈CM×1,w=[1,0,…,0]T∈CM×1;
g)迭代输出:不断重复步骤e和步骤f,直到达到最大迭代次数为止,将最终估计的辐射源位置作为最终定位结果。
所述建立幅相误差下运动单站接收的信号模型:假定一个运动的观测站沿已知运动轨迹运动,观测站为M阵元,通道幅相误差为Γ,Q个发射器的坐标为pq,q=1,…,Q,辐射的信号sq(t)为互不相关的窄带信号,观测站在位置坐标为pl=[xl,yl]T,l=1,…,L的观测间隙进行采样,是转置符号,将数据样本总时间分为J段,每段G点采样点,Bl(p)为第l个观测间隙下时域扩展的阵列流型矢量,各阵元接收的噪声是均值为零的加性高斯白噪声,第l个观测间隙的第j段G个采样数据为:
sl(j)=[sl,1(j),…,sl,Q(j)]T (2)
Bl(p)=[bl(p1),bl(p2),…,bl(pQ)] (4)
gl(pq)=[1,exp(i2πfl(pq)Ts),…,exp(i2πfl(pq)(G-1)Ts)]T (6)
其中xl(j),sl(j),wl(j)分别是第j段时间的接收信号、发射信号和噪声数据;其中ρm表示阵列中第m个阵元的幅度误差,φm表示阵列中第m 个阵元的相位误差,m=2,…,M,sl,q(j)表示第q个发射器的基带信号包络在第l个观测间隙内的第j个采样,T为采样周期,bl(pq)为第l个观测间隙的时域扩展后的阵列流型向量,代表克罗内克积,gl(pq)为第l个观测间隙的运动单站所产生的的多普勒相移向量,rl(jG+g)表示第l个观测间隙下第j段第g个采样信号;
nl(jG+g)∈CM×1为加性高斯白噪声向量,al(pq)表示在第l个观测间隙观测站对第q个发射器的辐射信号的理想导向矢量:
dl=[dl1,…,dlM] (9)
dl表示第m个阵元在第l个观测间隙内相对于参考阵元的位置,λq表示信号的波长,π表示圆周率。
在实际应用中,数据协方差矩阵由有限次采样数据的采样协方差矩阵代替:
式(13)为采样协方差矩阵。
本发明的有益效果在于本发明结合运动单站直接定位模型实现了幅相误差下的稳健直接定位算法,通过推导目标位置和幅相误差参数的闭式解实现了目标位置和幅相误差矢量的在线联合估计,大大提高了定位精度,提升了算法的稳健性。
附图说明
图1是基于幅相误差自校正的稳健单站直接定位方法流程图。
图2是运动单站的直接定位场景示意图。
图3是观测站与发射器位置的几何分布图。
图4是在信噪比SNR=10dB,采样点数JG=120,时域扩展因子G=3,观测站的通道幅度误差为g=[1,0.25,0.5,0.75,1],相位误差为φ=[0,0.2,0.2,0.2,0.2](rad)时采用子空间数据融合法的定位结果图,其中图4(a)为定位结果图,图4(b)为定位谱图。
图5是在信噪比SNR=10dB,采样点数JG=120,时域扩展因子G=3,观测站的通道幅度误差为g=[1,0.25,0.5,0.75,1],相位误差为φ=[0,0.2,0.2,0.2,0.2](rad)时采用本发明方法的定位结果图,其中图5(a)为定位结果图,图5(b)为定位谱图.
图6是采样点数JG=120,时域扩展因子G=3,观测站的通道幅度误差为 g=[1,0.25,0.5,0.75,1],相位误差为φ=[0,0.2,0.2,0.2,0.2](rad)本发明方法在不同信噪比下均方根误差变化图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明的实施例提供一种基于幅相误差自校正的稳健单站直接定位方法,图1为该方法流程图,图2是定位模型示意图,图3是该实例的观测站与发射器的几何分布图,本实例中观测站为5阵元阵列阵元间距为半波长,发射器数目Q=1,观测站的观测间隙L=4,每个观测间隙采样点数JG=120,时域扩展因子G=3,具体实施步骤如下:
步骤一:建立图2的运动单站直接定位模型,发射器位置为p0=[-1.5,0.5]T(km),发射器辐射的信号载频fc=200MHz,观测站以v=[2000,0]T(m/s)的运动速度从pstart=[-3,-3]T(km)运动到pend=[3,-3]T(km),共L=4个观测间隙。每个观测间隙以 fs=4000(采样/秒)的采样速率进行120点数据采样,将采样点划分为40段,每段数据3 个采样点,得到第l个观测间隙的第j段G个采样数据xl(j)。
步骤二:构造接收数据的采样协方差矩阵。由120个采样数据得到第l个观测间隙的采样协方差矩阵:
步骤五:给定幅相误差估计值估计发射器位置。子空间数据融合的直接定位代价函数如下式,搜索区间步长为10米。
步骤七:重复步骤五和步骤六,多次迭代直到达到迭代次数上限10次为止,把当前的发射器位置估计值作为最终定位结果。
图3是本发明实验的观测站与发射器坐标位置分布图,即发射器位置为 p0=[-1.5,0.5]T(km),观测站位置从pstart=[-3,-3]T(km)运动到pend=[3,-3]T(km);
图4是本发明实验采用未校正算法的定位结果图,未校正算法定位坐标为 [-1.55,0.46]T(km),定位谱峰增益为-9.715dB,其中图4(a)为定位结果图,图4(b) 为定位谱图;
图5是本发明实验采用迭代信号子空间算法的定位结果图,其中图5(a)为定位结果图,图5(b)为定位谱图,该校正算法定位坐标为[-1.55,0.47]T(km),定位谱峰增益为11.36dB,相比未校正算法,定位结果更加准确,谱峰峰值增益更高。
图6是本发明实验获得的上述两种算法比较的定位误差曲线,可以看出在低信噪比下稳健直接定位算法相比未校正算法具有更高的定位精度。
综上,幅相误差自校正的稳健单站直接定位方法在定位精度和分辨率要优于未校正的直接定位算法,大大提高了定位精度,提升了算法的稳健性。
Claims (3)
1.一种幅相误差自校正的稳健单站直接定位方法,其特征在于包括下述步骤:
a)建立幅相误差下运动单站接收的信号模型;
b)构造接收数据的协方差矩阵:
c)对采样协方差矩阵特征分解:
e)根据已知的幅相误差估计值估计发射器位置:子空间数据融合的直接定位代价函数如下式:
f)根据已知的发射器位置估计值估计幅相误差矢量:
式中
其中αl(pq)=diag(al(pq))∈CM×M,δ=vec(Γ)∈CM×1,w=[1,0,…,0]T∈CM×1;
g)迭代输出:不断重复步骤e和步骤f,直到达到最大迭代次数为止,将最终估计的辐射源位置作为最终定位结果。
2.根据权利要求1所述的幅相误差自校正的稳健单站直接定位方法,其特征在于:
所述建立幅相误差下运动单站接收的信号模型:假定一个运动的观测站沿已知运动轨迹运动,观测站为M阵元,通道幅相误差为Γ,Q个发射器的坐标为pq,q=1,…,Q,辐射的信号sq(t)为互不相关的窄带信号,观测站在位置坐标为pl=[xl,yl]T,l=1,…,L的观测间隙进行采样,[·]T是转置符号,将数据样本总时间分为J段,每段G点采样点,Bl(p)为第l个观测间隙下时域扩展的阵列流型矢量,各阵元接收的噪声是均值为零的加性高斯白噪声,第l个观测间隙的第j段G个采样数据为:
sl(j)=[sl,1(j),…,sl,Q(j)]T (2)
Bl(p)=[bl(p1),bl(p2),…,bl(pQ)] (4)
gl(pq)=[1,exp(i2πfl(pq)Ts),…,exp(i2πfl(pq)(G-1)Ts)]T (6)
其中xl(j),sl(j),wl(j)分别是第j段时间的接收信号、发射信号和噪声数据;其中ρm表示阵列中第m个阵元的幅度误差,φm表示阵列中第m个阵元的相位误差,m=2,…,M,sl,q(j)表示第q个发射器的基带信号包络在第l个观测间隙内的第j个采样,T为采样周期,bl(pq)为第l个观测间隙的时域扩展后的阵列流型向量,代表克罗内克积,gl(pq)为第l个观测间隙的运动单站所产生的的多普勒相移向量,rl(jG+g)表示第l个观测间隙下第j段第g个采样信号;
nl(jG+g)∈CM×1为加性高斯白噪声向量,al(pq)表示在第l个观测间隙观测站对第q个发射器的辐射信号的理想导向矢量:
dl=[dl1,…,dlM] (9)
dl表示第m个阵元在第l个观测间隙内相对于参考阵元的位置,λq表示信号的波长,π表示圆周率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110595981.9A CN113608164B (zh) | 2021-05-29 | 2021-05-29 | 一种幅相误差自校正的稳健单站直接定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110595981.9A CN113608164B (zh) | 2021-05-29 | 2021-05-29 | 一种幅相误差自校正的稳健单站直接定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113608164A true CN113608164A (zh) | 2021-11-05 |
CN113608164B CN113608164B (zh) | 2023-06-30 |
Family
ID=78303419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110595981.9A Active CN113608164B (zh) | 2021-05-29 | 2021-05-29 | 一种幅相误差自校正的稳健单站直接定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113608164B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018049595A1 (zh) * | 2016-09-14 | 2018-03-22 | 深圳大学 | 一种基于交替方向乘子法的稳健稀疏恢复stap方法及其系统 |
CN109975749A (zh) * | 2019-02-26 | 2019-07-05 | 中国人民解放军战略支援部队信息工程大学 | 一种校正源存在条件下的短波单站直接定位方法 |
-
2021
- 2021-05-29 CN CN202110595981.9A patent/CN113608164B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018049595A1 (zh) * | 2016-09-14 | 2018-03-22 | 深圳大学 | 一种基于交替方向乘子法的稳健稀疏恢复stap方法及其系统 |
CN109975749A (zh) * | 2019-02-26 | 2019-07-05 | 中国人民解放军战略支援部队信息工程大学 | 一种校正源存在条件下的短波单站直接定位方法 |
Non-Patent Citations (1)
Title |
---|
黄秀琼;郝克钢;卢铭迪;樊荣;万群;: "基于多次观测的机载天线阵列幅相误差自校正方法", 电讯技术, no. 01 * |
Also Published As
Publication number | Publication date |
---|---|
CN113608164B (zh) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110208735B (zh) | 一种基于稀疏贝叶斯学习的相干信号doa估计方法 | |
CN109490819B (zh) | 一种基于稀疏贝叶斯学习的离格波达方向估计方法 | |
CN110687500B (zh) | 智能天线获取无线信号到达角的识别和定位方法及系统 | |
US11982698B2 (en) | Joint denoising and delay estimation for the extraction of pulse-width of signals in RF interference | |
CN112130111B (zh) | 一种大规模均匀十字阵列中单快拍二维doa估计方法 | |
CN107966676B (zh) | 复杂噪声环境下阵列天线角度和信源个数联合估计方法 | |
CN111965596A (zh) | 一种基于联合参数估计的低复杂度单锚节点定位方法及装置 | |
Wang et al. | Iterative constrained weighted least squares estimator for TDOA and FDOA positioning of multiple disjoint sources in the presence of sensor position and velocity uncertainties | |
CN111965595A (zh) | 一种基于无人机的多非圆信源高精度直接定位方法 | |
CN111199280A (zh) | 短波信道模型误差存在下联合信号复包络和载波相位信息的多站目标源地理坐标估计方法 | |
CN109696651B (zh) | 一种基于m估计的低快拍数下波达方向估计方法 | |
CN107144811A (zh) | 一种单通道接收阵列信号的循环子空间测向方法 | |
CN112187315B (zh) | 一种大规模扩频通信数字阵列同时多用户快速角度估计方法及系统 | |
CN113608164A (zh) | 一种幅相误差自校正的稳健单站直接定位方法 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
CN115407266A (zh) | 一种基于互谱子空间正交性的直接定位方法 | |
CN113589223B (zh) | 基于互耦情况下嵌套阵列的测向方法 | |
CN112881971B (zh) | 电磁方向性互耦效应下相干干扰源测向方法 | |
CN113093098B (zh) | 基于lp范数补偿的轴向不一致矢量水听器阵列测向方法 | |
CN112415469B (zh) | 一种两维数字阵列雷达快速干扰测向方法 | |
CN115079090A (zh) | 基于esprit与加权降维搜索的多阵列非圆源快速定位方法 | |
CN100399730C (zh) | 一种阵列天线通道误差的盲估计方法 | |
CN108833320B (zh) | 一种基于Keystone变换的宽带DS/FH信号参数估计方法及系统 | |
CN113093093A (zh) | 基于线性阵列波达方向估计的车辆定位方法 | |
CN112579972A (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 |