CN103002502A - Positioning method and system in code division multiple access (CDMA) based on measurement report (MR) - Google Patents

Abstract

The invention discloses a positioning method and a system in CDMA based on an MR. The positioning method comprises obtaining a cell identifier CELL ID in report information in a CDMA product, and determining a cell coverage range; obtaining time of arrival (TOA) of signals of a mobile terminal report branch in MT tracking, and calculating the distance between the mobile terminal and a base station according to the TOA; and performing positioning analysis according to the determined cell coverage range and the distance between the mobile terminal and the base station. The positioning method is fast in calculation response speed, and the accuracy is higher than field intensity positioning accuracy and positioning accuracy on the basis of characteristic matching.

Description

A kind of method and system that in CDMA, realize the location based on measurement report

Technical field

The present invention relates to communication technical field, relate in particular to a kind of method and system that in CDMA, realize the location based on measurement report.

Background technology

MR(Measurement Report) be called measurement report, the cycle that is defined in the agreement in the communication process sends or Event triggered reports, and has comprised the radio environment information in online all user geographical position of living in.Doing the network optimization with the MR data can be near the data reflection of true wireless environment, the data of the perception of being close to the users most.Compare with traditional drive test DT, the positional information of measurement report is accurate not as drive test behind the location, but because measurement report can be got the information of all user's random time sections, its magnanimity statistical nature is that drive test is incomparable.

Network optimization analysis based on MR is the real-time status of utilizing the existing network cellphone subscriber, based on fast real magnanimity mobile phone measurement report, directly gather the front reverse link data of cellphone subscriber, under actual traffic model, be optimized, can find in time that covering problem, adjacent area polygamy leak and join problem, intuitively understand all sidedly traffic density, front reverse quality and disturbance regime.Measurement report is done the location, position, is all terminal locations to be got ready at map play up, and embodies the fine or not degree of evaluation index with the variation of different colours.

Based on the situation in the code division multiple access CDMA technology, existing mobile terminal locating method is a lot, but does not also have effectively realization based on the location implementing method of MR.

Summary of the invention

The present invention is based on the deficiency of existing location technology, a kind of method and system realize the location based on measurement report in CDMA are provided, can more high efficiencyly satisfy terminal positioning.

The invention provides a kind of method that in CDMA, realizes the location based on measurement report, comprise the steps:

Obtain the cell number CELLID in the reporting information in the CDMA product, determine cell coverage area;

In measurement report MR follows the tracks of, obtain the time of arrival (toa) TOA of mobile terminal reporting branch, calculate the distance of portable terminal and base station according to TOA;

Position analysis according to the cell coverage area of determining and the distance of portable terminal and base station.

The cell coverage area that described basis is determined and the distance of portable terminal and base station position in the analysis and comprise:

Judge the way situation that reports branch, position analysis according to the way situation that reports branch of judging.

The cell coverage area that described basis is determined and the distance of portable terminal and base station position to analyze and comprise:

When judging when only having one tunnel branch to report TOA, in conjunction with TOA and sector directions information, the method for positioning analyzing that use direction of arrival degree AOA combines with TOA analyzes the position of terminal.

The cell coverage area that described basis is determined and the distance of portable terminal and base station position to analyze and comprise:

When judgement has the TOA that two-way branch reports, use TOA positioning analysis method, analyze the position of terminal by two circle intersection algorithms.

The cell coverage area that described basis is determined and the distance of portable terminal and base station position to analyze and comprise:

When judgement had branch more than three tunnel to report TOA, the signal arrival time difference TDOA by the TOA of branch and reference radio channel TOA difference structure two-way branch used hyperbolic method to analyze the position of terminal.

Accordingly, the embodiment of the invention also provides a kind of system that realizes the location in CDMA based on measurement report, comprising:

The cell number acquisition module, the cell number CELLID in the reporting information that obtains the CDMA product determines cell coverage area;

The time of arrival (toa) acquisition module is used for following the tracks of the time of arrival (toa) TOA that obtains mobile terminal reporting branch at measurement report MR, calculates the distance of portable terminal and base station according to TOA;

The positioning analysis module is used for positioning analysis according to the cell coverage area of determining and the distance of portable terminal and base station.

Comprise in the described positioning analysis module:

Report the Bifurcation Analysis unit, be used for judging the way situation that reports branch, position analysis according to the way situation that reports branch of judging.

Described locating module comprises:

