JP2020141316A - Radio wave identification device, radio wave identification method, and radio wave environment visualization system - Google Patents

Radio wave identification device, radio wave identification method, and radio wave environment visualization system Download PDF

Info

Publication number
JP2020141316A
JP2020141316A JP2019036507A JP2019036507A JP2020141316A JP 2020141316 A JP2020141316 A JP 2020141316A JP 2019036507 A JP2019036507 A JP 2019036507A JP 2019036507 A JP2019036507 A JP 2019036507A JP 2020141316 A JP2020141316 A JP 2020141316A
Authority
JP
Japan
Prior art keywords
time
received signal
signal sequence
radio wave
packet
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2019036507A
Other languages
Japanese (ja)
Other versions
JP7135938B2 (en
Inventor
沢田 健介
Kensuke Sawada
健介 沢田
達也 菊月
Tatsuya Kikuzuki
達也 菊月
真 濱湊
Makoto Hamaminato
真 濱湊
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.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2019036507A priority Critical patent/JP7135938B2/en
Publication of JP2020141316A publication Critical patent/JP2020141316A/en
Application granted granted Critical
Publication of JP7135938B2 publication Critical patent/JP7135938B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

To compensate for individual differences in time information at each receiving node.SOLUTION: From among first and second received signal sequences received by first and second receiving nodes, while the scale of a time axis of the second received signal sequence is expanded and contracted at the time expansion/contraction rate, the correlation value with the first received signal sequence is calculated, the second received signal sequence is offset on the time axis by a time offset to determine the expansion/contraction rate and the time offset at which the correlation value becomes the maximum value, packets of the first and second received signal sequences are compared with each other by applying the time expansion/contraction rate and the time offset at which the correlation value becomes the maximum value to time information of the second received signal sequence, and packets having similar transmission start times of the packets, packet lengths, and signal strengths are associated with each other on a one-to-one basis, and when the maximum correlation value exceeds a threshold, the associated packet outputs a radio wave identification result determined to be a received signal sequence of radio waves received by the first and second receiving nodes from the same transmission source.SELECTED DRAWING: Figure 5

Description

本発明は、電波同定装置、電波同定方法及び電波環境可視化システムに関する。 The present invention relates to a radio wave identification device, a radio wave identification method, and a radio wave environment visualization system.

例えば、無線システムを新規に構築する際には、干渉や妨害等を既存の無線システムに与えたり、既存の無線システムから受けたりしないことを確認するために、電波の利用状況を事前に調査する。また、無線システムを運用中に、通信不達や不具合が発生した際にも、不具合の原因を特定するために電波の利用状況を調査する。 For example, when constructing a new wireless system, the usage status of radio waves is investigated in advance in order to confirm that interference and interference are not given to the existing wireless system or received from the existing wireless system. .. In addition, even if communication is not delivered or a problem occurs during operation of the wireless system, the usage status of radio waves will be investigated in order to identify the cause of the problem.

電波の利用状況を調査するために、電波の利用状況を可視化する電波環境可視化技術が提案されている。電波環境可視化技術を活用することで、場所毎の電波の受信信号強度(RSSI:Received Signal Strength Indicator)やトラフィックのレベルを例えば色で表現したヒートマップを作成することができる。また、送信局の位置情報やID情報をヒートマップ上に表示することで、周辺の電波環境の調査や不具合の原因特定に活用できる。 In order to investigate the usage status of radio waves, a radio wave environment visualization technology that visualizes the usage status of radio waves has been proposed. By utilizing the radio wave environment visualization technology, it is possible to create a heat map in which the received signal strength (RSSI: Received Signal Strength Indicator) and the traffic level of radio waves for each location are expressed in colors, for example. In addition, by displaying the location information and ID information of the transmitting station on the heat map, it can be used for investigating the surrounding radio wave environment and identifying the cause of the problem.

ヒートマップを作成するためには、空間的に離れた複数の場所に受信ノードを設置し、各受信ノードで受信された電波(例えば、パケット)の夫々が同一の送信局から送信されたものであるか否かを判定したり、送信局の位置を推定したりする。各受信ノードで受信された電波が同一の送信局から送信されたものであるか否かを判定する処理とは、例えば複数の異なる受信ノードで受信された同一の電波を同定する処理である。 In order to create a heat map, receiving nodes are installed in multiple spatially separated locations, and the radio waves (for example, packets) received by each receiving node are transmitted from the same transmitting station. It determines whether or not there is one, and estimates the position of the transmitting station. The process of determining whether or not the radio wave received by each receiving node is transmitted from the same transmitting station is, for example, a process of identifying the same radio wave received by a plurality of different receiving nodes.

しかし、複数の場所に個別に設置された各受信ノードの内部時計は、互いに独立しているため、複数の受信ノードの内部時計が受信電波に付与する時刻情報(例えば、タイムスタンプ)には差分が発生して一致しない。このため、複数の受信ノードで測定されたRSSIの時系列データを用いて複数の電波を関連付けることは困難である。 However, since the internal clocks of the receiving nodes individually installed in a plurality of locations are independent of each other, the time information (for example, time stamp) given to the received radio waves by the internal clocks of the plurality of receiving nodes is different. Occurs and does not match. Therefore, it is difficult to associate a plurality of radio waves using RSSI time series data measured by a plurality of receiving nodes.

時刻情報に差分が発生する仕組みを図1及び図2と共に説明する。図1は、電波環境可視化システムの一例を示す図である。図1に示す電波環境可視化システム1は、送信局2と、受信ノード3A,3Bと、サーバ4とを有する。この例では、TIMETXは送信局2における送信時刻情報、Dist#1は送信局2と受信ノード3Aとの間の距離、TIMERX#1は受信ノード3Aにおける受信時刻情報、TIMERX#2は受信ノード3Bにおける受信時刻情報、Dist#2は送信局2と受信ノード3Bとの間の距離を示す。また、TIME#1は、受信ノード3Aが出力した受信時刻情報、TIME#2は受信ノード3Bが出力した受信時刻情報を示す。ここで、光速をcで表すと、c(TIMERX#1−TIMETX)=Dist#1,c(TIMERX#2−TIMETX)=Dist#2なる関係が成り立つ。しかし、個々の受信ノード3A,3Bの処理速度等に差があることから、TIMERX#1≠TIME#1、且つ、TIMERX#2≠TIME#2である。 The mechanism by which the difference occurs in the time information will be described with reference to FIGS. 1 and 2. FIG. 1 is a diagram showing an example of a radio wave environment visualization system. The radio wave environment visualization system 1 shown in FIG. 1 has a transmitting station 2, receiving nodes 3A and 3B, and a server 4. In this example, TIME TX is the transmission time information at the transmitting station 2, Dist # 1 is the distance between the transmitting station 2 and the receiving node 3A, TIME RX # 1 is the receiving time information at the receiving node 3A, and TIME RX # 2 is. The reception time information in the reception node 3B, Dist # 2, indicates the distance between the transmission station 2 and the reception node 3B. Further, TIME # 1 indicates the reception time information output by the reception node 3A, and TIME # 2 indicates the reception time information output by the reception node 3B. Here, when the speed of light is expressed by c, the relationship c (TIME RX # 1 −TIME TX ) = Dist # 1 and c (TIME RX # 2 −TIME TX ) = Dist # 2 holds. However, since there is a difference in the processing speeds of the individual receiving nodes 3A and 3B, TIME RX # 1 ≠ TIME # 1 and TIME RX # 2 ≠ TIME # 2.

図2は、上記の如きタイムスタンプの差分が発生するメカニズムを説明する図である。図2は、正解受信信号と正解受信時刻TIMERX#1,TIMERX#2、受信ノード3Aの受信時刻情報(以下、「タイムスタンプ」とも言う)TIME#1、受信ノード3BのタイムスタンプTIME#2、及びサーバ4に届いた受信ノード3A,3Bの信号を示す。図2中、tは正解受信時刻、tは受信ノード3Aにおける時刻、tは受信ノード3Bにおける時刻、及びtはサーバ4における時刻を、夫々同じ任意単位の時間軸で示す。図2において、受信ノード3A,3Bの互いに独立した内部時計が決定するパケット受信時刻の決定結果には差分が発生する。このため、同一の信号に付与されたタイムスタンプTIME#1, TIME#2と、実際に当該信号が受信された時刻TIMERX#1, TIMERX#2との間に、例えば10−3ppmオーダの誤差が生じる。この結果、サーバ4に届いた信号の時刻情報に差分が発生する。そこで、このようなメカニズムで発生するタイムスタンプの差分を吸収して補償することが望ましい。 FIG. 2 is a diagram illustrating a mechanism for generating the above-mentioned time stamp difference. FIG. 2 shows the correct answer reception signal and the correct answer reception time TIME RX # 1, TIME RX # 2 , the reception time information of the reception node 3A (hereinafter, also referred to as “time stamp”) TIME # 1, and the time stamp TIME # of the reception node 3B. The signals of the receiving nodes 3A and 3B that arrived at 2 and the server 4 are shown. In FIG. 2, t indicates the correct answer reception time, t A indicates the time at the receiving node 3A, t B indicates the time at the receiving node 3B, and t S indicates the time at the server 4 on the same arbitrary unit time axis. In FIG. 2, a difference occurs in the determination result of the packet reception time determined by the internal clocks of the reception nodes 3A and 3B that are independent of each other. Therefore, there is an error of, for example, 10 -3 ppm between the time stamps TIME # 1 and TIME # 2 given to the same signal and the time TIMERX # 1 and TIMERX # 2 when the signal was actually received. Occurs. As a result, a difference occurs in the time information of the signal that arrives at the server 4. Therefore, it is desirable to absorb and compensate for the difference in time stamps generated by such a mechanism.

複数の時刻同期していない受信機において、基準となる送信機から送信された信号の受信フレームから生成したテンプレート波形データと、実際の受信波形データの相関値が最大になる点を求めるシステムが提案されている(例えば、特許文献1参照)。この提案されているシステムでは、実際の受信波形データの相関値が最大になるときのタイムスタンプを記録して、各受信機における受信時間差を求める。更に、受信機より上位の位置計算機で送信機の位置を推定し、複数の送受信機に対しても同様にして位置を推定する。しかし、この提案されているシステムは、送信機の測位を目的とした構成を有する。このため、各受信機で受信された信号同士を同定できるような高い精度で時刻情報を同期させる機能は有していない。つまり、各受信機から上位の位置計算機までの異なる距離に応じた異なる遅延が生じるため、位置計算機は、例えば各受信機で受信されたパケット同士を紐付けできるような高い精度でタイムスタンプの個体差を補償することはできない。 A system that finds the point where the correlation value between the template waveform data generated from the reception frame of the signal transmitted from the reference transmitter and the actual reception waveform data is maximized in multiple receivers that are not time-synchronized is proposed. (See, for example, Patent Document 1). In this proposed system, the time stamp when the correlation value of the actual received waveform data is maximized is recorded, and the reception time difference at each receiver is obtained. Further, the position of the transmitter is estimated by a position computer higher than the receiver, and the positions are estimated in the same manner for a plurality of transmitters and receivers. However, this proposed system has a configuration intended for transmitter positioning. Therefore, it does not have a function of synchronizing time information with high accuracy so that signals received by each receiver can be identified. In other words, since different delays occur according to different distances from each receiver to the higher-level position calculator, the position calculator is an individual time stamp with high accuracy that can link packets received by each receiver, for example. The difference cannot be compensated.

特開2014−085205号公報Japanese Unexamined Patent Publication No. 2014-0852005 特開平11−308658号公報JP-A-11-308658 特開2014−045283号公報JP-A-2014-045283

従来技術では、各受信ノードにおける時刻情報の個体差を補償することは難しい。 With the prior art, it is difficult to compensate for individual differences in time information at each receiving node.

そこで、1つの側面では、各受信ノードにおける時刻情報の個体差を補償することができる電波同定装置、電波同定方法及び電波環境可視化システムを提供することを目的とする。 Therefore, in one aspect, it is an object of the present invention to provide a radio wave identification device, a radio wave identification method, and a radio wave environment visualization system capable of compensating for individual differences in time information in each receiving node.

1つの案によれば、第1の受信ノードが受信した第1の受信信号系列と、第2の受信ノードが受信した第2の受信信号系列のうち、前記第2の受信信号系列の時間軸のスケールを時間伸縮率で伸縮しながら前記第1の受信信号系列との相関値を算出し、前記第2の受信信号系列を時間軸上で時間オフセットだけオフセットして前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する決定手段と、前記相関値が最大値になる前記時間伸縮率と前記時間オフセットを前記第2の受信信号系列の時刻情報に適用することで、前記第1の受信信号系列及び前記第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付ける紐付け手段と、前記最大値になる前記相関値が閾値を超えている場合に、紐付けられたパケットは前記第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する判定手段と、を備えた電波同定装置が提供される。 According to one proposal, the time axis of the second received signal sequence among the first received signal sequence received by the first receiving node and the second received signal sequence received by the second receiving node. The correlation value with the first received signal sequence is calculated while expanding and contracting the scale according to the time expansion / contraction rate, and the second received signal sequence is offset by the time offset on the time axis so that the correlation value becomes the maximum value. By applying the time expansion / contraction rate and the time offset for determining the time expansion / contraction rate and the time offset, and the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value to the time information of the second received signal series Each packet of the first received signal sequence and the second received signal sequence is compared with each other, and packets having similar transmission start times, packet lengths, and signal strengths of each packet are one-to-one. When the linking means to be linked with and the correlation value that becomes the maximum value exceeds the threshold value, the linked packet is the radio wave received by the first and second receiving nodes from the same source. Provided is a radio wave identification device including a determination means for outputting a radio wave identification result determined to be a received signal sequence.

一態様によれば、各受信ノードにおける時刻情報の個体差を補償することができる。 According to one aspect, individual differences in time information at each receiving node can be compensated.

電波環境可視化システムの一例を示す図である。It is a figure which shows an example of the radio wave environment visualization system. タイムスタンプの差分が発生するメカニズムを説明する図である。It is a figure explaining the mechanism which the difference of a time stamp occurs. 第1実施例における電波環境可視化システムの一例を示す図である。It is a figure which shows an example of the radio wave environment visualization system in 1st Example. パーソナルコンピュータのハードウェア構成の一例を示すブロック図である。It is a block diagram which shows an example of the hardware configuration of a personal computer. 第1実施例における電波同定処理の一例を説明するフローチャートである。It is a flowchart explaining an example of the radio wave identification process in 1st Example. ログデータ読込処理の一例を説明するフローチャートである。It is a flowchart explaining an example of a log data reading process. ログデータクリッピング処理の一例を説明するフローチャートである。It is a flowchart explaining an example of a log data clipping process. ログデータクリッピング処理を4つのケースについて説明する図である。It is a figure explaining the log data clipping process for four cases. ログデータ相関算出処理の一例を説明するフローチャート(その1)である。It is a flowchart (the 1) explaining an example of a log data correlation calculation process. ログデータ相関算出処理の一例を説明するフローチャート(その2)である。It is a flowchart (the 2) explaining an example of a log data correlation calculation process. 伸縮率とオフセット決定処理の一例を説明するフローチャートである。It is a flowchart explaining an example of expansion / contraction ratio and offset determination processing. タイムスタンプの個体差の補償を説明する図である。It is a figure explaining the compensation of the individual difference of a time stamp. パケット紐付け処理の一例を説明するフローチャートである。It is a flowchart explaining an example of a packet association processing. 正常なパケットの紐付け結果の一例を説明する図である。It is a figure explaining an example of the association result of a normal packet. 誤った結果を含むパケットの紐付け結果の一例を説明する図である。It is a figure explaining an example of the association result of the packet including an erroneous result. パケットの紐付け誤りの訂正を説明する図である。It is a figure explaining the correction of the packet association error. 電波同定結果出力処理の一例を説明するフローチャートである。It is a flowchart explaining an example of radio wave identification result output processing. 第2実施例における電波環境可視化システムの一例を示す図である。It is a figure which shows an example of the radio wave environment visualization system in 2nd Example. 第2実施例における電波同定処理の一例を説明するフローチャートである。It is a flowchart explaining an example of the radio wave identification process in 2nd Example. ID組付け結果出力処理の一例を説明するフローチャートである。It is a flowchart explaining an example of ID assembly result output processing. 第3実施例における電波環境可視化システムの一例を示す図である。It is a figure which shows an example of the radio wave environment visualization system in 3rd Example. 第3実施例における電波同定処理の一例を説明するフローチャートである。It is a flowchart explaining an example of radio wave identification processing in 3rd Example. 電波環境可視化処理の一例を説明するフローチャートである。It is a flowchart explaining an example of a radio wave environment visualization process. ヒートマップの一例を示す図である。It is a figure which shows an example of a heat map.

以下に、開示の電波同定装置、電波同定方法及び電波環境可視化システムの各実施例を図面と共に説明する。 Hereinafter, examples of the disclosed radio wave identification device, radio wave identification method, and radio wave environment visualization system will be described with drawings.

(第1実施例)
図3は、第1実施例における電波環境可視化システムの一例を示す図である。図3に示す電波環境可視化システム11−1は、送信局12と、受信モジュール13−1,13−2と、サーバ14と、パーソナルコンピュータ(PC:Personal Computer)15を有する。送信局12は、送信源の一例であり、信号系列のパケット等をアンテナを介して無線送信する周知の構成を有し、例えば周知の送信機により形成可能である。各受信モジュール13−1,13−2は、受信ノードの一例であり、送信局12から送信された信号系列のパケット等をアンテナを介して無線受信する周知の構成を有し、例えば内部時計を有する周知の受信機により形成可能である。この例では、受信モジュール13−1は、受信信号系列のパケットから測定した受信信号強度RSSIと、受信信号系列から検出した送信局12のID情報とを、時刻情報と共にサーバ14に送信する機能を備えた周知のWi−Fi(登録商標)モジュールである。Wi−Fiモジュールは、受信信号系列のパケットから周知の方法で受信信号強度RSSIを測定し、受信信号系列パケットから周知の方法で送信局12のID情報(以下、単に「ID」とも言う)を検出する。一方、受信モジュール13−2は、受信信号系列のパケットから周知の方法で測定した受信信号強度RSSIを、時刻情報と共にサーバ14に送信する機能を備えた周知のSDR(Software Defined Radio)モジュールである。
(First Example)
FIG. 3 is a diagram showing an example of the radio wave environment visualization system according to the first embodiment. The radio wave environment visualization system 11-1 shown in FIG. 3 includes a transmitting station 12, receiving modules 13-1 and 13-2, a server 14, and a personal computer (PC) 15. The transmission station 12 is an example of a transmission source, has a well-known configuration for wirelessly transmitting a packet of a signal series or the like via an antenna, and can be formed by, for example, a well-known transmitter. Each of the receiving modules 13-1 and 13-2 is an example of a receiving node, and has a well-known configuration of wirelessly receiving a packet or the like of a signal series transmitted from the transmitting station 12 via an antenna, for example, an internal clock. It can be formed by a well-known receiver. In this example, the receiving module 13-1 has a function of transmitting the received signal strength RSSI 1 measured from the packet of the received signal series and the ID information of the transmitting station 12 detected from the received signal series to the server 14 together with the time information. It is a well-known Wi-Fi (registered trademark) module provided with. The Wi-Fi module measures the received signal strength RSSI 1 from the received signal series packet by a well-known method, and the ID information of the transmitting station 12 (hereinafter, also simply referred to as “ID”) from the received signal series packet by a well-known method. Is detected. On the other hand, the receiving module 13-2 is a well-known SDR (Software Defined Radio) module having a function of transmitting the received signal strength RSSI 2 measured by a well-known method from the packet of the received signal series to the server 14 together with the time information. is there.

SDRモジュールは、Wi−Fi以外の規格の信号も受信して分析できる。このため、受信モジュールにWi−FiモジュールとSDRモジュールの2種類の受信モジュールを用いることで、例えばWi−Fiモジュールの受信信号に干渉や妨害がある場合、SDRモジュールの受信信号を用いて干渉や妨害の原因を特定できる。しかし、電波環境可視化システムが有する複数の受信モジュールは、異なる種類の受信モジュールに限定されず、同一種類の受信モジュールであっても良い。 The SDR module can also receive and analyze signals of standards other than Wi-Fi. Therefore, by using two types of receiving modules, a Wi-Fi module and an SDR module, for example, if there is interference or interference in the received signal of the Wi-Fi module, the reception signal of the SDR module may be used for interference or interference. The cause of the interference can be identified. However, the plurality of receiving modules included in the radio wave environment visualization system are not limited to different types of receiving modules, and may be the same type of receiving modules.

サーバ14は、受信モジュール13−1からの受信信号強度RSSI、ID、及び時刻情報と、受信モジュール13−2からの受信信号強度RSSI及び時刻情報とを受信可能な汎用コンピュータであり、サーバ14の構成は特に限定されない。この例では、サーバ14は、電波同定装置の一例を形成可能である。PC15は、例えばコマンド(Command)をサーバ14へ送信し、サーバ14から後述する電波同定結果、ID紐付け結果等を受信可能であれば、PC15の構成は特に限定されない。PC15が発行するコマンドは、後述する電波同定処理、電波環境可視化処理等を要求するコマンド等を含む。 The server 14 is a general-purpose computer capable of receiving the received signal strength RSSI 1 , ID, and time information from the receiving module 13-1, and the received signal strength RSSI 2 and time information from the receiving module 13-2. The configuration of 14 is not particularly limited. In this example, the server 14 can form an example of a radio wave identification device. The configuration of the PC 15 is not particularly limited as long as the PC 15 can transmit, for example, a command to the server 14 and receive the radio wave identification result, the ID association result, and the like described later from the server 14. The command issued by the PC 15 includes a command for requesting a radio wave identification process, a radio wave environment visualization process, and the like, which will be described later.

図4は、PCのハードウェア構成の一例を示すブロック図である。図4に示すPC15は、バス156により接続されたCPU(Central Processing Unit)151、メモリ152、インタフェース(I/F)153、入力装置154、及び表示装置155を有する。CPU151は、PC15全体の制御を司るプロセッサの一例である。CPU151は、メモリ152に格納されたプログラムを実行して、各種処理を実行可能である。後述するように、CPU151が実行する各種処理は、電波同定処理、電波環境可視化処理等を含む。メモリ152は、CPU151が実行するプログラム、及びログデータを含む各種データを記憶する記憶装置の一例である。記憶装置は、CD−ROM(Compact Disk-Read Only Memory)、DVD(Digital Versatile Disk)、USB(Universal Serial Bus)メモリ等の可搬型記録媒体、フラッシュメモリ等の半導体メモリ等であっても良い。また、記憶装置は、プログラムを記憶したコンピュータ読み取り可能な記憶媒体であっても良い。更に、記憶装置は、複数の記憶装置で形成されていても良く、この場合、PC15に外部接続された記憶装置を含んでも良い。 FIG. 4 is a block diagram showing an example of the hardware configuration of the PC. The PC 15 shown in FIG. 4 has a CPU (Central Processing Unit) 151, a memory 152, an interface (I / F) 153, an input device 154, and a display device 155 connected by a bus 156. The CPU 151 is an example of a processor that controls the entire PC 15. The CPU 151 can execute various processes by executing the program stored in the memory 152. As will be described later, various processes executed by the CPU 151 include radio wave identification processing, radio wave environment visualization processing, and the like. The memory 152 is an example of a storage device that stores various data including a program executed by the CPU 151 and log data. The storage device may be a portable recording medium such as a CD-ROM (Compact Disk-Read Only Memory), a DVD (Digital Versatile Disk), or a USB (Universal Serial Bus) memory, a semiconductor memory such as a flash memory, or the like. Further, the storage device may be a computer-readable storage medium that stores the program. Further, the storage device may be formed of a plurality of storage devices, and in this case, the storage device may include a storage device externally connected to the PC 15.

I/F153は、PC15とサーバ14との間の情報やコマンド等の授受を司る。入力装置154は、PC15の操作者が操作するキーボード等で形成可能であり、PC15に情報やコマンド等を入力するのに用いられる。表示装置155は、操作者に対してメッセージや処理結果等を表示するのに用いられる。表示装置155が表示する処理結果は、電波同定結果や、電波環境可視化結果の一例であるRSSIのヒートマップ等を含んでも良い。なお、入力装置154と表示装置155とが一体的に設けられたタッチパネルを用いても良い。 The I / F 153 controls the exchange of information, commands, and the like between the PC 15 and the server 14. The input device 154 can be formed by a keyboard or the like operated by the operator of the PC 15, and is used for inputting information, a command, or the like into the PC 15. The display device 155 is used to display a message, a processing result, or the like to the operator. The processing result displayed by the display device 155 may include a radio wave identification result, an RSSI heat map which is an example of a radio wave environment visualization result, and the like. A touch panel in which the input device 154 and the display device 155 are integrally provided may be used.

サーバ14の構成は、図4に示すPC15の構成と同様であっても良いため、サーバ14の構成の図示及び説明は省略する。サーバ14のCPUが実行する各種処理は、電波同定処理等を含む。 Since the configuration of the server 14 may be the same as the configuration of the PC 15 shown in FIG. 4, the illustration and description of the configuration of the server 14 will be omitted. Various processes executed by the CPU of the server 14 include radio wave identification processes and the like.

図5は、第1実施例における電波同定処理の一例を説明するフローチャートである。図5に示す電波同定処理は、PC15からの電波同定処理を要求するコマンドに応答して、サーバ14のCPUが、サーバ14のメモリが記憶するプログラムを実行することで実行可能である。PC15は、操作者による入力装置154の操作に応答して、電波同定処理を要求するコマンドを発行する。サーバ14のメモリは、ファイル、配列、変数等も格納する。 FIG. 5 is a flowchart illustrating an example of the radio wave identification process in the first embodiment. The radio wave identification process shown in FIG. 5 can be executed by the CPU of the server 14 executing a program stored in the memory of the server 14 in response to a command requesting the radio wave identification process from the PC 15. The PC 15 issues a command requesting the radio wave identification process in response to the operation of the input device 154 by the operator. The memory of the server 14 also stores files, arrays, variables, and the like.

図5において、ステップS1では、サーバ14のCPUが、受信モジュール13−1,13−2のログデータを読み込むログデータ読込処理を行う。サーバ14のCPUは、受信モジュール13−1から受信した受信信号強度RSSI、ID、及び時刻情報と、受信モジュール13−2から受信した受信信号強度RSSI及び時刻情報とを、夫々のログデータとしてサーバ14のメモリに格納する。 In FIG. 5, in step S1, the CPU of the server 14 performs a log data reading process for reading the log data of the receiving modules 13-1 and 13-2. The CPU of the server 14 records the received signal strength RSSI 1 , ID, and time information received from the receiving module 13-1 and the received signal strength RSSI 2 and time information received from the receiving module 13-2, respectively. Is stored in the memory of the server 14.

図6は、ログデータ読込処理の一例を説明するフローチャートである。図6において、ステップS11では、サーバ14のCPUが、サーバ14のメモリに格納された受信モジュール13−1のログデータファイルを開ける(以下「オープンする」とも言う)。ステップS12では、サーバ14のCPUが、受信モジュール13−1のログデータファイルを配列RS1に格納する。図6の右側上部に示すように、配列RS1には、例えば時、分、及び秒を表すタイムスタンプ(hh:mm:ss)、パケット長(μsec)、及びRSSI(dBm)が、各パケットのインデックスRS1-1〜RS1-M(Mは2以上の自然数)に対して格納されている。タイムスタンプは、受信モジュールの内部時計が周知の方法で受信電波に付与する時刻情報の一例である。ステップS13では、サーバ14のCPUが、配列RS1のパケット数Mを変数NumPktRS1に格納する。 FIG. 6 is a flowchart illustrating an example of log data reading processing. In FIG. 6, in step S11, the CPU of the server 14 opens the log data file of the receiving module 13-1 stored in the memory of the server 14 (hereinafter, also referred to as “open”). In step S12, the CPU of the server 14 stores the log data file of the receiving module 13-1 in the array RS1. As shown in the upper right part of FIG. 6, the array RS1 contains, for example, a time stamp (hh: mm: ss) representing hours, minutes, and seconds, a packet length (μsec), and RSSI (dBm) of each packet. It is stored for indexes RS1-1 to RS1-M (M is a natural number of 2 or more). The time stamp is an example of time information given to the received radio wave by the internal clock of the receiving module by a well-known method. In step S13, the CPU of the server 14 stores the number of packets M of the array RS1 in the variable NumPktRS1.

ステップS14では、サーバ14のCPUが、受信モジュール13−2のログデータファイルをオープンする。ステップS15では、サーバ14のCPUが、受信モジュール13−2のログデータファイルを配列RS2に格納する。図6の右側中部に示すように、配列RS2には、タイムスタンプ、パケット長、及びRSSIが、各インデックスRS2-1〜RS2-Mに対して格納されている。ステップS16では、サーバ14のCPUが、配列RS2のパケット数Mを変数NumPktRS2に格納する。 In step S14, the CPU of the server 14 opens the log data file of the receiving module 13-2. In step S15, the CPU of the server 14 stores the log data file of the receiving module 13-2 in the array RS2. As shown in the central part on the right side of FIG. 6, the time stamp, packet length, and RSSI are stored in the array RS2 for each index RS2-1 to RS2-M. In step S16, the CPU of the server 14 stores the number of packets M of the array RS2 in the variable NumPktRS2.

ステップS17では、サーバ14のCPUが、動作設定ファイルをオープンする。図6の右側下部に便宜上一部を示すように、動作設定ファイルは、RSSIの差を補償するためのRSSIのオフセットRSSI_Offset(dB)、時間軸のスケールを伸縮すると共に時間軸上のオフセットを加えるためのガード時間Guard_Time(sec)等を含む。この例では、動作設定ファイルは、時間オフセットtoffsetの最小値t_begin(sec)、時間オフセットtoffsetの最大値t_end(sec)、及びループ回数num_t(回)を含む。また、動作設定ファイルは、時間伸縮率(以下、単に「伸縮率」とも言う)kの最小値k_begin、伸縮率kの最大値k_end、ループ回数num_k(回)、RSSIのオフセットRSSI_Offset(dB)、ガード時間Guard_Time(sec)等を含む。動作設定ファイルは更に、パケット間相関算出用の過去のサーチウィンドウの時間長search_window_past(sec)、及びパケット間相関算出用の未来のサーチウィンドウの時間長search_window_future(sec)を含む。動作設定ファイルは、相関算出用の相関器のタイプを指定するタイプcorrelator_type、及びパケット間相関算出用の重み係数w1, w2, w3を含む。ループ回数num_tは、時間オフセットtoffsetをループ変数としたループ処理を実行する回数を示す。ループ回数num_kは、伸縮率kをループ変数としたループ処理を実行する回数を示す。過去のサーチウィンドウの時間長search_window_past(sec)は、過去のどこまで遡ってサーチ対象の範囲にするかを示す。未来のサーチウィンドウの時間長search_window_future(sec)は、未来のどこまでをサーチ対象の範囲にするかを示す。重み係数w1はタイムスタンプ用、重み係数w2はパケット長用、重み係数w3はRSSI用である。ステップS18では、サーバ14のCPUが、動作設定ファイルの情報を配列SetDに格納する。ステップS19では、サーバ14のCPUが、配列SetDの内容を動作設定ファイルのオフセットRSSI_Offset、ガード時間Guard_Time等のローカル変数に代入し、処理は図5に示す電波同定処理に戻ってステップS2の処理へ進む。 In step S17, the CPU of the server 14 opens the operation setting file. As shown in the lower right part of FIG. 6 for convenience, the operation setting file includes the RSSI offset RSSI_Offset (dB) for compensating for the difference in RSSI, the scale of the time axis is expanded and contracted, and the offset on the time axis is added. Guard time for Guard_Time (sec) etc. is included. In this example, the behavior configuration file contains the minimum time offset t offset t_begin (sec), the maximum time offset t offset t_end (sec), and the number of loops num_t (times). In addition, the operation setting file contains the minimum value k_begin of the time expansion / contraction rate (hereinafter, also simply referred to as "expansion / contraction rate") k, the maximum value k_end of the expansion / contraction rate k, the number of loops num_k (times), and the RSSI offset RSSI_Offset (dB). Includes guard time Guard_Time (sec), etc. The operation setting file further includes the time length search_window_past (sec) of the past search window for calculating the inter-packet correlation and the time length search_window_future (sec) of the future search window for calculating the inter-packet correlation. The operation setting file contains the type correlator_type that specifies the type of correlator for correlation calculation, and the weighting coefficients w1, w2, and w3 for interpacket correlation calculation. The number of loops num_t indicates the number of times to execute loop processing with the time offset t offset as a loop variable. The number of loops num_k indicates the number of times that loop processing is executed with the expansion / contraction rate k as a loop variable. The time length of the past search window search_window_past (sec) indicates how far back in the past the search target range is set. The time length of the future search window search_window_future (sec) indicates how far in the future the search target range should be. The weighting factor w1 is for the time stamp, the weighting factor w2 is for the packet length, and the weighting factor w3 is for RSSI. In step S18, the CPU of the server 14 stores the information of the operation setting file in the array SetD. In step S19, the CPU of the server 14 assigns the contents of the array SetD to local variables such as offset RSSI_Offset and guard time Guard_Time in the operation setting file, and the process returns to the radio wave identification process shown in FIG. 5 and proceeds to the process of step S2. move on.

このように、サーバ14のCPUは、ステップS1のログデータ読込処理を実行することで、サーバ14のメモリから受信モジュール13−1,13−2のログデータを読み込む読込手段として機能する。 In this way, the CPU of the server 14 functions as a reading means for reading the log data of the receiving modules 13-1 and 13-2 from the memory of the server 14 by executing the log data reading process of step S1.

図5において、ステップS2では、サーバ14のCPUが、受信モジュール13−1,13−2のログデータをガード時間Guard_Timeを含めてクリッピングしてタイムスタンプの頭出しをするログデータクリッピング処理を行う。ガード時間Guard_Timeは、後述するタイムスタンプの時間軸のスケールの伸縮と、時間軸上のオフセットを行う際のマージンとして用いられる。図7は、ログデータクリッピング処理の一例を説明するフローチャートである。図7において、ステップS21では、サーバ14のCPUが、配列RS1のデータと配列RS2のデータのクリッピングタイミングを抽出する。具体的には、例えば受信モジュール13−2に届いた先頭データと末尾データのタイムスタンプが、受信モジュール13−1に届いた先頭データと末尾データのタイムスタンプより常にガード時間Guard_Time以上広くなるようなクリッピングタイミングを抽出する。ステップS22では、サーバ14のCPUが、配列RS1のデータを対応するクリッピングタイミングでクリッピングして新たな配列NewRS1に格納する。ステップS23では、サーバ14のCPUが、配列RS2のデータを対応するクリッピングタイミングでクリッピングして新たな配列NewRS2に格納する。ステップS24では、サーバ14のCPUが、配列NewRS1のパケット数を変数NumPktNewRS1に格納する。ステップS25では、サーバ14のCPUが、配列NewRS2のパケット数を変数NumPktNewRS2に格納する。ステップS26では、サーバ14のCPUが、配列NewRS2のパケット開始時刻を変数NewRS2_Startに格納する。ステップS27では、サーバ14のCPUが、配列NewRS2のパケット終了時刻を変数NewRS2_Stopに格納する。ステップS28では、サーバ14のCPUが、変数NewRS2_Startと変数NewRS2_Stopをマージする。 In FIG. 5, in step S2, the CPU of the server 14 performs a log data clipping process of clipping the log data of the receiving modules 13-1 and 13-2 including the guard time Guard_Time to find the time stamp. The guard time Guard_Time is used as a margin for expanding and contracting the scale of the time axis of the time stamp, which will be described later, and offsetting on the time axis. FIG. 7 is a flowchart illustrating an example of log data clipping processing. In FIG. 7, in step S21, the CPU of the server 14 extracts the clipping timing of the data of the array RS1 and the data of the array RS2. Specifically, for example, the time stamps of the first data and the last data that arrived at the receiving module 13-2 are always wider than the time stamps of the first data and the last data that arrived at the receiving module 13-1 by the guard time Guard_Time or more. Extract the clipping timing. In step S22, the CPU of the server 14 clips the data of the array RS1 at the corresponding clipping timing and stores it in the new array NewRS1. In step S23, the CPU of the server 14 clips the data of the array RS2 at the corresponding clipping timing and stores it in the new array NewRS2. In step S24, the CPU of the server 14 stores the number of packets in the array NewRS1 in the variable NumPktNewRS1. In step S25, the CPU of the server 14 stores the number of packets in the array NewRS2 in the variable NumPktNewRS2. In step S26, the CPU of the server 14 stores the packet start time of the array NewRS2 in the variable NewRS2_Start. In step S27, the CPU of the server 14 stores the packet end time of the array NewRS2 in the variable NewRS2_Stop. In step S28, the CPU of the server 14 merges the variable NewRS2_Start and the variable NewRS2_Stop.

図8は、ログデータクリッピング処理を4つのケースについて説明する図である。図8では、受信信号系列のデータストリーム長をハッチングで示す。ログデータクリッピング処理は、例えば受信モジュール13−2に届いた先頭データと末尾データのタイムスタンプが、受信モジュール13−1に届いた先頭データと末尾データのタイムスタンプより常にガード時間Guard_Time(sec)以上広くなるようにデータを切り取る。図8では便宜上、配列RS1のデータのタイムスタンプをTSRS1で示し、配列RS2のデータのタイムスタンプをTSRS2で示す。 FIG. 8 is a diagram illustrating four cases of log data clipping processing. In FIG. 8, the data stream length of the received signal sequence is shown by hatching. In the log data clipping process, for example, the time stamps of the first data and the last data that arrived at the receiving module 13-2 are always guard time Guard_Time (sec) or more than the time stamps of the first data and the last data that arrived at the receiving module 13-1. Cut the data so that it is wide. In FIG. 8, for convenience, the time stamp of the data of the sequence RS1 is shown by TS RS1 , and the time stamp of the data of the sequence RS2 is shown by TS RS2 .

図8中、(a)に示すケースIでは、配列RS1の先頭データのタイムスタンプが、配列RS2の先頭データのタイムスタンプより遅く、配列RS1の末尾データのタイムスタンプが、配列RS2の末尾データのタイムスタンプより早い。このケースIでは、クリッピングされた配列RS1のタイムスタンプは、クリッピング前のタイムスタンプと同じである。一方、クリッピングされた配列RS2のタイムスタンプは、クリッピング前と比べると、先頭データのタイムスタンプが遅く、末尾データのタイムスタンプが早い。また、クリッピングされた配列RS2のタイムスタンプは、クリッピングされた配列RS1のタイムスタンプと比べると、先頭データがガード時間Guard_Timeだけ早く、末尾データがガード時間Guard_Timeだけ遅い。 In case I shown in FIG. 8A, the time stamp of the start data of the array RS1 is later than the time stamp of the start data of the array RS2, and the time stamp of the end data of the array RS1 is the time stamp of the end data of the array RS2. Faster than the time stamp. In this case I, the time stamp of the clipped array RS1 is the same as the time stamp before clipping. On the other hand, the time stamp of the clipped array RS2 has a slower time stamp of the first data and a faster time stamp of the last data than before clipping. In addition, the time stamp of the clipped array RS2 is earlier than the time stamp of the clipped array RS1 by the guard time Guard_Time for the first data and later by the guard time Guard_Time for the last data.

図8中、(b)に示すケースIIでは、配列RS1先頭データと末尾データのタイムスタンプが、配列RS2の先頭データと末尾データのタイムスタンプより夫々早い。このケースIIでは、クリッピングされた配列RS1のタイムスタンプは、クリッピング前と比べると、先頭データのタイムスタンプが遅い。一方、クリッピングされた配列RS2のタイムスタンプは、クリッピング前と比べると、末尾データのタイムスタンプが早い。また、クリッピングされた配列RS2のタイムスタンプは、クリッピングされた配列RS1のタイムスタンプと比べると、先頭データがガード時間Guard_Timeだけ早く、末尾データがガード時間Guard_Timeだけ遅い。 In case II shown in FIG. 8 (b), the time stamps of the start data and the end data of the array RS1 are earlier than the time stamps of the start data and the end data of the array RS2, respectively. In this case II, the time stamp of the clipped array RS1 has a slower time stamp of the first data than that before clipping. On the other hand, the time stamp of the clipped array RS2 has a faster time stamp of the last data than before clipping. In addition, the time stamp of the clipped array RS2 is earlier than the time stamp of the clipped array RS1 by the guard time Guard_Time for the first data and later by the guard time Guard_Time for the last data.

図8中、(c)に示すケースIIIでは、配列RS1先頭データと末尾データのタイムスタンプが、配列RS2の先頭データと末尾データのタイムスタンプより夫々遅い。このケースIIIでは、クリッピングされた配列RS1のタイムスタンプは、クリッピング前と比べると、先頭データのタイムスタンプが遅く、末尾データのタイムスタンプが早い。一方、クリッピングされた配列RS2のタイムスタンプは、クリッピング前と比べると、先頭データのタイムスタンプが遅い。また、クリッピングされた配列RS2のタイムスタンプは、クリッピングされた配列RS1のタイムスタンプと比べると、先頭データがガード時間Guard_Time以上早く、末尾データがガード時間Guard_Timeだけ遅い。 In case III shown in FIG. 8C, the time stamps of the start data and the end data of the array RS1 are later than the time stamps of the start data and the end data of the array RS2, respectively. In this case III, the time stamp of the clipped array RS1 has a later time stamp of the first data and a faster time stamp of the last data than before clipping. On the other hand, the time stamp of the clipped array RS2 has a slower time stamp of the first data than before clipping. In addition, the time stamp of the clipped array RS2 is earlier than the time stamp of the clipped array RS1 by the guard time Guard_Time or faster for the first data and later than the guard time Guard_Time for the last data.

図8中、(d)に示すケースIVでは、配列RS1の先頭データのタイムスタンプが、配列RS2の先頭データのタイムスタンプより早く、配列RS1の末尾データのタイムスタンプが、配列RS2の末尾データのタイムスタンプより遅い。このケースIVでは、クリッピングされた配列RS1のタイムスタンプは、クリッピング前と比べると、先頭データのタイムスタンプが遅く、末尾データのタイムスタンプが早い。一方、クリッピングされた配列RS2のタイムスタンプは、クリッピング前のタイムスタンプと同じである。また、クリッピングされた配列RS2のタイムスタンプは、クリッピングされた配列RS1のタイムスタンプと比べると、先頭データがガード時間Guard_Timeだけ早く、末尾データがガード時間Guard_Timeだけ遅い。 In case IV shown in FIG. 8 (d), the time stamp of the start data of the array RS1 is earlier than the time stamp of the start data of the array RS2, and the time stamp of the end data of the array RS1 is the time stamp of the end data of the array RS2. Slower than the time stamp. In this case IV, the time stamp of the clipped array RS1 has a later time stamp of the first data and a faster time stamp of the last data than before clipping. On the other hand, the time stamp of the clipped array RS2 is the same as the time stamp before clipping. In addition, the time stamp of the clipped array RS2 is earlier than the time stamp of the clipped array RS1 by the guard time Guard_Time for the first data and later by the guard time Guard_Time for the last data.

ステップS28の後、処理は図5に示す電波同定処理に戻る。具体的には、図5に示すステップS3,S4に対して、時間オフセットtoffsetをループ変数としたループ処理をループ回数num_tだけ実行する処理を、伸縮率kをループ変数としたループ処理をループ回数num_kだけ実行する処理へ進む。伸縮率kは、タイムスタンプの時間軸のスケールを伸縮する比率を表す。時間オフセットtoffsetは、タイムスタンプの絶対値に時間軸上のオフセット量を加減算する秒数を表す。時間オフセットtoffsetをループ変数としたループ処理は、ループステップパラメータt_step=(t_end-t_begin)/num_t)をt_begin〜t_endの範囲で変化させながらループ回数num_tだけ実行する。また、伸縮率kをループ変数としたループ処理は、ループステップパラメータk_step=(k_end-k_begin)/num_t)をk_begin〜k_endの範囲で変化させながらループ回数num_kだけ実行する。 After step S28, the process returns to the radio wave identification process shown in FIG. Specifically, for steps S3 and S4 shown in FIG. 5, a loop process in which the time offset t offset is set as a loop variable is executed for the number of loops num_t, and a loop process in which the expansion / contraction rate k is set as a loop variable is looped. Proceed to the process to execute the number of times num_k. The expansion / contraction ratio k represents the ratio of expansion / contraction of the scale on the time axis of the time stamp. The time offset t offset represents the number of seconds for adding or subtracting the offset amount on the time axis to the absolute value of the time stamp. Loop processing with the time offset t offset as the loop variable is executed for the number of loops num_t while changing the loop step parameter t_step = (t_end-t_begin) / num_t) in the range of t_begin to t_end. Further, the loop processing using the expansion / contraction rate k as the loop variable is executed by the number of loops num_k while changing the loop step parameter k_step = (k_end-k_begin) / num_t) in the range of k_begin to k_end.

