CN110275133B - Non-common-view pulse signal passive positioning method based on virtual arrival time difference - Google Patents
Non-common-view pulse signal passive positioning method based on virtual arrival time difference Download PDFInfo
- Publication number
- CN110275133B CN110275133B CN201910567061.9A CN201910567061A CN110275133B CN 110275133 B CN110275133 B CN 110275133B CN 201910567061 A CN201910567061 A CN 201910567061A CN 110275133 B CN110275133 B CN 110275133B
- Authority
- CN
- China
- Prior art keywords
- toa
- pulse
- lost
- arrival time
- curve
- 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 13
- 230000005855 radiation Effects 0.000 claims abstract description 18
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 230000000630 rising effect Effects 0.000 claims abstract description 4
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000011426 transformation method Methods 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
- 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/06—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to a non-common-view pulse signal passive positioning method based on virtual arrival time difference, which comprises the following steps: receiving pulse signals continuously emitted by a radiation source u by using M observation platforms with the same time reference, detecting the rising edge of the received pulse signals at each observation platform, and estimating to obtain a time of arrival (TOA); simultaneously acquiring the position information of each platform corresponding to each arrival time; carrying out plane transformation on the received TOA sequence to obtain a two-dimensional image; the TOA sequence emitted by the same radiation source can form a continuous curve, and due to the non-common condition, part of points of the curve in the image are lost; fitting a curve in the two-dimensional image to obtain a curve parameter; estimating to obtain the coordinates of the lost points, and restoring the TOA of the lost pulse according to the coordinates of the lost points; estimating to obtain a virtual arrival time difference according to the TOA of the lost pulse; the virtual arrival time difference and the arrival time difference are jointly used for positioning, and the radiation source position is obtained through calculation by adopting a method for solving a nonlinear equation.
Description
Technical Field
The invention relates to the field of passive positioning, in particular to a non-common-view pulse signal passive positioning method based on virtual arrival time difference, aiming at a non-cooperative radiation source for transmitting pulse signals.
Background
The passive positioning means that the observation station is in a silent state, passively receives signals of the target radiation source, and positioning of the target radiation source is completed by receiving the signals. In practical application, due to the reasons of non-cooperation of targets, obstruction, different antenna direction directions and the like, signal loss can be caused, and an observation station cannot receive signals sent by all radiation sources, namely non-common-view. Whereas traditional observations, such as time difference of arrival, require that at least two platforms receive the same signal. Therefore, the non-common view causes a decrease in the conventional observation amount, resulting in a decrease in the positioning accuracy.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a passive positioning method for non-common-view pulse signals based on virtual arrival time differences, which can effectively improve the positioning accuracy.
In order to achieve the purpose, the invention adopts the following technical scheme: a non-common-view pulse signal passive positioning method based on virtual arrival time difference comprises the following steps: 1) receiving pulse signals continuously emitted by a radiation source u by using M observation platforms with the same time reference, detecting the rising edge of the received pulse signals at each observation platform, and estimating to obtain a time of arrival (TOA); simultaneously acquiring the position information of each platform corresponding to each arrival time; 2) carrying out plane transformation on the received TOA sequence, wherein the width of the plane transformation is w, and obtaining a two-dimensional image; the TOA sequence emitted by the same radiation source can form a continuous curve, and due to the non-common condition, part of points of the curve in the image are lost; 3) fitting a curve in the two-dimensional image to obtain a curve parameter; 4) estimate the derived lossRestoring the TOA of the lost pulse according to the coordinates of the lost point; 5) from the TOA of the missing pulse, the virtual arrival time difference is estimated6) By virtual arrival time differencesPositioning observation quantity can be obtained by combining the obtained time difference of arrival TDOA; the position of the observation platform corresponding to the TOA can be obtained by a self-positioning module of the observation platform, and the radiation source position can be obtained by resolving by adopting a method of solving a nonlinear equation in combination with the positioning observed quantity.
Further, in the step 4), after obtaining the coefficient of the curve fitting function, according to the cumulative plane transformation, wherein the X axis represents the observation window number and takes a value of a continuous positive integer, if a pulse is lost, the X coordinate of the midpoint of the curve is discontinuous, so that the corresponding point of the lost pulse in the cumulative plane transformation is estimated, and then the TOA of the lost pulse is obtained.
in the formula, f (x)m,i,αm) Representing a fitting function of the estimated curve. (ii) a x is the number ofm,iTransforming the abscissa of the corresponding point for the cumulative plane; alpha is alphamIs an estimation equation of the fitting coefficients.
wherein m is1≠m2,m1And m2Indicating the number of the two observation platforms receiving the ith pulse.
Due to the adoption of the technical scheme, the invention has the following advantages: aiming at the problem that a plurality of observation stations receive pulse signals of a radiation source, the positioning accuracy is reduced due to non-common vision, and a one-dimensional TOA sequence can be converted into a two-dimensional image by combining a plane transformation method according to the front-back constraint relation of the TOA sequence. The two-dimensional image has the characteristics that the ordinate corresponds to time and the value is continuous, and the abscissa corresponds to the number of the observation window and the value is a continuous positive integer. Therefore, the invention estimates the coordinates of the lost points in the curve by a curve fitting method, and effectively improves the positioning precision.
Drawings
FIG. 1 is a schematic overall flow diagram of the present invention;
FIG. 2 is a schematic view of the scenario of the present invention;
FIG. 3 is a schematic diagram of a transmit time series and an arrival time series of the present invention;
fig. 4 is a schematic diagram of a planar transformation of a transmit time series and an arrival time series of the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and examples.
As shown in fig. 1, the present invention provides a non-common-view pulse signal passive positioning method based on virtual arrival time difference, which includes the following steps:
1) as shown in fig. 2, M observation platforms with the same time reference are used to receive pulse signals continuously emitted by a radiation source u, and the rising edge of the received pulse signal is detected at each observation platform to obtain a time sequence of arrival TOA by estimation; and simultaneously acquiring the position information of each platform corresponding to each arrival time.
As shown in FIG. 3, let TOT be the transmission time of the ith pulsei:
TOTi=TOT0+i·PRI
Wherein, i is 1, 2, and N is the total number of transmitted pulses; TOT0PRI is the pulse repetition interval, which is the firing time of the initial pulse.
The arrival time t of the ith pulse received by the mth observation platformm,iComprises the following steps:
tm,i=TOTi+Δtm,i=TOT0+i·PRI+Δtm,i
wherein M is 1, 2, M; Δ tm,iThe propagation delay from the target to the observation station m for the ith pulse can be expressed as:
wherein c represents the speed of light, sm,iIndicating the position of the m-th observation platform when the i-th pulse is received. r ism,iDenotes the radiation sources u and sm,iThe distance between them. I | · | purple wind2Representing the two-norm of.
2) Carrying out plane transformation on the received TOA sequence, wherein the width of the plane transformation is w, and obtaining a two-dimensional image; the TOA sequence emitted by the same radiation source can form a continuous curve (as shown in fig. 4), and due to the non-common condition, the points of the curve part in the image are lost;
each TOA corresponds to a point (x, y) on the two-dimensional plane as:
x=mod(tm,i,w)
y=tm,i-wx;
3) fitting a curve in the two-dimensional image to obtain a curve parameter;
for the pulse received by the observation platform m, the arrival time is tm,iThe coordinates of the point corresponding to the cumulative planar transform are (x)m,i,ym,i) Assuming that the fitting function of the curve in the cumulative planar transform is f (x, α)m),αmRepresenting the undetermined coefficients in the fitting function, then:
error of fit em,iComprises the following steps:
em,i=ym,i-f(xm,i,αm)
error of fit em,iOf vector form emComprises the following steps:
em=[em,1,em,2,...,em,N]T
the estimated equation for the fitting coefficient is:
wherein |2Representing a 2 norm. This is a least squares fitting method whose principle is to make the residuals em,iHas the smallest sum of squares, i.e.
The basis functions for fitting can be expressed as
f(xm,i,αm)=f(xm,i,α1)+f(xm,i,α2)+…+f(xm,i,αK)
Where K is the number of basis functions.
In the present embodiment, a polynomial basis function is adopted, which can be expressed as
f(x,α0)=α0
f(x,α1)=α1x1
f(x,α2)=α2x2
f(x,αK)=αKxK-1
Wherein alpha ism=[α0,α1,...,αK]。
Let ym=Hmαm+emIn which H ism=[1,x,x2,...,xK-1]。
Solving polynomial basis functions, namely linear least square problem, can obtain a fitting coefficient alphamIs composed of
4) Estimating to obtain the coordinates of the lost point, and then restoring the TOA of the lost pulse according to the coordinates of the lost point: after the coefficients of the curve fitting function are obtained, the X-axis represents the number of the observation window and takes the value of a continuous positive integer according to the cumulative plane transformation, and if pulse loss occurs, the X coordinate of the midpoint of the curve is discontinuous. Therefore, the corresponding point of the missing pulse in the cumulative plane transformation can be estimated, and the TOA of the missing pulse can be obtained.
in the formula, f (x)m,i,αm) Representing a fitting function of the estimated curve.
Wherein m is1≠m2,m1And m2Indicating the number of the two observation platforms receiving the ith pulse.
6) By virtual arrival time differencesAnd combining the obtained arrival time difference TDOA to obtain a positioning observation quantity. The time difference of arrival can be expressed as
Wherein m is3≠m4,m3And m4Indicating the number of the two observation platforms receiving the ith pulse.
The self-positioning module of the observation platform can obtain the position of the observation platform corresponding to the TOA, and the observation quantity is positioned in combinationThe radiation source position can be obtained by solving a nonlinear equation. After the virtual arrival time difference observed quantity is increased, the positioning precision can be effectively improved, especially under the condition that more pulses are lost.
The above embodiments are only for illustrating the present invention, and the steps may be changed, and on the basis of the technical solution of the present invention, the modification and equivalent changes of the individual steps according to the principle of the present invention should not be excluded from the protection scope of the present invention.
Claims (1)
1. A non-common-view pulse signal passive positioning method based on virtual arrival time difference is characterized by comprising the following steps:
1) receiving pulse signals continuously emitted by a radiation source u by using M observation platforms with the same time reference, detecting the rising edge of the received pulse signals at each observation platform, and estimating to obtain a time of arrival (TOA); simultaneously acquiring the position information of each platform corresponding to each arrival time;
2) carrying out plane transformation on the received TOA sequence, wherein the width of the plane transformation is w, and obtaining a two-dimensional image; the TOA sequence emitted by the same radiation source can form a continuous curve, and due to the non-common condition, part of points of the curve in the image are lost;
3) fitting a curve in the two-dimensional image to obtain a curve parameter;
4) estimating to obtain the coordinates of the lost points, and then restoring the TOA of the lost pulse according to the coordinates of the lost points;
6) By virtual arrival time differencesPositioning observation quantity can be obtained by combining the obtained time difference of arrival TDOA; the position of the observation platform corresponding to the TOA can be obtained by a self-positioning module of the observation platform, and the radiation source position can be obtained by resolving by adopting a method of solving a nonlinear equation in combination with the positioning observed quantity;
in the step 4), after the coefficient of the curve fitting function is obtained, according to the cumulative plane transformation, wherein the X axis represents the number of the observation window and takes the value as a continuous positive integer, if the pulse is lost, the X coordinate of the middle point of the curve is discontinuous, so that the corresponding point of the lost pulse in the cumulative plane transformation is estimated, and the TOA of the lost pulse is further obtained;
in the formula (I), the compound is shown in the specification,a fitting function representing the estimated curve; x is the number ofm,iTransforming the abscissa of the corresponding point for the cumulative plane;an estimation equation for the fitting coefficients;
wherein m is1≠m2,m1And m2Indicating the number of the two observation platforms receiving the ith pulse.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910567061.9A CN110275133B (en) | 2019-06-27 | 2019-06-27 | Non-common-view pulse signal passive positioning method based on virtual arrival time difference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910567061.9A CN110275133B (en) | 2019-06-27 | 2019-06-27 | Non-common-view pulse signal passive positioning method based on virtual arrival time difference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110275133A CN110275133A (en) | 2019-09-24 |
CN110275133B true CN110275133B (en) | 2021-06-29 |
Family
ID=67962400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910567061.9A Active CN110275133B (en) | 2019-06-27 | 2019-06-27 | Non-common-view pulse signal passive positioning method based on virtual arrival time difference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110275133B (en) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2359903A1 (en) * | 2008-11-21 | 2011-08-24 | Tohoku University | Signal processing device, signal processing method, signal processing program, computer-readable recording medium on which is recorded a signal processing program, and radiation therapy device |
CN102608573A (en) * | 2012-03-29 | 2012-07-25 | 清华大学 | Mutual-fuzzy-accumulation passive location method based on multiple observing points |
CN103645485A (en) * | 2013-10-28 | 2014-03-19 | 中国科学院国家授时中心 | Pseudorange differential method based on dual-satellite time difference and frequency difference passive positioning |
CN103713277A (en) * | 2013-12-19 | 2014-04-09 | 中国航天科工集团八五一一研究所 | Location information field-based radiation source localization algorithm |
WO2014086688A1 (en) * | 2012-12-07 | 2014-06-12 | Thales | Method for the passive localization of radar transmitters |
CN105353345A (en) * | 2015-10-27 | 2016-02-24 | 清华大学 | Single-channel radar signal passive positioning method based on higher-order observed quantities |
CN105572635A (en) * | 2016-01-25 | 2016-05-11 | 西安电子科技大学 | Single-station passive quick positioning method based on least square method |
CN106324563A (en) * | 2016-08-12 | 2017-01-11 | 中国电子进出口总公司 | Multipoint passive detection all-phase signal sorting and arrival time difference measurement system |
CN107526089A (en) * | 2017-08-25 | 2017-12-29 | 清华大学 | A kind of non-based on time delay second order difference regards radar signal passive location method altogether |
CN107870317A (en) * | 2017-11-10 | 2018-04-03 | 中国电子科技集团公司第四十研究所 | A kind of narrow band signal TDOA localization methods and device calculated based on the time difference |
CN108112071A (en) * | 2016-11-11 | 2018-06-01 | 中兴通讯股份有限公司 | Localization method, locating base station, location-server and alignment system |
CN109541531A (en) * | 2018-11-22 | 2019-03-29 | 中电科仪器仪表有限公司 | A method of reducing frequency difference influences time delay estimadon |
US10302758B2 (en) * | 2016-11-10 | 2019-05-28 | China University Of Mining & Technology, Beijing | Method and device for detecting discontinuous body with ground penetrating radar |
CN109917333A (en) * | 2019-01-31 | 2019-06-21 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Merge the passive location method of AOA observed quantity and TDOA observed quantity |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BR9917165A (en) * | 1999-01-08 | 2002-01-22 | Trueposition Inc | Method for use in a wireless location system and method for locating a transmitter |
CN103234564B (en) * | 2013-04-15 | 2015-11-11 | 山东联友通信科技发展有限公司 | Based on localization method and the positioning system of two-dimensional magnetic coding |
-
2019
- 2019-06-27 CN CN201910567061.9A patent/CN110275133B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2359903A1 (en) * | 2008-11-21 | 2011-08-24 | Tohoku University | Signal processing device, signal processing method, signal processing program, computer-readable recording medium on which is recorded a signal processing program, and radiation therapy device |
CN102608573A (en) * | 2012-03-29 | 2012-07-25 | 清华大学 | Mutual-fuzzy-accumulation passive location method based on multiple observing points |
WO2014086688A1 (en) * | 2012-12-07 | 2014-06-12 | Thales | Method for the passive localization of radar transmitters |
CN103645485A (en) * | 2013-10-28 | 2014-03-19 | 中国科学院国家授时中心 | Pseudorange differential method based on dual-satellite time difference and frequency difference passive positioning |
CN103713277A (en) * | 2013-12-19 | 2014-04-09 | 中国航天科工集团八五一一研究所 | Location information field-based radiation source localization algorithm |
CN105353345A (en) * | 2015-10-27 | 2016-02-24 | 清华大学 | Single-channel radar signal passive positioning method based on higher-order observed quantities |
CN105572635A (en) * | 2016-01-25 | 2016-05-11 | 西安电子科技大学 | Single-station passive quick positioning method based on least square method |
CN106324563A (en) * | 2016-08-12 | 2017-01-11 | 中国电子进出口总公司 | Multipoint passive detection all-phase signal sorting and arrival time difference measurement system |
US10302758B2 (en) * | 2016-11-10 | 2019-05-28 | China University Of Mining & Technology, Beijing | Method and device for detecting discontinuous body with ground penetrating radar |
CN108112071A (en) * | 2016-11-11 | 2018-06-01 | 中兴通讯股份有限公司 | Localization method, locating base station, location-server and alignment system |
CN107526089A (en) * | 2017-08-25 | 2017-12-29 | 清华大学 | A kind of non-based on time delay second order difference regards radar signal passive location method altogether |
CN107870317A (en) * | 2017-11-10 | 2018-04-03 | 中国电子科技集团公司第四十研究所 | A kind of narrow band signal TDOA localization methods and device calculated based on the time difference |
CN109541531A (en) * | 2018-11-22 | 2019-03-29 | 中电科仪器仪表有限公司 | A method of reducing frequency difference influences time delay estimadon |
CN109917333A (en) * | 2019-01-31 | 2019-06-21 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Merge the passive location method of AOA observed quantity and TDOA observed quantity |
Non-Patent Citations (5)
Title |
---|
Joint TDOA-DOA Localization Scheme for Passive Coherent Location Systems;Meng-Chang Hua;《8th IEEE, IET International Symposium on Communication Systems, Networks and Digital Signal Processing》;20121231;1-4 * |
Multi-Sensor Passive Localization Using Second Difference of Coherent Time Delays With Incomplete Measurements;SHANGYU ZHANG;《IEEE Access》;20190413;43167-43178 * |
TDOA定位中到达时间及时间差误差的统计模型;张志良;《重庆大学学报》;20060131;第29卷(第1期);85-88 * |
基于卡尔曼滤波的TDOA/AOA混合定位算法;张怡;《计算机工程与应用》;20151231;62-66 * |
运动观测站测时无源定位新方法研究;刘洋;《中国博士学位论文全文数据库 信息科技辑》;20171115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110275133A (en) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9407317B2 (en) | Differential ultra-wideband indoor positioning method | |
CN108717184B (en) | Error correction-based DOA and TOA combined single-station passive positioning method | |
CN109633581B (en) | Positioning method based on external radiation source TDOA/FDOA error correction | |
US9891306B2 (en) | Geolocating a remote emitter | |
KR101040705B1 (en) | The system and method for wireless positioning | |
KR101767924B1 (en) | Method and system for detecting location of multi-target | |
CN109870695B (en) | Deep sea bottom reflected sound-based non-cooperative target multi-node underwater sound positioning method | |
CN107526089B (en) | A kind of non-view radar signal passive location method altogether based on time delay second order difference | |
KR101789188B1 (en) | An underwater integrated navigation system for tracking underwater moving objects | |
CN106249197B (en) | The method for self-calibrating of receiver location error in a kind of multipoint location system | |
KR100977246B1 (en) | Apparatus and method for estmating positon using forward link angle of arrival | |
JP7326422B2 (en) | Positioning method and system for locating at least one object using wave-based signals | |
EP2866046B1 (en) | Self-locating mobile receiving device | |
CN112887901A (en) | Convex optimization target positioning method based on quantitative TOA measurement | |
JP2012083136A (en) | Bias error estimation apparatus, bias error estimation method and position estimation apparatus | |
Lui et al. | Range-based source localisation with pure reflector in presence of multipath propagation | |
KR101627419B1 (en) | Method for estmating location of mobile node and apparatus thereof | |
CN110471029B (en) | Single-station passive positioning method and device based on extended Kalman filtering | |
KR100882590B1 (en) | Device and method for measuring location | |
CN109856623B (en) | Target state estimation method for multi-radar linear flight path line | |
CN109752690B (en) | Method, system and device for eliminating NLOS (non-line of sight) positioned by unmanned aerial vehicle and storage medium | |
CN110275133B (en) | Non-common-view pulse signal passive positioning method based on virtual arrival time difference | |
RU2562616C1 (en) | Method of acquiring radio information and radio system therefor | |
CN112684411B (en) | Underwater target positioning method based on improved arrival frequency difference | |
US10156630B2 (en) | Method for passively locating a non-movable transmitter |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Huang Zhen Inventor after: Shi Lei Inventor after: Zhang Shangyu Inventor after: He Jiazhi Inventor after: Feng Xuefeng Inventor before: Huang Zhen Inventor before: Zhang Shangyu Inventor before: Feng Xuefeng Inventor before: He Jiazhi Inventor before: Shi Lei |