The first branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains one the tunnel analyzes when reporting TOA, wakes the first analyzing and positioning processing unit up and processes;

The first analyzing and positioning processing unit is used in conjunction with TOA and sector directions information, and the method for positioning analyzing that use direction of arrival degree AOA combines with TOA analyzes the position of terminal.

Described locating module comprises:

The second branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains is the two-way analysis when reporting TOA, wakes the second analyzing and positioning processing unit up and processes;

The second analyzing and positioning processing unit is used for using TOA positioning analysis method, analyzes the position of terminal by two circle intersection algorithms.

Described locating module comprises:

The 3rd branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains is branch more than three tunnel when reporting TOA, wakes the 3rd analyzing and positioning processing unit up and processes;

The 3rd analyzing and positioning processing unit is used for using TOA positioning analysis method, analyzes the position of terminal by two circle intersection algorithms.

By implementing the embodiment of the invention, the present invention mainly determines the terminal approximate range based on terminal to report reference cell or adjacent cell ID, based on reference pilot and activate the distance that time delay that level branch reports is determined terminal and known location residential quarter.According to the different number of branches of reporting message, select in conjunction with AOA, TOA, the signal parameter localization methods such as TDOA, thus finally determine terminal location.This localization method calculated response speed is fast, can obtain 50 meters to 300 meters positioning accuracy, and this precision is higher than the field intensity positioning accuracy, is higher than based on the characteristic matching positioning accuracy.

Description of drawings

Fig. 1 for the embodiment of the invention provide according to TOA location estimation MS position view;

Fig. 2 for the embodiment of the invention provide according to TDOA location estimation MS position view;

Fig. 3 for the embodiment of the invention provide according to AOA location estimation MS position view;

The method flow diagram of in CDMA, realizing the location based on measurement report that Fig. 4 provides for the embodiment of the invention;

The method particular flow sheet of in CDMA, realizing the location based on measurement report that Fig. 5 provides for the embodiment of the invention;

Two circles that Fig. 6 provides for the embodiment of the invention comprise the exemplary plot of situation;

Two round exemplary plot from situation that Fig. 7 provides for the embodiment of the invention;

Two exemplary plot of justifying tangent situations that Fig. 8 provides for the embodiment of the invention;

Two circles that Fig. 9 provides for the embodiment of the invention intersect the exemplary plot of situation;

Two circles that Figure 10 provides for the embodiment of the invention intersect another exemplary plot of situation;

Effective branch that Figure 11 provides for the embodiment of the invention is more than or equal to the schematic diagram of three situations;

The system configuration schematic diagram of in CDMA, realizing the location based on measurement report that Figure 12 provides for the embodiment of the invention;

The positioning analysis modular structure schematic diagram that Figure 13 provides for the embodiment of the invention.

Embodiment

Describe the wireless terminal schematic diagram that shows based on miniature literal in detail below in conjunction with embodiments of the invention.

Wireless location is exactly the real time information that self-position is provided to the mobile subscriber in the valid wireless coverage.All be the changing features of utilizing wireless signal in essence, such as propagation time, time difference, signal strength, phase place or incidence angle etc., by its certain geometrical principle calculated target positions.

Propagation delay time estimated result according to signal positions portable terminal, needs at least three base stations.The angle of utilizing signal to arrive positions, and needs at least two base stations.If utilize simultaneously propagation delay time and angle to position, then only need a base station to get final product.Since be subjected to the mobile telecommunication channel noise, the impact of the factors such as multipath transmisstion or non line of sight, the cellular radio navigation system is difficult to reach higher positioning accuracy, and the location coverage also is subject to the restriction of cell mobile communication systems covering range of field intensity.In the situation of Multipath Transmission, it is slower to utilize the base station that portable terminal is carried out the variation that location estimation comprised the time-delay of sense and time, and decline only affects the intensity of signal.That is to say that angle and time delays parameter are metastable.

A. play source honeycomb district COO (Cell Of Origin) location technology

The source honeycomb district location technology is to determine user's position according to the residing Cell Identity CI ID of travelling carriage.Travelling carriage just has in the database of system corresponding residential quarter ID number after the current area registration.As long as the center (position in local map) that system can arrange this cell base station and the covering radius of residential quarter are broadcast to all travelling carriages in the cell range, these travelling carriages just can know where oneself be in, Query Database can obtain positional information.The advantage of this targeting scheme is to need not network and mobile phone are made amendment, and the response time is short.But its positioning accuracy depends on the radius of residential quarter.

B. arrival time TOA (Time Of Arrival) location technology

Fig. 1 is according to TOA location estimation MS position, if mobile station MS and base station BS are in the visual scope, then therebetween apart from d _iFor: d _i=c Δ t _iWherein c is the light velocity, Δ t _iFor radio signal at MS and BS _iBetween propagation delay.

When mobile station MS is known with the propagation delay of three BS at least, just can calculate the position (seeing Fig. 1) of MS.Be not in the direct-view scope when between mobile station MS and BS, then can produce error to the estimation of MS position, namely the circular arc among the figure can not intersect at a point.So can adopt the distance between MS position and TOA circular arc is the method for mean square error minimum.If the TOA measurement result more than 3 is arranged, then precision can improve.

C. arrival time difference TDOA (Time Difference OfArrival) location technology

Fig. 2 is according to TDOA location estimation MS position, and this method is based on the time difference of the signal of receiving different base station BS.If MS receives neighbor bss ₁And BS ₂Time difference of signal be Δ t, if at this moment MS and BS are in the visual range, then the position of MS is on a hyperbola (seeing Fig. 2).d ₁-d ₂=c Δ t, wherein d ₁For MS arrives BS ₁Distance, d ₂For MS arrives BS ₂Distance, c is the light velocity, adopts three different base stations can measure two TDOA (step-out time), mobile radio station is positioned on the hyp intersection point that two TDOA determine.

Similar to the TOA method, the localization method of TDOA method can adopt mobile radio station MS to the algorithm of hyperbolic distance mean square error minimum, and condition is to have plural TDOA value can be used for calculating.

One of TDOA method and TOA method comparison advantage are: when calculating TDOA value, the error of calculation be identical to all base stations and itself and be zero, these errors comprise public multidiameter delay and synchronous error.This location technology has reduced the synchronous requirement to the time.Performance is relatively superior under the error environment, receives much concern in the location technology of cellular communication system.

D. angle arrives AOA (Angle Of Arrival) location technology

Fig. 3 is according to AOA location estimation MS position, the angle that this technology arrives according to signal, determine travelling carriage with respect to the angular relationship of base station, as long as measure a travelling carriage apart from the direction of arrival degree of two base stations, just can determine the position (seeing Fig. 3) of travelling carriage.The measurement of direction of arrival degree needs directional antenna, and shortcoming is to need complicated antenna, is prone to the location blind spot.

In the embodiment of the invention based on the localization method of MR.In the reporting information CELLID is arranged in the CDMA product at present, corresponding engineering parameter table can be determined cell coverage area substantially.In MR tracking or CHR tracking, can obtain the TOA that mobile phone reports branch, can calculate mobile phone to the distance of base station according to TOA, can be used for positioning analysis.

When only having one tunnel branch to report TOA, in conjunction with TOA and sector directions information, the positioning mode of using AOA to be combined with TOA provides the position of terminal.When the TOA that has two-way branch to report, use the TOA positioning mode, Xie Shuan circle intersection point is asked terminal location.Optimal solution is got in the AOA cooperation.In view of the TDOA localization method is better than because the TOA localization method, as the TOA that reports during more than or equal to three, by the TDOA of the TOA of branch with reference pilot TOA difference structure two-way branch, the position that the use hyperbolic method can unique definite terminal.

Concrete, Fig. 4 shows the method flow diagram of realizing the location in CDMA based on measurement report in the embodiment of the invention, comprises the steps:

S401: obtain the cell number CELLID in the reporting information in the CDMA product, determine cell coverage area;

S402: in measurement report MR follows the tracks of, obtain the time of arrival (toa) TOA of mobile terminal reporting branch, calculate the distance of portable terminal and base station according to TOA;

S403: position analysis according to the cell coverage area of determining and the distance of portable terminal and base station.

Concrete, Fig. 5 shows the method particular flow sheet of realizing the location in CDMA based on measurement report in the invention process example, comprises the steps:

S501: obtain the cell number CELLID in the reporting information in the CDMA product, determine cell coverage area;

S502: in measurement report MR follows the tracks of, obtain the time of arrival (toa) TOA of mobile terminal reporting branch, calculate the distance of portable terminal and base station according to TOA;

S503: judge the way situation that reports branch;

The follow-up main way situation that reports branch according to judging positions analysis, and it mainly embodies in S504 to S506 to some extent.

S504: when judging when only having one tunnel branch to report TOA, in conjunction with TOA and sector directions information, the method for positioning analyzing that use direction of arrival degree AOA combines with TOA analyzes the position of terminal.

Choose residential quarter corresponding to reference pilot TOA as the reference residential quarter, determine between terminal and reference cell apart from R=cTOA by the TOA value.

If the omnidirectional station utilizes the range information of reference cell position and terminal and reference cell to select a bit as positioning result.

If Directing Station according to sector directions angle, subtended angle and TOA distance, is done a circular arc, terminal location is randomly dispersed on this circular arc.

Reference cell longitude and latitude (Lat ₀, Long ₀), being located in the rectangular coordinate system, the reference cell position is the origin of coordinates (0,0), terminal and reference cell distance are R, the actual coordinate (x, y) of terminal, then: x ²+ y ²=R ²

If the omnidirectional station is got at random and is satisfied x ²+ y ²=R ²(x, the y) of formula is the positioning result that the longitude and latitude form can obtain terminal with Coordinate Conversion.

If Directing Station, establishing the sector directions angle is α, and subtended angle is β, gets at random and satisfies x ²+ y ²=R ²(x, the y) of formula, wherein $x\∈[R\cos (\mathrm{\α}-\frac{1}{2}\mathrm{\β}),R\cos (\mathrm{\α}+\frac{1}{2}\mathrm{\β})]$ And $y\∈[R\sin (\mathrm{\α}-\frac{1}{2}\mathrm{\β}),R\sin (\mathrm{\α}+\frac{1}{2}\mathrm{\β})].$ Be the positioning result that the longitude and latitude form can obtain terminal with Coordinate Conversion.

The processing of occuring simultaneously is got in 2 sectors of same BTS, consist of new logic district.New logic district angle of coverage diminishes, and radius is got 2 radiuses that the sector is less.The deflection of supposing two sectors is respectively: α ₁, α ₂Subtended angle is respectively: β ₁, β ₂, then the deflection of new logic district is:

Subtended angle is:

α wherein ₁+ α ₂∈ (0,2 π].Get the mean value of 2 TOA as new TOA.Location algorithm by an effective TOA positions.

S505: when judgement has the TOA that two-way branch reports, use TOA positioning analysis method, analyze the position of terminal by two circle intersection algorithms.

Get residential quarter corresponding to reference pilot TOA as the reference residential quarter.Find the solution the positioning result of MS by 2 circle intersection algorithms.

According to base station location, and the distance of the TOA between terminal and base station can be determined a circle.When existing two TOA to adhere to different base station separately, can ask terminal location by the intersection point of circle on the geometry.But exist two circles not have intersection point, an intersection point and three kinds of situations of two intersection points, the below will describe respectively.

If d ₁=cTOA ₁, d ₂=cTOA ₂, the longitude and latitude (Lat of two BTS of acquisition ₁, Long ₁), (Lat ₂, Long ₂) convert rectangular coordinate (0,0), (m, 0) to.Wherein m is two BTS distances.

1) do not have intersection point, two kinds of situations arranged, be the TOA measure error and cause, can only get an approximate result:

Fig. 6 shows the situation that two circles comprise, and two circles comprise:

If two base stations are the omnidirectional station, on two base station lines, get the mid point of the nearest line of two circles as positioning result.

（a）d ₁＞d ₂， $x=\frac{m+{\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">d</figure-callout>}}_1+d_2}{2}$

（b）d ₁＜d ₂， $x=\frac{m-d_1-d_2}{2}$

If it is Directing Station that a base station is only arranged, utilize the deflection of this Directing Station, subtended angle and Directing Station TOA position analysis according to an effective TOA;

If two base stations are Directing Station, utilize the separately deflection of Directing Station, subtended angle and Directing Station TOA position computational methods according to an effective TOA, choose at random on two camber lines a bit as positioning result.

It is two round from situation that Fig. 7 shows

If two base stations are the omnidirectional station, then on two base station lines, get the mid point of the nearest line of two circles as positioning result. $x=\frac{m+d_1-d_2}{2}$

If it is Directing Station that a base station is only arranged, utilize the deflection of this Directing Station, subtended angle and Directing Station TOA position analysis according to an effective TOA;

If two base stations are Directing Station, utilize the separately deflection of Directing Station, subtended angle and Directing Station TOA position analysis according to an effective TOA, choose at random on two camber lines a bit as positioning result.

2) Fig. 8 shows and only has an intersection point situation, and directly removing intersection point is the terminal positioning result;

x＝d ₁,y＝0

When 3) having two intersection points, as shown in Fig. 9 and Figure 10:

If the coordinate of two intersection points is respectively (x, y ₁), (x, y ₂), d is then arranged ₁ ²-x ²=d ₂ ²-(m-x) ²

$x=\frac{{{\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">d</figure-callout>}}_1}^2-{d_2}^2-m^2}{2m},$ ${\mathrm{<figure-callout\; id="1"\; label="y"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">y</figure-callout>}}_1=\sqrt{{{\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">d</figure-callout>}}_1}^2-x^2},$ $y_2=-\sqrt{{{\mathrm{<figure-callout\; id="1"\; label="d"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">d</figure-callout>}}_1}^2-x^2}$

Can judge by the sector directions angle of two base stations, the intersection point that is taken in the overlapping area of coverage is the terminal physical location, such as figure below, calculates two intersection points of A, B, branch is based upon on BTS1-SECTOR1 and the BTS2-SECTOR2, covers angle by the two and can determine that A is intersection point.If but two intersection points are all in covering angle, which can't be distinguished is real intersection point, gets the mid point of intersection point line as positioning result.

S506: when judgement had branch more than three tunnel to report TOA, the signal arrival time difference TDOA by the TOA of branch and reference radio channel TOA difference structure two-way branch used hyperbolic method to analyze the position of terminal.

If only obtain two TDOA, the two pairs of hyperbolas intersect and are bound to obtain significant solution so; But in cdma system, common pilot signal is only measured TOA.TOA by adjacent cell i _iThe TOA of measured value and reference cell ₁Measured value removes to construct TDOA _I1Come the locating terminal physical location.

When existing more than or equal to 3 effective TOA measured values, can construct more than or equal to 3 TDOA.As shown above, at this moment, three branch of a curve intersect may not can in a bit.And also there is certain error in each TDOA value of structure, and it is banded that curve is.So need to adopt optimized mode to solve an equation, this paper adopts weighted least-squares method.Principle TDOA value error is larger, and the weight that then gives this equation is less, thereby so that final MS positioning result is optimum.Make a concrete analysis of as follows:

Step 1 obtains initial solution:

In CDMA, the TOA of known adjacent cell i _iThe TOA of measured value and reference cell ₁Measured value just can obtain

${r_i}^0\&ap;r_i=c\&times;{\mathrm{TOA}}_i,$ Directly obtain initial solution:

At first be defined as follows:

Z _p＝[x,y] ^T

Z _a＝[Z _p ^T,r ₁] ^T＝[x，y，r ₁] ^T

The formula 5 of deriving in the following Chen Shi of the being algorithm rudimentary knowledge)

r _i,1 ²+2r _i,1r ₁＝-2x _i,1x-2y _i,1y+k _i-k ₁ i＝2,3,.M

Add reference cell formula r ₁=r ₁, as follows from writing following formula:

