JP2011058928A - System, device, method and program for estimating position - Google Patents

System, device, method and program for estimating position Download PDF

Info

Publication number
JP2011058928A
JP2011058928A JP2009208201A JP2009208201A JP2011058928A JP 2011058928 A JP2011058928 A JP 2011058928A JP 2009208201 A JP2009208201 A JP 2009208201A JP 2009208201 A JP2009208201 A JP 2009208201A JP 2011058928 A JP2011058928 A JP 2011058928A
Authority
JP
Japan
Prior art keywords
distance
mobile station
base stations
error
position estimation
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
JP2009208201A
Other languages
Japanese (ja)
Other versions
JP4957770B2 (en
Inventor
Kurato Maeno
蔵人 前野
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.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry Co Ltd
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 Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP2009208201A priority Critical patent/JP4957770B2/en
Publication of JP2011058928A publication Critical patent/JP2011058928A/en
Application granted granted Critical
Publication of JP4957770B2 publication Critical patent/JP4957770B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To estimate a position of a mobile station with high precision even in an NLOS state or in an environment where a direct wave and a reflected wave are mixedly received. <P>SOLUTION: The position estimating system 10 for estimating a position of a mobile station based on predetermined signals obtained by one-way or two-way communication between a first base station 201 to a fourth base station 204 and the mobile station 100 includes: a distance measuring part 310 for measuring, as a first distance, a distance between each of the plurality of the base stations and the mobile station for each base station based on the obtained predetermined signals; an error estimating part 320 for estimating a bias amount indicating an error amount contained in the first distance between one of the base stations and the mobile station measured by the distance measuring part 310 this time by using a bias amount estimated to be contained in the first distance between the other plurality of the base stations among the plurality of the base stations and the mobile station at the previous time or prior to the previous time; and a position estimating part 340 for estimating the position of the mobile station by using the first distance and the bias amount estimated this time. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、位置推定システム、位置推定装置、位置推定方法及びプログラムに関する。より詳細には、本発明は、移動局と複数の基地局で通信を行い、その通信結果である所定の信号に基づいて移動局の位置を推定する位置推定システム、位置推定装置、位置推定方法及びプログラムに関する。   The present invention relates to a position estimation system, a position estimation apparatus, a position estimation method, and a program. More specifically, the present invention relates to a position estimation system, a position estimation device, and a position estimation method for performing communication between a mobile station and a plurality of base stations and estimating the position of the mobile station based on a predetermined signal as a result of the communication. And the program.

従来から、移動局と複数の基地局で通信を行い、その通信結果である無線信号や受発信時刻などに基づいて移動局の位置を推定する測位システムが提案されている。たとえば、特許文献1〜3には、複数の無線基地局と移動局との間で無線信号を送受信し、移動局の位置を推定するシステムが提案されている。   Conventionally, a positioning system has been proposed in which a mobile station communicates with a plurality of base stations, and the position of the mobile station is estimated on the basis of a radio signal or a transmission / reception time as a communication result. For example, Patent Documents 1 to 3 propose systems for transmitting and receiving radio signals between a plurality of radio base stations and a mobile station and estimating the position of the mobile station.

無線信号としては、電波、音波、光、磁気などを用いたものがある。それらを用いて位置を求める方式のうち主要なものを大別すると、信号の受信強度を用いるRSSI(Received Signal Strength Indicator)方式(特許文献1参照)、信号の伝播時間を用いるTOA(Time Of Arrival)方式(特許文献2参照)、信号の伝播時間差を用いるTDOA(Time Difference Of Arrival)方式、信号の受信方位を用いるAOA(Angle Of Arrival)方式(特許文献3参照)の4方式が挙げられる。   Some wireless signals use radio waves, sound waves, light, magnetism, and the like. The major ones of the methods for obtaining the position using them are roughly classified into the RSSI (Received Signal Strength Indicator) method (refer to Patent Document 1) using the signal reception strength, and the TOA (Time Of Arrival) using the signal propagation time. ) Method (refer to Patent Document 2), TDOA (Time Difference Of Arrival) method using a signal propagation time difference, and AOA (Angle Of Arrival) method (refer to Patent Document 3) using a signal reception direction.

RSSI方式は、特許文献1に示したように、信号の受信強度から距離を推定し、そこから3辺測量で位置を推定する方式である。この方式は、高精度に距離を推定することが難しい一方で、厳密かつ高精度な時刻管理が不要で、様々な機器に造りこみやすいといった特徴がある。受信強度から距離を推定することが難しいため、複数の受信信号の強弱のパターンの学習と識別により位置を推定するものもあるが、これには付近に学習済のエリアがなければ位置を推定できないといった課題がある。   As shown in Patent Document 1, the RSSI method is a method of estimating a distance from a signal receiving intensity and then estimating a position by triangulation. This method is characterized in that it is difficult to estimate the distance with high accuracy, but it does not require strict and high-accuracy time management and can be easily built into various devices. Since it is difficult to estimate the distance from the received intensity, some positions are estimated by learning and discriminating patterns of multiple received signals, but this cannot be estimated unless there is a learned area nearby. There is a problem.

TOA方式は、信号の伝播時間から距離を計測し、そこから3辺測量で位置を推定する方式である。例えば、移動局から基地局に向かって信号を送信した際の信号の伝播時間をtとし(通信方向が逆の場合や、双方向通信のものもある)、信号の伝播速度をcとすると、基地局と移動局間の距離lは、l=c×tで求めることができる。無線信号の一例として電波を媒体としたものの場合、一般に空気中の電波の伝播速度は一定(約30万km/秒)のため、時計の精度が十分に高ければ安定した位置精度を実現しやすい。例えば移動局と基地局の全ての時計が同期している場合、伝播時間tの計測は1回の通信で行なえる。他にも同期なしで距離を求める方法もあるが、その場合には、移動局から発信した電波が基地局に到着し、更に基地局から移動局にて受信した電波を検出する必要があるため処理の負荷が高くなる。   The TOA method is a method of measuring a distance from a signal propagation time and estimating the position by triangulation from there. For example, let t be the signal propagation time when a signal is transmitted from the mobile station to the base station (when the communication direction is reversed or there are two-way communication), and the signal propagation speed is c. The distance l between the base station and the mobile station can be obtained by l = c × t. In the case of a radio signal as an example of a radio signal, the propagation speed of the radio wave in the air is generally constant (approximately 300,000 km / sec), so that stable position accuracy can be easily realized if the clock accuracy is sufficiently high. . For example, when all the clocks of the mobile station and the base station are synchronized, the propagation time t can be measured by one communication. There is another method for obtaining the distance without synchronization, but in that case, the radio wave transmitted from the mobile station arrives at the base station, and it is necessary to detect the radio wave received from the base station at the mobile station. The processing load increases.

例えばTOA方式で用いられる測距の一方式であるTWR(Two Way Ranging)方式では、移動局から基地局に信号を送出し、信号を受信した基地局は直ちに移動局に返信を行ない、移動局は、自局が信号を送出してから基地局からの信号を受信するまでの往復伝播時間から、基地局で折り返しに必要な処理時間を減算して2で割ることにより、片道にかかった伝播時間を求めることができる。いずれにしても、信号の伝播時間を計測するため、高精度な内部時計は必要である。たとえば電波を用いた場合、約1nsの時刻誤差で約30cmの距離誤差が生じる。   For example, in the TWR (Two Way Ranging) method, which is a ranging method used in the TOA method, a signal is transmitted from the mobile station to the base station, and the base station that has received the signal immediately returns a response to the mobile station. Is a one-way propagation by subtracting the processing time required for loopback at the base station and dividing by 2 from the round-trip propagation time from when the station sends the signal to when it receives the signal from the base station. You can ask for time. In any case, a high-accuracy internal clock is necessary to measure the signal propagation time. For example, when radio waves are used, a distance error of about 30 cm occurs with a time error of about 1 ns.

TDOA方式は、TOA方式と同様に信号の伝播時間をベースとしているが、基地局と移動局の距離を計測せずに、移動局からマルチキャスト発信した信号を同時に複数の基地局にて受信し、その受信時間差に基づいて位置を求める。2つの基地局で受信した信号の受信時刻差は、信号の伝播速度を掛け合わせると、移動局から見た距離差となる。2つの基地局から同一の距離差を発生し得る位置は双曲線を描くため、基地局3局以上で同時に信号を受信すると、複数の双曲線が得られ、それらの交点を求めることにより位置が決定する。この方式は、TOA方式と比較し、移動局が信号を受信する必要がない点、1回の測位に必要な信号の送信回数が少ない(通常、移動局からの1回の信号発信で位置が確定する)といった利点がある一方、基地局間で時刻を厳密に同期しなければならず、また基地局に内蔵する時計には必ず誤差があるため、時計の精度にあわせて逐次時刻同期を行なう必要がある。   The TDOA method is based on the propagation time of the signal as in the TOA method, but without the distance between the base station and the mobile station being measured, the signals transmitted from the mobile station by multicast transmission are simultaneously received by a plurality of base stations, The position is obtained based on the reception time difference. The reception time difference between the signals received by the two base stations is a distance difference seen from the mobile station when multiplied by the signal propagation speed. Since a position where the same distance difference can be generated from two base stations draws a hyperbola, when signals are simultaneously received by three or more base stations, a plurality of hyperbolas are obtained, and the positions are determined by obtaining the intersections thereof. . Compared with the TOA method, this method does not require the mobile station to receive a signal, and the number of signal transmissions required for one positioning is small (usually, the position is determined by one signal transmission from the mobile station). The time must be synchronized precisely between base stations, and there is always an error in the clock built in the base station, so time synchronization is performed according to the accuracy of the clock. There is a need.

AOA方式は、複数の基地局が移動局から信号を受信した際の信号の到達角度に基づき、三角測量を行なう方式である。信号の到達角度は、例えば電波の場合は、アレイアンテナ、音波の場合はマイクロフォンアレイなどを用い、MUSIC(MUltiple SIgnal Classification)法などを用いることで推定できる。距離的に非常に近いアンテナ(マイクロフォン)間の受信時刻差を用いるため、非常に高い精度で個々の受信信号をサンプリングする必要があるといった難しさがある。   The AOA method is a method of performing triangulation based on the arrival angle of signals when a plurality of base stations receive signals from mobile stations. The arrival angle of the signal can be estimated by using, for example, an array antenna in the case of radio waves, a microphone array in the case of sound waves, and using a MUSIC (Multiple Signal Classification) method. Since a reception time difference between antennas (microphones) that are very close to each other is used, it is difficult to sample individual received signals with very high accuracy.

このように、各方式には様々な状況により得手、不得手があり、すべてのケースにおいて常に最良となる方式は存在しない。特に、これら無線信号を用いる方式の共通の課題として、反射波による影響がある。   As described above, each method has advantages and disadvantages in various situations, and there is no method that is always the best in all cases. In particular, as a common problem of the systems using these radio signals, there is an influence of reflected waves.

