JP6685654B2

JP6685654B2 - Position specifying support device, position specifying method, and position specifying program

Info

Publication number: JP6685654B2
Application number: JP2015104780A
Authority: JP
Inventors: 坂野　鋭; 鋭坂野; 佑樹石黒; 康平小島
Original assignee: NTT Data Corp
Current assignee: NTT Data Corp
Priority date: 2015-05-22
Filing date: 2015-05-22
Publication date: 2020-04-22
Anticipated expiration: 2035-05-22
Also published as: JP2016217945A

Description

本発明の実施形態は、位置特定支援装置、位置特定方法、および位置特定プログラムに関する。 Embodiments of the present invention relate to a position specifying support device, a position specifying method, and a position specifying program.

Ｗｉ-ＦｉやＢｌｕｅｔｏｏｔｈ（登録商標）等の無線通信規格を用いた電波発信源からの電波強度を測定することで、スマートフォン等の端末装置の位置を特定する技術が知られている。このような従来の技術では、予め電波発信源からの電波強度を測定しておき、測定した電波強度と測定時の位置とを対応づけたデータを保持しておき、運用時において新たに取得した電波強度と既知のデータとの関係に基づいて、運用時において観測される電波の受信源である端末装置の位置を特定する。 A technique is known in which the position of a terminal device such as a smartphone is specified by measuring the radio field intensity from a radio wave transmission source that uses a wireless communication standard such as Wi-Fi or Bluetooth (registered trademark). In such a conventional technique, the radio field intensity from the radio wave transmission source is measured in advance, and the data that associates the measured radio field intensity with the position at the time of measurement is held and newly acquired during operation. The position of the terminal device, which is the reception source of radio waves observed during operation, is specified based on the relationship between the radio field intensity and known data.

しかしながら、従来の技術では、参照するデータ量が多い場合、或いは取得する電波強度のサンプル数が多い場合等において、計算コストが増大する傾向にある。この結果、計算時の処理負荷が増大する場合があった。 However, in the conventional technique, the calculation cost tends to increase when the amount of data to be referred to is large or when the number of samples of the radio field intensity to be acquired is large. As a result, the processing load at the time of calculation may increase.

特開２００６−６５７０３号公報JP, 2006-65703, A

本発明は、このような事情に鑑みてなされたもので、その目的は、処理負荷を軽減することができる位置特定支援装置、位置特定方法、および位置特定プログラムを提供することである。 The present invention has been made in view of such circumstances, and an object thereof is to provide a position identification support device, a position identification method, and a position identification program that can reduce the processing load.

上述した課題を解決するために、本発明の一態様は、発信源から送信される電波を受信する端末装置から、端末装置により複数の測定時刻において測定された電波の強度に基づく情報を取得する取得部と、空間相関法に基づいて、取得部により取得された複数の測定時刻に係る電波の強度の相関行列を解くことで、測定時刻における電波の強度の安定成分を示す固有ベクトルを抽出する抽出部と、抽出部により安定成分として抽出された固有ベクトルと、予め知られた端末装置の位置を示す座標ベクトルとを対応付けて対応情報として記憶部に記憶させる記憶制御部と、記憶部に記憶された対応情報と、取得部により取得された電波の強度に基づく情報とに基づいて、電波の強度に基づく情報の取得元である端末装置の位置を特定する特定部と、を備え、特定部は、測定時刻よりも後の運用時刻において取得部によって電波の強度を示す情報が取得された場合、記憶部に対応情報として記憶された複数の前記固有ベクトルのうち、固有値が大きい上位所定数の固有ベクトルを抽出し、抽出した所定数の固有ベクトルのそれぞれと、運用時刻における電波の強度を示す電波ベクトルとの内積を算出し、固有ベクトルと電波ベクトルとの内積の総和が最大となる座標ベクトルが示す位置を、運用時刻における端末装置の位置として特定することを特徴とする。 In order to solve the above problems, one embodiment of the present invention obtains information based on the intensity of a radio wave measured by a terminal device at a plurality of measurement times from a terminal device that receives a radio wave transmitted from a transmission source. Extraction to extract the eigenvector indicating the stable component of the radio wave intensity at the measurement time by solving the correlation matrix of the radio wave intensity at the plurality of measurement times acquired by the acquisition unit and the spatial correlation method based on the spatial correlation method. and parts, and eigenvectors extracted as a stable component by the extraction unit, is stored and a storage control unit which associates the coordinate vector representing the position of the previously-known terminal device is stored in the storage unit as the corresponding information in the storage unit Based on the correspondence information and the information based on the strength of the radio wave acquired by the acquisition unit, the specifying unit that specifies the position of the terminal device from which the information based on the strength of the radio wave is acquired. , The specifying unit, when the information indicating the intensity of the radio wave is acquired by the acquiring unit at the operation time after the measurement time, the eigenvalue among the plurality of eigenvectors stored as the correspondence information in the storage unit is A large predetermined number of eigenvectors are extracted, and the inner product of each of the extracted predetermined number of eigenvectors and the radio wave vector indicating the strength of the radio wave at the operation time is calculated, and the sum of the inner products of the eigenvector and the radio wave vector is maximized. The position indicated by the coordinate vector is specified as the position of the terminal device at the operation time .

また、本発明の一態様は、上記の位置特定支援装置において、前記記憶部に記憶された対応情報と、前記取得部により取得された電波の強度に基づく情報とに基づいて、前記電波の強度の取得元である端末装置の位置を特定する特定部をさらに備えることを特徴とする。 Further, according to an aspect of the present invention, in the above-described position identification support device, the strength of the radio wave is based on the correspondence information stored in the storage unit and the information based on the strength of the radio wave acquired by the acquisition unit. Is further provided with a specifying unit that specifies the position of the terminal device that is the acquisition source.

また、本発明の一態様は、上記の位置特定支援装置において、前記抽出部は、前記相関行列を解くことで、前記安定成分を示す固有ベクトルを抽出し、前記記憶制御部は、前記抽出部により抽出された固有ベクトルを、予め知られた前記端末装置の位置に対応付けて対応情報として前記記憶部に記憶させることを特徴とする。 Further, one aspect of the present invention is the location assistance device described above, the extraction unit, by solving the prior Symbol-correlation matrix, to extract the eigenvectors indicating the stable component, the storage control unit, said extraction The eigenvector extracted by the unit is stored in the storage unit as correspondence information in association with a previously known position of the terminal device.

また、本発明の一態様は、上記の位置特定支援装置において、前記特定部は、前記記憶部に記憶された対応情報に含まれる固有ベクトルと、前記取得部により取得された電波の強度に基づく情報に含まれる所定のベクトルとの類似度を算出し、前記所定のベクトルとの類似度が最大となる前記固有ベクトルと、前記記憶部に記憶された対応情報とに基づいて、前記電波の強度に基づく情報の取得元である端末装置の位置を特定することを特徴とする。 Further, according to one aspect of the present invention, in the above-described position identification support device, the identification unit is based on an eigenvector included in the correspondence information stored in the storage unit and information based on the intensity of the radio wave acquired by the acquisition unit. Is calculated based on the strength of the radio wave, based on the eigenvector having the maximum similarity to the predetermined vector and the correspondence information stored in the storage unit. It is characterized in that the position of the terminal device from which the information is acquired is specified.

また、本発明の一態様は、上記の位置特定支援装置において、前記取得部は、発信源から送信される電波を受信する端末装置から、前記端末装置により複数の時刻において測定された電波の強度に基づく情報から、前記端末装置により抽出された前記電波の強度における安定成分を示す情報を取得し、前記特定部は、前記記憶部に記憶された対応情報と、前記取得部により取得された安定成分を示す情報とに基づいて、前記電波の強度に基づく情報の取得元である端末装置の位置を特定することを特徴とする。 Further, according to one aspect of the present invention, in the above-described position identification support device, the acquisition unit is a radio wave intensity measured at a plurality of times by the terminal device from a terminal device receiving a radio wave transmitted from a transmission source. the information that is based on, and acquires information indicating the stable component in the strength before Symbol waves extracted by the terminal device, the specifying unit, the corresponding information stored in the storage unit, which is acquired by the acquisition unit The position of the terminal device from which the information based on the intensity of the radio wave is acquired is specified based on the information indicating the stable component.

また、本発明の他の態様は、発信源から送信される電波を受信する端末装置から、前記端末装置により複数の測定時刻において測定された前記電波の強度に基づく情報を取得し、空間相関法に基づいて、前記取得した複数の測定時刻に係る電波の強度の相関行列を解くことで、前記測定時刻における前記電波の強度の安定成分を示す固有ベクトルを抽出し、前記安定成分として抽出した前記固有ベクトルと、予め知られた前記端末装置の位置を示す座標ベクトルとを対応付けて対応情報として記憶部に記憶させせ、前記記憶部に記憶させた対応情報と、前記取得した電波の強度に基づく情報とに基づいて、前記電波の強度に基づく情報の取得元である端末装置の位置を特定し、前記測定時刻よりも後の運用時刻において前記電波の強度を示す情報を取得した場合、前記記憶部に前記対応情報として記憶させた複数の前記固有ベクトルのうち、固有値が大きい上位所定数の前記固有ベクトルを抽出し、前記抽出した所定数の固有ベクトルのそれぞれと、前記運用時刻における前記電波の強度を示す電波ベクトルとの内積を算出し、前記固有ベクトルと前記電波ベクトルとの内積の総和が最大となる前記座標ベクトルが示す位置を、前記運用時刻における前記端末装置の位置として特定することを特徴とする。 Another aspect of the present invention is to obtain information based on the intensity of the radio wave measured at a plurality of measurement times by the terminal device from a terminal device that receives the radio wave transmitted from the transmission source, and perform a spatial correlation method. based on, by solving the correlation matrix of the intensity of the radio wave according to a plurality of measurement times in which the acquired, the extracted eigenvectors indicating the stable component of the intensity of the radio wave in the measurement time, the eigenvectors extracted as the stabilizing component And the coordinate vector indicating the position of the terminal device known in advance are stored in the storage unit as the correspondence information in association with each other, and the correspondence information stored in the storage unit and the information based on the strength of the acquired radio wave. The position of the terminal device that is the acquisition source of the information based on the strength of the radio wave is specified based on, and the strength of the radio wave is indicated at the operation time after the measurement time. When the information is acquired, among the plurality of eigenvectors stored as the correspondence information in the storage unit, a predetermined upper number of eigenvectors having a large eigenvalue is extracted, and each of the extracted predetermined number of eigenvectors and the operation The inner product of the radio wave vector indicating the strength of the radio wave at the time is calculated, and the position indicated by the coordinate vector in which the sum of the inner products of the eigenvector and the radio wave vector is the maximum is the position of the terminal device at the operating time. It is characterized by specifying .

