CN108495271A - A kind of method of locating terminal suitable for the test of high-speed railway network - Google Patents
A kind of method of locating terminal suitable for the test of high-speed railway network Download PDFInfo
- Publication number
- CN108495271A CN108495271A CN201810194087.9A CN201810194087A CN108495271A CN 108495271 A CN108495271 A CN 108495271A CN 201810194087 A CN201810194087 A CN 201810194087A CN 108495271 A CN108495271 A CN 108495271A
- Authority
- CN
- China
- Prior art keywords
- base station
- test
- curve
- gps
- speed railway
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- 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 present invention relates to fields of communication technology, are a kind of method of locating terminal suitable for the test of high-speed railway network, terminal positioning is tested when being lacked for communications industry main line of communication network test GPS signal.The present invention utilizes the azimuths operator eNodeB and high ferro map, test terminal positioning is carried out when train penetration loss causes GPS to lack in network test, the discretization track longitude and latitude point in eNodeB coverage areas, using the ratio of track longitude and latitude points and test terminal RSR sampling numbers as the interpolation period, interpolation track longitude and latitude point when GPS is lacked, as the location information of test terminal.This method is in the case of different GPS miss rates, i.e., when whole or part way GPS is lacked, the position of energy assignment test terminal does not increase equipment operating cost, while improving vehicle location efficiency without changing existing positioning device configuration.
Description
Technical field
The present invention relates to fields of communication technology, are a kind of method of locating terminal suitable for the test of high-speed railway network, use
Terminal positioning is tested when communications industry main line of communication network test GPS signal lacks.
Background technology
Network test is the method for radio network optimization data acquisition, is the downlink signal to telecom operators' wireless network
The acquisition of data is step essential in radio network optimization.Network test uses test equipment combining geographic information, right
Air interface is comprehensively measured, to verify the parameter of wireless interface, analysis, optimization and positioning failure.It is wireless by accessing
Network, using business such as the voice of telecom operators, data and internets, real-time perception network wireless parameter, quality of service, and
Preserve entire test process.The information of network test record can be output in processing software to analyze network quality situation.
Traditional test tool (for example, century ancient cooking vessel profit Pilot Pioneer drive test software) is used in high ferro network test,
Testing tool contains GPS positioning module, test software, terminal and computer.Test software records test eventually by interval of the sampling period
The gathered data at end and GPS.Normally, test terminal is about spaced 1s and reports 1 Reference Signal Received Power RSRP information, GPS
Report 1 latitude and longitude information.It tests terminal and GPS gathers data passes through association analysis, the RSRP values of each location point can be obtained.
High ferro closure aluminium alloy carriage body has very strong penetration loss to electromagnetic wave, and GPS signal is made to generate deep fades.For example, CRH is arranged
Car body penetration loss is about 17dB, and northwest region anti-blown sand, resistance to high and cold CRH5G train carriage body penetration loss are up to 24dB.
Since high ferro network test is by factors shadows such as train carriage body penetration loss, tunnel masking, massif blocking and U-type groove multipath fadings
It rings, causes GPS search satellite difficult and location information is easy serious loss.In addition, the influence of atmospheric duct makes in high ferro carriage body
Faint GPS signal is irregular to be sought.Under GPS deletion conditions, then it can not determine the coverage area of each eNodeB, make
The planning of the addresses eNodeB, network structure regulation lose location information support.It is each to lead to not assessment for the missing of location information
The communication quality of test point causes extreme difficulties to Neighborhood Optimization, parameter adjustment.
In recent years, some are suggested by the test method of locating terminal of attachment device, generally comprise gyroscope and acceleration
Degree meter positioning.They are to be combined into inertial nevigation apparatus with GPS, realize test terminal positioning.But telecom operators' centralized purchasing
Third party's drive test software and equipment, only meet the basic demand of outdoor road test scene, high ferro test scene can not be solved
GPS lacks problem.As high ferro operating mileage increases severely year by year, on the basis of existing outdoor test software and test equipment, reach
It is telecom operators and test equipment, test software to the real-time requirement for obtaining high ferro special screne and testing terminal positional information
The problem of supplier's urgent need to resolve.
Invention content
In view of the above technical problems, the present invention provides a kind of terminal positioning sides suitable for the test of high-speed railway network
Method carries out test terminal positional information using the communication engineering parameter cell azimuths eNodeB and the discrete longitude and latitude point of high ferro map
Compensation positioning.In order to solve the problems, such as high ferro network test GPS missings, test of the present invention in existing telecom operators' centralized purchasing
In Equipment Foundations, do not increase any test instrumentation, provides a kind of localization method of test terminal, expanded special screne test
Range, while reducing telecom operators and the additional investment of telecommunications equipment provider.
The technical solution adopted by the present invention to solve the technical problems is:
A kind of method of locating terminal suitable for the test of high-speed railway network includes the following steps
S1, tracing point substitute into base station azimuth and list antenna main lobe oriented radial equation
In land communication base station be denoted as eNodeB1, eNodeB2,, vehicle line be equipped with multiple gesture determination points
It is denoted as A1、A2,, A1、A2,, horizontal plane latitude and longitude coordinates be denoted as (x respectively1、y1)、(x2、y2),, any base station eNodeB
With the line of gesture determination point with positive north orientation angle be antenna azimuth it is α, base station eNodeB eyeball level coordinates difference
For (x, y), it is f that eyeball level coordinates, which are substituted into eyeball azimuth antenna main lobe oriented radial equation,1(x)=tan
(90°-a)·(x-x1)+y1;
S2, fitting high-speed railway rail curve
Actual measurement track circuit curve is obtained by open map, track circuit curve latitude coordinates parameter is brought into curve and intends
Close Equation f (x)=p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7;Utilize the Curve Fitting of traditional MATLAB
The parameter value that Tool calculates the fitting of 6 rank multinomials is 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 title of natural constant, about
2.71828 p1To p7For the fitting coefficient of 6 rank multinomials;
S3, the intersecting point coordinate for solving antenna main lobe oriented radial and orbital curve equation
Base station azimuth lists antenna main lobe oriented radial equation and solves intersecting point coordinate with Fitting curve equation,
When solving for only 1, this solution is the intersecting point coordinate of antenna main lobe oriented radial and high-speed railway rail matched curve, note
For (x'1, y'1), it is expressed as (x', y')=solve (f1(x),f(x));
When there is N number of solution, disaggregation is { (x'1j, y'1j), j=1,2 ... N are expressed as { (x'n, y'n) (1≤n≤N), N
A solution value is brought intoIt calculates nearest apart from base station
Intersecting point coordinate;
S4, according to intersecting point coordinate calculation base station switching point
Coordinate (the x' of switching pointB, y'B), substitute into adjoining nodes coordinate
Calculate the coordinate of base station switching point;
S5, it is inserted into intersecting point coordinate by the period
The GPS sampling numbers of terminal between two base station switching points are m, and M step S4 is inserted between two base station switching points
Intersecting point coordinate, the interpolation period isM >=m corrects high ferro location information by being inserted into intersecting point coordinate, prevents GPS from believing
Number missing.
Two base station switching points is smaller than 1.6km in the step S4.
The GPS sampling numbers of terminal in the step S5 between two base station switching points are m, and two base station switching points interleave
Enter the intersecting point coordinate of M step S4, the interpolation period isM > m.
Beneficial effects of the present invention are:
A kind of method of locating terminal suitable for the test of high-speed railway network, the present invention utilize operator eNodeB the (the 4th
Base station in third generation mobile communication network LTE) azimuth and high ferro map, when train penetration loss causes GPS to lack in network test
Test terminal positioning is carried out when mistake.Include the following steps:According to the azimuths eNodeB, antenna main lobe oriented radial is calculated;With height
On the basis of iron map, it is fitted high-speed railway rail curvilinear equation;Resolve the intersection point of ray and curve, i.e. ray equation and curvilinear equation
Solution;The switching point of adjacent eNodeB is found out according to intersection point;Coverage areas of the eNodeB on high-speed railway rail is calculated, as two are cut
Track between changing a little;The discretization track longitude and latitude point in eNodeB coverage areas is counted with track longitude and latitude and is tested eventually
Hold the ratio of RSRP (Reference Signal Receiving Power, Reference Signal Received Power) sampling number as slotting
Phase on weekly duty, interpolation track longitude and latitude point when GPS is lacked, as the location information of test terminal.This method is lacked in different GPS
In the case of mistake rate, i.e., when whole or part way GPS is lacked, the position of equal energy assignment test terminal is existing without changing
Positioning device configures, and does not increase equipment operating cost, while improving vehicle location efficiency.
Description of the drawings
Fig. 1 is the method for locating terminal model of the present invention;
Fig. 2 is the method for locating terminal flow of the present invention;
Fig. 3 is specific 1 contrast effect of case embodiment of the present invention;
Fig. 4 is specific 2 contrast effect of case embodiment of the present invention;
Fig. 5 is specific 3 contrast effect of case embodiment of the present invention;
Fig. 6 is the high-speed railway rail curve matching single order contrast effect of the present invention;
Fig. 7 is the high-speed railway rail curve matching second order contrast effect of the present invention;
Fig. 8 is the three rank contrast effect of high-speed railway rail curve matching of the present invention;
Fig. 9 is the high-speed railway rail curve matching quadravalence contrast effect of the present invention;
Figure 10 is the five rank contrast effect of high-speed railway rail curve matching of the present invention;
Figure 11 is the six rank contrast effect of high-speed railway rail curve matching of the present invention.
Specific implementation method
The present invention and its effect are expanded on further below in conjunction with drawings and examples.
A kind of method of locating terminal suitable for the test of high-speed railway network includes the following steps
S1, tracing point substitute into base station azimuth and list antenna main lobe oriented radial equation
In land communication base station be denoted as eNodeB1, eNodeB2,, vehicle line be equipped with multiple gesture determination points
It is denoted as A1、A2,, A1、A2,, horizontal plane latitude and longitude coordinates be denoted as (x respectively1、y1)、(x2、y2),, any base station eNodeB
With the line of gesture determination point with positive north orientation angle be antenna azimuth it is α, base station eNodeB eyeball level coordinates difference
For (x, y), it is f that eyeball level coordinates, which are substituted into eyeball azimuth antenna main lobe oriented radial equation,1(x)=tan
(90°-a)·(x-x1)+y1;
S2, fitting high-speed railway rail curve
Actual measurement track circuit curve is obtained by open map, track circuit curve latitude coordinates parameter is brought into curve and intends
Close Equation f (x)=p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7;Utilize the Curve Fitting of traditional MATLAB
The parameter value that Tool calculates the fitting of 6 rank multinomials is 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 title of natural constant, about
2.71828 p1To p7For the fitting coefficient of 6 rank multinomials;
S3, the intersecting point coordinate for solving antenna main lobe oriented radial and orbital curve equation
Base station azimuth lists antenna main lobe oriented radial equation and solves intersecting point coordinate with Fitting curve equation,
When solving for only 1, this solution is the intersecting point coordinate of antenna main lobe oriented radial and high-speed railway rail matched curve, note
For (x'1, y'1), it is expressed as (x', y')=solve (f1(x),f(x));
When there is N number of solution, disaggregation is { (x'1j, y'1j), j=1,2 ... N are expressed as { (x'n, y'n) (1≤n≤N), N
A solution value is brought intoIt calculates apart from base station most
Close intersecting point coordinate;
S4, according to intersecting point coordinate calculation base station switching point
Coordinate (the x' of switching pointB, y'B), substitute into adjoining nodes coordinate
Calculate the coordinate of base station switching point;
S5, it is inserted into intersecting point coordinate by the period
The GPS sampling numbers of terminal between two base station switching points are m, and M step S4 is inserted between two base station switching points
Intersecting point coordinate, the interpolation period isM >=m corrects high ferro location information by being inserted into intersecting point coordinate, prevents GPS from believing
Number missing.
Two base station switching points is smaller than 1.6km in the step S4.
The GPS sampling numbers of terminal in the step S5 between two base station switching points are m, and two base station switching points interleave
Enter the intersecting point coordinate of M step S4, the interpolation period isM > m.
Include two base station eNodeBs in Fig. 2, in secondary line, corresponds to two slave antennas respectively, direct north is clockwise with two
The angle of slave antenna main lobe direction is respectively β 1 and β 2, i.e. the azimuth of eNodeB2-1 cells is β 1, eNodeB2-2 cells
Azimuth is β 2;Switching point between eNodeB indicates that the serving cell that test terminal occupies before and after this position switches;
A1, A2, A3 and A4 are the intersection point of antenna main lobe oriented radial and orbital curve;B1 is the switching point of eNodeB1 and eNodeB2,
Test terminal occupies eNodeB1 radio resources in B1 left rails, and eNodeB2 radio resources are occupied in B1 right rails;
B2 is the switching point of eNodeB2 and eNodeB3, and test terminal occupies eNodeB2 radio resources in B2 left rails, on the right sides B2
ENodeB3 radio resources are occupied on siding track.The track of eNodeB2 coverings be can be seen that between B1, B2 point, it thus can be true
Determine the coverage area of each eNodeB along high ferro.If testing GPS signal when terminal occupies eNodeB2 to lack, B1, B2 point can be used
Between track longitude and latitude compensation missing location information.
Fig. 3 is shown in based on the azimuthal localization method flow of cell in high ferro VoLTE tests, steps are as follows:
1) antenna main lobe oriented radial equation is listed according to azimuth;
2) it is fitted high-speed railway rail curvilinear equation;
3) intersecting point coordinate of antenna main lobe oriented radial and orbital curve equation is solved;
If 4) intersection point number be more than 1, ask each intersection point between eNodeB at a distance from and selected distance minimum intersection point;If intersection point
Number is equal to 1, then executes next step;
4) according to intersecting point coordinate, the switching point of adjacent eNodeB is solved;
5) coverage areas of the eNodeB on high-speed railway rail, the track between as two switching points are calculated;
7) the discretization track longitude and latitude point in eNodeB coverage areas is counted with track longitude and latitude and tests terminal RSRP
The ratio of sampling number is as the interpolation period, interpolation track longitude and latitude point when GPS is lacked, and believes as the position of test terminal
Breath.
S1:Calculate high-speed railway rail curvilinear equation
Longitude and latitude point is turned to by high ferro map is discrete, and the distance between longitude and latitude point value is 5m.Railroad track curve side
Journey can be obtained by approximating method.Using MATLAB Curve Fitting Tool to Zhangye western station-Jiuquan southern station track
Figure carries out 6 rank multinomial fittings, and obtaining orbital curve equation is
F (x)=p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
The parameter value that the fitting of 6 rank multinomials is calculated using the Curve Fitting Tool of MATLAB is 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 be natural constant title, about 2.71828.p1To p7For the fitting coefficient of 6 rank multinomials.
S2:Calculate antenna in cell main lobe direction ray equation
Enable the azimuth point of eNodeB1-1 cells, eNodeB2-1 cells, eNodeB2-2 cells and eNodeB3-1 cells
It Wei not α, β1, β2And γ.Assuming that the two-dimensional position coordinate of eNodeB1, eNodeB2 and eNodeB3 are respectively (x1,y1)、(x2,y2)
(x3,y3).The then antenna main lobe of eNodeB1-1 cells, eNodeB2-1 cells, eNodeB2-2 cells and eNodeB3-1 cells
Oriented radial equation is shown in formula (3)~formula (6)
f1(x)=tan (90 ° of-a) (x-x1)+y1 (3)
f2(x)=tan (90 ° of-β1)·(x-x2)+y2 (4)
f3(x)=tan (90 ° of-β2)·(x-x2)+y2 (5)
f4(x)=tan (90 ° of-γ) (x-x3)+y3 (6)
When having and only 1 solves, this solution is that antenna main lobe oriented radial and the intersection point of high-speed railway rail matched curve are sat
Mark, is denoted as (x'1, y'1), (x'2, y'2), (x'3, y'3) and (x'4, y'4), it is represented by
Consider eNodeB1-1 antenna in cell main lobe direction ray equations f1(x) with the more intersection points of Gauss curve fitting orbital curve Equation f (x)
Situation.Disaggregation can be written as { (x'1j, y'1j), j=1,2 ... N.Assuming that { (x'1n, y'1n) (1≤n≤N) distance eNodeB1 is most
Closely, it is represented byThen retain n-th of solution, (x'1n,
y'1n) be intersection point A1 coordinate.Similarly, the coordinate of intersection point A2, A3 and A4 can be found out.
Ideally, switching point is equal at a distance from two adjacent intersection points.I.e. that is, in fig. 2 eNodeB1 and
The distance of switching point B1 range intersections A1, A2 of eNodeB2 is equal, and the switching point B2 distances of eNodeB2 and eNodeB3 are handed over
The distance of point A3, A4 are equal.Therefore, according to formula (7), the coordinate (x' of switching point B1, B2 can be found outB1, y'B1) and (x'B2, y'B2)
Track between switching point B1, B2 is coverage areas of the eNodeB2 on high-speed railway rail.Similarly, it can obtain
The coverage area of all eNodeB along high ferro.
The GPS sampling numbers of terminal between two base station switching points are m, are inserted into M intersection point between two base station switching points and sit
Mark, test terminal have lacked GPS position information in eNodeB2 coverage areas, and it is RSRP sampled point numbers to enable m, and M is on track
Discrete longitude and latitude between switching point B1, B2 is counted out, then the interpolation period beTherefore, with the extraction of interpolation cycle T
Discrete longitude and latitude point between switching point B1, B2 can be used as location information of the test terminal in eNodeB2 coverage areas, to
Realize test terminal positioning when any sampling time GPS missings.
Embodiment 1
Assignment test terminal of the present invention occupies location information when eNodeB1.The coordinate of eNodeB1 be (99.346405,
39.323753), the azimuth of eNodeB1-1 cells is 105 °, and the azimuth of eNodeB1-2 cells is 265 °.Two antenna in cell
Main lobe direction ray equation is respectively:
f2(x)=tan (90 ° -105 °) (x-99.346405)+39.323753 and f3(x)=tan (90 ° -265 °)
(x-99.346405)+39.323753.Neighboring community's coordinate of the adjacent eNodeB of eNodeB1-1 cell southeastern directions is
(99.361522,39.322955), azimuth are 260 °, and antenna main lobe oriented radial equation is:f1(x)=tan (90 °-
260°)·(x-99.361522)+39.322955.The adjacent neighboring communities the eNodeB coordinate of eNodeB1-2 cell direction northwests is
(99.333510,39.325803), azimuth are 120 °, and antenna main lobe oriented radial equation is:f4(x)=tan (90 °-
120°)·(x-99.333510)+39.325803。
The 6 rank multinomial fit curve equations at Lanzhou-Xinjiang high ferro Zhangye western station-Jiuquan southern station are:
F (x)=p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
The parameter value that the fitting of 6 rank multinomials is calculated using the Curve Fitting Tool of MATLAB is 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 natureConstantTitle, about 2.71828.P1 to P 7 is the fitting coefficient of 6 rank multinomials.
f1(x)、f2(x)、f3(x) and f4(x) respectively with the intersection point of f (x) be A1 (99.357261,39.322252), A2
(99.348314,39.323287), A3 (99.345036,39.323681) and A4 (99.334996,39.324936).On track
With A1 (99.357261,39.322252) and the equidistant points of A2 (99.348314,39.323287) be B1 (99.352703,
39.322498), equidistant with A3 (99.345036,39.323681) and A4 (99.334996,39.324936) on track
Point is B2 (99.340070,39.324550).I.e. B1 (99.352703,39.322498) and B2 (99.340070,
39.324550) coverage area in high ferro that the track between is eNodeB1.Test terminal, which occupies eNodeB2,17
RSRP sampled points, B1 (99.352703,39.322498) and B2 (99.340070,39.324550) uniformly takes 17 on high-speed railway rail
A longitude and latitude point, i.e., longitude and latitude is listed as follows when occupying eNodeB1 using assignment test terminal of the present invention:
Serial number | Position longitude | Position latitude |
1 | 99.352183 | 39.322518 |
2 | 99.351473 | 39.322547 |
3 | 99.350763 | 39.322575 |
4 | 99.350054 | 39.322603 |
5 | 99.349344 | 39.322632 |
6 | 99.348634 | 39.322779 |
7 | 99.347924 | 39.322925 |
8 | 99.347215 | 39.323072 |
9 | 99.346505 | 39.323219 |
10 | 99.345795 | 39.323366 |
11 | 99.345085 | 39.323513 |
12 | 99.344376 | 39.323660 |
13 | 99.343666 | 39.323806 |
14 | 99.342956 | 39.323953 |
15 | 99.342247 | 39.324100 |
16 | 99.340922 | 39.324374 |
17 | 99.340259 | 39.324511 |
Embodiment 2
Assignment test terminal of the present invention occupies location information when eNodeB2.The coordinate of eNodeB2 be (99.169312,
39.359299), the azimuth of eNodeB2-1 cells is 90 °, and the azimuth of eNodeB2-2 cells is 290 °.Two antenna in cell
Main lobe direction ray equation is respectively:f2(x)=tan (90 ° -90 °) (x-99.169312)+39.359299 and f3(x)=
tan(90°-295°)·(x-99.169312)+39.359299.The adjacent eNodeB's of eNodeB2-1 cell southeastern directions is adjacent
Cell coordinate is (99.184611,39.359588), and azimuth is 265 °, and antenna main lobe oriented radial equation is:f1(x)=
tan(90°-265°)·(x-99.184611)+39.359588.The adjacent eNodeB of eNodeB2-2 cell direction northwests is adjacent small
Area's coordinate is (99.155389,39.361827), and azimuth is 90 °, and antenna main lobe oriented radial equation is:f4(x)=tan
(90°-90°)·(x-99.155389)+39.361827。
The 6 rank Gauss curve fitting curvilinear equations at Lanzhou-Xinjiang high ferro Zhangye western station-Jiuquan southern station are:
F (x)=p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
The parameter value that the fitting of 6 rank multinomials is calculated using the Curve Fitting Tool of MATLAB is 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 natureConstantTitle, about 2.71828.p1To p7For the fitting coefficient of 6 rank multinomials.
f1(x)、f2(x)、f3(x) and f4(x) respectively with the intersection point of f (x) be A1 (99.172832,39.35857), A2
(99.170402,39.359301), A3 (99.167242,39.360259) and A4 (99.162254,39.361802).On track
With A1 (99.172832,39.35857) and the equidistant points of A2 (99.170402,39.359301) be B1 (99.171959,
39.360699), equidistant with A3 (99.167242,39.360259) and A4 (99.162254,39.361802) on track
Point is B2 (99.164958,39.362086).I.e. B1 (99.171959,39.360699) and B2 (99.164958,
39.362086) coverage area in high ferro that the track between is eNodeB2.Test terminal, which occupies eNodeB2,9
RSRP sampled points, B1 (99.171959,39.360699) and B2 (99.164958,39.362086) uniformly takes 9 on high-speed railway rail
A longitude and latitude point, that is, when using localization method of the present invention, test terminal is listed as follows when occupying eNodeB2 longitudes and latitudes:
Serial number | Position longitude | Position latitude |
1 | 99.171724 | 39.360746 |
2 | 99.170967 | 39.360897 |
3 | 99.170210 | 39.361046 |
4 | 99.169453 | 39.361195 |
5 | 99.168696 | 39.361346 |
6 | 99.167939 | 39.361497 |
7 | 99.167182 | 39.361644 |
8 | 99.166424 | 39.361796 |
9 | 99.165667 | 39.361948 |
Embodiment 3
Assignment test terminal of the present invention occupies location information when eNodeB3.The coordinate of eNodeB3 be (98.976210,
39.443899), the azimuth of 1 cell is 160 °, and the azimuth of 2 cells is 290 °.Two antenna in cell main lobe direction ray equations
Respectively:f2(x)=tan (90 ° -160 °) (x-98.976210)+39.443899 and f3(x)=tan (90 ° -290 °)
(x-98.976210)+39.443899.Neighboring community's coordinate of the adjacent eNodeB of eNodeB3-1 cell southeastern directions is
(98.982611,39.39.435891), azimuth are 300 °, and antenna main lobe oriented radial equation is:f1(x)=tan (90 °-
300°)·(x-98.982611)+39.435891.The adjacent neighboring communities the eNodeB coordinate of eNodeB3-2 cell direction northwests is
(98.968178,39.449521), azimuth are 155 °, and antenna main lobe oriented radial equation is:f4(x)=tan (90 °-
155°)·(x-98.968178)+39.449521。
The 6 rank Gauss curve fitting curvilinear equations at Lanzhou-Xinjiang high ferro Zhangye western station-Jiuquan southern station are:
F (x)=p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7。
The parameter value that the fitting of 6 rank multinomials is calculated using the Curve Fitting Tool of MATLAB is 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 natureConstantTitle, about 2.71828.p1To p7For the fitting coefficient of 6 rank multinomials.
f1(x)、f2(x)、f3(x) and f4(x) respectively with the intersection point of f (x) be A1 (98.982611,39.435891), A2
(98.976975,39.44177), A3 (98.967195,39.447177) and A4 (98.970015,39.445607).On track
With A1 (98.982611,39.435891) and the equidistant points of A2 (98.976975,39.44177) be B1 (98.981471,
39.440866), equidistant with A3 (98.967195,39.447177) and A4 (98.970015,39.445607) on track
Point is B2 (98.967939,39.451956).I.e. B1 (98.981471,39.440866) and B2 (98.967939,
39.451956) coverage area in high ferro that the track between is eNodeB3.Test terminal, which occupies eNodeB2,9
RSRP sampled points, B1 (98.981471,39.440866) and B2 (98.967939,39.451956) uniformly takes 32 on high-speed railway rail
A longitude and latitude point, that is, when using localization method of the present invention, longitude and latitude is listed as follows when test terminal occupies eNodeB3:
Real case 1 is compared with traditional technology
Case 1, terminal occupy about 935 meters of the measuring distance of eNodeB1 signals, manually along about 55 meters of acquisitions of high-speed railway rail
1 GPS information, acquisition latitude and longitude information is as follows, contrast effect such as Fig. 4.
Real case 2 is compared with traditional technology
Case 2, terminal occupy about 495 meters of the measuring distance of eNodeB1 signals, manually along about 55 meters of acquisitions of high-speed railway rail
1 GPS information, acquisition latitude and longitude information is as follows, contrast effect such as Fig. 5.
Real case 3 is compared with traditional technology
Case 3, terminal occupy eNodeB1 signals about 1760 meters of measuring distance, manually along high-speed railway rail about 55 meters adopt
Collect 1 GPS information, acquisition latitude and longitude information is as follows, contrast effect such as Fig. 6.
Claims (3)
1. a kind of method of locating terminal suitable for the test of high-speed railway network, it is characterised in that include the following steps
S1, tracing point substitute into base station azimuth and list antenna main lobe oriented radial equation
In land communication base station be denoted as eNodeB1, eNodeB2,, vehicle line is denoted as equipped with multiple gesture determination points
A1、A2,, A1、A2,, horizontal plane latitude and longitude coordinates be denoted as (x respectively1、y1)、(x2、y2),, any base station eNodeB and rail
It is α that the line of mark measuring point and positive north orientation angle, which are antenna azimuth, base station eNodeB eyeball level coordinates be respectively (x,
Y), it is f eyeball level coordinates to be substituted into eyeball azimuth antenna main lobe oriented radial equation1(x)=tan (90 °-
a)·(x-x1)+y1;
S2, fitting high-speed railway rail curve
Actual measurement track circuit curve is obtained by open map, brings track circuit curve latitude coordinates parameter into curve matching side
Journey f (x)=p1*x6+p2*x5+p3*x4+p4*x3+p5*x2+p6*x+p7;Utilize the Curve Fitting Tool of traditional MATLAB
The parameter value for calculating the fitting of 6 rank multinomials is 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 natural constant, about 2.71828, p1To p7For
The fitting coefficient of 6 rank multinomials;
S3, the intersecting point coordinate for solving antenna main lobe oriented radial and orbital curve equation
Base station azimuth lists antenna main lobe oriented radial equation and solves intersecting point coordinate with Fitting curve equation,
When solving for only 1, this solution is the intersecting point coordinate of antenna main lobe oriented radial and high-speed railway rail matched curve, is denoted as
(x′1, y '1), it is expressed as (x', y')=solve (f1(x),f(x));
When there is N number of solution, disaggregation is { (x '1j, y '1j), j=1,2 ... N are expressed as { (x 'n, y 'n) (1≤n≤N), N number of solution
Value is brought intoThe intersection point nearest apart from base station is calculated to sit
Mark;
S4, according to intersecting point coordinate calculation base station switching point
Coordinate (the x ' of switching pointB, y 'B), substitute into adjoining nodes coordinate
Calculate the coordinate of base station switching point;
S5, it is inserted into intersecting point coordinate by the period
The GPS sampling numbers of terminal between two base station switching points are m, and the friendship of M step S4 is inserted between two base station switching points
Point coordinates, interpolation period beM >=m corrects high ferro location information by being inserted into intersecting point coordinate, GPS signal is prevented to lack
It loses.
2. a kind of method of locating terminal suitable for the test of high-speed railway network according to claim 1, it is characterised in that
Two base station switching points is smaller than 1.6km in the step S4.
3. a kind of method of locating terminal suitable for the test of high-speed railway network according to claim 1, it is characterised in that
The GPS sampling numbers of terminal in the step S5 between two base station switching points are m, and M step is inserted between two base station switching points
Suddenly the intersecting point coordinate of S4, interpolation period areM > 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 true CN108495271A (en) | 2018-09-04 |
CN108495271B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111288942A (en) * | 2018-12-07 | 2020-06-16 | 比亚迪股份有限公司 | Track transponder position measuring method and device and computer equipment |
CN111340908A (en) * | 2020-02-10 | 2020-06-26 | 兰州交通大学 | Track 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无线传播模型校正", 《通信与网络》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111288942A (en) * | 2018-12-07 | 2020-06-16 | 比亚迪股份有限公司 | Track transponder position measuring method and device and computer equipment |
CN111288942B (en) * | 2018-12-07 | 2021-09-03 | 比亚迪股份有限公司 | Track transponder position measuring method and device and computer equipment |
CN111340908A (en) * | 2020-02-10 | 2020-06-26 | 兰州交通大学 | Track electronic map generation method |
CN111340908B (en) * | 2020-02-10 | 2023-07-25 | 兰州交通大学 | Rail electronic map generation method |
Also Published As
Publication number | Publication date |
---|---|
CN108495271B (en) | 2020-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105611568B (en) | A method of LTE terminal position is accurately positioned based on MRO measurement report | |
CN102149192B (en) | Cellular network wireless positioning method based on cooperation of mobile stations | |
CN102272617B (en) | Interference detection, characterization and location in a wireless communications or broadcast system | |
CN104735705B (en) | A kind of base station covering and the check method and device of basic data | |
CN108207009B (en) | Rail network quality testing method and device | |
CN106535133A (en) | Indoor telephone traffic accurate location method based on machine learning in cellular network | |
CN101674523B (en) | Method and system for acquiring road information by using mobile communication network | |
CN109905892A (en) | Method for positioning user and device | |
CN105491586B (en) | Cell-site antenna orientation angles measurement method and system | |
CN111416676B (en) | High-speed rail railway crossing and merging section field strength prediction method based on ray tracking | |
CN105629080A (en) | Drawing method of electromagnetic distribution situation diagram or path electromagnetic distribution diagram | |
CN108495271A (en) | A kind of method of locating terminal suitable for the test of high-speed railway network | |
CN105501464A (en) | Method for testing performance of aviation airborne radio communication and navigation system during test flight | |
Jeong et al. | RSS-based LTE base station localization using single receiver in environment with unknown path-loss exponent | |
CN103323856A (en) | Method for locating non-cooperative wireless signal source based on high-rail three-satellite time-difference system | |
CN107607908A (en) | A kind of localization method based on signal time delay feature spatial and temporal distributions | |
CN104199072B (en) | Space flight thing is found and track tagging system and its method for realizing track mark | |
CN109996259A (en) | Aerial angle method of adjustment, device, equipment and medium | |
CN113347587B (en) | Subcarrier diversity precision ranging positioning method based on multicarrier signals | |
Jan et al. | Evaluation of positioning algorithms for wide area multilateration based alternative positioning navigation and timing (APNT) using 1090 MHz ADS-B signals | |
Šafář | Analysis, modelling and mitigation of cross-rate interference in enhanced loran | |
CN112363190A (en) | Different-rail time-sharing single-satellite time difference positioning method | |
CN108646218A (en) | A kind of lightning location net precision assessment method measured based on dynamic | |
CN108490393A (en) | A kind of positioning of marine mobile platform and air navigation aid based on lightning location net | |
CN107422340A (en) | A kind of multipoint location system receiving station positional error correction method |
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 |