CN112904274A - Multi-moving-object positioning method for improving TDOA/FDOA algorithm - Google Patents
Multi-moving-object positioning method for improving TDOA/FDOA algorithm Download PDFInfo
- Publication number
- CN112904274A CN112904274A CN202110082721.1A CN202110082721A CN112904274A CN 112904274 A CN112904274 A CN 112904274A CN 202110082721 A CN202110082721 A CN 202110082721A CN 112904274 A CN112904274 A CN 112904274A
- Authority
- CN
- China
- Prior art keywords
- target
- sensor
- positioning
- equation
- algorithm
- 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.)
- Pending
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/0205—Details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Radar, Positioning & Navigation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to a multi-moving-object positioning method for improving a TDOA/FDOA algorithm, which mainly solves the problems of poor adaptability of a traditional algorithm to noise errors in positioning of a plurality of dynamic objects and limited positioning accuracy. Firstly, setting K non-coincident moving targets in a three-dimensional space, deploying M sensors, and acquiring a measurement vector of target information; then, a TDOA/FDOA positioning equation system related to the position information of the target and the sensor is constructed, so that the speed and the position of the target are estimated; and constructing a global constraint condition according to the relation between the target and the additional variable, introducing a Lagrange multiplier technology, and solving a positioning equation by adopting a quasi-Newton BFGS iterative formula. Constraint conditions are considered comprehensively and systematically, and a Lagrange multiplier technology and a quasi-Newton BFGS iterative formula are introduced, so that the Hessian matrix is prevented from being calculated, the calculated amount is greatly reduced, and the positioning accuracy is improved.
Description
Technical Field
The invention relates to the technical field of wireless sensor networks, in particular to a multi-moving-object positioning method for improving a TDOA/FDOA algorithm.
Background
In real life, no matter in aspects of national defense and military, medical assistance, intelligent transportation, environmental monitoring, industry or agricultural control and the like, the WSN (wireless sensor network) is widely applied and plays an increasingly important role, while location-based services are a key core part of all applications, and a positioning target generally occurs in a noise environment, so that the measurement error is large, and the error is also large when the target positioning is realized by adopting a single TDOA (time difference of arrival) ranging technology. In order to improve the positioning accuracy of the algorithm, a hybrid positioning technology is usually adopted to estimate the position information of the target, and when the position information estimation is carried out for multiple moving targets at the same time, a positioning model combining TDOA and FDOA is usually required to estimate the speed and the position of the target, so that the positioning accuracy of the moving target is effectively improved. However, when a dynamic target is positioned, a large positioning deviation occurs in the sensor position error, and therefore the sensor position error becomes an important factor to be considered when estimating the target position.
In the prior art, when a multi-motion target is located by a TSWLS algorithm based on a TDOA/FDOA combined measurement technology and combined with a position error of a sensor, an experimental result shows that: the location estimate RMSE can only reach CRLB if the signal-to-noise ratio is high. Because the measurement error exists together with the sensor position error, the TSWLS algorithm has poor adaptability in terms of noise error, so that the positioning deviation is large and the accuracy is insufficient.
Further, the constrained total least square algorithm system proposed in the prior art considers all coefficient matrix noise in the positioning equation and solves the equation according to the newton iteration method, so as to obtain the position information of the target, and the experimental result shows that: compared with the TSWLS algorithm, the algorithm has better positioning performance and can reach CRLB under the condition of larger noise. The disadvantage is that the algorithm is based on TDOA ranging techniques to estimate the position information of a fixed object and is therefore not suitable for locating moving objects. Furthermore, by analyzing the main factors influencing the positioning performance of the TSWLS algorithm, the invention provides a moving object positioning algorithm based on TDOA/FDOA to correct the positioning error, and the positioning error can be reduced and the capability of adapting to the measurement noise can be improved without increasing the operation complexity.
Aiming at the problem of positioning of related information among multiple (group) targets, Yang Zeu, university of strategic support army information engineering, provides a positioning algorithm for coordinating target position and speed constraint. According to the relevance information of the position and the speed parameters between the targets, on the basis of an observation model combining Time Difference of Arrival (TDOA) and Frequency Difference of Arrival (FDOA), the relevance information of the distance and the speed between the targets is converted into inequality constraints for improving the positioning precision, and the estimation performance bound of the position and the speed of the targets containing the constraints is deduced. Compared with an unconstrained positioning algorithm, the provided algorithm has higher positioning precision and robustness to noise, and the performance of the algorithm can approach the positioning performance bound with constraint. The main problem is that a fixed multi-target positioning algorithm for collaborative constraint information is researched, but in actual positioning scenes and applications, the target position usually changes dynamically, so that the constraint condition may change dynamically within a specific range, and the maximum likelihood estimation algorithm has no advantage in dynamic target positioning.
Therefore, the position error and the measurement error of the sensor are analyzed in detail on the basis of the prior art, an Improved constrained total least square (ICTRS) algorithm is provided for estimating the position coordinates of a plurality of non-coincident moving targets, and a Lagrange multiplier is introduced into the algorithm and a BFGS formula of a quasi-Newton method is utilized to solve the CTLS target positioning problem. The experimental result shows that the algorithm has better positioning performance than the TSWLS algorithm, can greatly reduce the positioning deviation, has stronger adaptability to the measurement error and the sensor position error, and can reach CRLB even under the condition of larger noise error.
However, no matter the classic TSWLS, CTLS or ICTLS positioning algorithms, there is a fatal defect in positioning the target, which results in a large positioning deviation and a sharp drop in positioning performance. When the target is positioned on or close to any one coordinate axis of the set reference sensor, the positioning algorithm has larger positioning deviation, and the positioning performance is sharply reduced.
Disclosure of Invention
The invention aims to solve the technical problem of providing a multi-moving object positioning method for improving a TDOA/FDOA algorithm.
In order to solve the technical problems, the invention adopts the following technical scheme:
a multi-moving object positioning method for improving TDOA/FDOA algorithm comprises the following steps:
step 1, setting K non-coincident moving targets in a three-dimensional space, namely a target 1 and a target 2 … …, respectively, deploying M sensors in the three-dimensional space, namely a sensor 1 and a sensor 2 … …, taking the sensor 1 as a reference sensor, and constructing an FDOA measurement vector of speed information of a target i:
wherein, i is 1,2 … … K, j is 1,2 … … M, rj1,iFor the true value of the difference in distance from the object i to the sensor 1 and the sensor j,the true distance between target i and sensor j,is the true distance, Δ r, between the target i and the sensor 1j1,iThe difference value between the actual value and the measured value of the distance difference between the target i and the sensor 1 and the sensor j is an error value;
step 2: constructing a TDOA/FDOA system of positioning equations for object i and sensor location information, introducing additional variablesWherein u isiIs the position of the target i and,to transmitTrue position of sensor jAndsubstituted typeAnd squaring the two sides, and obtaining a TDOA equation set after arrangement:
and step 3: will be provided withAt the sensor noise position s1And performing first-order Taylor series expansion, neglecting higher-order terms with more than two orders, entering a TDOA equation set to obtain a new equation set, differentiating the new equation set with time to obtain an FDOA equation set, and synthesizing the new equation set and the FDOA equation set to obtain a positioning equation set in a matrix form: (a + Δ a) θ ═ b + Δ b; in the formula (I), the compound is shown in the specification,w=ui-s1,data matrixData vector Is r1,iIs determined by the estimated value of (c),is composed ofAn estimated value of (d);
and 4, step 4: taking into account the errors of all coefficient matrices in the equation, and incorporating additional variablesAndand w andconstraint condition between, get constraint condition thetaTΣ1θ=0,θTΣ2θ=0;
And 5: constructing a value function of a vector theta to be solved, and changing the total least square solution of the vector theta for solving the positioning equation set into a CTLS solution;
step 6: introducing Lagrange multiplier technique lambda1And λ2Solving the problem of solving the CTLS constraint optimization minimum value to obtain the final estimated value of the vector theta
θ=[(A-D)TWθA+λ1Σ1+λ2Σ2]-1(A-D)TWθb (18)
Substituting equation (18) into θ againTΣ1θ=0,θTΣ2θ is 0 to yield:
θTΣ1θ=0,θTΣ2θ=0 (19);
where θ is an estimate of θ;
and 7: definition of λ ═ λ1 λ2]T,f=[f1 f2]TSolving λ in equation (19) using a quasi-newton iterative equation:
λ(k+1)=λ(k)-α(k)H(k)g(k) (21)
in the formula, λ(k)K number of λIteration, α(k)In order to be the step size,H(k)and updating the symmetric matrix of the corresponding dimension by a BFGS formula as follows:
let parameter delta be (0,1), sigma be (0,0.5), Wθ=Q-1Precision of 10-3Taking an initial point λ(1)=[0,0]TInitial matrix H(1)Is a unit matrix I, k is 0;
and 8: computingIf g | | |(k)Stopping iteration if | | is less than or equal to epsilon, and outputting lambdakAs an approximate solution;
and step 9: calculating a search direction p(k)=-H(k)g(k)And is provided with mkIs satisfied with f (λ)(k))+δmp(k)≤f(λ(k))+σδmg(k )Tp(k)The smallest non-negative integer of (c);
step 10: updating H by BFGS formula of quasi-Newton method(k+1)Setting k to k +1, and turning to step 8;
step 11: lambda to be obtained finallykOptimal value substitutionThus, the optimal estimation value of the vector theta to be solved is obtained, wherein theta is the final estimation position information of the target i.
Further, the step 1 specifically includes the following steps:
The real coordinates of the sensor are represented,which is indicative of the true speed of the sensor,andrespectively representing the true position and velocity of the jth sensor [. ]]0Is []True value of [ ·]TIs []The transposition operation of (1);
the sensor with j equal to 1 is used as the reference sensor, so the real distance between the target i and the sensor j is
The actual and measured values of the difference in distance from the object i to the sensors 1 and j are respectively
The TDOA-form measurement vector of the target i-position information is recorded as
ri=[r21,i,r31,i,…,rM1,i]T=ri 0+Δri (3)
FDOA measurement of target velocity information:
further, the step 3 comprises the following steps:
will be provided withAt the sensor noise position s1And performing first-order Taylor series expansion on the position, neglecting higher-order terms above the second order, and obtaining after arrangement:
differentiating equation (7) with respect to time to obtain the FDOA equation set:
integrating (7) and (8) to obtain a positioning equation system in a matrix form:
(A+ΔA)θ=b+Δb (9)
Further, the step 4 comprises the following steps:
taking into account all of the equationsError in coefficient matrix, combined with additional variablesAndand w andand (4) estimating the target position information by adopting a CTLS positioning algorithm under the constraint condition. The cost function of the candidate vector θ can be configured as:
J(θ)=(Aθ-b)TWθ(Aθ-b) (10)
wherein J (θ) is a cost function of θ, WθIs a weighting matrix. Theta is an estimated value of theta, and theta is an optimized value of theta. Due to the introduced variablesAndand w andsince there is a certain constraint relationship, equation (11) needs to be solved by combining the constraint condition. The constraint conditions are obtained by integrating (11) and (4)
Transforming the formula (12) to obtain the formula (10), and arranging to obtain the constraint condition in a matrix form:
θTΣ1θ=0,θTΣ2θ=0 (13)
in the formula (I), the compound is shown in the specification,o is a 4 × 4 all-0 matrix, sigma11=Σ12=Σ22=diag{1,1,1,-1}。
Further, the step 5 comprises the following steps:
the overall least squares solution for vector θ becomes the CTLS solution:
S.t.θTΣ1θ=0,θTΣ2θ=0 (15)
2[(A-D)TWθA+λ1Σ1+λ2Σ2]θ-2(A-D)TWθb=0 (16)
further, the step 6 comprises the following steps:
equation (16) is a quadratic constraint non-convex optimization problem, and solving the problem comprises the steps of firstly solving partial differential of equation (16) on theta and then enabling equation (16) to be subjected to partial differentialAfter finishing, obtaining:
2[(A-D)TWθA+λ1Σ1+λ2Σ2]θ-2(A-D)TWθb=0 (17)
in the formula (I), the compound is shown in the specification,l-8 is the dimension of the vector theta,Q=E(nnT) Is an initial weighting matrixAnd n ═ Δ βT,ΔαT]TAs random noise, F1,…,FLA coefficient matrix which is a disturbance matrix delta A, delta b; thus, the final estimate of vector θ is
θ=[(A-D)TWθA+λ1Σ1+λ2Σ2]-1(A-D)TWθb (18)
Substituting equation (18) into constraint equation (13) again yields:
θTΣ1θ=0,θTΣ2θ=0 (19)
formula (19) relates to1And λ2By solving the system of equations to obtain λ1And λ2And substituting the position information into the formula (18) to obtain the final estimated position information of the target i.
After the technical scheme is adopted, compared with the prior art, the invention has the following advantages:
the invention provides a correction constraint total least square algorithm-ICTLS (information communication technology) aiming at the problems of low precision and large calculation amount of the traditional algorithm when a plurality of dynamic targets are positioned. Constraint conditions are considered comprehensively and systematically by the algorithm, and a Lagrange multiplier technology and a quasi-Newton BFGS iterative formula are introduced, so that the calculation of a Hessian matrix is avoided, the calculation amount is greatly reduced, and the positioning accuracy is improved.
The present invention will be described in detail below with reference to the accompanying drawings and examples.
Drawings
FIG. 1 is a diagram of the spatial location of a target and sensors according to the present invention;
FIG. 2 is a comparison of RMSE and CRLB for the next three positioning algorithms under the present scenario;
FIG. 3 is a comparison of RMSE and CRLB for the next three positioning algorithms according to the present invention scenario;
FIG. 4 is a scenario of the present invention (RMSE versus CRLB for the next three positioning algorithms);
FIG. 5 is a target and reference sensor profile of the present invention;
FIG. 6 is a block diagram of the optimization algorithm flow of the present invention;
FIG. 7 is a distribution of target and sensor positions in a particular scenario of the present invention;
FIG. 8 is a comparison of the performance of 2 algorithms in locating close range targets in accordance with the present invention;
FIG. 9 is a comparison of the performance of 2 algorithms of the present invention in locating a remote target;
FIG. 10 is a distribution of target and sensor positions in a stochastic scenario of the invention;
FIG. 11 is a comparison of the performance of the 2 algorithms when the targets and sensors of the present invention are randomly distributed;
FIG. 12 is a diagram of the position estimation MSE and bias versus sensor count in accordance with the present invention;
Detailed Description
The principles and features of this invention are described below in conjunction with the following drawings, which are set forth by way of illustration only and are not intended to limit the scope of the invention.
The working principle of the invention is as follows:
(1) setting K moving targets which are not overlapped mutually in a three-dimensional space, deploying M sensors, and constructing a measurement vector of target information:
(1.1) let the i-th target's true position and velocity be uiAndalso deploying M sensors in three-dimensional space, orderThe real coordinates of the sensor are represented,which is indicative of the true speed of the sensor,and respectively representing the true position and velocity of the jth sensor [. ]]0Is []True value of [ ·]TIs []The transpose operation of (1).
(1.2) using the sensor with j equal to 1 as the reference sensor, the real distance between the target i and the sensor j is
The actual and measured values of the difference in distance from the object i to the sensors 1 and j are respectively
The TDOA-form measurement vector of the target i-position information is recorded as
ri=[r21,i,r31,i,…,rM1,i]T=ri 0+Δri (3)
FDOA measurement of target velocity information:
(2) constructing a TDOA/FDOA location equation set for the object and sensor location information to estimate the object's velocity and position:
(2.1) introduction of additional variablesHandleAndsubstituting the formula (4), squaring the two sides, and finishing to obtainTo the TDOA system of equations:
(2.2) mixingAt the sensor noise position s1And performing first-order Taylor series expansion on the position, neglecting higher-order terms above the second order, and obtaining after arrangement:
differentiating equation (7) with respect to time to obtain the FDOA equation set:
(2.3) integrating the equations (7) and (8) to obtain a positioning equation system in a matrix form:
(A+ΔA)θ=b+Δb (9)
(3) constructing a global constraint condition according to the relation between the target and the additional variable, introducing a Lagrange multiplier, and solving a positioning equation by adopting a quasi-Newton BFGS iterative formula:
(3.1) taking into account the errors present in all the coefficient matrices in the equation, and incorporating additional variablesAndand w andand (4) estimating the target position information by adopting a CTLS positioning algorithm under the constraint condition. The cost function of the candidate vector θ can be configured as:
J(θ)=(Aθ-b)TWθ(Aθ-b) (10)
(3.2) solving the constraint condition of the comprehensive formula (11) and (4)
(3.3) the overall least squares solution of vector θ becomes the CTLS solution:
S.t.θTΣ1θ=0,θTΣ2θ=0 (14)
2[(A-D)TWθA+λ1Σ1+λ2Σ2]θ-2(A-D)TWθb=0 (15)
(3.5) definition of λ ═ λ1 λ2]T,f=[f1 f2]TAnd solving lambda by adopting a quasi-Newton iterative formula.
(4) In a specific scene, for example, when the target i is exactly located on a coordinate axis, a first step of primary estimation of the TSWLS algorithm is adopted to obtain a target position, and then the reference sensor and the rotating coordinate axis are reselected according to the estimated target position to correct the defects of the TSWLS algorithm:
(4.1) referring to the positioning scene and model establishment in the step (1), and expressing a measurement equation based on a TDOA algorithm as follows:
r=r0+Δr (16)
wherein r is ═ r2,k,r3,k,…,rM,k]T,Δr=[Δr2,k,Δr3,k,…,ΔrM,k]TThe error vector is measured randomly as a Gaussian whose mean is zero, and the covariance matrix is Qt。
Introducing additional variablesPlacing it inPerforms expansion of a first order Taylor series, and only retains Δ skLinear terms, which after sorting can yield:
(4.2) orderIn combination with the full pseudolinearity equation of equation (17), one can derive the equation for the vectorThe localization equation of (c):
(4.3) when vectorCan be regarded as the estimation error ofIs thatA random vector fluctuating up and down nearby, so that the estimated vectorCan be expressed as:
(4.4) define reference Sensors for further analysis why reselectWhen the coordinate axes are rotated for the first time, phi is 45 degrees, the coordinate axes are rotated clockwise with the z-axis as the rotation axis, and the rotation angle alpha is pi/4. Then the ith sensor located under the new coordinate system can be expressed as:
in the second rotation of the coordinate axes, let θ ≈ 55 ° to pass through the origin of coordinates and be perpendicular to the planeIs clockwise rotated with the rotation angle β being 55 pi/180, and the cosine direction of the rotation axis is n z0. The ith sensor can be expressed as:
Thus, the estimated position of the target under the new coordinate system is
Wherein pi is diag { sgn (u '(1: 3) -s')k) And sgn (×) is a sign function.
Resulting in a sudden increase in estimation error when the target is at or near the reference sensor. On the basis of keeping a closed-form solution, the optimization algorithm effectively reduces the position estimation error by reselecting the reference sensor and the rotating coordinate system, corrects the defects of the TSWLS algorithm and comprehensively improves the positioning performance of the TSWLS algorithm.
The effects of the invention can be further illustrated by simulation:
A. simulation conditions
In order to verify the feasibility of the ICTLS algorithm for positioning a plurality of dynamic targets, the invention designs a simulation experiment of three target positioning scenes. In order to compare and analyze the positioning performance of the ICTLS algorithm in detail, the simulation experiment in the text compares and analyzes the simulation results of the three positioning algorithms of TSWLS, CTLS, and ICTLS. Firstly, a single dynamic target is positioned, and the real position and speed of the target are respectivelyAnd② locating two different dynamic targets, except for targetIn addition to the true position and velocity of another objectAre respectively asAndpositioning two different static targets, assuming two targetsAndspeed is zero and selected close toSensor (2)Is a reference sensor. The distribution of the positions of the targets and sensors in space is shown in fig. 1. The experiment assumes that there are 6 observation sensors in three-dimensional space, and the specific positions and velocities of the sensors are shown in table 1. The simulation experiment of multi-target positioning is carried out by taking two targets as an example, and the simulation experiment can be extended to more targets. The real positions of two targets to be positioned are respectivelyAndthe two target speeds are respectivelyAndlet the position error of the sensor be a Gaussian random process with zero mean and the covariance matrix be Andrespectively, the position and velocity error variances of the sensor, anThe matrix R ═ diag [1,1,1,2,2,2,10,10,10,40,40,40,20,20,20,3,3,3]. The measurement errors based on TDOA and FDOA are Gaussian random vectors with zero mean value, and the covariance matrix is Andfor measuring time difference and frequency difference variance, respectivelyAnd J is a (M-1) × (M +1) square matrix with a major diagonal element of 1 and other elements of 0.5.
TABLE 1 location distribution of targets and sensors in space
Aiming at the positioning algorithm in a specific scene, the performance of the positioning target of the optimization algorithm and the TSWLS algorithm is compared through simulation experiments in two positioning scenes. Two positioning scenarios are: a special positioning scene, namely that the target is positioned at or close to a reference sensor; and secondly, randomly positioning scenes, namely randomly distributing the target and the sensor in a three-dimensional space.
Firstly, in a special positioning scene, 2 targets to be positioned are distributed in a monitoring area, and one target is a short-distance targetThe other is a distant targetThe two targets to be positioned estimate the target position information from the data measured by 6 sensors in the area, the real position coordinates of the sensors are shown in table 2, and the distribution of the target sensors in space is shown in fig. 7.
TABLE 2 true position of sensor
Positioning targets in a random scene, wherein 100 target sources to be positioned are randomly distributed in a set space, and the target sources are randomly distributed in a cube with the side length of 600 m. Likewise, 100 sensors are randomly distributed in a cube with a side length of 200m, and the centers of the two cubes coincide. And randomly selecting 6 sensors and 1 target as a positioning scene. The distribution of targets and sensors in space is shown in fig. 10.
B. Emulated content
In order to verify the feasibility of the ICTLS algorithm in locating a plurality of dynamic targets, 10000 Monte Carlo experiments are run on a PC, and the RMSE of the ICTLS algorithm and the TSWLS and CTLS methods for target position and speed estimation is compared. The RMSE of the target estimated position and velocity is defined as follows:
wherein: P-RMSE and V-RMSE represent the RMSE, u, of the estimated position and velocity of the target, respectivelyi(l) The estimation result of the first Monte Carlo experiment is shown.
5000 Monte Carlo experiments are run on a PC (personal computer) aiming at experiments under two positioning scenes in a specific scene, and the positioning performance of the algorithm is analyzed by comparing the MSE and the deviation of the position estimation of the two algorithms. The MSE and bias of the position estimate are defined as
In the formula ulBit ulAnd L is 5000. The position error variance of the sensor isAnd covariance matrixThe variance of the mean measurement error isHaving a covariance matrix ofAnd R is an (M-1) × (M-1) matrix with a main diagonal element of 1 and other elements of 0.5; it is assumed that the measurement error is independent of the position error of the sensor.
C. Simulation result
FIG. 2 shows the positioning performance of the three positioning algorithms TSWLS, CTLS and ICTLS when positioning a single dynamic target, and the comparison of the RMSE and CRLB of the target position velocity estimation bias. The curves in the figure show that: for the estimation of the target position, the method includesIn time, the RMSE of the position estimation biases of the three positioning algorithms can all reach CRLB. However, the three positioning algorithms of TSWLS, CTLS and ICTLS are respectively inThe RMSE of the temporal position estimation error is gradually higher than CRLB. The reason why the three positioning algorithms of TSWLS, CTLS and ICTLS have different thresholds is that the WLS result deviation of the first step of the TSWLS positioning algorithm is increased along with the increase of noise errors, so that the final positioning accuracy is poor and the robustness is weak;although the CTLS positioning algorithm integrally considers factors such as measurement errors and position errors of the sensors, the additional variables introduced are not consideredWith target position coordinatesThe relation between the CTLS and the CTLS cannot reach the optimal solution because the correlation between the measurement noise and the position error of the sensor is not considered; in the ICTLS positioning algorithm, the correlation between measurement noise and the position error of the sensor is considered on the basis of the CTLS algorithm, the introduced additional variable is subjected to approximate processing, the estimation error of the first step is reduced, meanwhile, the correlation between the introduced variable and the target position is considered in the positioning algorithm, and the optimal solution can be achieved under the constraint condition.
Therefore, the ICTLS positioning algorithm in the text is optimal in positioning dynamic target performance among the three positioning algorithms. (ii) target speed estimation aspect whenIn the process, the RMSE of the speed estimation of the three positioning algorithms of TSWLS, CTLS and ICTLS can reach CRLB, and the threshold values of the three positioning algorithms are respectively 2.5dBm2,5dBm2,7.5dBm2Therefore, the position estimation threshold values corresponding to the three positioning algorithms are respectively larger than the speed estimation threshold values. Therefore, the three positioning algorithms have a poorer noise adaptation capability for speed estimation when positioning a dynamic target than for position estimation, and the randomness of estimation errors is increased.
FIG. 3 depicts the RMSE and CRLB alignment of the three TSWLS, CTLS and ICTLS localization algorithms to the target position and velocity estimation errors when localizing a plurality of dynamic targets. The RMSE curves for the position and velocity estimate bias in FIG. 3, the upper curve being the positioning targetThe lower curve is a positioning targetRMSE curve of (2). The two dynamic targets are positioned according to the curves in the graphAndthe estimation of the ICTLS positioning algorithm no matter the position or the speed has higher precision and stronger robustness than other two algorithms, and even if the measurement error and the sensor position error are larger, the CRLB can also be achieved. In addition, TSWLS and CTLS positioning algorithms obtained from curves in the graph are used for positioning close-range targetsThe threshold value occurring in the estimation of the target position or velocity is larger than that of a distant targetThe threshold of (a) is small, which indicates that the performance of the TSWLS and CTLS positioning algorithms for positioning close-range targets is poor, mainly due to the fact that the noise interaction between close-range targets is large and the correlation of noise errors is not considered.
FIG. 4 depicts the results of RMSE versus CRLB comparisons of position estimation errors for three positioning algorithms TSWLS, CTLS, and ICTLS in locating two static targets. For long distance targetsWhen in useThe RMSE of all three positioning algorithms can reach CRLB even inThe RMSE of the three positioning algorithms increases along with the increase of the error noise power, but under the condition of the same noise power increment, the static state positioning algorithms have the same noise power incrementThe RMSE increment of the target position estimation error is smaller than that of the dynamic target position estimation error because the static target is positioned without parameters such as speed frequency and the like, so that the mutual influence among the participating operation variables is weakened, and the calculation efficiency and the positioning precision are improved. However, locating objects closer to the reference sensorThe thresholds for the TSWLS and CTLS positioning algorithms are significantly higher than the targetThe threshold of (a) is small, and the positioning performance of the positioning algorithm is poor when positioning the target close to the reference sensor.
Fig. 8 depicts the MSE and bias of 2 algorithms for close range target position estimation as a function of sensor position error variance under certain scenarios. As can be seen from the figures, the,while the threshold for the optimization algorithm in this context is-10 dBm2The optimization algorithm is improved by 30dBm compared with the threshold value of the TSWLS algorithm2And the robustness is stronger. In addition, whenThe time context optimization algorithm reduces the estimation bias by about 8dB compared to the TSWLS algorithm. Therefore, the optimization algorithm shows higher positioning accuracy and stronger robustness when positioning a close-range target.
Fig. 9 presents the MSE and bias of 2 algorithms for remote target position estimation as a function of sensor position error variance. As can be seen from the graph, whenIn time, the MSEs for the position estimation of the positioning algorithm in 2 can all reach the CRLB. The noise thresholds for the TSWLS algorithm and the optimization algorithm are-10 dBm, respectively2And-5 dBm2Obviously, the optimization algorithm has stronger robustness than the TSWLS algorithm. In addition, the optimization algorithm is also better thanThe TSWLS algorithm has a small deviation in the position estimate. In conclusion, when the remote target is positioned, the optimization algorithm also shows better positioning performance.
FIG. 11 depicts the MSE and bias of the 2 algorithms on the target position estimate as a function of the variance of the sensor position error when the target and sensor are randomly distributed. According to the curve in the figure, whenOrThe MSE for the target position estimate by the TSWLS algorithm suddenly deviates from CRLB and the deviation also increases dramatically, indicating that the target being located at this time may be just at or near the reference sensor, resulting in an increased estimation error. In contrast, the optimization algorithm herein can reach CRLB at the error threshold, and the optimization algorithm reduces the position estimation MSE by 4dBm and 7dBm respectively compared with the TSWLS algorithm at the error variance point of the two sensor positions, and also reduces the position estimation error by 7dBm and 6dBm respectively, mainly because the optimization algorithm corrects the defects in the TSWLS algorithm by reselecting the reference sensor and rotating the coordinate system. However, the optimization algorithm in the same positioning scenario locates the target a little longer than the TSWLS algorithm, because the optimization algorithm introduces an increase in the amount of computation due to the re-selection of the reference sensor and the rotating coordinate system. Comprehensively, when randomly distributed targets are positioned by using randomly distributed sensors, the overall performance of the optimization algorithm on positioning the targets is superior to that of the TSWLS algorithm.
Fig. 12 shows the average MSE and the variance of the target position estimate with 2 positioning algorithms as a function of the number of sensors, which increased gradually from 3 to 8 in the experiment. At the same time, the average position error variance of the sensor is set toThe settings of other parameters were the same as in the above experiment. The experimental result shows that the MSE and the deviation of the target position estimation of the 2 algorithms are reduced along with the increase of the number of the sensors, particularly the number of the sensors is 3The MSE and bias of the position estimate decrease sharply when the number of sensors increases to 5, but the change will be slower when the number of sensors increases from 6 to 8, which is why 6 sensors are set to locate the target in the above experiment.
And (3) synthesizing the simulation results and analysis, and providing an ICTLS positioning algorithm on the basis of the TSWLS and CTLS positioning algorithms aiming at the positioning problems of a plurality of non-coincident dynamic targets. The ICTLS positioning algorithm provided by the invention corrects two defect problems in the TSWLS algorithm: firstly, the deviation of the WLS estimation result increases along with the increase of noise errors; secondly, the calculation error caused by the nonlinear operation introduced in the WLS in the second step is larger, thereby reducing the precision of the positioning algorithm. The ICTLS algorithm provided in the text modifies the relation between the additional variables and the target position coordinates which are not considered to be introduced in the CTLS algorithm, and establishes global constraint conditions according to the relation between the target and the additional variables; meanwhile, a Lagrange multiplier technology is introduced to solve a positioning equation, calculation of a Hessian matrix is avoided through a BFGS iterative formula of a quasi-Newton method, and the calculated amount is reduced so as to accelerate the convergence speed. Simulation experiments show that under the condition of moderate measurement errors and sensor position errors, the ICTLS positioning algorithm has smaller RMSE than TSWLS and CTLS positioning algorithms, and shows higher positioning accuracy and stronger robustness.
Aiming at a specific scene, the TSWLS optimized positioning algorithm provided by the patent solves the hidden danger problem existing in various positioning algorithms, namely, estimation error is increased suddenly when a target is positioned at or close to a reference sensor. On the basis of keeping a closed-form solution, the optimization algorithm effectively reduces the position estimation error by reselecting the reference sensor and the rotating coordinate system, corrects the defects of the TSWLS algorithm and comprehensively improves the positioning performance of the TSWLS algorithm.
The foregoing is illustrative of the best mode of the invention and details not described herein are within the common general knowledge of a person of ordinary skill in the art. The scope of the present invention is defined by the appended claims, and any equivalent modifications based on the technical teaching of the present invention are also within the scope of the present invention.
Claims (6)
1. A multi-moving object positioning method for improving TDOA/FDOA algorithm is characterized by comprising the following steps:
step 1, setting K non-coincident moving targets in a three-dimensional space, namely a target 1 and a target 2 … …, respectively, deploying M sensors in the three-dimensional space, namely a sensor 1 and a sensor 2 … …, taking the sensor 1 as a reference sensor, and constructing an FDOA measurement vector of speed information of a target i:
wherein, i is 1,2 … … K, j is 1,2 … … M, rj1,iFor the true value of the difference in distance from the object i to the sensor 1 and the sensor j,the true distance between target i and sensor j,is the true distance, Δ r, between the target i and the sensor 1j1,iThe difference value between the actual value and the measured value of the distance difference between the target i and the sensor 1 and the sensor j is an error value;
step 2: constructing a TDOA/FDOA system of positioning equations for object i and sensor location information, introducing additional variablesWherein u isiIs the position of the target i and,for the true position of sensor j, handleAndsubstituted typeAnd squaring the two sides, and obtaining a TDOA equation set after arrangement:
and step 3: will be provided withAt the sensor noise position s1And performing first-order Taylor series expansion, neglecting higher-order terms with more than two orders, entering a TDOA equation set to obtain a new equation set, differentiating the new equation set with time to obtain an FDOA equation set, and synthesizing the new equation set and the FDOA equation set to obtain a positioning equation set in a matrix form: (a + Δ a) θ ═ b + Δ b; in the formula (I), the compound is shown in the specification,data matrixData vector Is r1,iIs determined by the estimated value of (c),is composed ofAn estimated value of (d);
and 4, step 4: taking into account the errors of all coefficient matrices in the equation, and incorporating additional variablesAndand w andconstraint condition between, get constraint condition thetaTΣ1θ=0,θTΣ2θ=0;
And 5: constructing a value function of a vector theta to be solved, and changing the total least square solution of the vector theta for solving the positioning equation set into a CTLS solution;
step 6: introducing Lagrange multiplier technique lambda1And λ2Solving the problem of solving the CTLS constraint optimization minimum value to obtain the final estimated value of the vector theta
θ=[(A-D)TWθA+λ1Σ1+λ2Σ2]-1(A-D)TWθb (18)
Substituting equation (18) into θ againTΣ1θ=0,θTΣ2θ is 0 to yield:
θTΣ1θ=0,θTΣ2θ=0 (19);
where θ is an estimate of θ;
and 7: definition of λ ═ λ1 λ2]T,f=[f1 f2]TSolving λ in equation (19) using a quasi-newton iterative equation:
λ(k+1)=λ(k)-α(k)H(k)g(k) (21)
in the formula, λ(k)For the k-th iteration of λ, α(k)In order to be the step size,H(k)is a corresponding dimensionThe symmetric matrix of the numbers is updated by the BFGS formula as follows:
let parameter delta be (0,1), sigma be (0,0.5), Wθ=Q-1Precision of 10-3Taking an initial point λ(1)=[0,0]TInitial matrix H(1)Is a unit matrix I, k is 0;
and 8: computingIf g | | |(k)Stopping iteration if | | is less than or equal to epsilon, and outputting lambdakAs an approximate solution;
and step 9: calculating a search direction p(k)=-H(k)g(k)And is provided with mkIs satisfied with f (λ)(k))+δmp(k)≤f(λ(k))+σδmg(k)Tp(k)The smallest non-negative integer of (c);
step 10: updating H by BFGS formula of quasi-Newton method(k+1)Setting k to k +1, and turning to step 8;
2. The method for multi-moving object location with improved TDOA/FDOA algorithm of claim 1, wherein said step 1 specifically comprises the steps of:
let the i-th target's true position and speed be recorded as uiAnd1,2, …, K, orderThe real coordinates of the sensor are represented,which is indicative of the true speed of the sensor,and(j 1, 2.. times.m) represents the true position and speed, [. times.]0Is []True value of [ ·]TIs []The transposition operation of (1);
the sensor with j equal to 1 is used as the reference sensor, so the real distance between the target i and the sensor j is
The actual and measured values of the difference in distance from the object i to the sensors 1 and j are respectively
The TDOA-form measurement vector of the target i-position information is recorded as
ri=[r21,i,r31,i,…,rM1,i]T=ri 0+Δri (3)
FDOA measurement of target velocity information:
3. the method for multi-moving object location with improved TDOA/FDOA algorithm as recited in claim 1, wherein said step 3 comprises the steps of:
will be provided withAt the sensor noise position s1And performing first-order Taylor series expansion on the position, neglecting higher-order terms above the second order, and obtaining after arrangement:
differentiating equation (7) with respect to time to obtain the FDOA equation set:
integrating (7) and (8) to obtain a positioning equation system in a matrix form:
(A+ΔA)θ=b+Δb (9)
4. The method for multi-moving object location with improved TDOA/FDOA algorithm as recited in claim 1, wherein said step 4 comprises the following:
taking into account the errors of all coefficient matrices in the equation, and incorporating additional variablesAndand w andestimating target position information by adopting a CTLS positioning algorithm under the constraint condition; the cost function of the candidate vector θ can be configured as:
J(θ)=(Aθ-b)TWθ(Aθ-b) (10)
wherein J (θ) is a cost function of θ, WθIs a weighting matrix; theta is an estimated value of theta, and theta is an optimized value of theta; due to the introduced variablesAndand w andhas a certain constraint relation, so that the constraint strip needs to be combinedSolving equation (11); the constraint conditions are obtained by integrating (11) and (4)
Transforming the formula (12) to obtain the formula (10), and arranging to obtain the constraint condition in a matrix form:
θTΣ1θ=0,θTΣ2θ=0 (13)
5. The method for multi-moving object location with improved TDOA/FDOA algorithm as recited in claim 1, wherein said step 5 comprises the following:
the overall least squares solution for vector θ becomes the CTLS solution:
S.t.θTΣ1θ=0,θTΣ2θ=0 (15)
2[(A-D)TWθA+λ1Σ1+λ2Σ2]θ-2(A-D)TWθb=0 (16)。
6. the method for multi-moving object location with improved TDOA/FDOA algorithm as recited in claim 1, wherein said step 6 comprises the following:
equation (16) is a quadratic constraint non-convex optimization problem, and solving the problem comprises the steps of firstly solving partial differential of equation (16) on theta and then enabling equation (16) to be subjected to partial differentialAfter finishing, obtaining:
2[(A-D)TWθA+λ1Σ1+λ2Σ2]θ-2(A-D)TWθb=0 (17)
in the formula (I), the compound is shown in the specification,l-8 is the dimension of the vector theta,is an initial weighting matrix and n ═ Δ βT,ΔαT]TAs random noise, F1,…,FLA coefficient matrix which is a disturbance matrix delta A, delta b; thus, the final estimate of vector θ is
θ=[(A-D)TWθA+λ1Σ1+λ2Σ2]-1(A-D)TWθb (18)
Substituting equation (18) into constraint equation (13) again yields:
θTΣ1θ=0,θTΣ2θ=0 (19)
formula (19) relates to1And λ2By solving the system of equations to obtain λ1And λ2And substituting the position information into the formula (18) to obtain the final estimated position information of the target i.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110082721.1A CN112904274A (en) | 2021-01-21 | 2021-01-21 | Multi-moving-object positioning method for improving TDOA/FDOA algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110082721.1A CN112904274A (en) | 2021-01-21 | 2021-01-21 | Multi-moving-object positioning method for improving TDOA/FDOA algorithm |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112904274A true CN112904274A (en) | 2021-06-04 |
Family
ID=76118085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110082721.1A Pending CN112904274A (en) | 2021-01-21 | 2021-01-21 | Multi-moving-object positioning method for improving TDOA/FDOA algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112904274A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113691934A (en) * | 2021-08-23 | 2021-11-23 | 武汉科技大学 | UWB indoor positioning method based on LOS credibility recognition |
CN114089273A (en) * | 2021-11-22 | 2022-02-25 | 电子科技大学 | GPS and UWB based motion platform positioning method |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5570099A (en) * | 1993-10-15 | 1996-10-29 | Loral Federal Systems Company | TDOA/FDOA technique for locating a transmitter |
CN1537395A (en) * | 2001-07-18 | 2004-10-13 | ��ʵ��λ��˾ | Improved method for estimating TDOA and FDOA in wireless location system |
CN102308228A (en) * | 2008-12-23 | 2012-01-04 | 泰勒斯公司 | Method for locating multiple rays of a source with or without AOA by multi-channel estimation of the TDOA and FDOA |
CN102331581A (en) * | 2011-05-27 | 2012-01-25 | 哈尔滨工业大学 | Rapid positioning method of binary TDOA/FDOA satellite-to-earth integration positioning system |
CN105116377A (en) * | 2015-08-18 | 2015-12-02 | 西安电子科技大学 | An FDOA positioning method based on an HHT instantaneous energy spectrum |
US20180149729A1 (en) * | 2015-06-17 | 2018-05-31 | Thales | Method and device for locating an electromagnetic emission source and system implementing such a method |
CN108732532A (en) * | 2018-04-19 | 2018-11-02 | 郑州航空工业管理学院 | A kind of generalized rank target is based on the polymorphic fusion and positioning methods of TDOA and FDOA |
CN109633581A (en) * | 2018-12-26 | 2019-04-16 | 杭州电子科技大学 | Based on the localization method under external sort algorithm TDOA/FDOA error correction |
CN111551895A (en) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | Motion radiation source TDOA and FDOA positioning method based on weighted multidimensional scale and Lagrange multiplier technology |
US20200386851A1 (en) * | 2019-06-05 | 2020-12-10 | Rincon Research Corporation | Multi-antenna detection, localization, and filtering of complex time-and-doppler-shifted signals |
CN112114296A (en) * | 2020-09-18 | 2020-12-22 | 王玉冰 | Parameter estimation method and system for unmanned aerial vehicle cooperative TDOA/FDOA composite positioning |
-
2021
- 2021-01-21 CN CN202110082721.1A patent/CN112904274A/en active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5570099A (en) * | 1993-10-15 | 1996-10-29 | Loral Federal Systems Company | TDOA/FDOA technique for locating a transmitter |
CN1537395A (en) * | 2001-07-18 | 2004-10-13 | ��ʵ��λ��˾ | Improved method for estimating TDOA and FDOA in wireless location system |
CN102308228A (en) * | 2008-12-23 | 2012-01-04 | 泰勒斯公司 | Method for locating multiple rays of a source with or without AOA by multi-channel estimation of the TDOA and FDOA |
CN102331581A (en) * | 2011-05-27 | 2012-01-25 | 哈尔滨工业大学 | Rapid positioning method of binary TDOA/FDOA satellite-to-earth integration positioning system |
US20180149729A1 (en) * | 2015-06-17 | 2018-05-31 | Thales | Method and device for locating an electromagnetic emission source and system implementing such a method |
CN105116377A (en) * | 2015-08-18 | 2015-12-02 | 西安电子科技大学 | An FDOA positioning method based on an HHT instantaneous energy spectrum |
CN108732532A (en) * | 2018-04-19 | 2018-11-02 | 郑州航空工业管理学院 | A kind of generalized rank target is based on the polymorphic fusion and positioning methods of TDOA and FDOA |
CN109633581A (en) * | 2018-12-26 | 2019-04-16 | 杭州电子科技大学 | Based on the localization method under external sort algorithm TDOA/FDOA error correction |
US20200386851A1 (en) * | 2019-06-05 | 2020-12-10 | Rincon Research Corporation | Multi-antenna detection, localization, and filtering of complex time-and-doppler-shifted signals |
CN111551895A (en) * | 2020-04-25 | 2020-08-18 | 中国人民解放军战略支援部队信息工程大学 | Motion radiation source TDOA and FDOA positioning method based on weighted multidimensional scale and Lagrange multiplier technology |
CN112114296A (en) * | 2020-09-18 | 2020-12-22 | 王玉冰 | Parameter estimation method and system for unmanned aerial vehicle cooperative TDOA/FDOA composite positioning |
Non-Patent Citations (3)
Title |
---|
闫雷兵: ""基于无线传感器网络的目标定位与跟踪技术研究"", 《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》 * |
闫雷兵等: ""改进的时差频差多目标定位算法"", 《北京邮电大学学报》 * |
陈少昌等: ""传感器位置误差条件下的约束总体最小二乘时差定位算法"", 《航空学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113691934A (en) * | 2021-08-23 | 2021-11-23 | 武汉科技大学 | UWB indoor positioning method based on LOS credibility recognition |
CN113691934B (en) * | 2021-08-23 | 2024-06-04 | 武汉科技大学 | UWB indoor positioning method based on LOS credibility identification |
CN114089273A (en) * | 2021-11-22 | 2022-02-25 | 电子科技大学 | GPS and UWB based motion platform positioning method |
CN114089273B (en) * | 2021-11-22 | 2023-05-26 | 电子科技大学 | GPS and UWB-based motion platform positioning method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109116349B (en) | Multi-sensor cooperative tracking joint optimization decision method | |
CN107703480B (en) | Mixed kernel function indoor positioning method based on machine learning | |
Yan et al. | An improved NLOS identification and mitigation approach for target tracking in wireless sensor networks | |
CN112904274A (en) | Multi-moving-object positioning method for improving TDOA/FDOA algorithm | |
Liu et al. | Smallest enclosing circle-based fingerprint clustering and modified-WKNN matching algorithm for indoor positioning | |
CN109756842A (en) | Wireless indoor location method and system based on attention mechanism | |
CN113794983B (en) | Multi-target indoor positioning method based on nonlinear geometric constraint optimization | |
Yu et al. | AI based location tracking in WiFi indoor positioning application | |
CN111157943A (en) | TOA-based sensor position error suppression method in asynchronous network | |
CN107576932B (en) | Alternating Kalman spatial registration method for coexistence of cooperative target and non-cooperative target | |
CN109862510B (en) | Convex region positioning method based on compressed sensing | |
Wang | Wireless network indoor positioning method using nonmetric multidimensional scaling and RSSI in the Internet of things environment | |
CN110145966A (en) | A kind of finite time convergence control Global robust sliding formwork method of guidance | |
US20240094343A1 (en) | Method, device, system, and storage medium for tracking moving target | |
CN116047495B (en) | State transformation fusion filtering tracking method for three-coordinate radar | |
Assayag et al. | Indoor positioning system using synthetic training and data fusion | |
CN109375160B (en) | Angle measurement error estimation method in pure-azimuth passive positioning | |
CN110649911B (en) | Distributed nonlinear Kalman filtering method based on alpha divergence | |
Li et al. | Coordinate mapping method of crowdsourced data for indoor fingerprint localization | |
Ou et al. | An improved node localization based on adaptive iterated unscented Kalman filter for WSN | |
Li et al. | Performance comparison and evaluation of indoor positioning technology based on machine learning algorithms | |
CN112800889A (en) | Target tracking method based on distributed matrix weighting and Gaussian filtering fusion | |
Zhang et al. | RFID localization system based on k-nearest neighbor algorithm and extreme learning machine algorithm with virtual reference tags | |
CN111983561A (en) | TDOA (time difference of arrival) positioning method for multiple unmanned aerial vehicle targets under receiver position error | |
Guo et al. | A novel density-based geolocation algorithm for a noncooperative radio emitter using power difference of arrival |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210604 |
|
RJ01 | Rejection of invention patent application after publication |