また、本発明の他の態様は、上記の位置特定支援装置のコンピュータに、発信源から送信される電波を受信する端末装置から、前記端末装置により複数の時刻において測定された電波の強度に基づく情報から、前記端末装置により抽出された前記測定時刻における前記電波の強度の安定成分を示す固有ベクトルを取得すること、前記安定成分として抽出した前記固有ベクトルと、予め知られた前記端末装置の位置を示す座標ベクトルとを対応付けて対応情報として記憶部に記憶させること、前記記憶部に記憶させた対応情報と、前記取得した固有ベクトルとに基づいて、前期固有ベクトルの取得元である端末装置の位置を特定すること、前記測定時刻よりも後の運用時刻において前記固有ベクトルを取得した場合、前記記憶部に前記対応情報として記憶させた複数の前記固有ベクトルのうち、固有値が大きい上位所定数の前記固有ベクトルを抽出すること、前記抽出した所定数の固有ベクトルのそれぞれと、前記運用時刻における前記固有ベクトルとの内積を算出すること、前記固有ベクトル同士の内積の総和が最大となる前記座標ベクトルが示す位置を、前記運用時刻における前記端末装置の位置として特定すること、を実行させるための位置特定プログラムである。

Another aspect of the present invention is based on the strength of radio waves measured by the terminal device at a plurality of times from a terminal device that receives radio waves transmitted from a transmission source to the computer of the position identification support device. Obtaining from the information an eigenvector indicating a stable component of the intensity of the radio wave at the measurement time extracted by the terminal device, indicating the eigenvector extracted as the stable component, and the position of the terminal device known in advance. be stored in the storage unit as the corresponding information in association with the coordinate vector, the corresponding information stored in the storage unit, the obtained based on the eigenvector, the location of the terminal device is a retrieve the previous year eigenvectors be specified when acquiring the eigenvectors in the operation time later than the measurement time, the corresponding information in the storage unit Extracting a predetermined upper number of eigenvectors having a large eigenvalue from among the plurality of stored eigenvectors, and calculating an inner product of each of the predetermined number of extracted eigenvectors and the eigenvector at the operating time. , A position specifying program for specifying , as the position of the terminal device at the operating time, the position indicated by the coordinate vector having the maximum sum of inner products of the eigenvectors .

以上説明したように、この発明によれば、処理負荷を軽減することができる位置特定支援装置、位置特定方法、および位置特定プログラムを提供できる。 As described above, according to the present invention, it is possible to provide a position specifying support device, a position specifying method, and a position specifying program that can reduce the processing load.

第１の実施形態における位置特定支援装置２００を含む位置特定システム１の一例を示す概略図である。It is a schematic diagram showing an example of position specific system 1 containing position specific supporting device 200 in a 1st embodiment. 第１の実施形態における端末装置１００の機能構成の一例を示す図である。It is a figure showing an example of functional composition of terminal unit 100 in a 1st embodiment. 第１の実施形態における位置特定支援装置２００の機能構成の一例を示す図である。It is a figure showing an example of functional composition of position specific assistance device 200 in a 1st embodiment. 第１の実施形態における位置特定支援装置２００の制御部２１０の事前処理の流れの一例を示すフローチャートである。It is a flow chart which shows an example of a flow of pre-processing of control part 210 of position specific assistance device 200 in a 1st embodiment. 事前処理の段階において、複数の発信源Ｓから端末装置１００が電波を受信する様子を上方から見た図である。It is the figure which looked at the terminal device 100 receiving a radio wave from a plurality of transmission sources S from the upper side at the stage of preprocessing. 第１の実施形態における端末装置１００の運用処理の流れの一例を示すフローチャートである。It is a flow chart which shows an example of a flow of operation processing of terminal unit 100 in a 1st embodiment. 第１の実施形態における位置特定支援装置２００の運用処理の流れの一例を示すフローチャートである。It is a flow chart which shows an example of a flow of operation processing of position specific assistance device 200 in a 1st embodiment. 第２の実施形態における端末装置１００の機能構成の一例を示す図である。It is a figure showing an example of functional composition of terminal unit 100 in a 2nd embodiment. 第２の実施形態における端末装置１００の運用処理の流れの一例を示すフローチャートである。It is a flowchart which shows an example of the flow of the operation process of the terminal device 100 in 2nd Embodiment.

以下、実施形態の位置特定支援装置、位置特定方法、および位置特定プログラムを、図面を参照して説明する。 Hereinafter, a position specifying support device, a position specifying method, and a position specifying program according to embodiments will be described with reference to the drawings.

（第１の実施形態）
図１は、第１の実施形態における位置特定支援装置２００を含む位置特定システム１の一例を示す概略図である。本実施形態における位置特定システム１において、位置特定支援装置２００は、端末装置１００が発信源Ｓ１からＳｎのうちいずれか１つ以上から受信した電波の信号強度ＲＳＳＩ（Received Signal Strength Indication）に基づいて、端末装置１００の位置を特定する。位置特定システム１は、発信源Ｓ１からＳｎ（ｎは任意の自然数）と、端末装置１００と、位置特定支援装置２００とを備える。 (First embodiment)
FIG. 1 is a schematic diagram showing an example of a position specifying system 1 including a position specifying support device 200 according to the first embodiment. In the position identification system 1 according to the present embodiment, the position identification support device 200 is based on the signal strength RSSI (Received Signal Strength Indication) of the radio wave that the terminal device 100 receives from any one or more of the transmission sources S1 to Sn. , The position of the terminal device 100 is specified. The position specifying system 1 includes sources S1 to Sn (n is an arbitrary natural number), a terminal device 100, and a position specifying support device 200.

発信源Ｓ１からＳｎは、例えば、Ｗｉ-ＦｉやＢｌｕｅｔｏｏｔｈ等の無線技術を用いて所定の電波を発信する。所定の電波は、例えばビーコンであり、複数の発信源Ｓｎをそれぞれ特定することが可能な識別情報を示す発信源ＩＤが含まれる。これによって、後述する端末装置１００は、受信した電波の発信元の発信源Ｓを特定することができる。 The transmission sources S1 to Sn emit predetermined radio waves using a wireless technology such as Wi-Fi or Bluetooth. The predetermined radio wave is, for example, a beacon, and includes a transmission source ID indicating identification information that can specify each of the plurality of transmission sources Sn. As a result, the terminal device 100 described later can specify the transmission source S that is the transmission source of the received radio wave.

発信源Ｓ１からＳｎは、例えば、デパートなどの建物内に設けられ、当該建物内に滞在する利用者が保有する端末装置１００に対して電波を発信する。なお、発信源Ｓ１からＳｎを特段区別しない場合、単に「発信源Ｓ」と記載する。 The transmission sources S1 to Sn are provided in, for example, a building such as a department store, and emit radio waves to the terminal device 100 owned by a user who stays in the building. In addition, when the sources S1 to Sn are not particularly distinguished, they are simply described as "source S".

端末装置１００は、例えば、ユーザが保有するスマートフォンや携帯電話等の電子デバイスである。図２は、第１の実施形態における端末装置１００の機能構成の一例を示す図である。第１の実施形態における端末装置１００は、アンテナ１０２と、変調／復調部１０４と、ＡＤ／ＤＡ変換部１０６と、入力部１０８と、出力部１１０と、制御部１３０と、記憶部１５０とを備える。 The terminal device 100 is, for example, an electronic device such as a smartphone or a mobile phone owned by the user. FIG. 2 is a diagram illustrating an example of a functional configuration of the terminal device 100 according to the first embodiment. The terminal device 100 according to the first embodiment includes an antenna 102, a modulation / demodulation unit 104, an AD / DA conversion unit 106, an input unit 108, an output unit 110, a control unit 130, and a storage unit 150. Prepare

アンテナ１０２は、例えば、発信源Ｓにより発信される電波を受信し、受信した電波に基づく信号を変調／復調部１０４に入力する。また、アンテナ１０２は、変調／復調部１０４により変調された信号に基づいた電波を外部装置に送信する。外部装置は、例えば、位置特定支援装置２００、無線基地局、ルータやプロキシサーバ等の通信中継装置等である。 The antenna 102 receives, for example, a radio wave transmitted from the transmission source S and inputs a signal based on the received radio wave to the modulation / demodulation unit 104. Further, the antenna 102 transmits a radio wave based on the signal modulated by the modulation / demodulation unit 104 to an external device. The external device is, for example, the position identification support device 200, a wireless base station, a communication relay device such as a router or a proxy server, or the like.

