CN111031474B - Method for predicting longitude and latitude of base station based on user MDT data - Google Patents

Method for predicting longitude and latitude of base station based on user MDT data Download PDF

Info

Publication number
CN111031474B
CN111031474B CN201911263121.4A CN201911263121A CN111031474B CN 111031474 B CN111031474 B CN 111031474B CN 201911263121 A CN201911263121 A CN 201911263121A CN 111031474 B CN111031474 B CN 111031474B
Authority
CN
China
Prior art keywords
latitude
user
base station
longitude
users
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911263121.4A
Other languages
Chinese (zh)
Other versions
CN111031474A (en
Inventor
王计斌
宋铁成
张程
张华�
胡静
余健
朱格苗
杜利涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Howso Technology Co ltd
Southeast University
Original Assignee
Nanjing Howso Technology Co ltd
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Howso Technology Co ltd, Southeast University filed Critical Nanjing Howso Technology Co ltd
Priority to CN201911263121.4A priority Critical patent/CN111031474B/en
Publication of CN111031474A publication Critical patent/CN111031474A/en
Application granted granted Critical
Publication of CN111031474B publication Critical patent/CN111031474B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02—Services making use of location information
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/08—Testing, supervising or monitoring using real traffic

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The invention relates to a method for predicting longitude and latitude of a base station based on user MDT data, which comprises the following steps: (1) collecting user MDT data; (2) removing longitude and latitude abnormal values of the user, (3) converting the longitude and latitude of the user into a two-dimensional plane coordinate; (4) calculating an initial iteration point (x)0,y0) (ii) a (5) Converting the time advance Ta of the user into a distance; (6) performing secondary screening according to the iteration points in the step (5) for iterative users; (7) calculating a matrix G and a vector b according to the initial point coordinates of a given base station; (8) then, a displacement vector delta x is calculated according to the matrix G and the vector b, if the length of the displacement vector ||Δx||<ε, turning to the step (10), otherwise, executing the step (9); (9) setting k to k +1, returning to the step (7), and circulating the step (7) and the step (8) until ||Δx||<ε; (10) obtaining an estimated value x of the coordinates of the base station, converting the two-dimensional plane coordinates x into the longitude and latitude, and obtaining a predicted value of the longitude and latitude of the base station
Figure DDA0002312116630000011

Description

Method for predicting longitude and latitude of base station based on user MDT data
Technical Field
The invention relates to the technical field of wireless communication, in particular to a method for predicting longitude and latitude of a base station based on user MDT data.
Background
In recent years, the latitude and longitude of a base station is the most basic parameter of the base station, and is also the reference position covered by different cells under the base station, which has an important role in network planning. The accuracy of the latitude and longitude information recorded by the base station working reference table not only depends on the precision of a measuring instrument, but also has the problem of data recording errors. In addition, the longitude and latitude recorded by the base station are planned positions before station building, the station building positions of part of the base stations are not consistent with the recorded positions, and the work parameter table is not updated. The problems occur, so that the longitude and latitude information of a large part of base stations in the working parameter data of the current base stations has problems. The base station is the core and carrier of mobile communication construction, and LTE site planning generally needs to satisfy utilization maximization, coverage maximization, benefit maximization and customer satisfaction maximization. Base stations are generally divided into two types, one is an omni-base station, i.e., a signal can be transmitted to 360 degrees. In another type, the directional base station transmits signals in a certain direction, and since 360 degrees of direction needs to be covered as required, the directional base station generally includes three base station antennas.
The traditional base station working parameters mainly depend on a manual measurement method, and are compared with historical records to confirm whether the working parameters have deviation or not. Due to the fact that the number of the base stations is too large, if the working parameters of all the base stations are checked, consumed manpower and material resources are too large. In addition, studies have shown that the time difference of arrival TDOA location principle can be used to predict the location of a base station using the time of arrival of signals at different users. However, this requires knowledge of the location of the user and the absolute time difference of arrival of the base station signal at different users.
The Minimization of Drive Test (MDT) technology is an automatic Drive test technology introduced in the LTE system by 3GPP at stage R10. And the base station issues related measurement configuration to the terminal according to the MDT measurement task configured by the network management, and the terminal performs measurement and reports measurement information when meeting measurement conditions. Compared with a conventional Measurement Report (MR), the MDT data includes, in addition to RSRP of the user, longitude and latitude of the user, Ta, and neighbor information. Although the MDT data is large in volume and easy to obtain, part of the MDT data has problems and needs to be discarded due to high-speed movement of the terminal and delay in data uploading. In addition, due to the complexity of the wireless propagation environment and the uneven distribution of users, the MDT data under different base stations is very different.
Therefore, it is necessary to develop a method for predicting the latitude and longitude of the base station based on the MDT data of the user, so as to improve the accuracy and stability of the prediction of the latitude and longitude of the base station and effectively reduce the error introduced by the coordinate conversion.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method for predicting the longitude and latitude of the base station based on user MDT data, so that the prediction precision and stability of the longitude and latitude of the base station are improved, and errors introduced by coordinate conversion are effectively reduced.
In order to solve the technical problems, the invention adopts the technical scheme that: the method for predicting the latitude and longitude of the base station based on the MDT data of the user comprises the following steps:
(1) collecting user MDT data including longitude theta and latitude of user
Figure GDA0002933055610000024
Time advance TaAnd a reference signal received strength, RSRP;
(2) removing longitude and latitude abnormal values of a user, carrying out first screening on longitude and latitude of the user, wherein the number of the users left after the first screening is M;
(3) converting the longitude and latitude of the M users obtained in the step (2) into a two-dimensional plane coordinate by adopting a mercator projection method;
(4) calculating an initial iteration point (x)0,y0)；
(5) Converting the time advance Ta of the user into a distance;
(6) performing secondary screening according to the iteration points in the step (5) for iterative users to obtain K users meeting the conditions;
(7) calculating a matrix G and a vector b according to the initial point coordinates of a given base station;
(8) calculating a displacement vector delta x according to the matrix G and the vector b, if the length of the displacement vector is | | delta x | | < epsilon, wherein epsilon is a converged positioning error, turning to the step (10), and otherwise, executing the step (9);
(9) setting k as k +1, returning to the step (7), and circulating the step (7) and the step (8) until | | delta x | < epsilon;
(10) obtaining an estimated value x of the coordinates of the base station, converting the two-dimensional plane coordinates x into longitude and latitude by adopting a reverse ink carto projection method, and obtaining the longitude and latitude of the current base station
Figure GDA0002933055610000021
As a preferred technical scheme of the present invention, the step (2) of screening the latitude and longitude of the user according to the 3-Sigma criterion specifically comprises: let the number of users under the base station be N, and the longitude of the i (i ═ 1, 2.., N) th user be θiLatitude is
Figure GDA0002933055610000022
Calculating the mean of all user longitudes
Figure GDA0002933055610000023
Standard deviation of
Figure GDA0002933055610000031
Mean value of all user latitudes
Figure GDA0002933055610000032
And standard deviation of
Figure GDA0002933055610000033
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure GDA0002933055610000034
and keeping the user i, otherwise, eliminating the user, wherein the number of the users left after screening is M.
As a preferred embodiment of the present invention, the steps described above(3) The specific process of converting the longitude and latitude of the M users obtained in the step (2) into a two-dimensional plane coordinate by adopting the ink card holder projection method comprises the following steps: taking the reference latitude as
Figure GDA0002933055610000035
The major radius of the earth is a, the minor radius is b, and the first eccentricity of the earth is
Figure GDA0002933055610000036
The latitude of the reference has a weft circle radius of
Figure GDA0002933055610000037
Equivalent latitude of ith user
Figure GDA0002933055610000038
Then the plane coordinate U of the ith useri(xi,yi) Is xi＝r0θi，yi＝r0qiWherein i ═ 1, 2.., M.
As a preferred technical scheme of the invention, the specific process of the step (4) is as follows: finding out that the time advance Ta is 0 and the RSRP is greater than P0If the number of users is less than n, selecting Ta is less than or equal to 1 and RSRP is greater than P1If the number of users is still less than n, selecting Ta less than or equal to 2 and RSRP greater than P2Until the number of users is more than or equal to n, taking the average value of M users meeting the conditions to calculate the initial point B (x) of the base station coordinate0,y0)。
As a preferable technical solution of the present invention, the step (5) is to select the user TaConverting into a distance r from the i (i ═ 1, 2.., M) th user to the base stationi＝78*Tai+39 of wherein TaiIndicating the timing advance of the ith user.
As a preferred technical solution of the present invention, the step (6) specifically comprises:
user i (i ═ 1, 2.., M) to base station initiation point B (x)0,y0) A distance of
Figure GDA0002933055610000039
Screening for the second time to screen out the product satisfying di-ri|≤d0Of d, wherein d0For the distance difference threshold, K users satisfying the condition are obtained.
As a preferred technical solution of the present invention, the step (7) specifically comprises:
s71 setting the initial point coordinate of the base station as
Figure GDA0002933055610000041
Setting k to 0;
s72 substituting the K users obtained in the step (6) to calculate vectors
Figure GDA0002933055610000042
Sum matrix
Figure GDA0002933055610000043
Obtaining a displacement vector Δ x ═ GTG)-1GTb。
As a preferred technical solution of the present invention, the step (8) specifically comprises: if the length of the displacement vector
Figure GDA0002933055610000044
And (5) showing that the coordinates of the base station converge to the given precision, turning to the step (10), and otherwise, executing the step (9).
As a preferred technical solution of the present invention, the step (9) specifically comprises: updating xk+1＝xk+ Δ x, i.e. x(k+1)＝x(k)+Δx，y(k+1)＝y(k)And + Δ y, setting k to k +1, returning to step (7), and looping step (7) and step (8).
As a preferable technical scheme of the invention, the step (10) comprises the following steps:
s101 the estimated value of the plane coordinate of the base station is
Figure GDA0002933055610000045
Taking initial point of latitude of base station
Figure GDA0002933055610000046
Setting k to 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure GDA0002933055610000051
Wherein
Figure GDA0002933055610000052
Where e and r0The first eccentricity of the earth and the reference latitude in the step (2) are respectively
Figure GDA0002933055610000053
The weft circle radius of (1);
s103 if
Figure GDA0002933055610000054
If δ is the allowable latitude error, which indicates that the latitude calculation has converged, step S104 is executed, otherwise, k is set to k +1, and the process returns to step S103;
s104 obtains the longitude forecast value of the base station as
Figure GDA0002933055610000055
Predicted value of latitude is
Figure GDA0002933055610000056
Compared with the prior art, the invention has the following advantages and beneficial effects:
1) because the influence of abnormal values is eliminated by adopting twice screening on the longitude and latitude of the user, the base station longitude and latitude prediction precision is high;
2) because the conversion of the longitude and latitude and the plane coordinate is carried out under the same reference latitude, the invention effectively reduces the error introduced by the coordinate conversion;
3) due to the fact that according to TaThe obtained initial point of the base station coordinate is very close to the actual coordinate, and the convergence rate of the base station coordinate prediction algorithm is high;
4) user latitude and longitude, time advance T in MDT data onlyaAnd referenceUnder the condition of signal receiving strength RSRP, the invention can quickly converge and calculate to obtain the predicted value of the longitude and latitude of the base station, and the method can realize that the prediction error is within 20 meters.
Drawings
The technical scheme of the invention is further described by combining the accompanying drawings as follows:
FIG. 1 is a schematic diagram of base station coordinate prediction of the prediction method of base station latitude and longitude based on user MDT data according to the present invention;
fig. 2 is a base station longitude and latitude prediction flow chart in the base station longitude and latitude prediction method based on the user MDT data of the present invention.
Detailed Description
For the purpose of enhancing the understanding of the present invention, the present invention will be described in further detail with reference to the accompanying drawings and examples, which are provided for the purpose of illustration only and are not intended to limit the scope of the present invention.
Example 1: the method for predicting the latitude and longitude of the base station based on the user MDT data, as shown in figure 2, comprises the following steps:
(1) collecting user MDT data including longitude theta, latitude of user
Figure GDA0002933055610000061
Time advance TaReference signal received strength, RSRP;
(2) removing longitude and latitude abnormal values of a user, carrying out first screening on longitude and latitude of the user, wherein the number of the users left after the first screening is M; the step (2) is to screen the latitude and longitude of the user according to a 3-Sigma criterion, and specifically comprises the following steps: let the number of users under the base station be N, and the longitude of the i (i ═ 1, 2.., N) th user be θiLatitude is
Figure GDA0002933055610000062
Calculating the mean of all user longitudes
Figure GDA0002933055610000063
Standard deviation of
Figure GDA0002933055610000064
Mean value of all user latitudes
Figure GDA0002933055610000065
And standard deviation of
Figure GDA0002933055610000066
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure GDA0002933055610000067
if not, the user i is removed, and the number of the users left after screening is M;
(3) converting the longitude and latitude of the M users obtained in the step (2) into a two-dimensional plane coordinate by adopting a mercator projection method; the specific process of converting the longitude and latitude of the M users obtained in the step (2) into a two-dimensional plane coordinate by adopting the Mocha projection method in the step (3) is as follows: taking the reference latitude as
Figure GDA0002933055610000068
The major radius of the earth is a, the minor radius is b, and the first eccentricity of the earth is
Figure GDA0002933055610000069
The latitude of the reference has a weft circle radius of
Figure GDA00029330556100000610
Equivalent latitude of ith user
Figure GDA00029330556100000611
Then the plane coordinate U of the ith useri(xi,yi) Is xi＝r0θi，yi＝r0qiWherein i 1, 2.. said, M;
(4) calculating an initial iteration point (x)0,y0): finding out that the time advance Ta is 0 and the RSRP is greater than P0If the number of users is less than n, selecting Ta is less than or equal to 1 and RSRP is greater than P1If the number of users is still less than n,then Ta ≦ 2 and RSRP greater than P is selected2Until the number of users is more than or equal to n, taking the average value of M users meeting the conditions to calculate the initial point B (x) of the base station coordinate0,y0)；
(5) Advancing the time T of the useraConverted to a distance, i.e. the distance r from the i (i ═ 1, 2.., M) th user to the base stationi＝78*Tai+39 of wherein TaiRepresenting the time advance of the ith user;
(6) and (5) performing secondary screening according to the iteration points in the step (5) for iterative users to obtain K users meeting the conditions: user i (i ═ 1, 2.. m.,) to the base station initiation point B (x)0,y0) A distance of
Figure GDA0002933055610000071
Screening for the second time to screen out the product satisfying di-ri|≤d0Of d, wherein d0Obtaining K users meeting the condition for the distance difference threshold;
(7) calculating a matrix G and a vector b according to the initial point coordinates of a given base station;
s71 setting the initial point coordinate of the base station as
Figure GDA0002933055610000072
Setting k to 0;
s72 substituting the K users obtained in the step (6) to calculate vectors
Figure GDA0002933055610000073
Sum matrix
Figure GDA0002933055610000074
Obtaining a displacement vector Δ x ═ GTG)-1GTb；
(8) Then, a displacement vector delta x is calculated according to the matrix G and the vector b, if the length of the displacement vector is | | delta x | | survival of the naked eye<E, turning to the step (10), otherwise, executing the step (9); the method specifically comprises the following steps: if the length of the displacement vector
Figure GDA0002933055610000075
Showing that the coordinates of the base station converge to the given precision, turning to the step (10), otherwise, executing the step (9);
(9) setting k to k +1, returning to the step (7), and circulating the step (7) and the step (8) until | Δ x | | survival rate is reached<Epsilon; the step (9) is specifically as follows: updating xk+1＝xk+ Δ x, i.e. x(k+1)＝x(k)+Δx，y(k+1)＝y(k)+ Δ y, set k to k +1, return to step (7), loop through steps (7) and (8) until | | | Δ x | | luminance<ε；
(10) Obtaining the value of the coordinates of the base station, converting the two-dimensional plane coordinates x into the longitude and latitude by adopting a reverse ink carto projection method, and obtaining the longitude and latitude of the current base station
Figure GDA0002933055610000081
S101 the estimated value of the plane coordinate of the base station is
Figure GDA0002933055610000082
Taking initial point of latitude of base station
Figure GDA0002933055610000083
Setting k to be 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure GDA0002933055610000084
Wherein
Figure GDA0002933055610000085
Where e and r0The first eccentricity of the earth and the reference latitude in the step (2) are respectively
Figure GDA0002933055610000086
The weft circle radius of (1);
s103 if
Figure GDA0002933055610000087
The latitude calculation is converged, step S104 is executed, otherwise, the step is returned toStep S103;
s104 obtains the longitude forecast value of the base station as
Figure GDA0002933055610000088
Predicted value of latitude is
Figure GDA0002933055610000089
The invention is an effective method for realizing the longitude and latitude prediction of a base station by utilizing the longitude and latitude of a user, the time advance Ta and the reference signal receiving strength RSRP on the basis that the MDT data of the user can be acquired under the base station, and a schematic diagram of the coordinate prediction of the base station is shown in figure 1;
the specific steps of predicting according to the base station coordinates shown in fig. 1 are:
(1) survey station longitude of base station A is thetabs118.641220 degrees, survey station latitude is
Figure GDA00029330556100000810
Inputting MDT data of all users in 3 cells 1,2 and 3 of the base station, wherein the MDT data mainly comprises longitude theta and latitude theta of the users
Figure GDA00029330556100000811
A timing advance Ta and a reference signal received strength RSRP;
(2) and screening the latitude and longitude of the user according to a 3-Sigma criterion, wherein the specific operation is as follows: the number of users in all cells under the base station is N32247, wherein the 1 st user has a longitude of 118.640823 ° and a latitude of 32247%
Figure GDA00029330556100000815
Longitude of an i (i) ═ 1, 2., N ═ 32247) th user is θiLatitude is
Figure GDA00029330556100000812
Calculating the mean of all user longitudes
Figure GDA00029330556100000813
Standard deviation of
Figure GDA00029330556100000814
Mean value of all user latitudes
Figure GDA0002933055610000091
And standard deviation of
Figure GDA0002933055610000092
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure GDA0002933055610000093
keeping the user i, otherwise, removing the user, wherein the number of the users left after screening is M-31867;
(3) converting the longitude and latitude of the 31867 users obtained in the step (2) into two-dimensional plane coordinates by the mercator projection method: reference latitude is taken
Figure GDA0002933055610000094
The long radius a of the earth is 6378137, the short radius b is 6356752, and the first eccentricity of the earth is
Figure GDA0002933055610000095
Radius of latitude of reference
Figure GDA0002933055610000096
The equivalent latitude of the i (i 1, 2.., M31867) th user
Figure GDA0002933055610000097
Then the plane coordinate of the ith user is xi＝r0θi＝5788686.8θi，yi＝r0qi＝5788686.8qiThe plane coordinates U of all users can be obtainedi(xi,yi) And take θi＝θbs＝118.641220°，
Figure GDA0002933055610000098
Equivalent weft of survey base station can be obtainedDegree qi＝qbs0.44767272, calculating the plane coordinates of the survey base station
Figure GDA0002933055610000099
ybs＝r0qbs＝5788686.8×0.44767272≈2591437.2；
(4) Finding Ta 0 and RSRP greater than P0-85dBm user, if the number of users is less than n-10, then select Ta ≦ 1 and RSRP greater than P1-95dBm user, if the number of users is still less than n 10, then pick Ta ≦ 2 and RSRP greater than P2Until the number of users is 10 or more, the average value of m 1135 users satisfying the condition is calculated as the initial point B (x) of the base station coordinate for the user of-110 dBm0,y0) Is taken immediately
Figure GDA00029330556100000910
(5) The user Ta is converted into a distance, i.e., the distance r from the i-th (i-1, 2.., M-31867) user to the base stationi＝78*Tai+39 of wherein TaiRepresenting the time advance of the ith user;
(6) user i (i 1, 2.., M31867) is located at the initial point of the base station by the distance of
Figure GDA0002933055610000101
Screening out the satisfying | di-ri|≤d050 user ═ 50
Figure GDA0002933055610000105
Distance threshold d0Taking 50 meters, so that 6447 users meeting the conditions can be obtained;
(7) the initial point coordinates of a given base station are:
x0＝(x(0),y(0))＝(x0,y0) (11986496.9,2591424.8) with k equal to 0;
substituting K into 6447 users in step (6), and calculating
Figure GDA0002933055610000102
Sum matrix
Figure GDA0002933055610000103
Obtain the displacement vector Δ x ═ (G)TG)-1GTb；
(8) If the length of the displacement vector
Figure GDA0002933055610000104
Showing that the coordinates of the base station converge to the given precision, turning to the step (10), otherwise, executing the step (9);
(9) updating xk+1＝xk+ Δ x, i.e. x(k+1)＝x(k)+Δx，y(k+1)＝y(k)+ Δ y, where k is k +1, and returning to step (7);
(10) the converged base station coordinates are predicted as x ═ x, y ═ 11986499.0,2591436.4, and the base station coordinates prediction error is:
Figure GDA0002933055610000111
converting the two-dimensional plane coordinate x into longitude and latitude by using a reverse ink cartoost projection method to obtain a predicted value of the longitude and latitude of the current base station as
Figure GDA0002933055610000112
S101 the estimated value of the plane coordinate of the base station is
Figure GDA0002933055610000113
Taking initial point of latitude of base station
Figure GDA0002933055610000114
Setting k to be 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure GDA0002933055610000115
Wherein
Figure GDA0002933055610000116
S103 if
Figure GDA0002933055610000117
Explaining the convergence of the latitude calculation, turning to step S104, otherwise, setting k to k +1, and returning to step S103;
s104 obtains the longitude forecast value of the base station as
Figure GDA0002933055610000118
Predicted value of latitude is
Figure GDA0002933055610000119
Unit is rad, then by
Figure GDA00029330556100001110
And
Figure GDA00029330556100001111
obtaining the longitude and latitude predicted value of the base station with the unit of degree
Figure GDA00029330556100001112
It is obvious to those skilled in the art that the present invention is not limited to the above embodiments, and it is within the scope of the present invention to adopt various insubstantial modifications of the method concept and technical scheme of the present invention, or to directly apply the concept and technical scheme of the present invention to other occasions without modification.

