(57)【特許請求の範囲】
【請求項１】
未知位置に配置された発信源から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサのそれぞれにおいて受信されており、前記複数のセンサのそれぞれから受信信号を取得する取得部と、
前記取得部において取得した受信信号に対してセンサ間における相互相関を導出するとともに、相互相関の絶対値の平均として予期される期待値を導出する第１導出部と、
前記第１導出部において導出した期待値から減算すべきレプリカであって、かつ未知位置の候補位置と、発信源の送信電力が１の電力スペクトルである規格化フィルタの候補と、伝搬学習プロセスで取得した伝送路特性とからレプリカを導出する第２導出部と、
前記第２導出部において導出したレプリカと前記第１導出部において導出した期待値との誤差を減算によって導出する第３導出部と、
前記第２導出部における候補位置を変更しながら、前記第２導出部と前記第３導出部とに対して処理を繰り返し実行させることによって導出された複数の誤差から、最小の誤差を選択する選択部と、
前記選択部において選択した誤差に対応した候補位置を、発信源が設置された位置として出力する出力部とを備え、
前記第２導出部は、規格化フィルタの候補から作成したレプリカ送信信号の自己相関と、候補位置および２つのセンサ間の伝送路特性の相互相関とを畳み込んでから、絶対値を計算した後に、平均として予期される期待値としてレプリカを導出することを特徴とする発信源推定装置。
【請求項２】
前記取得部は、センサに複数のアンテナが備えられている場合、アンテナごとの受信信号を取得し、
前記第１導出部は、異なったセンサに備えられたアンテナ間において、相互相関の絶対値の平均として予期される期待値を導出するとともに、同一のセンサに備えられたアンテナ間における受信信号の相互相関の期待値も導出し、
前記第２導出部は、異なったセンサに備えられたアンテナ間において、相互相関の絶対値の平均として予期される期待値を使用してレプリカを導出するとともに、同一のセンサに備えられたアンテナ間において、規格化フィルタの候補から作成したレプリカ送信信号の自己相関と、候補位置および２つのセンサ間の伝送路特性の相互相関とを畳み込んでから、絶対値を計算した後に、平均として予期される期待値としてレプリカを導出し、
前記第３導出部は、異なったセンサに備えられたアンテナ間に対する期待値とレプリカから誤差を導出するとともに、同一のセンサに備えられたアンテナ間における期待値とレプリカからも誤差を導出することを特徴とする請求項１に記載の発信源推定装置。
【請求項３】
前記選択部において選択した誤差に対応した候補位置を推定した後に当該候補位置における発信源の送信電力を推定する電力推定部をさらに備えることを特徴とする請求項１または２に記載の発信源推定装置。
【請求項４】
未知位置に配置された発信源から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサのそれぞれにおいて受信されており、前記複数のセンサのそれぞれから受信信号を取得するステップと、
取得した受信信号に対してセンサ間における相互相関を導出するとともに、相互相関の絶対値の平均として予期される期待値を導出するステップと、
導出した期待値から減算すべきレプリカであって、かつ未知位置の候補位置と、発信源の送信電力が１の電力スペクトルである規格化フィルタの候補と、伝搬学習プロセスで取得した伝送路特性とからレプリカを導出するステップと、
導出したレプリカと導出した期待値との誤差を減算によって導出するステップと、
候補位置を変更しながら、前記レプリカを導出するステップと前記誤差を導出するステップとに対して処理を繰り返し実行させることによって導出された複数の誤差から、最小の誤差を選択するステップと、
選択した誤差に対応した候補位置を、発信源が設置された位置として出力するステップとを備え、
前記レプリカを導出するステップは、規格化フィルタの候補から作成したレプリカ送信信号の自己相関と、候補位置および２つのセンサ間の伝送路特性の相互相関とを畳み込んでから、絶対値を計算した後に、平均として予期される期待値としてレプリカを導出することを特徴とする発信源推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、発信源推定技術に関し、未知の発信源を推定する発信源推定方法およびそれを利用した発信源推定装置に関する。
【背景技術】
【０００２】
近年の無線通信技術の急激な発展に伴って周波数資源が逼迫しているため、多数の無線システムが共存することができる電波の有効利用が求められている。そのため、一般に無線システムの利用にあたっては無線局免許が必要であるが、無線局免許を取得せず電波を利用する不法無線局が存在している場合があり、既存の無線システムに対して干渉を及ぼしてしまう。そのため日本の総務省ではＤｅｔｅｃｔ  Ｕｎｌｉｃｅｎｓｅｄ  Ｒａｄｉｏ  Ｓｔａｔｉｏｎｓ（ＤＥＵＲＡＳ）と呼ばれる電波監視システムを用いて不法無線局の探知を行っている。ＤＥＵＲＡＳでは複数のセンサを用いて不法無線局が送信している電波の電界強度および到来方向Ｄｉｒｅｃｔｉｏｎ  ｏｆ  Ａｒｒｉｖａｌ（ＤＯＡ）を測定することで不法無線局の位置を推定している（例えば、非特許文献１参照）。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】総務省、電波監視システム、［online］、インターネット＜URL：http://www.tele.soumu.go.jp/j/adm/monitoring/moni/type/deurasys/index.htm＞
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ＤＥＵＲＡＳ以外にも各センサの平均受信電力Ｒｅｃｅｉｖｅｄ  Ｓｉｇｎａｌ  Ｓｔｒｅｎｇｔｈ  Ｉｎｄｉｃａｔｏｒ（ＲＳＳＩ）を用いた方法やセンサ間の電波の到来遅延時間差Ｔｉｍｅ  Ｄｉｆｆｅｒｅｎｃｅ  ｏｆ  Ａｒｒｉｖａｌ（ＴＤＯＡ）、各センサへの電波の到来角方向ＤＯＡを用いた未知発信源の検出方法が提案されている。また、近年ではこれらの手法を複数センサを用いた手法への拡張が検討されている。現在では実際にのようなＲＳＳＩ、ＴＤＯＡ、ＤＯＡを利用した未知発信源の位置推定を実用化した製品も存在する。しかしながらこれらの手法は特に都市部においてマルチパスの影響から位置推定精度が劣化する。
【０００５】
本発明はこうした状況に鑑みてなされたものであり、その目的は、マルチパスの存在下においても推定精度の悪化を抑制する技術を提供することである。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、本発明のある態様の発信源推定装置は、未知位置に配置された発信源から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサのそれぞれにおいて受信されており、複数のセンサのそれぞれから受信信号を取得する取得部と、取得部において取得した受信信号に対してセンサ間における相互相関を導出するとともに、相互相関の絶対値の期待値を導出する第１導出部と、第１導出部において導出した期待値に対するレプリカであって、かつ未知位置が候補位置であると仮定した場合のレプリカを導出する第２導出部と、第２導出部において導出したレプリカと第１導出部において導出した期待値との誤差を導出する第３導出部と、第２導出部における候補位置を変更しながら、第２導出部と第３導出部とに対して処理を繰り返し実行させることによって導出された複数の誤差から、ひとつの誤差を選択する選択部と、選択部において選択した誤差に対応した候補位置を、発信源が設置された位置として出力する出力部とを備える。第２導出部は、候補位置から既知位置までの伝送路特性の相互相関であって、かつセンサ間における相互相関の絶対値の期待値と、発信源から送信された送信信号の自己相関とをもとにレプリカを導出する。
【０００７】
この態様によると、センサ間における相互相関の絶対値の期待値と、発信源から送信された送信信号の自己相関とをもとにレプリカを導出するので、ＲＳＳＩとＴＤＯＡを利用できる。
【０００８】
取得部は、センサに複数のアンテナが備えられている場合、アンテナごとの受信信号を取得し、第１導出部は、異なったセンサに備えられたアンテナ間において、相互相関の絶対値の期待値を導出するとともに、同一のセンサに備えられたアンテナ間における受信信号の相互相関の期待値も導出し、第２導出部は、異なったセンサに備えられたアンテナ間において、相互相関の絶対値の期待値を使用してレプリカを導出するとともに、同一のセンサに備えられたアンテナ間において、候補位置から既知位置までの伝送路特性の相互相関であって、かつアンテナ間における相互相関の期待値と、発信源から送信された送信信号の自己相関とをもとにレプリカを導出し、第３導出部は、異なったセンサに備えられたアンテナ間に対する期待値とレプリカから誤差を導出するとともに、同一のセンサに備えられたアンテナ間における期待値とレプリカからも誤差を導出してもよい。この場合、同一のセンサに備えられたアンテナ間における期待値とレプリカも導出するので、ＲＳＳＩとＴＤＯＡとに加えてＤＯＡも利用できる。
【０００９】
選択部において選択した誤差に対応した候補位置をもとに、発信源の送信電力を推定する電力推定部をさらに備えてもよい。この場合、候補位置を導出してから送信電力を推定するので、推定精度を向上できる。
【００１０】
本発明の別の態様は、発信源推定方法である。この方法は、未知位置に配置された発信源から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサのそれぞれにおいて受信されており、複数のセンサのそれぞれから受信信号を取得するステップと、取得した受信信号に対してセンサ間における相互相関を導出するとともに、相互相関の絶対値の期待値を導出するステップと、導出した期待値に対するレプリカであって、かつ未知位置が候補位置であると仮定した場合のレプリカを導出するステップと、導出したレプリカと導出した期待値との誤差を導出するステップと、候補位置を変更しながら、レプリカを導出するステップと誤差を導出するステップとに対して処理を繰り返し実行させることによって導出された複数の誤差から、ひとつの誤差を選択するステップと、選択した誤差に対応した候補位置を、発信源が設置された位置として出力するステップとを備える。レプリカを導出するステップは、候補位置から既知位置までの伝送路特性の相互相関であって、かつセンサ間における相互相関の絶対値の期待値と、発信源から送信された送信信号の自己相関とをもとにレプリカを導出する。
【００１１】
なお、以上の構成要素の任意の組合せ、本発明の表現を方法、装置、システム、記録媒体、コンピュータプログラムなどの間で変換したものもまた、本発明の態様として有効である。
【発明の効果】
【００１２】
本発明によれば、マルチパスの存在下においても推定精度の悪化を抑制できる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施例に係るＤＳＭシステムの構成を示す図である。
【図２】図１のＤＳＭシステムの処理手順を示す図である。
【図３】図１のＤＳＭシステムによる周波数利用状況マッピングを示す図である。
【図４】図１のＤＳＭシステムに対する数値解析の諸元を示す図である。
【図５】図１のＤＳＭシステムに対する数値解析におけるセンサ、未知発信源の配置を示す図である。
【図６】図１のＤＳＭシステムによる送信電力１０ｄＢｍでの位置推定特性評価結果を示す図である。
【図７】図１のＤＳＭシステムによる位置推定特性評価結果を示す図である。
【図８】図１のＤＳＭシステムによる伝搬路学習とＲＳＳＩを利用した送信電力推定結果を示す図である。
【図９】図９（ａ）−（ｂ）は、図１のＤＳＭシステムによるフィルタ推定結果を示す図である。
【図１０】図１のフュージョンセンタの構成を示す図である。
【発明を実施するための形態】
【００１４】
本発明を具体的に説明する前に、まず概要を述べる。本発明の実施例は、不法無線局や未知発信源の位置を推定する発信源推定装置に関する。発信源推定装置は、ＲＳＳＩ、ＴＤＯＡ、ＤＯＡの３つの情報をすべて活用する。さらに、本発明の実施例は、不法無線局や未知発信源の位置推定だけでなく、分散スペクトラムモニター（Ｄｉｓｔｒｉｂｕｔｅｄ  Ｓｐｅｃｔｒｕｍ  Ｍｏｎｉｔｏｒ，  ＤＳＭ）にも関する。分散スペクトラムモニターによって、既存システムも含めて、電波利用が希薄な周波数およびそのエリア（ホワイトスペース）の特定や、効率的にエナジーハーベスティングが可能な電波利用が盛んな周波数およびそのエリア（チャージングスペース）の特定も同時になされる。その結果、より高効率な周波数資源の活用が実現される
【００１５】
以下では、１．ＤＳＭの概要、２．ＤＳＭにおける複数センサを用いた未知発信源の推定手法、３．ＤＳＭを実現するための送信電力推定、４．特性評価、５．装置構成の順番で本実施例を説明する。なお、２．では、（ｉ）単一アンテナ、（ｉｉ）複数アンテナを用いた場合に分けて説明する。
【００１６】
１．ＤＳＭの概要
図１は、本発明の実施例に係るＤＳＭシステム１００の構成を示す。ＤＳＭシステム１００は、フュージョンセンタ１０、センサ１２と総称される第１センサ１２ａ、第２センサ１２ｂ、第３センサ１２ｃ、未知送信源１４と総称される第１未知送信源１４ａ、第２未知送信源１４ｂを含む。ここでは、３つのセンサ１２が示されているが、センサ１２の数はこれに限定されない。第１未知送信源１４ａあるいは第２未知送信源１４ｂから送信された信号は、複数のセンサ１２のそれぞれにおいて受信される。未知送信源１４が設置されている位置、信号の送信電力は未知である。各センサ１２は、受信信号をフュージョンセンタ１０に出力する。フュージョンセンタ１０は、各センサ１２からの受信電力をもとに、未知送信源１４が設置されている位置、信号の送信電力を推定する。推定の具体的な処理は後述する。
【００１７】
図２は、ＤＳＭシステム１００の処理手順を示す。ＤＳＭシステム１００は、前述のフュージョンセンタ１０、センサ１２に加えて既知情報データベース１６も含む。また、ＤＳＭシステム１００は、伝搬路学習プロセス、未知送信源推定プロセス、そしてスペクトラムマップ作成プロセスの３つからなる。以下では、各プロセスにおける具体的な処理を順を追って説明する。
【００１８】
まず、伝搬路学習プロセスでは、未知発信源推定に必要な、実際の環境における任意位置からセンサ１２間の伝搬路情報とその計算量の学習・取得を行う。しかしながら、任意位置すべてに既知発信源を設置し伝搬路推定を行うのは現実的に考えると非常に難しい。そこで地形データや建物データを用いてレイトレースシミュレーションによって、任意位置からセンサ１２までの伝搬路およびその統計量をモデリングする。しかしながら、実際の環境では人や車などの移動物や、事前データに存在しない障害物の存在などによってシミュレートした伝搬路とは異なる場合がある。そこでセルラシステムなどを既知発信源として活用して、既知位置からセンサ１２までの伝搬路を実際に測定し、その統計量を取得する。これを用いて先ほどのシミュレートした任意位置からの伝搬路を繰り返し補正し、実際の環境における任意位置とセンサ１２間の伝搬路情報およびその統計量を取得する。
【００１９】
次に、未知発信源推定プロセスでは、伝搬路学習プロセスで取得した伝搬路情報の統計量、未知送信源１４の送信信号パターン候補データベース、および未知送信源１４からの受信信号を用いて未知送信源１４の情報の推定を行う。各センサ１２もしくはアンテナ間における伝搬路情報の統計量から作成した相互相関と、送信信号の自己相関の畳み込み積分が、各センサ１２もしくはアンテナにおける受信信号の相互相関に等しくなるという特徴を用いて未知発信源を推定する。
【００２０】
最後に、スペクトラムマップ作成プロセスでは、推定した未知送信源１４の情報と、あらかじめ所有している既知発信源の情報を用いてスペクトラムマップの作成を行う。図３は、ＤＳＭシステム１００による周波数利用状況マッピングを示す。任意の位置や周波数における、受信電力を地形データ上にマッピングし、十分な電力を受けることのできるチャージングスペース２０２や利用されていない周波数領域ホワイトスペース２００を発見する。
【００２１】
２．ＤＳＭにおける複数センサを用いた未知発信源の推定手法
ここでは広帯域信号を想定し、伝搬路は伝搬遅延プロファイルが指数分布にしたがうマルチパスである統計モデルを用いる。また、各パスはレイリー分布にしたがい無相関である。なお、簡単のため未知送信源１４はひとつであるとし、最初のプロセスである伝搬路学習プロセスは学習済みであると仮定する。
【００２２】
（ｉ）単一アンテナの場合
ｉ番目の（ｉ＝１，・・・，Ｍ）センサ１２が任意の位置から受ける受信信号は、広帯域信号および伝搬環境が時々刻々と変化する場合、未知送信源１４の位置ψ０における動的伝搬チャネルと未知送信源１４の送信信号を用いて以下のように表される。
【数１】
0006032462.tif 000002
ここでｔとτｉは信号系列上の時間および遅延時間を表しており、
【数２】
0006032462.tif 000003
は伝搬路の変動を表す時間である。ｇ０は未知発信源の送信電力が１の電力スペクトル（以下、「規格化フィルタ」という）の形であり送信フィルタおよび変調方式に依存する。ｐ０は未知発信源の送信電力である。また、
【数３】
0006032462.tif 000004
は未知発信源の送信信号、
【数４】
0006032462.tif 000005
は未知発信源の位置ψ０からｉ番目のセンサ１２までの伝搬路応答、
【数５】
0006032462.tif 000006
は分散σ２の白色ガウス雑音である。
【００２３】
式（１）で表現されるｉ番目のセンサ１２の受信信号とｊ番目のセンサ１２の受信信号の時間ｔに関する相互相関は、Δｔをｉ番目のセンサ１２およびｊ番目のセンサ１２で観測される信号の伝搬遅延差分(ＴＤＯＡ)を表すものとして以下のように計算できる。ただし、本来ｉ=ｊは自己相関と呼ぶが以下では簡略化のために相互相関と表記する。また、式の簡潔化のためにψ０、ｇ０、ｐ０の記載を省略し、Ｒｘ［Ａ，Ｂ］を関数ＡとＢのｘに対する相関と置く。
【数６】
0006032462.tif 000007
【００２４】
ただし、Ｄｉｊはｉ=ｊのとき１、
【数７】
0006032462.tif 000008
のとき０である。
【数８】
0006032462.tif 000009
は統計的処理であるため
【数９】
0006032462.tif 000010
に依存しない。ここで、
【数１０】
0006032462.tif 000011
と置けば、式（２）はさらに以下のように変形できる。
【数１１】
0006032462.tif 000012
【００２５】
【数１２】
0006032462.tif 000013
はそれぞれ未知発信源の送信信号の時間に関する自己相関および未知発信源からｉ番目のセンサ１２までの伝搬路とｊ番目のセンサ１２までの伝搬路の相互相関を表している。すなわち各センサ１２における未知発信源からの受信信号同士の相互相関は、未知発信源から各センサ１２までの伝搬路同士の相互相関と、送信信号の自己相関との畳み込み積分で表すことができる。このとき送信信号の自己相関は、未知発信源の送信電力および規格化フィルタの形によって形状が決定する。また、伝搬路応答同士の相互相関は、
【数１３】
0006032462.tif 000014
のときにピーク値を持つことからｉ番目のセンサ１２までの伝搬遅延およびｊ番目のセンサ１２までの伝搬遅延の差分の情報を持つ関数となり、その形状はそれぞれのセンサまでの伝搬路の環境によって決定される。ここで、
【数１４】
0006032462.tif 000015
は
【数１５】
0006032462.tif 000016
に関して変化するため期待値を取る必要がある。
【００２６】
しかし、伝搬路応答の位相は一様分布するため、動的な環境で式（２）の期待値を取ると相互相関の値は０となり、位置に関する情報を抽出することができない。そのため、式（２）の絶対値を取ることで位相情報を棄却する。この処理によって
【数１６】
0006032462.tif 000017
に関して変化するのはｉ番目のセンサ１２とｊ番目のセンサ１２のそれぞれのチャネルの振幅に限られる。よって式（２）について絶対値を取り、
【数１７】
0006032462.tif 000018
に対する期待値を取ると以下のように表される。なお、Ｅｘ［・］をｘに対する期待値と置く。
【数１８】
0006032462.tif 000019
【００２７】
伝搬路応答の振幅は確率的に変動するため、式（４）からＲＳＳＩとＴＤＯＡ情報抽出が可能であることを示している。すなわち伝搬路の相互相関にはＴＤＯＡとＲＳＳＩの位置によって決定される２つの情報が含まれていることを示している。よって式（１）に示されるように受信信号は位置、送信電力および規格化フィルタにも依存するため、式（３）で表される相互相関および式（４）で示した絶対値の期待値を各センサの受信信号から作成し、さらに受信信号をセンサ中の最大平均受信電力によって規格化すると、以下の式で書き表すことができる。
【数１９】
0006032462.tif 000020
【００２８】
ここで、推定候補の場所から相互相関のレプリカを作成し、そのレプリカと
【数２０】
0006032462.tif 000021
の値の差分が最小となる場所を推定値とする。具体的に、位置候補ψ、規格化フィルタ候補ｇおよび伝搬学習プロセスで取得した伝搬路情報の統計量から作成したレプリカは式(３)から以下のように表される。
【数２１】
0006032462.tif 000022
【００２９】
【数２２】
0006032462.tif 000023
は推定規格化フィルタ候補ｇから作成したレプリカ送信信号
【数２３】
0006032462.tif 000024
の自己相関
【数２４】
0006032462.tif 000025
と、推定候補位置ψおよび伝搬路の学習プロセスで取得した動的伝搬路応答の統計量から作成したｉ番目のセンサ１２およびｊ番目のセンサ１２の伝搬路応答
【数２５】
0006032462.tif 000026
の相互相関
【数２６】
0006032462.tif 000027
を畳み込み、絶対値を取った後に
【数２７】
0006032462.tif 000028
に対する期待値を取ったものである。
【００３０】
よって、式（５）、式（８）を用いて、未知発信源の位置および規格化フィルタを以下の最尤推定法によって推定できる。
【数２８】
0006032462.tif 000029
式（９）において未知発信源の位置および規格化フィルタの２つのパラメータを変化させ、実際の受信信号から作成した
【数２９】
0006032462.tif 000030
と推定値のレプリカから作成した
【数３０】
0006032462.tif 000031
の差の二乗和が最小となったものがそれら２つのパラメータの推定値となる。また、真値同士の相互相関を引いた値は雑音となる。
【００３１】
（ｉｉ）複数アンテナの場合
ｉ番目のセンサ１２の第ｋアンテナ(ｋ=１，・・・，Ｎｉ)が任意の位置から受ける受信信号のレプリカは、広帯域信号、線形アレーアンテナおよび伝搬環境が時々刻々と動的に変化することを想定した場合、未知発信源の位置ψ０における動的伝搬路応答と未知発信源の送信信号を用いて以下のように表される。
【数３１】
0006032462.tif 000032
【００３２】
ここでｄｋは第１素子のアンテナを基準点としたときの、基準点から第ｋ素子の位置までの距離であり、λは到来波の波長を、θｉ（τｉ）は到来角度を表す。このとき同一センサ間の相互相関と異なるセンサにおけるアンテナ間の相互相関では期待値を取った場合の振る舞いが異なる。なお、以下では単一アンテナの場合と同様に式の簡潔のためにψ０、ｇ０、ｐ０の記載を省略する。
【００３３】
（ｉｉ−１）同一センサにおけるアンテナ間の相互相関
ｉ番目のセンサ１２に関して同一センサ間における第ｋアンテナと第ｌアンテナの相互相関は以下のように表される。
【数３２】
0006032462.tif 000033
ただし
【数３３】
0006032462.tif 000034
である。したがって、
【数３４】
0006032462.tif 000035
のときに相互相関のピークが現れることを示しており、
【数３５】
0006032462.tif 000036
に関する変化が相互相関に与える影響は、伝搬路応答の振幅と到来角度であるといえる。
【００３４】
式（１２）の相互相関において
【数３６】
0006032462.tif 000037
に関する期待値を取ると以下のように表される。
【数３７】
0006032462.tif 000038
【００３５】
よって、
【数３８】
0006032462.tif 000039
のときは値が０となり、
【数３９】
0006032462.tif 000040
のときは以下の式で表される。
【数４０】
0006032462.tif 000041
式（１４）から
【数４１】
0006032462.tif 000042
は同一センサの異なるアンテナ間における相互相関は動的な伝搬路応答において、その変化に対して期待値を取ったとしてもＲＳＳＩ、ＤＯＡの情報を抽出することが可能であることを示している。
【００３６】
（ｉｉ−２）異なるセンサにおけるアンテナ間の相互相関
ｉ番目のセンサ１２の第ｋアンテナとｊ番目のセンサ１２の第ｌアンテナとの相互相関は以下のように表される。
【数４２】
0006032462.tif 000043
【００３７】
よって、
【数４３】
0006032462.tif 000044
の変化が相互相関に与える影響は、伝搬路応答の振幅および位相と到来角度であるといえる。一般に動的伝搬路応答の位相は一様分布するため、式（１５）の期待値を取ると位置に関する情報を抽出することができない。そこで単一アンテナのときと同様に式（１５）について絶対値を取る。この処理によって
【数４４】
0006032462.tif 000045
に関して変化するのはｉ番目のセンサ１２およびｊ番目のセンサ１２のチャネルの振幅に限られる。よって式（１５）の絶対値について、
【数４５】
0006032462.tif 000046
に対する期待値を取ると以下のように表される。
【数４６】
0006032462.tif 000047
【００３８】
一般に、動的な伝搬路応答では振幅の確率分布は任意の形となるため、式（１６）はＲＳＳＩとＴＤＯＡ情報の抽出が可能であることを示している。以上から複数アンテナを用いた場合、伝搬路の相互相関にはＲＳＳＩ、ＴＤＯＡ、ＤＯＡの位置によって決定される３つの情報が含まれていることを示した。よって、式（１０）に示されるように受信信号は位置、送信電力および規格化フィルタにも依存するため、式（１２）、式（１５）で表される相互相関および式（１３）、式（１６）に示した絶対値の期待値を各センサの受信信号から作成し、さらに受信信号をセンサにおける最大平均受信電力によって規格化すると、以下の式に書き表すことができる。
【数４７】
0006032462.tif 000048
【００３９】
また、位置候補ψ、規格化フィルタ候補ｇおよび伝搬学習プロセスで取得した伝搬路情報の統計量から作成したレプリカは式（３）から以下のように表される。
【数４８】
0006032462.tif 000049
【数４９】
0006032462.tif 000050
は推定規格化フィルタ候補ｇから作成したレプリカ送信信号
【数５０】
0006032462.tif 000051
の自己相関
【数５１】
0006032462.tif 000052
と、推定候補位置ψおよび伝搬路の学習プロセスで取得した動的伝搬路応答の統計量から作成したｉ番目のセンサ１２の第ｋアンテナおよびｊ番目のセンサ１２の第ｌアンテナの伝搬路応答
【数５２】
0006032462.tif 000053
の相互相関を畳み込み、
【数５３】
0006032462.tif 000054
に対する期待値を取ったものである。
【００４０】
また、
【数５４】
0006032462.tif 000055
は畳み込みに対して絶対値を取った後に
【数５５】
0006032462.tif 000056
に対する期待値を取ったものである。よって、式（１７）、式（１８）、式（２１）を用いて、未知発信源の位置および規格化フィルタを以下の最尤推定法によって推定できる。
【数５６】
0006032462.tif 000057
式（２２）において未知発信源の位置および規格化フィルタの２つのパラメータを変化させ、実際の受信信号から作成した
【数５７】
0006032462.tif 000058
とレプリカから作成した
【数５８】
0006032462.tif 000059
との差の二乗和が最小となったものがそれら２つのパラメータの推定値となる。また、真値同士の相互相関を引いた値は雑音となる。
【００４１】
３．ＤＳＭを実現するための送信電力推定
式（９）および式（２２）を用いることで未知発信源の位置推定を行うことができるが、ＤＳＭのスペクトラムマップ作成のためには未知発信源の送信電力も知る必要がある。そこで、未知発信源の位置推定後に学習した伝搬路情報およびＲＳＳＩを用いることで、以下の最尤推定法から未知発信源の送信電力Ｐｔｘを推定する。
【数５９】
0006032462.tif 000060
【００４２】
ただし、
【数６０】
0006032462.tif 000061
であり、平均受信電力を表している。また、Ｐｔｘ，０は未知発信源の真の送信電力であり、
【数６１】
0006032462.tif 000062
は未知発信源の推定送信電力および学習した伝搬路情報から作成した推定平均受信電力である。式（２３）において未知発信源の送信電力をパラメータとして変化させ、実際の平均受信電力との差の和が最小となったものが送信電力パラメータの推定値となる。
【００４３】
４．特性評価
提案手法の未知発信源の位置、送信電力およびフィルタ推定精度を数値解析によって評価する。ここで、図４は、ＤＳＭシステム１００に対する数値解析の諸元を示し、図５は、ＤＳＭシステム１００に対する数値解析におけるセンサ１２、未知送信源１４の配置を示す。位置推定および送信電力推定では図４の諸元に加え未知発信源の送信電力を変化させ、各送信電力において推定を繰り返し行うことで評価を行った。フィルタ推定では図４の諸元に加え未知発信源の送信電力を−１０ｄＢｍとし、三角フィルタと矩形フィルタを用いて推定した。なおいずれの場合においても伝搬路は学習済みであるとしている。
【００４４】
図６は、ＤＳＭシステム１００による送信電力１０ｄＢｍでの位置推定特性評価結果を示す。これは、未知送信源１４の送信電力を１０ｄＢｍ、フィルタ既知、アンテナ２本を用いた場合の式（２２）の評価関数における誤差を示す。誤差は見やすくするためにログスケールで表記した。図６から、未知発信源の位置における評価関数の誤差が最小となっており、位置推定が正しく行われていることが確認できる。
【００４５】
図７は、ＤＳＭシステム１００による位置推定特性評価結果を示す。これは、未知送信源１４の位置推定について送信電力を−１０ｄＢｍから１５ｄＢｍまで５ｄＢｍごとに変化させた場合の、フィルタ既知の条件下において、ＲＳＳＩ、ＴＤＯＡそれぞれを用いた従来法と本実施例であるＤＳＭにおいてアンテナ１本とアンテナ２本を用いた場合の位置推定誤差の送信電力に対する標準偏差を示す。ここで、ＲＳＳＩ、ＴＤＯＡはそれぞれ本実施例と同様に伝搬路応答の相互相関を既知としている。いずれの送信電力においても本実施例はＲＳＳＩ、ＴＤＯＡを用いた従来法のどちらよりも推定誤差の標準偏差は小さく、推定精度がよいといえる。また、送信電力が−５ｄＢｍにおいても提案手法ではほぼ推定できている。
【００４６】
送信電力が−５ｄＢｍにおける各センサの平均受信ＳＮＲは図５の左のセンサから０．４、８．０、−２．６［ｄＢ］となっている。このことから受信電力が雑音電力より小さくなるような劣悪な受信環境においてもある程度の推定誤差を持って位置推定が可能であるといえる。さらに、アレーアンテナの素子数を２とした場合、アンテナ数が１の場合と大きな差はないものの推定誤差の標準偏差が小さくなっているといえる。したがって、さらに多くのアンテナをセンサに用いることでより精度の高い位置推定が可能であるといえる。
【００４７】
図８は、ＤＳＭシステム１００による伝搬路学習とＲＳＳＩを利用した送信電力推定結果を示す。これは、式（２３）を用いた未知送信源１４の送信電力推定について送信電力を−２０ｄＢｍから２０ｄＢｍまで５ｄＢｍごとに変化させた場合の、送信電力推定誤差の標準偏差を示す。ただし未知送信源１４の位置およびフィルタが既知であるとする。図８より、未知送信源１４の送信電力が０ｄＢｍより大きくなる場合は完全に推定できているが、これは送信電力が０ｄＢｍにおける各センサの平均受信ＳＮＲが５．４、１３．０、２．４［ｄＢ］となっており、各センサにおける平均受信ＳＮＲが０ｄＢを上回っているため推定できていると考えられる。未知発信源の送信電力が−５ｄＢｍのときに推定精度が劣化する理由としては、一部センサの平均受信ＳＮＲが０ｄＢを下回っているためであると考えられる。
【００４８】
図９（ａ）−（ｂ）は、ＤＳＭシステム１００によるフィルタ推定結果を示す。これは、送信電力を１０ｄＢｍ、アンテナ２本および未知送信源１４のフィルタを矩形フィルタおよび三角フィルタとした場合の式（２２）の評価関数における誤差をそれぞれ示す。未知送信源１４の位置における矩形フィルタを用いた場合と三角フィルタを用いた場合の誤差は矩形フィルタの誤差が最小となっているため、位置およびフィルタが正しく推定できているといえる。すなわち、未知発信源の位置とフィルタを同時に推定することが可能であるといえる。
【００４９】
５．装置構成
図１０は、フュージョンセンタ１０の構成を示す。フュージョンセンタ１０は、取得部２０、第１導出部２２、第２導出部２４、第３導出部２６、選択部２８、出力部３０、電力推定部３２を含む。まず、（ｉ）単一アンテナの場合について説明する。取得部２０は、未知位置に配置された未知送信源１４から送信された送信信号が、互いに異なった既知位置に配置された複数のセンサ１２のそれぞれにおいて受信されており、複数のセンサ１２のそれぞれから受信信号を取得する。受信信号は、式（１）のように示される。
【００５０】
第１導出部２２は、取得部２０において取得した受信信号に対してセンサ１２間における相互相関を導出するとともに、相互相関の絶対値の期待値を導出する。第１導出部２２での処理は、式（５）から式（７）に相当する。特に、相互相関の絶対値の期待値は、式（５）に相当する。
【００５１】
第２導出部２４は、第１導出部２２において導出した期待値に対するレプリカであって、かつ未知位置が候補位置であると仮定した場合のレプリカを導出する。具体的に説明すると、第２導出部２４は、候補位置から既知位置までの伝送路特性の相互相関であって、かつセンサ１２間における相互相関の絶対値の期待値と、未知送信源１４から送信された送信信号の自己相関とをもとにレプリカを導出する。レプリカは、式（８）に相当する。
【００５２】
第３導出部２６は、第２導出部２４において導出したレプリカと第１導出部２２において導出した期待値との誤差を導出する。選択部２８は、第２導出部２４における候補位置を変更しながら、第２導出部２４と第３導出部２６とに対して処理を繰り返し実行させることによって導出された複数の誤差から、ひとつの誤差を選択する。これら処理は、式（９）において、最小値を選択することに相当する。
【００５３】
出力部３０は、選択部２８において選択した誤差に対応した候補位置を、発信源が設置された位置として出力する。出力部３０から出力される位置は、式（９）において選択された２つのパラメータの推定値に相当する。電力推定部３２は、選択部２８において選択した誤差に対応した候補位置をもとに、発信源の送信電力を推定する。これは、式（２３）、式（２４）に相当する。
【００５４】
次に、（ｉｉ）複数アンテナの場合について説明する。取得部２０は、センサ１２に複数のアンテナが備えられている場合、アンテナごとの受信信号を取得する。受信信号は、式（１０）、式（１１）のように示される。第１導出部２２は、異なったセンサ１２に備えられたアンテナ間において、相互相関の絶対値の期待値を導出する。相互相関の絶対値の期待値は、式（１８）に相当する。また、第１導出部２２は、同一のセンサ１２に備えられたアンテナ間における受信信号の相互相関の期待値も導出する。相互相関の期待値は、式（１７）に相当する。
【００５５】
第２導出部２４は、異なったセンサ１２に備えられたアンテナ間において、相互相関の絶対値の期待値を使用してレプリカを導出する。レプリカは、式（２１）の下段に相当する。また、第２導出部２４は、同一のセンサ１２に備えられたアンテナ間において、候補位置から既知位置までの伝送路特性の相互相関であって、かつアンテナ間における相互相関の期待値と、未知送信源１４から送信された送信信号の自己相関とをもとにレプリカを導出する。レプリカは、式（２１）の上段に相当する。
【００５６】
第３導出部２６は、異なったセンサ１２に備えられたアンテナ間に対する期待値とレプリカから誤差を導出するとともに、同一のセンサ１２に備えられたアンテナ間における期待値とレプリカからも誤差を導出する。前者は、式（２２）の右辺の第２項に相当し、後者は、式（２２）の右辺の第１項に相当する。
【００５７】
この構成は、ハードウエア的には、任意のコンピュータのＣＰＵ、メモリ、その他のＬＳＩで実現でき、ソフトウエア的にはメモリにロードされたプログラムなどによって実現されるが、ここではそれらの連携によって実現される機能ブロックを描いている。したがって、これらの機能ブロックがハードウエアのみ、ソフトウエアのみ、またはそれらの組合せによっていろいろな形で実現できることは、当業者には理解されるところである。
【００５８】
本発明の実施例によれば、センサ間における相互相関の絶対値の期待値と、発信源から送信された送信信号の自己相関とをもとにレプリカを導出するので、ＲＳＳＩとＴＤＯＡを利用できる。また、ＲＳＳＩとＴＤＯＡが利用されるので、マルチパスの存在下においても推定精度の悪化を抑制できる。また、同一のセンサに備えられたアンテナ間における期待値とレプリカも導出するので、ＲＳＳＩとＴＤＯＡとに加えてＤＯＡも利用できる。また、ＲＳＳＩとＴＤＯＡとに加えてＤＯＡが利用されるので、マルチパスの存在下においても推定精度の悪化を抑制できる。また、未知送信源１４の位置を導出してから送信電力を推定するので、推定精度を向上できる。また、屋外における高精度な位置推定および波源の情報の推定が可能となり、違法電波発信源の発見できる。また、屋外における高精度な位置推定および波源の情報の推定が可能となり、ホワイトスペースやチャージングスペースを検出できる。
【００５９】
以上、本発明を実施例をもとに説明した。この実施例は例示であり、それらの各構成要素の組合せにいろいろな変形例が可能なこと、またそうした変形例も本発明の範囲にあることは当業者に理解されるところである。
【符号の説明】
【００６０】
１０  フュージョンセンタ、  １２  センサ、  １４  未知送信源、  １６  既知情報データベース、  ２０  取得部、  ２２  第１導出部、  ２４  第２導出部、  ２６  第３導出部、  ２８  選択部、  ３０  出力部、  ３２  電力推定部、  １００  ＤＳＭシステム。
【図１】
0006032462.tif 000063
【図２】
0006032462.tif 000064
【図３】
0006032462.tif 000065
【図４】
0006032462.tif 000066
【図５】
0006032462.tif 000067
【図７】
0006032462.tif 000068
【図８】
0006032462.tif 000069
【図１０】
0006032462.tif 000070
【図６】
0006032462.tif 000071
【図９】
0006032462.tif 000072