例えば特許文献2では、TOA方式で位置を推定するシステムにおいて、反射波の影響を低減する方式が提案されている。通常、基地局と移動局の間に見通し(LOS: Line of sight)があり、直線で結ばれる経路(図1のD〜D)を通る直接波で通信が行なわれる場合には、直線距離を求めることが可能である。しかしながら、現実には経路途中の構造物Obなどにより反射が発生することがあり、このような反射波を用いて距離を計測すると、直線距離よりも長い距離(図1のD’)が得られてしまう。これは、基地局と移動局を最短で結ぶ経路がお互いを直線で結ぶ直接波の通る経路であり、それ以外の反射波の通る経路はすべて直線よりも長くなるD’>Dの関係が成り立つためである。しかし、この関係を逆に利用することで、なるべく直線の経路に近いものを選ぶことができる。例えば1回の通信の中で直接波と反射波が混在する合成波の場合は、一番初めに到達した成分が直接波による成分であり、この成分により伝播時間を求めることができれば、直線距離を算出することができる。また例えば複数回、信号の伝播時間を計測でき、その中に反射波と直接波による計測値が混在している場合は、最小の計測値を選ぶことで直接波による計測結果を選ぶことができ、それに基づき直線距離が得ることができる。また、この最小値を選ぶという処理は、たとえ直接波による計測結果がなく、すべて反射波によるものであったとしても、直接波と比較して最も差の少ない反射波を選ぶことになり、距離計測の誤差を低減するような直線距離に近い値を選ぶことができる。通信の受信品質をあわせて用いると、多重反射波などの品質の悪い通信を排除でき、より直接波に近いものを選ぶことも可能かもしれない。長期間にわたって継続して計測を行なう場合、移動局の基地局から見た平均移動速度を算出することで、よりもっともらしい位置を算出することも可能である。 For example, Patent Document 2 proposes a method for reducing the influence of reflected waves in a system that estimates the position using the TOA method. Usually, there is a line of sight (LOS) between the base station and the mobile station, and when communication is performed by a direct wave passing through a straight line (D 1 to D 4 in FIG. 1), a straight line is used. It is possible to determine the distance. However, in reality, reflection may occur due to the structure Ob or the like along the route, and when a distance is measured using such a reflected wave, a distance longer than the linear distance (D 1 ′ in FIG. 1) is obtained. It will be. This is a path through which a direct wave connecting the base station and the mobile station in the shortest path connects each other with a straight line, and all other paths through which reflected waves pass are longer than the line D 1 ′> D 1 This is because However, by using this relationship in reverse, it is possible to select a route that is as close to a straight path as possible. For example, in the case of a composite wave in which a direct wave and a reflected wave are mixed in one communication, the component that arrives first is the component due to the direct wave, and if the propagation time can be obtained from this component, the linear distance Can be calculated. In addition, for example, when the signal propagation time can be measured multiple times, and the measurement value by the reflected wave and the direct wave is mixed, the measurement result by the direct wave can be selected by selecting the minimum measurement value. Based on this, a linear distance can be obtained. In addition, the process of selecting this minimum value is to select the reflected wave with the smallest difference compared to the direct wave, even if there is no measurement result of the direct wave and all are based on the reflected wave. A value close to a linear distance that reduces measurement errors can be selected. If the reception quality of communication is used together, poor quality communication such as multiple reflected waves can be eliminated, and it may be possible to select the one closer to the direct wave. When measurement is performed continuously over a long period of time, a more likely position can be calculated by calculating an average moving speed viewed from the base station of the mobile station.

特許文献3では、TDOA方式とAOA方式を組み合わせることで、位置精度を高めようとする技術が開示されている。AOA方式は、角度の分解能が一定の場合、基地局と移動局の距離が離れるに従い、位置の分解能が低下する。一方で基地局と距離が十分に近い場合は、高い位置分解能を有する。かかる特徴により、信号の受信時刻が早いものが、移動局に最も近い基地局であり、位置分解能が高いとし、高い重み付けを行なう。このように時刻により重みを変えて加重平均をとることで、よりもっともらしい位置の推定を行なう。また、特許文献3では、受信時刻を得るため、TDOA方式の測位もあわせて行い、それらを重み付けによって加重平均を取ることにより、もっともらしい位置を計測する方式が提案されている。   Patent Document 3 discloses a technique for increasing the position accuracy by combining the TDOA method and the AOA method. In the AOA scheme, when the angle resolution is constant, the position resolution decreases as the distance between the base station and the mobile station increases. On the other hand, when the distance from the base station is sufficiently close, it has high position resolution. Due to this feature, the base station with the earliest signal reception time is the base station closest to the mobile station, and the position resolution is high, and high weighting is performed. In this way, the most probable position is estimated by changing the weight according to the time and taking the weighted average. Further, Patent Document 3 proposes a method of measuring a plausible position by performing TDOA positioning in order to obtain a reception time and taking a weighted average by weighting them.

特開2009−85780号公報JP 2009-85780 A 特開2001−275148号公報JP 2001-275148 A 特開2007−13500号公報JP 2007-13500 A

しかし、上記の手法は、例えば基地局と移動局の間に見通しがなく(このような状態を以下、NLOS(Non Line of sight)と呼ぶ。)、直接波が届かずに反射波のみで伝播経路が構成される場合や、直接波による測定が十分に行なえない場合に、正しく位置を推定することは非常に難しい。   However, in the above method, for example, there is no line of sight between the base station and the mobile station (such a state is hereinafter referred to as NLOS (Non Line of Light)), and the direct wave does not reach and propagates only by the reflected wave. It is very difficult to estimate the position correctly when a path is constructed or when direct wave measurement cannot be performed sufficiently.

