CN108279411B - 一种基于mds的被动mimo时差定位方法 - Google Patents
一种基于mds的被动mimo时差定位方法 Download PDFInfo
- Publication number
- CN108279411B CN108279411B CN201810100693.XA CN201810100693A CN108279411B CN 108279411 B CN108279411 B CN 108279411B CN 201810100693 A CN201810100693 A CN 201810100693A CN 108279411 B CN108279411 B CN 108279411B
- Authority
- CN
- China
- Prior art keywords
- time difference
- positioning
- target
- vector
- measurement
- 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
- 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/06—Systems determining position data of a target
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于电子对抗技术领域,具体涉及一种基于MDS的被动MIMO时差定位方法。本发明通过在被动MIMO雷达定位系统下,将时差定位中的距离和作为一个纯虚数维度引入到三维位置空间,来构造一个含有纯虚维度的特殊四维空间,从而实现时差定位问题的多维标度分析,得到目标位置估计的加权最小二乘解,最后再通过解相关技术,进一步优化对目标的位置估计。本发明的方法实现了MIMO雷达系统下时差定位问题的多维标度分析,对于近场/远场辐射源目标均可实现精确定位,定位解算方法计算量小,具有良好的稳健性和准确性。
Description
技术领域
本发明属于电子对抗技术领域,具体涉及一种基于MDS的被动MIMO时差定位方法。
背景技术
MIMO(Multiple-Input Multiple-Output)雷达是多输入多输出雷达的简称。由于MIMO技术在移动通信技术的研究取得实际进展,在本世纪初,人们开始将MIMO技术拓展应用于雷达探测技术领域,用于目标检测,定位跟踪等。MIMO雷达系统较传统雷达在信号检测能力、参数估计精度、目标分辨率等方面有明显优势,尤其在抗雷达信号截获、强杂波条件下的弱目标检测及慢速目标检测跟踪等方面的性能较传统雷达有显著改善。
被动MIMO雷达系统是指利用已广泛存在的多基站外辐射源如:调频广播信号,电视信号,手机基站信号等作为照射源,在接收端接收经目标反射的信号,从而基于时差定位原理,获得距离和测量值,实现目标检测与定位跟踪。近年来,随着电子对抗技术的迅猛发展,传统雷达的生存能力和战场性能受到极大威胁,被动MIMO雷达系统因其隐蔽性好,定位构型灵活而受到各国军方的广泛关注。
多维标度(MDS:multidimensional scaling)是现代统计学中的重要分析方法,它是用来表示在低维空间中由任意两点的空间距离来刻画对象之间的相似度量。多维标度分析方法在现代信息技术中的无线传感器网络(WSN:wireless sensor network)节点定位、蜂窝无线通信网络移动台定位及位置校正等领域展现出良好的定位稳健性。因此,将多维标度分析方法用于被动MIMO雷达系统中进行目标定位,能够有效提升雷达定位系统的抗电磁打击能力及定位稳健性能,具有很强的工程实践价值。
发明内容
本发明的目的,就是针对上述问题,提出了一种基于MDS的被动MIMO时差定位方法,该方法通过在被动MIMO雷达定位系统下,将时差定位中的距离和作为一个纯虚数维度引入到三维位置空间,来构造一个含有纯虚维度的特殊四维空间,从而实现时差定位问题的多维标度分析,得到目标位置估计的加权最小二乘解,最后再通过解相关技术,进一步优化对目标的位置估计。
本发明所采用的技术方案为:
本发明的方法在MIMO雷达定位系统下,假定MIMO雷达系统有M个发射站和N个接收站,其站址坐标分别表示为k=1,2,...,M和l=1,2,...,N,所要定位的目标辐射源的位置为u=[x,y,z]T,定位模型如图1所示。在本系统中,真实值用 等表示,测量值用τkl,dkl等表示,测量噪声用等表示。
则时差测量值可表示为:
将其转化为距离和表示为:
则,记获得的距离和测量值向量为:
d=[d11 d12 ... d1N d21 ... dM1]T
故,系统测量方程可表示为:
本发明主要包括以下步骤:
a、构造一个以距离和为纯虚维度的特殊四维空间,建立时差定位问题的多维标度分析模型;
b、通过多维标度分析原理,建立线性方程,得到目标向量的最小二乘估计。
c、进一步,得到存在时差测量误差条件下的目标向量的加权最小二乘估计。
d、通过对目标向量进行解相关分析计算,获得目标辐射源的位置估计。
算法流程图如图2所示。
具体的,所述步骤a中,本发明基于以下原理:
根据上述距离和测量方程,可得:
则,建立多维标度分析中的内积矩阵B为:
B=ZZT
将B分块表示为:
其中具体表达式为:
上述Amn,Cmn,Dmn分别为矩阵A,C,D中位置下标为(m,n)的元素,下同。
进一步的,在所述步骤b中:
根据多维标度分析理论,可建立如下线性方程:
式中,F=Ψ(ΨTΨ)-1,其中:
定义:
则上述线性方程可重写为:
Gx=g
进而,得到目标向量的最小二乘解:
则最终目标向量的最小二乘解为:
如上所述,在时差定位系统中,存在时差测量误差,则需进一步实现存在时差测量误差条件下的多维标度分析,得到目标向量的加权最小二乘解。
更进一步的,在所述步骤c中:
测量误差向量ε可表示为:
其中:
其中:
则:
E(ε)=0
式中,Qn为测量噪声n的协方差矩阵,具体表达式为:
则,最终可得目标向量的加权最小二乘解为:
最后,在所述步骤d中:
由上述可知,系统测量方程可表示为:
距离和测量值向量为:d=[d11 d12 ... d1N d21 ... dM1]T
则通过求解目标函数极小值,可得目标位置解相关后的加权最小二乘解:
本发明的有益效果为,本发明的方法实现了MIMO雷达系统下时差定位问题的多维标度分析,对于近场/远场辐射源目标均可实现精确定位,定位解算方法计算量小,具有良好的稳健性和准确性。
附图说明
图1为MIMO雷达系统下的时差定位模型图;
图2为定位算法流程图;
图3为近场辐射源目标下定位算法的均方误差随观测误差的变化曲线;
图4为远场辐射源目标下定位算法的均方误差随观测误差的变化曲线;
具体实施方式
下面结合附图对上述基于MDS的被动MIMO时差定位方法进行仿真验证说明,首先对系统模型作如下合理假定:
1.为了便于仿真,假定MIMO雷达系统为3发3收模型,即由3个发射站和3个接收站构成MIMO雷达系统。
2.假定测量误差服从均值为零的高斯分布,且误差之间相互独立。
如上所述,假定3个发射站的位置坐标分别为t1=[4300 2500 50]T,t2=[-43002500 150]T,t3=[0 -5000 100]T,单位均为m,下同;3个接收站的位置坐标分别为r1=[05000 200]T,r2=[-4300 -2500 100]T,r3=[4300 -2500 150]T。
(1)近场辐射源目标时定位效果:
如图2所示,假定近场辐射源目标位置为u=[300 800 2000]T,采用上述的定位解算方法在时差测量误差变化下进行Monte Carlo仿真,并将解相关前后的定位均方误差进行了对比。
从图2中可以看出,本文提出的方法对于近场辐射源目标可以进行有效的定位解算,随着测量误差的增大,定位的均方误差平稳增加,表现出良好的稳健性和准确性。同时由图中对比可以看出,基于MDS的时差定位估计经过解相关之后定位效果又进一步得到提升。
(2)远场辐射源目标时定位效果:
如图3所示,假定近场辐射源目标位置为u=[4000 6000 2000]T,采用上述的定位解算方法在时差测量误差变化下进行Monte Carlo仿真,并将解相关前后的定位均方误差进行了对比。
从图3中可以看出,本文提出的方法对于远场辐射源目标可以进行有效的定位解算,随着测量误差的增大,定位的均方误差平稳增加表现出良好的稳健性和准确性。同时由图中对比可以看出,基于MDS的时差定位估计经过解相关之后定位效果又进一步得到提升。
如(1),(2)中所述,本文提出的方法对于近场/远场辐射源目标均可以实现有效的定位解算,同时针对近场/远场辐射源目标的定位效果来看,近场辐射源目标的定位效果要稍优于远场辐射源目标的定位效果。
Claims (3)
1.一种基于MDS的时差定位方法,该方法用于被动MIMO雷达系统,所述被动MIMO雷达系统有M个发射站和N个接收站,其站址坐标分别表示为:发射站为k=1,2,...,M和接收站为并设定所要定位的目标辐射源的位置为u=[x,y,z]T;将时差测量值表示为:
将其转化为距离和表示为:
记获得的距离和测量值向量为:
d=[d11 d12...d1N d21...dM1]T
则系统测量方程为:
其特征在于,所述时差定位方法包括以下步骤:
a、构造一个以距离和为纯虚维度的四维空间,建立时差定位问题的多维标度分析模型;具体方法为:
根据距离和测量方程,可得:
则,建立多维标度分析中的内积矩阵B为:
B=ZZT
将B分块表示为:
其中:
上述Amn,Cmn,Dmn分别为矩阵A,C,D中位置下标为(m,n)的元素;
b、通过多维标度分析原理,建立线性方程,得到目标向量的最小二乘估计;具体方法为:
根据多维标度分析理论,建立如下线性方程:
式中,F=Ψ(ΨTΨ)-1,其中:
定义:
则上述线性方程可重写为:
Gx=g
得到目标向量的最小二乘解:
则最终目标向量的最小二乘解为:
c、根据步骤b的结果,获取存在时差测量误差条件下的目标向量的加权最小二乘估计;
d、通过对目标向量进行解相关分析计算,获得目标辐射源的位置估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810100693.XA CN108279411B (zh) | 2018-02-01 | 2018-02-01 | 一种基于mds的被动mimo时差定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810100693.XA CN108279411B (zh) | 2018-02-01 | 2018-02-01 | 一种基于mds的被动mimo时差定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108279411A CN108279411A (zh) | 2018-07-13 |
CN108279411B true CN108279411B (zh) | 2020-04-14 |
Family
ID=62807403
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810100693.XA Active CN108279411B (zh) | 2018-02-01 | 2018-02-01 | 一种基于mds的被动mimo时差定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108279411B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111239718B (zh) * | 2020-01-17 | 2022-11-15 | 电子科技大学 | 一种基于单星辐射源的多基站目标定位方法 |
CN111352087B (zh) * | 2020-03-25 | 2022-05-24 | 电子科技大学 | 基于dbscan的被动mimo雷达多目标定位方法 |
CN113589269B (zh) * | 2021-07-14 | 2023-07-21 | 太原理工大学 | 基于线性融合的被动多基站雷达目标检测方法 |
CN113484854B (zh) * | 2021-07-21 | 2023-04-11 | 电子科技大学 | 一种外辐射源位置未知的目标定位方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100842817B1 (ko) * | 2007-01-17 | 2008-07-01 | 국방과학연구소 | Mdsa를 이용한 저 연산량 mimo 수신방식 |
CN102231912A (zh) * | 2011-07-29 | 2011-11-02 | 杭州电子科技大学 | 一种基于rssi测距的室内无线传感器网络定位方法 |
JP2013542640A (ja) * | 2010-09-20 | 2013-11-21 | インテル コーポレイション | 通信方法、装置及び製品 |
CN104902565A (zh) * | 2015-06-04 | 2015-09-09 | 杭州电子科技大学 | 一种分布式的无线传感器网络三维mds定位方法 |
CN105158730A (zh) * | 2015-08-14 | 2015-12-16 | 上海交通大学 | 基于mds子空间第4和第5特征向量的tdoa定位方法 |
CN105259533A (zh) * | 2015-10-28 | 2016-01-20 | 上海交通大学 | 基于多维标度法子空间分析的三阶段到达时间差定位方法 |
CN105891776A (zh) * | 2016-04-06 | 2016-08-24 | 上海交通大学 | 基于mds模型的直接法到达时间差定位方法 |
CN106255203A (zh) * | 2016-09-19 | 2016-12-21 | 哈尔滨工业大学 | 基于mds的终端rsrp差异补偿的定位方法 |
CN107426817A (zh) * | 2017-08-28 | 2017-12-01 | 西安电子科技大学 | 距离修正混沌粒子群多维标度无线传感器网络定位方法 |
-
2018
- 2018-02-01 CN CN201810100693.XA patent/CN108279411B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100842817B1 (ko) * | 2007-01-17 | 2008-07-01 | 국방과학연구소 | Mdsa를 이용한 저 연산량 mimo 수신방식 |
JP2013542640A (ja) * | 2010-09-20 | 2013-11-21 | インテル コーポレイション | 通信方法、装置及び製品 |
CN102231912A (zh) * | 2011-07-29 | 2011-11-02 | 杭州电子科技大学 | 一种基于rssi测距的室内无线传感器网络定位方法 |
CN104902565A (zh) * | 2015-06-04 | 2015-09-09 | 杭州电子科技大学 | 一种分布式的无线传感器网络三维mds定位方法 |
CN105158730A (zh) * | 2015-08-14 | 2015-12-16 | 上海交通大学 | 基于mds子空间第4和第5特征向量的tdoa定位方法 |
CN105259533A (zh) * | 2015-10-28 | 2016-01-20 | 上海交通大学 | 基于多维标度法子空间分析的三阶段到达时间差定位方法 |
CN105891776A (zh) * | 2016-04-06 | 2016-08-24 | 上海交通大学 | 基于mds模型的直接法到达时间差定位方法 |
CN106255203A (zh) * | 2016-09-19 | 2016-12-21 | 哈尔滨工业大学 | 基于mds的终端rsrp差异补偿的定位方法 |
CN107426817A (zh) * | 2017-08-28 | 2017-12-01 | 西安电子科技大学 | 距离修正混沌粒子群多维标度无线传感器网络定位方法 |
Non-Patent Citations (3)
Title |
---|
"一种多维尺度分析到达时间差定位算法";王琳等;《导航定位学报》;20141231;第2卷(第4期);第58-73页 * |
"观测站存在位置误差条件下基于MDS的多站时差定位算法";吴魏等;《信号处理》;20150731;第31卷(第7期);第770-776页 * |
20141231;H-W Wei等;《IET signal processing》;20081231;第2卷(第4期);第415-423页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108279411A (zh) | 2018-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108279411B (zh) | 一种基于mds的被动mimo时差定位方法 | |
Ye et al. | Neural-network-assisted UE localization using radio-channel fingerprints in LTE networks | |
Zhao et al. | Similarity analysis-based indoor localization algorithm with backscatter information of passive UHF RFID tags | |
CN104977558B (zh) | 一种基于贝叶斯压缩感知的分布源中心波达方向估计方法 | |
Noroozi et al. | Closed-form solution for elliptic localization in distributed MIMO radar systems with minimum number of sensors | |
Fang et al. | Calibration-free approaches for robust Wi-Fi positioning against device diversity: A performance comparison | |
CN109151714A (zh) | 一种三维鲁棒估计定位方法 | |
CN108387876A (zh) | 基于ctls的外辐射源雷达网双基距误差配准方法 | |
CN110954865A (zh) | 一种基于电离层信息的短波时差定位方法 | |
CN106353720A (zh) | 基于tdoa/groa的多站连续定位模型 | |
Chen et al. | Improved two-step weighted least squares algorithm for TDOA-based source localization | |
Khan et al. | Enhanced hybrid positioning in wireless networks I: AoA-ToA | |
CN107613458B (zh) | 一种tdoa条件下最优联合时间同步与定位的定位方法 | |
Song et al. | A novel source localization scheme based on unitary esprit and city electronic maps in urban environments | |
CN112255629A (zh) | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 | |
CN115150744A (zh) | 一种大型会议场馆室内信号干扰源定位方法 | |
CN107592654B (zh) | 一种基于压缩感知的同频多辐射源场强定位方法 | |
CN106908760B (zh) | 基于阵列自相关矩阵的单站无源定位方法 | |
CN111505566B (zh) | 一种特高频射频信号doa估计方法 | |
CN107144815B (zh) | 一种基于一维测向的三维定位方法 | |
Wang et al. | Robust Wi-Fi location fingerprinting against device diversity based on spatial mean normalization | |
CN116930963A (zh) | 一种基于无线通信系统的穿墙成像方法 | |
CN107682822A (zh) | 一种基于电磁场强度的压缩感知室外定位方法 | |
CN110536410B (zh) | 非视距环境下基于rss和tdoa测量的定位方法 | |
CN108445445B (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 |