変調／復調部１０４は、アンテナ１０２により受信された信号をダウンコンバートした後に復調し、復調した信号（アナログ信号）をＡＤ／ＤＡ変換部１０６に入力する。また、変調／復調部１０４は、自装置内で生成されＡＤ／ＤＡ変換部１０６によりアナログ信号に変換された信号を変調し、アップコンバートしてアンテナ１０２に出力する。 The modulation / demodulation unit 104 down-converts the signal received by the antenna 102, demodulates it, and inputs the demodulated signal (analog signal) to the AD / DA conversion unit 106. The modulation / demodulation unit 104 also modulates a signal generated in the device itself and converted into an analog signal by the AD / DA conversion unit 106, up-converts it, and outputs it to the antenna 102.

ＡＤ／ＤＡ変換部１０６は、変調／復調部１０４から入力されたアナログ信号をディジタル信号に変換し、変換したディジタル信号を制御部１３０に入力する。また、ＡＤ／ＤＡ変換部１０６は、制御部１３０から出力されるディジタル信号をアナログ信号に変換し、変換したアナログ信号を変調／復調部１０４に入力する。 The AD / DA conversion unit 106 converts the analog signal input from the modulation / demodulation unit 104 into a digital signal, and inputs the converted digital signal to the control unit 130. Further, the AD / DA conversion unit 106 converts the digital signal output from the control unit 130 into an analog signal, and inputs the converted analog signal to the modulation / demodulation unit 104.

入力部１０８は、ユーザの操作を受け付けて、操作に応じた信号を生成する。入力部１０８は、例えば、タッチパネル、キーボード、マウス等である。 The input unit 108 receives a user operation and generates a signal according to the operation. The input unit 108 is, for example, a touch panel, a keyboard, a mouse, or the like.

出力部１１０は、入力された信号に基づいた画像を表示する。出力部１１０は、液晶ディスプレイ、有機ＥＬ（Electroluminescence）ディスプレイ等である。なお、出力部１１０は、入力部１０８の機能を有するタッチパネルであってもよい。 The output unit 110 displays an image based on the input signal. The output unit 110 is a liquid crystal display, an organic EL (Electroluminescence) display, or the like. The output unit 110 may be a touch panel having the function of the input unit 108.

制御部１３０は、ＲＳＳＩデータ生成部１３２と、出力制御部１３４とを備える。上述した制御部１３０の機能部のうち一部または全部は、ＣＰＵ（Central Processing Unit）等のプロセッサが記憶部１５０に記憶されたプログラムを実行することにより機能するソフトウェア機能部である。プログラムは、例えば、ネットワークＮＷを介してアプリサーバからダウンロードされる。また、制御部１３０の機能部のうち一部または全部は、ＬＳＩ（Large Scale Integration）、ＡＳＩＣ（Application Specific Integrated Circuit）等のハードウェア機能部であってもよい。 The control unit 130 includes an RSSI data generation unit 132 and an output control unit 134. Some or all of the functional units of the control unit 130 described above are software functional units that function by a processor such as a CPU (Central Processing Unit) executing a program stored in the storage unit 150. The program is downloaded from the application server via the network NW, for example. Further, some or all of the functional units of the control unit 130 may be hardware functional units such as LSI (Large Scale Integration) and ASIC (Application Specific Integrated Circuit).

記憶部１５０は、例えば、ＲＯＭ（Read Only Memory）、フラッシュメモリ、ＳＤカード等の不揮発性の記憶媒体と、ＲＡＭ（Random Access Memory）、レジスタ等の揮発性の記憶媒体とを有する。記憶部１５０に記憶される情報は、プロセッサが実行するプログラムの他、発信源ＩＤ、後述する図６に示すフローチャートの処理を制御部１３０（プロセッサ）に実行させるためのアプリケーションプログラム等の情報を含む。 The storage unit 150 has, for example, a non-volatile storage medium such as a ROM (Read Only Memory), a flash memory, and an SD card, and a volatile storage medium such as a RAM (Random Access Memory) and a register. The information stored in the storage unit 150 includes a program executed by the processor, a source ID, and an application program for causing the control unit 130 (processor) to execute the process of the flowchart illustrated in FIG. .

上述した制御部１３０におけるＲＳＳＩデータ生成部１３２は、ＡＤ／ＤＡ変換部１０６により変換された信号、すなわち発信源Ｓから受信した電波に基づく信号を基に、発信源ＩＤと対応付けた電波の信号強度を表すデータを生成する。以下、発信源ＩＤと対応付けた電波の信号強度を表すデータを、「ＲＳＳＩデータ」と称する。例えば、ＲＳＳＩデータ生成部１３２は、発信源ＳごとにＲＳＳＩデータを生成する。 The RSSI data generation unit 132 in the control unit 130 described above uses the signal converted by the AD / DA conversion unit 106, that is, the signal of the radio wave associated with the transmission source ID based on the signal based on the radio wave received from the transmission source S. Generate data representing intensity. Hereinafter, the data indicating the signal strength of the radio wave associated with the transmission source ID will be referred to as “RSSI data”. For example, the RSSI data generation unit 132 generates RSSI data for each transmission source S.

より具体的には、ＲＳＳＩデータ生成部１３２は、発信源Ｓｋから発信された電波に基づく信号に対して、発信源Ｓｋの発信源ＩＤを対応付けたＲＳＳＩデータ（発信源毎データＲＳＳＩ_Ｓｋ）を生成する（ｋ＝１〜ｎ）。このようにして、ＲＳＳＩデータ生成部１３２は、発信源Ｓの数ｎに応じて、ＲＳＳＩデータを生成する。 More specifically, the RSSI data generation unit 132 generates RSSI data (source-specific data RSSI _Sk ) in which the source ID of the source Sk is associated with the signal based on the radio wave transmitted from the source Sk. Generate (k = 1 to n). In this way, the RSSI data generation unit 132 generates RSSI data according to the number n of the transmission sources S.

出力制御部１３４は、ＲＳＳＩデータ生成部１３２により生成された発信源毎データＲＳＳＩ_Ｓｋを、ネットワークＮＷを介して位置特定支援装置２００に送信するように、アンテナ１０２、変調／復調部１０４、およびＡＤ／ＤＡ変換部１０６を制御する。ネットワークＮＷは、ＬＡＮ（Local Area Network）やＷＡＮ（Wide Area Network）、シリアル通信線、無線基地局、通信中継装置等を含む。前述した発信源Ｓ１からＳｎも同様にネットワークＮＷに接続され得るため、端末装置１００がネットワークにアクセスする際の窓口装置は発信源Ｓ１からＳｎのいずれかであってよい。また、出力制御部１３４は、アンテナ１０２を用いた無線通信の代わりに、ＵＳＢ、ＰＣＩ等の図示しない通信インターフェースを用いた有線接続によって、発信源毎データＲＳＳＩ_Ｓｋを位置特定支援装置２００に送信するように制御してもよい。 The output control unit 134 causes the antenna 102, the modulation / demodulation unit 104, and the AD so that the transmission source data RSSI _Sk generated by the RSSI data generation unit 132 is transmitted to the position identification support device 200 via the network NW. The / DA converter 106 is controlled. The network NW includes a LAN (Local Area Network), a WAN (Wide Area Network), a serial communication line, a wireless base station, a communication relay device, and the like. Since the above-mentioned transmission sources S1 to Sn can be similarly connected to the network NW, the window device when the terminal device 100 accesses the network may be any of the transmission sources S1 to Sn. Further, the output control unit 134 transmits the data RSSI _Sk for each transmission source to the position identification support apparatus 200 by wire connection using a communication interface (not shown) such as USB or PCI, instead of wireless communication using the antenna 102. May be controlled as follows.

本実施形態における位置特定支援装置２００は、運用開始までに行う事前処理と、運用後に行う運用処理との２種類の処理パターンを有する。位置特定支援装置２００は、事前処理として、上述した端末装置１００から発信源毎データＲＳＳＩ_Ｓｋを予め取得しておき、取得した発信源毎データＲＳＳＩ_Ｓｋに基づいて、運用時に参照するためのデータベースを生成する。また、位置特定支援装置２００は、運用処理として、端末装置１００からネットワークＮＷを介して取得した発信源毎データＲＳＳＩ_Ｓｋに基づき、予め生成したデータベースを参照して、その時点における端末装置１００の位置を特定する。なお、事前処理時の端末装置１００と運用処理時の端末装置１００は異なってもよいし、同じであってもよい。 The position identification support device 200 according to the present embodiment has two types of processing patterns, that is, pre-processing performed before the start of operation and operation processing performed after the operation. As the pre-processing, the position identification support device 200 acquires in advance the source-specific data RSSI _Sk from the above-described terminal device 100, and based on the acquired source-specific data RSSI _Sk , a database to be referred to during operation. To generate. Further, the position identification support device 200 refers to a database generated in advance based on the data RSSI _Sk for each transmission source acquired from the terminal device 100 via the network NW as an operation process, and the position of the terminal device 100 at that time. Specify. The terminal device 100 at the time of preprocessing and the terminal device 100 at the time of operation processing may be different or the same.

以下、位置特定支援装置２００の構成図、およびフローチャートを参照して、位置特定支援装置２００における２種類の処理パターンを具体的に説明する。図３は、第１の実施形態における位置特定支援装置２００の機能構成の一例を示す図である。 Hereinafter, two types of processing patterns in the position identification support device 200 will be specifically described with reference to the configuration diagram of the position identification support device 200 and a flowchart. FIG. 3 is a diagram illustrating an example of a functional configuration of the position identification support device 200 according to the first embodiment.

本実施形態における位置特定支援装置２００は、通信インターフェース２０２と、入力部２０４と、表示部２０６と、制御部２１０と、記憶部２３０とを備える。通信インターフェース２０２は、ネットワークＮＷに接続するための通信インターフェースであり、例えば、ネットワークカード等を含む。 The position identification support device 200 in this embodiment includes a communication interface 202, an input unit 204, a display unit 206, a control unit 210, and a storage unit 230. The communication interface 202 is a communication interface for connecting to the network NW, and includes, for example, a network card.

入力部２０４は、ユーザの操作を受け付けて、操作に応じた情報を生成する。入力部２０４は、例えば、タッチパネル、キーボード、マウス等である。入力部２０４には、例えば、事前処理の段階において、発信源Ｓの発信源ＩＤや、発信源Ｓから事前に電波を受信した際の端末装置１００の位置等が、ユーザにより入力される。 The input unit 204 receives a user operation and generates information according to the operation. The input unit 204 is, for example, a touch panel, a keyboard, a mouse, or the like. In the input unit 204, for example, the source ID of the source S, the position of the terminal device 100 when a radio wave is previously received from the source S, and the like are input by the user in the pre-processing stage.

表示部２０６は、入力された信号に基づいた画像を表示する。表示部２０６は、液晶ディスプレイ、有機ＥＬディスプレイ等の表示装置である。なお、表示部２０６は、入力部２０４の機能を有するタッチパネルであってもよい。 The display unit 206 displays an image based on the input signal. The display unit 206 is a display device such as a liquid crystal display or an organic EL display. The display unit 206 may be a touch panel having the function of the input unit 204.

位置特定支援装置２００の制御部２１０は、取得部２１２と、相関行列演算部２１４と、固有ベクトル導出部２１６と、記憶制御部２１８と、類似度演算部２２０と、位置特定部２２２とを備える。上述した制御部２１０の機能部のうち一部または全部は、ＣＰＵ等のプロセッサが記憶部２３０に記憶されたプログラムを実行することにより機能するソフトウェア機能部である。また、制御部２１０の機能部のうち一部または全部は、ＬＳＩ、ＡＳＩＣ等のハードウェア機能部であってもよい。なお、上述した相関行列演算部２１４および固有ベクトル導出部２１６は、「抽出部」の一例である。また、類似度演算部２２０および位置特定部２２２は、「特定部」の一例である。 The control unit 210 of the position identification support device 200 includes an acquisition unit 212, a correlation matrix calculation unit 214, an eigenvector derivation unit 216, a storage control unit 218, a similarity calculation unit 220, and a position identification unit 222. Some or all of the functional units of the control unit 210 described above are software functional units that function when a processor such as a CPU executes a program stored in the storage unit 230. Further, some or all of the functional units of the control unit 210 may be hardware functional units such as LSI and ASIC. The correlation matrix calculation unit 214 and the eigenvector derivation unit 216 described above are examples of an “extraction unit”. Further, the similarity calculation unit 220 and the position specifying unit 222 are examples of a “specifying unit”.

位置特定支援装置２００の記憶部２３０は、例えば、ＲＯＭ、フラッシュメモリ、ＨＤＤ等の不揮発性の記憶媒体と、ＲＡＭ、レジスタ等の揮発性の記憶媒体とを有する。記憶部２３０に記憶される情報は、位置特定支援装置２００のプロセッサが実行するプログラムの他、後述するフィンガープリント、数式、ＲＳＳＩデータ（発信源毎データＲＳＳＩ_Ｓｋ）等の情報を含む。 The storage unit 230 of the position identification support device 200 has, for example, a non-volatile storage medium such as a ROM, a flash memory, and an HDD, and a volatile storage medium such as a RAM and a register. The information stored in the storage unit 230 includes a program executed by the processor of the position identification support device 200, as well as information such as fingerprints, mathematical expressions, and RSSI data (source-specific data RSSI _Sk ) described later.

（事前処理）
以下、制御部２１０が実行する事前処理をフローチャートに則して説明する。図４は、第１の実施形態における位置特定支援装置２００の制御部２１０の事前処理の流れの一例を示すフローチャートである。 (Pre-processing)
Hereinafter, the pre-processing executed by the control unit 210 will be described with reference to the flowchart. FIG. 4 is a flowchart showing an example of a pre-processing flow of the control unit 210 of the position identification support device 200 according to the first embodiment.

まず、取得部２１２は、通信インターフェース２０２に接続されるネットワークＮＷを介して端末装置１００から発信源毎データＲＳＳＩ_Ｓｋを取得し（ステップＳ１００）、相関行列演算部２１４に出力する。次に、取得部２１２は、ユーザにより入力された、発信源Ｓから電波を受信した際の端末装置１００の位置と、発信源Ｓの発信源ＩＤとを、入力部２０４から取得する（ステップＳ１０２）。 First, the acquisition unit 212 acquires the data RSSI _Sk for each transmission source from the terminal device 100 via the network NW connected to the communication interface 202 (step S100) and outputs it to the correlation matrix calculation unit 214. Next, the acquisition unit 212 acquires, from the input unit 204, the position of the terminal device 100 when the radio wave is received from the transmission source S and the transmission source ID of the transmission source S, which are input by the user (step S102). ).

次に、相関行列演算部２１４は、取得部２１２により取得された発信源毎データＲＳＳＩ_Ｓｋに基づいて、ＲＳＳＩデータ数ｎに対応するｎ次元のベクトルを生成する（ステップＳ１０４）。以下、ＲＳＳＩデータ数ｎに対応するｎ次元のベクトルを、「ＲＳＳＩベクトル→ｘ」と称し、ＲＳＳＩベクトル→ｘの一例を数式（１）に示す。なお、矢印→は、ベクトルを表すものとする。 Next, the correlation matrix calculation unit 214 generates an n-dimensional vector corresponding to the number of RSSI data n based on the source-specific data RSSI _Sk acquired by the acquisition unit 212 (step S104). Hereinafter, an n-dimensional vector corresponding to the number n of RSSI data is referred to as “RSSI vector → x”, and an example of the RSSI vector → x is shown in Expression (1). The arrow → indicates a vector.

数式（１）に表すように、ＲＳＳＩベクトル→ｘは、端末装置１００が各発信源Ｓｎから受信した電波強度をそれぞれ示す発信源毎データＲＳＳＩ_Ｓ１から発信源毎データＲＳＳＩ_Ｓｎを並べたものである。ここで、端末装置１００が存在し得る空間の座標を、数式（２）に示すように→ξにより表す。数式（２）に示すように空間座標→ξは、例えば、Ｘ軸上の座標（Ｘ成分）と、Ｘ軸に直行するＹ軸上の座標（Ｙ成分）とによって表される。 As represented by the mathematical expression (1), the RSSI vector → x is an array of source-source data RSSI _S1 to source-source data RSSI _Sn that indicate the radio field intensity received by the terminal device 100 from each source Sn. . Here, the coordinates of the space in which the terminal device 100 can exist are represented by → ξ as shown in Expression (2). As shown in Expression (2), the spatial coordinate → ξ is represented by, for example, a coordinate on the X axis (X component) and a coordinate on the Y axis (Y component) orthogonal to the X axis.

図５は、事前処理の段階において、複数の発信源Ｓから端末装置１００が電波を受信する様子を上方から見た図である。図５の例の場合、端末装置１００は、座標→ξにより示される位置で、発信源Ｓ１からＳｎのそれぞれから電波を所定時間（例えば１０秒間程度）受信する。端末装置１００は、例えば、所定時間として１０秒間の間に１回でも電波を受信することができた発信源Ｓを、着目対象に含めるようにしてｎを決定し、その上で着目対象に含めるように決定した発信源Ｓについて、電波を受信できない時刻に対応する発信源毎データＲＳＳＩ_Ｓｋの値を“０”にする。 FIG. 5 is a view from above showing how the terminal device 100 receives radio waves from a plurality of transmission sources S in the pre-processing stage. In the case of the example in FIG. 5, the terminal device 100 receives a radio wave from each of the transmission sources S1 to Sn for a predetermined time (for example, about 10 seconds) at a position indicated by coordinates → ξ. The terminal device 100 determines, for example, the source S that was able to receive the radio wave even once in 10 seconds as the predetermined time so as to be included in the attention target, and then determines n to be included in the attention target. With respect to the transmission source S thus determined, the value of the transmission source data RSSI _Sk corresponding to the time when radio waves cannot be received is set to "0".

本実施形態における端末装置１００は、ある時刻において、複数の発信源Ｓの全てから電波を受信しない場合と、いずれか１つ以上（全てを含む）の発信源Ｓから電波を受信する場合との２つのケースが存在する。例えば、端末装置１００が複数の発信源Ｓの全てから電波を受信しない場合、相関行列演算部２１４は、上述した数式（１）に示すＲＳＳＩベクトル→ｘ内の成分（要素）を、いずれも“０”としたベクトルを生成する。また、端末装置１００が複数の発信源Ｓｎのうち、いずれか１つ以上（全てを含む）から電波を受信する場合、相関行列演算部２１４は、上述した数式（１）に示すＲＳＳＩベクトル→ｘ内の成分（要素）のうち、電波を受信できた発信源Ｓに対応した成分（要素）を、受信電波の強度値に対応した数値に置き換えたベクトルを生成する。 The terminal device 100 according to the present embodiment may receive radio waves from all of the plurality of transmission sources S at a certain time, or may receive radio waves from any one or more (including all) transmission sources S. There are two cases. For example, when the terminal device 100 does not receive radio waves from all of the plurality of transmission sources S, the correlation matrix calculation unit 214 makes all of the components (elements) in the RSSI vector → x shown in the above mathematical expression (1) “. A vector with 0 "is generated. Further, when the terminal device 100 receives radio waves from any one or more (including all) of the plurality of transmission sources Sn, the correlation matrix calculation unit 214 causes the RSSI vector → x shown in the mathematical expression (1) described above. Among these components (elements), a vector is generated by replacing the component (element) corresponding to the transmission source S that was able to receive the radio wave with a numerical value corresponding to the intensity value of the received radio wave.

