CN111031474A - 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
CN111031474A
CN111031474A CN201911263121.4A CN201911263121A CN111031474A CN 111031474 A CN111031474 A CN 111031474A CN 201911263121 A CN201911263121 A CN 201911263121A CN 111031474 A CN111031474 A CN 111031474A
Authority
CN
China
Prior art keywords
latitude
base station
user
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.)
Granted
Application number
CN201911263121.4A
Other languages
Chinese (zh)
Other versions
CN111031474B (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 is | delta x | | survival of the fly ash<E, 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; (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 BDA0002312116610000024
Time advance TaAnd reference signal receptionA 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 less than epsilon, wherein epsilon is a positioning error of convergence, 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 BDA0002312116610000021
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 BDA0002312116610000022
Calculating the mean of all user longitudes
Figure BDA0002312116610000023
Standard deviation of
Figure BDA0002312116610000031
Mean value of all user latitudes
Figure BDA0002312116610000032
And standard deviation of
Figure BDA0002312116610000033
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure BDA0002312116610000034
and keeping the user i, otherwise, eliminating the user, wherein the number of the users left after screening is M.
As a preferred technical scheme of the present invention, the specific process of converting the M user longitudes and latitudes obtained in step (2) into two-dimensional plane coordinates by the mercator projection method in step (3) is as follows: taking the reference latitude as
Figure BDA0002312116610000035
The major radius of the earth is a, the minor radius is b, and the first eccentricity of the earth is
Figure BDA0002312116610000036
The latitude of the reference has a weft circle radius of
Figure BDA0002312116610000037
Equivalent latitude of ith user
Figure BDA0002312116610000038
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 larger than or equal toAt n, taking the average value of M users satisfying the condition 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 BDA0002312116610000039
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 base station as x0＝(x(0)，y(0))＝(x0，y0) Setting k to 0;
s72 substituting the K users obtained in the step (6) to calculate vectors
Figure BDA0002312116610000041
Sum matrix
Figure BDA0002312116610000042
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 BDA0002312116610000043
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).
Is an advantage of the present inventionSelecting a technical scheme, wherein the step (9) specifically comprises the following steps: 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, taking the base station latitude initial point as the predicted value x ═ x, y
Figure BDA0002312116610000044
Setting k to 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure BDA0002312116610000051
Wherein
Figure BDA0002312116610000052
Where e and r0The first eccentricity of the earth and the reference latitude in the step (2) are respectively
Figure BDA0002312116610000053
The weft circle radius of (1);
s103 if
Figure BDA0002312116610000054
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 BDA0002312116610000055
Predicted value of latitude is
Figure BDA0002312116610000056
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 onlyaUnder the condition of reference signal received strength RSRP, the method 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 BDA00023121166100000611
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 used for screening the longitude and latitude of the user according to the 3-Sigma criterionComprises 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 BDA0002312116610000061
Calculating the mean of all user longitudes
Figure BDA0002312116610000062
Standard deviation of
Figure BDA0002312116610000063
Mean value of all user latitudes
Figure BDA0002312116610000064
And standard deviation of
Figure BDA0002312116610000065
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure BDA0002312116610000066
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 BDA0002312116610000067
The major radius of the earth is a, the minor radius is b, and the first eccentricity of the earth is
Figure BDA0002312116610000068
The latitude of the reference has a weft circle radius of
Figure BDA0002312116610000069
Equivalent latitude of ith user
Figure BDA00023121166100000610
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, 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)；
(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) carrying out secondary screening according to the iteration points in the step (5), and obtaining K users meeting the conditions by using the users of the dry iteration: user i (i ═ 1,2,. M,) to base station initiation point B (x)0，y0) A distance of
Figure BDA0002312116610000071
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 base station as x0＝(x(0)，y(0))＝(x0，y0) Setting k to 0;
s72 substituting the K users obtained in the step (6) to calculate vectors
Figure BDA0002312116610000072
Sum matrix
Figure BDA0002312116610000073
Obtaining a displacement vector Δ x ═ GTG)-1GTb；
(8) Calculating a displacement vector delta x according to the matrix G and the vector b, if the length of the displacement vector is less than epsilon, turning to the step (10), and if not, executing the step (9); the method specifically comprises the following steps: if the length of the displacement vector
Figure BDA0002312116610000074
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 as k +1, returning to the step (7), and circulating the step (7) and the step (8) until | | delta x | | < 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, setting k to k +1, returning to step (7), and circulating step (7) and step (8) until | | | Δ x | | < epsilon;
(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 BDA0002312116610000081
S101, taking the base station latitude initial point as the predicted value x ═ x, y
Figure BDA0002312116610000082
Setting k to be 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure BDA0002312116610000083
Wherein
Figure BDA0002312116610000084
Where e and r0The first eccentricity of the earth and the reference latitude in the step (2) are respectively
Figure BDA0002312116610000085
The weft circle radius of (1);
s103 if
Figure BDA0002312116610000086
If the latitude calculation is converged, executing the step S104, otherwise, returning to the step S103;
s104 obtains the longitude forecast value of the base station as
Figure BDA0002312116610000087
Predicted value of latitude is
Figure BDA0002312116610000088
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 BDA0002312116610000089
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 BDA00023121166100000810
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, where the longitude of the 1 st user is 118.640823 °, the latitude is 24.980077 °, and the longitude of the i (i) th user is θ 1,2iLatitude is
Figure BDA00023121166100000811
Calculating the longitude of all usersMean value
Figure BDA00023121166100000812
Standard deviation of
Figure BDA00023121166100000813
Mean value of all user latitudes
Figure BDA0002312116610000091
And standard deviation of
Figure BDA0002312116610000092
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure BDA0002312116610000093
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 BDA0002312116610000094
The long radius a of the earth is 6378137, the short radius b is 6356752, and the first eccentricity of the earth is
Figure BDA0002312116610000095
Radius of latitude of reference
Figure BDA0002312116610000096
The equivalent latitude of the i (i 1, 2.., M31867) th user
Figure BDA0002312116610000097
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 BDA0002312116610000098
The equivalent latitude q of the survey base station can be obtainedi＝qbs0.44767272, calculating the plane coordinates of the survey base station
Figure BDA0002312116610000099
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 BDA00023121166100000910
(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 BDA0002312116610000101
Screening out the satisfying | di-ri|≤d050 user ═ 50
Figure BDA0002312116610000105
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), setting k equal to 0;
substituting K into 6447 users in step (6), and calculating
Figure BDA0002312116610000102
Sum matrix
Figure BDA0002312116610000103
Obtain the displacement vector Δ x ═ (G)TG)-1GTb；
(8) If the length of the displacement vector
Figure BDA0002312116610000104
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 coordinate prediction is x ═ x, y ═ 11986499.0, 2591436.4, and the base station coordinate prediction error is:
Figure BDA0002312116610000111
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 BDA0002312116610000112
S101, taking the base station latitude initial point as the predicted value of the base station plane coordinate of x ═ x, y ═ 11986499.0, 2591436.4
Figure BDA0002312116610000113
Setting k to be 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure BDA0002312116610000114
Wherein
Figure BDA0002312116610000115
S103 if
Figure BDA0002312116610000116
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 BDA0002312116610000117
Predicted value of latitude is
Figure BDA0002312116610000118
Unit is rad, then by
Figure BDA0002312116610000119
And
Figure BDA00023121166100001110
obtaining the longitude and latitude predicted value of the base station with the unit of degree
Figure BDA00023121166100001111
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 (10)

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 including longitude theta and latitude of user
Figure FDA0002312116600000011
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 FDA0002312116600000012
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 FDA0002312116600000013
Where i is 1,2, … …, N, the mean of all user longitudes is calculated
Figure FDA0002312116600000014
Standard deviation of
Figure FDA0002312116600000015
Mean value of all user latitudes
Figure FDA0002312116600000016
And standard deviation of
Figure FDA0002312116600000017
If the longitude and latitude of the ith user satisfy the thetai-μθ|＜3σθAnd
Figure FDA0002312116600000018
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 FDA0002312116600000021
The major radius of the earth is a, the minor radius is b, and the first eccentricity of the earth is
Figure FDA0002312116600000022
The latitude of the reference has a weft circle radius of
Figure FDA0002312116600000023
Equivalent latitude of ith user
Figure FDA0002312116600000024
Then the plane coordinate U of the ith useri(xi,yi) Is xi＝r0θi，yi＝r0qiWherein i ═ 1, 2.., M.
4. The method for predicting the latitude and longitude of the base station based on the MDT data of the user according to claim 3, wherein 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)。
5. The method of claim 4, 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.
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 (6) is specifically as follows:
user i, i ═ 1,2, … …, M to base station initiation point B (x)0,y0) A distance of
Figure FDA0002312116600000025
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.
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 (7) is specifically as follows:
s71 setting the initial point coordinate of base station as x0＝(x(0),y(0))＝(x0,y0) Setting k to 0;
s72 substituting the K users obtained in the step (6) to calculate vectors
Figure FDA0002312116600000031
Sum matrix
Figure FDA0002312116600000032
Obtaining a displacement vector Δ x ═ GTG)-1GTb。
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 (8) is specifically as follows: if the length of the displacement vector
Figure FDA0002312116600000033
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).
9. The method for predicting the latitude and longitude of the base station based on the MDT data of the user according to claim 8, 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).
10. The method of claim 9, wherein the step (10) comprises the steps of:
s101, taking the initial point of the latitude of the base station, wherein the estimated value of the plane coordinate of the base station is x (x, y)
Figure FDA0002312116600000034
Setting k to 0;
s102, calculating a base station latitude predicted value obtained after the kth iteration
Figure FDA0002312116600000035
Wherein
Figure FDA0002312116600000036
Where e and r0The first eccentricity of the earth and the reference latitude in the step (2) are respectively
Figure FDA0002312116600000037
The weft circle radius of (1);
s103 if
Figure FDA0002312116600000041
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 FDA0002312116600000042
Predicted value of latitude is
Figure FDA0002312116600000043
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 true CN111031474A (en) 2020-04-17
CN111031474B 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)