$\left\{\begin{array}{cc}{r}_1-{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_1=0& \\ \frac{1}{2}({r_{i,1}}^2-k_i+k_1)+x_{i,1}x+y_{i,1}y+r_{i,1}r=0& \end{array}\right.,i=\mathrm{2,3},\&CenterDot;\&CenterDot;\&CenterDot;M$

Get final product equationof structure: h-G _aZ _aThe concrete meaning of each symbol is as follows in=0 equation:

$0={\left[\begin{array}{c}0\\ 0\\ 0\\ ...\\ 0\end{array}\right]}_{M\&times;1}$ $h=\frac{1}{2}{\left[\begin{array}{c}{2r}_1\\ {{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_{\mathrm{2,1}}}^2-k_2+k_1\\ {{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_{\mathrm{3,1}}}^2-k_3+{\mathrm{<figure-callout\; id="1"\; label="k"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">k</figure-callout>}}_1\\ ...\\ {r_{M,1}}^2-k_M+k_1\end{array}\right]}_{M\&times;1}$ $G_a=-{\left[\begin{array}{ccc}0& 0& -1\\ x_{\mathrm{2,1}}& y_{\mathrm{2,1}}& r_{\mathrm{2,1}}\\ x_{\mathrm{3,1}}& y_{\mathrm{3,1}}& {\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_{\mathrm{3,1}}\\ ...& ...& ...\\ x_{M,1}& y_{M,1}& r_{M,1}\end{array}\right]}_{M\&times;3}$

Suppose Z _a ⁰Expression Z _aActual value, so

Wherein,

Expression equation h-G _aZ _a=0 residual vector, the purpose of any optimization solution are so that residual vector

Weighted mean square and the minimum of each component.According to classical weighted least-squares method, its weighting matrix ψ is residual vector Covariance matrix.

n _TOA=TOA ₁-TOA ₁ ⁰n _{I, 1}TDOA measure error between expression adjacent cell i and reference cell, n measure error vector, then:

Wherein: $\frac{1}{2}({{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_{\mathrm{2,1}}}^{0^2}-k_2+k_1)+x_{\mathrm{2,1}}{x}^0+y_{\mathrm{2,1}}{\mathrm{<figure-callout\; id="0"\; label="y"\; filenames="DEST\_PATH\_HDA00002732089200031.png,DEST\_PATH\_HDA00002732089200032.png"\; state="\{\{state\}\}">y</figure-callout>}}^0+{{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_{\mathrm{2,1}}}^0{{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_1}^0=0$

Again because:

And r ₂ ⁰=r _2,1 ⁰+ r ₁ ⁰

So:

Then can derive:

Wherein: $B={\left[\begin{array}{ccccc}1& & & & \\ & {r_2}^0& & & \\ & & {r_3}^0& & \\ & & & ...& \\ & & & & {r_M}^0\end{array}\right]}_{M\&times;M}$ $B={\left[\begin{array}{ccccc}1& & & & \\ & {\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_1+r_{\mathrm{2,1}}& & & \\ & & {\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_1+{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_{\mathrm{3,1}}& & \\ & & & ...& \\ & & & & {\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_1+r_{M,1}\end{array}\right]}_{M\&times;M}$

N is TDOA measure error vector

1) expression that can be similar to of measure error vector n becomes the normally distributed random variable for zero-mean;

2) Cn _{I, 1}Much smaller than Namely can ignore

So, have:

ψ is

Covariance matrix,

$Q=\left[\begin{array}{ccccc}{{\mathrm{<figure-callout\; id="1"\; label="w"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">w</figure-callout>}}_1}^2& & & & \\ & {w_2}^2& & & \\ & & {w_3}^2& & \\ & & & \&CenterDot;\&CenterDot;\&CenterDot;& \\ & & & & {w_M}^2\end{array}\right]$

W in the Q matrix ₁ ²It is the measurement variance of reference cell TOA; And w _i ², i 〉=2nd, TDOA _{I, 1}The measurement variance.

According to classical weighted least-squares method, derive slightly herein, can obtain Z _a=[X, y, R ₁] ^TSolution:

Z _a＝argmin{(h-G _aZ _a) ^Tψ ^-1(h-G _aZ _a)}＝(G _a ^Tψ ^-1G _a) ^-1G _a ^Tψ ^-1h

Estimation error:

Step 2 takes full advantage of restriction relation, calculates Z _aVector (x, y, r ₁)

Take full advantage of Z _aVector (x, y, r ₁) restriction relation between three unknown numbers:

$r_i=\sqrt{{(x-x_i)}^2+{(y-y_i)}^2},i=\mathrm{1,2,3}...M$

When i=1, have: r ₁ ²=(x-x ₁) ²+ (y-y ₁) ²

With (x-x ₁) ²With (y-y ₁) ²Consist of new unknown vector Z _a ^'(comprise very reference cell position (x of MS position (x, y) ₁, y ₁))

${Z_a}^{\&prime;}=\left[\begin{array}{c}{(x-x_1)}^2\\ {(y-y_1)}^2\end{array}\right]$ $h^{\&prime;}=\left[\begin{array}{c}{(Z_{a,1}-x_1)}^2\\ {(Z_{a,2}-y_1)}^2\\ {Z_{a,3}}^2\end{array}\right]$ ${G_a}^{\&prime;}=\left[\begin{array}{cc}1& 0\\ 0& 1\\ 1& 1\end{array}\right]$ $0=\left[\begin{array}{c}0\\ 0\\ 0\end{array}\right]$

Then: h '-G _a' Z _a'=0

Vector Z _a' actual value be Z _a' ⁰, construct so residual vector Obtain the covariance matrix of this equation

After, can obtain the weighted least-squares solution, as follows:

Z _a'＝argmin{(h′-G _a'Z _a') ^Tψ′ ^-1(h′-G _a'Z _a')}＝(G _a' ^Tψ′ ^-1G _a') ^-1G _a' ^Tψ′ ^-1h'cov(Z _a')＝(G _a′ ^Tψ' ^-1G _a') ^-1

$B^{\&prime;}=\left[\begin{array}{ccc}{x}^0-x_1& & \\ & y^0-y_1& \\ & & {{\mathrm{<figure-callout\; id="1"\; label="r"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">r</figure-callout>}}_1}^0\end{array}\right]\&ap;\left[\begin{array}{ccc}{Z}_{a,1}-x_1& & \\ & Z_{a,2}-y_1& \\ & & Z_{a,3}\end{array}\right]$

Step 3 obtains the terminal location result

As above obtain vector Z in the step _a' solution after, we just can be according to vector Z _a' definition obtain the solution of MS Location-Unknown number (x, y):

Because: ${Z_a}^{\&prime;}=\left[\begin{array}{c}{(x-x_1)}^2\\ {(y-y_1)}^2\end{array}\right]$ So $Z_p=\left[\begin{array}{c}x\\ y\end{array}\right]=\&PlusMinus;\sqrt{{Z_a}^{\&prime;}}+\left[\begin{array}{c}{x}_1\\ {\mathrm{<figure-callout\; id="1"\; label="y"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">y</figure-callout>}}_1\end{array}\right]$

Following formula may draw a plurality of $Z_p=\left[\begin{array}{c}x\\ y\end{array}\right]$ Solution, finally solve house get can according to whether fall within reference cell be covered as judge.

Step 4 obtains position error and estimates that proof is similar to step 1 herein, and proof slightly herein.

$Z_p=\left[\begin{array}{c}x\\ y\end{array}\right]=\&PlusMinus;\sqrt{{Z_a}^{\&prime;}}+\left[\begin{array}{c}{x}_1\\ {\mathrm{<figure-callout\; id="1"\; label="y"\; filenames="DEST\_PATH\_HDA00002732089200041.png"\; state="\{\{state\}\}">y</figure-callout>}}_1\end{array}\right]$

$\mathrm{cov}\left(Z_p\right)=\frac{1}{4}{B}^{\&prime;\&prime;-1}\mathrm{cov}\left({Z_a}^{\&prime;}\right)B^{\&prime;\&prime;-1},$ Wherein $B^{\&prime;\&prime;}=\left[\begin{array}{cc}(x^0-x_1)& 0\\ 0& (y^0-y_1)\end{array}\right]$

Accordingly, Figure 12 also shows the system configuration schematic diagram of realizing the location in CDMA based on measurement report in the embodiment of the invention, comprises as follows:

The cell number acquisition module, the cell number CELLID in the reporting information that obtains the CDMA product determines cell coverage area;

The time of arrival (toa) acquisition module is used for following the tracks of the time of arrival (toa) TOA that obtains mobile terminal reporting branch at measurement report MR, calculates the distance of portable terminal and base station according to TOA

The positioning analysis module is used for positioning analysis according to the cell coverage area of determining and the distance of portable terminal and base station.

Concrete Figure 13 also shows positioning analysis modular structure schematic diagram, comprises in the positioning analysis module here:

Report the Bifurcation Analysis unit, be used for judging the way situation that reports branch, position analysis according to the way situation that reports branch of judging.

And and/or

The first branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains one the tunnel analyzes when reporting TOA, wakes the first analyzing and positioning processing unit up and processes;

The first analyzing and positioning processing unit is used in conjunction with TOA and sector directions information, and the method for positioning analyzing that use direction of arrival degree AOA combines with TOA analyzes the position of terminal.

And and/or

The second branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains is the two-way analysis when reporting TOA, wakes the second analyzing and positioning processing unit up and processes;

The second analyzing and positioning processing unit is used for using TOA positioning analysis method, analyzes the position of terminal by two circle intersection algorithms.

And and/or

The 3rd branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains is branch more than three tunnel when reporting TOA, wakes the 3rd analyzing and positioning processing unit up and processes;

The 3rd analyzing and positioning processing unit is used for using TOA positioning analysis method, analyzes the position of terminal by two circle intersection algorithms.

The present invention mainly determines the terminal approximate range based on terminal to report reference cell or adjacent cell ID, based on reference pilot and activate the distance that time delay that level branch reports is determined terminal and known location residential quarter.According to the different number of branches of reporting message, select in conjunction with AOA, TOA, the signal parameter localization methods such as TDOA, thus finally determine terminal location.This localization method calculated response speed is fast, can obtain 50 meters to 300 meters positioning accuracy, and this precision is higher than the field intensity positioning accuracy, is higher than based on the characteristic matching positioning accuracy.

Above disclosed is a kind of preferred embodiment of the present invention only, certainly can not limit with this interest field of the present invention, and the equivalent variations of therefore doing according to claim of the present invention still belongs to the scope that the present invention is contained.

Claims (10)

1. a method that realizes the location in CDMA based on measurement report is characterized in that, comprises the steps:

Obtain the cell number CELLID in the reporting information in the CDMA product, determine cell coverage area;

In measurement report MR follows the tracks of, obtain the time of arrival (toa) TOA of mobile terminal reporting branch, calculate the distance of portable terminal and base station according to TOA;

Position analysis according to the cell coverage area of determining and the distance of portable terminal and base station.

2. the method that realizes the location in CDMA based on measurement report as claimed in claim 1 is characterized in that, the cell coverage area that described basis is determined and the distance of portable terminal and base station position in the analysis and comprise:

Judge the way situation that reports branch, position analysis according to the way situation that reports branch of judging.

3. the method that realizes the location in CDMA based on measurement report as claimed in claim 2 is characterized in that, the cell coverage area that described basis is determined and the distance of portable terminal and base station position to analyze and comprise:

When judging when only having one tunnel branch to report TOA, in conjunction with TOA and sector directions information, the method for positioning analyzing that use direction of arrival degree AOA combines with TOA analyzes the position of terminal.

4. the method that realizes the location in CDMA based on measurement report as claimed in claim 2 is characterized in that, the cell coverage area that described basis is determined and the distance of portable terminal and base station position to analyze and comprise:

When judgement has the TOA that two-way branch reports, use TOA positioning analysis method, analyze the position of terminal by two circle intersection algorithms.

5. the method that realizes the location in CDMA based on measurement report as claimed in claim 2 is characterized in that, the cell coverage area that described basis is determined and the distance of portable terminal and base station position to analyze and comprise:

When judgement had branch more than three tunnel to report TOA, the signal arrival time difference TDOA by the TOA of branch and reference radio channel TOA difference structure two-way branch used hyperbolic method to analyze the position of terminal.

6. a system that realizes the location in CDMA based on measurement report is characterized in that, comprising:

The cell number acquisition module, the cell number CELLID in the reporting information that obtains the CDMA product determines cell coverage area;

The time of arrival (toa) acquisition module is used for following the tracks of the time of arrival (toa) TOA that obtains mobile terminal reporting branch at measurement report MR, calculates the distance of portable terminal and base station according to TOA;

The positioning analysis module is used for positioning analysis according to the cell coverage area of determining and the distance of portable terminal and base station.

7. the system that realizes the location in CDMA based on measurement report as claimed in claim 6 is characterized in that, comprises in the described positioning analysis module:

Report the Bifurcation Analysis unit, be used for judging the way situation that reports branch, position analysis according to the way situation that reports branch of judging.

8. the system that realizes the location in CDMA based on measurement report as claimed in claim 7 is characterized in that described locating module comprises:

The first branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains one the tunnel analyzes when reporting TOA, wakes the first analyzing and positioning processing unit up and processes;

The first analyzing and positioning processing unit is used in conjunction with TOA and sector directions information, and the method for positioning analyzing that use direction of arrival degree AOA combines with TOA analyzes the position of terminal.

9. the system that realizes the location in CDMA based on measurement report as claimed in claim 7 is characterized in that described locating module comprises:

The second branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains is the two-way analysis when reporting TOA, wakes the second analyzing and positioning processing unit up and processes;

The second analyzing and positioning processing unit is used for using TOA positioning analysis method, analyzes the position of terminal by two circle intersection algorithms.

10. the system that realizes the location in CDMA based on measurement report as claimed in claim 7 is characterized in that described locating module comprises:

The 3rd branch outcome acquiring unit is used for from reporting the Bifurcation Analysis unit to obtain the result of judgement, if the result who obtains is branch more than three tunnel when reporting TOA, wakes the 3rd analyzing and positioning processing unit up and processes;

The 3rd analyzing and positioning processing unit is used for using TOA positioning analysis method, analyzes the position of terminal by two circle intersection algorithms.

Images