JP6032462B2 - 発信源推定方法およびそれを利用した発信源推定装置 - Google Patents
発信源推定方法およびそれを利用した発信源推定装置 Download PDFInfo
- Publication number
- JP6032462B2 JP6032462B2 JP2012078280A JP2012078280A JP6032462B2 JP 6032462 B2 JP6032462 B2 JP 6032462B2 JP 2012078280 A JP2012078280 A JP 2012078280A JP 2012078280 A JP2012078280 A JP 2012078280A JP 6032462 B2 JP6032462 B2 JP 6032462B2
- Authority
- JP
- Japan
- Prior art keywords
- replica
- deriving
- correlation
- transmission
- expected value
- 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
Links
- 238000000034 method Methods 0.000 title claims description 45
- 230000005540 biological transmission Effects 0.000 claims description 152
- 230000008569 process Effects 0.000 claims description 24
- 238000009795 derivation Methods 0.000 claims description 17
- 238000010606 normalization Methods 0.000 claims description 13
- 238000001228 spectrum Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 description 15
- 238000011156 evaluation Methods 0.000 description 10
- 230000004927 fusion Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000006866 deterioration Effects 0.000 description 4
- 238000007476 Maximum Likelihood Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 239000000470 constituent Substances 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Description
図1は、本発明の実施例に係るDSMシステム100の構成を示す。DSMシステム100は、フュージョンセンタ10、センサ12と総称される第1センサ12a、第2センサ12b、第3センサ12c、未知送信源14と総称される第1未知送信源14a、第2未知送信源14bを含む。ここでは、3つのセンサ12が示されているが、センサ12の数はこれに限定されない。第1未知送信源14aあるいは第2未知送信源14bから送信された信号は、複数のセンサ12のそれぞれにおいて受信される。未知送信源14が設置されている位置、信号の送信電力は未知である。各センサ12は、受信信号をフュージョンセンタ10に出力する。フュージョンセンタ10は、各センサ12からの受信電力をもとに、未知送信源14が設置されている位置、信号の送信電力を推定する。推定の具体的な処理は後述する。
ここでは広帯域信号を想定し、伝搬路は伝搬遅延プロファイルが指数分布にしたがうマルチパスである統計モデルを用いる。また、各パスはレイリー分布にしたがい無相関である。なお、簡単のため未知送信源14はひとつであるとし、最初のプロセスである伝搬路学習プロセスは学習済みであると仮定する。
i番目の(i=1,・・・,M)センサ12が任意の位置から受ける受信信号は、広帯域信号および伝搬環境が時々刻々と変化する場合、未知送信源14の位置ψ0における動的伝搬チャネルと未知送信源14の送信信号を用いて以下のように表される。
i番目のセンサ12の第kアンテナ(k=1,・・・,Ni)が任意の位置から受ける受信信号のレプリカは、広帯域信号、線形アレーアンテナおよび伝搬環境が時々刻々と動的に変化することを想定した場合、未知発信源の位置ψ0における動的伝搬路応答と未知発信源の送信信号を用いて以下のように表される。
i番目のセンサ12に関して同一センサ間における第kアンテナと第lアンテナの相互相関は以下のように表される。
i番目のセンサ12の第kアンテナとj番目のセンサ12の第lアンテナとの相互相関は以下のように表される。
式(9)および式(22)を用いることで未知発信源の位置推定を行うことができるが、DSMのスペクトラムマップ作成のためには未知発信源の送信電力も知る必要がある。そこで、未知発信源の位置推定後に学習した伝搬路情報およびRSSIを用いることで、以下の最尤推定法から未知発信源の送信電力Ptxを推定する。
提案手法の未知発信源の位置、送信電力およびフィルタ推定精度を数値解析によって評価する。ここで、図4は、DSMシステム100に対する数値解析の諸元を示し、図5は、DSMシステム100に対する数値解析におけるセンサ12、未知送信源14の配置を示す。位置推定および送信電力推定では図4の諸元に加え未知発信源の送信電力を変化させ、各送信電力において推定を繰り返し行うことで評価を行った。フィルタ推定では図4の諸元に加え未知発信源の送信電力を−10dBmとし、三角フィルタと矩形フィルタを用いて推定した。なおいずれの場合においても伝搬路は学習済みであるとしている。
図10は、フュージョンセンタ10の構成を示す。フュージョンセンタ10は、取得部20、第1導出部22、第2導出部24、第3導出部26、選択部28、出力部30、電力推定部32を含む。まず、(i)単一アンテナの場合について説明する。取得部20は、未知位置に配置された未知送信源14から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサ12のそれぞれにおいて受信されており、複数のセンサ12のそれぞれから受信信号を取得する。受信信号は、式(1)のように示される。
Claims (4)
- 未知位置に配置された発信源から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサのそれぞれにおいて受信されており、前記複数のセンサのそれぞれから受信信号を取得する取得部と、
前記取得部において取得した受信信号に対してセンサ間における相互相関を導出するとともに、相互相関の絶対値の平均として予期される期待値を導出する第1導出部と、
前記第1導出部において導出した期待値から減算すべきレプリカであって、かつ未知位置の候補位置と、発信源の送信電力が1の電力スペクトルである規格化フィルタの候補と、伝搬学習プロセスで取得した伝送路特性とからレプリカを導出する第2導出部と、
前記第2導出部において導出したレプリカと前記第1導出部において導出した期待値との誤差を減算によって導出する第3導出部と、
前記第2導出部における候補位置を変更しながら、前記第2導出部と前記第3導出部とに対して処理を繰り返し実行させることによって導出された複数の誤差から、最小の誤差を選択する選択部と、
前記選択部において選択した誤差に対応した候補位置を、発信源が設置された位置として出力する出力部とを備え、
前記第2導出部は、規格化フィルタの候補から作成したレプリカ送信信号の自己相関と、候補位置および2つのセンサ間の伝送路特性の相互相関とを畳み込んでから、絶対値を計算した後に、平均として予期される期待値としてレプリカを導出することを特徴とする発信源推定装置。 - 前記取得部は、センサに複数のアンテナが備えられている場合、アンテナごとの受信信号を取得し、
前記第1導出部は、異なったセンサに備えられたアンテナ間において、相互相関の絶対値の平均として予期される期待値を導出するとともに、同一のセンサに備えられたアンテナ間における受信信号の相互相関の期待値も導出し、
前記第2導出部は、異なったセンサに備えられたアンテナ間において、相互相関の絶対値の平均として予期される期待値を使用してレプリカを導出するとともに、同一のセンサに備えられたアンテナ間において、規格化フィルタの候補から作成したレプリカ送信信号の自己相関と、候補位置および2つのセンサ間の伝送路特性の相互相関とを畳み込んでから、絶対値を計算した後に、平均として予期される期待値としてレプリカを導出し、
前記第3導出部は、異なったセンサに備えられたアンテナ間に対する期待値とレプリカから誤差を導出するとともに、同一のセンサに備えられたアンテナ間における期待値とレプリカからも誤差を導出することを特徴とする請求項1に記載の発信源推定装置。 - 前記選択部において選択した誤差に対応した候補位置を推定した後に当該候補位置における発信源の送信電力を推定する電力推定部をさらに備えることを特徴とする請求項1または2に記載の発信源推定装置。
- 未知位置に配置された発信源から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサのそれぞれにおいて受信されており、前記複数のセンサのそれぞれから受信信号を取得するステップと、
取得した受信信号に対してセンサ間における相互相関を導出するとともに、相互相関の絶対値の平均として予期される期待値を導出するステップと、
導出した期待値から減算すべきレプリカであって、かつ未知位置の候補位置と、発信源の送信電力が1の電力スペクトルである規格化フィルタの候補と、伝搬学習プロセスで取得した伝送路特性とからレプリカを導出するステップと、
導出したレプリカと導出した期待値との誤差を減算によって導出するステップと、
候補位置を変更しながら、前記レプリカを導出するステップと前記誤差を導出するステップとに対して処理を繰り返し実行させることによって導出された複数の誤差から、最小の誤差を選択するステップと、
選択した誤差に対応した候補位置を、発信源が設置された位置として出力するステップとを備え、
前記レプリカを導出するステップは、規格化フィルタの候補から作成したレプリカ送信信号の自己相関と、候補位置および2つのセンサ間の伝送路特性の相互相関とを畳み込んでから、絶対値を計算した後に、平均として予期される期待値としてレプリカを導出することを特徴とする発信源推定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012078280A JP6032462B2 (ja) | 2012-03-29 | 2012-03-29 | 発信源推定方法およびそれを利用した発信源推定装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012078280A JP6032462B2 (ja) | 2012-03-29 | 2012-03-29 | 発信源推定方法およびそれを利用した発信源推定装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013205398A JP2013205398A (ja) | 2013-10-07 |
JP6032462B2 true JP6032462B2 (ja) | 2016-11-30 |
Family
ID=49524598
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012078280A Active JP6032462B2 (ja) | 2012-03-29 | 2012-03-29 | 発信源推定方法およびそれを利用した発信源推定装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6032462B2 (ja) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6331072B2 (ja) * | 2014-03-03 | 2018-05-30 | 株式会社国際電気通信基礎技術研究所 | ホワイトスペース検出装置、ホワイトスペース検出方法、及びプログラム |
JP6663098B2 (ja) * | 2015-03-12 | 2020-03-11 | 株式会社国際電気通信基礎技術研究所 | 波源位置選択装置、波源位置算出装置、波源位置選択方法、波源位置算出方法、及びプログラム |
JP6488808B2 (ja) * | 2015-03-24 | 2019-03-27 | 三菱電機株式会社 | 時間差測位システム |
JP6819797B2 (ja) | 2017-11-30 | 2021-01-27 | 日本電気株式会社 | 位置推定装置、位置推定方法とプログラム、並びに、位置推定システム |
JP6659054B2 (ja) * | 2018-05-17 | 2020-03-04 | Necプラットフォームズ株式会社 | 電波監視装置、電波監視システム、電波監視方法及びプログラム |
JP7180133B2 (ja) * | 2018-06-14 | 2022-11-30 | 株式会社Soken | 端末位置推定システム、端末位置推定装置、端末位置推定方法、及び制御プログラム |
JP7315349B2 (ja) * | 2019-03-25 | 2023-07-26 | 株式会社国際電気通信基礎技術研究所 | 位置推定装置、及び位置推定方法 |
CN109932677B (zh) * | 2019-04-08 | 2021-03-23 | 南京航空航天大学 | 基于无人机和云gis***的黑广播定位方法 |
WO2020208987A1 (ja) * | 2019-04-12 | 2020-10-15 | 日本電気株式会社 | 位置推定装置、位置推定方法及びプログラム |
WO2021048907A1 (ja) * | 2019-09-09 | 2021-03-18 | 日本電気株式会社 | 電波発信源位置推定システム |
JP7007023B2 (ja) * | 2020-01-23 | 2022-01-24 | Necプラットフォームズ株式会社 | 電波監視装置、電波監視システム、電波監視方法及びプログラム |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05288823A (ja) * | 1992-04-07 | 1993-11-05 | Kansai Tec:Kk | 微小放射電磁波の検出方法及びその放射源の位置検出方法並びにこれらを利用した微小放射電磁波検出装置 |
JP2917145B1 (ja) * | 1998-06-11 | 1999-07-12 | 宇宙開発事業団 | 複数衛星による電磁波発生源観測方法 |
JP4115379B2 (ja) * | 2003-08-28 | 2008-07-09 | 株式会社国際電気通信基礎技術研究所 | 送受信システム |
US7869810B2 (en) * | 2005-04-25 | 2011-01-11 | Agilent Technologies, Inc. | Method and system for computing and displaying location information from cross-correlation data |
JP2009175096A (ja) * | 2008-01-28 | 2009-08-06 | Mitsubishi Electric Corp | 信号源位置推定装置 |
JP2011179946A (ja) * | 2010-03-01 | 2011-09-15 | Doshisha | 位置推定方法および位置推定システム |
-
2012
- 2012-03-29 JP JP2012078280A patent/JP6032462B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2013205398A (ja) | 2013-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6032462B2 (ja) | 発信源推定方法およびそれを利用した発信源推定装置 | |
JP6032469B2 (ja) | 発信源推定方法およびそれを利用した発信源推定装置 | |
US11002825B2 (en) | Position determination system having a deconvolution decoder using a joint snr-time of arrival approach | |
JP6251930B2 (ja) | 位置推定システム | |
WO2019159965A1 (ja) | 電波環境推定方法及び電波環境推定装置 | |
WO2015177559A1 (en) | Method for localising scattering elements in a 3d environment | |
Chaudhari et al. | Spatial interpolation of cyclostationary test statistics in cognitive radio networks: Methods and field measurements | |
JP6396037B2 (ja) | データ解析装置及び方法 | |
KR101627419B1 (ko) | 이동 노드의 위치 추정 방법 및 그 장치 | |
Traßl et al. | Deriving an empirical channel model for wireless industrial indoor communications | |
KR101219913B1 (ko) | 이동 노드 위치 추정 방법 | |
JP5583170B2 (ja) | 散乱***置推定装置、散乱***置推定方法及びプログラム | |
Silva et al. | Towards non-line-of-sight ranging error mitigation in industrial wireless sensor networks | |
De Sousa et al. | Mobile station localization emitter in urban NLoS using multipath ray tracing fingerprints and machine learning | |
JP6331072B2 (ja) | ホワイトスペース検出装置、ホワイトスペース検出方法、及びプログラム | |
JP6747688B2 (ja) | 発信源推定方法およびそれを利用した発信源推定装置 | |
Haniz et al. | Localization of illegal radios utilizing cross-correlation of channel impulse response with interpolation in urban scenarios | |
US20220276333A1 (en) | Radio wave source position estimation system | |
Feng et al. | Sensor calibration using nonparametric statistical characterization of error models | |
Betta et al. | Experimental analysis of the frequency diversity to improve localization in WSNs | |
Venieris et al. | Preprocessing algorithm for source localisation in a multipath environment | |
Sen et al. | A node deployment mechanism accounting into received signal strength and frequency diversity for a wireless sensor network | |
KhanbashI et al. | Measurements and analysis of fingerprinting structures for wlan localization systems | |
Wang et al. | Convex hull based node selection NLoS mitigation for indoor localization | |
Chang et al. | Integrated cross-correlation and genetic algorithm to improve TDOA-based locating system in preventing radio interference |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150326 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20150326 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20160120 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160202 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160401 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160406 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160421 |
|
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: 20160913 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20161011 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161012 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20161011 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6032462 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |