Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP6830300B2 - センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム - Google Patents
[go: Go Back, main page]

JP6830300B2 - センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム - Google Patents

センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム Download PDF

Info

Publication number
JP6830300B2
JP6830300B2 JP2016228118A JP2016228118A JP6830300B2 JP 6830300 B2 JP6830300 B2 JP 6830300B2 JP 2016228118 A JP2016228118 A JP 2016228118A JP 2016228118 A JP2016228118 A JP 2016228118A JP 6830300 B2 JP6830300 B2 JP 6830300B2
Authority
JP
Japan
Prior art keywords
wave
fmcw
sensor information
interrogation
bandwidth
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2016228118A
Other languages
English (en)
Other versions
JP2018084513A (ja
Inventor
柴田 孝基
孝基 柴田
陽敏 藤原
陽敏 藤原
高橋 英紀
英紀 高橋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Japan Radio Co Ltd
Original Assignee
Japan Radio Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Japan Radio Co Ltd filed Critical Japan Radio Co Ltd
Priority to JP2016228118A priority Critical patent/JP6830300B2/ja
Publication of JP2018084513A publication Critical patent/JP2018084513A/ja
Application granted granted Critical
Publication of JP6830300B2 publication Critical patent/JP6830300B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

本開示は、センサ情報を生成するセンサ情報生成装置と、センサ情報を取得するセンサ情報取得装置と、を備えるセンサ情報生成・取得システムに関する。
センサ情報を生成するセンサ情報生成装置と、センサ情報を取得するセンサ情報取得装置と、を備えるセンサ情報生成・取得システムに関して、特許文献1等の先行技術文献が存在しており、特許文献1の内容及び課題について説明する。
送受信装置は、センサ情報を取得するセンサ情報取得装置に対応し、アンテナを介して質問波を無線応答装置へ送信する。無線応答装置は、センサ情報を生成するセンサ情報生成装置に対応し、アンテナを介して質問波を送受信装置から受信する。
無線応答装置は、インピーダンス変換部及びSAW(Surface Acoustic Wave)−IDタグを備える。インピーダンス変換部は、アンテナとSAW−IDタグとの間に備えられ、外部からの刺激に応じてインピーダンスを変化させる。
外部からの刺激がないときには、インピーダンス変換部は、アンテナとSAW−IDタグとの間でインピーダンスを整合させないため、SAW−IDタグは、質問波を入力されず、応答波を生成できず、無線応答装置は、アンテナを介して応答波を送受信装置へ送信できない。外部からの刺激があるときには、インピーダンス変換部は、アンテナとSAW−IDタグとの間でインピーダンスを整合させるため、SAW−IDタグは、質問波を入力され、応答波を生成でき、無線応答装置は、アンテナを介して応答波を送受信装置へ送信できる。送受信装置は、応答波の有無に基づいて、外部からの刺激の有無を判定する。
特開2007−129657号公報
特許文献1に記載されていないが、送受信装置は、質問波と応答波との間の時間領域での相互相関を算出してもよい。すると、送受信装置は、時間領域での相互相関のピーク時間に基づいて、ピーク時間に等しいSAW遅延時間を有する無線応答装置を識別することができる。そして、送受信装置は、時間領域での相互相関のピーク時間での振幅に基づいて、無線応答装置における外部からの刺激の有無を判定することができる。
しかし、実際には、送受信アンテナ間で、質問波が回り込む。ここで、回り込み波は、短距離の空間伝搬による減衰を受けるのみであるため、電力レベルが高い。一方で、応答波は、長距離の空間伝搬及びSAW伝搬による減衰を受けるため、電力レベルが低い。なお、応答波について、長距離の空間伝搬の時間は、SAW伝搬の時間よりはるかに短い。
すると、送受信装置は、質問波と受信波(応答波及び回り込み波を含む。)との間の時間領域での相互相関を算出することになる。よって、送受信装置は、時間領域での相互相関のピーク時間が、ピーク時間に等しいSAW遅延時間を有する無線応答装置によるものであるのか、質問波の送受信アンテナ間での回り込みによるものであるのか、識別することができない。そして、送受信装置は、時間領域での相互相関のピーク時間での振幅が、無線応答装置における外部からの刺激の有無によるものであるのか、質問波の送受信アンテナ間での回り込みによるものであるのか、識別することができない。
そこで、前記課題を解決するために、本開示は、センサ情報取得装置において、質問波の送受信アンテナ間での回り込みがあるとしても、センサ情報生成装置から受信した応答波からセンサ情報生成装置が生成したセンサ情報を高精度で抽出することを目的とする。
上記目的を達成するために、レーダシステムで用いられるFMCW(Frequency Modulated Continuous Wave)方式を、センサ情報生成・取得システムで用いることとした。FMCW方式を用いるにあたり、信号長Tを長くすることで、及び/又は、帯域制限範囲内で信号帯域幅ΔFを広くすることで、パルス圧縮比TΔFを上げて、ノイズ耐性を高めることができる。帯域制限をさらに施すことにより、電力レベルを抑え過ぎなくても、無線チャネルマスク規格を満たすことができる。
ここで、FMCW質問波のスペクトルは、信号帯域幅において、ほぼ平坦であるといえる。よって、FMCW質問波の自己相関は、信号帯域幅の逆数の自然数倍の時間において、0になるはずである。そして、FMCW質問波とFMCW回り込み波との間の相互相関は、信号帯域幅の逆数の自然数倍の時間において、0になるはずである。なぜならば、FMCW回り込み波の空間伝搬時間は、FMCW応答波の空間伝搬時間+意図的に施された遅延時間(例えば、SAW遅延時間)と比べて、ほぼ0であるといえるからである。
そこで、センサ情報生成装置が施す遅延の時間を、信号帯域幅の逆数の自然数倍の時間に等しくする。すると、FMCW質問波とFMCW受信波(FMCW応答波及びFMCW回り込み波を含む。)との間の相互相関は、センサ情報生成装置が施す遅延の時間において、FMCW回り込み波からの寄与に影響されず、FMCW応答波からの寄与のみに影響されるはずである。よって、センサ情報取得装置において、FMCW質問波の送受信アンテナ間での回り込みがあるとしても、センサ情報生成装置から受信したFMCW応答波から、センサ情報生成装置が生成したセンサ情報を高精度で抽出することができるはずである。
しかし、実際には、FMCW質問波のスペクトルは、信号帯域中心周辺においては、平坦なレベルを維持しているが、信号帯域端においては、平坦なレベルから減衰していることがある。すると、FMCW質問波の自己相関は、全信号帯域幅の逆数の自然数倍の時間において、0にならない。そして、FMCW質問波とFMCW回り込み波との間の相互相関は、全信号帯域幅の逆数の自然数倍の時間において、0にならない。
よって、センサ情報生成装置が施す遅延の時間を、全信号帯域幅の逆数の自然数倍の時間に等しくしても、上述の効果を奏することができない。つまり、FMCW質問波とFMCW受信波(FMCW応答波及びFMCW回り込み波を含む。)との間の相互相関は、センサ情報生成装置が施す遅延の時間において、FMCW回り込み波からの寄与にも影響される。そして、センサ情報取得装置において、FMCW質問波の送受信アンテナ間での回り込みがあるときには、センサ情報生成装置から受信したFMCW応答波から、センサ情報生成装置が生成したセンサ情報を高精度で抽出することができない。
そこで、本開示は、FMCW質問波のスペクトルのうち、平坦なレベルを維持している信号帯域中心周辺のスペクトルのみを切り出して、時間領域での相関関数を算出することとした。すると、FMCW質問波の自己相関は、平坦な帯域幅の逆数の自然数倍の時間において、0になる。そして、FMCW質問波とFMCW回り込み波との間の相互相関は、平坦な帯域幅の逆数の自然数倍の時間において、0になる。
よって、センサ情報生成装置が施す遅延の時間を、平坦な帯域幅の逆数の自然数倍の時間に等しくすれば、上述の効果を奏することができる。つまり、FMCW質問波とFMCW受信波(FMCW応答波及びFMCW回り込み波を含む。)との間の相互相関は、センサ情報生成装置が施す遅延の時間において、FMCW回り込み波からの寄与には影響されない。そして、センサ情報取得装置において、FMCW質問波の送受信アンテナ間での回り込みがあるとしても、センサ情報生成装置から受信したFMCW応答波から、センサ情報生成装置が生成したセンサ情報を高精度で抽出することができる。
なお、平坦な帯域幅は、全信号帯域幅より狭い。よって、センサ情報生成装置が施す遅延の時間の分解能は、本開示では見かけ上悪化しているように思われる。しかし、FMCW回り込み波からの寄与の影響は、本開示ではほぼ見られない。よって、センサ情報生成装置が施す遅延の時間の分解能は、本開示では実質的に向上することができる。
具体的には、本開示は、センサ情報をセンサ情報生成装置から取得するセンサ情報取得装置であって、帯域制限されたFMCW質問波を、前記センサ情報生成装置に送信する質問波送信部と、前記センサ情報による変調と、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間だけの遅延と、が前記FMCW質問波に施されたFMCW応答波を、前記センサ情報生成装置から受信する応答波受信部と、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅において、前記FMCW質問波と、前記FMCW応答波及びFMCW回り込み波を含むFMCW受信波と、の間の時間領域での相互相関を算出する相互相関算出部と、前記時間領域での相互相関の振幅がピークをなす0以外の時間において、及び/又は、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間において、前記時間領域での相互相関の振幅及び/又は位相を算出することにより、前記FMCW受信波のうち前記FMCW回り込み波を除く前記FMCW応答波から、振幅変調及び/又は位相変調に用いられた前記センサ情報を抽出するセンサ情報抽出部と、を備えることを特徴とするセンサ情報取得装置である。
また、本開示は、センサ情報をセンサ情報取得装置に提供するセンサ情報生成装置であって、帯域制限されたFMCW質問波を、前記センサ情報取得装置から受信する質問波受信部と、前記センサ情報による変調を、前記FMCW質問波に施す質問波変調部と、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間だけの遅延を、前記FMCW質問波に施す質問波遅延部と、前記センサ情報による変調と、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間だけの遅延と、が前記FMCW質問波に施されたFMCW応答波を、前記センサ情報取得装置に送信する応答波送信部と、を備えることを特徴とするセンサ情報生成装置である。
また、本開示は、以上に記載のセンサ情報取得装置と、以上に記載のセンサ情報生成装置と、を備えることを特徴とするセンサ情報生成・取得システムである。
この構成によれば、センサ情報取得装置が、時間領域での相関関数を算出する帯域幅を適切に設定するとともに、センサ情報生成装置が、意図的に施す遅延時間(例えば、SAW遅延時間)を適切に設定することにより、システム全体として、上述の効果を奏することができる。つまり、センサ情報取得装置において、FMCW質問波の送受信アンテナ間での回り込みがあるとしても、センサ情報生成装置から受信したFMCW応答波から、センサ情報生成装置が生成したセンサ情報を高精度で抽出することができる。
このように、本開示によれば、センサ情報取得装置において、質問波の送受信アンテナ間での回り込みがあるとしても、センサ情報生成装置から受信した応答波からセンサ情報生成装置が生成したセンサ情報を高精度で抽出することができる。
本開示のセンサ情報生成・取得システムの概要を示す図である。 本開示のセンサ情報取得装置の構成を示す図である。 本開示のセンサ情報生成装置の構成を示す図である。 本開示の帯域制限前のFMCW質問波を示す図である。 本開示の帯域制限後のFMCW質問波を示す図である。 本開示のFMCW質問波のスペクトルを示す図である。 第1の比較例の時間相関の算出方法を示す図である。 第1の比較例のFMCW質問波の自己相関の算出結果を示す図である。 第2の比較例の時間相関の算出方法を示す図である。 第2の比較例のFMCW質問波の自己相関の算出結果を示す図である。 第1の本開示の時間相関の算出方法を示す図である。 第1の本開示のFMCW質問波の自己相関の算出結果を示す図である。 第2の本開示の時間相関の算出方法を示す図である。 第2の本開示のFMCW質問波の自己相関の算出結果を示す図である。 本開示のFMCW受信波のスペクトルを示す図である。 本開示のFMCW質問波とFMCW受信波との間の相互相関の算出結果を示す図である。
添付の図面を参照して本開示の実施形態を説明する。以下に説明する実施形態は本開示の実施の例であり、本開示は以下の実施形態に制限されるものではない。
(本開示のセンサ情報生成・取得システムの概要)
本開示のセンサ情報生成・取得システムの概要を図1に示す。センサ情報生成・取得システムSは、センサ情報取得装置1及びセンサ情報生成装置2−1、2−2、2−3から構成される。センサ情報生成・取得システムSは、例えば、橋梁Bの振動を検出する。センサ情報取得装置1は、例えば、橋梁Bの路肩に設置される。各センサ情報生成装置2−1、2−2、2−3は、例えば、橋梁Bの各車線の直下に設置される。
センサ情報取得装置1は、センサ情報生成装置2−1、2−2、2−3に質問波を送信する。各センサ情報生成装置2−1、2−2、2−3は、センサ情報取得装置1から質問波を受信し、橋梁Bの各車線の振動に応じた変調及び各所定時間だけの遅延(例えば、SAW遅延)を施された質問波を、センサ情報取得装置1に各応答波として送信する。センサ情報取得装置1は、各センサ情報生成装置2−1、2−2、2−3から各応答波を受信し、橋梁Bの各車線の振動に応じた変調及び各所定時間だけの遅延(例えば、SAW遅延)を施された各応答波から、橋梁Bの各車線の振動の情報を抽出する。橋梁Bの各車線の振動の情報は、直接に、又は、インターネットI等を介して、管理者に通知される。
このように、各センサ情報生成装置2−1、2−2、2−3が、各所定時間だけの遅延(例えば、SAW遅延)を質問波に施すことにより、センサ情報取得装置1は、各センサ情報生成装置2−1、2−2、2−3を識別することができる。そして、各センサ情報生成装置2−1、2−2、2−3が、橋梁Bの各車線の振動に応じた変調を質問波に施すことにより、センサ情報取得装置1は、橋梁Bの各車線の振動の情報を応答波から抽出することができる。さらに、各センサ情報生成装置2−1、2−2、2−3は、質問波及び橋梁Bの各車線の振動に応じて、パッシブに応答するのみであり、無電源で動作することができ、橋梁Bの各車線の直下での設置や保守を容易にすることができる。
本開示では、レーダシステムで用いられるFMCW方式を、センサ情報生成・取得システムSで用いることとした。FMCW方式を用いるにあたり、信号長Tを長くすることで、及び/又は、帯域制限範囲内で信号帯域幅ΔFを広くすることで、パルス圧縮比TΔFを上げて、ノイズ耐性を高めることができる。帯域制限をさらに施すことにより、電力レベルを抑え過ぎなくても、無線チャネルマスク規格を満たすことができる。
本開示のセンサ情報取得装置の構成を図2に示す。センサ情報取得装置1は、質問波送信部11、応答波受信部12、相互相関算出部13及びセンサ情報抽出部14から構成される。相互相関算出部13は、FFT算出部131、帯域中心切出部132、複素共役部133、FFT算出部134、帯域中心切出部135、スペクトル乗算部136及びIFFT算出部137から構成される。センサ情報抽出部14は、ピーク検出部141、振幅/位相算出部142及び振動等情報抽出部143から構成される。
本開示のセンサ情報生成装置の構成を図3に示す。N台のセンサ情報生成装置2−1、・・・、2−Nが設置される。センサ情報生成装置2−1は、質問波受信部21−1、質問波変調部22−1、質問波遅延部23−1及び応答波送信部24−1から構成される。・・・センサ情報生成装置2−Nは、質問波受信部21−N、質問波変調部22−N、質問波遅延部23−N及び応答波送信部24−Nから構成される。
質問波送信部11は、帯域制限されたFMCW質問波を、センサ情報生成装置2−1、・・・、2−Nに送信する。質問波受信部21−1、・・・、21−Nは、帯域制限されたFMCW質問波を、センサ情報取得装置1から受信する。
各質問波変調部22−1、・・・、22−Nは、各センサ情報による変調を、FMCW質問波に施す。各センサ情報として、例えば、振動、圧力、温度及び湿度等の情報(より具体的には、橋梁Bの各車線の振動の情報)が挙げられる。各質問波遅延部23−1、・・・、23−Nは、各所定時間T、・・・、Tだけの遅延を、FMCW質問波に施す。各質問波遅延部として、例えば、電磁波を音波等に変換し、音波等を伝搬させ、音波等を電磁波に再変換する、音波等伝搬回路(より具体的には、SAW遅延回路)が挙げられる。
各応答波送信部24−1、・・・、24−Nは、各センサ情報による変調と、各所定時間T、・・・、Tだけの遅延と、がFMCW質問波に施された各FMCW応答波を、センサ情報取得装置1に送信する。応答波受信部12は、各センサ情報による変調と、各所定時間T、・・・、Tだけの遅延と、がFMCW質問波に施された各FMCW応答波を、各応答波送信部24−1、・・・、24−Nから受信する。
実際には、質問波送信部11と応答波受信部12と間で、FMCW質問波が回り込む。ここで、FMCW回り込み波は、短距離の空間伝搬による減衰を受けるのみであるため、電力レベルが高い。一方で、各FMCW応答波は、長距離の空間伝搬及び各質問波遅延部23−1、・・・、23−Nによる減衰を受けるため、電力レベルが低い。なお、各FMCW応答波について、図8、10、12、14、16で示すように、長距離の空間伝搬の時間は、各質問波遅延部23−1、・・・、23−Nの伝搬の時間よりはるかに短い。
相互相関算出部13は、FMCW質問波と、各FMCW応答波及びFMCW回り込み波を含むFMCW受信波と、の間の時間領域での相互相関を算出する。センサ情報抽出部14は、FMCW受信波のうちFMCW回り込み波を除く各FMCW応答波から、各変調に用いられた各センサ情報を抽出する。相互相関算出部13及びセンサ情報抽出部14については、以下の段落で比較例と本開示とを比較しながら詳述する。
以下の段落では、等価低域系を用いて説明する。本開示の帯域制限前のFMCW質問波を図4に示す。質問波送信部11は、1msecにわたり−500kHzから+500kHzまで、振幅を一定にしながら周波数を掃引することにより、図4に示したFMCW質問波を生成する。本開示の帯域制限後のFMCW質問波を図5に示す。質問波送信部11は、図4に示したFMCW質問波に対して、信号帯域外の電力を減らす窓関数を時間領域で乗算することにより、図5に示したFMCW質問波を生成する。信号帯域外の電力を減らす窓関数として、Cosine−Tukey窓関数、Hamming窓関数、Hanning窓関数及びBlackmann窓関数等が挙げられる。本開示のFMCW質問波のスペクトルを図6に示す。図6の縦軸は、FMCW質問波のスペクトルの電力レベルを対数表示dBで示す。FMCW質問波のスペクトルの電力レベルは、後述する数1、5に記載のU(k)に比例する。帯域制限後の信号帯域幅は、1MHzである。
図6に示したように、FMCW質問波のスペクトルは、信号帯域幅1MHzにおいて、信号帯域端を除いて、ほぼ平坦であるといえる。よって、FMCW質問波の自己相関は、信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecにおいて、0になるはずである。そして、FMCW質問波とFMCW回り込み波との間の相互相関は、信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecにおいて、0になるはずである。なぜならば、FMCW回り込み波の空間伝搬時間は、各FMCW応答波の各空間伝搬時間+意図的に施された各遅延時間(例えば、各SAW遅延時間)と比べて、ほぼ0であるといえるからである。
そこで、各センサ情報生成装置2−1、・・・、2−Nが施す各遅延の時間T、・・・、Tを、信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecに等しくする。すると、FMCW質問波とFMCW受信波(各FMCW応答波及びFMCW回り込み波を含む。)との間の相互相関は、各センサ情報生成装置2−1、・・・、2−Nが施す各遅延の時間T、・・・、Tにおいて、FMCW回り込み波からの寄与に影響されず、各FMCW応答波からの寄与のみに影響されるはずである。よって、センサ情報取得装置1において、FMCW質問波の質問波送信部11と応答波受信部12との間での回り込みがあるとしても、各センサ情報生成装置2−1、・・・、2−Nから受信した各FMCW応答波から、各センサ情報生成装置2−1、・・・、2−Nが生成した各センサ情報を高精度で抽出することができるはずである。しかし、実際には上述の効果を奏することができないことを、まず比較例として説明する。そして、上述の効果を奏することができるための方法を、次に本開示として説明する。
(比較例の時間相関の算出方法及びFMCW質問波の自己相関の算出結果)
比較例では、FMCW質問波のスペクトルのうち、帯域中心周辺で平坦である帯域幅とともに、信号帯域端で減衰する帯域幅において、時間領域での相関関数を算出する。
FFT算出部131は、数1に示すように、時間領域でのFMCW質問波u(n)を、周波数領域でのFMCW質問波U(k)にフーリエ変換する。ここで、nは、サンプリング周期を用いて正規化された時間を表わす。帯域中心切出部132は、図7、9に示すように、周波数領域でのFMCW質問波U(k)から、帯域中心周辺で平坦である帯域幅とともに、信号帯域端で減衰する帯域幅において、周波数領域でのFMCW質問波U(k)を切り出す。複素共役部133は、周波数領域でのFMCW質問波U(k)を、周波数領域でのFMCW質問波U (k)に複素共役化する。
Figure 0006830300
時間領域でのFMCW受信波x(n)は、時間領域での各FMCW応答波v(n)とFMCW回り込み波w(n)とを重畳したものである。FFT算出部134は、数2に示すように、時間領域でのFMCW受信波x(n)=v(n)+w(n)を、周波数領域でのFMCW受信波X(k)=V(k)+W(k)にフーリエ変換する。帯域中心切出部135は、周波数領域でのFMCW受信波X(k)=V(k)+W(k)から、帯域中心切出部132が周波数領域でのFMCW質問波U(k)を切り出す帯域幅と同様の帯域幅において、周波数領域でのFMCW受信波X(k)=V(k)+W(k)を切り出す。
Figure 0006830300
スペクトル乗算部136は、数3に示すように、周波数領域でのFMCW質問波U (k)と、周波数領域でのFMCW受信波X(k)=V(k)+W(k)と、を乗算して、周波数領域での相互相関Y(k)を算出する。IFFT算出部137は、数4に示すように、周波数領域での相互相関Y(k)を、時間領域での相互相関y(n)に逆フーリエ変換する。時間領域での相互相関y(n)は、yB,v(n)及びyB,w(n)を含む。yB,v(n)は、FMCW質問波と各FMCW応答波との間の時間領域での相互相関を表わす。yB,w(n)は、FMCW質問波とFMCW回り込み波との間の時間領域での相互相関を表わし、FMCW質問波の時間領域での自己相関に比例する。
Figure 0006830300
Figure 0006830300
ここで、FMCW質問波と各FMCW応答波との間の時間領域での相互相関yB,v(n)を算出して、各FMCW応答波の遅延時間及びセンサ情報を抽出したい。すると、FMCW質問波とFMCW回り込み波との間の時間領域での相互相関yB,w(n)を低減する必要がある。つまり、FMCW質問波の時間領域での自己相関を低減する必要がある。
第1の比較例の時間相関の算出方法を図7に示す。図7の縦軸は、FMCW質問波のスペクトルの電力レベルを対数表示dBで示す。FMCW質問波のスペクトルの電力レベルは、上述した数1に記載のU(k)に比例する。時間相関の算出に用いる帯域幅は、帯域中心周辺で平坦である帯域幅とともに、信号帯域端で減衰する帯域幅を含めて、1MHzである。サンプリング周波数は、4MHzである。
第1の比較例のFMCW質問波の自己相関の算出結果を図8に示す。図8の縦軸は、第1の比較例のFMCW質問波の自己相関の電力レベルを対数表示dBで示す。第1の比較例のFMCW質問波の自己相関の電力レベルは、上述した数4に記載のyB,w(n)に比例する。このように、FMCW質問波の自己相関は、全信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecにおいて、0にならない。つまり、FMCW質問波とFMCW回り込み波との間の相互相関yB,w(n)は、全信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecにおいて、0にならない。
第2の比較例の時間相関の算出方法を図9に示す。図9の縦軸は、FMCW質問波のスペクトルの電力レベルを対数表示dBで示す。FMCW質問波のスペクトルの電力レベルは、上述した数1に記載のU(k)に比例する。時間相関の算出に用いる帯域幅は、帯域中心周辺で平坦である帯域幅とともに、信号帯域端で減衰する帯域幅を含めて、1MHzである。サンプリング周波数も、1MHzである。
第2の比較例のFMCW質問波の自己相関の算出結果を図10に示す。図10の縦軸は、第2の比較例のFMCW質問波の自己相関の電力レベルを対数表示dBで示す。第2の比較例のFMCW質問波の自己相関の電力レベルは、上述した数4に記載のyB,w(n)に比例する。このように、FMCW質問波の自己相関は、全信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecにおいて、0にならない。つまり、FMCW質問波とFMCW回り込み波との間の相互相関yB,w(n)は、全信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecにおいて、0にならない。
よって、各センサ情報生成装置2−1、・・・、2−Nが施す各遅延の時間T、・・・、Tを、全信号帯域幅1MHzの逆数1μsecの自然数1、・・・、N倍の時間1、・・・、Nμsecに等しくしても、上述の効果を奏することができない。つまり、FMCW質問波とFMCW受信波との間の相互相関y(n)は、各センサ情報生成装置2−1、・・・、2−Nが施す各遅延の時間1、・・・、Nμsecにおいて、FMCW回り込み波からの寄与にも影響される。そして、センサ情報取得装置1において、FMCW質問波の質問波送信部11と応答波受信部12との間での回り込みがあるときには、各センサ情報生成装置2−1、・・・、2−Nから受信した各FMCW応答波から、各センサ情報生成装置2−1、・・・、2−Nが生成した各センサ情報を高精度で抽出することができない。
(本開示の時間相関の算出方法及びFMCW質問波の自己相関の算出結果)
本開示では、FMCW質問波のスペクトルのうち、信号帯域端で減衰する帯域幅を除いて、帯域中心周辺で平坦である帯域幅において、時間領域での相関関数を算出する。
FFT算出部131は、数5に示すように、時間領域でのFMCW質問波u(n)を、周波数領域でのFMCW質問波U(k)にフーリエ変換する。ここで、nは、サンプリング周期を用いて正規化された時間を表わす。帯域中心切出部132は、図11、13に示すように、周波数領域でのFMCW質問波U(k)から、信号帯域端で減衰する帯域幅を除いて、帯域中心周辺で平坦である帯域幅において、周波数領域でのFMCW質問波U(k)を切り出す。複素共役部133は、周波数領域でのFMCW質問波U(k)を、周波数領域でのFMCW質問波U (k)に複素共役化する。
Figure 0006830300
時間領域でのFMCW受信波x(n)は、時間領域での各FMCW応答波v(n)とFMCW回り込み波w(n)とを重畳したものである。FFT算出部134は、数6に示すように、時間領域でのFMCW受信波x(n)=v(n)+w(n)を、周波数領域でのFMCW受信波X(k)=V(k)+W(k)にフーリエ変換する。帯域中心切出部135は、周波数領域でのFMCW受信波X(k)=V(k)+W(k)から、帯域中心切出部132が周波数領域でのFMCW質問波U(k)を切り出す帯域幅と同様の帯域幅において、周波数領域でのFMCW受信波X(k)=V(k)+W(k)を切り出す。
Figure 0006830300
スペクトル乗算部136は、数7に示すように、周波数領域でのFMCW質問波U (k)と、周波数領域でのFMCW受信波X(k)=V(k)+W(k)と、を乗算して、周波数領域での相互相関Y(k)を算出する。IFFT算出部137は、数8に示すように、周波数領域での相互相関Y(k)を、時間領域での相互相関y(n)に逆フーリエ変換する。時間領域での相互相関y(n)は、yC,v(n)及びyC,w(n)を含む。yC,v(n)は、FMCW質問波と各FMCW応答波との間の時間領域での相互相関を表わす。yC,w(n)は、FMCW質問波とFMCW回り込み波との間の時間領域での相互相関を表わし、FMCW質問波の時間領域での自己相関に比例する。
Figure 0006830300
Figure 0006830300
ここで、FMCW質問波と各FMCW応答波との間の時間領域での相互相関yC,v(n)を算出して、各FMCW応答波の遅延時間及びセンサ情報を抽出したい。すると、FMCW質問波とFMCW回り込み波との間の時間領域での相互相関yC,w(n)を低減する必要がある。つまり、FMCW質問波の時間領域での自己相関を低減する必要がある。
第1の本開示の時間相関の算出方法を図11に示す。図11の縦軸は、FMCW質問波のスペクトルの電力レベルを対数表示dBで示す。FMCW質問波のスペクトルの電力レベルは、上述した数5に記載のU(k)に比例する。時間相関の算出に用いる帯域幅は、信号帯域端で減衰する帯域幅を除いて、帯域中心周辺で平坦である帯域幅を含めて、0.8MHzである。サンプリング周波数は、4MHzである。
第1の本開示のFMCW質問波の自己相関の算出結果を図12に示す。図12の縦軸は、第1の本開示のFMCW質問波の自己相関の電力レベルを対数表示dBで示す。第1の本開示のFMCW質問波の自己相関の電力レベルは、上述した数8に記載のyC,w(n)に比例する。このように、FMCW質問波の自己相関は、平坦な帯域幅0.8MHzの逆数1.25μsecの自然数1、2、・・・、N倍の時間1.25、2.50、・・・、1.25Nμsecにおいて、ほぼ0になる。つまり、FMCW質問波とFMCW回り込み波との間の相互相関yC,w(n)は、平坦な帯域幅0.8MHzの逆数1.25μsecの自然数1、2、・・・、N倍の時間1.25、2.50、・・・、1.25Nμsecにおいて、ほぼ0になる。
第2の本開示の時間相関の算出方法を図13に示す。図13の縦軸は、FMCW質問波のスペクトルの電力レベルを対数表示dBで示す。FMCW質問波のスペクトルの電力レベルは、上述した数5に記載のU(k)に比例する。時間相関の算出に用いる帯域幅は、信号帯域端で減衰する帯域幅を除いて、帯域中心周辺で平坦である帯域幅を含めて、0.8MHzである。サンプリング周波数も、0.8MHzである。
第2の本開示のFMCW質問波の自己相関の算出結果を図14に示す。図14の縦軸は、第2の本開示のFMCW質問波の自己相関の電力レベルを対数表示dBで示す。第2の本開示のFMCW質問波の自己相関の電力レベルは、上述した数8に記載のyC,w(n)に比例する。このように、FMCW質問波の自己相関は、平坦な帯域幅0.8MHzの逆数1.25μsecの自然数1、2、・・・、N倍の時間1.25、2.50、・・・、1.25Nμsecにおいて、ほぼ0になる。つまり、FMCW質問波とFMCW回り込み波との間の相互相関yC,w(n)は、平坦な帯域幅0.8MHzの逆数1.25μsecの自然数1、2、・・・、N倍の時間1.25、2.50、・・・、1.25Nμsecにおいて、ほぼ0になる。
よって、各センサ情報生成装置2−1、2−2、・・・、2−Nが施す各遅延の時間T、T、・・・、Tを、平坦な帯域幅0.8MHzの逆数1.25μsecの自然数1、2、・・・、N倍の時間1.25、2.50、・・・、1.25Nμsecに等しくすれば、上述の効果を奏することができる。つまり、FMCW質問波とFMCW受信波との間の相互相関y(n)は、各センサ情報生成装置2−1、2−2、・・・、2−Nが施す各遅延の時間1.25、2.50、・・・、1.25Nμsecにおいて、FMCW回り込み波からの寄与には影響されない。そして、センサ情報取得装置1において、FMCW質問波の質問波送信部11と応答波受信部12との間での回り込みがあるとしても、各センサ情報生成装置2−1、・・・、2−Nから受信した各FMCW応答波から、各センサ情報生成装置2−1、・・・、2−Nが生成した各センサ情報を高精度で抽出することができる。
なお、平坦な帯域幅0.8MHzは、全信号帯域幅1MHzより狭い。よって、各センサ情報生成装置2−1、・・・、2−Nが施す各遅延の時間T、・・・、Tの分解能は、本開示では見かけ上悪化しているように思われる(比較例での1μsec→本開示での1.25μsec)。しかし、FMCW回り込み波からの寄与の影響は、本開示ではほぼ見られない。よって、各センサ情報生成装置2−1、・・・、2−Nが施す各遅延の時間T、・・・、Tの分解能は、本開示では実質的に向上することができる。
(本開示のFMCW質問波とFMCW受信波との間の相互相関の算出結果)
本開示のFMCW受信波のスペクトルを図15に示す。図15の縦軸は、FMCW受信波のスペクトルの電力レベルを対数表示dBで示す。FMCW受信波のスペクトルの電力レベルは、上述した数6に記載のX(k)に比例する。時間領域でのFMCW受信波x(n)は、時間領域での各FMCW応答波v(n)とFMCW回り込み波w(n)とを重畳したものである。時間領域での各FMCW応答波v(n)として、センサ情報生成装置2−3が遅延時間T=3.75μsecを施したFMCW応答波のみを取り上げる。
本開示のFMCW質問波とFMCW受信波との間の相互相関の算出結果を図16に示す。図16の縦軸は、本開示のFMCW質問波とFMCW受信波との間の相互相関の電力レベルを対数表示dBで示す。本開示のFMCW質問波とFMCW受信波との間の相互相関の電力レベルは、上述した数8に記載のy(n)=yC,v(n)+yC,w(n)に比例する。時間相関の算出に用いる帯域幅は、FMCW質問波における信号帯域端で減衰する帯域幅を除いて、FMCW質問波における帯域中心周辺で平坦である帯域幅を含めて、0.8MHzである。サンプリング周波数も、0.8MHzである。
ピーク検出部141は、FMCW質問波とFMCW受信波との間の相互相関y(n)の振幅|y(n)|のピーク位置を検出する。一般的には、|y(n)|のピーク位置は、時間T=T、・・・Tに対応する正規化時間n=n、・・・nである。図16では、|y(n)|のピーク位置は、時間T=Tに対応する正規化時間n=nである。ただし、|y(n)|のピーク位置として、図12に示したような、FMCW質問波の自己相関の振幅のピーク位置を採用するべきではなく、つまり、時間TN−1<T<Tに対応する正規化時間nN−1<n<nを採用するべきではない。
次に、振幅/位相算出部142は、FMCW質問波とFMCW受信波との間の相互相関y(n)の振幅|y(n)|がピークをなす0以外の正規化時間n=n、・・・nにおいて、FMCW質問波とFMCW受信波との間の相互相関y(n)の振幅|y(n)|及び/又は位相θ(n)を算出する。又は、振幅/位相算出部142は、FMCW質問波のスペクトルU(k)が帯域中心周辺で平坦である帯域幅の逆数の自然数倍の正規化時間n=n、・・・nにおいて、FMCW質問波とFMCW受信波との間の相互相関y(n)の振幅|y(n)|及び/又は位相θ(n)を算出する。
振動等情報抽出部143は、FMCW質問波とFMCW受信波との間の相互相関y(n)の振幅|y(n)|及び/又は位相θ(n)(n=n、・・・n)に基づいて、FMCW受信波のうちFMCW回り込み波を除く各FMCW応答波から、各振幅変調及び/又は各位相変調に用いられた各センサ情報s(n)(n=n、・・・n)を抽出する。
そして、振動等情報抽出部143は、各々単数回抽出した各センサ情報s(n)(n=n、・・・n)を、各センサ情報生成装置2−1、・・・、2−Nが設置される橋梁Bの各車線の直下における加速度についての情報として抽出する。さらに、振動等情報抽出部143は、各々複数回抽出した各センサ情報s(n)(n=n、・・・n)に対するFFT結果を、各センサ情報生成装置2−1、・・・、2−Nが設置される橋梁Bの各車線の直下における振動についての情報として抽出する。各センサ情報生成装置2−1、・・・、2−Nが設置される橋梁Bの各車線の直下における加速度及び振動についての情報は、直接に、又は、インターネットI等を介して、管理者に通知される。
ここで、センサ情報取得装置1は、各FMCW応答波を振幅復調するにあたり、各FMCW応答波の振幅が、各センサ情報生成装置2−1、・・・、2−Nとの間の空間伝搬で減衰されたものであるのか、各センサ情報生成装置2−1、・・・、2−Nでの各センサ情報で生成されたものであるのか、区別することが望ましい。よって、各センサ情報生成装置2−1、・・・、2−Nは、図3で示した構成を設置するのみならず、各センサ情報と無関係なリファレンス素子を設置することが望ましい。そして、センサ情報取得装置1は、図2で示した構成からの各FMCW応答波と、各センサ情報と無関係なリファレンス素子からのリファレンス信号と、の差分を計測することが望ましい。
本開示のセンサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システムは、橋梁等の構造物の検査等、センサ情報を応答器から質問器へと無線通信する用途に適用することができる。
S:センサ情報生成・取得システム
B:橋梁
I:インターネット
1:センサ情報取得装置
2−1、2−2、2−3、2−N:センサ情報生成装置
11:質問波送信部
12:応答波受信部
13:相互相関算出部
14:センサ情報抽出部
21−1、21−N:質問波受信部
22−1、22−N:質問波変調部
23−1、23−N:質問波遅延部
24−1、24−N:応答波送信部
131:FFT算出部
132:帯域中心切出部
133:複素共役部
134:FFT算出部
135:帯域中心切出部
136:スペクトル乗算部
137:IFFT算出部
141:ピーク検出部
142:振幅/位相算出部
143:振動等情報抽出部

Claims (3)

  1. センサ情報をセンサ情報生成装置から取得するセンサ情報取得装置であって、
    帯域制限されたFMCW(Frequency Modulated Continuous Wave)質問波を、前記センサ情報生成装置に送信する質問波送信部と、
    前記センサ情報による変調と、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間だけの遅延と、が前記FMCW質問波に施されたFMCW応答波を、前記センサ情報生成装置から受信する応答波受信部と、
    前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅において、前記FMCW質問波と、前記FMCW応答波及びFMCW回り込み波を含むFMCW受信波と、の間の時間領域での相互相関を算出する相互相関算出部と、
    前記時間領域での相互相関の振幅がピークをなす0以外の時間において、及び/又は、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間において、前記時間領域での相互相関の振幅及び/又は位相を算出することにより、前記FMCW受信波のうち前記FMCW回り込み波を除く前記FMCW応答波から、振幅変調及び/又は位相変調に用いられた前記センサ情報を抽出するセンサ情報抽出部と、
    を備えることを特徴とするセンサ情報取得装置。
  2. センサ情報をセンサ情報取得装置に提供するセンサ情報生成装置であって、
    帯域制限されたFMCW(Frequency Modulated Continuous Wave)質問波を、前記センサ情報取得装置から受信する質問波受信部と、
    前記センサ情報による変調を、前記FMCW質問波に施す質問波変調部と、
    前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅における、前記FMCW質問波と、FMCW応答波及びFMCW回り込み波を含むFMCW受信波と、の間の時間領域での相互相関が、前記センサ情報生成装置が前記FMCW質問波に施す遅延の時間において、前記FMCW回り込み波からの寄与には影響されないように、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間だけの遅延を、前記FMCW質問波に施す質問波遅延部と、
    前記センサ情報による変調と、前記FMCW質問波のスペクトルが帯域中心周辺で平坦である帯域幅の逆数の自然数倍の時間だけの遅延と、が前記FMCW質問波に施された前記FMCW応答波を、前記センサ情報取得装置に送信する応答波送信部と、
    を備えることを特徴とするセンサ情報生成装置。
  3. 請求項1に記載のセンサ情報取得装置と、
    請求項2に記載のセンサ情報生成装置と、
    を備えることを特徴とするセンサ情報生成・取得システム。
JP2016228118A 2016-11-24 2016-11-24 センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム Active JP6830300B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016228118A JP6830300B2 (ja) 2016-11-24 2016-11-24 センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016228118A JP6830300B2 (ja) 2016-11-24 2016-11-24 センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム

Publications (2)

Publication Number Publication Date
JP2018084513A JP2018084513A (ja) 2018-05-31
JP6830300B2 true JP6830300B2 (ja) 2021-02-17

Family

ID=62236700

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016228118A Active JP6830300B2 (ja) 2016-11-24 2016-11-24 センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム

Country Status (1)

Country Link
JP (1) JP6830300B2 (ja)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09147284A (ja) * 1995-11-07 1997-06-06 Siemens Ag 表面音響波により作動する無線問合わせ装置
US7034660B2 (en) * 1999-02-26 2006-04-25 Sri International Sensor devices for structural health monitoring
JP2013162293A (ja) * 2012-02-03 2013-08-19 Mitsubishi Electric Corp 無線通信装置、無線送信装置および無線受信装置
EP3709730A1 (en) * 2013-01-16 2020-09-16 Interdigital Patent Holdings, Inc. Discovery signal generation and reception

Also Published As

Publication number Publication date
JP2018084513A (ja) 2018-05-31

Similar Documents

Publication Publication Date Title
DE602004028950D1 (de) Vorrichtung und verfahren zur mehrantennen-senderstrahlformung von breitbandigen paketierten drahtlosen kommunikationssignalen mit hoher datenrate
Han et al. Experimental demonstration of underwater acoustic communication using bionic signals
CN110178434A (zh) 传输探测参考信号的方法、终端设备和网络设备
CN103308892B (zh) 一种多扫频雷达信号的产生及压缩方法
CN107076834A (zh) 具有增大的多普勒能力的雷达操作
JP2017515362A5 (ja)
US10746863B2 (en) Target extraction system, target extraction method, information processing apparatus, and control method and control program of information processing apparatus
US10225126B2 (en) Synchronization method, receive end, and transmit end
JP2015522977A5 (ja)
RU2014149209A (ru) Устройство и способ применения специальных спектральных масок для передачи в суб-гигагерцовых диапазонах
Ivanov et al. Software-defined radio technology in the problem concerning with the successive sounding of HF ionospheric communication channels
CN109863773A (zh) 传输上行信号的方法、终端设备和网络侧设备
JP5249115B2 (ja) パルス圧縮装置および物標探知装置
CN101159482A (zh) 一种移动通信系统及其信号传输方法
CN109644347A (zh) 传输控制信道的方法、网络设备和终端设备
GB2589233A (en) System and method for analyzing an oilfield casing using an apparent thickness approach
JP5229251B2 (ja) 目標検出方法、パッシブレーダ装置及びレーダ装置
JP6830300B2 (ja) センサ情報取得装置、センサ情報生成装置及びセンサ情報生成・取得システム
CN106771845B (zh) 一种确定第一反射点距离射频单元的距离的方法和装置
JP3920194B2 (ja) 超音波測定装置
CN113176541B (zh) 一种自适应抗频谱弥散干扰方法及系统
JP2016180679A (ja) 時間差測位システム
CN110463080A (zh) 一种导频处理方法及装置、系统
JP6029287B2 (ja) 距離測定方法及びレーダー装置
JP6451382B2 (ja) 目標物検出装置、目標物検出方法、目標物検出プログラムおよび記憶媒体

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191120

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20201030

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20201104

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201224

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210122

R150 Certificate of patent or registration of utility model

Ref document number: 6830300

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150