このように、サーバ14のCPUは、ログデータクリッピング処理を実行することで、受信モジュール13−1,13−2のログデータをガード時間Guard_Timeを含めてクリッピングしてタイムスタンプの頭出しをする抽出手段として機能する。受信モジュール13−2が受信した第2の受信信号系列の先頭データと末尾データの時刻情報は、受信モジュール13−1が受信した第1の受信信号系列の先頭データと末尾データの時刻情報より常にガード時間Guard_Time以上広くなる。 In this way, the CPU of the server 14 performs the log data clipping process to clip the log data of the receiving modules 13-1 and 13-2 including the guard time Guard_Time to cue the time stamp. Functions as a means. The time information of the first data and the last data of the second received signal sequence received by the receiving module 13-2 is always more than the time information of the first data and the last data of the first received signal sequence received by the receiving module 13-1. Guard time is wider than Guard_Time.

図5において、ステップS3では、サーバ14のCPUが、ログデータの相関を算出するログデータ相関算出処理を行う。図9及び図10は、ログデータ相関算出処理の一例を説明するフローチャートである。図9において、ステップS301では、サーバ14のCPUが、ループステップパラメータk_step, t_stepを算出する。ステップS302では、サーバ14のCPUが、ループカウンタnum_loopをリセットする。ステップS303では、サーバ14のCPUが、k,toffsetループ中で相関値corr(k,toffset)の最大値を保持する最大相関値保持変数sigma_correlation_maxをリセットする。k,toffsetループは、時間オフセットtoffsetをループ変数としたループ処理をループ回数num_tだけ実行するループ処理を、伸縮率kをループ変数としたループ処理をループ回数num_kだけ実行する処理を含む。相関値corr(k,toffset)は、ある[k,toffset]の組み合わせにおける配列と配列(例えば、配列RS1と配列RS2)の相関値の和である。ステップS304では、サーバ14のCPUが、最大値保持変数k_at_max, t_at_maxをリセットする。最大値保持変数k_at_maxは、相関値corr(k,toffset)が最大値を取るときの伸縮率kの値を保持する変数である。最大値保持変数t_at_maxは、相関値corr(k,toffset)が最大値を取るときの時間オフセットtoffsetの値を保持する変数である。 In FIG. 5, in step S3, the CPU of the server 14 performs a log data correlation calculation process for calculating the log data correlation. 9 and 10 are flowcharts illustrating an example of log data correlation calculation processing. In FIG. 9, in step S301, the CPU of the server 14 calculates the loop step parameters k_step and t_step. In step S302, the CPU of the server 14 resets the loop counter num_loop. In step S303, the CPU of the server 14 resets the maximum correlation value holding variable sigma_correlation_max that holds the maximum value of the correlation value corr (k, t offset ) in the k, t offset loop. The k, t offset loop includes a loop process that executes a loop process with the time offset t offset as a loop variable for the number of loops num_t, and a process for executing a loop process with the expansion / contraction rate k as a loop variable for the number of loops num_k. The correlation value corr (k, t offset ) is the sum of the correlation values of the array and the array (for example, array RS1 and array RS2) in a certain [k, t offset ] combination. In step S304, the CPU of the server 14 resets the maximum value holding variables k_at_max and t_at_max. The maximum value holding variable k_at_max is a variable that holds the value of the expansion / contraction rate k when the correlation value corr (k, t offset ) takes the maximum value. The maximum value holding variable t_at_max is a variable that holds the value of the time offset t offset when the correlation value corr (k, t offset ) takes the maximum value.

