CN1499873A - Method for eveluating position - Google Patents
Method for eveluating position Download PDFInfo
- Publication number
- CN1499873A CN1499873A CNA021483396A CN02148339A CN1499873A CN 1499873 A CN1499873 A CN 1499873A CN A021483396 A CNA021483396 A CN A021483396A CN 02148339 A CN02148339 A CN 02148339A CN 1499873 A CN1499873 A CN 1499873A
- Authority
- CN
- China
- Prior art keywords
- sigma
- nlos
- error
- nlos error
- location estimation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
First, combining recognized result of NLOS in location estimation, using mean of NLOS error determined by probability density function and distribution parameter of error, zero mean rectification is carried out for NLOS error in TDOA measurement. Then, combining recognized result of NLOS, using variance determined by probability density function and distribution parameter of error, weighting matrix in method of location estimation is adjusted so as to restrain NLOS error in preliminary. Finally, averaging estimated results in multiple positions realizes restrain of NLOS error further. The above scheme is applicable to LOS environment as well as applicable to NLOS environment and mixed environment between LOS and NLOS. The method restrains NLOS error in realtime effectively without need of follow-up measuring mobile station for long time.
Description
Technical field
The present invention relates to the location estimation method in radio communication and the radio position finding radio directional bearing field, particularly the location estimation method in the cellular mobile station location.
Background technology
Traditional application in wireless location technology, GPS location as radar location, radiocoustic position finding, vast and bare area, non-visual (NLOS) propagation path is not ubiquitous phenomenon, and the localization method that produces in these fields also all has been based upon visual (LOS) propagation path and has existed on the basis.But in the cellular mobile station navigation system, because blocking or the fluctuating of landform of surface structures, the NLOS of signal propagates becomes a kind of universal phenomenon, and this NLOS error of being introduced by the non-visual propagation path relative time delay of LOS propagation path (promptly with respect to) can cause position estimation accuracy significantly to reduce.
In the location estimation method of existing poor based on the time of advent (TDOA), basically all be under the prerequisite that has line of sight between hypothesis radiation source and the transducer (location receiver), to construct, under this hypothesis prerequisite, the error signal that is input to location estimation method only derives from the TDOA measure error of transducer, and this TDOA measure error is the zero-mean Gaussian Profile.In numerous location estimation methods based on TDOA, Y.T.Chan (Chen Shi) proposition (referring to, a kind of simple and effective position estimator " A simple and efficient estimator forhyperbolic location " that is used for hyperbolic fix, IEEE Trans Signal processing, vol.42, no.8, Aug.1994, pp.1905-1915) typical.The estimator of this method realizes it being to estimate to come the approximate maximal possibility estimation of realizing with weighted linear least square (twice use), theory analysis shows, under the less situation of TDOA measure error, this approximate processing can guarantee that the location estimation performance reaches higher level.In addition, this method can be applicable to that promptly TDOA measures the situation that number equals radiation source coordinate number, is applicable to that also TDOA measures the situation of number greater than radiation source coordinate number; Both be applicable to the situation that linear sensor is arranged, also be applicable to the situation that transducer is arranged arbitrarily.On function, can say that the Chen Shi method has realized the unification of existing methodological function is better than (being not less than at least) other TDOA location estimation method again on performance, its operand significantly is lower than platform labor Series Method, and the performance of estimated accuracy and Tai Lao Series Method is suitable.Above-mentioned Chen Shi method is the same with existing other location estimation methods based on TDOA, shortcoming is: when having the NLOS error, the measure error that is input to the Chen Shi method no longer satisfies the such hypothesis of zero-mean Gaussian Profile, at this moment, the Chen Shi method no longer is the location estimation method of optimizing, and that is to say that the Chen Shi method does not possess the NLOS error and suppresses ability.Therefore, existing location estimation method based on TDOA is because the influence of NLOS error causes position estimation accuracy not high.
Summary of the invention
The object of the present invention is to provide the higher location estimation method of a kind of position estimation accuracy, use this method to suppress of the influence of NLOS error with adaptive mode to positioning accuracy.
For achieving the above object, location estimation method provided by the invention comprises:
Step 1: when carrying out location estimation, utilize the average of non-visual (NLOS) error of obtaining in real time that the NLOS error that comprises in poor (TDOA) measuring amount time of advent is carried out the zero-mean rectification;
Step 2: utilize the variance of the NLOS error of obtaining in real time that the element in the weighting matrix in the weighted least-squares location estimation is carried out the self adaptation adjustment;
Step 3: the weighted least-squares that positions the position is estimated, obtains this location estimation value;
Step 4: the location estimation value of repeatedly obtaining is averaged, get location estimation value to the end.
Described step 1 further comprises:
Step 11: from least one TDOA measuring amount, identify the TDOA measuring amount that comprises the NLOS error respectively;
Step 12: the average of NLOS error is estimated in the TDOA measuring amount that above-mentioned steps 11 is identified;
Step 13: the average of the NLOS error that above-mentioned steps 12 is obtained is carried out the zero-mean of NLOS error and is corrected.
Described step 11 realizes by radiation source being carried out NLOS identification.
Above-mentioned radiation source is carried out NLOS identification, pick out the sample coefficient of dispersion of one group of most powerful path one group of power time delay that utilization is obtained from radiation source distributes and realize NLOS identification, or utilize between the footpath of the single or multiple power time delay of obtaining from radiation source on distributing difference power to realize NLOS identification, or utilize between the footpath of the single or multiple power time delay of obtaining from radiation source on distributing amplitude difference to realize NLOS identification.
Described step 12 further comprises:
Step 121: the distributed constant that from one group of power time delay of the channel of the TDOA measuring amount correspondence that comprises the NLOS error distributes, obtains the NLOS error under discrete form and the conitnuous forms;
Step 122: utilize the distribution form of NLOS error in NLOS error profile parameter and the TDOA measuring amount, calculate the average of NLOS error;
Step 121 is obtained the distributed constant of the NLOS error under the discrete form according to following method:
In the formula: p
iBe NLOS error profile parameter under the discrete form; m
kBe k (k ∈ 1,2 ... the N) number in detected footpath in the individual scattering object statistic window; W is the width of scattering object statistic window, and unit is a chip; N is for obtaining a p
iThe number that distributes of the power time delay that estimated value adopted; α is the sampling number that carries out in the chip.
Step 121 also can be obtained the distributed constant of the NLOS error under the discrete form according to following method:
In the formula: p
iBe NLOS error profile parameter under the discrete shape; s
kBe k (k ∈ 1,2 ... N) the detected number that surpasses the sampling point of detection threshold in the individual scattering object statistic window; W is the width of scattering object statistic window, and unit is a sampling point; N is for obtaining a p
iThe number that distributes of the power time delay that estimated value adopted.
Step 121 is obtained the distributed constant of the NLOS error under the conitnuous forms according to following method:
In the formula, θ
iBe the distributed constant of the NLOS error under the conitnuous forms, T is systematic sampling sampling point blanking time, p
iBe NLOS error profile parameter under the discrete form.
Described step 2 further comprises:
Step 21: the form of determining weighting matrix;
Step 22: determine the element that needs are adjusted in the weighting matrix according to the TOA measuring amount that comprises the NLOS error of step 11;
Step 23: need in the weighting matrix to determine the element value adjusted.
Step 21 is determined the weighting matrix of following form:
Q
(r)=Q
(l)+Q
(n)
Wherein,
Q
(l)Be TDOA measure error N=[n under the LOS environment
2,1n
3,1N
M, 1]
TCovariance matrix, σ
I, j 2Covariance for TDOA measure error under the LOS environment;
Q
(n)Be the NLOS error n after the zero-mean rectification
I, jThat constitute is M-1 n dimensional vector n N
(n)=[n
(n) 2,1n
(n) 3,1N
(n) M, 1]
TCovariance matrix.
Described step 23 further comprises:
Step 231: determine that each need adjust the value of element in the oblique variance matrix of TDOA measure error under visual in the weighting matrix (LOS) environment: be similar to by the statistics of under the LOS channel circumstance system being carried out the TDOA measure error and obtain, perhaps obtain by system emulation is approximate;
Step 232: determine that each need adjust the value of element in the oblique variance matrix of the multidimensional vector that the NLOS error after zero-mean in the weighting matrix is corrected constitutes: the probability density function by the NLOS error in the TOA measuring amount that comprises the NLOS error obtains.
Because location estimation method of the present invention can suppress the NLOS error, make this method widen the suitable environment of existing TDOA location estimation method, make it not to be only applicable to the LOS environment, and be applicable to NLOS environment and LOS and NLOS mixed environment mutually; Simultaneously,, can suppress the NLOS error in real time, need not carry out tracking measurement to travelling carriage for a long time because NLOS error mean and variance that the present invention uses can obtain in real time; Therefore, the invention solves because the influence of NLOS error causes the not high problem of position estimation accuracy.
Description of drawings
Fig. 1 is the embodiment flow chart of the method for the invention;
Fig. 2 is the NLOS error zero-mean correcting process flow chart that Fig. 1 adopts;
Fig. 3 is that the self adaptation that Fig. 1 adopts is adjusted the weighting matrix process flow diagram.
Embodiment
In order to suppress NLOS error in the TDOA measuring amount to the influence of position estimation accuracy, the basic ideas that the present invention adopts are: under the prerequisite of the average of obtaining the NLOS error and variance, the average of using the NLOS error is the stochastic variable of NLOS error rectification as zero-mean, use the variance structure weighted least-squares of the NLOS error weighting matrix in estimating tentatively to suppress NLOS error (average of NLOS error at this moment has been zero) then to the position estimation effect, last again according to the zero mean characteristic of the NLOS error after correcting, by multiple averaging, further suppress the NLOS error to the position estimated result.
The present invention is described in further detail below in conjunction with accompanying drawing.
Fig. 1 is the embodiment flow chart of the method for the invention.Four basic steps that shown in Figure 1 is with mobile position estimation embodiment:
Step 101 utilizes the average of the NLOS error of being obtained by the NLOS error profile parameter of obtaining in real time that the NLOS error that comprises in the TDOA measuring amount is carried out the zero-mean rectification.This step is made up of substep 201,202,203 as shown in Figure 2.
Step 201 identifies the TDOA measuring amount TDOA that comprises the NLOS error
(m) i, j, TDOA
(m) i, jExpression is through the zero-mean of NLOS error TDOA measuring amount that correct, original, and subscript i, j represent that this TDOA measuring amount is that time of advent between i radiation source (as the base station in the cellular network) and j the radiation source is poor;
Wherein, TDOA
(los) i, j (0)It is the desirable TDOA value that does not comprise any error; n
(n) i, jBe the residual error of the NLOS error after correcting through zero-mean, its average is zero; μ
(n) i, jIt is the average of the NLOS error that goes out according to NLOS error profile calculation of parameter; n
I, jBeing the TDOA time delay evaluated error of system under the LOS environment, is the stochastic variable (not considering that clock drift is to the influence of TDOA measuring amount between each radiation source) of zero-mean normal distribution.In step 201, judge TDOA
(m) i, jThe method that whether comprises the NLOS error can be by carrying out the method realization that NLOS discerns to i radiation source (as the base station in the cellular network) and j radiation source respectively.The method of i radiation source being carried out NLOS identification has multiple, as, can utilize one group of power time delay obtaining from i radiation source distribute (as, the pseudo noise code of utilizing this radiation source to transmit is extracted) in pick out the sample coefficient of dispersion of one group of most powerful path (each power time delay distribute on pick out a most powerful path) size realize NLOS identification, can utilize also that power (or amplitude) difference realizes that NLOS discerns between the footpath of the single or multiple power time delay of obtaining from i radiation source on distributing.
Step 202 is carried out the estimation of the average of NLOS error in the TDOA measuring amount.NLOS recognition result according to step 201, at first, from distributing, one group of power time delay of the channel of the TDOA measuring amount correspondence that comprises the NLOS error obtains the distributed constant of NLOS error under discrete form and the conitnuous forms, during specific implementation, both can adopt following formula (1), and also can adopt following formula (2) to finish NLOS error profile parameter p under the discrete form
iAnd p
jEstimation (for LOS channel ch
i, its NLOS error profile parameter p
iBe 1).
In the formula: p
iBe NLOS error profile parameter under the discrete form; m
kBe k (k ∈ 1,2 ... the N) number in the interior detected footpath of individual scattering object statistic window (from the intercepting that distributes of k power time delay, the starting point of scattering object statistic window can be certain position after the first footpath, also can comprise first footpath, first footpath); W is the width of scattering object statistic window, and unit is a chip, and usually, the value of W is between 1~10 chip; N is for obtaining a p
iThe number that distributes of the power time delay that estimated value adopted, usually, the value of N is between 1~10, used N power time delay distribution is to carry out N Multipath searching to obtain in the regular hour interval; α is the sampling number that carries out in the chip, and usually, α is value between 1~32, and the α value is exactly a number of samples that the footpath comprises.
In the formula: p
iBe NLOS error profile parameter under the discrete form; s
kBe k (k ∈ 1,2 ... the N) number of detected sampling point above detection threshold in the individual scattering object statistic window (from the intercepting that distributes of k power time delay, the starting point of scattering object statistic window can be certain position after the first footpath, also can comprise first footpath, first footpath); W is the width of scattering object statistic window, and unit is a sampling point, and usually, the value of W is within 40 sampling points; N is for obtaining a p
iThe number that distributes of the power time delay that estimated value adopted, usually, the value of N is between 1~10, used N power time delay distribution is to carry out N Multipath searching to obtain in the regular hour interval.
Then, utilize NLOS error profile parameter p under the above-mentioned discrete form
iAnd p
jAnd the distribution form of NLOS error, calculate the average of NLOS error.The average of the NLOS error of discrete form can be utilized p
i, p
jNLOS error delta with the TDOA of the discrete form of following formula (3) expression
(s) ijProbability density function directly obtain.
In the formula, δ
(s) ijBe i and j base station correspondence be the NLOS margin of error of the TDOA of unit with the sampling point number, unit is the sampling point number, δ
(s) ij ∈ (...-3 ,-2 ,-1,0,1,2,3 ...), this number of samples and sampling point product at interval is exactly NLOS error (a time dimension); p
i, p
jBe respectively the NLOS margin of error δ of the TOA measurement of i and j base station
(s) i, δ
(s) jDistributed constant; δ
(s) ij=δ
(s) i-δ
(s) j
The average of the NLOS error of conitnuous forms can be utilized p
i, p
jObtain the distributed constant θ of the NLOS error in the TDOA measuring amount of conitnuous forms with following formula (4)
iAnd θ
j, utilize θ then
i, θ
jNLOS error delta with the TDOA of the conitnuous forms of following formula (5) expression
I, jProbability density function obtain the average of NLOS error.
In the formula, θ
iBe the distributed constant of the NLOS error under the conitnuous forms, T is systematic sampling sampling point blanking time, p
iBe NLOS error profile parameter under the discrete form.
Obtain θ with reference to above-mentioned formula (4)
j Work as δ
I, j<0 is δ
i<δ
jThe time (5a)
Step 203 is utilized the average μ of the NLOS error that step 202 obtains
(n) i, j, (6) carry out the zero-mean rectification of NLOS error according to the following equation.
In the formula, TDOA
I, j (nlos_miti)The TDOA measuring amount after NLOS error zero-mean is corrected is carried out in expression.
Step 102 utilizes the variance of the NLOS error of being obtained in real time by NLOS error profile parameter that the element in the weighting matrix in the weighted least-squares location estimation is carried out the self adaptation adjustment.This step comprises substep shown in Figure 3 301,302 and 303.
Step 301 is determined the form of weighting matrix, the weighting matrix that this example adopts following formula (7) to estimate as the least square location of mobile station.
Q
(r)=Q
(l)+Q
(n)
Wherein,
Q
(l)Be TDOA measure error N=[n under the LOS environment
2,1n
3,1N
M, 1]
TCovariance matrix, N is the M-1 n dimensional vector n, the average of N is zero, the covariance matrix Q of N
(l)Be (M-1) * (M-1) dimension symmetrical matrix.σ
I, j 2Covariance for TDOA measure error under the LOS environment.
Q
(n)Be the NLOS error n after the zero-mean rectification
I, jThat constitute is M-1 n dimensional vector n N
(n)=[n
(n) 2,1n
(n) 3,1N
(n) M, 1]
TCovariance matrix, N
(n)Average be zero, Q
(n)Be (M-1) * (M-1) dimension symmetrical matrix;
Step 302 need in the weighting matrix to determine the element adjusted, and this step is directly utilized the result of step step 201, as long as constitute TDOA measuring amount TDOA
(m) i, jI radiation source (as the base station in the cellular network) and j radiation source among have one to comprise the NLOS error, will adjust element (weight coefficient) corresponding in the weighting matrix.As, through NLOS identification, determine that the 1st radiation source (corresponding channel) is the LOS channel, then, Q in the formula (7)
(n)In comprise σ
(n) 1 2Element, all to adjust (this moment σ
(n) 1 2Get null value);
Step 303 need in the weighting matrix to determine the value of the element adjusted, and this substep is again by two parts content:
1, at first determines matrix Q
(l)In the value of each element, the σ in the promptly above-mentioned formula (8)
I, j 2Or σ
i 2Value.σ
I, j 2And σ
i 2The pass be
As long as determined σ
I, j 2And σ
i 2In one, just can determine another one, as, can be similar to by the statistics of under the LOS channel circumstance system being carried out the TDOA measure error and obtain σ
I, j 2, also can be similar to and obtain σ by system emulation
I, j 2, then, if necessary, just can utilize
Determine σ
i 2
2, determine matrix Q
(n)In the value of each element.Q
(n)Have formula (9a) and (9b) two kinds of forms, corresponding two kinds of different definite Q
(n)The method of element value.
For Q
(n)The form of formula (9a), only need ask for the variance of the NLOS error in the TOA measuring amount, the p that this can utilize step 202 to ask for
iNLOS error delta with the TOA of following formula (10) expression
(s) iProbability density function obtain σ
(n) i 2, the also p that can utilize step 202 to ask for
iAnd the NLOS error delta of the TOA of formula (4), following formula (11) expression
iProbability density function obtain σ
(n) i 2
f
δ(δ
(s) iδ is worked as in)=0
(s) iGet (0,1,2 ...) during in addition value (10)
δ in the formula
(s) iThe number of samples of expression systematic sampling, this number of samples and sampling point product at interval is exactly NLOS error (a time dimension), p
iDistributed constant for geometric distributions.
f
δ(δ
iδ is worked as in)=0
iWhen getting other value (11)
For Q
(n)The form of formula (9b) calculate σ except needs
(n) i 2Outside, also to calculate σ
(n) i, j 2σ
(n) i, j 2σ
(n) i 2Calculating can be according to Q
(n)The discussion of formula (9a) carry out; σ
(n) i, j 2Calculating just need to utilize the p that asks for by step 202
i, p
jAnd the probability density function of the NLOS error of the discrete form in the TDOA measuring amount of formula (3) expression, perhaps use the p that asks for by step 202
i, p
jp
j, obtain θ by formula (4)
iAnd θ
jAfter, utilize the probability density function of the NLOS error of the conitnuous forms in the TDOA measuring amount of formula (5) expression to calculate σ again
(n) i, j 2
Step 103 is weighted the least square location estimation, obtains one group of estimated value that has tentatively suppressed the location of mobile station of NLOS error effect.This step at first is transformed into the distance relation of following equation (12) expression the Hyperbolic Equation group of formula (13) expression,
In the formula:
(x
i, y
i) be the position coordinates parameter of known base station, (x y) is the location of mobile station coordinate parameter of unknown (waiting to find the solution).r
iBe the distance of travelling carriage to i base station.According to the propagation velocity of signal between travelling carriage and base station is c, can obtain following one group of Hyperbolic Equation:
r
i,1=cd
i,1=r
i-r
1 (13)
In the formula: d
I, 1The expression travelling carriage arrives time delay poor of first base station to the time delay of i base station and travelling carriage; r
I, 1The expression travelling carriage to i base station apart from r
iWith travelling carriage to first base station apart from r
1Poor.
Solving equation (13) just can obtain travelling carriage position coordinates (x, y).The basic step of solving equation (13) is: 1) introduce intermediate variable d
1=(x-x
1)
2+ (y-y
1)
22) equation (13) twice use weighted least-squares are estimated.
When utilizing weighted least-squares to estimate to find the solution this Nonlinear System of Equations twice, the Q that use step 102 to determine
(r)Replace the estimation error covariance matrix of equation group (13), the location estimation the when radiation source (base station or satellite, base station and satellite) that so just can obtain four or more participates in the location.What obtain here is the weighted least-squares estimated value that has tentatively suppressed the location of mobile station of NLOS error effect.
Step 104 averages by one group of location of mobile station value that step 103 is obtained, and further suppresses the influence of NLOS error to position estimation accuracy.The specific implementation step is: at first by step 103, in short time interval, obtain the location estimation value that a plurality of (greater than) export on independent TDOA based measurement, then, these location estimation values (coordinate components) are averaged, obtain a position coordinates.This step utilization is zero characteristics according to the average that formula (6) carries out the NLOS error of the TDOA measured value after zero-mean is corrected, average by a plurality of location estimation results that step 103 is obtained, further suppress the NLOS error and correct the influence of residual error to position estimation accuracy, variance of estimaion error dips.
Claims (11)
1, a kind of location estimation method comprises:
Step 1: when carrying out location estimation, utilize the average of non-visual (NLOS) error of obtaining in real time that the NLOS error that comprises in poor (TDOA) measuring amount time of advent is carried out the zero-mean rectification;
Step 2: utilize the variance of the NLOS error of obtaining in real time that the element in the weighting matrix in the weighted least-squares location estimation is carried out the self adaptation adjustment;
Step 3: the weighted least-squares that positions the position is estimated, obtains this location estimation value;
Step 4: the location estimation value of repeatedly obtaining is averaged, get location estimation value to the end.
2, location estimation method according to claim 1 is characterized in that, described step 1 further comprises:
Step 11: from least one TDOA measuring amount, identify the TDOA measuring amount that comprises the NLOS error respectively;
Step 12: the average of NLOS error is estimated in the TDOA measuring amount that above-mentioned steps 11 is identified;
Step 13: the average of the NLOS error that above-mentioned steps 12 is obtained is carried out the zero-mean of NLOS error and is corrected.
3, location estimation method according to claim 2 is characterized in that, described step 11 realizes by radiation source being carried out NLOS identification.
4, location estimation method according to claim 3, it is characterized in that: described radiation source is carried out NLOS identification, pick out the sample coefficient of dispersion of one group of most powerful path one group of power time delay that utilization is obtained from radiation source distributes and realize NLOS identification, or utilize between the footpath of the single or multiple power time delay of obtaining from radiation source on distributing difference power to realize NLOS identification, or utilize between the footpath of the single or multiple power time delay of obtaining from radiation source on distributing amplitude difference to realize NLOS identification.
5, location estimation method according to claim 2 is characterized in that, described step 12 further comprises:
Step 121: the distributed constant that from one group of power time delay of the channel of the TDOA measuring amount correspondence that comprises the NLOS error distributes, obtains the NLOS error under discrete form and the conitnuous forms;
Step 122: utilize the distribution form of NLOS error in NLOS error profile parameter and the TDOA measuring amount, calculate the average of NLOS error;
6, location estimation method according to claim 5 is characterized in that, step 121 is obtained the distributed constant of the NLOS error under the discrete form according to following method:
In the formula: p
iBe NLOS error profile parameter under the discrete form; m
kBe k (k ∈ 1,2 ... the N) number in detected footpath in the individual scattering object statistic window; W is the width of scattering object statistic window, and unit is a chip; N is for obtaining a p
iThe number that distributes of the power time delay that estimated value adopted; α is the sampling number that carries out in the chip.
7, location estimation method according to claim 5 is characterized in that, step 121 is obtained the distributed constant of the NLOS error under the discrete form according to following method:
In the formula: p
iBe NLOS error profile parameter under the discrete shape; s
kBe k (k ∈ 1,2 ... N) the detected number that surpasses the sampling point of detection threshold in the individual scattering object statistic window; W is the width of scattering object statistic window, and unit is a sampling point; N is for obtaining a p
iThe number that distributes of the power time delay that estimated value adopted.
According to claim 6 or 7 described location estimation methods, it is characterized in that 8, step 121 is obtained the distributed constant of the NLOS error under the conitnuous forms according to following method:
In the formula, θ
iBe the distributed constant of the NLOS error under the conitnuous forms, T is systematic sampling sampling point blanking time, p
iBe NLOS error profile parameter under the discrete form.
9, location estimation method according to claim 2 is characterized in that, described step 2 further comprises:
Step 21: the form of determining weighting matrix;
Step 22: determine the element that needs are adjusted in the weighting matrix according to the TOA measuring amount that comprises the NLOS error of step 11;
Step 23: need in the weighting matrix to determine the element value adjusted.
10, location estimation method according to claim 9 is characterized in that, step 21 is determined the weighting matrix of following form:
Q
(r)=Q
(l)+Q
(n)
Wherein,
Q
(l)Be TDOA measure error N=[n under the LOS environment
2,1n
3,1N
M, 1]
TCovariance matrix, σ
I, j 2Covariance for TDOA measure error under the LOS environment;
Q
(n)Be the NLOS error n after the zero-mean rectification
I, jThat constitute is M-1 n dimensional vector n N
(n)=[n
(n) 2,1n
(n) 3,1N
(n) M, 1]
TCovariance matrix.
11, location estimation method according to claim 9 is characterized in that, described step 23 further comprises:
Step 231: determine that each need adjust the value of element in the oblique variance matrix of TDOA measure error under visual in the weighting matrix (LOS) environment: be similar to by the statistics of under the LOS channel circumstance system being carried out the TDOA measure error and obtain, perhaps obtain by system emulation is approximate;
Step 232: determine that each need adjust the value of element in the oblique variance matrix of the multidimensional vector that the NLOS error after zero-mean in the weighting matrix is corrected constitutes: the probability density function by the NLOS error in the TOA measuring amount that comprises the NLOS error obtains.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02148339 CN1243453C (en) | 2002-11-08 | 2002-11-08 | Method for eveluating position |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02148339 CN1243453C (en) | 2002-11-08 | 2002-11-08 | Method for eveluating position |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1499873A true CN1499873A (en) | 2004-05-26 |
CN1243453C CN1243453C (en) | 2006-02-22 |
Family
ID=34233134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 02148339 Expired - Fee Related CN1243453C (en) | 2002-11-08 | 2002-11-08 | Method for eveluating position |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1243453C (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102170658A (en) * | 2011-04-28 | 2011-08-31 | 北京交通大学 | Geometric positioning improvement method under NLOS (non-line-of-sight) environment |
CN101449306B (en) * | 2006-03-21 | 2011-11-02 | 天米公司 | Private, auditable vehicle positioning system and on-board unit for same |
CN103105778A (en) * | 2013-02-05 | 2013-05-15 | 中原工学院 | Estimation method for industrial process simulation mathematical model parameters |
CN103869345A (en) * | 2014-01-01 | 2014-06-18 | 深圳北斗国芯科技有限公司 | United time synchronizing and positioning method for Beidou satellite navigation positioning system |
CN107241797A (en) * | 2017-06-02 | 2017-10-10 | 重庆邮电大学 | Mono-station location method based on scattering object information under NLOS environment |
CN108375754A (en) * | 2018-02-10 | 2018-08-07 | 沈阳理工大学 | Node positioning method based on mobile node original state and mobile status in WSN |
CN113973262A (en) * | 2020-07-24 | 2022-01-25 | 华为技术有限公司 | Method for single anchor point positioning under multipath and communication device |
-
2002
- 2002-11-08 CN CN 02148339 patent/CN1243453C/en not_active Expired - Fee Related
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101449306B (en) * | 2006-03-21 | 2011-11-02 | 天米公司 | Private, auditable vehicle positioning system and on-board unit for same |
CN102170658A (en) * | 2011-04-28 | 2011-08-31 | 北京交通大学 | Geometric positioning improvement method under NLOS (non-line-of-sight) environment |
CN102170658B (en) * | 2011-04-28 | 2013-10-30 | 北京交通大学 | Geometric positioning improvement method under NLOS (non-line-of-sight) environment |
CN103105778A (en) * | 2013-02-05 | 2013-05-15 | 中原工学院 | Estimation method for industrial process simulation mathematical model parameters |
CN103105778B (en) * | 2013-02-05 | 2015-07-15 | 中原工学院 | Estimation method for industrial process simulation mathematical model parameters |
CN103869345A (en) * | 2014-01-01 | 2014-06-18 | 深圳北斗国芯科技有限公司 | United time synchronizing and positioning method for Beidou satellite navigation positioning system |
CN103869345B (en) * | 2014-01-01 | 2016-03-30 | 深圳北斗国芯科技有限公司 | For associating time synchronized and the localization method of Beidou satellite navigation and positioning system |
CN107241797A (en) * | 2017-06-02 | 2017-10-10 | 重庆邮电大学 | Mono-station location method based on scattering object information under NLOS environment |
CN107241797B (en) * | 2017-06-02 | 2019-06-14 | 重庆邮电大学 | Based on the mono-station location method of scatterer information under NLOS environment |
CN108375754A (en) * | 2018-02-10 | 2018-08-07 | 沈阳理工大学 | Node positioning method based on mobile node original state and mobile status in WSN |
CN113973262A (en) * | 2020-07-24 | 2022-01-25 | 华为技术有限公司 | Method for single anchor point positioning under multipath and communication device |
CN113973262B (en) * | 2020-07-24 | 2023-05-02 | 华为技术有限公司 | Method for positioning single anchor point under multipath and communication device |
Also Published As
Publication number | Publication date |
---|---|
CN1243453C (en) | 2006-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1235429C (en) | Method for estimating position | |
KR100791179B1 (en) | Location calculating method and location estimating apparatus | |
CN109975749B (en) | Short-wave single-station direct positioning method under condition of existence of correction source | |
CN1582600A (en) | System and method for identification of transmitters with limited information | |
CN1131644C (en) | Method and system for comparing measured radio frequency signal propagation characteristics in wireless communication system | |
CN1488076A (en) | Method and apparatus for location estimation | |
CN102480743B (en) | Positioning method of cross-cell coverage, apparatus and system thereof | |
CN1541338A (en) | Improvements to tracking systems | |
CN1243453C (en) | Method for eveluating position | |
US20030054845A1 (en) | Enhanced time of arrival estimation using reduced complexity optimal processing | |
CN109975842A (en) | A kind of blind catching method of Big Dipper satellite signal high-precision based on wavelet transformation | |
CN1470885A (en) | Angle evaluating method for restraining multi-path influence | |
EP1256814A3 (en) | GPS method and apparatus, navigation system, and computer program | |
CN1788211A (en) | A method for the location of mobile terminals | |
US8188920B2 (en) | Location measurement acquisition optimization with Monte Carlo simulation | |
CN111199281A (en) | Short wave single-station direct positioning deviation compensation method based on geographical coordinate airspace position spectrum | |
CN1235430C (en) | Method for eveluating position | |
US10241188B2 (en) | Method and apparatus for obtaining time of arrival TOA when mobile terminal is located | |
CN1256853C (en) | Method for obtaining mean and variance of timedelay error for non visibility paths | |
US10151823B2 (en) | Method for passive approximate localization using frequency modulation and software defined radio | |
US8063824B2 (en) | Magnitude-squared coherence filtering in a time difference of arrival (TDOA) or group delay interferometric geolocation system | |
CN1885746A (en) | Doppler frequency detector and doppler frequency estimation method | |
US8436773B2 (en) | Method for leveraging diversity for enhanced location determination | |
KR20160044244A (en) | A positioning apparatus and a positioning method | |
CN1592434A (en) | Method of estimating arrival time additing delay error |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20060222 Termination date: 20191108 |