Claims (9)

1. A method for predicting the latitude and longitude of a base station based on user MDT data is characterized by comprising the following steps:
(1) collecting user MDT data, packetsIncluding longitude theta and latitude of user
Figure FDA0002933055600000011
A timing advance Ta and a reference signal received strength RSRP;
(2) removing longitude and latitude abnormal values of a user, carrying out first screening on longitude and latitude of the user, wherein the number of the users left after the first screening is M;
(3) converting the longitude and latitude of the M users obtained in the step (2) into a two-dimensional plane coordinate by adopting a mercator projection method;
(4) calculating an initial iteration point (x)0,y0)；
(5) Converting the time advance Ta of the user into a distance;
(6) performing secondary screening according to the iteration points in the step (4) for iterative users to obtain K users meeting the conditions;
(7) calculating a matrix G and a vector b according to the initial point coordinates of a given base station;
(8) calculating a displacement vector delta x according to the matrix G and the vector b, if the length of the displacement vector is | | delta x | | < epsilon, wherein epsilon is an allowable positioning error, turning to the step (10), and otherwise, executing the step (9);
(9) setting k as k +1, returning to the step (7), and circulating the step (7) and the step (8) until | | delta x | < epsilon;
(10) obtaining an estimated value x of the coordinates of the base station, converting the two-dimensional plane coordinates x into longitude and latitude by adopting a reverse ink carto projection method, and obtaining the longitude and latitude of the current base station
Figure FDA0002933055600000012
The specific process of the step (4) is as follows: finding out that the time advance Ta is 0 and the RSRP is greater than P0If the number of users is less than n, selecting Ta is less than or equal to 1 and RSRP is greater than P1If the number of users is still less than n, selecting Ta less than or equal to 2 and RSRP greater than P2Until the number of users is more than or equal to n, taking the average value of M users meeting the conditions to calculate the initial point B (x) of the base station coordinate0,y0)。
2. The method for predicting the latitude and longitude of the base station based on the user MDT data according to claim 1, wherein the step (2) is to screen the latitude and longitude of the user according to the 3-Sigma criterion, and specifically comprises: let the number of users under the base station be N and the longitude of the ith user be thetaiLatitude is
Figure FDA0002933055600000013
Where i is 1,2, … …, N, the mean of all user longitudes is calculated
Figure FDA0002933055600000021
Standard deviation of
Figure FDA0002933055600000022
Mean value of all user latitudes
Figure FDA0002933055600000023
And standard deviation of
Figure FDA0002933055600000024
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure FDA0002933055600000025
and keeping the user i, otherwise, eliminating the user, wherein the number of the users left after screening is M.
3. The method for predicting the latitude and longitude of the base station based on the MDT data of the user according to claim 2, wherein the step (3) of converting the latitude and longitude of the M users obtained in the step (2) into the two-dimensional plane coordinates by the mercator projection method comprises the following specific steps: taking the reference latitude as
Figure FDA0002933055600000026
The major radius of the earth is a, the minor radius is b, and the first eccentricity of the earth is
Figure FDA0002933055600000027
The latitude of the reference has a weft circle radius of
Figure FDA0002933055600000028
Equivalent latitude of ith user
Figure FDA0002933055600000029
Then the plane coordinate U of the ith useri(xi,yi) Is xi＝r0θi，yi＝r0qiWherein i ═ 1, 2.., M.
4. The method of claim 3, wherein step (5) predicts the user T based on the base station latitude and longitude of the user MDT dataaConverting into a distance r from the ith user to the base stationi＝78*Tai+39, where i ═ 1,2, … …, M, TaiIndicating the timing advance of the ith user.
5. The method for predicting the latitude and longitude of the base station based on the MDT data of the user according to claim 4, wherein the step (6) is specifically as follows:
user i, i ═ 1,2, … …, M to base station initiation point B (x)0,y0) A distance of
Figure FDA00029330556000000210
Screening for the second time to screen out the product satisfying di-ri|≤d0Of d, wherein d0For the distance difference threshold, K users satisfying the condition are obtained.
6. The method for predicting the latitude and longitude of the base station based on the MDT data of the user according to claim 5, wherein the step (7) is specifically as follows:
s71 setting the initial point coordinate of the base station as
Figure FDA00029330556000000211
Setting k to 0;
s72 substituting the K users obtained in the step (6) to calculate vectors
Figure FDA0002933055600000031
Sum matrix
Figure FDA0002933055600000032
Obtaining a displacement vector Δ x ═ GTG)-1GTb。
7. The method for predicting the latitude and longitude of the base station based on the MDT data of the user according to claim 6, wherein the step (8) is specifically as follows: if the length of the displacement vector
Figure FDA0002933055600000033
And (5) showing that the coordinates of the base station converge to the given precision, turning to the step (10), and otherwise, executing the step (9).
8. The method for predicting the latitude and longitude of the base station based on the MDT data of the user according to claim 7, wherein the step (9) is specifically as follows: updating xk+1＝xk+ Δ x, i.e. x(k+1)＝x(k)+Δx，y(k+1)＝y(k)And + Δ y, setting k to k +1, returning to step (7), and looping step (7) and step (8).
9. The method of claim 8, wherein the step (10) comprises the steps of:
s101 the estimated value of the plane coordinate of the base station is
Figure FDA0002933055600000034
Taking initial point of latitude of base station
Figure FDA0002933055600000035
Setting k to 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure FDA0002933055600000036
Where e and r0The first eccentricity of the earth and the reference latitude in the step (2) are respectively
Figure FDA0002933055600000037
The weft circle radius of (1);
s103 if
Figure FDA0002933055600000041
If δ is the allowable latitude error, which indicates that the latitude calculation has converged, step S104 is executed, otherwise, k is set to k +1, and the process returns to step S103;
s104 obtains the longitude forecast value of the base station as
Figure FDA0002933055600000042
Predicted value of latitude is
Figure FDA0002933055600000043
CN201911263121.4A 2019-12-11 2019-12-11 Method for predicting longitude and latitude of base station based on user MDT data Active CN111031474B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911263121.4A CN111031474B (en) 2019-12-11 2019-12-11 Method for predicting longitude and latitude of base station based on user MDT data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911263121.4A CN111031474B (en) 2019-12-11 2019-12-11 Method for predicting longitude and latitude of base station based on user MDT data