例えば特許文献2の方法では、NLOSの状態にあったとしても、時間的な経過により移動局の移動や環境の変化などで、NLOSの状態が一時的にでも解消し、直接波が到達すると仮定している。例えば屋外で移動局が道路上を継続的に移動しており、ビルなどの建築物の影響で一時的にNLOS状態となっている場合は、時間的経過によってLOS状態(Line of sight:見通しがある状態)とNLOS状態が交互に遷移し、長期的にみれば直接波による距離計測が可能な場合も多くある。このような場合では、文献2でもある程度正しく位置を推定可能である。しかし、たとえばオフィス空間などの屋内環境で人や物体を対象とした位置検知を行う場合では、室内の壁や棚などの構造物によりNLOS状態になると、検知対象が継続的に移動しない場合が多いため、なかなかNLOS状態が解消しないといった課題がある。特許文献2の方法では、このような課題に対処できない。たとえば電波の場合、電波暗室でない屋内など周囲に電波を良く反射する壁面のある環境などでは、直接波と反射波が複雑に混在することになるが、このような場合に電波の強度や雑音などの問題で、直接波を取り逃がす場合がよく起こり得る。このような場合、直接波を受信できる頻度が問題となり、頻度が低い場合には、十分な精度で位置測位が出来ないといった問題も起こる。   For example, in the method of Patent Document 2, even if the NLOS state is present, it is assumed that the NLOS state is temporarily canceled due to the movement of the mobile station or the environment changes with time, and a direct wave arrives. is doing. For example, when the mobile station is continuously moving on the road outdoors and is temporarily in the NLOS state due to the influence of buildings such as buildings, the LOS state (Line of light: outlook There are many cases where distance measurement by direct waves is possible in the long term. In such a case, the position can be estimated correctly to some extent even in Document 2. However, for example, when position detection is performed on a person or object in an indoor environment such as an office space, the detection target often does not move continuously when the NLOS state is caused by a structure such as an indoor wall or shelf. Therefore, there is a problem that the NLOS state is not easily resolved. The method of Patent Document 2 cannot cope with such a problem. For example, in the case of radio waves, direct waves and reflected waves are mixedly mixed in environments where there are wall surfaces that reflect radio waves well, such as indoors that are not anechoic chambers. It is often the case that a direct wave is missed by this problem. In such a case, the frequency with which direct waves can be received becomes a problem. When the frequency is low, there is a problem in that position measurement cannot be performed with sufficient accuracy.

また、特許文献3の方法では、特にAOA方式の分解能が高い基地局近くに移動局がある場合に反射波を拾ってしまうと、移動局が遠方にある場合よりも信号の入射角が大きく変わり、位置精度を大きく低下させる原因となる。   In addition, in the method of Patent Document 3, if a reflected wave is picked up when there is a mobile station near a base station with a high resolution of the AOA method, the incident angle of the signal changes more than when the mobile station is far away. This causes a significant decrease in positional accuracy.

上記問題に鑑み、本発明は、上記のようなNLOSの状態の場合や直接波と反射波が混在して受信するような環境下であっても、高い精度で移動局の位置を推定することを目的とする。   In view of the above problems, the present invention estimates the position of a mobile station with high accuracy even in the above-described NLOS state or in an environment where a direct wave and a reflected wave are mixedly received. With the goal.

上記課題を解決するために、本発明のある観点によれば、位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき該移動局の位置を推定する位置推定システムであって、前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測部と、今回、前記距離計測部により計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する誤差推定部と、前記距離計測部により計測された前記第1の距離、及び前記誤差推定部により今回推定されたバイアス量を用いて前記移動局の位置を推定する位置推定部と、を備える位置推定システムが提供される。   In order to solve the above problems, according to an aspect of the present invention, based on a predetermined signal acquired by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target. A position estimation system for estimating a position of the mobile station, wherein each base station measures a distance between each of the plurality of base stations and the mobile station as a first distance based on the acquired predetermined signal A distance measurement unit, and a bias amount indicating an error amount included in the first distance between one base station and the mobile station among the plurality of base stations measured by the distance measurement unit, Or an error estimating unit that estimates using a bias amount estimated to be included in the first distance between the plurality of base stations and the mobile station among the plurality of base stations before the previous time, and The first measured by the distance measuring unit Distance, and the position estimation system and a position estimating section for estimating the position of the mobile station using the current estimated bias amount by the error estimation unit is provided.

前記距離計測部は、前記所定の信号から各基地局と前記移動局間の信号の伝搬時間を計測し、計測された伝搬時間に所与の伝搬速度を乗算することによって前記第1の距離をそれぞれ算出してもよい。   The distance measuring unit measures a propagation time of a signal between each base station and the mobile station from the predetermined signal, and multiplies the measured propagation time by a given propagation speed to obtain the first distance. Each may be calculated.

前記誤差推定部は、前記第1の距離に基づき各基地局及び前記移動局間の確からしい距離を推定するための複数の状態ベクトルの密度分布を推定し、前記推定された状態ベクトルの密度分布に基づき前記各基地局及び前記移動局間の確からしい距離として第2の距離を算出し、前記位置推定部は、前記第2の距離に基づき前記移動局の位置を推定してもよい。状態ベクトルは、バイアス量や、第2の距離である推定距離を要素に持つ1次元以上のベクトルであってよい。また、複数の電波の伝播路が存在する場合にも精度良くバイアス量と第2の距離を推定するため、状態ベクトルにバイアス量と推定距離の2要素をもった2次元以上の多次元ベクトルとしても良い。   The error estimation unit estimates a density distribution of a plurality of state vectors for estimating a probable distance between each base station and the mobile station based on the first distance, and the density distribution of the estimated state vectors The second distance may be calculated as a probable distance between each of the base stations and the mobile station, and the position estimation unit may estimate the position of the mobile station based on the second distance. The state vector may be a one-dimensional or more vector having elements of the bias amount and the estimated distance that is the second distance. Further, in order to accurately estimate the bias amount and the second distance even when there are a plurality of radio wave propagation paths, the state vector is a two-dimensional or more multidimensional vector having two elements of the bias amount and the estimated distance. Also good.

前記誤差推定部は、前記推定された状態ベクトルの密度分布に基づきバイアス量を求め、求められたバイアス量と前記第1の距離とから前記第2の距離を算出してもよい。   The error estimation unit may obtain a bias amount based on the estimated density distribution of the state vector, and calculate the second distance from the obtained bias amount and the first distance.

前記誤差推定部は、前記状態ベクトルの密度分布の分散値と所与の閾値とを比較することにより、前記状態ベクトルの密度分布が収束していると判断した場合、前記求められたバイアス量と前記第1の距離とから前記第2の距離を算出してもよい。   When the error estimation unit determines that the density distribution of the state vector has converged by comparing the variance value of the density distribution of the state vector with a given threshold value, The second distance may be calculated from the first distance.

前記誤差推定部は、前記状態ベクトルの密度分布の分散値と所与の閾値とを比較することにより、前記状態ベクトルの密度分布が収束していないと判断した場合、前記第1の距離をそのまま前記第2の距離としてもよい。   When the error estimation unit determines that the density distribution of the state vector is not converged by comparing the variance value of the density distribution of the state vector and a given threshold, the first distance is used as it is. The second distance may be used.

前記距離計測部は、所定時間毎に前記所定の信号を取得し、取得した所定の信号に基づき所定時間毎に前記第1の距離を計測し、前記誤差推定部は、前記第1の距離が計測される都度、前記状態ベクトル密度分布を更新し、更新された状態ベクトル密度分布に基づき前記第2の距離を算出してもよい。   The distance measurement unit acquires the predetermined signal every predetermined time, measures the first distance every predetermined time based on the acquired predetermined signal, and the error estimation unit determines that the first distance is Each time measurement is performed, the state vector density distribution may be updated, and the second distance may be calculated based on the updated state vector density distribution.

また、上記課題を解決するために、本発明の別の観点によれば、位置が特定される複数の基地局と一方向又は双方向通信することにより取得した所定の信号に基づき自局の位置を推定する移動局であって、前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと自局との距離を第1の距離として基地局毎に計測する距離計測部と、今回、前記距離計測部により計測された前記複数の基地局のうちの一の基地局と前記自局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記自局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する誤差推定部と、前記距離計測部により計測された前記第1の距離及び前記誤差推定部により今回推定されたバイアス量を用いて前記自局の位置を推定する位置推定部と、を備える移動局が提供される。   In order to solve the above problem, according to another aspect of the present invention, the position of the local station is based on a predetermined signal acquired by one-way or two-way communication with a plurality of base stations whose positions are specified. A distance measuring unit that measures the distance between each of the plurality of base stations and the own station as a first distance based on the acquired predetermined signal for each base station, and this time, A bias amount indicating an error amount included in the first distance between one base station of the plurality of base stations measured by the distance measurement unit and the own station is set to the previous or previous time. An error estimation unit that estimates using a bias amount estimated to be included in the first distance between the plurality of other base stations and the own station, and the distance measurement unit Current estimation by the first distance and the error estimation unit Using said bias amount mobile station and a position estimating section for estimating the position of the own station is provided.

また、上記課題を解決するために、本発明の別の観点によれば、位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき該移動局の位置を推定する位置推定方法であって、前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測ステップと、今回、前記距離計測ステップにて計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する誤差推定ステップと、前記距離計測ステップにて算出された前記第1の距離及び前記誤差推定ステップにて今回推定されたバイアス量を用いて前記移動局の位置を推定する位置推定ステップと、を含む位置推定方法が提供される。   In order to solve the above problem, according to another aspect of the present invention, a predetermined acquired by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target. A position estimation method for estimating a position of the mobile station based on a signal, based on the acquired predetermined signal, for each base station with a distance between each of the plurality of base stations and the mobile station as a first distance And a bias indicating an error amount included in the first distance between one of the plurality of base stations measured in the distance measurement step and the mobile station. An error estimation that estimates the amount using the bias amount estimated to be included in the first distance between the mobile station and a plurality of other base stations of the plurality of base stations before or before the previous time Step and the distance measuring step Position estimating method comprising: a position estimation step of estimating a location of the mobile station using the current estimated bias amount at the calculated first distance and the error estimation step Te is provided.

また、上記課題を解決するために、本発明の別の観点によれば、位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測処理と、今回、前記距離計測処理により計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する位置推定処理と、前記距離計測処理により算出された前記第1の距離及び前記位置推定処理により今回推定されたバイアス量を用いて前記移動局の位置を推定する位置推定処理と、をコンピュータに実行させるためのプログラムが提供される。   In order to solve the above problem, according to another aspect of the present invention, a predetermined acquired by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target. A program for causing a computer to perform a process of estimating the position of the mobile station based on a signal, based on the acquired predetermined signal, a distance between each of the plurality of base stations and the mobile station is a first A distance measurement process for measuring each base station as a distance, and this time, included in the first distance between one of the plurality of base stations and the mobile station measured by the distance measurement process A bias amount indicating an error amount is used as the bias amount estimated to be included in the first distance between the mobile station and a plurality of other base stations among the plurality of base stations before or before the previous time. Position estimation processing A program for causing a computer to execute a position estimation process for estimating the position of the mobile station using the first distance calculated by the distance measurement process and the bias amount currently estimated by the position estimation process. Is provided.

以上説明したように本発明によれば、NLOS状態の場合や直接波と反射波が混在して受信するような環境下であっても、高い精度で移動局の位置を推定することができる。   As described above, according to the present invention, it is possible to estimate the position of the mobile station with high accuracy even in an NLOS state or in an environment where a direct wave and a reflected wave are mixedly received.

本発明の一実施形態に係る位置推定システムの全体構成図である。1 is an overall configuration diagram of a position estimation system according to an embodiment of the present invention. 同実施形態に係る位置推定装置の内部構成図である。It is an internal block diagram of the position estimation apparatus which concerns on the same embodiment. 同実施形態に係る推定処理を示したフローチャートである。It is the flowchart which showed the estimation process which concerns on the same embodiment. 同実施形態に係る位置推定方法を説明するための図である。It is a figure for demonstrating the position estimation method which concerns on the embodiment. 同実施形態に係る位置推定方法を説明するための図である。It is a figure for demonstrating the position estimation method which concerns on the embodiment. 同実施形態に係る位置推定方法を説明するための図である。It is a figure for demonstrating the position estimation method which concerns on the embodiment. 同実施形態に係る状態ベクトルの収束過程を概念的に示した図である。It is the figure which showed notionally the convergence process of the state vector based on the embodiment.

以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。   Exemplary embodiments of the present invention will be described below in detail with reference to the accompanying drawings. In addition, in this specification and drawing, about the component which has the substantially same function structure, duplication description is abbreviate | omitted by attaching | subjecting the same code | symbol.

[位置推定システムの全体構成]
まず、本発明の一実施形態に係る位置推定システムの全体構成について図1を参照しながら説明する。本実施形態に係る位置推定システム10は、移動局100、第1の基地局201、第2の基地局202、第3の基地局203及び第4の基地局204を有している。位置推定システム10は、移動局100の位置を推定する。第1の基地局201〜第4の基地局204は、位置が既知、もしくは位置推定演算時にその位置が取得でき、移動局100との一方向又は双方向通信が可能である。位置推定システム10は、第1の基地局201〜第4の基地局204と移動局100との通信の結果、送信又は受信により取得した所定の信号に基づき移動局100の位置を推定するようになっている。
[Overall configuration of position estimation system]
First, an overall configuration of a position estimation system according to an embodiment of the present invention will be described with reference to FIG. The position estimation system 10 according to the present embodiment includes a mobile station 100, a first base station 201, a second base station 202, a third base station 203, and a fourth base station 204. The position estimation system 10 estimates the position of the mobile station 100. The first base station 201 to the fourth base station 204 have known positions or can acquire the positions at the time of position estimation calculation, and can perform one-way or two-way communication with the mobile station 100. The position estimation system 10 estimates the position of the mobile station 100 based on a predetermined signal acquired by transmission or reception as a result of communication between the first base station 201 to the fourth base station 204 and the mobile station 100. It has become.

第1の基地局201〜第4の基地局204は、例えば屋内の無線LANの基地局や携帯電話の基地局等と同様に設置されている。また、第1の基地局201〜第4の基地局204は、測位信号にデータを付加することで、通信網としての共用も可能である。また、第1の基地局201〜第4の基地局204の位置を特定する手段があれば、基地局を固定する必要はなく、移動局100と同様に可動式の形態としても良い。   The first base station 201 to the fourth base station 204 are installed in the same manner as, for example, indoor wireless LAN base stations and mobile phone base stations. In addition, the first base station 201 to the fourth base station 204 can be shared as a communication network by adding data to the positioning signal. Further, if there is a means for specifying the positions of the first base station 201 to the fourth base station 204, it is not necessary to fix the base station, and the mobile station 100 may be movable.

移動局100としては、PC(Personal Computer)やセンサーネットワークのノード、携帯電話、PDA(Personal Digital Assistant)、アクティブタグなど、様々な形態への適用が可能である。ここでは、移動局を1局、基地局を4局の例をあげて説明するが、移動局は複数でも良く、まだ基地局も2局以上であれば何局であっても良い。通常、3辺測量による位置の推定には3局以上の基地局が必要だが、電波以外に位置に関連する制約条件(たとえば人がいる領域がある程度決まっている場合など)がある場合、基地局が2局であっても位置を推定することができる場合がある。例えば2つの基地局からの距離がそれぞれxメートルとyメートルとなるエリアは、空間内には環状に存在するが、そのうちオフィスフロア(たとえば床面から1m上部にある平面)やオフィス内で人が移動可能な範囲などの制約条件との交点を求めると、1点が得られる場合がある。このような場合は、2局からでも位置検知が可能である。   The mobile station 100 can be applied to various forms such as a PC (Personal Computer), a sensor network node, a mobile phone, a PDA (Personal Digital Assistant), and an active tag. Here, an example in which there is one mobile station and four base stations will be described. However, there may be a plurality of mobile stations, and any number of base stations may be used as long as there are two or more base stations. Normally, three or more base stations are required for position estimation by triangulation. However, if there are restrictions related to position (for example, the area where people are present is determined to some extent) in addition to radio waves, the base station Even if there are two stations, the position may be estimated. For example, an area where the distance from two base stations is x meters and y meters, respectively, exists in a ring shape in the space, but there are people on the office floor (for example, a plane 1 m above the floor) or in the office. If an intersection point with a constraint condition such as a movable range is obtained, one point may be obtained. In such a case, position detection is possible even from two stations.

[位置推定装置の内部構成]
図2に、本実施形態に係る位置推定装置300の内部構成(ブロック)図を示す。位置推定装置300は、基地局に内蔵されていてもよいし、移動局に内蔵されていてもよい。基地局や移動局と別体で存在していてもよい。本実施形態では、位置推定装置300は、基地局に内蔵されている。位置推定装置300は、距離計測部310、誤差推定部320、誤差記憶部330及び位置推定部340を有している。
[Internal configuration of position estimation device]
FIG. 2 shows an internal configuration (block) diagram of the position estimation apparatus 300 according to the present embodiment. The position estimation apparatus 300 may be built in the base station or may be built in the mobile station. It may exist separately from the base station or mobile station. In the present embodiment, the position estimation device 300 is built in the base station. The position estimation apparatus 300 includes a distance measurement unit 310, an error estimation unit 320, an error storage unit 330, and a position estimation unit 340.

距離計測部310は、第1の測距部311、第2の測距部312、第3の測距部313及び第4の測距部314を有している。第1の測距部311〜第4の測距部314は、移動局100と第1の基地局201〜第4の基地局204との通信により取得した所定の信号に基づき、第1の基地局201〜第4の基地局204と移動局100との観測距離を基地局毎に計測する。すなわち、第1の測距部311〜第4の測距部314は、第1の基地局201〜第4の基地局204から、移動局100との通信結果(上記所定の信号に相当)を取得し、移動局100とそれぞれの基地局との観測距離を出力する。たとえば第1の測距部311は、移動局100と第1の基地局201との通信結果から観測距離(図1のD’)を算出する。同様に、第2の測距手段302、第3の測距部313及び第4の測距部314は、移動局100との通信結果から観測距離D、D、Dを算出する。 The distance measuring unit 310 includes a first distance measuring unit 311, a second distance measuring unit 312, a third distance measuring unit 313, and a fourth distance measuring unit 314. The first ranging unit 311 to the fourth ranging unit 314 are configured to transmit the first base station based on a predetermined signal acquired by communication between the mobile station 100 and the first base station 201 to the fourth base station 204. The observation distance between the station 201 to the fourth base station 204 and the mobile station 100 is measured for each base station. That is, the first ranging unit 311 to the fourth ranging unit 314 send communication results (corresponding to the predetermined signal) from the first base station 201 to the fourth base station 204 to the mobile station 100. Obtaining and outputting the observation distance between the mobile station 100 and each base station. For example, the first distance measuring unit 311 calculates the observation distance (D 1 ′ in FIG. 1) from the communication result between the mobile station 100 and the first base station 201. Similarly, the second distance measuring means 302, the third distance measuring unit 313, and the fourth distance measuring unit 314 calculate observation distances D 2 , D 3 , and D 4 from the communication result with the mobile station 100.

なお、各測距部で距離を求めるために用いた通信結果は、必ずしも基地局から得る必要はなく、移動局側で収集したのち、特定の基地局へ複数基地局分をまとめて渡してもよいし、別の通信路で測距部へ渡しても良い。ここの通信路やプロトコルについては、いかような方法でも可能である。また、移動局が複数局ある場合は、それぞれの移動局毎に距離の収集と測位演算を行なうが、誤差が十分に収束していない移動局については、近傍の十分に誤差が収束した移動局の誤差情報を流用するなどの方法をとることも可能である。   Note that the communication results used to determine the distance in each distance measuring unit do not necessarily have to be obtained from the base station, and after collecting on the mobile station side, a plurality of base stations can be collectively delivered to a specific base station. Alternatively, it may be passed to the distance measuring unit via another communication path. Any communication path or protocol can be used here. Also, when there are multiple mobile stations, distance collection and positioning calculation are performed for each mobile station. For mobile stations where the error is not sufficiently converged, a mobile station in which the error has converged sufficiently in the vicinity. It is also possible to take a method such as diverting the error information.

第1の測距部311〜第4の測距部314は、前記所定の信号(通信結果)から第1の基地局201〜第4の基地局204及び移動局100間の信号の伝搬時間を計測し、計測された伝搬時間に所与の伝搬速度を乗算することによって観測距離をそれぞれ算出する。算出された第1の基地局201〜第4の基地局204と移動局100との間の観測距離は、何らの障害物がない場合、観測距離D、D、D、Dとなる。しかしながら、本実施形態では、経路途中の構造物Obなどにより反射が発生し、このような反射波を用いて距離が計測されるため、直線距離よりも長い観測距離D’が算出される。 The first ranging unit 311 to the fourth ranging unit 314 calculate the signal propagation time between the first base station 201 to the fourth base station 204 and the mobile station 100 from the predetermined signal (communication result). Each observation distance is calculated by measuring and multiplying the measured propagation time by a given propagation velocity. The calculated observation distances between the first base station 201 to the fourth base station 204 and the mobile station 100 are the observation distances D 1 , D 2 , D 3 , D 4 when there are no obstacles. Become. However, in the present embodiment, reflection occurs due to the structure Ob or the like along the path, and the distance is measured using such a reflected wave, so that an observation distance D 1 ′ longer than the linear distance is calculated.

なお、第1の測距部311〜第4の測距部314は、所定時間毎に前記所定の信号を取得し、取得した所定の信号に基づき所定時間毎に観測距離を計測することを繰り返す。各基地局及び移動局間の観測距離は、距離計測部310により計測される第1の距離に相当する。   Note that the first distance measuring unit 311 to the fourth distance measuring unit 314 repeatedly acquire the predetermined signal every predetermined time and measure the observation distance every predetermined time based on the acquired predetermined signal. . The observation distance between each base station and mobile station corresponds to the first distance measured by the distance measurement unit 310.

誤差推定部320は、第1の誤差推定部321、第2の誤差推定部322、第3の誤差推定部323及び第4の誤差推定部324を有している。第1の誤差推定部321〜第4の誤差推定部324は、第1の測距部311〜第4の測距部314により計測された観測距離と、過去の履歴などから推定した誤差の分布から観測誤差を予測し、距離を推定する。   The error estimator 320 includes a first error estimator 321, a second error estimator 322, a third error estimator 323, and a fourth error estimator 324. The first error estimation unit 321 to the fourth error estimation unit 324 are distributions of errors estimated from the observation distances measured by the first distance measurement unit 311 to the fourth distance measurement unit 314, past history, and the like. Predict the observation error and estimate the distance.

具体的には、第1の誤差推定部321〜第4の誤差推定部324は、今回、第1の測距部311〜第4の測距部314により計測された第1の基地局201〜第4の基地局204のうちの一の基地局と移動局100との観測距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に第1の基地局201〜第4の基地局204のうちの他の複数の基地局と移動局100との観測距離に含まれると推定されたバイアス量を用いて推定する。   Specifically, the first error estimation unit 321 to the fourth error estimation unit 324 are the first base stations 201 to 201 measured by the first distance measurement unit 311 to the fourth distance measurement unit 314 this time. The bias amount indicating the error amount included in the observation distance between one base station of the fourth base stations 204 and the mobile station 100 is set to the first base station 201 to the fourth base station before or before the previous time. The estimation is performed using the bias amount estimated to be included in the observation distance between the plurality of other base stations 204 and the mobile station 100.

バイアス量bは、図1の反射経路D’に示したように、直線距離Dに対して観測距離に含まれる誤差値(伸び量)である。たとえば第1の基地局201と移動局100との間に電波を通さない障害物Obがあると、電波は、たとえば建物Brの壁を反射して基地局まで到達する。このときの観測距離はD’となり、第1の基地局201及び移動局100間の真の距離よりバイアス量bだけ長くなる。 The bias amount b is an error value (elongation amount) included in the observation distance with respect to the linear distance D 1 as shown in the reflection path D 1 ′ of FIG. For example, if there is an obstacle Ob that does not transmit radio waves between the first base station 201 and the mobile station 100, the radio waves, for example, reflect the wall of the building Br and reach the base station. The observation distance at this time is D 1 ′, which is longer than the true distance between the first base station 201 and the mobile station 100 by the bias amount b.

第1の誤差推定部321〜第4の誤差推定部324は、状態ベクトルの密度分布に基づきバイアス量を求め、求められたバイアス量と観測距離とから各基地局と移動局との誤差量を取り除いた各基地局と移動局との確からしい距離を算出する。より具体的には、第1の誤差推定部321〜第4の誤差推定部324は、観測距離に基づき各基地局及び移動局間の確からしい距離を推定するための複数の状態ベクトルの密度分布を推定し、前記推定された状態ベクトルの密度分布に基づき前記各基地局及び前記移動局間の確からしい距離を算出する。なお、第1の距離から各基地局と移動局との誤差量を取り除いた(各基地局と移動局との確からしい)距離は、第2の距離に相当する。   The first error estimator 321 to the fourth error estimator 324 obtain a bias amount based on the density distribution of the state vector, and calculate the error amount between each base station and the mobile station from the obtained bias amount and the observation distance. A probable distance between each removed base station and the mobile station is calculated. More specifically, the first error estimator 321 to the fourth error estimator 324 are density distributions of a plurality of state vectors for estimating a probable distance between each base station and mobile station based on the observation distance. And a probable distance between each base station and the mobile station is calculated based on the estimated density distribution of the state vector. The distance obtained by removing the error amount between each base station and the mobile station from the first distance (probably between each base station and the mobile station) corresponds to the second distance.

第1の誤差推定部321〜第4の誤差推定部324は、状態ベクトルの密度分布の分散値と所与の閾値とを比較することにより、状態ベクトルの密度分布が収束していると判断した場合、求められたバイアス量と観測距離とから第2の距離を算出する。一方、第1の誤差推定部321〜第4の誤差推定部324は、状態ベクトルの密度分布の分散値と所与の閾値とを比較することにより、状態ベクトルの密度分布が収束していないと判断した場合、観測距離をそのまま前記第2の距離とする。   The first error estimator 321 to the fourth error estimator 324 determine that the density distribution of the state vector has converged by comparing the variance value of the density distribution of the state vector with a given threshold value. In this case, the second distance is calculated from the obtained bias amount and the observation distance. On the other hand, the first error estimator 321 to the fourth error estimator 324 compare the distribution value of the density distribution of the state vector with a given threshold value, so that the density distribution of the state vector has not converged. If it is determined, the observation distance is set as the second distance as it is.

第1の誤差推定部321〜第4の誤差推定部324は、前記第1の距離が計測される都度、上記の方法で状態ベクトル密度分布を更新し、更新された状態ベクトル密度分布に基づき第2の距離を算出する。   Each time the first distance is measured, the first error estimation unit 321 to the fourth error estimation unit 324 update the state vector density distribution by the above method, and the first error estimation unit 321 to the fourth error estimation unit 324 update the state vector density distribution based on the updated state vector density distribution. The distance of 2 is calculated.

誤差記憶部330は、第1の誤差記憶部331、第2の誤差記憶部332、第3の誤差記憶部333及び第4の誤差記憶部334を有している。第1の誤差記憶部331〜第4の誤差記憶部334は、たとえばバイアス量等の過去の誤差の推定値や誤差を推定するための状態ベクトルを記憶する。   The error storage unit 330 includes a first error storage unit 331, a second error storage unit 332, a third error storage unit 333, and a fourth error storage unit 334. The first error storage unit 331 to the fourth error storage unit 334 store estimated values of past errors such as bias amounts and state vectors for estimating errors.

よって、今回、第1の測距部311〜第4の測距部314のいずれかにより計測される一の基地局と移動局との観測距離に含まれるバイアス量は、第1の誤差記憶部331〜第4の誤差記憶部334に記憶された、前回又は前回より前に他の複数の基地局と移動局との観測距離に含まれると推定されたバイアス量を用いて推定されることとなる。   Therefore, this time, the bias amount included in the observation distance between the one base station and the mobile station measured by any one of the first distance measurement unit 311 to the fourth distance measurement unit 314 is the first error storage unit. 331 to the fourth error storage unit 334 stored in the previous time or estimated using the bias amount estimated to be included in the observation distance between the plurality of other base stations and the mobile station before the previous time; Become.

位置推定部340は、第1の誤差推定部321〜第4の誤差推定部324が推定する誤差を考慮して移動局100の位置を計算する。より詳しくは、位置推定部340は、第1の測距部311〜第4の測距部314により算出された各基地局及び移動局間の観測距離、及び第1の誤差推定部321〜第4の誤差推定部324のいずれかにより今回推定されたバイアス量を用いて移動局100の位置を推定する。移動局100の位置を推定する際には、算出された第2の距離を用いてもよい。   The position estimation unit 340 calculates the position of the mobile station 100 in consideration of errors estimated by the first error estimation unit 321 to the fourth error estimation unit 324. More specifically, the position estimation unit 340 includes the observation distances between the base stations and the mobile stations calculated by the first ranging unit 311 to the fourth ranging unit 314, and the first error estimation units 321 to 321. The position of the mobile station 100 is estimated using the bias amount estimated this time by any one of the four error estimation units 324. When estimating the position of the mobile station 100, the calculated second distance may be used.

なお、上記位置推定装置300の各部への指令は、専用の制御デバイスあるいはプログラムを実行する図示しないCPUにより実行される。各プログラムや誤差記憶部330が記憶する状態ベクトルは、図示しないROMや不揮発性メモリに予め記憶されている。CPUが、これらのメモリから各プログラムを読み出し実行することにより、距離計測部310、誤差推定部320、位置推定部340の各機能が実現される。   A command to each part of the position estimation apparatus 300 is executed by a dedicated control device or a CPU (not shown) that executes a program. Each program and the state vector stored in the error storage unit 330 are stored in advance in a ROM or non-volatile memory (not shown). The CPU reads out and executes each program from these memories, thereby realizing the functions of the distance measuring unit 310, the error estimating unit 320, and the position estimating unit 340.

[位置推定装置の動作]
次に、本実施形態に係る位置推定装置300において誤差を推定する動作について、図3の推定処理を示したフローチャートを参照しながら説明する。ここでは、第1の測距部311により測定された観測距離から第1の誤差推定部321により推定された誤差(バイアス量)を用いて第1の基地局201と移動局100との最も確からしい距離を推定し、移動局100の位置を推定する。
[Operation of position estimation device]
Next, an operation of estimating an error in the position estimation apparatus 300 according to the present embodiment will be described with reference to a flowchart illustrating the estimation process in FIG. Here, the most probable relationship between the first base station 201 and the mobile station 100 using the error (bias amount) estimated by the first error estimating unit 321 from the observation distance measured by the first distance measuring unit 311. The estimated distance is estimated, and the position of the mobile station 100 is estimated.

推定処理は、状態ベクトル初期化ステップ(S301)、距離計測ステップ(S302)、誤差推定ステップ(尤度算出ステップS303、バイアス量推定ステップS304)、出力ステップ(S305)の順に実行される。   The estimation process is executed in the order of a state vector initialization step (S301), a distance measurement step (S302), an error estimation step (likelihood calculation step S303, bias amount estimation step S304), and an output step (S305).

(状態ベクトル初期化ステップ)
状態ベクトル初期化ステップS301のでは、状態ベクトルが初期状態として取り得る値の範囲にほぼ一様に分布する擬似乱数値や、平均a,分散σのガウス分布に従う乱数などを与えることができる。例えば状態ベクトルの初期状態の一例が、図7(a)に図示されている。この状態ベクトルは、第1の誤差記憶部331に記憶される。
(State vector initialization step)
In the state vector initialization step S301, it is possible to give a pseudo-random value distributed almost uniformly in a range of values that the state vector can take as an initial state, a random number according to a Gaussian distribution with an average a, and a variance σ 2 . For example, an example of the initial state of the state vector is shown in FIG. This state vector is stored in the first error storage unit 331.

状態ベクトルの密度分布は、例えば図7(a)〜図7(g)に粒子として示された多数の状態ベクトルを用いて、前の誤差状態からの予測と現在の観測距離とから現在の誤差状態を推定するために使われる。第1の基地局201と移動局100との観測距離に含まれると推定される誤差(バイアス量)は、一つの状態ベクトルの密度分布で現される。第2の基地局202と移動局100との観測距離に含まれると推定される誤差も、他の一つの状態ベクトルの密度分布で現される。よって、4局の基地局と移動局との間の4つの観測距離に含まれると推定される誤差は、4つの状態ベクトルによってそれぞれ現される。そして、後述するように、今回、いずれかの基地局と移動局との観測距離に含まれると推定される誤差(バイアス量)は、いずれかの基地局以外の3局の基地局と移動局の観測距離に含まれると推定される誤差を示すための3つの状態ベクトルを用いて推定される。状態ベクトルの具体的な用い方については後述する。   For example, the density distribution of the state vector is obtained by using a large number of state vectors shown as particles in FIGS. 7A to 7G, and using the current error from the prediction from the previous error state and the current observation distance. Used to estimate the state. An error (bias amount) estimated to be included in the observation distance between the first base station 201 and the mobile station 100 is expressed by a density distribution of one state vector. The error estimated to be included in the observation distance between the second base station 202 and the mobile station 100 is also expressed by the density distribution of another one state vector. Therefore, errors estimated to be included in the four observation distances between the four base stations and the mobile station are represented by four state vectors, respectively. As will be described later, the error (bias amount) estimated to be included in the observation distance between any one of the base stations and the mobile station this time is the base station of the three stations other than any of the base stations and the mobile station. Are estimated using three state vectors for indicating an error estimated to be included in the observation distance. Specific usage of the state vector will be described later.

(距離計測ステップ)
距離計測ステップS301のでは、距離計測部310が、第1の基地局201〜第4の基地局204と移動局100との送信又は受信により取得した所定の信号に基づき、第1の基地局201〜第4の基地局204と移動局100との観測距離を基地局毎に計測する。以下に観測距離の計測方法について例を挙げて説明する。
(Distance measurement step)
In the distance measurement step S301, the distance measurement unit 310 uses the first base station 201 based on a predetermined signal acquired by transmission or reception between the first base station 201 to the fourth base station 204 and the mobile station 100. The observation distance between the fourth base station 204 and the mobile station 100 is measured for each base station. Below, an example is given and demonstrated about the measuring method of observation distance.

前述したように、第1の基地局201〜第4の基地局204は、それぞれ移動局100と通信を行なう。第1の測距部311〜第4の測距部314は、第1の基地局201〜第4の基地局204のそれぞれと移動局100との通信結果から、観測距離を算出する。   As described above, the first base station 201 to the fourth base station 204 communicate with the mobile station 100, respectively. The first distance measuring unit 311 to the fourth distance measuring unit 314 calculate the observation distance from the communication result between each of the first base station 201 to the fourth base station 204 and the mobile station 100.

第1の測距部311は、第1の基地局201と移動局100との通信結果から、第1の基地局201と移動局100との間の観測距離を求める。TOA方式の場合は、例えば基地局と移動局の間で信号を送受信し、信号の伝播に要した時間に信号の伝播速度をかけることで、観測距離を求める。   The first ranging unit 311 obtains the observation distance between the first base station 201 and the mobile station 100 from the communication result between the first base station 201 and the mobile station 100. In the case of the TOA system, for example, an observation distance is obtained by transmitting / receiving a signal between a base station and a mobile station and multiplying the time required for signal propagation by the signal propagation speed.

たとえば、第1の基地局201からの信号の発信時間をTS、第1の基地局201から移動局100に信号が到達した後、移動局100から信号を発信し、その信号を第1の基地局201が受信した時間をTR、移動局100で信号を受信してから信号を発信するまでの時間をtとし、信号の伝播速度をcとすると、第1の測距部311では、第1の基地局201と移動局100の観測距離Dを、以下のように求める。 For example, the transmission time of the signal from the first base station 201 is TS 1. After the signal arrives at the mobile station 100 from the first base station 201, the signal is transmitted from the mobile station 100 and the signal is transmitted to the first base station 201. The first distance measurement unit 311 assumes that the time received by the base station 201 is TR 1 , the time from when a signal is received by the mobile station 100 to the time when the signal is transmitted is t m , and the signal propagation speed is c. the observation distance D 1 of the mobile station 100 and the first base station 201, determined as follows.

Figure 2011058928
Figure 2011058928

移動局100と第1の基地局201〜第4の基地局204が時刻同期を行なっている場合、移動局100が信号を発信した時刻をTM、第1の基地局201が信号を受信した時刻をTBとすると、第1の基地局201と移動局100の観測距離Dは、以下のように求めることもできる。 When the mobile station 100 and the first base station 201 to the fourth base station 204 are performing time synchronization, the time when the mobile station 100 transmits a signal is TM 1 , and the first base station 201 receives the signal. Assuming that the time is TB 1 , the observation distance D 1 between the first base station 201 and the mobile station 100 can also be obtained as follows.

Figure 2011058928
Figure 2011058928

(観測誤差と誤差推定)
移動局100と第2の基地局202〜第4の基地局204との観測距離D〜Dについても、第2の測距部312〜第4の測距部314により同様に求めることが出来る。また、RSSI方式の場合は、たとえば基地局と移動局が通信を行なった際の電波強度を得て、電波強度と距離との関係から観測距離を求めることができる。
(Observation error and error estimation)
The observation distances D 2 to D 4 between the mobile station 100 and the second base station 202 to the fourth base station 204 can be similarly obtained by the second distance measuring unit 312 to the fourth distance measuring unit 314. I can do it. In the case of the RSSI system, for example, the observation distance can be obtained from the relationship between the radio field intensity and the distance by obtaining the radio field intensity when the base station and the mobile station communicate.

第1の測距部311〜第4の測距部314で得られる観測距離D〜Dは、観測値としての距離であり、実際の真の距離ではなく、様々な要因に基づく観測誤差e〜eを含む。真の距離をG〜Gとすると、以下の関係が成り立つ。 The observation distances D 1 to D 4 obtained by the first distance measurement unit 311 to the fourth distance measurement unit 314 are distances as observation values, and are not actual true distances but observation errors based on various factors. Includes e 1 to e 4 . When the true distance is G 1 to G 4 , the following relationship is established.

Figure 2011058928
ただし、i=1〜4
Figure 2011058928
However, i = 1 to 4

例えばTOA方式の場合は、信号の伝播時間を測定する際のクロックの分解能、揺らぎやクロック自身の誤差のほか、アンテナ(音波の場合はマイクロフォンなど)と信号を検波し受信時間を付与する部分までの伝送路内の伝播時間、検波処理時間等に誤差要因が存在しうる。TWR方式の場合は、通信の折り返しに必要な処理時間の推定誤差などの要因が発生しうる。また電波や音波等を利用するため、壁の反射などによる反射波の影響や多重反射と直接波の合成、壁などによる直接波の遮断や減衰などが起こりえるほか、電波環境・音波環境としての環境ノイズ重畳による検波ミスや誤検波なども発生する。このうち、位置や距離を測定するといった目的からは、反射波の影響が非常に大きく、反射波を捉えてしまうと、反射経路の推定が難しいため直接波の伝播経路(例えば図1のD)より延びた反射波の伝播経路(例えば図1のD’)を、直接波の経路と誤って観測してしまい、大きな観測誤差が発生することになる。 For example, in the case of the TOA system, in addition to the clock resolution when measuring the signal propagation time, fluctuations and errors of the clock itself, up to the part that detects the antenna (such as a microphone in the case of sound waves) and gives the reception time There may be an error factor in the propagation time, detection processing time, etc. in the transmission path. In the case of the TWR method, factors such as an estimation error in processing time necessary for communication return may occur. In addition, because of the use of radio waves and sound waves, the effects of reflected waves due to wall reflections, multiple reflections and direct wave synthesis, direct wave blockage and attenuation due to walls, etc. can occur, Misdetection and false detection due to environmental noise superimposition also occur. Among these, for the purpose of measuring the position and distance, the influence of the reflected wave is very large, and if the reflected wave is caught, it is difficult to estimate the reflected path, so that the direct wave propagation path (for example, D 1 in FIG. 1). ), The propagation path of the reflected wave extending further (for example, D 1 ′ in FIG. 1) is mistakenly observed as the direct wave path, resulting in a large observation error.