また、端末装置１００が、所定時間において同一の座標→ξの位置で発信源Ｓからの電波を測定し続けるため、取得部２１２は、座標→ξにおいて受信された複数の時系列の発信源毎データＲＳＳＩ_Ｓｋを端末装置１００から取得し、取得した複数の時系列の発信源毎データＲＳＳＩ_Ｓｋを相関行列演算部２１４に出力する。相関行列演算部２１４は、取得部２１２により出力された複数の時系列の発信源毎データＲＳＳＩ_Ｓｋを用いて、ＲＳＳＩベクトル→ｘを所定数ｎ_→ξ生成する。すなわち、相関行列演算部２１４は、サンプル数ｎ_→ξ個のＲＳＳＩベクトル→ｘを生成する。相関行列演算部２１４により生成された複数のＲＳＳＩベクトル→ｘのうち、例えば、ｉ回目に生成されたＲＳＳＩベクトル→ｘ_→ξｉを数式（３）に示し、ｊ回目に生成されたＲＳＳＩベクトル→ｘ_→ξｊを数式（４）に示す。 Further, since the terminal device 100 continues to measure the radio wave from the transmission source S at the same coordinate → ξ position for a predetermined time, the acquisition unit 212 determines that each of the plurality of time-series transmission sources received at the coordinate → ξ. The data RSSI _Sk is acquired from the terminal device 100, and the acquired plurality of time-series data per source RSSI _Sk is output to the correlation matrix calculation unit 214. The correlation matrix calculation unit 214 uses the plurality of time-series data per source RSSI _Sk output by the acquisition unit 212 to generate a predetermined number n _{→ ξ of} RSSI vector → x. That is, the correlation matrix calculation unit 214 generates the number of samples n _{→ ξ} RSSI vectors → x. Of the plurality of RSSI vectors → x generated by the correlation matrix calculation unit 214, for example, the RSSI vector → x _{→ ξi} generated at the i-th time is shown in Equation (3), and the RSSI vector → x generated at the j-th time. _{→ ξj} is shown in Equation (4).

次に、相関行列演算部２１４は、上述したように所定数ｎ_→ξ個生成したＲＳＳＩベクトル→ｘに対して、以下の数式（５）に示す自己相関行列Γ_→ξを計算する。 Next, the correlation matrix calculation unit 214 calculates the autocorrelation matrix Γ _{→ ξ} shown in the following mathematical expression (5) for the RSSI vector → x generated by the predetermined number n _{→ ξ} as described above.

数式（５）中に示すｎ_→ξは上述した位置→ξにおけるＲＳＳＩベクトル→ｘのサンプル数を表す。相関行列演算部２１４は、数式（５）に示すように、ｎ_→ξ個のＲＳＳＩベクトル→ｘから、例えば、１回目に生成したＲＳＳＩベクトル→ｘ_→ξ１と、このＲＳＳＩベクトル→ｘ_→ξ１の転置行列とを乗算し、２回目以降も同様に生成した各ＲＳＳＩベクトル→ｘに対してもその転置行列を乗算する。相関行列演算部２１４は、このような自己相関処理をｎ_→ξ個の全てのＲＳＳＩベクトル→ｘに対して行い、得られた乗算値を全て加算し、ｎ_→ξ値で除算して平均を取る。このようにして得られた行列は、Γ_→ξ＝Γ_→ξ ^Ｔを満たす対称行列となる。 N _{→ ξ} shown in the mathematical expression (5) represents the number of samples of the RSSI vector → x at the position → ξ described above. The correlation matrix calculation unit 214 calculates, as shown in Expression (5), from the n _{→ ξ} RSSI vectors → x, for example, the RSSI vector → x _{→ ξ1} generated at the first time and the RSSI vector → x _{→ ξ1} The transposed matrix is multiplied, and each RSSI vector → x similarly generated after the second time is also multiplied by the transposed matrix. The correlation matrix calculation unit 214 performs such autocorrelation processing on all n _{→ ξ} RSSI vectors → x, adds all the obtained multiplication values, divides by n _{→ ξ} values, and calculates the average. take. The matrix thus obtained is a symmetric matrix that satisfies Γ _{→ ξ} = Γ _{→ ξ} ^T.

固有ベクトル導出部２１６は、相関行列演算部２１４により導出された自己相関行列に基づいて、固有値問題を解き、固有ベクトル→ψ_→ξｋを導出する（ステップＳ１０６）。 The eigenvector deriving unit 216 solves the eigenvalue problem based on the autocorrelation matrix derived by the correlation matrix computing unit 214 and derives the eigenvector → ψ _{→ ξk} (step S106).

例えば、端末装置１００により受信される電波には、電波強度が常に安定しているものと、不安定なものとが存在している。この場合、大きな成分は、座標→ξに位置する端末装置１００が、所定時間内において安定的に電波を受信することができていた発信源Ｓの集合を表している。すなわち、固有ベクトル導出部２１６は、相関行列演算部２１４により導出された自己相関行列から、所定時間内において端末装置１００により受信された電波の信号強度の変化が小さく安定した成分を抽出する。固有ベクトル導出部２１６は、数式（６）に示す固有方程式を解くことで、座標→ξに対応した固有ベクトル→ψ_→ξｋを算出する。 For example, there are radio waves received by the terminal device 100 that are always stable and unstable. In this case, the large component represents the set of transmission sources S that the terminal device 100 located at the coordinates → ξ could stably receive the radio wave within the predetermined time. That is, the eigenvector deriving unit 216 extracts, from the autocorrelation matrix derived by the correlation matrix computing unit 214, a stable component with a small change in the signal strength of the radio wave received by the terminal device 100 within a predetermined time. The eigenvector deriving unit 216 calculates the eigenvector → ψ _{→ ξk} corresponding to the coordinate → ξ by solving the _{eigenequation} shown in Expression (6).

次に、記憶制御部２１８は、固有ベクトル導出部２１６により算出された固有ベクトル→ψ_→ξｋを、座標→ξに対応させて記憶部２３０に記憶させる（ステップＳ１０８）。記憶制御部２１８は、記憶部２３０に座標→ξと対応させて記憶させた固有ベクトル→ψ_→ξｋの数が所定数に達したか否かを判定し（ステップＳ１１０）、所定数未満の場合に、上述した事前処理を繰り返し行うように、取得部２１２、相関行列演算部２１４、および固有ベクトル導出部２１６を制御する。制御部２１０は、記憶部２３０に座標→ξと対応させて記憶させた固有ベクトル→ψ_→ξｋの数が所定数に達した場合に、図４に示すフローチャートの処理を終了する。 Next, the storage control unit 218 stores the eigenvector → ψ _{→ ξk} calculated by the eigenvector deriving unit 216 in the storage unit 230 in _association with the coordinate → ξ (step S108). The storage control unit 218 determines whether or not the number of eigenvectors → ψ _{→ ξk} stored in the storage unit 230 in association with the coordinate → ξ reaches a predetermined number (step S110). The acquisition unit 212, the correlation matrix calculation unit 214, and the eigenvector derivation unit 216 are controlled so as to repeatedly perform the above-described preprocessing. When the number of eigenvectors → ψ _{→ ξk} stored in the storage unit 230 in association with the coordinate → ξ reaches a predetermined number, the control unit 210 ends the process of the flowchart shown in FIG. 4.

このように、位置特定支援装置２００の制御部２１０は事前処理を行い、端末装置１００が発信源Ｓからの電波を測定した際の位置に応じて、所定数の固有ベクトル→ψ_→ξｋを導出し、導出した固有ベクトル→ψ_→ξｋを、ユーザによって入力部２０４に入力された端末装置１００の位置を示す座標→ξに対応付けて記憶部２３０に記憶させることにより、運用時において参照するデータベースを生成することができる。以下、事前処理を行って構築したデータベースを、「フィンガープリント」と称する。なお、フィンガープリントは、「対応情報」の一例である。 In this way, the control unit 210 of the position identification support device 200 performs the _{preprocessing} and derives a predetermined number of eigenvectors → ψ _{→ ξk} according to the position when the terminal device 100 measures the radio wave from the transmission source S. , The derived eigenvector → ψ _{→ ξk} is stored in the storage unit 230 in association with the coordinate → ξ indicating the position of the terminal device 100 input by the user in the input unit 204, thereby generating a database to be referred to during operation. can do. Hereinafter, the database constructed by performing the pre-processing is referred to as "fingerprint". The fingerprint is an example of “correspondence information”.

（運用処理）
以下、運用処理における、端末装置１００側の処理と、位置特定支援装置２００側の処理とを、それぞれ別々のフローチャートに則して説明する。まず、運用処理における端末装置１００側の処理について説明する。図６は、第１の実施形態における端末装置１００の運用処理の流れの一例を示すフローチャートである。 (Operation processing)
Hereinafter, the processing on the side of the terminal device 100 and the processing on the side of the position identification supporting device 200 in the operation processing will be described respectively according to different flowcharts. First, the processing on the terminal device 100 side in the operation processing will be described. FIG. 6 is a flowchart showing an example of the flow of operation processing of the terminal device 100 according to the first embodiment.