Publications (2)

Publication Number Publication Date
CN111031474A CN111031474A (en) 2020-04-17
CN111031474B true CN111031474B (en) 2021-05-07

Family

ID=70205617

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911263121.4A Active CN111031474B (en) 2019-12-11 2019-12-11 Method for predicting longitude and latitude of base station based on user MDT data

Country Status (1)

Country Link
CN (1) CN111031474B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114095968B (en) * 2020-08-24 2024-03-26 中国电信股份有限公司 Time advance prediction method, device and communication system
CN112911502B (en) * 2021-01-20 2022-07-01 东南大学 Base station longitude and latitude estimation method based on MDT technology and RSRP ranging
CN115529605B (en) * 2021-06-24 2024-07-02 中国移动通信集团广西有限公司 Method and device for identifying cell configuration repeater and electronic equipment
CN114374990B (en) * 2021-12-23 2024-03-22 北京红山信息科技研究院有限公司 Method for improving precision between MR and base station based on TA calculation

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002033699A (en) * 2000-07-18 2002-01-31 Pioneer Electronic Corp Mobile communication system
CN101986747A (en) * 2010-10-14 2011-03-16 北京拓明科技有限公司 Mobile terminal positioning method
CN103702415A (en) * 2013-12-16 2014-04-02 深圳市博瑞得科技有限公司 Mobile terminal positioning method and device as well as base station abnormal examination method and device
CN103858462A (en) * 2013-11-22 2014-06-11 华为技术有限公司 Positioning method and base station
CN104602336A (en) * 2014-12-25 2015-05-06 大连楼兰科技股份有限公司 Base station location method based on AT instructions in GSM network
CN108306699A (en) * 2018-02-08 2018-07-20 南京华苏科技有限公司 A kind of antenna feeder optimization method estimated based on gain
CN109302714A (en) * 2018-12-07 2019-02-01 南京华苏科技有限公司 Realize that base station location is studied and judged and area covered knows method for distinguishing based on user data

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106899988A (en) * 2016-06-03 2017-06-27 中国移动通信集团设计院有限公司 A kind of road measuring method and device
CN110351655B (en) * 2019-07-10 2021-06-18 博睿泰克科技（宁波）有限公司 Indoor positioning method and system based on signal multipath propagation measurement

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002033699A (en) * 2000-07-18 2002-01-31 Pioneer Electronic Corp Mobile communication system
CN101986747A (en) * 2010-10-14 2011-03-16 北京拓明科技有限公司 Mobile terminal positioning method
CN103858462A (en) * 2013-11-22 2014-06-11 华为技术有限公司 Positioning method and base station
CN103702415A (en) * 2013-12-16 2014-04-02 深圳市博瑞得科技有限公司 Mobile terminal positioning method and device as well as base station abnormal examination method and device
CN104602336A (en) * 2014-12-25 2015-05-06 大连楼兰科技股份有限公司 Base station location method based on AT instructions in GSM network
CN108306699A (en) * 2018-02-08 2018-07-20 南京华苏科技有限公司 A kind of antenna feeder optimization method estimated based on gain
CN109302714A (en) * 2018-12-07 2019-02-01 南京华苏科技有限公司 Realize that base station location is studied and judged and area covered knows method for distinguishing based on user data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GSM网络的规划和优化;高文波;《科技情报开发与经济》;20061231;全文 *

