CN102056224A - Method for positioning mobile terminal in cellular network - Google Patents

Method for positioning mobile terminal in cellular network Download PDF

Info

Publication number
CN102056224A
CN102056224A CN2010106061782A CN201010606178A CN102056224A CN 102056224 A CN102056224 A CN 102056224A CN 2010106061782 A CN2010106061782 A CN 2010106061782A CN 201010606178 A CN201010606178 A CN 201010606178A CN 102056224 A CN102056224 A CN 102056224A
Authority
CN
China
Prior art keywords
sub
district
group
base station
measurement report
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
CN2010106061782A
Other languages
Chinese (zh)
Other versions
CN102056224B (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN 201010606178 priority Critical patent/CN102056224B/en
Publication of CN102056224A publication Critical patent/CN102056224A/en
Application granted granted Critical
Publication of CN102056224B publication Critical patent/CN102056224B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention provides a method for positioning a mobile terminal in a cellular network. The method comprises the following steps of: generating measurement reports during communication of a terminal side and a network side and transmitting the measurement reports through an ABIS interface or an Iub interface of the network side; calculating an approximate geographic position of a target terminal by an approximate method according to geographic position, reception level and type of a cell included by each measurement report in the group of measurement report data records for the measurement report; calculating positioning error of the approximate geographic position; correcting the positioning error by iterative computation to acquire a group of geographic position data; and arithmetically averaging longitude data and latitude data of the group of geographic position data to acquire the final target terminal geographic position. The method can more accurately position the mobile terminal than the conventional method taking the cell as the terminal position, has low complexity, and can be applied to the second generation and the third generation mobile communication cellular networks such as global system for mobile communication (GSM), code division multiple access (CDMA), CDMA2000, wideband code division multiple access (WCDMA), time division-synchronization code division multiple access (TD-SCDMA) and the like.

Description

The method of localisation of mobile terminals in a kind of cellular network
Technical field
The present invention relates to moving communicating field, relate in particular to the method for localisation of mobile terminals in the cellular network.
Background technology
Development along with mobile communication technology, location-based service extensively is penetrated in people's the productive life, utilizes the terminal positioning technology, the relevant information of can inquiring user required locality, as neighbouring which hotel, market, gas station are arranged, and the position of local weather situation and certain company etc.; When people run into the emergency traffic accident in the other places, rob, when sick and wounded, as long as beat emergency relief telephone with mobile phone, cellular network is when delivering to rescue center with this urgent call, can send terminal use's positional information and user's voice information to rescue center in the lump, with fast reaction, the success rate that improves rescue and fight crime; Adorn a moving carriage-table on each automobile or the every row train, communications and transportation unit can determine the position of vehicle, realizes management of vehicles and scheduling easily; By the long-time statistics of great amount of terminals positional information, can also provide significant data in addition for rational deployment, structure and the raising wireless resource utility efficiency of honeycomb.
The portable terminal location technology is divided into satellite system location and cellular network location two big classes at present, the utilization of satellite system location technology runs on the multi-satellite of trapped orbit and the Time Calculation end coordinates that satellite-signal arrives receiver, accuracy is higher, but terminal can't receive signal when the user is positioned at city or shelter, and positioning chip and algorithm need by the terminal support, if use the terminal that chip is not installed then can't enjoy location-based service; Simultaneously except user self, other personnel are difficult to learn user's position, and when the user was in dangerous situation or awkward situation and is difficult to send Location Request, board service personnel can't the positive location user, and these all provide difficulty for the popularization of location-based service.And the cellular network location technology is provided by Mobile Network Operator, its wide coverage, can locate any zone that mobile network signals covers, need not the user changes or upgrading terminals, can be by the disaster-stricken user of network management personnel's positive location, and locating information is often combined closely with data such as the intensity of base station signal, quality, frequency, user's scramblers, can be applied in access network preferably and optimize the field.
The angle of arrival of signal of communication between cellular network location technology general using base station and the terminal, the time of advent, the time of advent parameter acquiring terminal such as difference and received power the position, yet traditional is lower as the localization method precision of terminal location with the sub-district, and amount of information is less, is difficult to be applied to access network and optimizes the field.
Summary of the invention
The invention provides the method for localisation of mobile terminals in a kind of cellular network, utilize end side with the network service process in produce and measurement report consumer positioning positioning accuracy by network side ABIS interface or Iub interface transmission is higher than traditional approximate location method, and complexity is low, can be applied in the second generations such as GSM, CDMA, CDMA2000, WCDMA, TD-SCDMA and the 3G (Third Generation) Moblie cellular network.
For solving the problems of the technologies described above, the present invention proposes the method for localisation of mobile terminals in a kind of cellular network, it is characterized in that, may further comprise the steps:
Step 1, utilize the cellular network interface snooping equipment to intercept interface between base station and the base station controller, to obtain group measurement reporting data of corresponding target terminal;
Step 2, obtain the engineering parameter table from network operation maintenance centre;
Step 3, each sub-district that each the bar measurement report in these group measurement reporting data is comprised, in the engineering parameter table, carry out match search, with the geographical position that obtains this sub-district, affiliated base station, judge whether this sub-district corresponding parameters is abandoned according to the geographical position of this sub-district;
Step 4, administer and determined the type of this sub-district according to base station under this sub-district by the number of cells of this measurement report record;
Step 5 utilizes approximation method to calculate the general geographic location at target terminal place;
Step 6 is to each sub-district that each the bar measurement report in these group measurement reporting data comprises, the distance between the geographical position of calculation plot and the general geographic location of target terminal;
Step 7, to each the bar measurement report in these group measurement reporting data, calculate the distance between each sub-district, the arithmetic average statistical expection and the mean square deviation of distance between the calculation plot, calculate the arithmetic average statistical expection and the mean square deviation of the corresponding incoming level in each sub-district, calculate the position error of described general geographic location according to the number of above-mentioned statistical index and sub-district;
Step 8, to each the bar measurement report in these group measurement reporting data, affiliated base station title according to its sub-district that comprises, all microzonations are divided into several groups, it is all inequality that sub-district in each group all belongs to the base station title that same base station is administered and each group is corresponding, according to incoming level, type, the antenna directional angle of all sub-districts in each group, calculate an angle value of affiliated base station, this group sub-district correspondence;
Step 9, to each the bar measurement report in these group measurement reporting data, according to the distance between the general geographic location of the general geographic location of the received signal level of its each sub-district that comprises, target terminal, each sub-district that it comprises and target terminal, the position error of approximation method, some groups, the corresponding angle value of each group of division, carry out iterative computation, obtain one group of target terminal geographic position data;
Step 10 at one group of terminal geographic position data that step 9 obtains, is carried out the arithmetic average statistical average of longitude data and latitude data respectively, obtains final objective terminal geographic position.
Further, described step 3 specifically may further comprise the steps:
Step a, the Serving cell that comprises for each bar measurement report data record, in the engineering parameter table, carry out match search, matching process with the geographical position that obtains this Serving cell, affiliated base station is: search the data line record that has comprised the same cells sign in the engineering parameter table, as the engineering parameter data record of this Serving cell correspondence;
Step b, each adjacent area that comprises for each bar measurement report data record, in the engineering parameter table, carry out match search, matching process with the geographical position that obtains this sub-district, affiliated base station is: for gsm system or cdma system, search the data record row that comprises identical broadcasts control channel number and base station identity code in the engineering parameter table; For CDMA2000, WCDMA or TD-SCDMA system, in the engineering parameter table, search the data record row that comprises identical common signal channel center frequency point and scrambler;
Step c, each adjacent area that comprises for each bar measurement report data record, when the record line number order of the coupling that searches greater than 1 the time, geographical position according to the sub-district in the geographical position of the Serving cell of this measurement report data record and many data record in the engineering parameter table, matching, distance between difference calculation services sub-district and the data record respective cell, that data line record of a selected distance value minimum then is as the engineering parameter data record of this adjacent area correspondence;
Steps d, each adjacent area that comprises for each bar measurement report data record, according to the geographical position of the Serving cell of this measurement report data record and the geographical position of this adjacent area, the distance between calculation services sub-district and the adjacent area, when this distance greater than f dWhen rice, then in thereafter step S07, the sub-district described in the S08, get rid of this adjacent area; f dBe empirical coefficient, value is a 2000-3500 rice.Described f dBe preferably 2640 meters.
Described step 4 specifically comprises: write down each adjacent area that comprises for each bar measurement report data, when the record line number order of the coupling that searches greater than 1 the time, geographical position according to the sub-district in the geographical position of the Serving cell of this measurement report data record and many data record in the engineering parameter table, matching, distance between difference calculation services sub-district and the data record respective cell, that data line record of a selected distance value minimum then is as the engineering parameter data record of this adjacent area correspondence.
The approximation method of described step 5 comprises: to each the bar measurement report data record in this group measurement reporting, its all subdistrict positions that comprise are carried out the statistical average of longitude and latitude respectively, perhaps the received signal level data of sub-district are carried out the weighting statistical average of longitude and latitude as weight coefficient respectively to its all subdistrict positions that comprise, then the one group of geographic position data that calculates is carried out the statistical average of longitude data and latitude data respectively.In the above-mentioned steps, for CDMA2000, WCDMA and TD-SCDMA system, the received signal level of any one sub-district is P-CCPCH received signal code power+110 of this sub-district.
In the described step 7, the position error that calculates approximation method specifically may further comprise the steps:
s = 10 f s × E ( d ) f smd × D ( d ) f sdd × E ( R ) f smr × D ( R ) f sdr × N f sn - - - ( 1 )
In the formula (1), E (d) represents the arithmetic average statistical expection of distance between each sub-district, and d represents distance between each sub-district; D (d) represents the arithmetic average statistics mean square deviation of distance between each sub-district; E (R) represents the arithmetic average statistical expection of the corresponding incoming level in each sub-district, and R represents the corresponding incoming level in each sub-district; D (R) represents the arithmetic average statistics mean square deviation of the corresponding incoming level in each sub-district; N represents the number of sub-district; f s, f Smd, f Sdd, f Smd, f Smr, f SdrAnd f SnAll be empirical coefficient, have following value respectively:
f s=0.8±0.2、f smd=0.8±0.2、f sdd=0.1±0.025、f smr=-0.1±0.025、f sdr=-0.3±0.075、f sn=-0.75±0.18。
Preferably, described f s=0.8, f Smd=0.8, f Sdd=0.1, f Smr=-0.1, f Sdr=-0.3, f Sn=-0.75.
Described step 8 specifically may further comprise the steps:
To each group sub-district that each the bar measurement report in these group measurement reporting data comprises, when this group sub-district belonged to the 1st class sub-district, this group only comprised a sub-district, and the antenna directional angle of this sub-district is the angle value of affiliated base station, this group sub-district correspondence.
When this group sub-district belongs to the 2nd class sub-district, this group comprises two sub-districts, if these two sub-districts are respectively sub-district 1 and sub-district 2, when the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 during, adopt following formula (2) to calculate the angle value of base station correspondence under this group sub-district less than 180 °:
Φ=(θ 12)/2+s(θ 12)×f Φ×(ε 12) (2)
Wherein, θ 1Represent the antenna directional angle of sub-district 1; θ 2Represent the antenna directional angle of sub-district 2; ε 1Represent the received signal level of sub-district 1; ε 2Represent the received signal level of sub-district 2; f ΦBe empirical coefficient, its value is 3.35 ± 0.85.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 during, adopt following formula (3) to calculate the angle value of base station correspondence under this group sub-district greater than 180 °:
Φ=(360+θ 12)/2+s(θ 12)×f Φ×(ε 21) (3)
Wherein, each variable to gush justice identical with the connotation of variable of the same name in the formula (2).
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 equals 180 °, then should choose the angle value of the antenna directional angle of a bigger sub-district of corresponding received power as base station correspondence under this group sub-district.
When this group sub-district belongs to the 3rd class sub-district, if this group be numbered 1, group 1 comprises 3 sub-districts, establishes these 3 sub-districts and is respectively sub-district 3, sub-district 4 and sub-district 5, and the incoming level of sub-district 3 is greater than the incoming level of sub-district 4, and the incoming level of sub-district 4 is greater than the incoming level of sub-district 5; Then sub-district 3 and sub-district 4 are considered as the 2nd class sub-district, and are combined into an interim group, that establishes this interim group is numbered 2, calculates the angle value of group 2 correspondences; To organize 2 then and be considered as an interim sub-district, establishing this interim sub-district is sub-district 6, and the incoming level of sub-district 6 is the mean value of the incoming level of sub-district 3 and sub-district 4, and the antenna directional angle of sub-district 6 is the angle value of group 2 correspondences.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 during, adopt following formula (4) to calculate the angle value of base station correspondence under the sub-district of group 1 less than 180 °:
Φ=(2×θ 43)/3+4s(θ 34)×f Φ×(ε 34)/3 (4)
Wherein, θ 3Represent the antenna directional angle of sub-district 5; θ 4Represent the antenna directional angle of sub-district 6; ε 3Represent the received signal level of sub-district 5; ε 4Represent the received signal level of sub-district 6.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 during, can adopt following formula (5) to calculate the angle value of base station correspondence under the sub-district of group 1 greater than 180 °:
Φ = [ 2 × θ 4 + ( θ 3 + 360 ) ] / 3 - 4 × f Φ × ( ϵ 4 - ϵ 3 ) / 3 , θ 4 > θ 3 [ 2 × ( θ 4 + 360 ) + θ 3 ] / 3 + 4 × f Φ × ( ϵ 4 - ϵ 3 ) / 3 , θ 3 > θ 4 - - - ( 5 )
Wherein, the connotation of the variable of the same name in the connotation of each variable and the formula (2) is identical.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 equals 180 °, then should choose the angle value of the antenna directional angle of a bigger sub-district of corresponding received power as base station correspondence under the sub-district of group 1.
Wherein, in formula (2), formula (3) and formula (4), s (θ 12) be sign function, its implication can be represented with following formula (6):
s ( &alpha; ) = 1 , &alpha; > 0 0 , &alpha; = 0 - 1 , &alpha; < 0 - - - ( 6 )
Described step 9 specifically may further comprise the steps:
If the E point is the general geographic location of target terminal, A, B, 4 difference of C, D are represented the geographical position of the affiliated base station, sub-district in pairing each group of this measurement report successively, a, b, four rays of c, d are represented the angle value of the above-mentioned base station correspondence that obtains respectively successively in step S11, the angular range that the dotted line representative is certain; During iterative computation, being the center of circle with the A point at first, is radius with AE, and the E point is rotated to an angle to the direction near ray a, arrives the F point; Then be the center of circle, be radius with BF, the F point is rotated to an angle to the direction near ray b, arrive the G point with the B point; Then be the center of circle, be radius with CG, the G point is rotated to an angle to the direction near ray c, arrive the H point with the C point; Then be the center of circle, be radius with DH, the G point is rotated to an angle to the direction near ray d, arrive the I point with the D point; Then, be the center of circle with A again, be radius rotation with AI; Before once rotating arbitrarily, if point to be rotated, for example, E point, F point, G point or H point, the orientation of pivot is in beyond the angular range of two dotted lines that pivot next time connects relatively next time, is not that the center of circle is rotated with this center then, and skips, whether the geographical position with next base station is the center of circle, rejudge and rotate.
Total number of revolutions is between 10 to 30 times, and the amplitude of each rotation satisfies following formula (7):
Figure BDA0000040574180000061
Wherein,
Figure BDA0000040574180000062
Be the center of circle to be rotated (under any one group of sub-district that comprises the geographical position of base station) with point (the interim result of iterative computation, preliminary examination value are the general geographic location of target terminal) to be rotated between distance; R iIt is the assembly average of organizing the received signal level of the sub-district that comprises; S is the position error of the approximation method that obtains in step S10; N is the number of organizing in this measurement report; Clip is the azimuth in point the to be rotated center of circle for the treatment of rotation mutually and the angle between the rolling target ray; f v, f Vd, f Vr, f VsAnd f VnAll be empirical coefficient, they generally have following value:
f v=-6.7±1.67、f vd=-1.5±0.37、f vr=3.35±0.8、f vs=1.85±0.46、f vn=-1±0.25。
Preferably, described f v=-6.7, f Vd=-1.5, f Vr=3.35, f Vs=1.85, f Vn=-1.
Complexity of the present invention is low, and positioning accuracy is moderate, under the sub-district mean radius is situation about 550 meters, field test according to all kinds of complex-terrains such as road, square, the foot of the hill, lakeside and building sides, average position error can reach about 250 meters, suitable general to positioning accuracy request, but the bigger occasion of deal with data amount.
Description of drawings
Accompanying drawing described herein has been used to be provided further understanding of the present invention, constitutes the application's a part, and illustrative examples of the present invention and explanation thereof are used to explain the present invention, do not constitute improper qualification of the present invention.Wherein:
Fig. 1 is the flow chart according to the method for localisation of mobile terminals in the cellular network of the embodiment of the invention;
Fig. 2 be basis of the present invention roughly terminal location progressively rotate process schematic diagram with the more accurate terminal location of final acquisition.
Embodiment
In the following description, provided some detail so that fully understanding to a plurality of embodiment of the present invention.But, those skilled in the art do not have these details can realize the present invention as can be known yet.In other cases, there are not detailed expression or description and the relevant known configurations of equipment such as computer, transmission medium, measuring equipment, recording equipment, to avoid unnecessary interference description to the embodiment of the invention.
Unless content needs, otherwise in specification and claims, word " comprises " and change the meaning of " comprising " of being understood as that evolution, promptly " includes, but are not limited to ".
About the description of " sub-district ", " Serving cell " and " adjacent area " notion, based on the general knowledge of this area, the three has different implications in the literary composition, and " sub-district " extension of conception has included " Serving cell " and " adjacent area ".
The present invention proposes the method for localisation of mobile terminals in a kind of cellular network, improved positioning accuracy, and can be applied in the cellular network of multiple standard.As shown in Figure 1, the method according to localisation of mobile terminals in the cellular network of the embodiment of the invention may further comprise the steps:
Step S01 utilizes the cellular network interface snooping equipment, as Signaling Analyzer, intercepts the interface between base station and the base station controller, and to obtain a group measurement reporting of corresponding target terminal, the data format of measurement report as shown in Figure 2.The data format of " measurement report that MS reports " among Fig. 2 as shown in Figure 3.One group measurement reporting comprises the measurement report data record of the constant duration record of some (being generally less than 60); For GSM or cdma system, the data record of every measurement report comprises following parameter: timestamp, Serving cell sign, Serving cell received signal level, the Broadcast Control Channel of several (maximum 6, minimum 0) adjacent areas number, base station identity code, received signal level.
CDMA2000, WCDMA or TD-SCDMA system for 3G, every the measurement report data record comprises following parameter: timestamp, the Serving cell sign, the P-CCPCH received signal code power of Serving cell, common signal channel center frequency point, scrambler, the P-CCPCH received signal code power of several (maximum 6, minimum 0) adjacent areas.
For GSM or cdma system, the interface between base station and the base station controller is meant the ABIS interface; For the CDMA2000 system, the interface between base station and the base station controller is meant the ABIS interface; For WCDMA or TD-SCDMA system, the interface between base station and the base station controller is meant Iub interface.
Step S02 obtains the engineering parameter table from network operation maintenance centre;
The engineering parameter table comprises the several rows data record, and each line data writes down the configuration parameter of a sub-district in unique corresponding objective network.
For GSM or cdma system, each line data record has all comprised following configuration parameter: geographical position, cell ID, affiliated basic station number, adjacent area Broadcast Control Channel number, neighbor base station identification code;
For CDMA2000, WCDMA or the TD-SCDMA system of 3G, each line data record has all comprised following configuration parameter: geographical position, cell ID, affiliated basic station number, adjacent area common signal channel center frequency point, adjacent area scrambler.
Step S03, write down the Serving cell that comprises for each bar measurement report data, in the engineering parameter table, carry out match search, with the geographical position that obtains this Serving cell, affiliated base station, matching process is: search the data line record that has comprised the same cells sign in the engineering parameter table, as the engineering parameter data record of this Serving cell correspondence.
Step S04, write down each adjacent area that comprises for each bar measurement report data, in the engineering parameter table, carry out match search, with the geographical position that obtains this sub-district, affiliated base station, matching process is: for gsm system or cdma system, search the data record row that comprises identical broadcasts control channel number and base station identity code in the engineering parameter table; For CDMA2000, WCDMA or TD-SCDMA system, in the engineering parameter table, search the data record row that comprises identical common signal channel center frequency point and scrambler.
Step S05, write down each adjacent area that comprises for each bar measurement report data, when the record line number order of the coupling that searches greater than 1 the time, geographical position according to the sub-district in the geographical position of the Serving cell of this measurement report data record and many data record in the engineering parameter table, matching, distance between difference calculation services sub-district and the data record respective cell, that data line record of a selected distance value minimum then is as the engineering parameter data record of this adjacent area correspondence.
Step S06, each adjacent area that comprises for each bar measurement report data record, according to the geographical position of the Serving cell of this measurement report data record and the geographical position of this adjacent area, the distance between calculation services sub-district and the adjacent area, when this distance greater than f dWhen rice, then in thereafter step S07, the sub-district described in the S08, get rid of this adjacent area; f dBe empirical coefficient, general value is a 2000-3500 rice, is preferably 2640 meters.
Step S07, for each bar measurement report data record, judging whether each sub-district has with this sub-district belongs to the administration of same base station is documented in other sub-district in this measurement report data record equally, if do not have other sub-district and this sub-district when belonging to the administration of same base station, this sub-district belongs to the 1st class sub-district; When there being and only having other sub-district and this sub-district when belonging to the administration of same base station in this measurement report, this sub-district belongs to the 2nd class sub-district; When two other sub-district and this sub-districts when belonging to the administration of same base station are arranged in this measurement report, this sub-district belongs to the 3rd class sub-district.
Step S08, utilize approximation method to obtain the general geographic location of target terminal, specifically comprise: to each the bar measurement report data record in this group measurement reporting, its all subdistrict positions that comprise are carried out the statistical average of longitude and latitude respectively, perhaps the received signal level data of sub-district are carried out the weighting statistical average of longitude and latitude as weight coefficient respectively to its all subdistrict positions that comprise, then the one group of geographic position data that calculates is carried out the statistical average of longitude data and latitude data respectively.In the above-mentioned steps, for CDMA2000, WCDMA and TD-SCDMA system, the received signal level of any one sub-district is P-CCPCH received signal code power+110 of this sub-district.
Step S09 to each sub-district that each the bar measurement report in these group measurement reporting data comprises, calculates the distance between the general geographic location of the geographical position of this sub-district and target terminal.
Step S10, to each the bar measurement report in these group measurement reporting data, calculate the distance between each sub-district, calculate the arithmetic average statistical expection and the mean square deviation of distance between each sub-district, calculate the arithmetic average statistical expection and the mean square deviation of the corresponding incoming level in each sub-district, according to the number of above-mentioned statistical index and sub-district, calculate the position error estimated value of approximation method according to following formula (1), unit is a rice:
s = 10 f s &times; E ( d ) f smd &times; D ( d ) f sdd &times; E ( R ) f smr &times; D ( R ) f sdr &times; N f sn - - - ( 1 )
In the formula (1), E (d) represents the arithmetic average statistical expection of distance between each sub-district, and d represents distance between each sub-district; D (d) represents the arithmetic average statistics mean square deviation of distance between each sub-district; E (R) represents the arithmetic average statistical expection of the corresponding incoming level in each sub-district, and R represents the corresponding incoming level in each sub-district; D (R) represents the arithmetic average statistics mean square deviation of the corresponding incoming level in each sub-district; N represents the number of sub-district; f s, f Smd, f Sdd, f Smd, f Smr, f SdrAnd f SnAll be empirical coefficient, have following value respectively:
f s=0.8±0.2、f smd=0.8±0.2、f sdd=0.1±0.025、f smr=-0.1±0.025、f sdr=-0.3±0.075、f sn=-0.75±0.18。
Step S11, to each the bar measurement report in these group measurement reporting data, affiliated base station title according to its sub-district that comprises, all microzonations are divided into several groups, it is all inequality that sub-district in each group all belongs to the base station title that same base station is administered and each group is corresponding, according to incoming level, type, the antenna directional angle of all sub-districts in each group, calculate an angle value of affiliated base station, this group sub-district correspondence;
Wherein, to each group sub-district that each the bar measurement report in these group measurement reporting data comprises, when this group sub-district belonged to the 1st class sub-district, this group only comprised a sub-district, and the antenna directional angle of this sub-district is the angle value of affiliated base station, this group sub-district correspondence.
When this group sub-district belongs to the 2nd class sub-district, this group comprises two sub-districts, if these two sub-districts are respectively sub-district 1 and sub-district 2, when the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 during, adopt following formula (2) to calculate the angle value of base station correspondence under this group sub-district less than 180 °:
Φ=(θ 12)/2+s(θ 12)×f Φ×(ε 12) (2)
Wherein, θ 1Represent the antenna directional angle of sub-district 1; θ 2Represent the antenna directional angle of sub-district 2; ε 1Represent the received signal level of sub-district 1; ε 2Represent the received signal level of sub-district 2; f ΦBe empirical coefficient, its value is 3.35 ± 0.85.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 during, adopt following formula (3) to calculate the angle value of base station correspondence under this group sub-district greater than 180 °:
Φ=(360+θ 12)/2+s(θ 12)×f Φ×(ε 21) (3)
Wherein, the connotation of the variable of the same name in the connotation of each variable and the formula (2) is identical.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 equals 180 °, then should choose the angle value of the antenna directional angle of a bigger sub-district of corresponding received power as base station correspondence under this group sub-district.
When this group sub-district belongs to the 3rd class sub-district, if this group be numbered 1, group 1 comprises 3 sub-districts, establishes these 3 sub-districts and is respectively sub-district 3, sub-district 4 and sub-district 5, and the incoming level of sub-district 3 is greater than the incoming level of sub-district 4, and the incoming level of sub-district 4 is greater than the incoming level of sub-district 5; Then sub-district 3 and sub-district 4 are considered as the 2nd class sub-district, and are combined into an interim group, that establishes this interim group is numbered 2, calculates the angle value of group 2 correspondences; To organize 2 then and be considered as an interim sub-district, establishing this interim sub-district is sub-district 6, and the incoming level of sub-district 6 is the mean value of the incoming level of sub-district 3 and sub-district 4, and the antenna directional angle of sub-district 6 is the angle value of group 2 correspondences.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 during, adopt following formula (4) to calculate the angle value of base station correspondence under the sub-district of group 1 less than 180 °:
Φ=(2×θ 43)/3+4s(θ 34)×f Φ×(ε 34)/3 (4)
Wherein, θ 3Represent the antenna directional angle of sub-district 5; θ 4Represent the antenna directional angle of sub-district 6; ε 3Represent the received signal level of sub-district 5; ε 4Represent the received signal level of sub-district 6.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 during, can adopt following formula (5) to calculate the angle value of base station correspondence under the sub-district of group 1 greater than 180 °:
&Phi; = [ 2 &times; &theta; 4 + ( &theta; 3 + 360 ) ] / 3 - 4 &times; f &Phi; &times; ( &epsiv; 4 - &epsiv; 3 ) / 3 , &theta; 4 > &theta; 3 [ 2 &times; ( &theta; 4 + 360 ) + &theta; 3 ] / 3 + 4 &times; f &Phi; &times; ( &epsiv; 4 - &epsiv; 3 ) / 3 , &theta; 3 > &theta; 4 - - - ( 5 )
Wherein, the connotation of the variable of the same name in the connotation of each variable and the formula (2) is identical.
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 equals 180 °, then should choose the angle value of the antenna directional angle of a bigger sub-district of corresponding received power as base station correspondence under the sub-district of group 1.
Wherein, in formula (2), formula (3) and formula (4), s (θ 12) be sign function, its implication can be represented with following formula (6):
s ( &alpha; ) = 1 , &alpha; > 0 0 , &alpha; = 0 - 1 , &alpha; < 0 - - - ( 6 )
Step S12, to each the bar measurement report in these group measurement reporting data, according to the distance between the general geographic location of the general geographic location of the received signal level of its each sub-district that comprises, target terminal, each sub-district that it comprises and target terminal, the position error of approximation method, some groups, the corresponding angle value of each group of division, carry out iterative computation, obtain one group of terminal geographic position data.
The method of above-mentioned iterative computation can be represented with Fig. 2, wherein, if the E point is the general geographic location of target terminal, A, B, 4 difference of C, D are represented the geographical position of the affiliated base station, sub-district in pairing each group of this measurement report successively, a, b, four rays of c, d are represented the angle value of the above-mentioned base station correspondence that obtains respectively successively in step S11, the angular range that the dotted line representative is certain; During iterative computation, being the center of circle with the A point at first, is radius with AE, and the E point is rotated to an angle to the direction near ray a, arrives the F point; Then be the center of circle, be radius with BF, the F point is rotated to an angle to the direction near ray b, arrive the G point with the B point; Then be the center of circle, be radius with CG, the G point is rotated to an angle to the direction near ray c, arrive the H point with the C point; Then be the center of circle, be radius with DH, the G point is rotated to an angle to the direction near ray d, arrive the I point with the D point; Then, be the center of circle with A again, be radius rotation with AI; Before once rotating arbitrarily, if point to be rotated, for example, E point, F point, G point or H point, the orientation of pivot is in beyond the angular range of two dotted lines that pivot next time connects relatively next time, is not that the center of circle is rotated with this center then, and skips, whether the geographical position with next base station is the center of circle, rejudge and rotate.
Total number of revolutions is between 10 to 30 times, and the amplitude of each rotation satisfies following formula (7):
Figure BDA0000040574180000121
Wherein,
Figure BDA0000040574180000122
Be the center of circle to be rotated (under any one group of sub-district that comprises the geographical position of base station) with point (the interim result of iterative computation, preliminary examination value are the general geographic location of target terminal) to be rotated between distance; R iIt is the assembly average of organizing the received signal level of the sub-district that comprises; S is the position error of the approximation method that obtains in step S10; N is the number of organizing in this measurement report; Clip is the azimuth in point the to be rotated center of circle for the treatment of rotation mutually and the angle between the rolling target ray; f v, f Vd, f Vr, f VsAnd f VnAll be empirical coefficient, they generally have following value:
f v=-6.7±1.67、f vd=-1.5±0.37、f vr=3.35±0.8、f vs=1.85±0.46、f vn=-1±0.25。
Step S13 at one group of terminal geographic position data that step S12 obtains, carries out the arithmetic statistical average of longitude data and latitude data respectively, finally obtains the geographical position of target terminal.
It should be noted last that, above embodiment is only unrestricted in order to technical scheme of the present invention to be described, although the present invention is had been described in detail with reference to preferred embodiment, those of ordinary skill in the art is to be understood that, can make amendment or be equal to replacement technical scheme of the present invention, and not breaking away from the spirit and scope of technical solution of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.

Claims (10)

1. the method for localisation of mobile terminals in the cellular network is characterized in that, may further comprise the steps:
Step 1, utilize the cellular network interface snooping equipment to intercept interface between base station and the base station controller, to obtain group measurement reporting data of corresponding target terminal;
Step 2, obtain the engineering parameter table from network operation maintenance centre;
Step 3, each sub-district that each the bar measurement report in these group measurement reporting data is comprised, in the engineering parameter table, carry out match search, with the geographical position that obtains this sub-district, affiliated base station, judge whether this sub-district corresponding parameters is abandoned according to the geographical position of this sub-district;
Step 4, administer and determined the type of this sub-district according to base station under this sub-district by the number of cells of this measurement report record;
Step 5 utilizes approximation method to calculate the general geographic location at target terminal place;
Step 6 is to each sub-district that each the bar measurement report in these group measurement reporting data comprises, the distance between the geographical position of calculation plot and the general geographic location of target terminal;
Step 7, to each the bar measurement report in these group measurement reporting data, calculate the distance between each sub-district, the arithmetic average statistical expection and the mean square deviation of distance between the calculation plot, calculate the arithmetic average statistical expection and the mean square deviation of the corresponding incoming level in each sub-district, calculate the position error of described general geographic location according to the number of above-mentioned statistical index and sub-district;
Step 8, to each the bar measurement report in these group measurement reporting data, affiliated base station title according to its sub-district that comprises, all microzonations are divided into several groups, it is all inequality that sub-district in each group all belongs to the base station title that same base station is administered and each group is corresponding, according to incoming level, type, the antenna directional angle of all sub-districts in each group, calculate an angle value of affiliated base station, this group sub-district correspondence;
Step 9, to each the bar measurement report in these group measurement reporting data, according to the distance between the general geographic location of the general geographic location of the received signal level of its each sub-district that comprises, target terminal, each sub-district that it comprises and target terminal, the position error of approximation method, some groups, the corresponding angle value of each group of division, carry out iterative computation, obtain one group of target terminal geographic position data;
Step 10 at one group of terminal geographic position data that step 9 obtains, is carried out the arithmetic average statistical average of longitude data and latitude data respectively, obtains final objective terminal geographic position.
2. the method for localisation of mobile terminals is characterized in that in the cellular network according to claim 1, and described step 3 specifically may further comprise the steps:
Step a, the Serving cell that comprises for each bar measurement report data record, in the engineering parameter table, carry out match search, matching process with the geographical position that obtains this Serving cell, affiliated base station is: search the data line record that has comprised the same cells sign in the engineering parameter table, as the engineering parameter data record of this Serving cell correspondence;
Step b, each adjacent area that comprises for each bar measurement report data record, in the engineering parameter table, carry out match search, matching process with the geographical position that obtains this sub-district, affiliated base station is: for gsm system or cdma system, search the data record row that comprises identical broadcasts control channel number and base station identity code in the engineering parameter table; For CDMA2000, WCDMA or TD-SCDMA system, in the engineering parameter table, search the data record row that comprises identical common signal channel center frequency point and scrambler;
Step c, each adjacent area that comprises for each bar measurement report data record, when the record line number order of the coupling that searches greater than 1 the time, geographical position according to the sub-district in the geographical position of the Serving cell of this measurement report data record and many data record in the engineering parameter table, matching, distance between difference calculation services sub-district and the data record respective cell, that data line record of a selected distance value minimum then is as the engineering parameter data record of this adjacent area correspondence;
Steps d, each adjacent area that comprises for each bar measurement report data record, according to the geographical position of the Serving cell of this measurement report data record and the geographical position of this adjacent area, the distance between calculation services sub-district and the adjacent area, when this distance greater than f dWhen rice, then in thereafter step S07, the sub-district described in the S08, get rid of this adjacent area; f dBe empirical coefficient, value is a 2000-3500 rice.
3. the method for localisation of mobile terminals is characterized in that in the cellular network according to claim 2, described f dBe preferably 2640 meters.
4. according to the method for localisation of mobile terminals in claim 2 or the 3 described cellular networks, it is characterized in that, described step 4 specifically comprises: write down each adjacent area that comprises for each bar measurement report data, when the record line number order of the coupling that searches greater than 1 the time, geographical position according to the sub-district in the geographical position of the Serving cell of this measurement report data record and many data record in the engineering parameter table, matching, distance between difference calculation services sub-district and the data record respective cell, that data line record of a selected distance value minimum then is as the engineering parameter data record of this adjacent area correspondence.
5. the method for localisation of mobile terminals in the cellular network according to claim 4, it is characterized in that, the approximation method of described step 5 comprises: to each the bar measurement report data record in this group measurement reporting, its all subdistrict positions that comprise are carried out the statistical average of longitude and latitude respectively, perhaps the received signal level data of sub-district are carried out the weighting statistical average of longitude and latitude as weight coefficient respectively to its all subdistrict positions that comprise, then the one group of geographic position data that calculates is carried out the statistical average of longitude data and latitude data respectively; For CDMA2000, WCDMA and TD-SCDMA system, the received signal level of any one sub-district is P-CCPCH received signal code power+110 of this sub-district.
6. the method for localisation of mobile terminals is characterized in that in the cellular network according to claim 5, in the described step 7, adopts following formula (1) to calculate the position error of approximation method:
s = 10 f s &times; E ( d ) f smd &times; D ( d ) f sdd &times; E ( R ) f smr &times; D ( R ) f sdr &times; N f sn - - - ( 1 )
In the formula (1), E (d) represents the arithmetic average statistical expection of distance between each sub-district, and d represents distance between each sub-district; D (d) represents the arithmetic average statistics mean square deviation of distance between each sub-district; E (R) represents the arithmetic average statistical expection of the corresponding incoming level in each sub-district, and R represents the corresponding incoming level in each sub-district; D (R) represents the arithmetic average statistics mean square deviation of the corresponding incoming level in each sub-district; N represents the number of sub-district; f s, f Smd, f Sdd, f Smd, f Smr, f SdrAnd f SnAll be empirical coefficient, have following value respectively:
f s=0.8±0.2、f smd=0.8±0.2、f sdd=0.1±0.025、f smr=-0.1±0.025、f sdr=-0.3±0.075、f sn=-0.75±0.18。
7. the method for localisation of mobile terminals is characterized in that in the cellular network according to claim 6, and described step 8 specifically may further comprise the steps:
To each group sub-district that each the bar measurement report in these group measurement reporting data comprises, when this group sub-district belonged to the 1st class sub-district, this group only comprised a sub-district, and the antenna directional angle of this sub-district is the angle value of affiliated base station, this group sub-district correspondence;
When this group sub-district belongs to the 2nd class sub-district, this group comprises two sub-districts, if these two sub-districts are respectively sub-district 1 and sub-district 2, when the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 during, adopt following formula (2) to calculate the angle value of base station correspondence under this group sub-district less than 180 °:
Φ=(θ 12)/2+s(θ 12)×f Φ×(ε 12) (2)
Wherein, θ 1Represent the antenna directional angle of sub-district 1; θ 2Represent the antenna directional angle of sub-district 2; ε 1Represent the received signal level of sub-district 1; ε 2Represent the received signal level of sub-district 2; f ΦBe empirical coefficient, its value is 3.35 ± 0.85;
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 during, adopt following formula (3) to calculate the angle value of base station correspondence under this group sub-district greater than 180 °:
Φ=(360+θ 12)/2+s(θ 12)×f Φ×(ε 21)(3)
Wherein, the connotation of the variable of the same name in the connotation of each variable and the formula (2) is identical;
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 1 and sub-district 2 equals 180 °, then should choose the angle value of the antenna directional angle of a bigger sub-district of corresponding received power as base station correspondence under this group sub-district;
When this group sub-district belongs to the 3rd class sub-district, if this group be numbered 1, group 1 comprises 3 sub-districts, establishes these 3 sub-districts and is respectively sub-district 3, sub-district 4 and sub-district 5, and the incoming level of sub-district 3 is greater than the incoming level of sub-district 4, and the incoming level of sub-district 4 is greater than the incoming level of sub-district 5; Then sub-district 3 and sub-district 4 are considered as the 2nd class sub-district, and are combined into an interim group, that establishes this interim group is numbered 2, calculates the angle value of group 2 correspondences; To organize 2 then and be considered as an interim sub-district, establishing this interim sub-district is sub-district 6, and the incoming level of sub-district 6 is the mean value of the incoming level of sub-district 3 and sub-district 4, and the antenna directional angle of sub-district 6 is the angle value of group 2 correspondences;
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 during, adopt following formula (4) to calculate the angle value of base station correspondence under the sub-district of group 1 less than 180 °:
Φ=(2×θ 43)/3+4s(θ 34)×f Φ×(ε 34)/3(4)
Wherein, θ 3Represent the antenna directional angle of sub-district 5; θ 4Represent the antenna directional angle of sub-district 6; ε 3Represent the received signal level of sub-district 5; ε 4Represent the received signal level of sub-district 6;
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 during, can adopt following formula (5) to calculate the angle value of base station correspondence under the sub-district of group 1 greater than 180 °:
&Phi; = [ 2 &times; &theta; 4 + ( &theta; 3 + 360 ) ] / 3 - 4 &times; f &Phi; &times; ( &epsiv; 4 - &epsiv; 3 ) / 3 , &theta; 4 > &theta; 3 [ 2 &times; ( &theta; 4 + 360 ) + &theta; 3 ] / 3 + 4 &times; f &Phi; &times; ( &epsiv; 4 - &epsiv; 3 ) / 3 , &theta; 3 > &theta; 4 - - - ( 5 )
Wherein, the connotation of the variable of the same name in the connotation of each variable and the formula (2) is identical;
When the absolute value of the difference of the antenna directional angle of the antenna directional angle of sub-district 5 and sub-district 6 equals 180 °, then should choose the angle value of the antenna directional angle of a bigger sub-district of corresponding received power as base station correspondence under the sub-district of group 1;
Wherein, in formula (2), formula (3) and formula (4), s (θ 12) be sign function, its implication can be represented with following formula (6):
s ( &alpha; ) = 1 , &alpha; > 0 0 , &alpha; = 0 - 1 , &alpha; < 0 - - - ( 6 )
8. the method for localisation of mobile terminals is characterized in that in the cellular network according to claim 7, and described step 9 specifically may further comprise the steps:
If the E point is the general geographic location of target terminal, A, B, 4 difference of C, D are represented the geographical position of the affiliated base station, sub-district in pairing each group of this measurement report successively, a, b, four rays of c, d are represented the angle value of the above-mentioned base station correspondence that obtains respectively successively in step S11, the angular range that the dotted line representative is certain; During iterative computation, being the center of circle with the A point at first, is radius with AE, and the E point is rotated to an angle to the direction near ray a, arrives the F point; Then be the center of circle, be radius with BF, the F point is rotated to an angle to the direction near ray b, arrive the G point with the B point; Then be the center of circle, be radius with CG, the G point is rotated to an angle to the direction near ray c, arrive the H point with the C point; Then be the center of circle, be radius with DH, the G point is rotated to an angle to the direction near ray d, arrive the I point with the D point; Then, be the center of circle with A again, be radius rotation with AI; Before once rotating arbitrarily, if point to be rotated, for example, E point, F point, G point or H point, the orientation of pivot is in beyond the angular range of two dotted lines that pivot next time connects relatively next time, is not that the center of circle is rotated with this center then, and skips, whether the geographical position with next base station is the center of circle, rejudge and rotate;
Total number of revolutions is between 10 to 30 times, and the amplitude of each rotation satisfies following formula (7):
Wherein,
Figure FDA0000040574170000053
Be the center of circle to be rotated (under any one group of sub-district that comprises the geographical position of base station) with point (the interim result of iterative computation, preliminary examination value are the general geographic location of target terminal) to be rotated between distance; R iIt is the assembly average of organizing the received signal level of the sub-district that comprises; S is the position error of the approximation method that obtains in step S10; N is the number of organizing in this measurement report; Clip is the azimuth in point the to be rotated center of circle for the treatment of rotation mutually and the angle between the rolling target ray; f v, f Vd, f Vr, f VsAnd f VnAll be empirical coefficient, they generally have following value:
f v=-6.7±1.67、f vd=-1.5±0.37、f vr=3.35±0.8、f vs=1.85±0.46、f vn=-1±0.25。
9. the method for localisation of mobile terminals is characterized in that in the cellular network according to claim 6, described f s=0.8, f Smd=0.8, f Sdd=0.1, f Smr=-0.1, f Sdr=-0.3, f Sn=-0.75.
10. the method for localisation of mobile terminals is characterized in that in the cellular network according to claim 8, described f v=-6.7, f Vd=-1.5, f Vr=3.35, f Vs=1.85, f Vn=-1.
CN 201010606178 2010-12-24 2010-12-24 Method for positioning mobile terminal in cellular network Expired - Fee Related CN102056224B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010606178 CN102056224B (en) 2010-12-24 2010-12-24 Method for positioning mobile terminal in cellular network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010606178 CN102056224B (en) 2010-12-24 2010-12-24 Method for positioning mobile terminal in cellular network

Publications (2)

Publication Number Publication Date
CN102056224A true CN102056224A (en) 2011-05-11
CN102056224B CN102056224B (en) 2013-03-20

Family

ID=43960021

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010606178 Expired - Fee Related CN102056224B (en) 2010-12-24 2010-12-24 Method for positioning mobile terminal in cellular network

Country Status (1)

Country Link
CN (1) CN102056224B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012119427A1 (en) * 2011-08-30 2012-09-13 华为技术有限公司 Cellular network location method, service platform, location mobile station and system
CN103024673A (en) * 2013-01-14 2013-04-03 北京中创信测科技股份有限公司 Terminal positioning method and device for GSM (global system for mobile communication) networks
CN103747525B (en) * 2014-01-26 2017-10-24 林汉忠 A kind of compound alignment system based on mobile communications network
CN108375363A (en) * 2017-12-05 2018-08-07 中国移动通信集团福建有限公司 Antenna bearingt angular deflection check method, device, equipment and medium
CN108882292A (en) * 2017-05-15 2018-11-23 中兴通讯股份有限公司 Measurement report localization method and device
CN110012416A (en) * 2017-12-29 2019-07-12 中国移动通信集团浙江有限公司 A kind of localization method and device of user terminal
CN114363805A (en) * 2021-12-29 2022-04-15 北京红山信息科技研究院有限公司 Improved triangulation location method for correcting back points
CN114363805B (en) * 2021-12-29 2024-04-30 北京红山信息科技研究院有限公司 Triangular positioning improvement method for correcting back point

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060240844A1 (en) * 2002-11-18 2006-10-26 Polaris Wireless Inc. Computationally-Efficient Estimation of the Location of a Wireless Terminal Based on Pattern Matching
CN101248680A (en) * 2005-07-05 2008-08-20 凯锐恩Iq公司 Rule based data collection and management in a wireless communications network

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060240844A1 (en) * 2002-11-18 2006-10-26 Polaris Wireless Inc. Computationally-Efficient Estimation of the Location of a Wireless Terminal Based on Pattern Matching
CN101248680A (en) * 2005-07-05 2008-08-20 凯锐恩Iq公司 Rule based data collection and management in a wireless communications network

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012119427A1 (en) * 2011-08-30 2012-09-13 华为技术有限公司 Cellular network location method, service platform, location mobile station and system
CN103024673A (en) * 2013-01-14 2013-04-03 北京中创信测科技股份有限公司 Terminal positioning method and device for GSM (global system for mobile communication) networks
CN103747525B (en) * 2014-01-26 2017-10-24 林汉忠 A kind of compound alignment system based on mobile communications network
CN108882292A (en) * 2017-05-15 2018-11-23 中兴通讯股份有限公司 Measurement report localization method and device
CN108375363A (en) * 2017-12-05 2018-08-07 中国移动通信集团福建有限公司 Antenna bearingt angular deflection check method, device, equipment and medium
CN108375363B (en) * 2017-12-05 2021-06-04 中国移动通信集团福建有限公司 Antenna azimuth deflection checking method, device, equipment and medium
CN110012416A (en) * 2017-12-29 2019-07-12 中国移动通信集团浙江有限公司 A kind of localization method and device of user terminal
CN110012416B (en) * 2017-12-29 2021-03-02 中国移动通信集团浙江有限公司 User terminal positioning method and device
CN114363805A (en) * 2021-12-29 2022-04-15 北京红山信息科技研究院有限公司 Improved triangulation location method for correcting back points
CN114363805B (en) * 2021-12-29 2024-04-30 北京红山信息科技研究院有限公司 Triangular positioning improvement method for correcting back point

Also Published As

Publication number Publication date
CN102056224B (en) 2013-03-20

Similar Documents

Publication Publication Date Title
CN102056224B (en) Method for positioning mobile terminal in cellular network
CN1320793C (en) Device and method for judging position of mobile communication terminal
CN1158890C (en) Method and system for determining position of mobile radio terminals
CA2660910C (en) Embedded wireless location validation benchmarking systems and methods
JP4855340B2 (en) Area map construction system that creates area maps based on communication quality information
US8233912B2 (en) Mobile station localization apparatus, systems, and methods
US8483714B2 (en) Method and apparatus for identification of local beacon systems
CN101631349B (en) Method, device and wireless operation maintenance center for positioning terminal
CN1909731B (en) Method for collecting position data of base station
CN1852591B (en) Method for positioning mobile terminal
CN1663299A (en) Method and system for determining the speed and distance of a mobile unit
WO1998035524A2 (en) Method and system for locating a mobile station within a mobile telecommunications network
JP2001512940A (en) System and method using elliptical search area coverage for mobile terminal location
CN103024895A (en) Positioning method for positioning by combining GPS (global positioning system) with base station
CN102740350A (en) Indoor positioning method, device and system
CN103874009A (en) Positioning infomration indicating method, signal sending device, access point, terminal and system
EP0885543B1 (en) Location method and system
CN102014429A (en) Method for fast positioning of mobile terminal in cellular network
CN101635985A (en) Positioning method, positioning device and recording medium
CN104168646A (en) Positioning method and apparatus of wireless access equipment
CN100512514C (en) Positioning method applied in multi carrier frequency system
CN103874029B (en) Message method, fake satellite positioning system and the relevant apparatus of pseudolite systems
WO2011097182A1 (en) Technique for effectively communication location information in a wireless communication service
US20150207860A1 (en) Information distribution system, information distribution server, information distribution method, and program
CN1554546A (en) Vehicle carried communication alarm positioning method and system based on mobile network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130320

Termination date: 20131224