端末装置１００の制御部１３０は、記憶部１５０に予め記憶されたアプリケーションプログラムを実行し（ステップＳ２００）、例えば、以下の処理をバックグラウンドで実施させる。次に、制御部１３０は、アンテナ１０２を介して発信源Ｓから電波を受信したか否かを判定する（ステップＳ２０２）。 The control unit 130 of the terminal device 100 executes the application program stored in advance in the storage unit 150 (step S200) and, for example, causes the following processing to be performed in the background. Next, the control unit 130 determines whether or not a radio wave has been received from the transmission source S via the antenna 102 (step S202).

ＲＳＳＩデータ生成部１３２は、アンテナ１０２を介して発信源Ｓから電波を受信した場合、発信源Ｓから受信した電波に基づく信号を基に、ＲＳＳＩデータを生成する（ステップＳ２０４）。制御部１３０は、アンテナ１０２を介して発信源Ｓから電波を受信していない場合、ステップＳ２０４の処理を行わずに、ステップＳ２０６の処理を行う。 When the radio wave is received from the transmission source S via the antenna 102, the RSSI data generation unit 132 generates the RSSI data based on the signal based on the radio wave received from the transmission source S (step S204). When the radio wave is not received from the transmission source S via the antenna 102, the control unit 130 performs the process of step S206 without performing the process of step S204.

次に、制御部１３０は、発信源Ｓからの電波を測定し始めてから、所定時間が経過したか否かを判定する（ステップＳ２０６）。制御部１３０は、所定時間が経過しない場合、ステップＳ２０２の処理に戻る。出力制御部１３４は、所定時間が経過した場合、アンテナ１０２、変調／復調部１０４、およびＡＤ／ＤＡ変換部１０６を制御して、当該所定時間において蓄積した複数のＲＳＳＩデータを、位置特定支援装置２００に出力する（ステップＳ２０８）。なお、出力制御部１３４は、例えば、所定時間において蓄積した複数のＲＳＳＩデータの平均値を、位置特定支援装置２００に出力するようにしてもよい。これによって、制御部１３０は、本フローチャートの処理を終了する。 Next, the control unit 130 determines whether or not a predetermined time has elapsed since the measurement of the radio wave from the transmission source S was started (step S206). When the predetermined time has not elapsed, the control unit 130 returns to the process of step S202. When a predetermined time has elapsed, the output control unit 134 controls the antenna 102, the modulation / demodulation unit 104, and the AD / DA conversion unit 106 so that the plurality of RSSI data accumulated during the predetermined time is stored in the position identification support device. It outputs to 200 (step S208). The output control unit 134 may output the average value of a plurality of RSSI data accumulated in a predetermined time to the position identification support device 200, for example. As a result, the control unit 130 ends the processing of this flowchart.

次に、運用過程における位置特定支援装置２００側の処理について説明する。図７は、第１の実施形態における位置特定支援装置２００の運用処理の流れの一例を示すフローチャートである。 Next, a process on the side of the position identification supporting device 200 in the operation process will be described. FIG. 7 is a flowchart showing an example of the flow of operation processing of the position identification support device 200 in the first embodiment.

取得部２１２は、端末装置１００により出力された複数のＲＳＳＩデータを、ネットワークＮＷを介して取得する（ステップＳ３００）。次に、相関行列演算部２１４は、取得部２１２により取得された複数のＲＳＳＩデータに基づいて、ＲＳＳＩベクトル→ｘを生成する（ステップＳ３０２）。 The acquisition unit 212 acquires a plurality of RSSI data output by the terminal device 100 via the network NW (step S300). Next, the correlation matrix calculation unit 214 generates the RSSI vector → x based on the plurality of RSSI data acquired by the acquisition unit 212 (step S302).

次に、類似度演算部２２０は、記憶部２３０に記憶されているフィンガープリントと、相関行列演算部２１４により生成されたＲＳＳＩベクトル→ｘとに基づいて、位置→ξに対応する類似度Ｃ_ξを算出する。数式（７）は、類似度Ｃの算出式の一例である。例えば、類似度演算部２２０は、上述した数式（７）に示すように、正規化したＲＳＳＩベクトル→ｘを、フィンガープリントとして記憶されている固有ベクトル→ψ_→ξｋに対して射影して類似度Ｃ_ξを算出する。なお、類似度演算部２２０は、例えば、ＲＳＳＩベクトル→ｘと固有ベクトル→ψ_→ξｋとの二乗誤差を用いて類似度Ｃ_ξを算出してもよい。 Next, the similarity calculation unit 220, based on the fingerprint stored in the storage unit 230 and the RSSI vector → x generated by the correlation matrix calculation unit 214, the similarity C _ξ corresponding to the position → _ξ To calculate. Formula (7) is an example of a formula for calculating the similarity C. For example, the similarity calculation unit 220 _projects the normalized RSSI vector → x onto the eigenvector → ψ _{→ ξk} stored as the fingerprint to _calculate the similarity C, as shown in the above-described mathematical expression (7). Calculate _ξ . Note that the similarity calculation unit 220 may calculate the similarity C _ξ , for example, by using the squared error of the RSSI vector → x and the eigenvector → ψ _{→ ξk} .

数式（７）に示すように、類似度演算部２２０は、相関行列演算部２１４により生成されたＲＳＳＩベクトル→ｘを正規化し、正規化したＲＳＳＩベクトル→ｘと、フィンガープリントとして記憶されている固有ベクトル→ψ_→ξｋとの内積を算出する。類似度演算部２２０は、複数の固有ベクトル→ψ_→ξｋから、他の固有ベクトル→ψ_→ξｋとの相関が高いｄ個の固有ベクトル→ψ_→ξｋを抽出する。例えば、類似度演算部２２０は、相関行列演算部２１４により導出された自己相関行列に対する固有値λの大きいものから順に所定数ｄ個の固有ベクトル→ψ_→ξｋを抽出し、抽出したｄ個の固有ベクトル→ψ_→ξｋのそれぞれに対して、上述したＲＳＳＩベクトル→ｘとの内積を算出し、算出した全ての内積の総和を位置→ξに対応する類似度Ｃ_ξとして導出する。 As shown in Expression (7), the similarity calculation unit 220 normalizes the RSSI vector → x generated by the correlation matrix calculation unit 214 and normalizes the RSSI vector → x, and the eigenvector stored as the fingerprint. Calculate the inner product of → ψ _{→ ξk} . Similarity calculating unit 220, a plurality of eigenvectors _→ ψ → ξk, extracts the other eigenvectors _→ ψ → ξk correlation is high d eigenvectors of _→ ψ → ξk. For example, the similarity calculation unit 220 extracts a predetermined number d of eigenvectors → ψ _{→ ξk} in order from the largest eigenvalue λ of the autocorrelation matrix derived by the correlation matrix calculation unit 214, and extracts the extracted d eigenvectors → For each of ψ _{→ ξk} , the inner product of the above-mentioned RSSI vector → x is calculated, and the sum of all the calculated inner products is derived as the similarity C _ξ corresponding to the position → ξ.

次に、位置特定部２２２は、類似度演算部２２０により算出された類似度Ｃ_ξに基づいて、取得部２１２の取得元である端末装置１００の位置を特定する（ステップＳ３０６）。具体的には、位置特定部２２２は、観測データであるＲＳＳＩベクトル→ｘとの類似度Ｃ_ξが最大となる座標→ξが示す位置を、取得部２１２の取得元である端末装置１００の位置として特定する。位置特定部２２２は、通信インターフェース２０２を介して、特定した端末装置１００の位置を、例えば、地図等の画像情報に加工して取得部２１２の取得元である端末装置１００に出力する（ステップＳ３０８）。これによって、制御部１３０は、本フローチャートの処理を終了する。 Next, the position specifying unit 222 specifies the position of the terminal device 100 that is the acquisition source of the acquisition unit 212 based on the similarity C _ξ calculated by the similarity calculation unit 220 (step S306). Specifically, the position specifying unit 222 determines the position of the terminal device 100, which is the acquisition source of the acquisition unit 212, as the position indicated by the coordinate → ξ that has the maximum similarity C _ξ with the RSSI vector → x that is the observation data. Specify as. The position identifying unit 222 processes the identified position of the terminal device 100 into image information such as a map through the communication interface 202 and outputs the image information to the terminal device 100 that is the acquisition source of the acquisition unit 212 (step S308). ). As a result, the control unit 130 ends the processing of this flowchart.