Also Published As

Publication number Publication date
CN111031474A (en) 2020-04-17

Similar Documents

Publication Publication Date Title
CN111031474B (en) Method for predicting longitude and latitude of base station based on user MDT data
JP5032586B2 (en) Method and configuration for improved cell identification and cell positioning
KR101729873B1 (en) Network coverage planning method and apparatus of evolution communication system
EP2283683B1 (en) Location services based on positioned wireless measurement reports
US6631263B1 (en) Cell hand-off border identification using MS positioning and signal strength values
KR101418663B1 (en) Mobile communication system, base station, and cell coverage control method
JP4981068B2 (en) Positioning by path loss polygon
Han et al. Automatic precision control positioning for wireless sensor network
KR20160121160A (en) Apparatus and method for positioning in wireless communication system
JP2005123662A (en) Position detecting system and apparatus thereof
EP2602634A1 (en) Mobile geolocation
EP3336578A1 (en) Device and method for estimating position of terminal in wireless communication system
CN101646201A (en) Method, device and system for determining terminal position
JP2022517223A (en) Positioning method and positioning device using beam information
JP2021536154A (en) Methods implemented on network equipment and methods implemented on terminal equipment
CN102647725B (en) Method and device for determining neighbor cell priority in neighbor cell planning
EP2865214A1 (en) Method and network node for enabling position determination of a user equipment measurement
CN108226860A (en) Ultra wide band mixing dimension localization method and alignment system based on RSS
CN114466448A (en) Positioning method, positioning device and processor-readable storage medium
CN105491586A (en) Method and system for measuring azimuth angle of cell base station antenna
CN102547836B (en) Method, device and system for determining motion trail of mobile terminal
US20240155537A1 (en) Methods for positioning a wireless device, a related wireless node and a related location network node
CN103686998A (en) Switching position positioning method and device based on AOA
EP2934052A1 (en) Method, base station, system and computer-readable medium to localize a user device
CN107113766A (en) Exchange method, the apparatus and system of positional parameter

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant