<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>5G networks three-dimensional indoor positioning based on improved Chan algorithm</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Jiazhi Han</institution>
          ,
          <addr-line>Fuxing Yang, Ke Han(Member, IEEE), Zhongliang Deng(Senior Member, IEEE), Jingwen Jiang, Jiawen Zhang, Lingjie Shi, Xiao Fu</addr-line>
          ,
          <institution>Jiawei Fu Beijing University of Posts and Telecommunications Beijing</institution>
          ,
          <country country="CN">China</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>In the 3D solution of indoor positioning based on 5G networks, the traditional 3D Chan algorithm cannot solve the problem when all base stations have the same height. Therefore, an improved Chan algorithm is proposed in this paper. Through the improvement of the Chan algorithm operation matrix, the 3D TDOA value and the height value of the barometric height measurement are used to solve the more accurate fixed-point 3D coordinates; And using the benchmark base station sorting selection and the second weighted least squares error absolute value symbol of the Chan algorithm to achieve the effect of reducing the positioning error and expanding the positioning area. According to the actual physical environment and base station distribution, using the improved Chan algorithm, the position coordinates of any point in the three-dimensional space within the coverage of the base station signal can be located. Experiments show that, when all base stations are at the same height, using the solution system proposed in this paper and the improved Chan algorithm, the positioning accuracy is significantly improved compared with the solution result of the traditional two-dimensional Chan algorithm. The positioning result with deviation less than one meter accounted for 84.36%.</p>
      </abstract>
      <kwd-group>
        <kwd>Indoor Positioning</kwd>
        <kwd>Chan algorithm</kwd>
        <kwd>base station at the same height</kwd>
        <kwd>Solution system</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>
        The arrival of 5G has brought indoor positioning to a new level, such as direct
positioning using large-scale Multiple-Input Multiple-Output [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]; Multi-path Angle of
Arrival and Time of Arrival indoor 3D positioning based on millimeter wave [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. The use
of 5G related technology for indoor positioning has become the development trend of
the future indoor positioning.
      </p>
      <p>
        In the positioning technology based on 5G signals, the solution methods mainly include
TOA and TDOA. The positioning solution method studied in this paper is the TDOA
algorithm which can eliminate the time synchronization by using the difference [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ].
      </p>
      <p>Among the TDOA solution methods, Chan algorithm is a non-recursive algorithm.
Although Chan algorithm can improve the accuracy, it cannot locate the position of the
mobile terminal through the base station of the same height. The existing improved</p>
      <p>
        Chan algorithm methods include residual weighted Chan algorithm improvement [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ],
Non-Line-Of-Sight (NLOS) condition by adding attenuation factor to smooth [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ].
      </p>
      <p>The traditional 3D Chan algorithm requires that the base stations are not all on the
same plane. If the height is the same, the calculation matrix in the Chan algorithm will
be equal to 0, and the solution cannot be solved. This paper proposes an improved
method for the Chan algorithm to locate the point outside the range of the base station，
and solve the problem when base station at the same height, and proposes a complete
indoor 3D spatial positioning system.
2
2.1</p>
    </sec>
    <sec id="sec-2">
      <title>System</title>
    </sec>
    <sec id="sec-3">
      <title>Model</title>
      <sec id="sec-3-1">
        <title>Lab Environment</title>
        <p>This section mainly studies the indoor positioning solution of any point in an
underground garage under LOS conditions. In Fig. 1, The five-pointed star is the distribution
position of the base station, and each base station is 3 meters high and in the same
height plane.
2.2</p>
        <p>Firstly, the received base station signal is processed, and the reference base station
is used to select the reduction error to calculate TDOA value and RSS value. If the
tobe-positioned point is located on the vertical line (surface), the improved RSS
positioning method is used to estimate the initial point, Finally, combined with the Taylor
algorithm, the base station can solve the problem of the position of any mobile terminal
in the 3D space in the signal coverage under the same height.</p>
      </sec>
      <sec id="sec-3-2">
        <title>2.3 Base Station Selection</title>
        <p>After the second weighted least squares estimation of Chan algorithm, the obtained
TDOA measurement error is a square root. For the judgment of positive and negative
of this error, the first least squares measurement error discrimination Method is
proposed. The specific method will be explained in detail later, and a small amplitude
flipping problem caused by the reference base station being too close to the mobile terminal
will occur. Therefore, according to the chip value size of the signal sent by each base
station, the base station with the second chip value is selected as the reference base
station, which is not the closest to the mobile terminal.</p>
      </sec>
      <sec id="sec-3-3">
        <title>2.4 Chan Algorithm Improvement</title>
        <sec id="sec-3-3-1">
          <title>1)3D TDOA value solving 2D Chan algorithm improvement</title>
          <p>When the base stations are of the same height, the height value measured by the air
pressure altimetry technique is used as the height of the mobile terminal. In the
improvement of the first least squares calculation of the Chan algorithm, the calculation
method of the 2D Chan algorithm is used, and one of the arithmetic formulas and the
two matrices are improved. The distance calculation formula of the base station and the
mobile terminal is improved to Equation 1.</p>
          <p>" = $(" − )) + (" − )) + (" − ))</p>
          <p>H is the height of the mobile terminal measured by the air pressure, and (", ", ")
is the coordinate of the ith base station. Since the TDOA is a 3D value, the height
difference obtained by the air pressure measurement is added when calculating the
distance. The distance from the base station to the mobile terminal of the conventional 2D
Chan algorithm in the Ga and h matrices in the first least squares formula of Equation
2 is replaced by the 3D space distance of Equation 1.</p>
          <p>0 = (0 2 Q0 )4 50 2 4 5ℎ (2)</p>
          <p>Where Q is an identity matrix of n-1 dimensions, n is the number of base stations, n
is greater than or equal to 4.</p>
          <p>The first least squares of the improved Chan algorithm is solved. However, the error
caused by the height difference of the three-dimensional TDOA is not introduced, and
the solution result is more accurate. The result obtained is the formula 3.</p>
          <p>0 = (, , r)2 (3)</p>
          <p>In Equation 3, x and y are the obtained x-axis and y-axis coordinates of the mobile
terminal, respectively, and r is the distance from the point to be measured to the
reference base station.</p>
          <p>The above formula is the improvement of the Chan algorithm for solving the 2D
coordinates of the 3D TDOA value, and does not introduce the error caused by the
height difference. In the second weighted least squares calculation, the original 3D
Chan algorithm solving matrix is changed, and the second 3D weighted least squares
of the Chan algorithm is improved. The second 3D weighted least squares formula of
the improved Chan algorithm uses a 3D least squares formula, as shown in Equation 4:
05 = (05 2 505)4 505 2 54 5ℎ5 (4)</p>
          <p>The Ga1 matrix and the Q1 matrix are shown in Equation 5 and Equation 6: matrix
and the Q1 matrix are shown in Equation 5 and Equation 6:
(1)
5 = 4Bcov(0 )B
cov(0 ) = (0 ?4 50 ?)4 5 (6)</p>
          <p>The third column of the B matrix in Equation 5 are improved, and the height value
obtained by the barometric height measurement is introduced, and the matrix inversion
is prevented. The singular value cannot be solved. The specific improved formula is
shown in Equation 7, Equation 8 and Equation 9.</p>
          <p>(5)
( − 5))
0</p>
        </sec>
        <sec id="sec-3-3-2">
          <title>2）Improved Chan algorithm to expand the locator range</title>
          <p>In the original Chan algorithm, only the point coordinates in the area enclosed by
each base station can be located, so the utility is not high. In response to this problem,
this paper proposes the first least squares measurement error discriminant method.
After the first weighted least squares solution of the Chan algorithm, Δx, Δy, and Δz with
positive and negative values are obtained, and the positive and negative values of Δx,
Δy, and Δz are obtained according to the first time, respectively, and the corresponding
correspondence is obtained for the second time. By this method, it is possible to locate
any mobile terminal location within the coverage of the base station signal.</p>
          <p>When the location of the mobile terminal is on the vertical line (surface) in the base
station, the position of the point cannot be located. The RSS algorithm uses the
attenuation characteristics of the signal strength in a uniform medium to roughly estimate the
position of the mobile terminal.</p>
          <p>O|QR5RT"|I4 SU
d = 10</p>
          <p>In Equation 13, d is the distance from the base station to the mobile terminal, rssi is
the received signal strength, A is the signal strength when the transmitting end is
different from the receiving end by 1 meter, and n is the environmental attenuation factor.</p>
          <p>The purpose of the solution here is to give the estimated initial positioning point in
the correct region to perform the iterative optimization solution of the Taylor algorithm.
First, the strength of each base station positioning signal is sorted, and the weighted
(13)
foot of the two base stations is selected according to the field strength, and the vertical
plane is made.</p>
          <p>) = 10WQRR"X4 QRR"YZ</p>
          <p>5T
5</p>
          <p>Then, the geometric center of gravity is calculated according to the space enclosed
by the weighted vertical plane of each of the three base stations, and the initial
positioning point estimated by the RSSI positioning of the three base stations is calculated.</p>
          <p>a "
[\ = ] " _"` (15)</p>
          <p>"b5 "</p>
          <p>In Equation 15, N is the number of RSSI positioning results, and µ is a weighting
factor positively correlated with the signal strength of the base station group, and x, y,
and z are 3D of each group of positioning estimation results. Finally, the positioning
result of each group of base stations is weighted according to the signal strength of the
base station in the group to obtain an optimal initial positioning point.
(14)
3
3.1</p>
        </sec>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>Simulation And Experiment</title>
      <sec id="sec-4-1">
        <title>Simulation Comparison</title>
        <sec id="sec-4-1-1">
          <title>1)Benchmark selection comparison</title>
          <p>In order to prevent other Chan algorithm improvements from affecting the
comparison results, this paper simulates the base stations with different heights. Fig. 2 is a
simulation result of the mobile terminal near the base station, and the triangle represents
the location of the base station. A total of five base stations are set, the circle is the
actual position of the mobile terminal, x is the position of the improved positioning
point. The lower graph in Fig. 2 is an enlarged view of the above graph in the vicinity
of the mobile terminal. The absolute values of the x and z coordinate differences of the
two positioning results are less than 0.01 m, which can be ignored, and the y-coordinate
values differ greatly. When the base station and the mobile terminal is very close, the
decision method that is not sorted by the base station may cause the final result
deviation to be excessive, and the improved decision can avoid this problem and improve the
positioning accuracy.
2)Improved comparison of Chan algorithm for the same height base station</p>
          <p>When all base stations are of the same height, the improved Chan algorithm uses the
standard TDOA value and the test point height value to solve.</p>
          <p>Fig. 3 Comparison of the improved algorithm of the same base station Chan
algorithm</p>
          <p>The triangle represents the location of the base station. All base stations are 3 meters
high. The circle represents the real position of the point to be located. x is the improved
Chan algorithm positioning result in this paper. The upper picture is a
three-dimensional space picture, and the lower picture is a top view of the above three-dimensional
space picture. Since the traditional Chan algorithm cannot solve the TDOA error caused
by the height difference, the positioning result will generate a large offset, However,
the improved Chan algorithm in this paper will also bring the height into the calculation,
which solves the problem that the singular value of the Chan algorithm matrix cannot
be calculated when all the base stations are the same, and the settlement result is
basically the same as the real position of the to-be-positioned point when using the standard
TDOA solution.</p>
        </sec>
        <sec id="sec-4-1-2">
          <title>3)Expanding the positioning range</title>
          <p>In order to prevent other Chan algorithm improvements from affecting the
comparison results, this paper simulates base stations with different heights. The only
difference is the second weighting of the 3D Chan algorithm. The absolute value of each
coordinate axis error obtained by least squares is handled differently.</p>
          <p>Fig. 4 Improved comparison of mobile terminals outside the location of the Chan
algorithm</p>
          <p>When the discriminant condition is not improved, the Chan algorithm locates a point
outside the area enclosed by the base station, and there is a phenomenon that the error
is positive or negative, which leads to double error and cannot accurately locate the
outer point of the area. However, the improved Chan algorithm has no such defect.
3.2</p>
        </sec>
      </sec>
      <sec id="sec-4-2">
        <title>Experimental Test Comparison</title>
        <p>This section collected 50 fixed-point TDOA data randomly and collected 200,000
data for each fixed point to solve, where the environment is basically unchanged and
the device is stationary. Due to space limitations, this paper randomly selects a set of
50 fixed points in the field to carry out the traditional 2D Chan algorithm solution and
the improved Chan algorithm solution for comparative analysis.</p>
        <p>Fig. 5-6 Traditional 2D Chan algorithm solution result vs Improved Chan
algorithm solution result</p>
        <p>Fig. 5 and Fig. 6 show the results of the traditional 2D Chan algorithm solution and
the improved Chan algorithm solution. Multipath interference processing is performed
for the physical environment, and only the LOS signal is taken to reduce the influence
of multipath interference. The triangle in the figure is the location of the base station,
the blue asterisk is the positioning solution result, and the red five-pointed star is the
real position of the test point. Compared with the results of the traditional 2D Chan
algorithm in Fig. 5, the results of the modified Chan algorithm in Fig. 6 are significantly
improved. Moreover, when the solution result exceeds the range enclosed by the base
station, a phenomenon of inversion with respect to the connection boundary of the base
station occurs. As shown in Fig. 5, the positioning result dot is flipped at y = 5.58m.
Through the improvement of the calculation matrix, the same height base station can
also correctly locate the point coordinates of the 3D space.</p>
        <p>Fig. 7 is a CDF diagram of the error distances of the two algorithms. The blue line
is the improved Chan algorithm and the red line is the traditional 2D Chan algorithm.
According to the test results, more than 84.63% of the test results reached the m-level
positioning, which is a qualitative leap compared with the results of the traditional 2D
Chan algorithm.</p>
      </sec>
    </sec>
    <sec id="sec-5">
      <title>Conclusion</title>
      <p>The improved Chan algorithm solves the problem that the same height base station
cannot use the Chan algorithm to accurately locate the mobile terminal position. In the
second weighted least squares of the Chan algorithm, the discriminant of the absolute
value of the error is added, and the positioning range is raised from the point coordinates
in the range enclosed by the base stations to all points in the coverage of the base station
signal. Through the improvement of the above Chan algorithm, the problem that the
same height base station is effectively solved, the positioning range is expanded, the
positioning error is reduced, and the Chan algorithm is more suitable for engineering
applications. The innovative RSS algorithm is used to compensate for the defect that
the Chan algorithm cannot locate the mobile terminal position on the vertical line,
which makes this set of TDOA solution system more widely applicable and can be
applied to various experimental physics.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <given-names>N.</given-names>
            <surname>Garcia</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H.</given-names>
            <surname>Wymeersch</surname>
          </string-name>
          ,
          <string-name>
            <given-names>E. G.</given-names>
            <surname>Larsson</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. M.</given-names>
            <surname>Haimovich</surname>
          </string-name>
          and
          <string-name>
            <given-names>M.</given-names>
            <surname>Coulon</surname>
          </string-name>
          ,
          <article-title>"Direct Localization for Massive MIMO,"</article-title>
          <source>in IEEE Transactions on Signal Processing</source>
          , vol.
          <volume>65</volume>
          , no.
          <issue>10</issue>
          , pp.
          <fpage>2475</fpage>
          -
          <lpage>2487</lpage>
          ,
          <fpage>15</fpage>
          <lpage>May15</lpage>
          ,
          <year>2017</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <given-names>Y.</given-names>
            <surname>Jia</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H.</given-names>
            <surname>Tian</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Fan</surname>
          </string-name>
          and
          <string-name>
            <given-names>B.</given-names>
            <surname>Liu</surname>
          </string-name>
          ,
          <article-title>"Motion Feature and Millimeter Wave Multi-path AoAToA Based 3D Indoor Positioning,"</article-title>
          <source>2018 IEEE 29th Annual International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC)</source>
          , Bologna,
          <year>2018</year>
          , pp.
          <fpage>1</fpage>
          -
          <lpage>7</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Baoping</surname>
            <given-names>Q</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gang</surname>
            <given-names>S</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Zhanrong</surname>
            <given-names>J</given-names>
          </string-name>
          .
          <article-title>An improved TDOA location algorithm</article-title>
          in LOS environment[C].
          <source>Intelligent Human-Machine Systems and Cybernetics (IHMSC)</source>
          ,
          <year>2012</year>
          4th International Conference on. IEEE,
          <year>2012</year>
          ,
          <fpage>P70</fpage>
          -
          <lpage>73</lpage>
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <given-names>Y.</given-names>
            <surname>Zhong</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Liu</surname>
          </string-name>
          and
          <string-name>
            <given-names>X.</given-names>
            <surname>Luo</surname>
          </string-name>
          ,
          <article-title>"</article-title>
          <source>Indoor UWB Location Based on Residual Weighted Chan Algorithm," 2018 7th International Conference on Digital Home (ICDH)</source>
          , Guilin, China,
          <year>2018</year>
          , pp.
          <fpage>274</fpage>
          -
          <lpage>279</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <given-names>S.</given-names>
            <surname>Bao</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Yan</surname>
          </string-name>
          and
          <string-name>
            <given-names>Y.</given-names>
            <surname>Bao</surname>
          </string-name>
          ,
          <article-title>"Three-dimensional indoor positioning based on wireless mobile communication network,"</article-title>
          <source>2017 First International Conference on Electronics Instrumentation &amp; Information Systems (EIIS)</source>
          ,
          <year>Harbin</year>
          ,
          <year>2017</year>
          , pp.
          <fpage>1</fpage>
          -
          <lpage>4</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>