以上説明した第１の実施形態の位置特定支援装置２００によれば、発信源Ｓから送信される電波を受信する端末装置１００から、複数の時刻において測定された電波の強度に基づく発信源毎データＲＳＳＩ_Ｓｋを取得する取得部と、取得部により取得された発信源毎データＲＳＳＩ_Ｓｋから、ＲＳＳＩベクトル→ｘを生成し、生成したＲＳＳＩベクトル→ｘに対して自己相関行列Γ_→ξを導出する相関行列演算部２１４と、相関行列演算部２１４により導出された自己相関行列から共起した成分を抽出し、抽出した成分を固有値λとした未知数の固有ベクトル→ψ_→ξｋを導出する固有ベクトル導出部２１６と、固有ベクトル導出部２１６により導出された固有ベクトル→ψ_→ξｋを、座標→ξに対応させて記憶部２３０に記憶させる記憶制御部２１８とを備えることにより、複数の測定データをまとめて処理することができる。この結果、第１の実施形態の位置特定支援装置２００は、処理負荷を軽減することができる。 According to the position identification support device 200 of the first exemplary embodiment described above, the data for each transmission source from the terminal device 100 receiving the radio wave transmitted from the transmission source S based on the intensity of the radio wave measured at a plurality of times. Correlation that generates an RSSI vector → x from the acquisition unit that acquires the RSSI _Sk and the data RSSI _{Sk for} each transmission source acquired by the acquisition unit, and that derives the autocorrelation matrix Γ _{→ ξ} for the generated RSSI vector → x A matrix calculation unit 214, and an eigenvector derivation unit 216 that extracts co-occurring components from the autocorrelation matrix derived by the correlation matrix calculation unit 214 and derives an unknown eigenvector → ψ _{→ ξk} with the extracted component as the eigenvalue λ. the eigenvectors _→ ψ → ξk derived by eigenvectors deriving unit 216 in the storage unit 230 in correspondence with the coordinates → xi] storage controller By providing a 218, can be processed together a plurality of measurement data. As a result, the position identifying support device 200 of the first embodiment can reduce the processing load.

また、第１の実施形態の位置特定支援装置２００によれば、自己相関行列から共起した成分（安定した成分）を抽出することにより、次元数ｎのＲＳＳＩベクトル→ｘのデータ量に比して位置特定に用いる際のデータ量を少なくすることができる。この結果、第１の実施形態の位置特定支援装置２００は、処理負荷を更に軽減することができる。 In addition, according to the position identification support device 200 of the first exemplary embodiment, the co-occurring component (stable component) is extracted from the autocorrelation matrix, so that the RSSI vector of the dimension number n is compared with the data amount of x. It is possible to reduce the amount of data used for position identification. As a result, the position identification support device 200 of the first embodiment can further reduce the processing load.

また、第１の実施形態の位置特定支援装置２００によれば、複数（所定数ｎ_→ξ）のＲＳＳＩベクトル→ｘと、フィンガープリントの固有ベクトル→ψ_→ξｋとに基づいて類似度Ｃを算出することにより、精度よく端末装置１００の位置を特定することができる。 Further, according to the position identifying support device 200 of the first embodiment, the similarity C is calculated based on a plurality (predetermined number n _{→ ξ} ) of RSSI vectors → x and fingerprint eigenvectors → ψ _{→ ξk.} As a result, the position of the terminal device 100 can be accurately specified.

（第２の実施形態）
以下、第２の実施形態における位置特定システム１について説明する。ここでは、第１の実施形態との相違点として、端末装置１００から位置特定支援装置２００に対して送信される情報が、固有ベクトル→ψ_→ξｋである場合について説明する。以下、上述した第１の実施形態と共通する機能等についての説明は省略する。 (Second embodiment)
Hereinafter, the position specifying system 1 according to the second embodiment will be described. Here, as a difference from the first embodiment, a case where the information transmitted from the terminal device 100 to the position identification support device 200 is an eigenvector → ψ _{→ ξk} will be described. Hereinafter, description of functions and the like common to the above-described first embodiment will be omitted.

図８は、第２の実施形態における端末装置１００の機能構成の一例を示す図である。図８に示すように、第２の実施形態における端末装置１００は、相関行列演算部１３６と、固有ベクトル導出部１３８とをさらに備える。相関行列演算部１３６と固有ベクトル導出部１３８とは、第１の実施形態における位置特定支援装置２００の相関行列演算部２１４と固有ベクトル導出部２１６とのそれぞれに相当する機能部である。 FIG. 8 is a diagram illustrating an example of a functional configuration of the terminal device 100 according to the second embodiment. As shown in FIG. 8, the terminal device 100 according to the second embodiment further includes a correlation matrix calculation unit 136 and an eigenvector deriving unit 138. The correlation matrix calculation unit 136 and the eigenvector derivation unit 138 are functional units corresponding to the correlation matrix calculation unit 214 and the eigenvector derivation unit 216 of the position identifying support device 200 according to the first embodiment.

（運用処理）
以下、運用過程における、第２の実施形態における端末装置１００側の処理を、ローチャートに則して説明する。図９は、第２の実施形態における端末装置１００の運用処理の流れの一例を示すフローチャートである。 (Operation processing)
Hereinafter, the process on the terminal device 100 side in the second embodiment in the operation process will be described according to a row chart. FIG. 9 is a flowchart showing an example of the flow of operation processing of the terminal device 100 according to the second embodiment.

端末装置１００の制御部１３０は、記憶部１５０に予め記憶されたアプリケーションプログラムを実行し（ステップＳ４００）、例えば、以下の処理をバックグラウンドで実施させる。次に、制御部１３０は、アンテナ１０２を介して発信源Ｓから電波を受信したか否かを判定する（ステップＳ４０２）。 The control unit 130 of the terminal device 100 executes the application program stored in advance in the storage unit 150 (step S400), and executes the following process in the background, for example. Next, the control unit 130 determines whether a radio wave is received from the transmission source S via the antenna 102 (step S402).

ＲＳＳＩデータ生成部１３２は、アンテナ１０２を介して発信源Ｓから電波を受信した場合、発信源Ｓから受信した電波に基づく信号を基に、発信源毎データＲＳＳＩ_Ｓｋを生成する（ステップＳ４０４）。制御部１３０は、アンテナ１０２を介して発信源Ｓから電波を受信していない場合、ステップＳ４０４の処理を行わずに、ステップＳ４０６の処理を行う。 When the radio wave is received from the transmission source S via the antenna 102, the RSSI data generation unit 132 generates the transmission source data RSSI _Sk based on the signal based on the radio wave received from the transmission source S (step S404). When the radio wave is not received from the transmission source S via the antenna 102, the control unit 130 performs the process of step S406 without performing the process of step S404.

次に、制御部１３０は、発信源Ｓからの電波を測定し始めてから、所定時間が経過したか否かを判定する（ステップＳ４０６）。制御部１３０は、所定時間が経過しない場合、ステップＳ４０２の処理に戻る。相関行列演算部１３６は、所定時間が経過した場合、生成した発信源毎データＲＳＳＩ_Ｓｋに基づいて、ＲＳＳＩベクトル→ｘを生成する（ステップＳ４０８）。 Next, the control unit 130 determines whether or not a predetermined time has elapsed since the measurement of the radio wave from the transmission source S was started (step S406). When the predetermined time has not elapsed, the control unit 130 returns to the process of step S402. When the predetermined time has elapsed, the correlation matrix calculation unit 136 generates the RSSI vector → x based on the generated source-specific data RSSI _Sk (step S408).

次に、相関行列演算部１３６は、生成したＲＳＳＩベクトル→ｘの数が所定数ｎ_→ξに達したか否かを判定する（ステップ４１０）。制御部１３０は、ＲＳＳＩベクトル→ｘの数が所定数ｎ_→ξに達しない場合、ステップＳ４０２の処理に戻る。固有ベクトル導出部１３８は、ＲＳＳＩベクトル→ｘの数が所定数ｎ_→ξに達した場合、所定数ｎ_→ξのＲＳＳＩベクトル→ｘに対して、上述した数式（５）に示す自己相関関数式Γ_→ξを適用して、自己相関行列を導出する。 Next, the correlation matrix calculation unit 136 determines whether or not the number of generated RSSI vectors → x has reached a predetermined number n _{→ ξ} (step 410). When the number of RSSI vectors → x does not reach the predetermined number n _{→ ξ} , the control unit 130 returns to the process of step S402. When the number of RSSI vectors → x reaches a predetermined number n _{→ ξ} , the eigenvector deriving unit 138 _responds to the predetermined number n _{→ ξ} of RSSI vectors → x by the autocorrelation function formula Γ shown in the above equation (5). _{→ ξ} is applied to derive the autocorrelation matrix.

次に、固有ベクトル導出部１３８は、相関行列演算部１３６により導出された自己相関行列に基づいて、固有値問題を解き、固有ベクトル→ψ_→ξｋ＃を導出する（ステップＳ４１２）。次に、出力制御部１３４は、アンテナ１０２、変調／復調部１０４、およびＡＤ／ＤＡ変換部１０６を制御して、固有ベクトル導出部１３８により導出された固有ベクトル→ψ_→ξｋ＃を、位置特定支援装置２００に出力する（ステップＳ４１４）。これによって、制御部１３０は、本フローチャートの処理を終了する。 Next, the eigenvector deriving unit 138 solves the eigenvalue problem based on the autocorrelation matrix derived by the correlation matrix computing unit 136, and derives the eigenvector → ψ _{→ ξk} # (step S412). Next, the output control unit 134 controls the antenna 102, the modulation / demodulation unit 104, and the AD / DA conversion unit 106 so that the eigenvector → ψ _{→ ξk} # derived by the eigenvector deriving unit 138 is _transferred to the position identification support device. It outputs to 200 (step S414). As a result, the control unit 130 ends the processing of this flowchart.

以上説明した第２の実施形態の端末装置１００によれば、上述した処理によって、運用時において生成した固有ベクトルψ_→ξｋ＃を位置特定支援装置２００に送信することにより、生成した固有ベクトルψ_→ξｋ＃とフィンガープリントして予め記憶されている固有ベクトル→ψ_→ξｋとの類似度Ｃを位置特定支援装置２００に算出させることができる。この結果、第２の実施形態の端末装置１００は、第１の実施形態と同様の効果を奏する上に、第１の実施形態に比して、通信負荷や位置特定支援装置２００側の処理負担を軽減することができる。 According to the second embodiment of the terminal apparatus 100 described above, the processing described above, by transmitting the eigenvectors _{ψ →} ξk _# generated during operation to location assistance unit 200, generates eigenvectors _{ψ →} ξk # It is possible to cause the position identification support apparatus 200 to calculate the similarity C of the eigenvector → ψ _{→ ξk} that is fingerprinted and stored in advance. As a result, the terminal device 100 of the second embodiment has the same effects as those of the first embodiment, and in addition to the communication load of the first embodiment, the processing load on the position identification support device 200 side. Can be reduced.