反射波による影響では、NLOS環境では常に反射波を捉えてしまい、非常に大きくて量がある程度固定的なバイアス的な誤差を生じうる。また、アンテナと検波回路までの伝播路の伝播時間は、アンテナの付け替えやケーブル長の影響により、量としては小さいもののバイアス的な誤差を発生しうる。このようなバイアス的な誤差は、観測誤差の中ではある程度固定的な成分となるため、バイアス量bと定義し、これと真の観測雑音e’との合成を観測誤差eと考えることができる。 Due to the influence of the reflected wave, the reflected wave is always captured in the NLOS environment, and a very large bias error with a fixed amount can occur. Further, although the propagation time of the propagation path from the antenna to the detection circuit is small in amount due to the change of the antenna and the influence of the cable length, a bias error may occur. Such bias specific errors, since the somewhat fixed components in the observation error is defined as the amount of bias b i, consider the synthesis of this and the true measurement noise e 'i the observation error e i be able to.

Figure 2011058928
Figure 2011058928

NLOS環境で、非常に大きなバイアス量が生じた場合の補償を考えると、バイアス量と観測誤差の関係はb>>e’となることがあり、このような場合は、観測誤差とバイアス量の関係を以下の式のようにみなすこともできる。 Considering compensation when a very large amount of bias occurs in the NLOS environment, the relationship between the amount of bias and the observation error may be b i >> e ′ i . In such a case, the observation error and the bias The relationship of quantity can also be regarded as the following equation.

