JP7420289B2 - Method and device for visualizing accuracy of position estimation by wireless communication - Google Patents

Method and device for visualizing accuracy of position estimation by wireless communication Download PDF

Info

Publication number
JP7420289B2
JP7420289B2 JP2022576909A JP2022576909A JP7420289B2 JP 7420289 B2 JP7420289 B2 JP 7420289B2 JP 2022576909 A JP2022576909 A JP 2022576909A JP 2022576909 A JP2022576909 A JP 2022576909A JP 7420289 B2 JP7420289 B2 JP 7420289B2
Authority
JP
Japan
Prior art keywords
accuracy
position estimation
calculating
factors
wireless communication
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.)
Active
Application number
JP2022576909A
Other languages
Japanese (ja)
Other versions
JPWO2022157936A1 (en
Inventor
匡史 岩渕
智明 小川
友規 村上
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
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 Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Publication of JPWO2022157936A1 publication Critical patent/JPWO2022157936A1/ja
Application granted granted Critical
Publication of JP7420289B2 publication Critical patent/JP7420289B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Mobile Radio Communication Systems (AREA)

Description

本開示は、無線通信による位置推定の精度を可視化する方法及び装置に関する。 The present disclosure relates to a method and apparatus for visualizing the accuracy of position estimation by wireless communication.

将来の社会基盤として、サイバーフィジカル融合と呼ばれる概念が議論されている。サイバーフィジカル融合は、IoT(Internet-of-Things)技術やAI(Artificial Intelligent)を活用して、実空間の状態把握と解析とによって新たな高付加価値サービスを提供可能にするものである。 A concept called cyber-physical fusion is being discussed as a future social infrastructure. Cyber-physical fusion makes it possible to provide new high-value-added services by utilizing IoT (Internet-of-Things) technology and AI (Artificial Intelligent) to understand and analyze the state of real space.

実空間の状態把握や解析を行うにあたっては、収集される情報が「どこで」「どのような状況で」得られたものであるかが重要であり、位置情報の重要性が増している。無線通信により位置情報を取得する測位技術は以前より多く検討されている。最近においても、例えば非特許文献1に示されるように、The Third Partnership Project (3GPP) Release
16やIEEE 802.11az等の標準規格で議論され、さらなる高精度化に向けて技術的にも注目され続けている。
When grasping and analyzing the state of real space, it is important to know where and under what circumstances the collected information was obtained, and the importance of location information is increasing. BACKGROUND OF THE INVENTION Positioning technology that acquires location information through wireless communication has been widely studied. Recently, as shown in Non-Patent Document 1, The Third Partnership Project (3GPP) Release
It has been discussed in standards such as IEEE 802.16 and IEEE 802.11az, and continues to receive attention from a technical perspective as it aims to further improve accuracy.

これまでに様々な位置推定方法が提案されている。例えば、受信電力を用いる方法、電波の到来角を用いる方法、電波の到来時間やラウンドトリップ時間を用いる方法、機械学習を用いて位置情報に紐づけられたチャネル情報から推定する方法等がある。ただし、いずれの方法も一長一短であるため、組み合わせて利用されることもある。例えば、一般に、受信電力よりも電波の到来時間のほうが変動は少なく、精度がよいとされる。しかし、電波の到来時間を測定するためには、位置推定用の信号を受信する全ての受信点の受信装置が精度よく同期されている必要がある等、実装上の課題がある。 Various position estimation methods have been proposed so far. For example, there are methods that use received power, methods that use the angle of arrival of radio waves, methods that use arrival time and round trip time of radio waves, and methods that use machine learning to estimate from channel information linked to location information. However, each method has advantages and disadvantages, so they may be used in combination. For example, it is generally said that the time of arrival of radio waves has less variation and is more accurate than the received power. However, in order to measure the arrival time of radio waves, there are implementation issues, such as the need for receiving devices at all receiving points that receive position estimation signals to be synchronized with high accuracy.

また、いずれの方法においても、測位を行う環境によって位置推定の精度は異なる。例えば、見通し外環境では電波は反射を繰り返しながら到来するため、電波の到来時間は長くなる傾向にある。そのため、電波の到来時間やラウンドトリップ時間を用いて距離を算出する位置推定方法の場合、送信点と受信点間の距離は実際の距離よりも長く推定されてしまう。また、見通し外でなくとも、ビル街など周辺に電波を反射する構造物が多い環境においては、直接波よりも反射波の電力が大きくなることにより、推定誤差を生じる場合も多い。位置情報の重要性が増すにつれ、位置推定の精度を把握することも重要な要素となりつつある。 Furthermore, in either method, the accuracy of position estimation varies depending on the environment in which positioning is performed. For example, in a non-line-of-sight environment, radio waves arrive while being repeatedly reflected, so the arrival time of radio waves tends to be longer. Therefore, in the case of a position estimation method that calculates distance using the arrival time or round trip time of radio waves, the distance between the transmission point and the reception point is estimated to be longer than the actual distance. Furthermore, even if it is not out of line of sight, in an environment where there are many structures that reflect radio waves in the surrounding area, such as in a building area, estimation errors often occur because the power of the reflected waves is greater than that of the direct waves. As the importance of location information increases, understanding the accuracy of location estimation is also becoming an important factor.

R. Keating, M. Saily, J. Hulkkonen, and J. Karjalainen, “Overview of Positioning in 5G New Radio”, 2019 16th International Symposium on Wireless Communication Systems, Aug. 2019.R. Keating, M. Saily, J. Hulkkonen, and J. Karjalainen, “Overview of Positioning in 5G New Radio”, 2019 16th International Symposium on Wireless Communication Systems, Aug. 2019.

位置情報を用いたデータ解析や、位置情報を用いるサービスの提供者の立場からすると、位置推定の精度はサービスの信頼性に直結する非常に重要な要素である。そのため、サービス開始にあたっては、サービスエリアにおける位置推定の精度を把握したいニーズがある。 From the standpoint of a provider of data analysis using location information or a service using location information, the accuracy of location estimation is a very important element that is directly linked to the reliability of the service. Therefore, when starting a service, there is a need to understand the accuracy of position estimation in the service area.

また上述の通り、位置推定方法によって位置推定の精度は異なる。加えて、事前にサービスエリアにおける位置推定の精度を把握することができれば、サービス開発においていずれの位置推定方法を採用すべきかの判断に活用することもできるようになる。なお、位置推定精度については、開発された位置推定システムの仕様値として知ることができる。しかし、それらはあくまでも目安であり、実際どれほどの精度が達成できるかは使用してみないとわからない。また、使用環境によって推定精度に差が生じるため、仕様値の情報だけでは不十分である。 Furthermore, as described above, the accuracy of position estimation varies depending on the position estimation method. In addition, if the accuracy of position estimation in a service area can be ascertained in advance, it will be possible to use this information to determine which position estimation method should be adopted during service development. Note that the position estimation accuracy can be known as the specification value of the developed position estimation system. However, these are just guidelines, and you won't know how much accuracy you can actually achieve until you use it. Furthermore, because estimation accuracy varies depending on the usage environment, information on specification values alone is insufficient.

本開示は、上記事情に着目してなされたもので、サービスエリア特有の無線通信の環境条件を考慮したうえで、位置推定の精度を可視化することができる技術を提供することを目的とする。 The present disclosure has been made in view of the above circumstances, and aims to provide a technology that can visualize the accuracy of position estimation while taking into account wireless communication environmental conditions specific to a service area.

本開示は、上記目的を達成するため、無線通信による位置推定の精度を可視化する方法を提供する。本方法の第1のステップは、無線通信の環境条件に基づき位置推定の精度に影響する互いに独立な複数の因子についてそれぞれの数値を算出するステップである。本方法の第2のステップは、複数の因子のそれぞれの基準値を用いて複数の因子のそれぞれの数値を無次元量に変換するステップである。本方法の第3のステップは、複数の因子のそれぞれの無次元量を説明変数とし位置推定の精度を目的変数とする一次多項式を用いて位置推定の精度を計算するステップである。本方法の第4のステップは、位置推定の精度の計算結果を表示するステップである。 In order to achieve the above object, the present disclosure provides a method for visualizing the accuracy of position estimation by wireless communication. The first step of the method is to calculate respective numerical values for a plurality of mutually independent factors that influence the accuracy of position estimation based on the environmental conditions of wireless communication. The second step of the method is to convert the numerical value of each of the plurality of factors into a dimensionless quantity using the reference value of each of the plurality of factors. The third step of this method is to calculate the accuracy of position estimation using a first-order polynomial in which the dimensionless quantities of each of the plurality of factors are used as explanatory variables and the accuracy of position estimation is used as an objective variable. The fourth step of the method is to display the calculation result of the accuracy of the position estimation.

また、本開示は、上記目的を達成するため、無線通信による位置推定の精度を可視化する装置を提供する。本装置は、1つ又は複数のプログラムを記憶した1つ又は複数のメモリと、1つ又は複数のメモリと結合された1つ又は複数のプロセッサとを備える。1つ又は複数のプロセッサは、1つ又は複数のプログラムの実行時、少なくとも第1乃至第4の処理を実行する。第1の処理は、無線通信の環境条件に基づき位置推定の精度に影響する互いに独立な複数の因子についてそれぞれの数値を算出する処理である。第2の処理は、複数の因子のそれぞれの基準値を用いて複数の因子のそれぞれの数値を無次元量に変換する処理である。第3の処理は、複数の因子のそれぞれの無次元量を説明変数とし位置推定の精度を目的変数とする一次多項式を用いて位置推定の精度を計算する処理である。第4の処理は、位置推定の精度の計算結果を表示する処理である。 Furthermore, in order to achieve the above object, the present disclosure provides a device that visualizes the accuracy of position estimation by wireless communication. The apparatus includes one or more memories storing one or more programs and one or more processors coupled to the one or more memories. One or more processors execute at least the first to fourth processes when executing one or more programs. The first process is a process of calculating respective numerical values for a plurality of mutually independent factors that influence the accuracy of position estimation based on the environmental conditions of wireless communication. The second process is a process of converting each numerical value of a plurality of factors into a dimensionless quantity using a reference value of each of the plurality of factors. The third process is a process of calculating the accuracy of position estimation using a first-order polynomial in which the dimensionless quantities of each of a plurality of factors are used as explanatory variables and the accuracy of position estimation is used as an objective variable. The fourth process is a process of displaying the calculation result of the position estimation accuracy.

本開示に係る方法及び装置によれば、位置推定の精度に影響する複数の因子のそれぞれの無次元量を説明変数とし位置推定の精度を目的変数とする一次多項式を用いて位置推定の精度を計算し表示することで、位置推定の精度を可視化することができる。 According to the method and device according to the present disclosure, the accuracy of position estimation is calculated using a first-order polynomial whose explanatory variables are dimensionless quantities of each of a plurality of factors that affect the accuracy of position estimation, and the accuracy of position estimation is the objective variable. By calculating and displaying, the accuracy of position estimation can be visualized.

無線通信による位置推定方法の一例を説明するための図である。FIG. 2 is a diagram for explaining an example of a position estimation method using wireless communication. 本開示の第1実施形態の位置推定精度可視化装置を説明するためのブロック図である。FIG. 1 is a block diagram for explaining a position estimation accuracy visualization device according to a first embodiment of the present disclosure. 入力部に入力される情報の一例を説明するための図である。FIG. 3 is a diagram for explaining an example of information input to an input unit. 伝搬計算部による電波のパスの計算例を示す図である。FIG. 6 is a diagram illustrating an example of calculation of a radio wave path by a propagation calculation unit. 伝搬計算部による各パスの受信電力及び到来時間の計算例を示す図である。FIG. 6 is a diagram illustrating an example of calculation of received power and arrival time of each path by a propagation calculation unit. 精度推定部の機能を説明するためのブロック図である。It is a block diagram for explaining the function of an accuracy estimating part. 精度推定部によるパラメータF5の計算の手順を説明するためのフローチャートである。12 is a flowchart for explaining a procedure for calculating a parameter F5 by an accuracy estimator. パラメータF5の計算で行われる処理を具体的に説明するための図である。FIG. 6 is a diagram for specifically explaining the process performed in calculating the parameter F5. 表示部による表示例を示す図である。It is a figure showing the example of a display by a display part. 本開示の第2実施形態の位置推定精度可視化装置を説明するためのブロック図である。FIG. 2 is a block diagram for explaining a position estimation accuracy visualization device according to a second embodiment of the present disclosure. 本開示の各実施形態の位置推定精度可視化装置のハードウェア構成の例を示すブロック図である。FIG. 2 is a block diagram illustrating an example of a hardware configuration of a position estimation accuracy visualization device according to each embodiment of the present disclosure.

1.無線通信による位置推定の方法の概要
図1は、無線通信による位置推定方法の一例を説明するための図である。無線通信による位置推定の対象は、例えば、スマートフォン等の携帯端末や自動車等の移動体を含む移動局2である。無線通信は移動局2と位置を固定された1以上の基地局4との間で行われる。移動局2の位置の推定には、移動局2との間で同時に通信が可能な複数の基地局4が用いられる。
1. Outline of position estimation method using wireless communication FIG. 1 is a diagram for explaining an example of a position estimation method using wireless communication. The target of position estimation by wireless communication is, for example, a mobile station 2 including a mobile terminal such as a smartphone or a mobile object such as a car. Wireless communication is performed between the mobile station 2 and one or more base stations 4 whose locations are fixed. To estimate the position of the mobile station 2, a plurality of base stations 4 that can simultaneously communicate with the mobile station 2 are used.

位置推定では、移動局2は位置推定用の信号を送信する送信点となる。以下、移動局2を送信点2と表記する。位置推定用の信号としては、例えば、SRS(Sounding Reference Signal)が用いられる。移動局2は位置推定用の信号を含む電波を全方位に放射する
。移動局2の周囲に存在する複数の基地局4は、位置推定用の信号を含む電波を受信する受信点となる。以下、基地局4を受信点4と表記する。なお、送信点2と受信点4とはMIMOシステムを構成する。ただし、送信点2と受信点4のそれぞれにおけるアンテナの数及び種類に限定は無い。
In position estimation, the mobile station 2 becomes a transmission point that transmits a signal for position estimation. Hereinafter, the mobile station 2 will be referred to as a transmission point 2. As a signal for position estimation, for example, an SRS (Sounding Reference Signal) is used. The mobile station 2 emits radio waves containing signals for position estimation in all directions. A plurality of base stations 4 existing around the mobile station 2 serve as reception points that receive radio waves including signals for position estimation. Hereinafter, the base station 4 will be referred to as a receiving point 4. Note that the transmission point 2 and the reception point 4 constitute a MIMO system. However, there is no limit to the number and type of antennas at each of the transmitting point 2 and the receiving point 4.

移動局2の位置推定は様々な方法で行うことができる。一つの例として、送信点から受信点までの電波の到来時間を用いる方法が挙げられる。この方法では、送信点2から同時刻に送信された電波を受信した複数(例えば3つ)の受信点4の間で電波の到来時間が比較される。具体的には、各受信点4から電波の到来時間に関する情報が位置推定システム6に送られる。位置推定システム6は、例えば、通信ネットワーク上に配置されたサーバである。位置推定システム6は、各受信点4における電波の到来時間と地図情報とに基づいて地図上における送信点2の現在位置を推定する。 Estimating the position of the mobile station 2 can be performed in various ways. One example is a method that uses the arrival time of radio waves from a transmission point to a reception point. In this method, the arrival times of radio waves are compared between a plurality of (for example, three) reception points 4 that received radio waves transmitted from the transmission point 2 at the same time. Specifically, information regarding the arrival time of radio waves is sent from each reception point 4 to the position estimation system 6. The position estimation system 6 is, for example, a server located on a communication network. The position estimation system 6 estimates the current position of the transmission point 2 on the map based on the arrival time of radio waves at each reception point 4 and map information.

ここに例示された位置推定方法を含めて、これまでに提案されているどの位置推定方法にも推定誤差は含まれる。例えば、ここで例示された位置推定方法は、到来時間を用いる位置推定方法であるが、受信電力を用いる位置推定方法や機械学習による位置推定方法などでも推定誤差は含まれる。本開示の各実施形態の位置推定の精度を可視化する方法及び装置は、位置推定の精度を可視化することによって、評価対象の位置推定方法のアプリケーションへの適用の可否等を容易に判断できるようにしたものである。評価対象となる位置推定方法には限定は無い。 All position estimation methods that have been proposed so far, including the position estimation method exemplified here, include estimation errors. For example, although the position estimation method exemplified here is a position estimation method using time of arrival, estimation errors are also included in position estimation methods using received power, position estimation methods based on machine learning, and the like. The method and device for visualizing the accuracy of position estimation according to each embodiment of the present disclosure makes it possible to easily determine whether or not a position estimation method to be evaluated is applicable to an application by visualizing the accuracy of position estimation. This is what I did. There are no limitations on the position estimation method to be evaluated.

2.第1実施形態の位置推定精度可視化装置
図2は、第1実施形態の位置推定精度可視化装置10を説明するためのブロック図である。位置推定精度可視化装置10は、位置推定の精度を数値化して表示することによって可視化する装置であって、この装置10によって位置推定の精度を可視化する方法が実行される。位置推定精度可視化装置10は、その機能別に、入力部11、伝搬計算部12、精度推定部13、及び表示部14から構成される。以下、位置推定精度可視化装置10を構成する各部の機能について図を参照しながら説明する。
2. Position estimation accuracy visualization device according to the first embodiment FIG. 2 is a block diagram for explaining the position estimation accuracy visualization device 10 according to the first embodiment. The position estimation accuracy visualization device 10 is a device that visualizes the accuracy of position estimation by converting it into numerical values and displaying it, and this device 10 executes a method of visualizing the accuracy of position estimation. The position estimation accuracy visualization device 10 includes an input section 11, a propagation calculation section 12, an accuracy estimation section 13, and a display section 14 according to their functions. Hereinafter, the functions of each part constituting the position estimation accuracy visualization device 10 will be explained with reference to the drawings.

入力部11は、無線通信の環境条件に関する情報の入力を受け付ける。図3は、入力部11に入力される情報の一例を説明するための図である。無線通信の環境条件には、位置推定精度の評価を行うエリアの情報が含まれる。エリアの情報には、図3に示すように、地図、構造物、およびそれらのレイアウトが含まれる。また、無線通信の環境条件には、図3に示すように、送信点2及び受信点4のレイアウトが含まれる。また、位置推定に用いる無線通信システムの無線パラメータ、位置推定方法、及び関連するパラメータも無線通信の環境条件に含まれる。 The input unit 11 receives input of information regarding environmental conditions of wireless communication. FIG. 3 is a diagram for explaining an example of information input to the input unit 11. The environmental conditions for wireless communication include information on the area in which position estimation accuracy is evaluated. As shown in FIG. 3, the area information includes maps, structures, and their layouts. Furthermore, the environmental conditions for wireless communication include the layout of the transmitting point 2 and the receiving point 4, as shown in FIG. Furthermore, the wireless parameters of the wireless communication system used for position estimation, the position estimation method, and related parameters are also included in the wireless communication environmental conditions.

伝搬計算部12は、入力部11に入力された情報を用いて送信点から受信点への電波の伝搬路の計算を行う。伝搬路の計算方法として、例えばレイトレース法を用いることができる。レイトレース法を用いることで、地図及び構造物に対応した伝搬路を計算することができる。以下の説明では、伝搬路の計算方法にレイトレース法が用いられるものとする。ただし、レイトレース法はあくまでも伝搬路の計算方法の一例であり、他の方法により伝搬路を計算してもよい。 The propagation calculation unit 12 uses the information input to the input unit 11 to calculate the propagation path of radio waves from the transmission point to the reception point. For example, a ray tracing method can be used as a method for calculating the propagation path. By using the ray tracing method, propagation paths corresponding to maps and structures can be calculated. In the following description, it is assumed that the ray tracing method is used as the propagation path calculation method. However, the ray tracing method is just one example of a method for calculating the propagation path, and the propagation path may be calculated using other methods.

図4は、伝搬計算部12による電波のパスの計算例を示す図である。具体的には、図3に示すように地図上に構造物、送信点2、及び受信点4が配置される場合おいて、送信点2から受信点4までの電波のパスをレイトレース法により計算した例が示されている。レイトレース法では、送信点2から放射される電波は光線と見なされ、周囲の構造物での反射・回折・透過を経て受信点4に到達する光線の軌跡が求められる。光線の軌跡を求める方法は、例えば、イメージング法でもよいしローンチング法でもよい。光線の軌跡が送信点2から受信点4までの電波のパスであり、図4に示す例ではP1からP5までの5つのパスが計算されている。このうちパスP1は、送信点2から受信点4に直接到達する直接波のパスである。 FIG. 4 is a diagram showing an example of calculation of a radio wave path by the propagation calculation unit 12. Specifically, when structures, transmitting points 2, and receiving points 4 are arranged on a map as shown in Figure 3, the radio wave path from transmitting point 2 to receiving point 4 is determined by the ray tracing method. A calculated example is shown. In the ray tracing method, radio waves emitted from the transmitting point 2 are regarded as light rays, and the trajectory of the light rays reaching the receiving point 4 after being reflected, diffracted, and transmitted by surrounding structures is determined. The method for determining the trajectory of the light ray may be, for example, an imaging method or a launching method. The trajectory of the light beam is the path of the radio wave from the transmission point 2 to the reception point 4, and in the example shown in FIG. 4, five paths from P1 to P5 are calculated. Among these, the path P1 is a direct wave path that directly reaches the receiving point 4 from the transmitting point 2.

伝搬計算部12は、評価地点ごとに、例えば、通信距離、パス数、各パスの受信電力、各パスの到来時間(または遅延)、直接波を受信できる受信点数、チャネル応答を計算する。図4に示す例では、通信距離は送信点2から受信点4までの距離であり、パス数はP1からP5までの5つである。各パスの受信電力は、受信点4での電波の振幅及び位相等に基づいて公知の方法で計算される。各パスの到来時間はパスの経路長に基づいて公知の方法で計算される。図5は、伝搬計算部12による各パスの受信電力及び到来時間の計算例を示す図である。また、チャネル応答は受信点4及び送信点2の各アンテナの構成と受信点4での各パスの電波の振幅及び位相等に基づいて公知の方法で計算される。 The propagation calculation unit 12 calculates, for each evaluation point, for example, the communication distance, the number of paths, the received power of each path, the arrival time (or delay) of each path, the number of reception points that can receive direct waves, and the channel response. In the example shown in FIG. 4, the communication distance is the distance from the transmission point 2 to the reception point 4, and the number of paths is five from P1 to P5. The received power of each path is calculated by a known method based on the amplitude and phase of the radio wave at the receiving point 4. The arrival time of each path is calculated in a known manner based on the path length of the path. FIG. 5 is a diagram showing an example of calculation of the received power and arrival time of each path by the propagation calculation unit 12. Further, the channel response is calculated by a known method based on the configuration of each antenna at the receiving point 4 and the transmitting point 2, and the amplitude and phase of the radio waves of each path at the receiving point 4.

精度推定部13は、伝搬計算部12から出力された情報を用い、各評価地点における位置推定精度を算出する。位置推定精度は、例えば、以下に式(1)で示す一次多項式を用いて算出される。この一次多項式において位置推定精度は目的変数であり、F1、F2、F3、F4、及びF5は説明変数である。a、b、c、d、及びeは重み係数である。
(位置推定精度)=a*F1+b*F2+c*F3+d*F4+e*F5 ・・・式(1)
The accuracy estimation unit 13 uses the information output from the propagation calculation unit 12 to calculate the position estimation accuracy at each evaluation point. The position estimation accuracy is calculated using, for example, a first-order polynomial shown in equation (1) below. In this first-order polynomial, the position estimation accuracy is the objective variable, and F1, F2, F3, F4, and F5 are explanatory variables. a, b, c, d, and e are weighting factors.
(Position estimation accuracy) = a*F1+b*F2+c*F3+d*F4+e*F5...Formula (1)

上記の一次多項式におけるF1、F2、F3、F4、及びF5は、位置推定精度に影響する互いに独立な5つの因子に関係するパラメータである。F1は、通信距離に基づいて計算されるパラメータである。F2は、パス数に基づいて計算されるパラメータである。F3は、各パスの受信電力と到来時間(または遅延)とに基づいて計算されるパラメータである。F4は、直接波を受信できる受信点数に基づいて計算されるパラメータである。そして、F5は、チャネル応答に基づいて計算されるパラメータである。 F1, F2, F3, F4, and F5 in the above first-order polynomial are parameters related to five mutually independent factors that affect position estimation accuracy. F1 is a parameter calculated based on the communication distance. F2 is a parameter calculated based on the number of passes. F3 is a parameter calculated based on the received power and arrival time (or delay) of each path. F4 is a parameter calculated based on the number of reception points that can receive direct waves. And F5 is a parameter calculated based on the channel response.

図6は、精度推定部13の機能を説明するためのブロック図である。上記の各パラメータF1-F5を計算するため、精度推定部13は、図6に示すようにF1計算部131、
F2計算部132、F3計算部133、F4計算部134、及びF5計算部135を備える。
FIG. 6 is a block diagram for explaining the functions of the accuracy estimation section 13. In order to calculate each of the above parameters F1 to F5, the accuracy estimation unit 13 includes an F1 calculation unit 131,
It includes an F2 calculation section 132, an F3 calculation section 133, an F4 calculation section 134, and an F5 calculation section 135.

F1計算部131は、通信距離が大きいほど位置推定精度が劣化しやすいという特性を考慮してパラメータF1を計算する。F1計算部131には、予め通信距離の基準値が設定されている。通信距離の基準値とは、要求する位置推定精度が得られると見込まれる値であり、ユーザが任意に設定できる。例えば、位置推定誤差10mを許容できるとしたとき、位置推定誤差10mが満足できる通信距離が50mであれば、伝搬計算部12の出力である通信距離の基準値は50mとなる。F1計算部131は、例えば、伝搬計算部12から出力された通信距離の計算値に対する通信距離の基準値の比(基準値/計算値)をパラメータF1として算出する。パラメータF1は、位置推定精度に影響する因子である通信距離をその基準値で無次元化した無次元量である。 The F1 calculation unit 131 calculates the parameter F1 in consideration of the characteristic that the longer the communication distance is, the more likely the position estimation accuracy is to deteriorate. A reference value of communication distance is set in advance in the F1 calculation unit 131. The communication distance reference value is a value that is expected to provide the required position estimation accuracy, and can be set arbitrarily by the user. For example, if it is assumed that a position estimation error of 10 m is acceptable, and the communication distance that satisfies the position estimation error of 10 m is 50 m, the reference value of the communication distance that is the output of the propagation calculation unit 12 is 50 m. The F1 calculation unit 131 calculates, for example, the ratio of the communication distance reference value to the communication distance calculation value output from the propagation calculation unit 12 (reference value/calculation value) as the parameter F1. The parameter F1 is a dimensionless quantity obtained by making the communication distance, which is a factor that influences position estimation accuracy, dimensionless using its reference value.

F2計算部132は、パス数が大きいほど位置推定精度が劣化しやすいという特性を考慮してパラメータF2を計算する。F2計算部132には、予めパス数の基準値が設定されている。パス数の基準値とは、要求する位置推定精度が得られると見込まれる値であり、ユーザが任意に設定できる。F2計算部132は、例えば、伝搬計算部12から出力されたパス数の計算値に対するパス数の基準値の比(基準値/計算値)をパラメータF2として算出する。パラメータF2は、位置推定精度に影響する因子であるパス数をその基準値で無次元化した無次元量である。 The F2 calculation unit 132 calculates the parameter F2 in consideration of the characteristic that the position estimation accuracy is more likely to deteriorate as the number of passes increases. A reference value for the number of passes is set in advance in the F2 calculation unit 132. The reference value for the number of passes is a value that is expected to provide the required position estimation accuracy, and can be set arbitrarily by the user. The F2 calculating unit 132 calculates, for example, the ratio of the reference value of the number of paths to the calculated value of the number of paths output from the propagation calculating unit 12 (reference value/calculated value) as the parameter F2. The parameter F2 is a dimensionless quantity obtained by making the number of passes, which is a factor that affects position estimation accuracy, dimensionless using its reference value.

F3計算部133は、複数のパスの受信電力が同程度であるとき、パスの到来時間(遅延)が大きいほど位置推定精度が劣化しやすいという特性を考慮してパラメータF3を計算する。具体的には、F3計算部133は、各パスのうち受信電力が同程度となるパスを抽出する。そして、同程度の受信電力を有するパスのうちで到来時間差の最大値(最大到来時間差)を計算する。図5に示す例では、パスP1、P2、及びP3が同程度の受信電力を有するパスであり、パスP1の到来時間とパスP3の到来時間との差が最大到来時間差である。F3計算部133には、予め最大到来時間差の基準値が設定されている。最大到来時間差の基準値とは、要求する位置推定精度が得られると見込まれる値であり、ユーザが任意に設定できる。F3計算部133は、例えば、上記のように計算された最大到来時間差の計算値に対する最大到来時間差の基準値の比(基準値/計算値)をパラメータF3として算出する。パラメータF3は、位置推定精度に影響する因子である最大到来時間差をその基準値で無次元化した無次元量である。 The F3 calculation unit 133 calculates the parameter F3 in consideration of the characteristic that when the received power of a plurality of paths is about the same, the longer the arrival time (delay) of a path is, the more likely the position estimation accuracy is to deteriorate. Specifically, the F3 calculation unit 133 extracts paths with the same received power from among the paths. Then, the maximum value of the arrival time difference (maximum arrival time difference) among the paths having the same received power is calculated. In the example shown in FIG. 5, paths P1, P2, and P3 have the same received power, and the difference between the arrival time of path P1 and the arrival time of path P3 is the maximum arrival time difference. A reference value for the maximum arrival time difference is set in advance in the F3 calculation unit 133. The reference value of the maximum arrival time difference is a value that is expected to provide the required position estimation accuracy, and can be set arbitrarily by the user. The F3 calculating unit 133 calculates, for example, the ratio (reference value/calculated value) of the reference value of the maximum arrival time difference to the calculated value of the maximum arrival time difference calculated as described above (reference value/calculated value) as the parameter F3. The parameter F3 is a dimensionless quantity obtained by making the maximum arrival time difference, which is a factor that affects position estimation accuracy, dimensionless using its reference value.

F4計算部134は、直接波を受信できる受信点数が少ないほど位置推定精度が劣化しやすいという特性を考慮してパラメータF4を計算する。F4計算部134には、予め受信点数の基準値が設定されている。受信点数の基準値とは、要求する位置推定精度が得られると見込まれる値であり、ユーザが任意に設定できる。F4計算部134は、例えば、伝搬計算部12から出力された受信点数の計算値に対する受信点数の基準値の比(基準値/計算値)をパラメータF4として算出する。パラメータF4は、位置推定精度に影響する因子である受信点数をその基準値で無次元化した無次元量である。 The F4 calculation unit 134 calculates the parameter F4 in consideration of the characteristic that the position estimation accuracy is more likely to deteriorate as the number of reception points that can receive direct waves is smaller. A reference value for the number of reception points is set in advance in the F4 calculation unit 134. The reference value of the number of reception points is a value that is expected to provide the required position estimation accuracy, and can be set arbitrarily by the user. The F4 calculation unit 134 calculates, for example, the ratio of the reference value of the number of reception points to the calculated value of the number of reception points output from the propagation calculation unit 12 (reference value/calculation value) as the parameter F4. The parameter F4 is a dimensionless quantity obtained by making the number of reception points, which is a factor that influences position estimation accuracy, dimensionless using its reference value.

次に、F5計算部135によるチャネル応答に基づいたパラメータF5の計算について説明する。チャネル応答の値が位置推定精度に影響するのは、例えば、機械学習を用いて位置情報に紐づけられたチャネル情報から位置推定が行われる場合である。相関の高いチャネル応答となる地点が複数存在した場合、位置推定に誤差が生じやすい。そこで、F5計算部135は、図7にフローチャートで示す手順にてパラメータF5を計算する。 Next, the calculation of the parameter F5 based on the channel response by the F5 calculation unit 135 will be explained. The value of the channel response affects the position estimation accuracy, for example, when position estimation is performed from channel information linked to position information using machine learning. If there are multiple points with highly correlated channel responses, errors are likely to occur in position estimation. Therefore, the F5 calculation unit 135 calculates the parameter F5 according to the procedure shown in the flowchart in FIG.

ステップS1では、F5計算部135は、位置推定の精度を求めたい場所と比較場所1地点の計2地点を選択する。ステップS2では、F5計算部135は、各地点で得られているチャネル応答を用いて、チャネル相関値を計算する。F5計算部135は、ステップS1とステップS2の処理を繰り返す。すなわち、ステップS1において異なる比較場所を選択し、ステップS2において選択した比較場所についてチャネル相関値を計算する。全ての比較場所についてのチャネル相関値の計算が完了すると、F5計算部135は、次にステップS3の処理を行う。 In step S1, the F5 calculation unit 135 selects a total of two points: a place where the accuracy of position estimation is desired and one comparison point. In step S2, the F5 calculation unit 135 calculates a channel correlation value using the channel response obtained at each point. The F5 calculation unit 135 repeats the processing of step S1 and step S2. That is, a different comparison location is selected in step S1, and a channel correlation value is calculated for the selected comparison location in step S2. When the calculation of channel correlation values for all comparison locations is completed, the F5 calculation unit 135 next performs the process of step S3.

ステップS3では、F5計算部135は、全ての比較場所の中からチャネル相関値が予め設定した閾値以上である場所を抽出する。続いて、ステップS4では、F5計算部135は、ステップS3で抽出された全ての場所について、位置推定の精度を求めたい場所からの距離を計算する。そして、計算された距離の中から最大距離を求める。 In step S3, the F5 calculation unit 135 extracts locations where the channel correlation value is greater than or equal to a preset threshold value from among all comparison locations. Subsequently, in step S4, the F5 calculation unit 135 calculates the distance from the location for which the accuracy of position estimation is desired for all the locations extracted in step S3. Then, the maximum distance is determined from among the calculated distances.

ここで、図8はステップS3及びS4で行われる処理を具体的に説明するための図である。図8には、地図上における位置推定の精度を求めたい場所TLと、全ての比較場所L1-L14の位置が示されている。図8に示す例では、比較場所L1-L14のうちチャネル相関値が閾値以上である場所は、L2、L5、L8、及びL11の計5地点である。そして、位置推定の精度を求めたい場所TLから各場所L2、L5、L8、L11までの距離D2、D5、D8、D11のうち最大距離は、場所L5までの距離D5である。この例の場合、ステップS4では、距離D5が最大距離として求められる。 Here, FIG. 8 is a diagram for specifically explaining the processing performed in steps S3 and S4. FIG. 8 shows a location TL on the map for which the accuracy of position estimation is desired and the locations of all comparison locations L1 to L14. In the example shown in FIG. 8, among the comparison locations L1 to L14, there are a total of five locations where the channel correlation value is greater than or equal to the threshold value: L2, L5, L8, and L11. Of the distances D2, D5, D8, and D11 from the location TL for which the accuracy of position estimation is desired to each location L2, L5, L8, and L11, the maximum distance is the distance D5 to the location L5. In this example, the distance D5 is determined as the maximum distance in step S4.

ステップS5では、F5計算部135は、許容できる最大距離を基準値として設定する。基準値の設定は毎回行われてもよいし、予め設定された基準値が用いられてもよい。ステップF6では、F5計算部135は、例えば、上記のように計算された最大距離の計算値に対する最大距離の基準値の比(基準値/計算値)をパラメータF5して算出する。パラメータF5は、位置推定精度に影響する因子である類似性の高い場所までの最大距離をその基準値で無次元化した無次元量である。 In step S5, the F5 calculation unit 135 sets the maximum allowable distance as a reference value. The reference value may be set every time, or a preset reference value may be used. In step F6, the F5 calculation unit 135 calculates, for example, the ratio of the reference value of the maximum distance to the calculated value of the maximum distance calculated as described above (reference value/calculated value) as a parameter F5. The parameter F5 is a dimensionless quantity obtained by making the maximum distance to a highly similar place, which is a factor that affects position estimation accuracy, dimensionless using its reference value.

再び図6に戻り、精度推定部13の機能について説明を続ける。精度推定部13は、位置推定精度計算部130を備える。F1計算部131、F2計算部132、F3計算部133、F4計算部134、及びF5計算部135で計算された各パラメータF1、F2、F3、F4、F5は、位置推定精度計算部130に入力される。位置推定精度計算部130は、上記の式(1)で表される一次多項式を用いて位置推定精度を算出する。 Returning to FIG. 6 again, the description of the function of the accuracy estimation unit 13 will be continued. The accuracy estimation unit 13 includes a position estimation accuracy calculation unit 130. The parameters F1, F2, F3, F4, and F5 calculated by the F1 calculation unit 131, F2 calculation unit 132, F3 calculation unit 133, F4 calculation unit 134, and F5 calculation unit 135 are input to the position estimation accuracy calculation unit 130. be done. The position estimation accuracy calculation unit 130 calculates the position estimation accuracy using the first-order polynomial expressed by the above equation (1).

なお、一次多項式を構成する重み係数a、b、c、d、eはユーザが任意に設定することができる。一次多項式で表される関係は採用する位置推定方法により異なるため、位置推定方法に応じて異なる値を重み係数a、b、c、d、eに設定してもよい。また、想定する無線環境において支配的になる特定誤差の要因が異なる場合には、その要因が多く反映されるように重み係数a、b、c、d、eを設定してもよい。位置推定方法により考慮しない因子に対応する重み係数には0が設定される。例えば、チャネル応答を用いない場合はe=0とされる。 Note that the weighting coefficients a, b, c, d, and e constituting the first-order polynomial can be arbitrarily set by the user. Since the relationship expressed by the first-order polynomial differs depending on the position estimation method employed, different values may be set for the weighting coefficients a, b, c, d, and e depending on the position estimation method. Further, if the factors of the specific error that are dominant in the assumed radio environment are different, the weighting coefficients a, b, c, d, and e may be set so as to reflect many of the factors. Weighting coefficients corresponding to factors not considered by the position estimation method are set to 0. For example, if no channel response is used, e=0.

以上のように、位置推定精度可視化装置10は、伝搬計算部12と精度推定部13とにおいて、無線通信の環境条件に基づき位置推定精度に影響する互いに独立な複数の因子についてそれぞれの数値を算出する処理を行う。位置推定精度可視化装置10は、次に、精度推定部13において、複数の因子のそれぞれの基準値を用いて複数の因子のそれぞれの数値を無次元量に変換する処理を行う。そして、位置推定精度可視化装置10は、精度推定部13において、複数の因子のそれぞれの無次元量を説明変数とし位置推定精度を目的変数とする一次多項式を用いて位置推定精度を計算する処理を行う。 As described above, the position estimation accuracy visualization device 10 uses the propagation calculation unit 12 and the accuracy estimation unit 13 to calculate respective numerical values for a plurality of mutually independent factors that affect position estimation accuracy based on the environmental conditions of wireless communication. Perform the processing to do. Next, in the position estimation accuracy visualization device 10, the accuracy estimation unit 13 performs a process of converting each numerical value of the plurality of factors into a dimensionless quantity using the reference value of each of the plurality of factors. Then, the position estimation accuracy visualization device 10 causes the accuracy estimation unit 13 to perform a process of calculating the position estimation accuracy using a first-order polynomial in which the dimensionless quantities of the plurality of factors are explanatory variables and the position estimation accuracy is the objective variable. conduct.

再び図2に戻り、位置推定精度可視化装置10の機能について説明を続ける。精度推定部13で計算された位置推定精度は表示部14に出力される。表示部14は、精度推定部13から出力された位置推定精度の計算結果をグラフまたは地図として表示する。例えば、入力部11にて入力した地図情報を参照して地図を表示するとともに、地図上の対応する位置の位置推定精度を表示する。表示方法としては、精度を数値として表示してもよいし、色で表示してもよい。 Returning to FIG. 2 again, the description of the functions of the position estimation accuracy visualization device 10 will be continued. The position estimation accuracy calculated by the accuracy estimation section 13 is output to the display section 14. The display unit 14 displays the calculation result of the position estimation accuracy output from the accuracy estimation unit 13 as a graph or a map. For example, a map is displayed with reference to map information input through the input unit 11, and the position estimation accuracy of a corresponding position on the map is displayed. As a display method, accuracy may be displayed as a numerical value or may be displayed in color.

図9は、表示部14による表示例を示す図である。図9に示すように、構造物が配置された地図上において、位置推定の精度が高いエリアと、位置推定精度が低いエリアとを色分けして表示するようにしてもよい。例えば、位置推定精度が高ければ赤、低ければ青のようエリア毎に色分けして表示してもよい。また、位置推定精度が高いか低いかで2色に色分けするのではなく、位置推定精度の値に応じて多段階にエリア毎に色分けするようにしてもよい。 FIG. 9 is a diagram showing an example of display on the display unit 14. As shown in FIG. 9, on a map where structures are placed, areas with high position estimation accuracy and areas with low position estimation accuracy may be displayed in different colors. For example, the area may be displayed in different colors, such as red if the position estimation accuracy is high and blue if it is low. Moreover, instead of color-coding into two colors depending on whether the position estimation accuracy is high or low, it is also possible to color-code each area in multiple stages according to the value of the position estimation accuracy.

3.第2実施形態の位置推定精度可視化装置
上記の式(1)で表される一次多項式の重み係数a、b、c、d、eの設定において、予め測定した位置推定の実測値を活用してもよい。重み係数の設定に位置推定の実測値を活用するため、第2実施形態の位置推定精度可視化装置10は、図10に示すように構成される。第2実施形態の位置推定精度可視化装置10は、その機能別に、入力部11、伝搬計算部12、精度推定部13、表示部14、誤差計算部15、及び重み推定部16から構成される。入力部11、伝搬計算部12、精度推定部13、及び表示部14は第1実施形態と共通であるので、以下では誤差計算部15と重み推定部16について説明する。
3. Position estimation accuracy visualization device of second embodiment In setting the weighting coefficients a, b, c, d, and e of the first-order polynomial expressed by the above equation (1), the actual values of position estimation measured in advance are used. Good too. In order to utilize actual measured values of position estimation for setting weighting coefficients, the position estimation accuracy visualization device 10 of the second embodiment is configured as shown in FIG. 10. The position estimation accuracy visualization device 10 of the second embodiment includes an input section 11, a propagation calculation section 12, an accuracy estimation section 13, a display section 14, an error calculation section 15, and a weight estimation section 16 according to their functions. Since the input unit 11, propagation calculation unit 12, accuracy estimation unit 13, and display unit 14 are the same as those in the first embodiment, the error calculation unit 15 and weight estimation unit 16 will be explained below.

誤差計算部15は、実測により得た位置推定情報(実測値)と、それに対応する正しい位置情報(正解値)の入力を受け付ける。誤差計算部15は、入力された実測値と正解値から誤差を計算する。例えば、実測値と正解値との間の距離が誤差として計算される。誤差計算部15は、誤差の計算結果を重み推定部16に出力する。 The error calculation unit 15 receives input of position estimation information (actual measurement value) obtained through actual measurement and correct position information (correct value) corresponding thereto. The error calculation unit 15 calculates an error from the input actual measurement value and correct value. For example, the distance between the measured value and the correct value is calculated as the error. The error calculation unit 15 outputs the error calculation result to the weight estimation unit 16.

重み推定部16は、誤差計算部15から出力された様々な場所における誤差の計算結果を用いて、実測により得られた位置推定精度に近づくように、上記の式(1)で表される一次多項式の重み係数a、b、c、d、eを調整する。その際、上記の式(1)の左辺の位置推定精度は以下の式(2)によって計算してもよい。なお、式(2)の右辺における許容誤差は、要求する位置推定精度が得られると見込まれる値であり、ユーザが任意に設定できる。
(位置推定精度)=(許容誤差)/(誤差計算部から出力された誤差) ・・・式(2)
The weight estimation unit 16 uses the calculation results of errors at various locations output from the error calculation unit 15 to calculate the linear estimation accuracy expressed by the above equation (1) so as to approach the position estimation accuracy obtained by actual measurement. Adjust the weighting coefficients a, b, c, d, and e of the polynomial. In this case, the position estimation accuracy on the left side of the above equation (1) may be calculated using the following equation (2). Note that the allowable error on the right side of equation (2) is a value that is expected to provide the required position estimation accuracy, and can be set arbitrarily by the user.
(Position estimation accuracy) = (Tolerance error) / (Error output from the error calculation unit) ...Formula (2)

このようにして重み係数a、b、c、d、eが算出されることにより、用いる位置推定方法及び位置推定を行う場所に対応した重み係数a、b、c、d、eを得ることができる。得られた重み係数a、b、c、d、eは精度推定部13に出力され、式(1)で表される一次多項式に適用される。 By calculating the weighting coefficients a, b, c, d, and e in this way, it is possible to obtain the weighting coefficients a, b, c, d, and e corresponding to the position estimation method used and the location where the position estimation is performed. can. The obtained weighting coefficients a, b, c, d, and e are output to the accuracy estimating section 13 and applied to the first-order polynomial expressed by equation (1).

第2実施形態においては、予め位置推定を行う場所において位置推定の実測値の測定が行われ、その測定結果が重み係数a、b、c、d、eの設定に反映される。その結果、第2実施形態によれば、位置推定の実測値は重み係数を算出するために用いられ、位置推定精度は伝搬計算部において算出された結果が用いられるため、位置推定の実測値を測定していない場所を含めて精度を確認することができる。 In the second embodiment, actual position estimation values are measured in advance at a location where position estimation is performed, and the measurement results are reflected in the settings of weighting coefficients a, b, c, d, and e. As a result, according to the second embodiment, the actual measured value of position estimation is used to calculate the weighting coefficient, and the result calculated in the propagation calculation unit is used for the position estimation accuracy. Accuracy can be confirmed including locations where measurements have not been made.

4.位置推定精度可視化装置のハードウェア構成
本開示の各実施形態の位置推定精度可視化装置10のハードウェア構成について図11を参照して説明する。図11は、位置推定精度可視化装置10のハードウェア構成の例を示すブロック図である。
4. Hardware Configuration of Position Estimation Accuracy Visualization Device The hardware configuration of the position estimation accuracy visualization device 10 of each embodiment of the present disclosure will be described with reference to FIG. 11. FIG. 11 is a block diagram showing an example of the hardware configuration of the position estimation accuracy visualization device 10.

図2及び図10に示す位置推定精度可視化装置10の各機能は、例えば、コンピュータ100と入力装置105と表示装置106とにより実現される。コンピュータ100は、1又は複数のプロセッサ101(以下、単にプロセッサ101と呼ぶ)と、プロセッサ101に結合された1又は複数のメモリ102(以下、単にメモリ102と呼ぶ)とを含んでいる。メモリ102は主記憶装置と補助記憶装置とを含む。メモリ102には、プロセッサ101で実行可能な1又は複数のプログラム103(以下、単にプログラム103と呼ぶ)とそれに関連する種々の情報とが記憶されている。 Each function of the position estimation accuracy visualization device 10 shown in FIGS. 2 and 10 is realized by, for example, the computer 100, the input device 105, and the display device 106. Computer 100 includes one or more processors 101 (hereinafter simply referred to as processors 101) and one or more memories 102 (hereinafter simply referred to as memory 102) coupled to processors 101. Memory 102 includes a main storage device and an auxiliary storage device. The memory 102 stores one or more programs 103 (hereinafter simply referred to as programs 103) executable by the processor 101 and various information related thereto.

プロセッサ101がプログラム103を実行することにより、プロセッサ101による各種処理が実現される。図2に示す入力部11、伝搬計算部12、精度推定部13、及び表示部14の各機能は、プログラム103がプロセッサ101により実行されることで実現される。図10示す入力部11、伝搬計算部12、精度推定部13、表示部14、誤差計算部15、及び重み推定部16の各機能も、プログラム103がプロセッサ101により実行されることで実現される。プログラム103は、主記憶装置に記憶されることもできるし、補助記憶装置であるコンピュータ読み取り可能な記録媒体に記憶されることもできる。 When the processor 101 executes the program 103, various processing by the processor 101 is realized. The functions of the input section 11, propagation calculation section 12, precision estimation section 13, and display section 14 shown in FIG. 2 are realized by the program 103 being executed by the processor 101. The functions of the input section 11, propagation calculation section 12, accuracy estimation section 13, display section 14, error calculation section 15, and weight estimation section 16 shown in FIG. 10 are also realized by the program 103 being executed by the processor 101. . The program 103 can be stored in the main storage device or in a computer-readable recording medium that is an auxiliary storage device.

入力装置105は、キーボード、マウス等の入力デバイスを含む。また、入力装置105の機能には、USBメモリ等を介してデータをコンピュータ100に入力する機能、L
AN経由でデータをコンピュータ100に入力する機能、或いは、インターネットからデータをダウンロードする機能も含まれる。表示装置106は、ユーザによって視認可能な表示画面を有するディスプレイである。
The input device 105 includes input devices such as a keyboard and a mouse. The functions of the input device 105 include a function of inputting data into the computer 100 via a USB memory, etc.
It also includes a function to input data into the computer 100 via the AN or a function to download data from the Internet. The display device 106 is a display having a display screen that can be viewed by the user.

5.その他
上記実施形態は、本開示の要旨を逸脱しない範囲で種々に変形して実施することができる。すなわち、上記実施形態において各要素の個数、数量、量、範囲などの数に言及されている場合、特に明示した場合や原理的に明らかにその数に特定される場合を除いて、その言及した数に、本開示に係る技術が限定されるものではない。また、上記実施形態において説明する構造等は、特に明示した場合や明らかに原理的にそれに特定される場合を除いて、本開示に係る技術に必ずしも必須のものではない。
5. Others The embodiments described above can be modified and implemented in various ways without departing from the gist of the present disclosure. In other words, in the above embodiments, when the number, amount, amount, range, etc. of each element is mentioned, unless it is specifically specified or it is clearly specified by the number in principle, the mentioned The technology according to the present disclosure is not limited to the number. Furthermore, the structures and the like described in the above embodiments are not necessarily essential to the technology according to the present disclosure, unless specifically specified or clearly specified in principle.

2 送信点(移動局)
4 受信点(基地局)
6 位置推定システム
10 位置推定精度可視化装置
11 入力部
12 伝搬計算部
13 精度推定部
14 表示部
15 誤差計算部
16 重み推定部
100 コンピュータ
101 プロセッサ
102 メモリ
103 プログラム
105 入力装置
106 表示装置
P1、P2、P3、P4、P5 パス
2 Transmission point (mobile station)
4 Receiving point (base station)
6 Position estimation system 10 Position estimation accuracy visualization device 11 Input section 12 Propagation calculation section 13 Accuracy estimation section 14 Display section 15 Error calculation section 16 Weight estimation section 100 Computer 101 Processor 102 Memory 103 Program 105 Input device 106 Display devices P1, P2, P3, P4, P5 pass

Claims (8)

無線通信による位置推定の精度を可視化する方法であって、
前記無線通信の環境条件に基づき前記精度に影響する互いに独立な複数の因子についてそれぞれの数値を算出するステップと、
前記複数の因子のそれぞれの基準値を用いて前記複数の因子のそれぞれの数値を無次元量に変換するステップと、
前記複数の因子のそれぞれの無次元量を説明変数とし前記精度を目的変数とする一次多項式を用いて前記精度を計算するステップと、
前記精度の計算結果を表示するステップと、
を含むことを特徴とする方法。
A method for visualizing the accuracy of position estimation by wireless communication, the method comprising:
calculating respective numerical values for a plurality of mutually independent factors that influence the accuracy based on the environmental conditions of the wireless communication;
converting each numerical value of the plurality of factors into a dimensionless quantity using the reference value of each of the plurality of factors;
calculating the accuracy using a first-order polynomial with dimensionless quantities of each of the plurality of factors as explanatory variables and the accuracy as an objective variable;
Displaying the calculation result of the accuracy;
A method characterized by comprising:
請求項1に記載の方法において、
前記精度を計算する前記ステップは、前記位置推定が行われるエリア毎に前記精度を計算するステップを含み、
前記計算結果を表示する前記ステップは、地図上の前記エリアに対応する位置に前記計算結果を表示するステップを含む
ことを特徴とする方法。
The method according to claim 1,
The step of calculating the accuracy includes calculating the accuracy for each area in which the position estimation is performed,
The method characterized in that the step of displaying the calculation result includes the step of displaying the calculation result at a position corresponding to the area on a map.
請求項1又は2に記載の方法において、
前記位置推定の実測値と正解値との誤差に基づいて前記一次多項式の重み係数を調整するステップをさらに含む
ことを特徴とする方法。
The method according to claim 1 or 2,
The method further comprises the step of adjusting a weighting coefficient of the first-order polynomial based on an error between the actual measured value of the position estimation and the correct value.
請求項1乃至3の何れか1項に記載の方法において、
前記複数の因子は電波の伝搬路に依存する因子を含み、
前記数値を算出するステップは、前記無線通信の前記環境条件に基づき前記伝搬路を計算するステップを含む
ことを特徴とする方法。
The method according to any one of claims 1 to 3,
The plurality of factors include factors that depend on the propagation path of the radio wave,
A method characterized in that the step of calculating the numerical value includes the step of calculating the propagation path based on the environmental conditions of the wireless communication.
無線通信による位置推定の精度を可視化する装置であって、
1つ又は複数のプログラムを記憶した1つ又は複数のメモリと、
前記1つ又は複数のメモリと結合された1つ又は複数のプロセッサと、を備え、
前記1つ又は複数のプロセッサは、前記1つ又は複数のプログラムの実行時、
前記無線通信の環境条件に基づき前記精度に影響する互いに独立な複数の因子についてそれぞれの数値を算出する処理と、
前記複数の因子のそれぞれの基準値を用いて前記複数の因子のそれぞれの数値を無次元量に変換する処理と、
前記複数の因子のそれぞれの無次元量を説明変数とし前記精度を目的変数とする一次多項式を用いて前記精度を計算する処理と、
前記精度の計算結果を表示する処理と、を実行する
を含むことを特徴とする装置。
A device that visualizes the accuracy of position estimation by wireless communication,
one or more memories storing one or more programs;
one or more processors coupled with the one or more memories,
When the one or more processors execute the one or more programs,
a process of calculating respective numerical values for a plurality of mutually independent factors that influence the accuracy based on the environmental conditions of the wireless communication;
a process of converting each numerical value of the plurality of factors into a dimensionless quantity using a reference value of each of the plurality of factors;
a process of calculating the accuracy using a first-order polynomial with dimensionless quantities of each of the plurality of factors as explanatory variables and the accuracy as an objective variable;
An apparatus characterized in that the apparatus comprises the steps of: displaying the calculation result of the accuracy;
請求項5に記載の装置において、
前記精度を計算する前記処理は、前記位置推定が行われるエリア毎に前記精度を計算する処理を含み、
前記計算結果を表示する前記処理、地図上の前記エリアに対応する位置に前記計算結果を表示する処理を含む
ことを特徴とする装置。
The apparatus according to claim 5,
The process of calculating the accuracy includes a process of calculating the accuracy for each area in which the position estimation is performed,
An apparatus characterized by comprising: the process of displaying the calculation result; and the process of displaying the calculation result at a position corresponding to the area on a map.
請求項5又は6に記載の装置において、
前記1つ又は複数のプロセッサは、前記1つ又は複数のプログラムの実行時、
前記位置推定の実測値と正解値との誤差に基づいて前記一次多項式の重み係数を調整する処理をさらに実行する
ことを特徴とする装置。
The device according to claim 5 or 6,
When the one or more processors execute the one or more programs,
An apparatus characterized in that the apparatus further executes a process of adjusting a weighting coefficient of the first-order polynomial based on an error between an actual measured value of the position estimation and a correct value.
請求項5乃至7の何れか1項に記載の装置において、
前記複数の因子は電波の伝搬路に依存する因子を含み、
前記数値を算出する処理は、前記無線通信の前記環境条件に基づき前記伝搬路を計算する処理を含む
ことを特徴とする装置。
The device according to any one of claims 5 to 7,
The plurality of factors include factors that depend on the propagation path of the radio wave,
An apparatus characterized in that the process of calculating the numerical value includes a process of calculating the propagation path based on the environmental conditions of the wireless communication.
JP2022576909A 2021-01-22 2021-01-22 Method and device for visualizing accuracy of position estimation by wireless communication Active JP7420289B2 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2021/002284 WO2022157936A1 (en) 2021-01-22 2021-01-22 Method and device for visualizing accuracy of position estimation using wireless communication

Publications (2)

Publication Number Publication Date
JPWO2022157936A1 JPWO2022157936A1 (en) 2022-07-28
JP7420289B2 true JP7420289B2 (en) 2024-01-23

Family

ID=82549570

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022576909A Active JP7420289B2 (en) 2021-01-22 2021-01-22 Method and device for visualizing accuracy of position estimation by wireless communication

Country Status (2)

Country Link
JP (1) JP7420289B2 (en)
WO (1) WO2022157936A1 (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010509571A (en) 2006-11-07 2010-03-25 スカイフック ワイヤレス,インク. System and method for estimating positioning error in a positioning system based on wireless LAN
JP2010130220A (en) 2008-11-26 2010-06-10 Hitachi-Ge Nuclear Energy Ltd System for assisting distribution of radio communication apparatus
JP2010186393A (en) 2009-02-13 2010-08-26 Nissan Motor Co Ltd Tracking controller and tracking control method
JP2010230501A (en) 2009-03-27 2010-10-14 Kddi Corp Position information acquiring device, position information acquisition system, position information acquisition method, and program
US20170074659A1 (en) 2015-09-16 2017-03-16 Here Global B.V. Method and apparatus for providing a location data error map
JP2019203823A (en) 2018-05-24 2019-11-28 株式会社デンソー Travel plan generator, travel plan generation method, and control program
WO2020066153A1 (en) 2018-09-26 2020-04-02 Necソリューションイノベータ株式会社 Information processing device, information processing system, information processing method, and non-temporary computer-readable medium

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010509571A (en) 2006-11-07 2010-03-25 スカイフック ワイヤレス,インク. System and method for estimating positioning error in a positioning system based on wireless LAN
JP2010130220A (en) 2008-11-26 2010-06-10 Hitachi-Ge Nuclear Energy Ltd System for assisting distribution of radio communication apparatus
JP2010186393A (en) 2009-02-13 2010-08-26 Nissan Motor Co Ltd Tracking controller and tracking control method
JP2010230501A (en) 2009-03-27 2010-10-14 Kddi Corp Position information acquiring device, position information acquisition system, position information acquisition method, and program
US20170074659A1 (en) 2015-09-16 2017-03-16 Here Global B.V. Method and apparatus for providing a location data error map
JP2019203823A (en) 2018-05-24 2019-11-28 株式会社デンソー Travel plan generator, travel plan generation method, and control program
WO2020066153A1 (en) 2018-09-26 2020-04-02 Necソリューションイノベータ株式会社 Information processing device, information processing system, information processing method, and non-temporary computer-readable medium

Also Published As

Publication number Publication date
WO2022157936A1 (en) 2022-07-28
JPWO2022157936A1 (en) 2022-07-28

Similar Documents

Publication Publication Date Title
KR102235544B1 (en) Channel delay determination method, positioning method, and related apparatus
US8509819B2 (en) Information processing apparatus and correction method
JP4207081B2 (en) Radio wave propagation characteristic estimation system, method and program
US11277214B2 (en) Radio environment estimation method and radio environment estimation apparatus
KR20160133716A (en) Positioning environment analysis apparatus, positioning performance projection method and system of terminal using the same
US10715263B2 (en) Information processing apparatus, method and non-transitory computer-readable storage medium
JP5699544B2 (en) Radio wave propagation characteristic estimation system, radio wave propagation characteristic estimation method, and computer program
JP2013003084A (en) Positioning device and positioning method
CN112584420B (en) Signal measurement method, terminal and network side equipment
KR101219913B1 (en) Location estimation method for mobile node
JP7420289B2 (en) Method and device for visualizing accuracy of position estimation by wireless communication
JP2019074422A (en) Wireless receiver, method for wireless reception, and wireless system
US20220394653A1 (en) Method and apparatus of positioning for accomodating wireless-environment change
Ji et al. Channel estimation using spherical-wave model for indoor LoS and obstructed LoS scenarios
US20220377699A1 (en) Access point (ap) placement using fine time measurement (ftm)
JP7082879B2 (en) Interference power estimation device and interference power estimation method
JP7130943B2 (en) Radio wave environment estimation device and radio wave environment estimation method
CN110708702A (en) Method and device for determining signal propagation model
JP2021170738A (en) Radio wave intensity estimation device, position estimation system, radio wave intensity estimation method and program
JP6607190B2 (en) POSITION ESTIMATION DEVICE, POSITION ESTIMATION SYSTEM, METHOD, AND RECORDING MEDIUM
JP7272050B2 (en) Position estimation device, position estimation method and position estimation program
CN112703415B (en) Network positioning RF planner
ALMĂJANU et al. Radio coverage analysis for mobile communication networks using ics telecom
Sapumohotti et al. Wilocsim: Simulation testbed for wlan location fingerprinting systems
KR102577520B1 (en) Method and apparatus of positioning for accommodating wireless-environment change

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230511

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20231212

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20231225

R150 Certificate of patent or registration of utility model

Ref document number: 7420289

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150