CN108495271B - Terminal positioning method suitable for high-speed railway network test - Google Patents
Terminal positioning method suitable for high-speed railway network test Download PDFInfo
- Publication number
- CN108495271B CN108495271B CN201810194087.9A CN201810194087A CN108495271B CN 108495271 B CN108495271 B CN 108495271B CN 201810194087 A CN201810194087 A CN 201810194087A CN 108495271 B CN108495271 B CN 108495271B
- Authority
- CN
- China
- Prior art keywords
- base station
- coordinates
- points
- track
- gps
- 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
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
- H04W4/029—Location-based management or tracking services
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/02—Arrangements for optimising operational condition
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/30—Services specially adapted for particular environments, situations or purposes
- H04W4/40—Services specially adapted for particular environments, situations or purposes for vehicles, e.g. vehicle-to-pedestrians [V2P]
- H04W4/42—Services specially adapted for particular environments, situations or purposes for vehicles, e.g. vehicle-to-pedestrians [V2P] for mass transport vehicles, e.g. buses, trains or aircraft
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/006—Locating users or terminals or network equipment for network management purposes, e.g. mobility management with additional information processing, e.g. for direction or speed determination
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to the technical field of communication, in particular to a terminal positioning method suitable for high-speed railway network testing, which is used for positioning a test terminal when a traffic trunk network test GPS signal is lost in the communication industry. The method comprises the steps of utilizing an eNodeB azimuth angle of an operator and a high-speed rail map, positioning a test terminal when GPS is lost due to train penetration loss in a network test, discretizing track longitude and latitude points in an eNodeB coverage range, taking the ratio of the track longitude and latitude points to RSR sampling points of the test terminal as an interpolation period, and interpolating the track longitude and latitude points when the GPS is lost as position information of the test terminal. The method can position the position of the test terminal under the condition of different GPS loss rates, namely when the GPS is lost in the whole course or partial road section, the configuration of the existing positioning equipment is not required to be changed, the running cost of the equipment is not increased, and meanwhile, the vehicle positioning efficiency is improved.
Description
Technical Field
The invention relates to the technical field of communication, in particular to a terminal positioning method suitable for high-speed railway network testing, which is used for positioning a test terminal when a traffic trunk network test GPS signal is lost in the communication industry.
Background
The network test is a method for acquiring optimized data of a wireless network, is used for acquiring downlink signal data of a wireless network of a telecommunication operator, and is an essential step in the optimization of the wireless network. Network testing uses test equipment in combination with geographic information to make comprehensive measurements of air interfaces to verify parameters of the radio interfaces, analyze, optimize and locate faults. By accessing a wireless network, the wireless parameters and the service quality of the network are sensed in real time by using the services of the telecom operator, such as voice, data, internet and the like, and the whole testing process is stored. The information of the network test record can be output to processing software for analyzing the network quality condition.
In the high-speed rail network test, a traditional test tool (such as the Pilot Pioneer road test software) is used, and the test tool comprises a GPS positioning module, test software, a terminal and a computer. And the test software records the collected data of the test terminal and the GPS at intervals of a sampling period. Generally, the test terminal reports Reference Signal Received Power (RSRP) information for 1 time at an interval of about 1s, and the GPS reports longitude and latitude information for 1 time. And the RSRP value of each position point can be obtained by the correlation analysis of the test terminal and the GPS collected data. The high-speed rail closed aluminum alloy carriage body has strong penetration loss to electromagnetic waves, so that GPS signals are seriously attenuated. For example, the CRH train carriage penetration loss is about 17dB, and the wind-sand resistant and high-cold resistant CRH5G train carriage penetration loss in the northwest region is as high as 24 dB. Because the high-speed rail network test is influenced by factors such as train carriage penetration loss, tunnel shielding, mountain blocking, U-shaped groove multipath fading and the like, the GPS is difficult to search satellites and the position information is easy to seriously lose. In addition, the influence of the atmospheric waveguide makes weak GPS signals in the high-speed railway carriage irregularly available. Under the condition of GPS loss, the coverage of each eNodeB cannot be determined, so that the eNodeB address planning and network structure adjustment lose the support of position information. The lack of the position information leads to the failure of evaluating the communication quality of each test point, and great difficulty is caused to the optimization and parameter adjustment of the adjacent cells.
In recent years, some test terminal positioning methods with additional means have been proposed, including gyroscope and accelerometer positioning in general. The inertial navigation device is combined with the GPS to realize the positioning of the test terminal. However, third-party drive test software and equipment purchased by telecom operators in a centralized manner only meet basic requirements of outdoor road test scenes, and the problem of GPS loss of high-speed rail test scenes cannot be solved. With the rapid increase of the operating mileage of the high-speed rail year by year, the requirement of acquiring the position information of the test terminal of the special scene of the high-speed rail in real time is met on the basis of the existing outdoor test software and test equipment, and the problem to be solved by telecommunication operators, test equipment and test software suppliers is urgent.
Disclosure of Invention
Aiming at the technical problem, the invention provides a terminal positioning method suitable for high-speed railway network testing, which utilizes communication eNodeB engineering parameter cell azimuth and high-speed rail map discrete longitude and latitude points to perform test terminal position information compensation positioning. In order to solve the problem of GPS loss in the high-speed rail network test, the invention provides a positioning method of a test terminal on the basis of the test equipment intensively purchased by the existing telecom operator without adding any test instrument, thereby expanding the test range of a special scene and simultaneously reducing the additional investment of the telecom operator and a communication equipment supplier.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a terminal positioning method suitable for high-speed railway network test comprises the following steps
S1, substituting the track points into the azimuth angle of the base station to list the ray equation of the main lobe direction of the antenna
The communication base station in the circuit area is marked as eNodeB1，eNodeB2，…eNodeBN-1The vehicle route is provided with a plurality of track measuring points marked as A1，A2，…An，A1，A2，…ANRespectively record the longitude and latitude coordinates of the horizontal plane as (x)1，y1)、(x2，y2)，…(xN，yN) The included angle between the connecting line of any base station eNodeB and the track measuring point and the positive north direction is α, the horizontal plane coordinates of the actual measuring points of the base station eNodeB are (x, y), the horizontal plane coordinates of the actual measuring points are substituted into the antenna main lobe direction ray equation of the actual measuring point azimuth angle to be f1(x)＝tan(90°-a)·(x-x1)+y1；
S2 fitting the high-speed rail curve
Obtaining an actually measured track line curve through an open map, substituting longitude and latitude coordinate parameters of the track line curve into a curve fitting equation f (x) p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7(ii) a Calculating the parameter value of 6-order polynomial fitting as p by using the Curve fit tool of the traditional MATLAB1＝0.001312；p2＝-0.3225；p3＝-34.46；p4＝1.979e+04；p5＝-2.612e+06；p6＝1.492e+08；p7-3.224e + 09; e is the name of a natural constant, about 2.71828, p1To p7Fitting coefficients that are 6 th order polynomials;
s3, solving intersection point coordinates of antenna main lobe direction rays and orbit curve equation
The base station azimuth angle lists the coordinates of the intersection point of the antenna main lobe direction ray equation and the curve fitting equation solution,
when only 1 solution exists, the solution is the coordinate of the intersection point of the main lobe direction ray of the antenna and the high-speed rail track fitting curve and is marked as (x'1，y′1) Expressed as:
(x′1，y′1)＝solve(f1(x)，f(x))；
when there are N solutions, the solution set is { (x'1j，y′1j)}，j＝1,2,…N，
N solution substitutions
Calculating the coordinate of the intersection point closest to the base station;
s4, calculating the switching point of the base station according to the intersection point coordinates
Coordinate (x ') of switching point'Bn，y′Bn) Substituting eNodeBn into the coordinates of the neighboring intersection
And calculating the coordinates of a base station switching point, wherein (x'B1，y′B1) Is the first switching point coordinate, (x'1n，y′1n) And (x'2n，y′2n) Coordinates of two adjacent base stations of the first switching point are respectively;
s5, inserting intersection point coordinates according to period
The GPS sampling frequency of the terminal between the two base station switching points is M, M intersection point coordinates of the step S4 are inserted between the two base station switching points, and the interpolation period is
M is more than or equal to M, through insertingAnd the coordinates of the intersection points are used for correcting the position information of the high-speed rail, so that the GPS signals are prevented from being lost.
The distance between the two base station switching points in the step S4 is less than 1.6 km.
The GPS sampling frequency of the terminal between the two base station switching points in step S5 is M, M intersection coordinates of step S4 are inserted between the two base station switching points, and the interpolation period is M
The invention has the beneficial effects that:
the invention relates to a terminal positioning method suitable for high-speed railway network testing, which utilizes an azimuth angle of an operator eNodeB (a base station in a fourth generation mobile communication network LTE) and a high-speed rail map to position a test terminal when a GPS is lost due to train penetration loss in network testing. The method comprises the following steps: calculating the main lobe direction ray of the antenna according to the azimuth angle of the eNodeB; fitting a curve equation of the high-speed rail by taking the high-speed rail map as a reference; solving the intersection point of the ray and the curve, namely solving the ray equation and the curve equation; solving the switching point of the adjacent eNodeB according to the intersection point; calculating the coverage range of the eNodeB on the high-speed rail, namely the track between the two switching points; discretizing track longitude and latitude points in the coverage range of eNodeB, taking the ratio of the track longitude and latitude points to the RSRP (Reference Signal Receiving Power) sampling points of the test terminal as an interpolation period, and interpolating the track longitude and latitude points when the GPS is absent, and taking the interpolated track longitude and latitude points as the position information of the test terminal. The method can position the position of the test terminal under the condition of different GPS loss rates, namely when the GPS is lost in the whole course or partial road section, the configuration of the existing positioning equipment is not required to be changed, the running cost of the equipment is not increased, and meanwhile, the vehicle positioning efficiency is improved.
Drawings
FIG. 1 is a model of a terminal location method of the present invention;
fig. 2 is a flow chart of a terminal positioning method of the present invention;
FIG. 3 shows the contrast effect of embodiment 1 of the present invention;
FIG. 4 shows the comparative effect of example 2 of the present invention;
FIG. 5 shows the comparison effect of embodiment 3 of the present invention;
FIG. 6 is a first order comparison of the high speed rail curve fitting of the present invention;
FIG. 7 is a graph of a second order comparison of the high speed rail curve fit of the present invention;
FIG. 8 is a third order comparison effect of curve fitting for high speed rail in accordance with the present invention;
FIG. 9 is a fourth order comparison of the high speed rail curve fit of the present invention;
FIG. 10 is a graph of the fifth order contrast effect of the high rail curve fitting of the present invention;
FIG. 11 is a graph of the sixth order contrast effect of the high rail curve fitting of the present invention.
Detailed description of the invention
The invention and its effects will be further elucidated with reference to the drawings and the embodiments.
A terminal positioning method suitable for high-speed railway network test comprises the following steps
S1, substituting the track points into the azimuth angle of the base station to list the ray equation of the main lobe direction of the antenna
The communication base station in the circuit area is marked as eNodeB1，eNodeB2，…eNodeBN-1The vehicle route is provided with a plurality of track measuring points marked as A1，A2，…An，A1，A2，…ANRespectively record the longitude and latitude coordinates of the horizontal plane as (x)1，y1)、(x2，y2) And (c), an included angle between a connecting line of any base station eNodeB and a track measuring point and a positive north direction is α, the horizontal plane coordinates of actual measuring points of the base station eNodeB are (x, y) respectively, and the real measuring point horizontal plane coordinates are substituted into a main lobe direction ray equation f of the antenna at the actual measuring point azimuth angle1(x)＝tan(90°-a)·(x-x1)+y1；
S2 fitting the high-speed rail curve
Obtaining an actually measured track line curve through an open map, substituting longitude and latitude coordinate parameters of the track line curve into a curve fitting equation f (x) p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7(ii) a Calculating the parameter value of 6-order polynomial fitting as p by using the Curve fit tool of the traditional MATLAB1＝0.001312；p2＝-0.3225；p3＝-34.46；p4＝1.979e+04；p5＝-2.612e+06；p6＝1.492e+08；p7-3.224e + 09; e is the name of a natural constant, about 2.71828, p1To p7Fitting coefficients that are 6 th order polynomials;
s3, solving intersection point coordinates of antenna main lobe direction rays and orbit curve equation
The base station azimuth angle lists the coordinates of the intersection point of the antenna main lobe direction ray equation and the curve fitting equation solution,
when only 1 solution exists, the solution is the coordinate of the intersection point of the main lobe direction ray of the antenna and the high-speed rail track fitting curve and is marked as (x'1，y′1) Is represented by (x'1，y′1)＝solve(f1(x)，f(x))；
When there are N solutions, the solution set is { (x'1j，y′1j)}，j＝1,2,…N，
N solution substitutions
Calculating the coordinate of the intersection point closest to the base station;
s4, calculating the switching point of the base station according to the intersection point coordinates
Coordinate (x ') of switching point'Bn，y′Bn) Substituting eNodeBn into the coordinates of the neighboring intersection
And calculating the coordinates of a base station switching point, wherein (x'B1，y′B1) Is the first switching point coordinate, (x'1n，y′1n) And (x'2n，y′2n) Coordinates of two adjacent base stations of the first switching point are respectively;
s5, inserting intersection point coordinates according to period
The GPS sampling frequency of the terminal between the two base station switching points is M, M intersection point coordinates of the step S4 are inserted between the two base station switching points, and the interpolation period is
M is larger than or equal to M, and the position information of the high-speed rail is corrected by inserting intersection point coordinates to prevent GPS signal loss.
The distance between the two base station switching points in the step S4 is less than 1.6 km.
The GPS sampling frequency of the terminal between the two base station switching points in step S5 is M, M intersection coordinates of step S4 are inserted between the two base station switching points, and the interpolation period is M
M＞m。
In fig. 1, the secondary circuit includes two base stations eNodeB corresponding to two pairs of antennas, respectively, and the included angles between the north-positive direction clockwise and the main lobe directions of the two pairs of antennas are β 1 and β 2, respectively, that is, the azimuth angle of the eNodeB2-1 cell is β 1, and the azimuth angle of the eNodeB2-2 cell is β 2; the switching point between eNodeBs indicates that the service cells occupied by the test terminal before and after the position are switched; a1, A2, A3 and A4 are the intersection points of the ray in the main lobe direction of the antenna and the track curve; b1 is a switching point of eNodeB1 and eNodeB2, the test terminal occupies eNodeB1 wireless resources on the left track of B1 and occupies eNodeB2 wireless resources on the right track of B1; b2 is the handover point between eNodeB2 and eNodeB3, and the test terminal occupies eNodeB2 radio resources on the left track of B2 and occupies eNodeB3 radio resources on the right track of B2. It can be seen that the eNodeB2 covers a trajectory between points B1, B2, from which the coverage of each eNodeB along the high-speed rail can be determined. If the GPS signal is missing when the test terminal occupies eNodeB2, the missing position information can be compensated by the track longitude and latitude between points B1 and B2.
The flow of the positioning method based on the cell azimuth angle in the high-speed rail VoLTE test is shown in fig. 2, and the steps are as follows:
1) listing an antenna main lobe direction ray equation according to the azimuth angle;
2) fitting a curve equation of the high-speed rail;
3) solving the coordinates of the intersection point of the main lobe direction ray of the antenna and the orbit curve equation;
4) if the number of the intersection points is larger than 1, the distance between each intersection point and the eNodeB is calculated, and the intersection point with the minimum distance is selected; if the number of the intersection points is equal to 1, executing the next step;
5) solving a switching point of an adjacent eNodeB according to the intersection point coordinates;
6) calculating the coverage range of the eNodeB on the high-speed rail, namely the track between the two switching points;
7) discretizing track longitude and latitude points in the coverage range of the eNodeB, taking the ratio of the track longitude and latitude points to the RSRP sampling points of the test terminal as an interpolation period, and interpolating the track longitude and latitude points when the GPS is absent, and taking the interpolated track longitude and latitude points as the position information of the test terminal.
S1: calculating curve equation of high-speed rail
And discretizing the high-speed rail map into longitude and latitude points, wherein the distance between the longitude and latitude points is 5 m. The railway track curve equation can be obtained by a fitting method. Using the Curve fixing Tool of MATLAB to perform 6-order polynomial Fitting on the track map of Zhangye West station-Jiuquan south station to obtain a track Curve equation of
f(x)＝p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
Calculating a parameter value of 6-order polynomial Fitting by using the Curve Fitting Tool of MATLAB
p1＝0.001312；p2＝-0.3225；p3＝-34.46；p4＝1.979e+04；p5＝-2.612e+06；p6＝1.492e+08；p7＝-3.224e+09。
e is the name of a natural constant, about 2.71828. p is a radical of1To p7Fitting coefficients of 6 th order polynomial.
S2: calculating the equation of the ray of the main lobe direction of the cell antenna
Let the azimuth angles of eNodeB1-1, 2-1, 2-2 and 3-1 cells be α1，β2And gamma. Suppose the two-dimensional position coordinates of eNodeB1, eNodeB2, and eNodeB3 are (x) respectively1,y1)、(x2,y2) And (x)3,y3). Then the equations for the antenna main lobe direction rays of eNodeB1-1 cell, eNodeB2-1 cell, eNodeB2-2 cell and eNodeB3-1 cell are found in equations (3) -6)
f1(x)＝tan(90°-a)·(x-x1)+y1(3)
f2(x)＝tan(90°-β1)·(x-x2)+y2(4)
f3(x)＝tan(90°-β2)·(x-x2)+y2(5)
f4(x)＝tan(90°-γ)·(x-x3)+y3(6)
When there are 1 and only 1 solutions, the solution is the coordinate of the intersection point of the antenna main lobe direction ray and the high-speed rail track fitting curve, and is marked as (x'1，y′1)，(x′2，y′2)，(x′3，y′3) And (x'4，y′4) Can be represented as
Consider the eNodeB1-1 cell antenna main lobe directional ray equation f1(x) Fitting to gaussian the case of multiple intersections of the orbital curve equation f (x). Deconstructed can be written as { (x'1j，y′1j) J ═ 1,2, … N. Suppose { (x'1n，y′1n) Is closest to eNodeB1 (1 ≦ N ≦ N), which may be expressed as
Then the nth solution is retained, (x'1n，y′1n) I.e., the coordinates of intersection a 1. Similarly, the coordinates of the intersection points a2, A3, and a4 can be found.
Ideally, the switching point is equidistant from two adjacent intersections. That is, in fig. 2, the distances from the intersection points a1 and a2 are equal for the switching point B1 of the eNodeB1 and the eNodeB2, and the distances from the intersection points A3 and a4 are equal for the switching point B2 of the eNodeB2 and the eNodeB 3. Thus, according to formula (7)The coordinates (x 'of the switching points B1 and B2 can be obtained'B1，y′B1) And (x'B2，y′B2)
The track between the switch points B1, B2 is the coverage of the eNodeB2 on the high-speed rail. Similarly, the coverage of all enodebs along the high-speed rail can be obtained.
The GPS sampling frequency of the terminal between two base station switching points is M, M intersection point coordinates are inserted between the two base station switching points, the test terminal lacks GPS position information in the coverage range of eNodeB2, the M is the number of RSRP sampling points, the M is the number of discrete longitude and latitude points between the switching points B1 and B2 on the track, and the interpolation period is M
Therefore, the discrete longitude and latitude points between the switching points B1 and B2 extracted in the interpolation period T can be used as the position information of the test terminal in the coverage of the eNodeB2, so that the test terminal can be positioned when the GPS is absent at any sampling time.
Example 1
The present invention locates the location information when the test terminal occupies eNodeB 1. The coordinates of the eNodeB1 are (99.346405, 39.323753), the azimuth of the eNodeB1-1 cell is 105 °, and the azimuth of the eNodeB1-2 cell is 265 °. The equations of the main lobe direction rays of the two cell antennas are respectively as follows:
f2(x) Tan (90 ° -105 °) · (x-99.346405) +39.323753 and f3(x) Tan (90 ° -265 °) · (x-99.346405) + 39.323753. The neighboring cell coordinates of the eNodeB1-1 cell southeast direction neighboring eNodeB are (99.361522,39.322955), the azimuth angle is 260 °, and the antenna main lobe direction ray equation is: f. of1(x) Tan (90 ° -260 °) · (x-99.361522) + 39.322955. The coordinate of the neighboring eNodeB in the northwest direction of the eNodeB1-2 cell is (99.333510,39.325803), the azimuth angle is 120 °, and the ray equation of the antenna main lobe direction is: f. of4(x)＝tan(90°-120°)·(x-99.333510)+39.325803。
The 6 th order polynomial fitting curve equation of the Lanxin high-iron Zhangye West station-the Jiuquan south station is as follows:
f(x)＝p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
calculating the parameter value of 6-order polynomial Fitting as p 1-0.001312 by using the Curve Fitting Tool of MATLAB; p2 ═ -0.3225; p3 ═ -34.46; p4 ═ 1.979e + 04; p5 ═ -2.612e + 06; p6 ═ 1.492e + 08; p7 ═ 3.224e + 09.
e is the name of a natural constant, about 2.71828. P1 to P7 are fitting coefficients of 6 th order polynomials.
f1(x)、f2(x)、f3(x) And f4(x) The intersection points with f (x) are a1(99.357261, 39.322252), a2(99.348314, 39.323287), A3(99.345036, 39.323681) and a4(99.334996, 39.324936), respectively. Points on the track equidistant from a1(99.357261, 39.322252) and a2(99.348314, 39.323287) are B1(99.352703, 39.322498), and points on the track equidistant from A3(99.345036, 39.323681) and a4(99.334996, 39.324936) are B2(99.340070, 39.324550). I.e., the track between B1(99.352703, 39.322498) and B2(99.340070, 39.324550), is the coverage of eNodeB1 on high-speed rail. The test terminal occupies 17 RSRP sampling points of eNodeB2, 17 longitude and latitude points are uniformly taken from B1(99.352703, 39.322498) and B2(99.340070, 39.324550) on a high-speed rail, namely, when the test terminal occupies eNodeB1, the longitude and latitude list is as follows:
example 2
The present invention locates the location information when the test terminal occupies eNodeB 2. The coordinates of eNodeB2 are (99.169312, 39.359299), the azimuth of eNodeB2-1 cell is 90 °, and the azimuth of eNodeB2-2 cell is 290 °. Two smallThe equations of the main lobe direction rays of the regional antenna are respectively as follows: f. of2(x) Tan (90 ° -90 °) · (x-99.169312) +39.359299 and f3(x) Tan (90 ° -295 °) · (x-99.169312) + 39.359299. The neighboring cell coordinates of eNodeB2-1 cell southeast direction neighboring eNodeB are (99.184611,39.359588), the azimuth angle is 265 °, the antenna main lobe direction ray equation is: f. of1(x) Tan (90 ° -265 °) · (x-99.184611) + 39.359588. The coordinate of the neighboring eNodeB in the northwest direction of the eNodeB2-2 cell is (99.155389,39.361827), the azimuth angle is 90 °, and the ray equation of the antenna main lobe direction is: f. of4(x)＝tan(90°-90°)·(x-99.155389)+39.361827。
The 6 th order Gaussian fitting curve equation of Lanxin Ganzhang Xita-Jiuquan Nanta is as follows:
f(x)＝p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
calculating a parameter value of 6-order polynomial Fitting by using the Curve Fitting Tool of MATLAB
p1＝0.001312；p2＝-0.3225；p3＝-34.46；p4＝1.979e+04；p5＝-2.612e+06；p6＝1.492e+08；p7＝-3.224e+09。
e is the name of a natural constant, about 2.71828. p is a radical of1To p7Fitting coefficients of 6 th order polynomial.
f1(x)、f2(x)、f3(x) And f4(x) The intersection points with f (x) are a1(99.172832, 39.35857), a2(99.170402, 39.359301), A3(99.167242, 39.360259) and a4(99.162254, 39.361802), respectively. Points on the track equidistant from a1(99.172832, 39.35857) and a2(99.170402, 39.359301) are B1(99.171959, 39.360699), and points on the track equidistant from A3(99.167242, 39.360259) and a4(99.162254, 39.361802) are B2(99.164958, 39.362086). I.e., the track between B1(99.171959, 39.360699) and B2(99.164958, 39.362086), is the coverage of eNodeB2 on high-speed rail. The test terminal occupies eNodeB2 with 9 RSRP sampling points, B1(99.171959, 39) on the high-speed rail track360699) and B2(99.164958, 39.362086) take 9 longitude and latitude points uniformly, that is, when the positioning method of the invention is used, the list of the test terminal occupying eNodeB2 longitude and latitude is as follows:
example 3
The present invention locates the location information when the test terminal occupies eNodeB 3. The coordinates of eNodeB3 are (98.976210, 39.443899), the azimuth of 1 cell is 160 °, and the azimuth of 2 cell is 290 °. The equations of the main lobe direction rays of the two cell antennas are respectively as follows: f. of2(x) Tan (90 ° -160 °) · (x-98.976210) +39.443899 and f3(x) Tan (90 ° -290 °) · (x-98.976210) + 39.443899. The neighboring cell coordinates of the eNodeB3-1 cell southeast direction neighboring eNodeB are (98.982611,39.39.435891), the azimuth angle is 300 °, and the antenna main lobe direction ray equation is: f. of1(x) Tan (90 ° -300 °) · (x-98.982611) + 39.435891. The coordinate of the neighboring eNodeB in the northwest direction of the eNodeB3-2 cell is (98.968178,39.449521), the azimuth angle is 155 °, and the ray equation of the antenna main lobe direction is: f. of4(x)＝tan(90°-155°)·(x-98.968178)+39.449521。
The 6 th order Gaussian fitting curve equation of Lanxin Ganzhang Xita-Jiuquan Nanta is as follows:
f(x)＝p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
calculating a parameter value of 6-order polynomial Fitting by using the Curve Fitting Tool of MATLAB
p1＝0.001312；p2＝-0.3225；p3＝-34.46；p4＝1.979e+04；p5＝-2.612e+06；p6＝1.492e+08；p7＝-3.224e+09。
e is the name of a natural constantAnd is about 2.71828. p is a radical of1To p7Fitting coefficients of 6 th order polynomial.
f1(x)、f2(x)、f3(x) And f4(x) The intersection points with f (x) are a1(98.982611, 39.435891), a2(98.976975, 39.44177), A3(98.967195, 39.447177) and a4(98.970015, 39.445607), respectively. Points on the track equidistant from a1(98.982611, 39.435891) and a2(98.976975, 39.44177) are B1(98.981471, 39.440866), and points on the track equidistant from A3(98.967195, 39.447177) and a4(98.970015, 39.445607) are B2(98.967939, 39.451956). I.e., the track between B1(98.981471, 39.440866) and B2(98.967939, 39.451956), is the coverage of eNodeB3 on high-speed rail. The test terminal occupies eNodeB2 and has 9 RSRP sampling points, and 32 longitude and latitude points are uniformly taken from B1(98.981471, 39.440866) and B2(98.967939, 39.451956) on a high-speed rail, namely when the positioning method is used, the longitude and latitude list when the test terminal occupies eNodeB3 is as follows:
comparison of practical case 1 with conventional technique
In case 1, the terminal occupies a test distance of eNodeB1 signal of about 935 m, and manually collects 1 piece of GPS information along about 55 m of the high-speed rail track, and collects longitude and latitude information as follows, and the comparison effect is as shown in fig. 3.
| Serial number | GPS longitude | GPS latitude |
| 1 | 99.352261 | 39.322591 |
| 2 | 99.351404 | 39.322537 |
| 3 | 99.350716 | 39.322497 |
| 4 | 99.350146 | 39.322630 |
| 5 | 99.349346 | 39.322669 |
| 6 | 99.348698 | 39.322795 |
| 7 | 99.347902 | 39.322854 |
| 8 | 99.347287 | 39.323102 |
| 9 | 99.346521 | 39.323259 |
| 10 | 99.345738 | 39.323299 |
| 11 | 99.345149 | 39.323594 |
| 12 | 99.344452 | 39.323732 |
| 13 | 99.343663 | 39.323784 |
| 14 | 99.342969 | 39.323974 |
| 15 | 99.342186 | 39.324067 |
| 16 | 99.340836 | 39.324284 |
| 17 | 99.340325 | 39.324578 |
Practical case 2 is compared with the conventional technology
In case 2, the test distance of the terminal occupying eNodeB1 signal is about 495 m, 1 piece of GPS information is manually acquired along about 55 m of the high-speed rail track, and the longitude and latitude information is acquired as follows, and the comparison effect is as shown in fig. 4.
| Serial number | GPS longitude | GPS latitude |
| 1 | 99.171768 | 39.360816 |
| 2 | 99.170904 | 39.360884 |
| 3 | 99.170294 | 39.361074 |
| 4 | 99.169408 | 39.361191 |
| 5 | 99.168753 | 39.361386 |
| 6 | 99.168018 | 39.361581 |
| 7 | 99.167169 | 39.361599 |
| 8 | 99.166451 | 39.361862 |
| 9 | 99.165712 | 39.361964 |
Practical case 3 is compared with the conventional technology
In case 3, the test distance of the terminal occupying eNodeB1 signal is about 1760 m, 1 piece of GPS information is manually acquired along about 55 m of the high-speed rail track, and the longitude and latitude information is acquired as follows, and the comparison effect is as shown in fig. 5.
Claims (3)
1. A terminal positioning method suitable for high-speed railway network test is characterized by comprising the following steps
S1, substituting the track points into the azimuth angle of the base station to list the ray equation of the main lobe direction of the antenna
The communication base station in the circuit area is marked as eNodeB1，eNodeB2，…eNodeBN-1The vehicle route is provided with a plurality of track measuring points marked as A1，A2，…AN，A1，A2，…ANRespectively record the longitude and latitude coordinates of the horizontal plane as (x)1，y1)、(x2，y2)，…(xN，yN) The included angle between the connecting line of any base station eNodeB and the track measuring point and the positive north direction is α, the horizontal plane coordinates of the actual measuring points of the base station eNodeB are (x, y), the horizontal plane coordinates of the actual measuring points are substituted into the antenna main lobe direction ray equation of the actual measuring point azimuth angle to be f1(x)＝tan(90°-a)·(x-x1)+y1；
S2 fitting the high-speed rail curve
Obtaining an actually measured track line curve through an open map, substituting longitude and latitude coordinate parameters of the track line curve into a curve fitting equation f (x) p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7(ii) a Calculating the parameter value of 6-order polynomial Fitting as p by using the Curve Fitting Tool of the traditional MATLAB1＝0.001312；p2＝-0.3225；p3＝-34.46；p4＝1.979e+04；p5＝-2.612e+06；p6＝1.492e+08；p7-3.224e + 09; e is the name of a natural constant, about 2.71828, p1To p7Fitting coefficients that are 6 th order polynomials;
s3, solving intersection point coordinates of antenna main lobe direction rays and orbit curve equation
The base station azimuth angle lists the coordinates of the intersection point of the antenna main lobe direction ray equation and the curve fitting equation solution,
when only 1 solution exists, the solution is the coordinate of the intersection point of the main lobe direction ray of the antenna and the high-speed rail track fitting curve and is marked as (x'1，y′1) Is shown as
(x′1，y′1)＝solve(f1(x)，f(x))；
When there are N solutions, the solution set is { (x'1j，y′1j)}，j＝1,2,…N，
N solution substitutions
Calculating the coordinate of the intersection point closest to the base station;
s4, calculating the switching point of the base station according to the intersection point coordinates
Coordinate (x ') of switching point'Bn，y′Bn) Substituting eNodeBn into the coordinates of the neighboring intersection
And calculating the coordinates of a base station switching point, wherein (x'B1，y′B1) Is the first switching point coordinate, (x'1n，y′1n) And (x'2n，y′2n) Coordinates of two adjacent base stations of the first switching point are respectively;
s5, inserting intersection point coordinates according to period
The GPS sampling frequency of the terminal between the two base station switching points is M, M intersection point coordinates of the step S4 are inserted between the two base station switching points, and the interpolation period is
M is larger than or equal to M, and the position information of the high-speed rail is corrected by inserting intersection point coordinates to prevent GPS signal loss.
2. The method as claimed in claim 1, wherein the distance between the two base station switching points in step S4 is less than 1.6 km.
3. The method as claimed in claim 1, wherein the GPS sampling time of the terminal between the two base station switch points in step S5 is M, M intersection coordinates of step S4 are inserted between the two base station switch points, and the interpolation period is M
M＞m。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810194087.9A CN108495271B (en) | 2018-03-09 | 2018-03-09 | Terminal positioning method suitable for high-speed railway network test |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810194087.9A CN108495271B (en) | 2018-03-09 | 2018-03-09 | Terminal positioning method suitable for high-speed railway network test |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108495271A CN108495271A (en) | 2018-09-04 |
| CN108495271B true CN108495271B (en) | 2020-09-25 |
Family
ID=63338464
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810194087.9A Active CN108495271B (en) | 2018-03-09 | 2018-03-09 | Terminal positioning method suitable for high-speed railway network test |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108495271B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111288942B (en) * | 2018-12-07 | 2021-09-03 | 比亚迪股份有限公司 | Track transponder position measuring method and device and computer equipment |
| CN111340908B (en) * | 2020-02-10 | 2023-07-25 | 兰州交通大学 | Rail electronic map generation method |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103523056A (en) * | 2013-10-14 | 2014-01-22 | 兰州交通大学 | Train positioning method based on railway mobile communication network communication base stations |
| WO2015165682A1 (en) * | 2014-05-02 | 2015-11-05 | Mclaren Applied Technologies Limited | Location based connections |
| CN105657819A (en) * | 2014-11-12 | 2016-06-08 | 北京神州泰岳软件股份有限公司 | Railway mobile network test longitude and latitude obtaining method and device |
| CN105871448A (en) * | 2016-03-30 | 2016-08-17 | 中科凯普（天津）卫星导航通信技术有限公司 | High-dynamic transformation method of broadband special network of high-speed rail |
| CN106504524A (en) * | 2016-09-14 | 2017-03-15 | 杭州诚智天扬科技有限公司 | A kind of method that express highway section is divided based on mobile signaling protocol dynamic |
| CN106535114A (en) * | 2016-09-29 | 2017-03-22 | 中国普天信息产业北京通信规划设计院 | Method and system for positioning terminal on high-speed railway and the monitoring method and system for high-speed railway network |
| CN106792943A (en) * | 2016-12-30 | 2017-05-31 | 广东海格怡创科技有限公司 | High ferro private network cell switches point detecting method and system |
| CN106878946A (en) * | 2017-02-16 | 2017-06-20 | 山东浪潮商用系统有限公司 | High ferro user identification method based on movement locus |
-
2018
- 2018-03-09 CN CN201810194087.9A patent/CN108495271B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103523056A (en) * | 2013-10-14 | 2014-01-22 | 兰州交通大学 | Train positioning method based on railway mobile communication network communication base stations |
| WO2015165682A1 (en) * | 2014-05-02 | 2015-11-05 | Mclaren Applied Technologies Limited | Location based connections |
| CN105657819A (en) * | 2014-11-12 | 2016-06-08 | 北京神州泰岳软件股份有限公司 | Railway mobile network test longitude and latitude obtaining method and device |
| CN105871448A (en) * | 2016-03-30 | 2016-08-17 | 中科凯普（天津）卫星导航通信技术有限公司 | High-dynamic transformation method of broadband special network of high-speed rail |
| CN106504524A (en) * | 2016-09-14 | 2017-03-15 | 杭州诚智天扬科技有限公司 | A kind of method that express highway section is divided based on mobile signaling protocol dynamic |
| CN106535114A (en) * | 2016-09-29 | 2017-03-22 | 中国普天信息产业北京通信规划设计院 | Method and system for positioning terminal on high-speed railway and the monitoring method and system for high-speed railway network |
| CN106792943A (en) * | 2016-12-30 | 2017-05-31 | 广东海格怡创科技有限公司 | High ferro private network cell switches point detecting method and system |
| CN106878946A (en) * | 2017-02-16 | 2017-06-20 | 山东浪潮商用系统有限公司 | High ferro user identification method based on movement locus |
Non-Patent Citations (2)
| Title |
|---|
| 高速铁路的GSM-R无线传播模型校正;肖蒙，刘佳佳等;《通信与网络》;20120601;全文 * |
| 高铁上终端的精准定位研究;王嘉兴，李果，卫钰;《信息通信》;20161227;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108495271A (en) | 2018-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101060689B (en) | A method and equipment for planning the communication system network | |
| US11800373B2 (en) | System and method for 3D propagation modelling for planning of a radio network | |
| CN108243451B (en) | Method and system for positioning in TD-LTE grid based on MR data | |
| CN102083090B (en) | Method and device for positioning interference source | |
| CN108207009B (en) | Rail network quality testing method and device | |
| CN102272617A (en) | Interference detection, characterization and location in a wireless communications or broadcast system | |
| CN101674523B (en) | Method and system for acquiring road information by using mobile communication network | |
| Spirito et al. | Preliminary experimental results of a GSM mobile phones positioning system based on timing advance | |
| CN106231621A (en) | A kind of many scene adaptives optimization method of propagation model in FDD LTE system | |
| CN106535133A (en) | Indoor telephone traffic accurate location method based on machine learning in cellular network | |
| CN108495271B (en) | Terminal positioning method suitable for high-speed railway network test | |
| Ibhaze et al. | Developing smart cities through optimal wireless mobile network | |
| Kan et al. | A dual-channel location estimation system for providing location services based on the GPS and GSM networks | |
| CN111491255B (en) | Backfill longitude and latitude positioning and expanding method, system and storage medium | |
| US20050003842A1 (en) | Network survey in radio telecommunications network | |
| CN108093414B (en) | Method and device for evaluating coverage effectiveness of cell | |
| CN110582058A (en) | wireless positioning method based on negative power summation ranging model | |
| CN106455052A (en) | Positioning method and device | |
| Tripkovic et al. | Benchmarking of mobile communications in high-speed scenarios: Active vs. passive modifications in high-speed trains | |
| CN108632740B (en) | Positioning method and device of user equipment | |
| Neuland et al. | Influence of positioning error on X-map estimation in LTE | |
| Tripkovic et al. | Measuring the Effects of AoA on Vehicle Penetration Loss in Cellular Networks | |
| Waadt et al. | An overview of positioning systems and technologies | |
| Khalife et al. | Characterization of sector clock biases in cellular CDMA systems | |
| Li et al. | Research on Differences of Prediction Models Related to Land Mobile Communication |
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 |