Figure 2011058928
Figure 2011058928

またクロック分解能に起因する誤差やハードウェア上の熱雑音等、ある程度ランダム状、白色雑音状の誤差が発生しうるが、これはe’に分類可能である。これらのように、複数の要因による誤差が合成された観測誤差が発生する。 In addition, errors such as errors due to clock resolution and thermal noise on hardware may occur to some extent, such as random noise and white noise, which can be classified as e ′ i . As described above, an observation error is generated by combining errors due to a plurality of factors.

前述したとおり、NLOS環境で、非常に大きなバイアス的な誤差が発生した場合、バイアス量b>>e’となることがあり、このような場合に、バイアス量を推定できれば、距離精度や位置精度を飛躍的に向上させることができる。たとえばIR−UWB(Impulse Radio Ultra wideband)無線を利用したTOA/TDOA型の位置検知デバイスでは、反射波の影響などが少ないLOS環境においては、一般に距離や位置の検知精度が50cm程度といわれている。NLOS環境では、基地局の設置間隔を20m間隔とすると、基地局と移動局間の観測距離(第1の距離)に数m〜十数mレベルのバイアス量が発生する場合があり、このバイアス量の影響で位置の検知精度が劣化する。通常、位置検知を行なう場合は、位置が分かっている複数の基地局を利用し、移動局と複数の基地局との通信結果から位置を推定する。また予め設定された条件に基づき、時間的経過を伴った複数回の通信およびそれに伴う位置の推定を行い、移動局の位置の更新を続けるなどの逐次的な方法で位置の推定を連続する。例えば、1秒間隔で基地局と移動局の通信を連続して行ない、1秒間隔で通信結果が得られる都度、移動局の位置を更新する、あるいは、電波条件が悪く通信が失敗しうる環境では、連続的に通信を試行し、通信が成功する都度、位置を更新する方法などを用いる。または、移動局に振動検知センサ等をつけ、移動時に発生する振動を検知する都度通信を行い、移動局の位置を更新するなどの方法を用いてもよい。
例えば、他の複数の基地局と移動局との距離が正しく推定できていると仮定すれば、複数の基地局との通信結果から得られる観測距離(第1の距離)から、個々の基地局と移動局間の直線距離やNLOS等に起因するバイアス量を推定することができる。一般に、基地局がLOS状態かNLOS状態かを判別することはできないが、1つの基地局の観測距離(第1の距離)に含まれるバイアス量は、他の基地局と移動局との間で推定されるバイアス量に基づく推定距離を用いればおおよそ推定できる。事前知識無く1回の推定処理で正しいバイアス量を推定するのは困難だが、逐次的に複数回、観測距離(第1の距離)が得られる都度、この処理を逐次的に繰り返すことで、個々の基地局毎に推定したバイアス量を、時間の経過に伴いより正しいバイアス量に収束させることが可能であり、その結果として、最終的な距離や位置の精度をLOS環境に近づけることができる。よって、以下の誤差推定ステップではバイアス量を推定する。
As described above, when a very large bias error occurs in the NLOS environment, the bias amount b i >> e ′ i may be obtained. In such a case, if the bias amount can be estimated, the distance accuracy or Position accuracy can be improved dramatically. For example, in a TOA / TDOA type position detection device using IR-UWB (Impulse Radio Ultra wideband) radio, it is generally said that the detection accuracy of distance and position is about 50 cm in an LOS environment where the influence of reflected waves is small. . In an NLOS environment, if the base station is set at an interval of 20 m, a bias amount of several m to several tens of m may occur in the observation distance (first distance) between the base station and the mobile station. The position detection accuracy deteriorates due to the amount. Usually, when position detection is performed, a plurality of base stations whose positions are known are used, and the position is estimated from communication results between the mobile station and the plurality of base stations. Further, based on a preset condition, the position estimation is continued by a sequential method such as performing a plurality of times of communication with the passage of time and estimating the position associated therewith and continuing to update the position of the mobile station. For example, an environment in which communication between the base station and the mobile station is continuously performed at intervals of 1 second, and the position of the mobile station is updated every time a communication result is obtained at intervals of 1 second, or communication may fail due to poor radio wave conditions. Then, a method is used in which communication is continuously attempted and the position is updated each time communication is successful. Alternatively, a method may be used in which a vibration detection sensor or the like is attached to the mobile station, communication is performed each time vibration generated during movement is detected, and the position of the mobile station is updated.
For example, if it is assumed that the distances between the plurality of other base stations and the mobile station can be correctly estimated, each base station can be determined from the observation distance (first distance) obtained from the communication result with the plurality of base stations. The amount of bias due to the linear distance between the mobile station and NLOS can be estimated. In general, it is impossible to determine whether a base station is in the LOS state or the NLOS state, but the bias amount included in the observation distance (first distance) of one base station is between the other base station and the mobile station. The estimated distance can be roughly estimated by using the estimated distance based on the estimated bias amount. Although it is difficult to estimate the correct amount of bias with a single estimation process without prior knowledge, each time an observation distance (first distance) is obtained multiple times in succession, it is possible to The bias amount estimated for each base station can be converged to a more accurate bias amount as time passes, and as a result, the final distance and position accuracy can be made closer to the LOS environment. Therefore, the bias amount is estimated in the following error estimation step.

(誤差推定ステップ)
誤差推定ステップS303では、誤差推定部320が、今回、距離計測ステップにより計測された第1の基地局201〜第4の基地局204のうちの一の基地局と移動局100との観測距離に含まれる誤差量をバイアス量として推定する。例えば第1の基地局201と移動局100との観測距離に含まれるバイアス量は、前回又は前回より前に第2の基地局202〜第4の基地局204(すなわち、第1の基地局201以外の基地局)と移動局100との観測距離に含まれると推定されたバイアス量を用いて推定される。
(Error estimation step)
In the error estimation step S303, the error estimation unit 320 sets the observation distance between one base station of the first base station 201 to the fourth base station 204 and the mobile station 100 that has been measured by the distance measurement step this time. The included error amount is estimated as a bias amount. For example, the amount of bias included in the observation distance between the first base station 201 and the mobile station 100 is the second base station 202 to the fourth base station 204 (that is, the first base station 201 before or last time). It is estimated using the bias amount estimated to be included in the observation distance between the mobile station 100 and the other base station).

第1の誤差推定部321〜第4の誤差推定部324では、第1の測距部311〜第4の測距部314から得られた観測距離D〜Dから推定距離g〜gおよび推定されるバイアス量b〜bを推定する。ここでは、第1の誤差推定部321について記載するが、第2の誤差推定部322〜第4の誤差推定部324についても同様な処理でバイアス量を推定することが可能である。 In the first error estimation unit 321 to the fourth error estimation unit 324, the estimated distances g 1 to g from the observation distances D 1 to D 4 obtained from the first distance measurement unit 311 to the fourth distance measurement unit 314. 4 and the estimated bias amounts b 1 to b 4 are estimated. Although the first error estimator 321 is described here, the second error estimator 322 to the fourth error estimator 324 can also estimate the bias amount by the same process.

第1の誤差推定部321では、N個の状態ベクトルを用いてバイアス量bを推定し補正距離Gを求める。Nは、値が大きくなるほど精度良く観測誤差を推定できる可能性が高まるが、一方で処理量がNに比例して高くなるため、10〜数千程度の範囲の値が良く用いられる。 The first error estimator 321 estimates the bias amount b 1 using N state vectors to obtain the correction distance G 1 . As N increases, the possibility that the observation error can be accurately estimated increases. On the other hand, since the processing amount increases in proportion to N, a value in the range of about 10 to several thousand is often used.

N個の状態ベクトルを、以下の式(6)と定義する。

Figure 2011058928
N state vectors are defined as the following equation (6).
Figure 2011058928

ここでtは観測時刻を示し、t=0は初期状態、t=1は、観測開始後最初に得られた観測値に基づく状態ベクトルである。また1個の状態ベクトルの一例を、たとえば以下のように表現できる。

Figure 2011058928
ただし、i=1,…,N Here, t indicates an observation time, t = 0 is an initial state, and t = 1 is a state vector based on an observation value obtained first after the start of observation. An example of one state vector can be expressed as follows, for example.
Figure 2011058928
However, i = 1, ..., N


Figure 2011058928

Figure 2011058928

前述したように、状態ベクトル初期化ステップ(S301)では、初期状態として、取り得る値の範囲にほぼ一様に分布する擬似乱数値や、平均a,分散σのガウス分布に従う乱数などを状態ベクトルに与えた(図7(a)参照)。この状態ベクトルは、第1の誤差記憶部331に記憶される。 As described above, in the state vector initialization step (S301), as an initial state, a pseudo-random value distributed almost uniformly in a range of possible values, a random number according to a Gaussian distribution with an average a, and a variance σ 2 , etc. The vector was given (see FIG. 7A). This state vector is stored in the first error storage unit 331.

誤差推定ステップでは、第1の誤差記憶部331で保持している時刻tの状態ベクトルの分布から、時刻(t+1)の状態ベクトルを予測する。予測した時刻(t+1)の状態ベクトルを、以下のように定義する(図7(b),図7(e)参照)。   In the error estimation step, the state vector at time (t + 1) is predicted from the state vector distribution at time t held in the first error storage unit 331. The state vector at the predicted time (t + 1) is defined as follows (see FIGS. 7B and 7E).

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

また、Pt,xyz(g,b)は、例えば次式で表すことができる。

Figure 2011058928
Further, P t, xyz (g, b) can be expressed by the following equation, for example.
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

以上のようにして、図3の推定処理中のステップS302〜ステップS305が所定時間経過毎に実行される。すなわち、時刻tで位置推定処理を実行した後、次の時刻t+1で新しい観測距離Dt+1,1が得られると、距離計測ステップS302に戻り、ステップS302〜ステップS305の処理を繰り返す。 As described above, step S302 to step S305 in the estimation process of FIG. 3 are executed every predetermined time. That is, after the position estimation process is executed at time t, when a new observation distance D t + 1,1 is obtained at the next time t + 1, the process returns to the distance measurement step S302, and the processes of steps S302 to S305 are repeated.

(位置推定ステップ)
位置推定部340は、距離計測部310により算出された複数の観測距離、及び誤差推定部320により今回推定されたバイアス量を用いて移動局100の位置を推定する。
(Position estimation step)
The position estimation unit 340 estimates the position of the mobile station 100 using the plurality of observation distances calculated by the distance measurement unit 310 and the bias amount currently estimated by the error estimation unit 320.

Figure 2011058928
Figure 2011058928

Figure 2011058928
Figure 2011058928

そこで、位置推定部340は、たとえば、第1の基地局201〜第4の基地局204のうちの2局の組み合わせ(1,4),(2,4),(3,4),(2,3),(1,3),(1,2)について、2局それぞれからの距離の交点を2か所(1点で交わる場合は1か所、交わらない場合は、基地局間を直線で結ぶ直線上で、それぞれの基地局からの推定距離の比で求めた1点)求め、それらの点の重心を移動局100の位置と推定してもよい。また、位置推定部340は、2局それぞれの交点2か所を通る線分を求め、3局〜4局の基地局との距離で得られる2本の線分の交点を位置とする方法(図6参照)により移動局の位置を推定してもよい。また、位置推定部340は、各基地局からの距離の確率分布である各基地局に対応する状態確率モデルの状態ベクトルに含まれる推定距離gの分布を密度関数に置き換え、基地局4局分合成し、もっとも密度の高い位置を移動局100の推定位置としても良い。   Therefore, the position estimation unit 340, for example, combines (1,4), (2,4), (3,4), (2) of two stations out of the first base station 201 to the fourth base station 204. , 3), (1, 3), (1, 2), two intersections of distances from each of the two stations (one if one point intersects, otherwise a straight line between the base stations) (1 point obtained by the ratio of the estimated distances from the respective base stations) on the straight line connecting the two points, and the center of gravity of these points may be estimated as the position of the mobile station 100. In addition, the position estimation unit 340 obtains a line segment passing through two intersections of two stations, and sets a position of the intersection of two line segments obtained by the distance from the base station of the third station to the fourth station ( The position of the mobile station may be estimated by referring to FIG. In addition, the position estimation unit 340 replaces the distribution of the estimated distance g included in the state vector of the state probability model corresponding to each base station, which is the probability distribution of the distance from each base station, with a density function. The position with the highest density may be used as the estimated position of the mobile station 100.

Figure 2011058928
Figure 2011058928

なお、位置推定方法は、以上に説明した方法以外の既存の方法を用いることができる。例えば、特開2009−47556号公報に記載のように、TOA(到来時間)方式とRSS(受信信号強度)方式とにより位置を推定する方法を利用してもよい。この方法では、推定距離と波発生源から各座標までの実際の距離との差の二乗平均平方根を各座標において算出して、二乗平均平方根が最小となる座標を選択することにより位置を推定する。   As the position estimation method, an existing method other than the method described above can be used. For example, as described in Japanese Patent Application Laid-Open No. 2009-47556, a method of estimating a position by a TOA (arrival time) method and an RSS (received signal strength) method may be used. In this method, the root mean square of the difference between the estimated distance and the actual distance from the wave source to each coordinate is calculated for each coordinate, and the position is estimated by selecting the coordinate with the smallest mean square root. .

また、特表2004−530358号公報に記載のように、TOA(距離)とAOA(角度)を組み合わせて位置を推定するようにしてもよい。また、特開2008−249419号公報に記載のように、あらかじめ測定エリアをメッシュで分割(障害物の周辺は大きく)し、対応座標を割り当てておく方法や、特開2006−3187号公報に記載された位置推定方法を用いることもできる。さらに、特表2009−515201号公報に記載されたTOA方式等を利用することも可能である。   Further, as described in JP-T-2004-530358, the position may be estimated by combining TOA (distance) and AOA (angle). Further, as described in Japanese Patent Application Laid-Open No. 2008-249419, a method in which a measurement area is divided in advance with a mesh (the periphery of an obstacle is large) and corresponding coordinates are assigned, or Japanese Patent Application Laid-Open No. 2006-3187 is described. It is also possible to use the estimated position estimation method. Furthermore, it is also possible to use the TOA method described in JP-T-2009-515201.

通常、1個の移動局に対する位置検知は、複数の基地局との通信の結果に基づき遂行されるため、一部の基地局との通信結果がNLOS状態で劣化し、観測距離が実際の直線距離と異なった場合でも、他の複数の基地局との通信結果を用いると、直線距離と比較したNLOSによる観測距離の伸び量(バイアス量)を含む誤差分布を推定できる場合がほとんどである。また、同一基地局との観測距離においても、長期的な分布をみると、移動速度が高速でないという条件下では直線距離を推定できる可能性は非常に高い。そのため、複数の基地局との通信結果から得られる観測距離(第1の距離)から、個々の基地局と移動局間の直線距離やNLOS等に起因するバイアス量を推定することができる。これにより、前記第1の距離及び今回推定されたバイアス量を用いて高い精度で移動局の位置を推定することができる。   Usually, the position detection for one mobile station is performed based on the results of communication with a plurality of base stations, so the communication results with some base stations deteriorate in the NLOS state, and the observation distance is an actual straight line. Even when the distance is different from the distance, in most cases, an error distribution including the amount of observation distance extension (bias amount) by NLOS compared to the linear distance can be estimated by using the communication results with other base stations. In addition, even in the observation distance to the same base station, when looking at the long-term distribution, it is very likely that the linear distance can be estimated under the condition that the moving speed is not high. Therefore, it is possible to estimate a bias amount due to a linear distance between each base station and a mobile station, NLOS, or the like, from an observation distance (first distance) obtained from communication results with a plurality of base stations. Thereby, the position of the mobile station can be estimated with high accuracy using the first distance and the bias amount estimated this time.

以上に説明したように、本実施形態に係る位置推定システム10によれば、見通しのないNLOS環境や、反射波の発生しやすい環境に対する耐性が得られる。また、NLOS環境かどうかを事前に設定する必要もなく、NLOS環境の場合に、NLOS環境の伝播経路の特性等を事前にモデル化して設定する必要もなく、伝播経路が未知である場合でも測位が可能なほか、動的に環境に追随することができるため、環境変化についても、従来方式よりも優れた精度で位置を検知可能である。   As described above, according to the position estimation system 10 according to the present embodiment, resistance to an NLOS environment where there is no line of sight or an environment where reflected waves are likely to be generated can be obtained. In addition, it is not necessary to set in advance whether or not the NLOS environment, and in the case of the NLOS environment, it is not necessary to model and set the characteristics of the propagation path of the NLOS environment in advance, and positioning is possible even when the propagation path is unknown. In addition to being able to follow the environment dynamically, it is possible to detect the position of the environment change with higher accuracy than the conventional method.

通常、基地局を多数配置すれば、確率的に信頼性の高い情報を収集する頻度も高くなり、位置精度を向上しやすいが、本発明を適用することで、基地局の設置数を増やさずに位置精度を高めることが可能となり、設置コストを削減することも可能である。   Usually, if a large number of base stations are arranged, the frequency of collecting highly reliable information is probabilistic and the position accuracy is easily improved, but by applying the present invention, the number of installed base stations is not increased. In addition, it is possible to improve the position accuracy and to reduce the installation cost.

上述した位置推定システムは、次のように応用することができる。例えば、通常、移動局は複数あることがあるが、このような場合、移動局は移動局間でユニークなIDを発信することで、基地局は移動局毎にデータを収集することができる。また、基地局と移動局の通信上の関係は、逆(移動局が基地局との距離を収集する)でも良い。さらに、基地局の位置は、位置を知る方法があれば固定である必要はない。   The position estimation system described above can be applied as follows. For example, normally, there may be a plurality of mobile stations. In such a case, the base station can collect data for each mobile station by transmitting a unique ID between the mobile stations. The communication relationship between the base station and the mobile station may be reversed (the mobile station collects the distance from the base station). Furthermore, the position of the base station does not need to be fixed if there is a method for knowing the position.

なお、上記実施形態において、各部の動作は互いに関連しており、互いの関連を考慮しながら、一連の動作及び一連の処理として置き換えることができる。これにより、位置推定システム及び位置推定装置の実施形態を、位置推定方法の実施形態とすることができる。また、位置推定システム及び位置推定装置の実施形態を、これらが有する機能をコンピュータに実現させるためのプログラムの実施形態とすることができる。   In the above embodiment, the operations of the respective units are related to each other, and can be replaced as a series of operations and a series of processes in consideration of the relationship between each other. Thereby, embodiment of a position estimation system and a position estimation apparatus can be made into embodiment of a position estimation method. Further, the embodiment of the position estimation system and the position estimation apparatus can be an embodiment of a program for causing a computer to realize the functions of these.

これにより、位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき該移動局の位置を推定する位置推定方法であって、前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測ステップと、今回、前記距離計測ステップにより計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する誤差推定ステップと、前記距離計測ステップにより算出された前記第1の距離、及び前記誤差推定ステップにより今回推定されたバイアス量を用いて前記移動局の位置を推定する位置推定ステップと、を含む位置推定方法を提供することができる。   Thus, a position estimation method for estimating the position of a mobile station based on a predetermined signal acquired by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target. , Based on the acquired predetermined signal, a distance measuring step of measuring each base station as a first distance with a distance between each of the plurality of base stations and the mobile station, and this time, measured by the distance measuring step In addition, a bias amount indicating an error amount included in the first distance between one base station of the plurality of base stations and the mobile station is set to be the other of the plurality of base stations before or before the previous time. An error estimation step using a bias amount estimated to be included in the first distance between the plurality of base stations and the mobile station, the first distance calculated by the distance measurement step, and The error It is possible to provide a position estimating method including the position estimation step of estimating a location of the mobile station using the estimated amount of bias current by a constant step.

また、これにより、位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により得られた所定の信号に基づき該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、前記所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測処理と、今回、前記距離計測処理により計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差値であるバイアス量を、前回又は前回より前に推定された前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する位置推定処理と、前記距離計測処理により算出された前記第1の距離、及び前記位置推定処理により今回又は今回以前に推定された1又は2以上の前記バイアス量を用いて前記移動局の位置を推定する位置推定処理と、をコンピュータに実行させるためのプログラムを提供することができる。   In addition, the computer performs processing for estimating the position of the mobile station based on a predetermined signal obtained by one-way or two-way communication between the plurality of base stations whose positions are specified and the mobile station that is the position estimation target. A distance measurement process for measuring, for each base station, a distance between each of the plurality of base stations and the mobile station as a first distance based on the predetermined signal; A bias amount, which is an error value included in the first distance between one of the plurality of base stations measured by the distance measurement processing and the mobile station, is estimated before or before the previous time. In addition, a position estimation process for estimating using a bias amount estimated to be included in the first distance between the plurality of base stations and the mobile station among the plurality of base stations, and the distance measurement process Before calculated A first distance, and a position estimation process for estimating the position of the mobile station using one or two or more of the bias amounts estimated at or before this time by the position estimation process. A program can be provided.

以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。   The preferred embodiments of the present invention have been described in detail above with reference to the accompanying drawings, but the present invention is not limited to such examples. It is obvious that a person having ordinary knowledge in the technical field to which the present invention pertains can come up with various changes or modifications within the scope of the technical idea described in the claims. Of course, it is understood that these also belong to the technical scope of the present invention.

例えば、上記説明では、図3の推定処理中に移動局の位置を推定する処理は含まれていなかった。しかし、移動局の位置推定処理は、図3の推定処理中の出力ステップS305の後に位置推定ステップを設けることにより、繰り返し実行されるようにしてもよい。   For example, in the above description, the process of estimating the position of the mobile station is not included in the estimation process of FIG. However, the position estimation process of the mobile station may be repeatedly executed by providing a position estimation step after the output step S305 in the estimation process of FIG.

10 位置推定システム
100 移動局
201 第1の基地局
202 第2の基地局
203 第3の基地局
204 第4の基地局
310 距離計測部
311 第1の測距部
312 第2の測距部
313 第3の測距部
314 第4の測距部
320 誤差推定部
321 第1の誤差推定部
322 第2の誤差推定部
323 第3の誤差推定部
324 第4の誤差推定部
330 誤差記憶部
331 第1の誤差記憶部
332 第2の誤差記憶部
333 第3の誤差記憶部
334 第4の誤差記憶部
340 位置推定部
DESCRIPTION OF SYMBOLS 10 Position estimation system 100 Mobile station 201 1st base station 202 2nd base station 203 3rd base station 204 4th base station 310 Distance measurement part 311 1st ranging part 312 2nd ranging part 313 Third distance measuring section 314 Fourth distance measuring section 320 Error estimating section 321 First error estimating section 322 Second error estimating section 323 Third error estimating section 324 Fourth error estimating section 330 Error storage section 331 First error storage unit 332 Second error storage unit 333 Third error storage unit 334 Fourth error storage unit 340 Position estimation unit

Claims (12)

位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき該移動局の位置を推定する位置推定システムであって、
前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測部と、
今回、前記距離計測部により計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する誤差推定部と、
前記距離計測部により計測された前記第1の距離、及び前記誤差推定部により今回推定されたバイアス量を用いて前記移動局の位置を推定する位置推定部と、を備える位置推定システム。
A position estimation system that estimates a position of a mobile station based on a predetermined signal acquired by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target,
Based on the acquired predetermined signal, a distance measuring unit that measures the distance between each of the plurality of base stations and the mobile station as a first distance for each base station;
This time, a bias amount indicating an error amount included in the first distance between one of the plurality of base stations measured by the distance measuring unit and the mobile station is set to be the previous time or before the previous time. An error estimation unit that estimates using a bias amount estimated to be included in the first distance between the plurality of base stations and the mobile station among the plurality of base stations;
A position estimation system comprising: a position estimation unit configured to estimate the position of the mobile station using the first distance measured by the distance measurement unit and the bias amount currently estimated by the error estimation unit.
前記距離計測部は、前記所定の信号から各基地局と前記移動局間の信号の伝搬時間を計測し、計測された伝搬時間に所与の伝搬速度を乗算することによって前記第1の距離をそれぞれ算出する請求項1に記載の位置推定システム。   The distance measuring unit measures a propagation time of a signal between each base station and the mobile station from the predetermined signal, and multiplies the measured propagation time by a given propagation speed to obtain the first distance. The position estimation system according to claim 1, wherein the position estimation system calculates each. 前記誤差推定部は、前記第1の距離に基づき各基地局及び前記移動局間の確からしい距離を推定するための複数の状態ベクトルの密度分布を推定し、前記推定された状態ベクトルの密度分布に基づき前記各基地局及び前記移動局間の確からしい距離として第2の距離を算出し、
前記位置推定部は、前記第2の距離に基づき前記移動局の位置を推定する請求項1又は2に記載の位置推定システム。
The error estimation unit estimates a density distribution of a plurality of state vectors for estimating a probable distance between each base station and the mobile station based on the first distance, and the density distribution of the estimated state vectors A second distance as a probable distance between each base station and the mobile station based on
The position estimation system according to claim 1, wherein the position estimation unit estimates the position of the mobile station based on the second distance.
前記状態ベクトルは、推定バイアス量と推定距離を含む2次元以上の多次元ベクトルである請求項3に記載の位置推定システム。   The position estimation system according to claim 3, wherein the state vector is a two-dimensional or more multidimensional vector including an estimated bias amount and an estimated distance. 前記距離計測部は、時刻の経過と共に逐次的に第1の距離計測を複数回行ない、前記誤差推定部は前記距離計測部によって前記第1の距離計測が行なわれる度にバイアス量を推定し、前記位置推定部は前記誤差推定部によってバイアス量の推定が行なわれる度に前記移動体の位置を推定する請求項3又は4に記載の位置推定システム。   The distance measurement unit sequentially performs a first distance measurement a plurality of times as time passes, and the error estimation unit estimates a bias amount each time the first distance measurement is performed by the distance measurement unit, The position estimation system according to claim 3 or 4, wherein the position estimation unit estimates the position of the moving body every time the bias estimation is performed by the error estimation unit. 前記誤差推定部は、前記推定された状態ベクトルの密度分布に基づきバイアス量を求め、求められたバイアス量と前記第1の距離とから前記第2の距離を算出する請求項3〜5のいずれかに記載の位置推定システム。   The error estimation unit obtains a bias amount based on the estimated density distribution of the state vector, and calculates the second distance from the obtained bias amount and the first distance. Crab position estimation system. 前記誤差推定部は、前記状態ベクトルの密度分布の分散値と所与の閾値とを比較することにより、前記状態ベクトルの密度分布が収束していると判断した場合、前記求められたバイアス量と前記第1の距離とから前記第2の距離を算出する請求項6に記載の位置推定システム。   When the error estimation unit determines that the density distribution of the state vector has converged by comparing the variance value of the density distribution of the state vector with a given threshold value, The position estimation system according to claim 6, wherein the second distance is calculated from the first distance. 前記誤差推定部は、前記状態ベクトルの密度分布の分散値と所与の閾値とを比較することにより、前記状態ベクトルの密度分布が収束していないと判断した場合、前記第1の距離をそのまま前記第2の距離とする請求項6又は7のいずれかに記載の位置推定システム。   When the error estimation unit determines that the density distribution of the state vector is not converged by comparing the variance value of the density distribution of the state vector and a given threshold, the first distance is used as it is. The position estimation system according to claim 6, wherein the second distance is set as the second distance. 前記距離計測部は、所定時間毎に前記所定の信号を取得し、取得した所定の信号に基づき所定時間毎に前記第1の距離を計測し、
前記誤差推定部は、前記第1の距離が計測される都度、前記状態ベクトル密度分布を更新し、更新された状態ベクトル密度分布に基づき前記第2の距離を算出する請求項1〜8のいずれか一項に記載の位置推定システム。
The distance measuring unit acquires the predetermined signal every predetermined time, measures the first distance every predetermined time based on the acquired predetermined signal,
The error estimation unit updates the state vector density distribution each time the first distance is measured, and calculates the second distance based on the updated state vector density distribution. The position estimation system according to claim 1.
位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき該移動局の位置を推定する位置推定装置であって、
前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測部と、
今回、前記距離計測部により計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する誤差推定部と、
前記距離計測部により計測された前記第1の距離、及び前記誤差推定部により今回推定されたバイアス量を用いて前記移動局の位置を推定する位置推定部と、を備える位置推定装置。
A position estimation apparatus that estimates a position of a mobile station based on a predetermined signal acquired by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target,
Based on the acquired predetermined signal, a distance measuring unit that measures the distance between each of the plurality of base stations and the mobile station as a first distance for each base station;
This time, a bias amount indicating an error amount included in the first distance between one of the plurality of base stations measured by the distance measuring unit and the mobile station is set to be the previous time or before the previous time. An error estimation unit that estimates using a bias amount estimated to be included in the first distance between the plurality of base stations and the mobile station among the plurality of base stations;
A position estimation device comprising: a position estimation unit configured to estimate the position of the mobile station using the first distance measured by the distance measurement unit and the bias amount currently estimated by the error estimation unit.
位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき該移動局の位置を推定する位置推定方法であって、
前記取得した所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測ステップと、
今回、前記距離計測ステップにて計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差量を示すバイアス量を、前回又は前回より前に前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する誤差推定ステップと、
前記距離計測ステップにて算出された前記第1の距離、及び前記誤差推定ステップにて今回推定されたバイアス量を用いて前記移動局の位置を推定する位置推定ステップと、を含む位置推定方法。
A position estimation method for estimating a position of a mobile station based on a predetermined signal acquired by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target,
Based on the acquired predetermined signal, a distance measuring step for measuring each of the plurality of base stations and the mobile station as a first distance for each base station;
This time, a bias amount indicating an error amount included in the first distance between one of the plurality of base stations measured in the distance measuring step and the mobile station is set to be the previous time or before the previous time. An error estimation step for estimating using a bias amount estimated to be included in the first distance between the plurality of base stations and the mobile station among the plurality of base stations;
A position estimation method comprising: a position estimation step of estimating the position of the mobile station using the first distance calculated in the distance measurement step and the bias amount currently estimated in the error estimation step.
位置が特定される複数の基地局と位置推定対象である移動局との一方向又は双方向通信により得られた所定の信号に基づき該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、
前記所定の信号に基づき、前記複数の基地局のそれぞれと前記移動局との距離を第1の距離として基地局毎に計測する距離計測処理と、
今回、前記距離計測処理により計測された前記複数の基地局のうちの一の基地局と前記移動局との前記第1の距離に含まれる誤差値であるバイアス量を、前回又は前回より前に推定された前記複数の基地局のうちの他の複数の基地局と前記移動局との前記第1の距離に含まれると推定されたバイアス量を用いて推定する位置推定処理と、
前記距離計測処理により算出された前記第1の距離、及び前記位置推定処理により今回又は今回以前に推定された1又は2以上の前記バイアス量を用いて前記移動局の位置を推定する位置推定処理と、をコンピュータに実行させるためのプログラム。
To cause a computer to perform processing for estimating a position of a mobile station based on a predetermined signal obtained by one-way or two-way communication between a plurality of base stations whose positions are specified and a mobile station that is a position estimation target A program,
Based on the predetermined signal, distance measurement processing that measures the distance between each of the plurality of base stations and the mobile station as a first distance for each base station;
This time, a bias amount that is an error value included in the first distance between one of the plurality of base stations measured by the distance measurement process and the mobile station is set to be the previous time or before the previous time. A position estimation process for estimating using a bias amount estimated to be included in the first distance between the plurality of other base stations and the mobile station among the plurality of estimated base stations;
Position estimation processing for estimating the position of the mobile station using the first distance calculated by the distance measurement processing and one or more bias amounts estimated this time or before this time by the position estimation processing. And a program for causing a computer to execute.
JP2009208201A 2009-09-09 2009-09-09 POSITION ESTIMATION SYSTEM, POSITION ESTIMATION DEVICE, POSITION ESTIMATION METHOD, AND PROGRAM Active JP4957770B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009208201A JP4957770B2 (en) 2009-09-09 2009-09-09 POSITION ESTIMATION SYSTEM, POSITION ESTIMATION DEVICE, POSITION ESTIMATION METHOD, AND PROGRAM

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009208201A JP4957770B2 (en) 2009-09-09 2009-09-09 POSITION ESTIMATION SYSTEM, POSITION ESTIMATION DEVICE, POSITION ESTIMATION METHOD, AND PROGRAM

Publications (2)

Publication Number Publication Date
JP2011058928A true JP2011058928A (en) 2011-03-24
JP4957770B2 JP4957770B2 (en) 2012-06-20

Family

ID=43946738

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009208201A Active JP4957770B2 (en) 2009-09-09 2009-09-09 POSITION ESTIMATION SYSTEM, POSITION ESTIMATION DEVICE, POSITION ESTIMATION METHOD, AND PROGRAM

Country Status (1)

Country Link
JP (1) JP4957770B2 (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011215001A (en) * 2010-03-31 2011-10-27 Oki Electric Industry Co Ltd Position estimation device and method of estimating position
KR101231799B1 (en) 2011-12-02 2013-02-08 충남대학교산학협력단 System for calculating self position information
CN104093207A (en) * 2014-07-28 2014-10-08 湘潭大学 Wireless positioning method based on distance-measuring error model correction
JP2014534427A (en) * 2011-10-17 2014-12-18 コミサリア ア レネルジィ アトミーク エ オ ゼネ ルジイ アルテアナティーフCommissariata L’Energie Atomique Et Aux Energies Alternatives Distance estimation method based on RSS measurements using a receiver with limited sensitivity
US9253750B2 (en) 2013-09-27 2016-02-02 Fujitsu Limited Location model updating apparatus and location estimating method
WO2016034832A1 (en) * 2014-09-04 2016-03-10 Khalifa University Of Science, Technology And Research Methods and devices for bias estimation and correction
JP2017207366A (en) * 2016-05-18 2017-11-24 千代田化工建設株式会社 Monitoring target management system
CN107861099A (en) * 2017-11-09 2018-03-30 特斯联(北京)科技有限公司 It is a kind of to determine method of the mobile device relative to the position of channel barrier equipment
US9945953B2 (en) 2014-09-04 2018-04-17 Khalifa University Of Science, Technology And Research Methods and devices for bias estimation and correction
CN111801588A (en) * 2018-03-26 2020-10-20 阿尔卑斯阿尔派株式会社 Position estimation device, position estimation system, position estimation method, and program
CN113625323A (en) * 2021-06-18 2021-11-09 北京千方科技股份有限公司 Vehicle real-time positioning system, method and medium based on vehicle-road cooperation and vehicle
CN116761253A (en) * 2023-08-17 2023-09-15 湘江实验室 UWB weighted positioning method based on triangular area
JP7384137B2 (en) 2020-09-22 2023-11-21 株式会社Soken positioning device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001298764A (en) * 2000-04-11 2001-10-26 Hitachi Kokusai Electric Inc Position calculation method
JP2007208962A (en) * 2005-12-07 2007-08-16 Ekahau Oy Location determining technique
JP2008128726A (en) * 2006-11-17 2008-06-05 Yokohama National Univ Positioning system, device and method using particle filter

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001298764A (en) * 2000-04-11 2001-10-26 Hitachi Kokusai Electric Inc Position calculation method
JP2007208962A (en) * 2005-12-07 2007-08-16 Ekahau Oy Location determining technique
JP2008128726A (en) * 2006-11-17 2008-06-05 Yokohama National Univ Positioning system, device and method using particle filter

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011215001A (en) * 2010-03-31 2011-10-27 Oki Electric Industry Co Ltd Position estimation device and method of estimating position
JP2014534427A (en) * 2011-10-17 2014-12-18 コミサリア ア レネルジィ アトミーク エ オ ゼネ ルジイ アルテアナティーフCommissariata L’Energie Atomique Et Aux Energies Alternatives Distance estimation method based on RSS measurements using a receiver with limited sensitivity
KR101231799B1 (en) 2011-12-02 2013-02-08 충남대학교산학협력단 System for calculating self position information
US9253750B2 (en) 2013-09-27 2016-02-02 Fujitsu Limited Location model updating apparatus and location estimating method
CN104093207A (en) * 2014-07-28 2014-10-08 湘潭大学 Wireless positioning method based on distance-measuring error model correction
CN104093207B (en) * 2014-07-28 2017-09-15 湘潭大学 A kind of wireless location method based on range error Modifying model
US9945953B2 (en) 2014-09-04 2018-04-17 Khalifa University Of Science, Technology And Research Methods and devices for bias estimation and correction
WO2016034832A1 (en) * 2014-09-04 2016-03-10 Khalifa University Of Science, Technology And Research Methods and devices for bias estimation and correction
JP2017207366A (en) * 2016-05-18 2017-11-24 千代田化工建設株式会社 Monitoring target management system
CN107861099A (en) * 2017-11-09 2018-03-30 特斯联(北京)科技有限公司 It is a kind of to determine method of the mobile device relative to the position of channel barrier equipment
CN111801588A (en) * 2018-03-26 2020-10-20 阿尔卑斯阿尔派株式会社 Position estimation device, position estimation system, position estimation method, and program
EP3779497A4 (en) * 2018-03-26 2021-12-22 Alps Alpine Co., Ltd. Position estimation device, position estimation system, position estimation method, and program
CN111801588B (en) * 2018-03-26 2024-01-23 阿尔卑斯阿尔派株式会社 Position estimation device, position estimation system, position estimation method, and program
JP7384137B2 (en) 2020-09-22 2023-11-21 株式会社Soken positioning device
CN113625323A (en) * 2021-06-18 2021-11-09 北京千方科技股份有限公司 Vehicle real-time positioning system, method and medium based on vehicle-road cooperation and vehicle
CN113625323B (en) * 2021-06-18 2022-12-30 北京千方科技股份有限公司 Vehicle real-time positioning system, method and medium based on vehicle-road cooperation and vehicle
CN116761253A (en) * 2023-08-17 2023-09-15 湘江实验室 UWB weighted positioning method based on triangular area
CN116761253B (en) * 2023-08-17 2023-10-20 湘江实验室 UWB weighted positioning method based on triangular area

Also Published As

Publication number Publication date
JP4957770B2 (en) 2012-06-20

Similar Documents

Publication Publication Date Title
JP4957770B2 (en) POSITION ESTIMATION SYSTEM, POSITION ESTIMATION DEVICE, POSITION ESTIMATION METHOD, AND PROGRAM
JP5278365B2 (en) Position estimation apparatus and position estimation method
JP2011214920A (en) Position estimating device, position estimation method, and program
KR101040705B1 (en) The system and method for wireless positioning
JP4968960B2 (en) Adaptive target positioning method, device, system, and target tracking method using the system by integrating TOA (arrival time) and RSS (received signal strength)
KR100671283B1 (en) System and method for asynchronous wireless positioning by ordered transmission
CN103402258B (en) Wi-Fi (Wireless Fidelity)-based indoor positioning system and method
EP2095151B1 (en) A method for localizing remote devices, using acoustical and electromagnetic waves
JP2007013500A (en) Radio terminal position estimating system, position estimating method for radio terminal position estimating system, and data processor
US10567918B2 (en) Radio-location method for locating a target device contained within a region of space
JP2008527394A (en) System and method for positioning using multipath signals
KR101600190B1 (en) Indoor positioning apparatus considering environmental parameters and method thereof
JP2010525300A (en) A method for locating a transmitter using reverse ray tracking.
JP2012083136A (en) Bias error estimation apparatus, bias error estimation method and position estimation apparatus
KR100882590B1 (en) Device and method for measuring location
KR101901407B1 (en) Apparatus and method for determining location
JP4724545B2 (en) Position detection method and position detection system
US20230266425A1 (en) Estimating a Location of a User Equipment
KR100897413B1 (en) Device and method for perceiving location of target
Alba et al. Uncertainty-boosted radiomap-based indoor positioning with RSSI fingerprints
KR100789420B1 (en) Method for accuracy enhancement for uwb indoor non-los location
KR101174199B1 (en) RFID-Based Localization System and Method in Mobile Environment
CN112352167B (en) Positioning device, positioning method and positioning system
Liu et al. Physical measurements
Johansson et al. Estimation of Orientation in a Dual-Tag Ultra Wideband Indoor Positioning System

Legal Events

Date Code Title Description
A977 Report on retrieval

Effective date: 20110810

Free format text: JAPANESE INTERMEDIATE CODE: A971007

A131 Notification of reasons for refusal

Effective date: 20110816

Free format text: JAPANESE INTERMEDIATE CODE: A131

A521 Written amendment

Effective date: 20111012

Free format text: JAPANESE INTERMEDIATE CODE: A523

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Effective date: 20120221

Free format text: JAPANESE INTERMEDIATE CODE: A01

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Effective date: 20120305

Free format text: JAPANESE INTERMEDIATE CODE: A61

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150330

Year of fee payment: 3

R150 Certificate of patent (=grant) or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150