次に、ステップS305〜S314の処理、ステップS315〜S322のループ変数配列RS12のループ処理をNumRS12分だけ実行する処理、及びステップS323〜S326の処理を、伸縮率kと時間オフセットtoffsetのk,toffsetループ分だけ実行する処理へ進む。 Next, the processing of steps S305 to S314, the processing of executing the loop processing of the loop variable array RS12 of steps S315 to S322 for NumRS12 minutes, and the processing of steps S323 to S326 are performed with the expansion / contraction ratio k and the time offset t offset k, The process proceeds to execute only the t offset loop.

図9において、ステップS305では、サーバ14のCPUが、k,toffsetループ内の伸縮率kを算出する。ステップS306では、サーバ14のCPUが、k,toffsetループ内の時間オフセットtoffsetを算出する。ステップS307では、サーバ14のCPUが、ループ変数配列RS12を宣言して初期化する。ステップS308では、サーバ14のCPUが、新たな配列NewRS1に伸縮率kと時間オフセットtoffsetを適用する。ステップS309では、サーバ14のCPUが、配列NewRS1のパケット開始時刻を変数NewRS1_Startに格納し、配列NewRS1のパケット終了時刻を変数NewRS1_Stopに格納する。ステップS310では、サーバ14のCPUが、変数NewRS1_Start, NewRS1_Stopをクリアする。ステップS311では、サーバ14のCPUが、変数NewRS1と変数NewRS2をマージしてループ変数配列RS12に格納する。ステップS312では、サーバ14のCPUが、ループ変数配列RS12をタイムスタンプでソートする。ステップS313では、サーバ14のCPUが、ループ変数配列RS12のパケット数を変数NumRS12に格納する。ステップS314では、サーバ14のCPUが、ループ変数配列RS12のループのループ内変数をリセットし、処理は図10に示すステップS315へ進む。 In FIG. 9, in step S305, the CPU of the server 14 calculates the expansion / contraction rate k in the k, t offset loop. In step S306, the CPU of the server 14 calculates the time offset t offset in the k, t offset loop. In step S307, the CPU of the server 14 declares and initializes the loop variable array RS12. In step S308, the CPU of the server 14 applies the stretch ratio k and the time offset t offset to the new array NewRS1. In step S309, the CPU of the server 14 stores the packet start time of the array NewRS1 in the variable NewRS1_Start, and stores the packet end time of the array NewRS1 in the variable NewRS1_Stop. In step S310, the CPU of the server 14 clears the variables NewRS1_Start and NewRS1_Stop. In step S311 the CPU of the server 14 merges the variable NewRS1 and the variable NewRS2 and stores them in the loop variable array RS12. In step S312, the CPU of the server 14 sorts the loop variable array RS12 by the time stamp. In step S313, the CPU of the server 14 stores the number of packets in the loop variable array RS12 in the variable NumRS12. In step S314, the CPU of the server 14 resets the in-loop variables of the loop of the loop variable array RS12, and the process proceeds to step S315 shown in FIG.

図10において、ステップS315では、サーバ14のCPUが、ループ変数配列RS12の値をループ変数配列RS12のループのループ内変数に格納する。ステップS316では、サーバ14のCPUが、ループ変数配列RS12の各パケットのパケット長を変数time_dutyに格納する。ステップS317では、サーバ14のCPUが、最大相関値保持変数sigma_correlation_maxをリセットする。ステップS318では、サーバ14のCPUが、ループ変数配列RS12の各パケットに配列RS1又は配列RS2のフラグを立て、ループ変数配列RS12の各パケットに、各受信モジュール13−1,13−2毎のインデックスを付与する。フラグは、変数NewRS1と変数NewRS2をマージしたループ変数配列RS12の各パケットが、配列RS1又は配列RS2であることを示す。ステップS319では、サーバ14のCPUが、ループ変数配列RS12の結果を変数RS12Recordに記録する。ステップS320では、サーバ14のCPUが、イベントがあるタイムスタンプの区切り毎に、配列NewRS1と配列NewRS2のRSSIの差分に基づく相関値を算出する。ここで、イベントとは、2つの受信信号系列のデータの遷移、即ち、立ち上がり又は立ち下がりのことを言う。ステップS321では、サーバ14のCPUが、前段で算出した相関値を積算して変数SigmaCorrに格納する。ステップS322では、サーバ14のCPUが、変数SigmaCorrを変数RS12CorRecordに記録する。ステップS323では、サーバ14のCPUが、ループ変数配列RS12をクリアする。ステップS324では、サーバ14のCPUが、変数SigmaCorrの最大値が閾値を超えていると変数SigmaCorrMaxに記録する。最大値が閾値以下である変数SigmaCorrを除外することで、最大値であっても低い相関値を排除できるので、例えばノイズ成分等の相関性が判定されてしまうことを防止可能である。ステップS325では、サーバ14のCPUが、変数SigmaCorrMaxとなるときの最大相関値保持変数k_at_max, t_at_maxを配列CorrRecordに格納する。ステップS326では、サーバ14のCPUが、k,toffsetループの回数を表すループ変数NumLoopをインクリメントする。 In FIG. 10, in step S315, the CPU of the server 14 stores the value of the loop variable array RS12 in the loop in-loop variable of the loop variable array RS12. In step S316, the CPU of the server 14 stores the packet length of each packet of the loop variable array RS12 in the variable time_duty. In step S317, the CPU of the server 14 resets the maximum correlation value holding variable sigma_correlation_max. In step S318, the CPU of the server 14 sets the flag of the array RS1 or the array RS2 for each packet of the loop variable array RS12, and indexes each of the receiving modules 13-1 and 13-2 for each packet of the loop variable array RS12. Is given. The flag indicates that each packet in the loop variable array RS12, which merges the variable NewRS1 and the variable NewRS2, is array RS1 or array RS2. In step S319, the CPU of the server 14 records the result of the loop variable array RS12 in the variable RS12Record. In step S320, the CPU of the server 14 calculates a correlation value based on the difference between the RSSIs of the array NewRS1 and the array NewRS2 for each time stamp delimiter with an event. Here, the event means a transition of data of two received signal sequences, that is, a rising edge or a falling edge. In step S321, the CPU of the server 14 integrates the correlation values calculated in the previous stage and stores them in the variable SigmaCorr. In step S322, the CPU of the server 14 records the variable SigmaCorr in the variable RS12CorRecord. In step S323, the CPU of the server 14 clears the loop variable array RS12. In step S324, the CPU of the server 14 records in the variable SigmaCorrMax that the maximum value of the variable SigmaCorr exceeds the threshold value. By excluding the variable SigmaCorr whose maximum value is equal to or less than the threshold value, it is possible to exclude a low correlation value even if it is the maximum value, so that it is possible to prevent the correlation of noise components and the like from being determined, for example. In step S325, the CPU of the server 14 stores the maximum correlation value holding variables k_at_max and t_at_max when the variables are SigmaCorrMax in the array CorrRecord. In step S326, the CPU of the server 14 increments the loop variable NumLoop, which represents the number of k, t offset loops.

これにより、図9に示すステップS305〜S314の処理、図10に示すステップS315〜S322の処理をループ変数配列RS12のループ分だけ実行する処理、及び図10に示すステップS323〜S326の処理を、伸縮率kと時間オフセットtoffsetのk,toffsetループ分だけ実行した後、処理は図5に示す電波同定処理に戻ってステップS4の処理へ進む。 As a result, the processing of steps S305 to S314 shown in FIG. 9, the processing of executing the processing of steps S315 to S322 shown in FIG. 10 for the loop of the loop variable array RS12, and the processing of steps S323 to S326 shown in FIG. 10 are performed. After executing the k and t offset loops of the expansion / contraction rate k and the time offset t offset , the process returns to the radio wave identification process shown in FIG. 5 and proceeds to the process of step S4.

図5において、ステップS4では、サーバ14のCPUが、相関値が閾値を超える最大値になる伸縮率と時間オフセットを決定する、伸縮率とオフセット決定処理を行う。図11は、伸縮率とオフセット決定処理の一例を説明するフローチャートである。図11において、ステップS41では、サーバ14のCPUが、配列CorrRecordに格納されている最大相関値保持変数k_at_max, t_at_maxを読み出す。ステップS42では、サーバ14のCPUが、k,toffsetループ終了時の最大相関値保持変数k_at_maxを最適な伸縮率kと決定する。また、ステップS42では、サーバ14のCPUが、k,toffsetループ終了時の最大相関値保持変数t_at_maxを最適な時間オフセットtoffsetと決定する。ステップS42の後、処理は図5に示す電波同定処理に戻る。 In FIG. 5, in step S4, the CPU of the server 14 performs the expansion / contraction rate and the offset determination process for determining the expansion / contraction rate and the time offset at which the correlation value becomes the maximum value exceeding the threshold value. FIG. 11 is a flowchart illustrating an example of expansion / contraction ratio and offset determination processing. In FIG. 11, in step S41, the CPU of the server 14 reads out the maximum correlation value holding variables k_at_max and t_at_max stored in the array CorrRecord. In step S42, the CPU of the server 14 determines the maximum correlation value holding variable k_at_max at the end of the k, t offset loop as the optimum expansion / contraction rate k. Further, in step S42, the CPU of the server 14 determines the maximum correlation value holding variable t_at_max at the end of the k, t offset loop as the optimum time offset t offset . After step S42, the process returns to the radio wave identification process shown in FIG.

これにより、図9乃至図11に示すステップS3,S4の処理に対して、時間オフセットtoffsetをループ変数としたループ処理をループ回数num_tだけ実行する処理を、伸縮率kをループ変数としたループ処理をループ回数num_kだけ実行する。その後、処理は図5に示すステップS5の処理へ進む。 As a result, with respect to the processes of steps S3 and S4 shown in FIGS. 9 to 11, the process of executing the loop process with the time offset t offset as the loop variable for the number of loops num_t is performed with the expansion / contraction rate k as the loop variable. Execute the process for the number of loops num_k. After that, the process proceeds to the process of step S5 shown in FIG.

このように、サーバ14のCPUは、ステップS3,S4の処理を実行することで、相関値が最大値になる時間伸縮率k及び前記時間オフセットtoffsetを決定する決定手段として機能する。決定手段は、受信モジュール13−1,13−2が受信した第1及び第2の受信信号系列うち、第2の受信信号系列の時間軸のスケールを時間伸縮率kで伸縮しながら第1の受信信号系列との相関値を算出する。また、決定手段は、第2の受信信号系列を時間軸上で時間オフセットtoffsetだけオフセットして相関値が最大値になる時間伸縮率k及び前記時間オフセットtoffsetを決定する。更に、決定手段は、上記抽出手段が抽出した時系列データから相関値が最大値になる時間伸縮率k及び前記時間オフセットtoffsetを決定する。 In this way, the CPU of the server 14 functions as a determination means for determining the time expansion / contraction rate k at which the correlation value becomes the maximum value and the time offset t offset by executing the processes of steps S3 and S4. The determining means is the first of the first and second received signal sequences received by the receiving modules 13-1 and 13-2, while expanding and contracting the scale of the time axis of the second received signal sequence with the time expansion / contraction rate k. Calculate the correlation value with the received signal series. Further, the determining means offsets the second received signal sequence by the time offset t offset on the time axis to determine the time expansion / contraction rate k at which the correlation value becomes the maximum value and the time offset t offset . Further, the determination means determines the time expansion / contraction rate k at which the correlation value becomes the maximum value and the time offset t offset from the time series data extracted by the extraction means.

従って、この例では、各受信モジュール13−1,13−2におけるタイムスタンプの個体差を次のように補償できる。先ず、サーバ14が、受信モジュール13−1,13−2から受信した受信信号系列の特徴が類似するように、一方の受信モジュールからの受信信号系列の時間軸のスケールを伸縮しながら相関を算出する。受信信号系列の特徴とは、受信信号系列に含まれるパケットの送信開始時刻、パケット長、RSSI等である。同時に、当該一方の受信モジュールからの受信信号系列を時間軸上でオフセットしながら相関を算出する。受信信号系列の時間軸上のオフセットとは、時間軸のスケールを伸縮した受信信号系列を当該時間軸上で時間オフセットだけ少なくとも一方向(正方向又は負方向)へスライドさせることを意味する。これにより、一方の受信モジュールからの受信信号系列の特徴が他方の受信モジュールからの受信信号系列の特徴に最も類似する時間軸の伸縮率kと時間オフセットtoffsetを決定できる。つまり、受信モジュール13−1,13−2から受信した2つの受信信号系列の相関値が閾値を超える最大値になる時間軸の伸縮率kと時間オフセットtoffsetを決定し、これら2つの受信信号系列を同期させることができる。 Therefore, in this example, the individual difference of the time stamp in each of the receiving modules 13-1 and 13-2 can be compensated as follows. First, the server 14 calculates the correlation while expanding and contracting the scale of the time axis of the received signal series from one of the receiving modules so that the characteristics of the received signal series received from the receiving modules 13-1 and 13-2 are similar. To do. The characteristics of the received signal sequence are the transmission start time, packet length, RSSI, and the like of the packet included in the received signal sequence. At the same time, the correlation is calculated while offsetting the received signal sequence from the one receiving module on the time axis. The offset on the time axis of the received signal sequence means that the received signal sequence with the scale of the time axis expanded or contracted is slid on the time axis by the time offset in at least one direction (positive direction or negative direction). Thereby, the expansion / contraction rate k and the time offset t offset on the time axis whose characteristics of the received signal sequence from one receiving module are most similar to the characteristics of the received signal sequence from the other receiving module can be determined. That is, the expansion / contraction rate k and the time offset t offset of the time axis at which the correlation value of the two received signal sequences received from the receiving modules 13-1 and 13-2 becomes the maximum value exceeding the threshold value are determined, and these two received signals are received. The series can be synchronized.

図12は、タイムスタンプの個体差の補償を説明する図である。図12に示す処理は、サーバ14のCPUが実行する図5に示すステップS3,S4に相当する手順ST1〜ST3を含む。手順ST1では、受信モジュール13−1と受信モジュール13−2夫々の受信信号系列の相関値を、例えば時間ts毎に算出する。例えば、時間ts=30(sec)である。手順ST2では、相関値が0でないと、受信モジュール13−2のタイムスタンプの時間軸(幅)に伸縮率kを乗算し、更に時間オフセットtoffsetを加算して、そのときの相関値を再算出する。ここで、伸縮率kは0ではない正の実数であり、時間オフセットtoffsetは実数である。手順ST3では、相関値が閾値を超える最大値になる伸縮率kと時間オフセットtoffsetを求める。具体的には、手順ST3では、伸縮率kと時間オフセットtoffsetを変えながら前段の手順ST2の処理を繰り返す。 FIG. 12 is a diagram illustrating compensation for individual differences in time stamps. The process shown in FIG. 12 includes procedures ST1 to ST3 corresponding to steps S3 and S4 shown in FIG. 5 executed by the CPU of the server 14. In the procedure ST1, the correlation value of the received signal series of the receiving module 13-1 and the receiving module 13-2 is calculated, for example, every time t s . For example, time t s = 30 (sec). In step ST2, if the correlation value is not 0, the time axis (width) of the time stamp of the receiving module 13-2 is multiplied by the expansion / contraction rate k, the time offset t offset is added, and the correlation value at that time is re-established. calculate. Here, the expansion / contraction rate k is a positive real number that is not 0, and the time offset t offset is a real number. In the procedure ST3, the expansion / contraction rate k and the time offset t offset at which the correlation value becomes the maximum value exceeding the threshold value are obtained. Specifically, in the procedure ST3, the process of the previous procedure ST2 is repeated while changing the expansion / contraction ratio k and the time offset t offset .

この結果、サーバ14に届いた受信モジュール13−2の受信信号系列は、図12に下部に示すように、サーバ14に届いた受信モジュール13−1の受信信号系列との相関値が閾値を超える最大値になるように、タイムスタンプの個体差が補正される。具体的には、受信モジュール13−1の受信信号系列の時間軸のスケールが伸縮率kで伸縮され、時間軸上で時間オフセットtoffsetだけ少なくとも一方向(正方向又は負方向)へスライドされる。図12中、時間軸上の正方向は右方向であり、負方向は左方向である。なお、図12は模式図であるため、各受信信号系列の振幅は、RSSIを同じスケールで表すものではない。 As a result, the correlation value of the received signal sequence of the receiving module 13-2 that arrived at the server 14 with the received signal sequence of the receiving module 13-1 that arrived at the server 14 exceeds the threshold value, as shown in the lower part of FIG. The individual difference of the time stamp is corrected so that it becomes the maximum value. Specifically, the scale of the received signal series of the receiving module 13-1 is expanded and contracted by the expansion / contraction rate k, and is slid in at least one direction (positive direction or negative direction) by the time offset t offset on the time axis. .. In FIG. 12, the positive direction on the time axis is the right direction, and the negative direction is the left direction. Since FIG. 12 is a schematic diagram, the amplitude of each received signal sequence does not represent RSSI on the same scale.

このようにして、伸縮率kに従って複数の受信モジュールが受信した受信信号系列に付与するタイムスタンプの時間軸のスケールを伸縮しながら相関値を算出し、同時に、時間軸上の時間オフセットtoffsetに従って受信信号系列をオフセットしながら相関値を算出する。これにより、相関値が閾値を超える最大値になる時間軸のスケールの伸縮率kと時間軸上の時間オフセットtoffsetを決定することができる。 In this way, the correlation value is calculated while expanding and contracting the scale of the time axis of the time stamps given to the received signal series received by the plurality of receiving modules according to the expansion / contraction rate k, and at the same time, according to the time offset t offset on the time axis. The correlation value is calculated while offsetting the received signal sequence. Thereby, the expansion / contraction rate k of the scale on the time axis and the time offset t offset on the time axis can be determined so that the correlation value becomes the maximum value exceeding the threshold value.

パケットの紐付け処理は、上記の如く算出した相関値が閾値を超える最大値になる伸縮率kと時間オフセットtoffsetを受信モジュール13−2のタイムスタンプに適用することで、各受信モジュール13−1,13−2が受信した同じパケット同士を紐付ける。具体的には、受信モジュール13−1の受信信号系列及び受信モジュール13−2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及びRSSIがいずれも類似するパケット同士を1対1で紐付ける。また、紐付け候補のパケットの受信順が異なっていたり、複数のパケットに重複紐付けが発生していたりするパケットの紐付けの誤りがある場合、後述する揺さぶりを行い、パケットの順番は入れ替わらないというルールを適用してパケットの紐付けの誤りを訂正する。 In the packet linking process, the expansion / contraction rate k and the time offset t offset at which the correlation value calculated as described above becomes the maximum value exceeding the threshold value are applied to the time stamp of the receiving module 13-2, so that each receiving module 13- The same packets received by 1, 13-2 are linked to each other. Specifically, each packet of the received signal sequence of the receiving module 13-1 and the received signal sequence of the receiving module 13-2 is compared with each other, and the transmission start time, packet length, and RSSI of each packet are compared. Associate similar packets on a one-to-one basis. In addition, if the reception order of the packet candidates for linking is different, or if there is an error in the linking of packets, such as multiple packets being duplicated, the packet order will be changed by shaking as described later. Correct the error in packet association by applying the rule that there is no packet.

図5において、ステップS5では、サーバ14のCPUが、パケットを紐付けるパケット紐付け処理を行う。図13は、パケット紐付け処理の一例を説明するフローチャートである。図13に示すパケット紐付け処理は、ステップS51〜S53の処理と、ステップS54〜S59の処理を配列NewRS2のループ分実行する処理とを、配列NewRS1のループ分実行する。 In FIG. 5, in step S5, the CPU of the server 14 performs a packet linking process for linking packets. FIG. 13 is a flowchart illustrating an example of the packet linking process. In the packet linking process shown in FIG. 13, the processes of steps S51 to S53 and the processes of steps S54 to S59 are executed for the loops of the array NewRS2 for the loops of the array NewRS1.

図13において、ステップS51では、サーバ14のCPUが、配列NewRS1の値を配列NewRS1のループのループ内変数に格納する。ステップS52では、サーバ14のCPUが、伸縮率kと時間オフセットtoffsetの夫々の最適値を配列NewRS1に適用して変数NewRS1Fixに格納する。ステップS53では、サーバ14のCPUが、相関値の最大値を格納する変数MaxCorrをリセットする。 In FIG. 13, in step S51, the CPU of the server 14 stores the value of the array NewRS1 in the in-loop variable of the loop of the array NewRS1. In step S52, the CPU of the server 14 applies the optimum values of the expansion / contraction rate k and the time offset t offset to the array NewRS1 and stores them in the variable NewRS1Fix. In step S53, the CPU of the server 14 resets the variable MaxCorr that stores the maximum value of the correlation value.

ステップS54では、サーバ14のCPUが、配列NewRS2の値を配列NewRS2のループのループ内変数に格納する。ステップS55では、サーバ14のCPUが、変数NewRS1Fixと配列NewRS2のパケットの差分情報を使って相関値を算出して変数Corrに格納する。この際、時間軸のスケールを伸縮した一方の受信信号系列を当該時間軸上でオフセットし、2つの受信信号系列を紐付けた後、揺さぶりにより再度紐付けをする処理も行う。揺さぶりとは、受信信号系列を時間軸上で一定範囲だけ両方向(正方向と負方向)へ振って(ずらして)微調整することで、パケットの紐付けの誤りを訂正する処理である。ステップS57では、サーバ14のCPUが、変数Corrの最大値を変数MaxCorrに記録し、変数Corrが最大値になるときの配列NewRS2のインデックスを変数NewRS2Indexに記録する。ステップS58では、サーバ14のCPUが、変数Corrの値を変数MemCorrに格納する。ステップS59では、変数MaxCorrのときの変数NewRS1Index, NewRS2Index(即ち、紐付き結果)を配列Identifiedに格納する。配列Identifiedは、後述する変数NewRS1Index, NewRS2Index、送信局12のソースアドレス等を含む紐付け結果を示す。 In step S54, the CPU of the server 14 stores the value of the array NewRS2 in the variable in the loop of the array NewRS2. In step S55, the CPU of the server 14 calculates the correlation value using the difference information of the packet of the variable NewRS1Fix and the packet of the array NewRS2 and stores it in the variable Corr. At this time, one of the received signal sequences whose scale on the time axis is expanded or contracted is offset on the time axis, the two received signal sequences are linked, and then the process of linking again by shaking is also performed. Shaking is a process of correcting an error in packet association by swinging (shifting) the received signal sequence in both directions (positive direction and negative direction) by a certain range on the time axis to make fine adjustments. In step S57, the CPU of the server 14 records the maximum value of the variable Corr in the variable MaxCorr, and records the index of the array NewRS2 when the variable Corr becomes the maximum value in the variable NewRS2Index. In step S58, the CPU of the server 14 stores the value of the variable Corr in the variable MemCorr. In step S59, the variables NewRS1Index and NewRS2Index (that is, the associated result) when the variable MaxCorr is stored are stored in the array Identified. The array Identified shows the association result including the variables NewRS1Index, NewRS2Index, the source address of the transmitting station 12, and the like, which will be described later.

このように、サーバ14のCPUは、ステップS5のパケット紐付け処理を実行することで、パケットを紐付ける紐付け手段として機能する。紐付け手段は、相関値が最大値になる時間伸縮率kと時間オフセットtoffsetを第2の受信信号系列の時刻情報に適用することで、第1の受信信号系列及び第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較する。また、紐付け手段は、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付ける。紐付け手段は、紐付け候補のパケットの受信順が異なっている場合、又は、複数のパケットに重複紐付けが発生している場合に、揺さぶりを行い、パケットの順番は入れ替わらないというルールを適用してパケットの紐付けの誤りを訂正する。揺さぶりは、紐付け候補のパケットを時間軸上で一定範囲だけ振って微調整することで、パケットの紐付けの誤りを訂正する。 In this way, the CPU of the server 14 functions as a linking means for linking packets by executing the packet linking process in step S5. The associating means applies the time expansion / contraction rate k at which the correlation value becomes the maximum value and the time offset t offset to the time information of the second received signal sequence, so that the first received signal sequence and the second received signal sequence Compare each packet of each of them. Further, the linking means links packets having similar transmission start times, packet lengths, and signal strengths of each packet on a one-to-one basis. The associating means shakes when the receiving order of the packets of the associating candidate is different, or when multiple packets are duplicated, and the order of the packets is not changed. Apply to correct packet association errors. Shaking corrects an error in packet linking by shaking the packet of the linking candidate by a certain range on the time axis and making fine adjustments.

図14は、正常なパケットの紐付け結果の一例を説明する図である。図14は、タイムスタンプの個体差を補償した後の受信モジュール13−1と受信モジュール13−2夫々の受信信号系列を示す。各受信信号系列中、パケットは四角形で表し、図14は、受信モジュール13−1と受信モジュール13−2の夫々が持つパケットを送信開始時刻、パケット長、及びRSSIにより比較し、比較した全ての要素が類似するパケット同士の紐付けに成功した例を示す。正しく紐付けられたパケットは、破線の矢印で示す。 FIG. 14 is a diagram illustrating an example of a normal packet association result. FIG. 14 shows the received signal sequences of the receiving module 13-1 and the receiving module 13-2 after compensating for the individual difference of the time stamp. In each received signal sequence, the packets are represented by rectangles, and FIG. 14 compares the packets of the receiving module 13-1 and the receiving module 13-2 by the transmission start time, the packet length, and RSSI, and all the comparisons are made. An example of successfully linking packets with similar elements is shown. Correctly linked packets are indicated by dashed arrows.

図15は、誤った結果を含むパケットの紐付け結果の一例を説明する図である。図15は、受信モジュール13−1の配置、又は、受信感度等の都合により、受信モジュール13−1においてパケットの受信欠損が発生した場合に、パケットの紐付けの誤りが生じた受信信号系列を示す。受信モジュール13−1で幾つかのパケットが破線の四角形で示すように欠損しているため、比較した全ての要素が類似するパケットが変化してしまい、その結果、パケットの紐付けの誤りが生じた例を示す。正しく紐付けられたパケットは、破線の矢印で示す。また、誤って紐付けられたパケットは、一点鎖線の矢印で示す。 FIG. 15 is a diagram illustrating an example of a packet association result including an erroneous result. FIG. 15 shows a received signal sequence in which an error in packet association occurs when a packet reception defect occurs in the receiving module 13-1 due to the arrangement of the receiving module 13-1 or the reception sensitivity. Shown. Since some packets are missing in the receiving module 13-1 as shown by the broken line rectangle, the packets having similar elements to all the compared elements are changed, and as a result, the packet association error occurs. An example is shown. Correctly linked packets are indicated by dashed arrows. In addition, packets that are erroneously linked are indicated by arrows on the alternate long and short dash line.

図16は、パケットの紐付け誤りの訂正を説明する図である。図16では、図15で示すパケットの紐付けの誤りを訂正する例を示す。図16の(a)に示すように、紐付け候補のパケットの受信順が異なっている場合、又は、複数のパケットに重複紐付けが発生している場合に、パケットの順番は入れ替わらないというルールを適用してパケットの紐付けの誤りを訂正する。具体的には、図16の(b)に示すように、受信モジュール13−2の紐付け候補のパケットの位置を時間軸上で揺さぶる。そして、図16の(c)に示すように、受信順の入れ替えや、重複紐付けが発生しない条件下で最も相関が高い、即ち、相関値が閾値を超える最大値になるパケット同士を紐付ける。 FIG. 16 is a diagram illustrating correction of packet association error. FIG. 16 shows an example of correcting an error in associating the packets shown in FIG. As shown in FIG. 16A, when the reception order of the packets of the linkage candidate is different, or when duplicate association occurs in a plurality of packets, the order of the packets is not changed. Apply the rule to correct the packet association error. Specifically, as shown in FIG. 16B, the position of the packet of the association candidate of the receiving module 13-2 is shaken on the time axis. Then, as shown in FIG. 16 (c), packets having the highest correlation under the condition that the reception order is changed and duplicate association does not occur, that is, the packets having the maximum correlation value exceeding the threshold value are associated with each other. ..

上記の例では、受信モジュール13−1においてパケットの受信欠損が発生した場合について、パケットの紐付け誤りの訂正を説明したが、受信モジュール13−1においてノイズ等に起因する余分なパケットを受信する場合もある。この場合、ノイズ等に起因する余分なパケットは、周知の方法で検出できるので、検出した余分なパケットを除外して、パケットの紐付けの誤りを上記と同様に訂正できる。更に、両方の受信モジュール13−1,13−2においてパケットの受信欠損が発生した場合、或いは、ノイズ等に起因する余分なパケットを受信した場合にも、パケットの紐付けの誤りを上記と同様に訂正できる。 In the above example, correction of packet association error has been described for the case where packet reception loss occurs in the receiving module 13-1, but the receiving module 13-1 receives an extra packet due to noise or the like. In some cases. In this case, since the extra packet caused by noise or the like can be detected by a well-known method, the detected extra packet can be excluded and the packet association error can be corrected in the same manner as described above. Further, when a packet reception defect occurs in both reception modules 13-1 and 13-2, or when an extra packet due to noise or the like is received, the packet association error is the same as described above. Can be corrected to.

各受信モジュールの処理時間の差分により、パケットの紐付け時に関連付けられた受信パケットの受信順が入れ替わったとしても、従来技術ではそのような受信パケットの受信順の入れ替わりを検出することができない。 Even if the reception order of the received packets associated with the packet association is changed due to the difference in the processing time of each receiving module, the conventional technology cannot detect such a change in the receiving order of the received packets.

これに対し、本実施例によれば、時間軸のスケールを伸縮した一方の受信信号系列を当該時間軸上でオフセット(スライド)し、2つの受信信号系列を紐付けた後、揺さぶりにより再度紐付けをする処理を行う。揺さぶりを行うことで、受信順の入れ替えや、重複紐付けが発生しない条件下で最も相関が高いパケット同士を紐付けることができるので、パケットの紐付けの誤りを訂正できる。 On the other hand, according to this embodiment, one of the received signal sequences whose scale on the time axis is expanded or contracted is offset (sliding) on the time axis, the two received signal sequences are linked, and then the strings are linked again by shaking. Perform the process of attaching. By shaking, it is possible to change the reception order and link the packets with the highest correlation under the condition that duplicate linking does not occur, so that it is possible to correct the packet linking error.

図13において、ステップS51〜S53の処理と、ステップS54〜S59の処理を配列NewRS2のループ分実行する処理とを、配列NewRS1のループ分実行した後、処理は図5に示す電波同定処理に戻ってステップS6の処理へ進む。 In FIG. 13, after the processing of steps S51 to S53 and the processing of steps S54 to S59 are executed for the loop of the array NewRS2 for the loop of the array NewRS1, the processing returns to the radio wave identification process shown in FIG. Then, the process proceeds to step S6.

図5において、ステップS6では、サーバ14のCPUが、電波同定結果をPC15に出力する電波同定結果出力処理を行う。図17は、電波同定結果出力処理の一例を説明するフローチャートである。図17において、ステップS61では、サーバ14のCPUが、配列Identifiedの値を読み込む。ステップS62では、サーバ14のCPUが、配列Identifiedの同一行に記載された変数NewRS1Index, NewRS2Indexが紐付けられたパケットであると判定して電波同定結果をPC15に出力し、処理は図5に示す電波同定処理に戻り、電波同定処理は終了する。PC15では、CPU151がサーバ14からの電波同定結果を例えば表示装置155に表示する。 In FIG. 5, in step S6, the CPU of the server 14 performs radio wave identification result output processing for outputting the radio wave identification result to the PC 15. FIG. 17 is a flowchart illustrating an example of radio wave identification result output processing. In FIG. 17, in step S61, the CPU of the server 14 reads the value of the array Identified. In step S62, the CPU of the server 14 determines that the variables NewRS1Index and NewRS2Index described in the same line of the array Identified are associated packets, outputs the radio wave identification result to the PC 15, and the process is shown in FIG. The process returns to the radio wave identification process, and the radio wave identification process ends. In the PC 15, the CPU 151 displays the radio wave identification result from the server 14 on, for example, the display device 155.

このように、サーバ14のCPUは、ステップS6の電波同定結果出力処理を実行することで、電波同定結果をPC15に出力する判定手段として機能する。判定手段は、最大値になる相関値が閾値を超えている場合に、紐付けられたパケットは第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する。 In this way, the CPU of the server 14 functions as a determination means for outputting the radio wave identification result to the PC 15 by executing the radio wave identification result output process in step S6. When the correlation value that becomes the maximum value exceeds the threshold value, the determination means determines that the associated packet is a reception signal sequence of radio waves received by the first and second receiving nodes from the same transmission source. Output the radio wave identification result.

第1実施例によれば、各受信モジュールにおけるタイムスタンプの個体差を補償することができる。 According to the first embodiment, it is possible to compensate for individual differences in time stamps in each receiving module.

(第2実施例)
図18は、第2実施例における電波環境可視化システムの一例を示す図である。図18中、図3と同一部分には同一符号を付し、その説明は省略する。図18に示す電波環境可視化システム11−2は、送信局12と、受信モジュール13−1,13−2と、サーバ14と、PC15を有する。受信モジュール13−1,13−2のうち、少なくとも一方は、送信局12のIDをサーバ14に送信する機能を有する。第2実施例では、サーバ14は、PC15からのID紐付けを含む電波同定処理を要求するコマンドに応答して、受信モジュール13−1,13−2が測定した信号受信強度RSSI,RSSIを、送信局12のIDと関連付ける、即ち、紐付ける。PC15は、操作者による入力装置154の操作に応答して、ID紐付けを含む電波同定処理を要求するコマンドを発行する。送信局12のIDは、この例では受信モジュール13−1が検出する。また、サーバ14は、電波同定結果及びID紐付け結果をPC15へ出力する。
(Second Example)
FIG. 18 is a diagram showing an example of the radio wave environment visualization system in the second embodiment. In FIG. 18, the same parts as those in FIG. 3 are designated by the same reference numerals, and the description thereof will be omitted. The radio wave environment visualization system 11-2 shown in FIG. 18 includes a transmitting station 12, receiving modules 13-1 and 13-2, a server 14, and a PC 15. At least one of the receiving modules 13-1 and 13-2 has a function of transmitting the ID of the transmitting station 12 to the server 14. In the second embodiment, the server 14 responds to the command requesting the radio wave identification process including the ID association from the PC 15, and the signal reception intensities RSSI 1 and RSSI 2 measured by the reception modules 13-1 and 13-2. Is associated with, that is, associated with the ID of the transmitting station 12. In response to the operation of the input device 154 by the operator, the PC 15 issues a command requesting radio wave identification processing including ID association. The ID of the transmitting station 12 is detected by the receiving module 13-1 in this example. Further, the server 14 outputs the radio wave identification result and the ID association result to the PC 15.

図19は、第2実施例における電波同定処理の一例を説明するフローチャートである。図19中、図5と同一ステップには同一符号を付し、その説明は省略する。この例では、サーバ14は、PC15からのID紐付けを含む電波同定処理を要求するコマンドを受信するので、ステップS6の後、処理はステップS7へ進む。なお、電波同定処理を要求するコマンドがID紐付けを含むか否かを判定するステップを設け、判定結果がYESである場合のみ処理がステップS7へ進み、判定結果がNOである場合はステップS7は実行せずにステップS6の後に処理を終了しても良い。図19において、ステップS7では、サーバ14のCPUが、受信モジュール13−1,13−2から受信した信号受信強度RSSI,RSSIと、送信局12のIDとのID紐付け結果をPC15に出力するID紐付け結果出力処理を行い、処理は終了する。 FIG. 19 is a flowchart illustrating an example of the radio wave identification process in the second embodiment. In FIG. 19, the same steps as those in FIG. 5 are designated by the same reference numerals, and the description thereof will be omitted. In this example, since the server 14 receives the command requesting the radio wave identification process including the ID association from the PC 15, the process proceeds to step S7 after step S6. A step of determining whether or not the command requesting the radio wave identification process includes ID association is provided, and the process proceeds to step S7 only when the determination result is YES, and step S7 when the determination result is NO. May end the process after step S6 without executing. In FIG. 19, in step S7, the CPU of the server 14 transfers the ID association result of the signal reception intensities RSSI 1 and RSSI 2 received from the reception modules 13-1 and 13-2 to the ID of the transmission station 12 to the PC 15. The ID association result output process to be output is performed, and the process ends.

図20は、ID組付け結果出力処理の一例を説明するフローチャートである。図20において、ステップS71では、サーバ14のCPUが、配列Identifiedの値を読み込む。ステップS72では、サーバ14のCPUが、配列Identifiedの同一行に記載された変数NewRS1Index, NewRS2Indexが紐付けられたパケットであると判定する。ステップS73では、サーバ14のCPUが、配列Identifiedの変数NewRS2Indexと同一行に記載された、送信局12のソースアドレスの一例であるMAC(Media Access Control)情報が紐付けられたIDであると判定してID紐付け結果をPC15へ出力する。ステップS73の後、処理は図5に示す電波同定処理に戻り、電波同定処理は終了する。送信局12のMAC情報は、送信源アドレス(又は、ソースアドレス)情報の一例である。PC15では、CPU151がサーバ14からの電波同定結果及びID紐付け結果を例えば表示装置155に表示する。 FIG. 20 is a flowchart illustrating an example of ID assembly result output processing. In FIG. 20, in step S71, the CPU of the server 14 reads the value of the array Identified. In step S72, the CPU of the server 14 determines that the packet is associated with the variables NewRS1Index and NewRS2Index described in the same line of the array Identified. In step S73, the CPU of the server 14 determines that the ID is associated with MAC (Media Access Control) information, which is an example of the source address of the transmitting station 12, described in the same line as the variable NewRS2Index of the array Identified. Then, the ID association result is output to the PC 15. After step S73, the process returns to the radio wave identification process shown in FIG. 5, and the radio wave identification process ends. The MAC information of the transmitting station 12 is an example of the source address (or source address) information. In the PC 15, the CPU 151 displays the radio wave identification result from the server 14 and the ID association result on, for example, the display device 155.

このように、サーバ14のCPUは、ステップS7のID紐付け結果出力処理を実行することで、受信モジュール13−1,13−2から受信した信号受信強度RSSI,RSSIと、送信局12のIDとのID紐付け結果を出力する出力手段として機能する。出力手段は、電波同定結果と、第1及び第2の受信信号系列の受信強度を受信モジュール13−1又は受信モジュール13−2で検出した送信局12のID情報と紐付けたID紐付け結果とを出力する。 In this way, the CPU of the server 14 executes the ID association result output process in step S7 to receive the signal reception strengths RSSI 1 and RSSI 2 and the transmission station 12 received from the reception modules 13-1 and 13-2. It functions as an output means for outputting the ID association result with the ID of. The output means is an ID linking result in which the radio wave identification result and the ID information of the transmitting station 12 detected by the receiving module 13-1 or the receiving module 13-2 are linked to the reception strength of the first and second received signal sequences. And output.

第2実施例によれば、各受信モジュールにおけるタイムスタンプの個体差を補償することができると共に、場所毎の電波の受信信号強度を送信源アドレス情報と紐付けることができる。 According to the second embodiment, it is possible to compensate for the individual difference of the time stamp in each receiving module, and it is possible to associate the received signal strength of the radio wave for each place with the source address information.

(第3実施例)
ところで、図5に示す第1実施例の電波同定処理を、各受信モジュール対について繰り返すことで、第3実施例のように、3台以上の受信モジュールが設けられている場合にも上記の電波同定処理を拡張できる。タイムスタンプの個体差を補償する処理と、パケットの紐付け処理は、2台の受信モジュール間で行なう処理である。これらの処理と同様の処理を異なる組合せの2台の受信モジュール間で繰り返すことにより、3台以上の受信モジュールのタイムスタンプの個体差を補償し、パケットを紐付けることができる。
(Third Example)
By the way, by repeating the radio wave identification process of the first embodiment shown in FIG. 5 for each receiving module pair, the above radio wave is provided even when three or more receiving modules are provided as in the third embodiment. The identification process can be extended. The process of compensating for individual differences in time stamps and the process of linking packets are processes performed between two receiving modules. By repeating the same processing as these processing between two receiving modules of different combinations, it is possible to compensate for individual differences in the time stamps of three or more receiving modules and associate packets.

図21は、第3実施例における電波環境可視化システムの一例を示す図である。図21中、図3及び図18と同一部分には同一符号を付し、その説明は省略する。図21に示す電波環境可視化システム11−3は、送信局12と、受信モジュール13−1〜13−n(nは3以上の自然数)と、サーバ14と、PC15を有する。この例では、受信モジュール13−1以外の受信モジュール13−2〜13−nは、SRDモジュールである。第3実施例では、サーバ14は、PC15からの電波同定処理を要求するコマンドがID紐付けを含む場合には、一対の受信モジュールが測定した信号受信強度と、ある受信モジュールが検出した送信局12のIDとを関連付ける、即ち、紐付ける。この例では、ある受信モジュールとは、受信モジュール13−1である。そして、サーバ14は、このようなID紐付けを、例えば隣り合う一対の受信モジュール等の各受信モジュール対が測定した信号受信強度について行い、最終的には全ての受信モジュールについて、電波同定結果及び/又はID紐付け結果をPC15へ出力する。つまり、サーバ14は、電波同定結果及びID紐付け結果のうち少なくとも一方をPC15へ出力することができる。 FIG. 21 is a diagram showing an example of the radio wave environment visualization system according to the third embodiment. In FIG. 21, the same parts as those in FIGS. 3 and 18 are designated by the same reference numerals, and the description thereof will be omitted. The radio wave environment visualization system 11-3 shown in FIG. 21 includes a transmitting station 12, receiving modules 13-1 to 13-n (n is a natural number of 3 or more), a server 14, and a PC 15. In this example, the receiving modules 13-2 to 13-n other than the receiving module 13-1 are SRD modules. In the third embodiment, when the command requesting the radio wave identification process from the PC 15 includes ID association, the server 14 has the signal reception strength measured by the pair of reception modules and the transmission station detected by a certain reception module. It associates with, that is, associates with 12 IDs. In this example, a receiving module is a receiving module 13-1. Then, the server 14 performs such ID association for the signal reception strength measured by each reception module pair such as a pair of adjacent reception modules, and finally, the radio wave identification result and the radio wave identification result for all the reception modules. / Or the ID association result is output to the PC15. That is, the server 14 can output at least one of the radio wave identification result and the ID association result to the PC 15.

図22は、第3実施例における電波同定処理の一例を説明するフローチャートである。図22において、ステップS1−1では、サーバ14のCPUが、一対の受信モジュール13−1,13−2のログデータを読み込み、例えば図5に示すステップS2〜S4と同様の処理を行う。ステップS7−1では、サーバ14のCPUが、受信モジュール13−1から受信した信号受信強度RSSIと受信モジュール13−2から受信した信号受信強度RSSIの、送信局12のIDとのID紐付け結果を出力するID紐付け結果出力処理を行う。ステップS7−1の処理は、例えば図5に示すステップS5の処理と同様である。図22において、例えばステップS7−1に記載された「1&2」は、信号受信強度RSSI1,RSSIの、送信局12のIDとのID紐付け結果を出力するID紐付け結果出力処理であることを示す。同様に、例えば後述するステップS7−n−1に記載された「n−1&n」は、信号受信強度RSSIn-1,RSSInの、送信局12のIDとのID紐付け結果を出力するID紐付け結果出力処理であることを示す。ID紐付け結果出力処理は、例えば図20と共に説明したID紐付け結果出力処理と同様で良い。 FIG. 22 is a flowchart illustrating an example of the radio wave identification process in the third embodiment. In FIG. 22, in step S1-1, the CPU of the server 14 reads the log data of the pair of receiving modules 13-1 and 13-2, and performs the same processing as in steps S2 to S4 shown in FIG. 5, for example. In step S7-1, the CPU of the server 14 associates the signal reception strength RSSI 1 received from the reception module 13-1 with the signal reception strength RSSI 2 received from the reception module 13-2 with the ID of the transmission station 12. The ID association result output processing for outputting the attachment result is performed. The process of step S7-1 is, for example, the same as the process of step S5 shown in FIG. In FIG. 22, for example, “1 & 2” described in step S7-1 is an ID association result output process for outputting the ID association result of the signal reception strength RSSI 1 and RSSI 2 with the ID of the transmitting station 12. Show that. Similarly, for example, “n-1 & n” described in step S7-n-1 described later is an ID that outputs the ID association result of the signal reception intensities RSSI n-1 and RSSI n with the ID of the transmitting station 12. Indicates that the linking result output processing is performed. The ID association result output process may be the same as the ID association result output process described with reference to FIG. 20, for example.

次に、ステップS1−2では、サーバ14のCPUが、一対の受信モジュール13−2,13−3のログデータを読み込み、例えば図5に示すステップS2〜S4と同様の処理を行う。ステップS7−2では、サーバ14のCPUが、受信モジュール13−2から受信した信号受信強度RSSIと受信モジュール13−3から受信した信号受信強度RSSIの、送信局12のIDとのID紐付け結果を出力するID紐付け結果出力処理を行う。ステップS7−2の処理は、例えば図5に示すステップS5の処理と同様である。以下、上記の如き処理を繰り返し、ステップS1−n−1では、サーバ14のCPUが、一対の受信モジュール13−n−1,13−nのログデータを読み込み、例えば図5に示すステップS2〜S4と同様の処理を行う。ステップS7−n−1では、サーバ14のCPUが、受信モジュール13−n−1から受信した信号受信強度RSSIn−1と受信モジュール13−nから受信した信号受信強度RSSIの、送信局12のIDとのID紐付け結果を出力するID紐付け結果出力処理を行う。ステップS7−n−1の処理は、例えば図5に示すステップS5の処理と同様である。次に、ステップS7−nでは、サーバ14のCPUが、各対の受信モジュールが測定した信号受信強度について、例えば図5に示すステップS6と同様の処理を行い求めた電波同定結果をマージする。また、ステップS7−nでは、サーバ14のCPUが、各対の受信モジュールが測定した信号受信強度について求めたID紐付け結果をマージする。更に、ステップS7−nでは、サーバ14のCPUが、マージした電波同定結果及び/又はマージしたID紐付け結果とを、全ての受信モジュール13−1〜13−nの電波同定結果及び/又はID紐付け結果としてPC15へ出力する。 Next, in step S1-2, the CPU of the server 14 reads the log data of the pair of receiving modules 13-2 and 13-3, and performs the same processing as in steps S2 to S4 shown in FIG. 5, for example. In step S7-2, the CPU of the server 14 links the signal reception strength RSSI 2 received from the reception module 13-2 and the signal reception strength RSSI 3 received from the reception module 13-3 with the ID of the transmission station 12. The ID association result output processing for outputting the attachment result is performed. The process of step S7-2 is the same as the process of step S5 shown in FIG. 5, for example. Hereinafter, the above processing is repeated, and in step S1-n-1, the CPU of the server 14 reads the log data of the pair of receiving modules 13-n-1, 13-n, and for example, steps S2 to 2 shown in FIG. The same process as in S4 is performed. In step S7-n-1, the CPU of the server 14 transmits the signal reception strength RSSI n-1 received from the reception module 13-n -1 and the signal reception strength RSSI n received from the reception module 13-n. The ID association result output process for outputting the ID association result with the ID of is performed. The process of step S7-n-1 is, for example, the same as the process of step S5 shown in FIG. Next, in step S7-n, the CPU of the server 14 merges the radio wave identification results obtained by performing the same processing as in step S6 shown in FIG. 5 for the signal reception strength measured by each pair of reception modules. Further, in step S7-n, the CPU of the server 14 merges the ID association results obtained for the signal reception strength measured by each pair of reception modules. Further, in step S7-n, the CPU of the server 14 transfers the merged radio wave identification result and / or the merged ID association result to the radio wave identification results and / or IDs of all the receiving modules 13-1 to 13-n. Output to PC15 as the linking result.

このように、サーバ14のCPUは、図22に示す処理を実行することで、上記決定手段及び上記紐付け手段が、3台以上の受信モジュール13−1〜13−nが受信した受信信号系列のうち、各受信モジュール対が受信した受信信号系列に対して第1及び第2の受信信号系列に対する処理を繰り返す。また、上記判定手段が、各受信ノード対に対する電波同定結果をマージする。 In this way, the CPU of the server 14 executes the process shown in FIG. 22, and the receiving signal sequence received by the three or more receiving modules 13-1 to 13-n by the determining means and the associating means. Among the received signal sequences received by each receiving module pair, the processing for the first and second received signal sequences is repeated. In addition, the determination means merges the radio wave identification results for each receiving node pair.

第3実施例によれば、3台以上の受信モジュールが設けられている場合にも電波同定処理を拡張することができる。また、必要に応じて、各受信モジュールにおけるタイムスタンプの個体差を補償したり、場所毎の電波の受信信号強度を送信源アドレス情報と紐付けることができる。更に、複数のID紐付け結果をマージするので、ID紐付け結果の精度を向上することができる。 According to the third embodiment, the radio wave identification process can be extended even when three or more receiving modules are provided. Further, if necessary, it is possible to compensate for individual differences in time stamps in each receiving module, and to associate the received signal strength of radio waves for each location with the source address information. Further, since a plurality of ID association results are merged, the accuracy of the ID association results can be improved.

(第4実施例)
次に、第4実施例について説明する。第4実施例では、上記第1乃至第3実施例のいずれかで出力した電波同定結果に基づいて、電波環境可視化処理を行う。
(Fourth Example)
Next, a fourth embodiment will be described. In the fourth embodiment, the radio wave environment visualization process is performed based on the radio wave identification result output in any of the first to third embodiments.

図23は、電波環境可視化処理の一例を説明するフローチャートである。図23に示す電波環境可視化処理は、例えばPC15のCPU151が、メモリ152が記憶するプログラムを実行することで実行可能である。図23に示す電波環境可視化処理を開始するには、例えばPC15の操作者が入力装置154を操作して電波環境可視化処理を要求するコマンドを発行する。 FIG. 23 is a flowchart illustrating an example of the radio wave environment visualization process. The radio wave environment visualization process shown in FIG. 23 can be executed, for example, by the CPU 151 of the PC 15 executing a program stored in the memory 152. To start the radio wave environment visualization process shown in FIG. 23, for example, the operator of the PC 15 operates the input device 154 to issue a command requesting the radio wave environment visualization process.

図23において、ステップS91では、PC15のCPU151が、電波環境可視化処理を要求するコマンドに応答して電波同定処理を要求するコマンドをサーバ14に送信する。また、ステップS91では、PC15のCPU151が、送信した電波同定処理を要求するコマンドに応答した、サーバ14からの電波同定結果出力を受信する。ステップS92では、PC15のCPU151が、電波同定結果出力に基づいて、受信モジュールの場所毎の電波の受信信号強度(RSSI)や、トラフィックのレベルを色で表現したヒートマップを作成する。ヒートマップの作成自体は、例えば周知の方法を採用可能である。また、送信局12のIDを出力する受信モジュール(例えば、受信モジュール13−1)が存在すれば、例えばサーバ14が受信した当該IDをヒートマップ上に表示する等して、周辺の電波環境の調査や不具合の原因特定に利用することができる。ステップS93では、PC15のCPU151が、例えば図24に示す如きヒートマップをPC15の表示装置155に表示し、処理は終了する。PC15のCPU141は、必要に応じてヒートマップ等をPC15のメモリ152に格納する。 In FIG. 23, in step S91, the CPU 151 of the PC 15 transmits a command requesting the radio wave identification process to the server 14 in response to the command requesting the radio wave environment visualization process. Further, in step S91, the CPU 151 of the PC 15 receives the radio wave identification result output from the server 14 in response to the transmitted command requesting the radio wave identification process. In step S92, the CPU 151 of the PC 15 creates a heat map in which the received signal strength (RSSI) of the radio wave for each location of the receiving module and the traffic level are expressed in color based on the radio wave identification result output. For the heat map creation itself, for example, a well-known method can be adopted. Further, if there is a receiving module (for example, receiving module 13-1) that outputs the ID of the transmitting station 12, for example, the ID received by the server 14 is displayed on the heat map to display the surrounding radio wave environment. It can be used for investigation and identification of the cause of defects. In step S93, the CPU 151 of the PC 15 displays a heat map as shown in FIG. 24 on the display device 155 of the PC 15, and the process ends. The CPU 141 of the PC 15 stores a heat map or the like in the memory 152 of the PC 15 as needed.

図24は、ヒートマップの一例を示す図である。図24中、縦軸はXY座標系におけるY方向の距離(m)を示し、横軸はXY座標系におけるX方向の距離(m)を示す。図24は、RSSIのヒートマップを示し、全てのアクセスポイント(AP:Access Point)からのRSSIレベルの最大値(dBm)を、各受信モジュールについてハーフトーンで表す。図24において、RSSIレベルは、ハーフトーンが明るい程低く、暗い程高いことを示す。この例では、受信モジュールは、Wi−FiモジュールWi-Fi001〜Wi-Fi005と、BLE(Bluetooth (登録商標)Low Energy)モジュールBLE001〜BLE004と、BT(Bluetooth(登録商標))モジュールBT001〜BT003を含む。 FIG. 24 is a diagram showing an example of a heat map. In FIG. 24, the vertical axis shows the distance (m) in the Y direction in the XY coordinate system, and the horizontal axis shows the distance (m) in the X direction in the XY coordinate system. FIG. 24 shows an RSSI heat map showing the maximum value (dBm) of the RSSI level from all access points (APs) in halftone for each receiving module. In FIG. 24, the RSSI level is lower as the halftone is brighter and higher as the halftone is darker. In this example, the receiving modules include Wi-Fi modules Wi-Fi001 to Wi-Fi005, BLE (Bluetooth (registered trademark) Low Energy) modules BLE001 to BLE004, and BT (Bluetooth (registered trademark)) modules BT001 to BT003. Including.

第4実施例によれば、複数の受信モジュールが受信した同一の送信源から送信された電波を同定できるので、各受信モジュールにおける電波受信の時間差から各受信モジュールから送信源までの距離を推定することができる。このため、電界強度分布のヒートマップの描画等の電波環境可視化に、送信源の位置の推定結果の表示を含めることもでき、ヒートマップを用いて電波環境の調査や、不具合の原因特定を迅速に行なえる。 According to the fourth embodiment, since the radio waves transmitted from the same source received by the plurality of receiving modules can be identified, the distance from each receiving module to the transmitting source is estimated from the time difference of radio wave reception in each receiving module. be able to. For this reason, it is possible to include the display of the estimation result of the position of the transmission source in the visualization of the radio wave environment such as drawing the heat map of the electric field strength distribution, and the heat map can be used to quickly investigate the radio wave environment and identify the cause of the problem. You can go to.

このように、上記第1乃至第3実施例における電波同定装置は、電波環境可視化技術に適用可能である。このため、上記第1乃至第3実施例における電波同定装置は、例えば工場等に設置された無線センサネットワークが動作不良を起こした際に、電波環境可視化技術に適用して原因を迅速に究明する場合等に好適である。 As described above, the radio wave identification device in the first to third embodiments can be applied to the radio wave environment visualization technique. Therefore, the radio wave identification device in the first to third embodiments is applied to the radio wave environment visualization technology to quickly investigate the cause when, for example, a wireless sensor network installed in a factory or the like malfunctions. It is suitable for cases and the like.

(変形例)
なお、変形例では、上記の各実施例において、サーバ14は、PC15の処理の少なくとも一部、即ち、PC15の処理の一部又は全部を実行しても良い。サーバ14がPC15の全部の処理を実行する場合、PC15は省略可能である。同様に、上記の各実施例において、PC15は、サーバ14の処理の少なくとも一部、即ち、サーバ14の処理の一部又は全部を実行しても良い。PC15がサーバ14の全部の処理を実行する場合、サーバ14は省略可能である。つまり、上記の各実施例において、サーバ14とPC15を、サーバ14とPC15の処理を実行する単一のコンピュータで置き換えても良い。
(Modification example)
In the modified example, in each of the above embodiments, the server 14 may execute at least a part of the processing of the PC 15, that is, a part or all of the processing of the PC 15. If the server 14 executes all the processing of the PC 15, the PC 15 can be omitted. Similarly, in each of the above embodiments, the PC 15 may execute at least a part of the processing of the server 14, that is, a part or all of the processing of the server 14. When the PC 15 executes all the processing of the server 14, the server 14 can be omitted. That is, in each of the above embodiments, the server 14 and the PC 15 may be replaced with a single computer that executes the processing of the server 14 and the PC 15.

なお、上記の各実施例に付されている第1乃至第4なる連番は、好ましい実施例の優先順位を表すものではない。また、当業者であれば、数々の変形や改変が可能であることは明らかである。 The first to fourth serial numbers attached to each of the above embodiments do not represent the priority of the preferred embodiments. It is also clear that those skilled in the art can make numerous modifications and modifications.

以上の実施例を含む実施形態に関し、更に以下の付記を開示する。
(付記1)
第1の受信ノードが受信した第1の受信信号系列と、第2の受信ノードが受信した第2の受信信号系列のうち、前記第2の受信信号系列の時間軸のスケールを時間伸縮率で伸縮しながら前記第1の受信信号系列との相関値を算出し、前記第2の受信信号系列を時間軸上で時間オフセットだけオフセットして前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する決定手段と、
前記相関値が最大値になる前記時間伸縮率と前記時間オフセットを前記第2の受信信号系列の時刻情報に適用することで、前記第1の受信信号系列及び前記第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付ける紐付け手段と、
前記最大値になる前記相関値が閾値を超えている場合に、紐付けられたパケットは前記第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する判定手段と、
を備えたことを特徴とする、電波同定装置。
(付記2)
前記第2の受信信号系列の先頭データと末尾データの時刻情報が、前記第1の受信信号系列の先頭データと末尾データの時刻情報より常にガード時間以上広くなるようなクリッピングタイミングで前記第1の受信信号系列及び前記第2の受信信号系列をクリッピングする抽出手段を更に備え、
前記決定手段は、前記抽出手段が抽出した受信信号系列から前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する
ことを特徴とする、付記1記載の電波同定装置。
(付記3)
前記紐付け手段は、紐付け候補のパケットの受信順が異なっている場合、又は、複数のパケットに重複紐付けが発生している場合に、揺さぶりを行い、パケットの順番は入れ替わらないというルールを適用してパケットの紐付けの誤りを訂正し、
前記揺さぶりは、前記紐付け候補のパケットを時間軸上で一定範囲だけ振って微調整することで、パケットの紐付けの誤りを訂正する
ことを特徴とする、付記1又は2記載の電波同定装置。
(付記4)
前記電波同定結果と、前記第1の受信信号系列及び前記第2の受信信号系列の受信強度を前記第1の受信ノード又は前記第2の受信ノードで検出した前記送信源のID情報と紐付けたID紐付け結果とを出力する出力手段を更に備えたことを特徴とする、付記1乃至3のいずれか1項記載の電波同定装置。
(付記5)
前記決定手段及び前記紐付け手段は、3台以上の受信ノードが受信した受信信号系列のうち、各受信ノード対が受信した受信信号系列に対して前記第1及び第2の受信信号系列に対する処理を繰り返し、
前記判定手段は、各受信ノード対に対する電波同定結果をマージすることを特徴とする、付記1乃至4のいずれか1項記載の電波同定装置。
(付記6)
第1の受信ノードが受信した第1の受信信号系列と、第2の受信ノードが受信した第2の受信信号系列のうち、前記第2の受信信号系列の時間軸のスケールを時間伸縮率で伸縮しながら前記第1の受信信号系列との相関値を算出し、前記第2の受信信号系列を時間軸上で時間オフセットだけオフセットして前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定し、
前記相関値が最大値になる前記時間伸縮率と前記時間オフセットを前記第2の受信信号系列の時刻情報に適用することで、前記第1の受信信号系列及び前記第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付け、
前記最大値になる前記相関値が閾値を超えている場合に、紐付けられたパケットは前記第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する、
処理をコンピュータが実行することを特徴とする、電波同定方法。
(付記7)
前記第2の受信信号系列の先頭データと末尾データの時刻情報が、前記第1の受信信号系列の先頭データと末尾データの時刻情報より常にガード時間以上広くなるようなクリッピングタイミングで受信信号系列を抽出する、
処理を前記コンピュータが更に実行し、
前記決定する処理は、前記抽出する処理が抽出した受信信号系列から前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する
ことを特徴とする、付記6記載の電波同定方法。
(付記8)
前記紐付ける処理は、紐付け候補のパケットの受信順が異なっている場合、又は、複数のパケットに重複紐付けが発生している場合に、揺さぶりを行い、パケットの順番は入れ替わらないというルールを適用してパケットの紐付けの誤りを訂正し、
前記揺さぶりは、前記紐付け候補のパケットを時間軸上で一定範囲だけ振って微調整することで、パケットの紐付けの誤りを訂正する
ことを特徴とする、付記6又は7記載の電波同定方法。
(付記9)
前記電波同定結果と、前記第1の受信信号系列及び前記第2の受信信号系列の受信強度を前記第1の受信ノード又は前記第2の受信ノードで検出した前記送信源のID情報と紐付けたID紐付け結果とを出力する、
処理を前記コンピュータに更に実行させることを特徴とする、付記6乃至8のいずれか1項記載の電波同定方法。
(付記10)
前記決定する処理及び前記紐付ける処理は、3台以上の受信ノードが受信した受信信号系列のうち、各受信ノード対が受信した受信信号系列に対して前記第1及び第2の受信信号系列に対する処理を繰り返し、
前記判定する処理は、各受信ノード対に対する電波同定結果をマージすることを特徴とする、付記6乃至9のいずれか1項記載の電波同定方法。
(付記11)
受信した第1の受信信号系列を出力する第1の受信ノードと、
前記第1の受信ノードとは異なる場所に設置され、前記第1の受信ノードの内部時計とは独立した内部時計を有し、受信した第2の受信信号系列を出力する第2の受信ノードと、
前記第1の受信信号系列及び前記第2の受信信号系列に基づいて電波同定結果を出力する電波同定装置と、
前記電波同定結果に基づいて、前記第1及び第2の受信ノードの場所毎の電波の受信信号強度とトラフィックのレベルを色で表現したヒートマップを作成する装置と、
を備え、
前記電波同定装置は、
前記第2の受信信号系列の時間軸のスケールを時間伸縮率で伸縮しながら前記第1の受信信号系列との相関値を算出し、前記第2の受信信号系列を時間軸上で時間オフセットだけオフセットして前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する決定手段と、
前記相関値が最大値になる前記時間伸縮率と前記時間オフセットを前記第2の受信信号系列の時刻情報に適用することで、前記第1の受信信号系列及び前記第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付ける紐付け手段と、
前記最大値になる前記相関値が閾値を超えている場合に、紐付けられたパケットは前記第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する判定手段と、
を有することを特徴とする、電波環境可視化システム。
(付記12)
前記ヒートマップを作成する前記装置は、前記送信源の位置情報及びID情報を前記ヒートマップ上に表示する手段を有することを特徴とする、付記11記載の電波環境可視化システム。
(付記13)
前記電波同定装置及び前記ヒートマップを作成する前記装置は、単一のコンピュータで形成されていることを特徴とする、付記11又は12記載の電波環境可視化システム。
The following additional notes will be further disclosed with respect to the embodiments including the above embodiments.
(Appendix 1)
Of the first received signal sequence received by the first receiving node and the second received signal sequence received by the second receiving node, the scale of the time axis of the second received signal sequence is scaled by the time expansion / contraction rate. While expanding and contracting, the correlation value with the first received signal sequence is calculated, and the second received signal sequence is offset by a time offset on the time axis to maximize the correlation value. The means of determining the time offset and
By applying the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value to the time information of the second received signal sequence, the first received signal sequence and the second received signal sequence are respectively. A means of associating packets having similar transmission start times, packet lengths, and signal strengths with each other on a one-to-one basis by comparing each packet of each packet.
When the correlation value that becomes the maximum value exceeds the threshold value, it is determined that the associated packet is a reception signal sequence of radio waves received by the first and second receiving nodes from the same transmission source. Judgment means to output the radio wave identification result and
A radio wave identification device characterized by being equipped with.
(Appendix 2)
The first data at a clipping timing such that the time information of the head data and the tail data of the second received signal sequence is always wider than the time information of the head data and the tail data of the first received signal series by a guard time or more. Further provided with an extraction means for clipping the received signal sequence and the second received signal sequence.
The radio wave identification device according to Appendix 1, wherein the determination means determines the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value from the received signal sequence extracted by the extraction means.
(Appendix 3)
The rule that the linking means shakes when the reception order of the packet of the linking candidate is different, or when duplicate linking occurs in a plurality of packets, and the order of the packets is not changed. Correct the error of packet linking by applying
The radio wave identification device according to Appendix 1 or 2, wherein the shaking corrects an error in packet association by shaking the packet of the association candidate by a certain range on the time axis and making fine adjustments. ..
(Appendix 4)
The radio wave identification result is associated with the reception strength of the first reception signal sequence and the second reception signal sequence with the ID information of the transmission source detected by the first reception node or the second reception node. The radio wave identification device according to any one of Supplementary note 1 to 3, further comprising an output means for outputting the ID association result.
(Appendix 5)
The determining means and the associating means process the first and second received signal sequences for the received signal sequences received by each receiving node pair among the received signal sequences received by the three or more receiving nodes. Repeat,
The radio wave identification device according to any one of Supplementary note 1 to 4, wherein the determination means merges the radio wave identification results for each receiving node pair.
(Appendix 6)
Of the first received signal sequence received by the first receiving node and the second received signal sequence received by the second receiving node, the scale of the time axis of the second received signal sequence is scaled by the time expansion / contraction rate. While expanding and contracting, the correlation value with the first received signal sequence is calculated, and the second received signal sequence is offset by a time offset on the time axis to maximize the correlation value. Determine the time offset,
By applying the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value to the time information of the second received signal sequence, the first received signal sequence and the second received signal sequence are respectively. Compare each packet that each packet has, and associate packets with similar transmission start time, packet length, and signal strength of each packet on a one-to-one basis.
When the correlation value that becomes the maximum value exceeds the threshold value, it is determined that the associated packet is a reception signal sequence of radio waves received by the first and second receiving nodes from the same transmission source. Output the radio wave identification result,
A radio wave identification method characterized in that a computer executes processing.
(Appendix 7)
The received signal sequence is set at a clipping timing such that the time information of the start data and the end data of the second received signal sequence is always wider than the time information of the start data and the end data of the first received signal sequence by a guard time or more. Extract,
The computer further executes the process,
The radio wave identification method according to Appendix 6, wherein the determination process determines the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value from the received signal sequence extracted by the extraction process.
(Appendix 8)
The rule that the linking process shakes when the receiving order of the packet of the linking candidate is different or when duplicate linking occurs in a plurality of packets, and the order of the packets is not changed. Correct the error of packet linking by applying
The radio wave identification method according to Appendix 6 or 7, wherein the shaking corrects an error in packet association by shaking the packet of the association candidate by a certain range on the time axis and making fine adjustments. ..
(Appendix 9)
The radio wave identification result is associated with the reception strength of the first reception signal sequence and the second reception signal sequence with the ID information of the transmission source detected by the first reception node or the second reception node. Outputs the ID link result
The radio wave identification method according to any one of Supplementary note 6 to 8, wherein the processing is further executed by the computer.
(Appendix 10)
The process of determining and the process of associating the first and second received signal sequences with respect to the received signal sequence received by each receiving node pair among the received signal sequences received by three or more receiving nodes. Repeat the process,
The radio wave identification method according to any one of Supplementary note 6 to 9, wherein the determination process merges the radio wave identification results for each receiving node pair.
(Appendix 11)
A first receiving node that outputs the first received signal sequence received, and
A second receiving node that is installed at a location different from the first receiving node, has an internal clock independent of the internal clock of the first receiving node, and outputs a second received signal sequence received. ,
A radio wave identification device that outputs a radio wave identification result based on the first received signal sequence and the second received signal sequence, and
Based on the radio wave identification result, a device for creating a heat map in which the received signal strength of the radio wave and the traffic level for each location of the first and second receiving nodes are expressed in color, and
With
The radio wave identification device is
The correlation value with the first received signal sequence is calculated while expanding and contracting the scale of the time axis of the second received signal sequence with the time expansion / contraction rate, and the second received signal sequence is offset only by the time on the time axis. A determination means for determining the time expansion / contraction rate and the time offset at which the correlation value is offset to the maximum value, and
By applying the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value to the time information of the second received signal sequence, the first received signal sequence and the second received signal sequence are respectively. A means of associating packets having similar transmission start times, packet lengths, and signal strengths with each other on a one-to-one basis by comparing each packet of each packet.
When the correlation value that becomes the maximum value exceeds the threshold value, it is determined that the associated packet is a reception signal sequence of radio waves received by the first and second receiving nodes from the same transmission source. Judgment means to output the radio wave identification result and
A radio wave environment visualization system characterized by having.
(Appendix 12)
The radio wave environment visualization system according to Appendix 11, wherein the device for creating the heat map has means for displaying the position information and ID information of the transmission source on the heat map.
(Appendix 13)
The radio wave environment visualization system according to Appendix 11 or 12, wherein the radio wave identification device and the device for creating the heat map are formed by a single computer.

以上、開示の電波同定装置、電波同定方法及び電波環境可視化システムを実施例により説明したが、本発明は上記実施例に限定されるものではなく、本発明の範囲内で種々の変形及び改良が可能であることは言うまでもない。 Although the disclosed radio wave identification device, radio wave identification method and radio wave environment visualization system have been described above by way of examples, the present invention is not limited to the above examples, and various modifications and improvements can be made within the scope of the present invention. It goes without saying that it is possible.

11−1〜11−3 電波環境可視化システム
12 送信局
13−1〜13−n 受信モジュール
14 サーバ
15 PC
151 CPU
152 メモリ
153 I/F
154 入力装置
155 表示装置
156 バス
11-1 to 11-3 Radio environment visualization system 12 Transmission station 13-1 to 13-n Reception module 14 Server 15 PC
151 CPU
152 Memory 153 I / F
154 Input device 155 Display device 156 Bus

Claims (5)

第1の受信ノードが受信した第1の受信信号系列と、第2の受信ノードが受信した第2の受信信号系列のうち、前記第2の受信信号系列の時間軸のスケールを時間伸縮率で伸縮しながら前記第1の受信信号系列との相関値を算出し、前記第2の受信信号系列を時間軸上で時間オフセットだけオフセットして前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する決定手段と、
前記相関値が最大値になる前記時間伸縮率と前記時間オフセットを前記第2の受信信号系列の時刻情報に適用することで、前記第1の受信信号系列及び前記第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付ける紐付け手段と、
前記最大値になる前記相関値が閾値を超えている場合に、紐付けられたパケットは前記第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する判定手段と、
を備えたことを特徴とする、電波同定装置。
Of the first received signal sequence received by the first receiving node and the second received signal sequence received by the second receiving node, the scale of the time axis of the second received signal sequence is scaled by the time expansion / contraction rate. While expanding and contracting, the correlation value with the first received signal sequence is calculated, and the second received signal sequence is offset by a time offset on the time axis to maximize the correlation value. The means of determining the time offset and
By applying the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value to the time information of the second received signal sequence, the first received signal sequence and the second received signal sequence are respectively. A means of associating packets having similar transmission start times, packet lengths, and signal strengths with each other on a one-to-one basis by comparing each packet of each packet.
When the correlation value that becomes the maximum value exceeds the threshold value, it is determined that the associated packet is a reception signal sequence of radio waves received by the first and second receiving nodes from the same transmission source. Judgment means to output the radio wave identification result and
A radio wave identification device characterized by being equipped with.
前記第2の受信信号系列の先頭データと末尾データの時刻情報が、前記第1の受信信号系列の先頭データと末尾データの時刻情報より常にガード時間以上広くなるようなクリッピングタイミングで前記第1の受信信号系列及び前記第2の受信信号系列をクリッピングする抽出手段を更に備え、
前記決定手段は、前記抽出手段が抽出した受信信号系列から前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する
ことを特徴とする、請求項1記載の電波同定装置。
The first data at a clipping timing such that the time information of the head data and the tail data of the second received signal sequence is always wider than the time information of the head data and the tail data of the first received signal series by a guard time or more. Further provided with an extraction means for clipping the received signal sequence and the second received signal sequence.
The radio wave identification device according to claim 1, wherein the determination means determines the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value from the received signal sequence extracted by the extraction means.
前記紐付け手段は、紐付け候補のパケットの受信順が異なっている場合、又は、複数のパケットに重複紐付けが発生している場合に、揺さぶりを行い、パケットの順番は入れ替わらないというルールを適用してパケットの紐付けの誤りを訂正し、
前記揺さぶりは、前記紐付け候補のパケットを時間軸上で一定範囲だけ振って微調整することで、パケットの紐付けの誤りを訂正する
ことを特徴とする、請求項1又は2記載の電波同定装置。
The rule that the linking means shakes when the reception order of the packets of the linking candidate is different, or when duplicate linking occurs in a plurality of packets, and the order of the packets is not changed. Correct the error of packet linking by applying
The radio wave identification according to claim 1 or 2, wherein the shaking corrects an error in packet association by shaking the packet of the association candidate by a certain range on the time axis and making fine adjustments. apparatus.
第1の受信ノードが受信した第1の受信信号系列と、第2の受信ノードが受信した第2の受信信号系列のうち、前記第2の受信信号系列の時間軸のスケールを時間伸縮率で伸縮しながら前記第1の受信信号系列との相関値を算出し、前記第2の受信信号系列を時間軸上で時間オフセットだけオフセットして前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定し、
前記相関値が最大値になる前記時間伸縮率と前記時間オフセットを前記第2の受信信号系列の時刻情報に適用することで、前記第1の受信信号系列及び前記第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付け、
前記最大値になる前記相関値が閾値を超えている場合に、紐付けられたパケットは前記第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する、
処理をコンピュータが実行することを特徴とする、電波同定方法。
Of the first received signal sequence received by the first receiving node and the second received signal sequence received by the second receiving node, the scale of the time axis of the second received signal sequence is scaled by the time expansion / contraction rate. While expanding and contracting, the correlation value with the first received signal sequence is calculated, and the second received signal sequence is offset by a time offset on the time axis to maximize the correlation value. Determine the time offset,
By applying the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value to the time information of the second received signal sequence, the first received signal sequence and the second received signal sequence are respectively. Compare each packet that each packet has, and associate packets with similar transmission start time, packet length, and signal strength of each packet on a one-to-one basis.
When the correlation value that becomes the maximum value exceeds the threshold value, it is determined that the associated packet is a reception signal sequence of radio waves received by the first and second receiving nodes from the same transmission source. Output the radio wave identification result,
A radio wave identification method characterized in that a computer executes processing.
受信した第1の受信信号系列を出力する第1の受信ノードと、
前記第1の受信ノードとは異なる場所に設置され、前記第1の受信ノードの内部時計とは独立した内部時計を有し、受信した第2の受信信号系列を出力する第2の受信ノードと、
前記第1の受信信号系列及び前記第2の受信信号系列に基づいて電波同定結果を出力する電波同定装置と、
前記電波同定結果に基づいて、前記第1及び第2の受信ノードの場所毎の電波の受信信号強度とトラフィックのレベルを色で表現したヒートマップを作成する装置と、
を備え、
前記電波同定装置は、
前記第2の受信信号系列の時間軸のスケールを時間伸縮率で伸縮しながら前記第1の受信信号系列との相関値を算出し、前記第2の受信信号系列を時間軸上で時間オフセットだけオフセットして前記相関値が最大値になる前記時間伸縮率及び前記時間オフセットを決定する決定手段と、
前記相関値が最大値になる前記時間伸縮率と前記時間オフセットを前記第2の受信信号系列の時刻情報に適用することで、前記第1の受信信号系列及び前記第2の受信信号系列の夫々が持つ1つ1つのパケット同士を比較し、各パケットの送信開始時刻、パケット長、及び信号強度が類似するパケット同士を1対1で紐付ける紐付け手段と、
前記最大値になる前記相関値が閾値を超えている場合に、紐付けられたパケットは前記第1及び第2の受信ノードが同一の送信源から受信した電波の受信信号系列であると判定した電波同定結果を出力する判定手段と、
を有することを特徴とする、電波環境可視化システム。
A first receiving node that outputs the first received signal sequence received, and
A second receiving node that is installed at a location different from the first receiving node, has an internal clock independent of the internal clock of the first receiving node, and outputs a second received signal sequence received. ,
A radio wave identification device that outputs a radio wave identification result based on the first received signal sequence and the second received signal sequence, and
Based on the radio wave identification result, a device for creating a heat map in which the received signal strength of the radio wave and the traffic level for each location of the first and second receiving nodes are expressed in color, and
With
The radio wave identification device is
The correlation value with the first received signal sequence is calculated while expanding and contracting the scale of the time axis of the second received signal sequence with the time expansion / contraction rate, and the second received signal sequence is offset only by the time on the time axis. A determination means for determining the time expansion / contraction rate and the time offset at which the correlation value is offset to the maximum value, and
By applying the time expansion / contraction rate and the time offset at which the correlation value becomes the maximum value to the time information of the second received signal sequence, the first received signal sequence and the second received signal sequence are respectively. A means of associating packets having similar transmission start times, packet lengths, and signal strengths with each other on a one-to-one basis by comparing each packet of each packet.
When the correlation value that becomes the maximum value exceeds the threshold value, it is determined that the associated packet is a reception signal sequence of radio waves received by the first and second receiving nodes from the same transmission source. Judgment means to output the radio wave identification result and
A radio wave environment visualization system characterized by having.
JP2019036507A 2019-02-28 2019-02-28 Radio wave identification device, radio wave identification method and radio wave environment visualization system Active JP7135938B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019036507A JP7135938B2 (en) 2019-02-28 2019-02-28 Radio wave identification device, radio wave identification method and radio wave environment visualization system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019036507A JP7135938B2 (en) 2019-02-28 2019-02-28 Radio wave identification device, radio wave identification method and radio wave environment visualization system

Publications (2)

Publication Number Publication Date
JP2020141316A true JP2020141316A (en) 2020-09-03
JP7135938B2 JP7135938B2 (en) 2022-09-13

Family

ID=72280738

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019036507A Active JP7135938B2 (en) 2019-02-28 2019-02-28 Radio wave identification device, radio wave identification method and radio wave environment visualization system

Country Status (1)

Country Link
JP (1) JP7135938B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022107403A1 (en) * 2020-11-18 2022-05-27 株式会社日立ソリューションズ Method, device, and program for determining difference between data sequences

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001201562A (en) * 2000-01-18 2001-07-27 Mitsubishi Electric Corp Radio detector
JP2008306240A (en) * 2007-06-05 2008-12-18 Kddi Corp Area map build up system for making area map based on communication quality information
JP2014041453A (en) * 2012-08-22 2014-03-06 Yokogawa Electric Corp Data similarity calculation method and data similarity calculation apparatus
JP2014231995A (en) * 2013-05-28 2014-12-11 有限会社アイ・アール・ティー Radar link device
JP2018029219A (en) * 2016-08-15 2018-02-22 富士通株式会社 Packet analysis program, packet analysis method and packet analyzer

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001201562A (en) * 2000-01-18 2001-07-27 Mitsubishi Electric Corp Radio detector
JP2008306240A (en) * 2007-06-05 2008-12-18 Kddi Corp Area map build up system for making area map based on communication quality information
JP2014041453A (en) * 2012-08-22 2014-03-06 Yokogawa Electric Corp Data similarity calculation method and data similarity calculation apparatus
JP2014231995A (en) * 2013-05-28 2014-12-11 有限会社アイ・アール・ティー Radar link device
JP2018029219A (en) * 2016-08-15 2018-02-22 富士通株式会社 Packet analysis program, packet analysis method and packet analyzer

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022107403A1 (en) * 2020-11-18 2022-05-27 株式会社日立ソリューションズ Method, device, and program for determining difference between data sequences
JP7438086B2 (en) 2020-11-18 2024-02-26 株式会社日立ソリューションズ Methods, devices and programs for determining differences between data series

Also Published As

Publication number Publication date
JP7135938B2 (en) 2022-09-13

Similar Documents

Publication Publication Date Title
CN110971295B (en) Optical fiber nonlinear noise monitoring device and optical receiver
JP6932494B2 (en) Systems and methods for applying aggregated cable test results data, as well as cloud-based computer servers
CN112072621B (en) Differential protection method for power distribution network line based on 5G communication and dynamic mode matching
CN110518990A (en) Calibration method, system and the computer readable storage medium of multiple antennas WiFi product
JP2020141316A (en) Radio wave identification device, radio wave identification method, and radio wave environment visualization system
JP3694281B2 (en) Multipath fading signal level calibration method, correction method, and mobile radio test system for sensitivity test of mobile radio
CN104518840B (en) A kind of method of wavelength calibration, optical network unit and optical line terminal
JP5532399B2 (en) Method and system for determining a dependency between device parameters and signal parameters of a mobile phone
JP3583730B2 (en) Wireless communication system and wireless transmission device
CN116527454A (en) Channel estimation method and device based on speech pilot frequency
CN108737964B (en) UWB indoor positioning system and positioning method thereof
US8412180B2 (en) Method for generating a signal train which can be played back on a tester for testing a mobile radio
JP5554800B2 (en) Fading simulator, mobile communication terminal test system, and fading simulation method
JP2019169877A (en) Monitoring system, monitoring method, and computer program
CN115913404A (en) Antenna calibration method and device and remote radio frequency unit
JP5539289B2 (en) Wireless communication device
JP7478302B1 (en) Wireless communication system, receiving device, and wireless communication method
CN101136678B (en) Location data processing method and device used for mobile station locating service
CN117292703B (en) Sound source positioning method and device for transformer equipment, electronic equipment and storage medium
CN116208247B (en) Signal processing method and system for underground optical fiber wireless transmission
WO2020203498A1 (en) Wireless communication characteristic evaluation method and wireless communication characteristic evaluation device
JP2007505532A (en) Synchronous monitoring method for transmitter of common wave network
CN115859481B (en) Simulation verification method and system for flight simulator
US20240036227A1 (en) Object detection system and object detection method
JP2000358273A (en) Base station unit and method for setting offset of pn code

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20211109

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220707

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: 20220802

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220815

R150 Certificate of patent or registration of utility model

Ref document number: 7135938

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150