Cited By (4)

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

Citations (9)

* 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
CN106899988A (en) * 2016-06-03 2017-06-27 中国移动通信集团设计院有限公司 A kind of road measuring method and device
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
CN110351655A (en) * 2019-07-10 2019-10-18 博睿泰克科技（宁波）有限公司 A kind of indoor orientation method and system based on the measurement of signal multipath transmisstion

Patent Citations (9)

* 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
CN106899988A (en) * 2016-06-03 2017-06-27 中国移动通信集团设计院有限公司 A kind of road measuring method and device
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
CN110351655A (en) * 2019-07-10 2019-10-18 博睿泰克科技（宁波）有限公司 A kind of indoor orientation method and system based on the measurement of signal multipath transmisstion

Non-Patent Citations (1)

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

Cited By (6)

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

Also Published As

Publication number Publication date
CN111031474B (en) 2021-05-07

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
JP4100320B2 (en) Position detection system and apparatus
KR101729873B1 (en) Network coverage planning method and apparatus of evolution communication system
JP4981068B2 (en) Positioning by path loss polygon
CN102869038B (en) Measurement compensation method and device in positioning of base station, server and system
KR102139792B1 (en) Apparatus and method for adjusting position of RF infrastructure
KR20160121160A (en) Apparatus and method for positioning in wireless communication system
US9807560B2 (en) Method and network node for enabling position determination of a user equipment measurement
CN109141427B (en) EKF positioning method based on distance and angle probability model under non-line-of-sight environment
CN108226860B (en) RSS (received signal strength) -based ultra-wideband mixed dimension positioning method and positioning system
JPWO2011114372A1 (en) Mobile communication system, base station, cell coverage control method
WO2009129344A1 (en) Location services based on positioned wireless measurement reports
CN114466448B (en) Positioning method, positioning device and processor-readable storage medium
JP2022517223A (en) Positioning method and positioning device using beam information
JP2021536154A (en) Methods implemented on network equipment and methods implemented on terminal equipment
CN108737959A (en) 5G super-intensive method for switching network based on Static Precise Point Positioning
WO2024087612A1 (en) Method and apparatus for positioning
CN102547836B (en) Method, device and system for determining motion trail of mobile terminal
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
CN105451263A (en) WCDMA network cell soft handover band identification method, equipment and system
CN110068794A (en) A kind of AOA positioning and optimizing compensation method
CN107113766A (en) Exchange method, the apparatus and system of positional parameter
CN105828342B (en) Method and device for confirming neighbor relation

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