CN107422302A - A kind of method for improving cellphone subscriber's displacement estimated accuracy - Google Patents
A kind of method for improving cellphone subscriber's displacement estimated accuracy Download PDFInfo
- Publication number
- CN107422302A CN107422302A CN201710351667.XA CN201710351667A CN107422302A CN 107422302 A CN107422302 A CN 107422302A CN 201710351667 A CN201710351667 A CN 201710351667A CN 107422302 A CN107422302 A CN 107422302A
- Authority
- CN
- China
- Prior art keywords
- subsequence
- base station
- displacement
- deviation
- track
- 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
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/14—Determining absolute distances from a plurality of spaced points of known location
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Mobile Radio Communication Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Navigation (AREA)
Abstract
The invention discloses a kind of method for improving cellphone subscriber's displacement estimated accuracy, belongs to displacement detection field, and the present invention determines to need all base station subsequences rectified a deviation first;Then estimate an advance track for current subsequence, determine the beginning and end of the advance track;Then replace current subsequence with the beginning and end of track and judge whether that also subsequence is untreated;Return and continue with when also subsequence is untreated, finally according to the accumulation displacement of the sequence estimation user after replacement;The present invention compared with prior art, can be directed to linear motion and carry out more accurate user's displacement estimation with the estimation precision of raising user's displacement of high degree.
Description
Technical field
The invention belongs to displacement detection field, and cellphone subscriber's displacement accuracy of detection is improved more particularly to one kind
Method.
Background technology
Telecom operators or user mobile phone can record the base station data passed through during user movement, and these data can be helped
Us are helped to estimate the displacement of user.When user advances along road, the base station of record is typically distributed on the both sides of road.
Now, if directly describing the movement locus of user using the line of base station location, the estimation of displacement will produce larger
Deviation.
The content of the invention
In view of the drawbacks described above of prior art, the technical problems to be solved by the invention are to provide a kind of mobile phone that improves and used
The method of family displacement estimated accuracy.
To achieve the above object, the invention provides a kind of method for improving cellphone subscriber's displacement estimated accuracy, bag
Include following steps:
Step 1, determine to need all base station subsequences rectified a deviation;
Step 2, it is that current subsequence estimates an advance track;
Step 3, the beginning and end for determining the advance track;
Step 4, the current subsequence of beginning and end replacement with track;
Step 5, judge whether that also subsequence is untreated;Returned when also subsequence is untreated and perform step 2, otherwise
Perform step 6;
Step 6, the accumulation displacement according to the sequence estimation user after replacement;
All base station subsequences for determining to need to rectify a deviation are carried out according to the following steps:
Step 101. determines the starting point and end point of subsequence;Starting point is the pending 1st in all base station sequences
Individual base station, end point are the 2nd pending base station;
Step 102. calculates the air line distance and Cumulative Distance of subsequence;
If S=(s1, s2..., sn) it is base station sequence, calculate l (S)=d (s1, sn) obtain S air line distance l (S);
Wherein d (a, b) represents the Euclidean distance of a points and b points, if a and b coordinate is respectively (x1, y1),(x2, y2), then
Have:
CalculateObtain S accumulation displacement a (S);The n is positive integer;
Step 103. determines to need the base station subsequence rectified a deviation;
If the Cumulative Distance of base station subsequence is more than threshold value T1 and Cumulative Distance and the ratio of direct range are more than threshold value
T2, then the subsequence need to rectify a deviation;
If step 104. also base station is untreated, step 101 is returned to, is otherwise terminated.
Base station subsequence (the s that the step 2 is rectified a deviation to needs1, s2..., sn), user is obtained using linear regression algorithm
Movement locus straight line.
The starting point that projects as track of the step 3 using base station subsequence starting point on the movement locus straight line,
Terminal using projection of the end point on the movement locus straight line as track.
The beneficial effects of the invention are as follows:The present invention can be with the estimation precision of raising user's displacement of high degree, phase
Than prior art, linear motion can be directed to and carry out more accurate user's displacement estimation.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of the embodiment of the invention.
Fig. 2 is the schematic diagram using linear regression algorithm estimation user movement track.
Fig. 3 is to determine the schematic diagram of track beginning and end.
Fig. 4 is that user's displacement estimates schematic diagram in the prior art.
Fig. 5 is scatter diagram when calculating movement locus straight line.
Embodiment
The invention will be further described with reference to the accompanying drawings and examples:
As shown in figure 4, when user advances along road, the base station of record is typically distributed on the both sides of road.Existing skill
Art directly describes the movement locus of user using the line of base station location, and the estimation of user's displacement as shown in Figure 4 produces
Relatively large deviation.
As shown in Figure 1, Figure 2 and Figure 3, the invention discloses a kind of side for improving cellphone subscriber's displacement estimated accuracy
Method, comprise the following steps:
Step 1, determine to need all base station subsequences rectified a deviation;
Step 2, it is that current subsequence estimates an advance track;
Step 3, the beginning and end for determining the advance track;
Step 4, the current subsequence of beginning and end replacement with track;
Step 5, judge whether that also subsequence is untreated;Returned when also subsequence is untreated and perform step 2, otherwise
Perform step 6;
Step 6, according to the sequence after replacement, calculateCome estimate the accumulation movement of user away from
From;
All base station subsequences for determining to need to rectify a deviation are carried out according to the following steps:
Step 101. determines the starting point and end point of subsequence;Starting point is the pending 1st in all base station sequences
Individual base station, end point are the 2nd pending base station;
Step 102. calculates the air line distance and Cumulative Distance of subsequence;
If S=(s1, s2..., sn) it is base station sequence, calculate l (S)=d (s1, sn) obtain S air line distance l (S);
Wherein d (a, b) represents the Euclidean distance of a points and b points, if a and b coordinate is respectively (x1, y1),(x2, y2), then
Have:
CalculateObtain S accumulation displacement a (S);The n is positive integer;
Step 103. determines to need the base station subsequence rectified a deviation;
If the Cumulative Distance of base station subsequence is more than threshold value T1 and Cumulative Distance and the ratio of direct range are more than threshold value
T2, then the subsequence need to rectify a deviation;If in general, a base station subsequence needed to rectify a deviation, the sequence and next base
The sequence formed of standing may be also required to rectify a deviation.Therefore, when finding that needs the subsequence rectified a deviation, attempt subsequent base stations not
It is disconnected to be added in subsequence, until Cumulative Distance and the ratio of direct range are less than or equal to threshold value T2.Now by last
Base station removes from subsequence, and marks the subsequence to need the base station subsequence rectified a deviation.
If step 104. also base station is untreated, step 101 is returned to, is otherwise terminated.
Base station subsequence (the s that the step 2 is rectified a deviation to needs1, s2..., sn), user is obtained using linear regression algorithm
Movement locus straight line.
As shown in figure 5, in statistics, linear regression (Linear Regression) is to utilize referred to as linear regression side
A kind of regression analysis that relation is modeled between the least square function pair one or more independent variable and dependent variable of journey.
According to the step of above, we have drawn the base station subsequence (s that needs are rectified a deviation1, s2..., sn), Si(i=1,
2 ... n) are determined by longitude and latitude, therefore SiCan be by coordinate pair (xi,yi) represent.Each point is put into plane right-angle coordinate,
The coordinate system represents its specific physical location.We can obtain Fig. 5, and this figure is referred to as scatter diagram.From figure 5 it can be seen that
Data point falls near straight line substantially, and this teaches that variable X and Y relation substantially can be regarded as linear relationship, i.e.,
Correlation between them can be described with linear relationship.But due to being not that all data points entirely fall within one directly
On line, therefore X and Y relation does not have the definite degree for arriving and Y value can uniquely being determined by an X value.In a practical situation,
The particular location of user can be influenced by many factors.If we will study X and Y relation, linear fit can be made
Formula (1) is regression equation, and a and b is from undetermined constant, referred to as regression coefficient.Theoretically, formula (1) has
Infinite more piece, the task of regression analysis are to obtain its optimal linear fit.Here, we ask a, b from least square method
Value.
Principle of least square method:If using regression equationIt is calculated | γ | i values (i=1,2 ...
N) it is referred to as regressand value, then actual measured value yiWith regressand valueBetween there is deviation, this deviation is called residual error, remembered
For ei(i=1,2 ... n).Then the extent of deviation of measured value and regression straight line is measured with residual sum of squares (RSS).Residual sum of squares (RSS)
It is defined as:
Determine that a and b value make it that Q (a, b) is minimum, i.e., with the regression straight line of least square methodIt is in institute
Have and minimum one of value residual sum of squares (RSS) Q is predicted in straight line.By formula (2) it is known that Q is on a, b quadratic function,
So what its minimum value was constantly present.It can be seen from the method for extreme value is sought in differential, Q (a, b)
Obtaining minimum value should meet
By formula (2) and above-mentioned condition, then
Formula (4) is referred to as normal equation group.Solving this equation group can obtain
Wherein
LxyReferred to as xy covariance sum, LxxReferred to as x difference of two squares sum.
Arrive this, it has been determined that a, b value, also determined that straight lineAlso the movement locus of cellphone subscriber has been determined that.
In the present embodiment, the step 3 using projection of the base station subsequence starting point on the movement locus straight line as
The starting point of track, the terminal using projection of the end point on the movement locus straight line as track.
Preferred embodiment of the invention described in detail above.It should be appreciated that one of ordinary skill in the art without
Creative work can is needed to make many modifications and variations according to the design of the present invention.Therefore, all technologies in the art
Personnel are available by logical analysis, reasoning, or a limited experiment on the basis of existing technology under this invention's idea
Technical scheme, all should be in the protection domain being defined in the patent claims.
Claims (3)
- A kind of 1. method for improving cellphone subscriber's displacement estimated accuracy, it is characterised in that comprise the following steps:Step 1, determine to need all base station subsequences rectified a deviation;Step 2, it is that current subsequence estimates an advance track;Step 3, the beginning and end for determining the advance track;Step 4, the current subsequence of beginning and end replacement with track;Step 5, judge whether that also subsequence is untreated;Returned when also subsequence is untreated and perform step 2, otherwise performed Step 6;Step 6, the accumulation displacement according to the sequence estimation user after replacement;All base station subsequences for determining to need to rectify a deviation are carried out according to the following steps:Step 101. determines the starting point and end point of subsequence;Starting point is the 1st base pending in all base station sequences Stand, end point is the 2nd pending base station;Step 102. calculates the air line distance and Cumulative Distance of subsequence;If S=(s1, s2..., sn) it is base station sequence, calculate l (S)=d (s1, sn) obtain S air line distance l (S);Wherein d (a, b) represents the Euclidean distance of a points and b points, if a and b coordinate is respectively (x1, y1),(x2, y2), then have:CalculateObtain S accumulation displacement a (S);The n is positive integer;Step 103. determines to need the base station subsequence rectified a deviation;If the Cumulative Distance of base station subsequence is more than threshold value T1 and Cumulative Distance and the ratio of direct range are more than threshold value T2, The subsequence needs to rectify a deviation;If step 104. also base station is untreated, step 101 is returned to, is otherwise terminated.
- 2. a kind of method for improving cellphone subscriber's displacement estimated accuracy as claimed in claim 1, it is characterized in that:The step The rapid 2 couples base station subsequence (s for needing to rectify a deviation1, s2..., sn), the movement locus that user is obtained using linear regression algorithm is straight Line.
- 3. a kind of method for improving cellphone subscriber's displacement estimated accuracy as claimed in claim 2, it is characterized in that:The step Rapid 3 starting point projected as track using base station subsequence starting point on the movement locus straight line, by end point described Terminal of the projection as track on movement locus straight line.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710351667.XA CN107422302B (en) | 2017-05-18 | 2017-05-18 | Method for improving mobile distance estimation precision of mobile phone user |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710351667.XA CN107422302B (en) | 2017-05-18 | 2017-05-18 | Method for improving mobile distance estimation precision of mobile phone user |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107422302A true CN107422302A (en) | 2017-12-01 |
CN107422302B CN107422302B (en) | 2020-12-01 |
Family
ID=60425245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710351667.XA Expired - Fee Related CN107422302B (en) | 2017-05-18 | 2017-05-18 | Method for improving mobile distance estimation precision of mobile phone user |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107422302B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454060A (en) * | 2023-12-22 | 2024-01-26 | 季华实验室 | Linear motion device position measuring method and device, electronic equipment and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020080064A1 (en) * | 2000-12-26 | 2002-06-27 | Lg Electronics Inc. | System for determining position of mobile communication terminal and method thereof |
CN101600257A (en) * | 2009-06-30 | 2009-12-09 | 中兴通讯股份有限公司 | The WIFI local system of supervisory control system, method for supervising and supervisory control system |
CN104828698A (en) * | 2015-04-15 | 2015-08-12 | 华东师范大学 | Ad Hoc network-based crane automatic cruising system of Beidou positioning system, and method thereof |
CN105743750A (en) * | 2016-03-23 | 2016-07-06 | 北京绿源普惠科技有限公司 | Adaptive sensing method for smart home |
CN106014058A (en) * | 2016-06-30 | 2016-10-12 | 李小司 | Deviation-correcting method for railless retractable door based on GPS positioning technology and deviation-correcting device for railless retractable door based on GPS positioning technology |
-
2017
- 2017-05-18 CN CN201710351667.XA patent/CN107422302B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020080064A1 (en) * | 2000-12-26 | 2002-06-27 | Lg Electronics Inc. | System for determining position of mobile communication terminal and method thereof |
CN101600257A (en) * | 2009-06-30 | 2009-12-09 | 中兴通讯股份有限公司 | The WIFI local system of supervisory control system, method for supervising and supervisory control system |
CN104828698A (en) * | 2015-04-15 | 2015-08-12 | 华东师范大学 | Ad Hoc network-based crane automatic cruising system of Beidou positioning system, and method thereof |
CN105743750A (en) * | 2016-03-23 | 2016-07-06 | 北京绿源普惠科技有限公司 | Adaptive sensing method for smart home |
CN106014058A (en) * | 2016-06-30 | 2016-10-12 | 李小司 | Deviation-correcting method for railless retractable door based on GPS positioning technology and deviation-correcting device for railless retractable door based on GPS positioning technology |
Non-Patent Citations (1)
Title |
---|
徐娟: ""蜂窝无线网络中手机定位算法的研究"", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454060A (en) * | 2023-12-22 | 2024-01-26 | 季华实验室 | Linear motion device position measuring method and device, electronic equipment and storage medium |
CN117454060B (en) * | 2023-12-22 | 2024-03-19 | 季华实验室 | Linear motion device position measuring method and device, electronic equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN107422302B (en) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10492022B2 (en) | System and method for robust and accurate RSSI based location estimation | |
CN105704652B (en) | Fingerprint base acquisition and optimization method in a kind of positioning of WLAN/ bluetooth | |
CN103533650A (en) | Cosine-similarity-based indoor positioning method capable of improving positioning precision | |
CN103916821A (en) | Floor distinguishing method based on RSSI difference between floors | |
CN108307301A (en) | Indoor orientation method based on RSSI rangings and track similitude | |
CN109151707B (en) | Sight distance/non-sight distance path identification method in moving state | |
CN105242239A (en) | Indoor subarea positioning method based on crowdsourcing fingerprint clustering and matching | |
CN103471586B (en) | The terminal combinations localization method that a kind of sensor is auxiliary and device | |
CN104080165A (en) | Indoor wireless sensor network positioning method based on TDOA | |
CN103379441A (en) | Indoor positioning method based on region segmentation and curve fitting | |
CN100415042C (en) | Method of evaluating a location of a mobile station within a cellular telecommunication network | |
CN114594495A (en) | Positioning deviation rectifying method and device | |
TW201329485A (en) | Wireless communication positioning method | |
CN103888979A (en) | Indoor positioning method based on wireless local area network | |
CN107422302A (en) | A kind of method for improving cellphone subscriber's displacement estimated accuracy | |
CN101201398B (en) | Method for positioning indirect wave using auxiliary measuring angle of adjacent station ranging | |
CN105116373B (en) | Target IP region city-class positioning algorithm based on indirect time delay | |
KR101593679B1 (en) | Method and device for estimating indoor location | |
CN106125070A (en) | A kind of nanoLOC range measurement exceptional value removing method | |
He et al. | A hybrid data fusion based cooperative localization approach for cellular networks | |
CN110944293B (en) | Radio search positioning method based on path attenuation and Kalman filtering fusion | |
CN107889211A (en) | A kind of method and device of positioning | |
CN107228667A (en) | A kind of improved Kalman filter device indoor positioning tracking for merging cartographic information | |
CN106063317B (en) | The localization method and device of traffic hotspots | |
CN112800889A (en) | Target tracking method based on distributed matrix weighting and Gaussian filtering fusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
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: 20201201 Termination date: 20210518 |