（他の実施形態）
以下、他の実施形態（変形例）について説明する。
上述した第１、２の実施形態における位置特定システム１は、測定した発信源毎データＲＳＳＩ_Ｓｋに基づく情報（ＲＳＳＩベクトル→ｘ）を、予め記憶させておいたフィンガープリントの固有ベクトルψ_→ξｋにより張られる空間に線形変換（射影）することによって、端末装置１００の位置を特定するものとして説明したがこれに限られない。例えば、位置特定システム１は、ｋｅｒｎｅｌＰＣＡ等の非線形変換を用いて固有ベクトルψ_→ξｋを導出する主成分分析手法を適用してもよい。 (Other embodiments)
Hereinafter, other embodiments (modifications) will be described.
The position identifying system 1 in the first and second embodiments described above extends the information (RSSI vector → x) based on the measured data RSSI _Sk for each transmission source by the eigenvector ψ _{→ ξk} of the fingerprint stored in advance. Although the position of the terminal device 100 is specified by performing a linear conversion (projection) into the defined space, the present invention is not limited to this. For example, the position identifying system 1 may apply a principal component analysis method for deriving the eigenvector ψ _{→ ξk} using a nonlinear transformation such as kernelPCA.

上述した実施形態における端末装置１００、または位置特定支援装置２００の処理をコンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＦＰＧＡ（Field Programmable Gate Array）等のプログラマブルロジックデバイスを用いて実現されるものであってもよい。 The processing of the terminal device 100 or the position identification support device 200 in the above-described embodiment may be realized by a computer. In that case, a program for realizing this function may be recorded in a computer-readable recording medium, and the program recorded in this recording medium may be read by a computer system and executed. The “computer system” mentioned here includes an OS and hardware such as peripheral devices. Further, the “computer-readable recording medium” refers to a portable medium such as a flexible disk, a magneto-optical disk, a ROM, a CD-ROM, or a storage device such as a hard disk built in a computer system. Further, the "computer-readable recording medium" means to hold a program dynamically for a short time like a communication line when transmitting the program through a network such as the Internet or a communication line such as a telephone line. In this case, a volatile memory inside a computer system that serves as a server or a client in that case may hold a program for a certain period of time. Further, the program may be for realizing a part of the functions described above, and may be a program that can realize the functions described above in combination with a program already recorded in a computer system, It may be realized using a programmable logic device such as FPGA (Field Programmable Gate Array).

以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。 Although the embodiment of the present invention has been described in detail above with reference to the drawings, the specific configuration is not limited to this embodiment, and includes a design and the like within a range not departing from the gist of the present invention.

１…位置特定システム、１００…端末装置、１０２…アンテナ、１０４…変調／復調部、１０６…ＡＤ／ＤＡ変換部、１０８…入力部、１１０…出力部、１３０…制御部、１３２…データ生成部、１３４…出力制御部、１３６…相関行列演算部、１３８…固有ベクトル導出部、１５０…記憶部、２００…位置特定支援装置、２０２…通信インターフェース、２０４…入力部、２０６…表示部、２１０…制御部、２１２…取得部、２１４…相関行列演算部、２１６…固有ベクトル導出部、２１８…記憶制御部、２２０…類似度演算部、２２２…特定部、２３０…記憶部、ＮＷ…ネットワーク DESCRIPTION OF SYMBOLS 1 ... Position identification system, 100 ... Terminal device, 102 ... Antenna, 104 ... Modulation / demodulation section, 106 ... AD / DA conversion section, 108 ... Input section, 110 ... Output section, 130 ... Control section, 132 ... Data generation section , 134 ... Output control unit, 136 ... Correlation matrix calculation unit, 138 ... Eigen vector derivation unit, 150 ... Storage unit, 200 ... Position identification support device, 202 ... Communication interface, 204 ... Input unit, 206 ... Display unit, 210 ... Control Part, 212 ... Acquisition part, 214 ... Correlation matrix calculation part, 216 ... Eigenvector derivation part, 218 ... Storage control part, 220 ... Similarity degree calculation part, 222 ... Specification part, 230 ... Storage part, NW ... Network

Claims

From a terminal device that receives radio waves transmitted from a transmission source, an acquisition unit that acquires information based on the intensity of the radio waves measured at a plurality of measurement times by the terminal device,
Based on the spatial correlation method, by extracting the correlation matrix of the intensity of the radio waves at a plurality of measurement times acquired by the acquisition unit, an extraction unit for extracting an eigenvector indicating a stable component of the intensity of the radio waves at the measurement time. When,
A storage control unit that stores the eigenvector extracted as the stable component by the extraction unit and the coordinate vector indicating the position of the terminal device known in advance in association with each other in the storage unit as correspondence information,
An identification unit that identifies the position of the terminal device that is the acquisition source of the information based on the intensity of the radio wave, based on the correspondence information stored in the storage unit and the information based on the intensity of the radio wave acquired by the acquisition unit. And
The specifying unit is
When information indicating the strength of the radio wave is acquired by the acquisition unit at an operation time after the measurement time, a higher predetermined value with a larger eigenvalue among the plurality of eigenvectors stored as the correspondence information in the storage unit. Extract a number of said eigenvectors,
Calculate an inner product of each of the extracted predetermined number of eigenvectors and a radio wave vector indicating the strength of the radio wave at the operating time,
The position indicated by the coordinate vector where the sum of the inner products of the eigenvector and the radio wave vector is maximum is specified as the position of the terminal device at the operation time.
Position specific support device.

An identification unit that identifies the position of the terminal device that is the acquisition source of the information based on the intensity of the radio wave, based on the correspondence information stored in the storage unit and the information based on the intensity of the radio wave acquired by the acquisition unit. Further comprising,
The position identification support device according to claim 1.

The extraction unit extracts the eigenvector indicating the stable component by solving the correlation matrix,
The storage control unit stores the eigenvector extracted by the extraction unit in the storage unit as correspondence information in association with a previously known position of the terminal device,
The position identification support device according to claim 2.

The specifying unit calculates the similarity between the eigenvector included in the correspondence information stored in the storage unit and a predetermined vector included in the information based on the intensity of the radio wave acquired by the acquiring unit, Based on the eigenvector having the maximum degree of similarity with the vector and the correspondence information stored in the storage unit, the position of the terminal device that is the acquisition source of the information based on the intensity of the radio wave is specified.
The position identification support device according to claim 3.

From the information based on the intensity of the radio wave measured by the terminal device at a plurality of times, from the terminal device that receives the radio wave transmitted from the transmission source, the acquisition unit determines the intensity of the radio wave extracted by the terminal device. Acquire information indicating stable components,
The specifying unit, based on the correspondence information stored in the storage unit and the information indicating the stable component acquired by the acquisition unit, determines the position of the terminal device that is the acquisition source of the information based on the intensity of the radio wave. Identify,
The position identification support device according to claim 2.

From a terminal device that receives radio waves transmitted from a transmission source, obtain information based on the intensity of the radio waves measured at a plurality of measurement times by the terminal device,
Based on the spatial correlation method, by solving the correlation matrix of the intensity of the radio wave at the plurality of acquired measurement time, to extract an eigenvector indicating a stable component of the intensity of the radio wave at the measurement time ,
The eigenvector extracted as the stable component and the coordinate vector indicating the position of the terminal device known in advance are associated with each other and stored in the storage unit as correspondence information ,
Based on the correspondence information stored in the storage unit and the information based on the acquired intensity of the radio wave, the position of the terminal device that is the acquisition source of the information based on the intensity of the radio wave is specified,
When the information indicating the strength of the radio wave is acquired at the operation time after the measurement time, among a plurality of the eigenvectors stored as the correspondence information in the storage unit, the eigenvalues having a large upper predetermined number of the eigenvectors are selected. Extract and
Calculate an inner product of each of the extracted predetermined number of eigenvectors and a radio wave vector indicating the strength of the radio wave at the operating time,
The position indicated by the coordinate vector where the sum of the inner products of the eigenvector and the radio wave vector is maximum is specified as the position of the terminal device at the operation time.
Location method.

A computer of the position identification support device according to claim 1,
From the terminal device receiving the radio wave transmitted from the transmission source, based on the information based on the intensity of the radio wave measured by the terminal device at a plurality of times, the stability of the radio wave intensity at the measurement time extracted by the terminal device obtaining an eigenvector the ingredients,
Storing the eigenvector extracted as the stable component and the coordinate vector indicating the position of the terminal device known in advance in the storage unit as correspondence information.
The corresponding information stored in the storage unit, based on the eigenvectors the acquired, identifying the location of the terminal device is a retrieve the previous period eigenvectors,
When the eigenvector is acquired at an operating time after the measurement time, of the plurality of eigenvectors stored as the correspondence information in the storage unit, extracting a predetermined upper number of eigenvectors having a large eigenvalue,
Calculating an inner product of each of the extracted predetermined number of eigenvectors and the eigenvector at the operating time,
Specifying the position indicated by the coordinate vector having the maximum sum of inner products of the eigenvectors as the position of the terminal device at the operating time,
Positioning program for executing .