JP2019148577A - 波源位置推定装置、コンピュータに実行させるためのプログラム、およびプログラムを記録したコンピュータ読み取り可能な記録媒体 - Google Patents

Abstract

【課題】計算量を抑制して正確に波源位置を推定可能な波源位置推定装置を提供する。【解決手段】波源位置推定装置１は、受信手段１２と、推定手段１４とを備える。推定手段１４は、受信手段１２によって受信された複数の実端末装置の複数の実端末情報に含まれる複数の位置情報と複数の受信電力とに基づいて受信電力によって重み付けされた複数の実端末装置の重心を波源位置として推定し、その推定した波源位置に誤差があると判定したとき、少なくとも１つの仮想端末装置の位置を示す位置情報と受信電力とを追加して複数の実端末装置および少なくとも１つの仮想端末装置の位置情報および受信電力に基づいて受信電力によって重み付けされた複数の実端末装置および少なくとも１つの仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する。【選択図】図２

Description

この発明は、波源位置推定装置、コンピュータに実行させるためのプログラム、およびプログラムを記録したコンピュータ読み取り可能な記録媒体に関するものである。

周波数の２次利用は、１次利用者への影響が無い範囲で１次利用者が利用する周波数を利用可能である。そこで、周波数の利用効率を高めるために、周波数の２次利用の検討が進められている。

周波数の２次利用の可否を判定するために、周波数の１次利用の有無を検出する必要がある。そして、周波数の１次利用の有無を検出するためには、波源位置を推定する必要がある。

従来、波源位置を推定する方法として、非特許文献１に記載の方法が知られている。非特許文献１に記載の方法では、ＴＤｏＡ（Ｔｉｍｅ Ｄｉｆｆｅｒｅｎｔｉａｌ ｏｆ Ａｒｒｉｖａｌ）を用いて波源位置を推定する。

より具体的には、長距離センサーと呼ばれる高精度なセンサーを用いて波源からの電波を受信し、複数のセンサーへの電波の到来時間差から波源を推定する。

即ち、電波が波源から２つのセンサーＳ_Ａ，Ｓ_Ｂへ到来するときの到来時間差｜ｔ_ａ−ｔ_ｂ｜を算出する。そして、到来時間差｜ｔ_ａ−ｔ_ｂ｜を用いて、電波が波源からセンサーＳ_Ａ，Ｓ_Ｂへ到来するときの伝搬距離差｜Ｌ_ａ−Ｌ_ｂ｜を算出する。

同様にして、電波が波源からセンサーＳ_Ｂ，Ｓ_Ｃへ到来するときの伝搬距離差｜Ｌ_ｂ−Ｌ_ｃ｜と、電波が波源からセンサーＳ_Ｃ，Ｓ_Ａへ到来するときの伝搬距離差｜Ｌ_ｃ−Ｌ_ａ｜とを算出する。

そうすると、３つの双曲線｜Ｌ_ａ−Ｌ_ｂ｜，｜Ｌ_ｂ−Ｌ_ｃ｜，｜Ｌ_ｃ−Ｌ_ａ｜を描画し、３つの双曲線の交点を波源位置として推定する。

堀端 研志、菅野 一生、古川 玲、長谷川 晃朗、武内 良男，"異種センサーによる局所的ホワイトスペース推定に関する要素技術の一検討，"２０１４年 電子情報通信学会総合大会，通信講演論文集１，Ｂ−１７−２８，ｐ６０７．

しかし、非特許文献１に記載の推定方法では、センサー間の同期を取る必要があり、高機能なシステムが必要となる問題がある。また、波源位置推定のために計算量が大きくなるという問題がある。

モニターするエリアおよび周波数範囲が広くなると、上記の問題は、より顕著になる。

そこで、この発明の実施の形態によれば、計算量を抑制して正確に波源位置を推定可能な波源位置推定装置を提供する。

また、この発明の実施の形態によれば、計算量を抑制して波源位置の正確な推定をコンピュータに実行させるためのプログラムを提供する。

更に、この発明の実施の形態によれば、計算量を抑制して波源位置の正確な推定をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体を提供する。

（構成１）
この発明の実施の形態によれば、波源位置推定装置は、受信手段と、推定手段とを備える。受信手段は、移動端末または静止端末からなる実端末装置の位置情報と実端末装置における受信電力とを含む実端末情報を複数の実端末装置から受信する。推定手段は、複数の実端末装置の複数の実端末情報に含まれる複数の位置情報と複数の受信電力とに基づいて受信電力によって重み付けされた複数の実端末装置の重心を波源位置として推定し、その推定した波源位置に誤差があると判定したとき、少なくとも１つの仮想端末装置の位置を示す位置情報と受信電力とを追加して複数の実端末装置および少なくとも１つの仮想端末装置の位置情報および受信電力に基づいて受信電力によって重み付けされた複数の実端末装置および少なくとも１つの仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する。

（構成２）
構成１において、推定手段は、電力分布の推定結果を示す電力分布推定結果に基づいて少なくとも１つの仮想端末装置の位置情報と受信電力とを追加し、その追加した少なくとも１つの仮想端末装置の位置情報と受信電力と、複数の実端末装置の複数の位置情報と複数の受信電力とに基づいて推定処理を収束条件が満たされるまで繰り返し実行する。

（構成３）
構成２において、推定手段は、少なくとも１つの仮想端末装置をランダムに配置し、その配置した位置を示す位置情報を少なくとも１つの仮想端末装置の位置情報とし、電力分布推定結果に基づいて、少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を少なくとも１つの仮想端末装置の受信電力とする。

（構成４）
前構成２において、推定手段は、電力分布推定結果によって示される等電力線上に少なくとも１つの仮想端末装置を配置し、等電力線上の電力値を少なくとも１つの仮想端末装置の受信電力とし、少なくとも１つの仮想端末装置を配置した位置を示す位置情報を少なくとも１つの仮想端末装置の位置情報とする。

（構成５）
構成２において、推定手段は、複数の実端末装置のうち、隣接する実端末装置間の距離を求め、その求めた距離が第１のしきい値以下である実端末装置が配置されている領域以外の領域に少なくとも１つの仮想端末装置を配置し、その配置した位置を示す位置情報を少なくとも１つの仮想端末装置の位置情報とし、電力分布推定結果に基づいて、少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を少なくとも１つの仮想端末装置の受信電力とする。

（構成６）
構成２から構成５のいずれかにおいて、推定手段は、電力分布推定結果によって示される等電力線上の電力値と実端末装置における受信電力との電力差が第２のしきい値よりも大きいか否かを複数の実端末装置の全てについて判定し、電力差が第２のしきい値よりも大きい実端末装置を除外して推定処理を収束条件が満たされるまで繰り返し実行する。

（構成７）
構成１から構成６のいずれかにおいて、推定手段は、複数の実端末装置および少なくとも１つの仮想端末装置の位置の平均が第３のしきい値以下であるとき、収束条件を満たすと判定し、推定処理を完了する。

（構成８）
構成１から構成６のいずれかにおいて、推定手段は、推定処理の実行回数が所定回数に達したとき、収束条件を満たすと判定し、推定処理を完了する。

（構成９）
構成１から構成６のいずれかにおいて、推定手段は、推定処理を２回以上実行したとき、推定処理をｋ−１（ｋは２以上の整数）回実行したときに推定された波源位置と、推定処理をｋ回実行したときに推定された波源位置との誤差が第４のしきい値以下であるとき、収束条件を満たすと判定し、推定処理を完了する。

（構成１０）
構成１から構成９のいずれかにおいて、波源位置推定装置は、通信推定手段を更に備える。通信推定手段は、１次利用者および２次利用者の通信方式または変調方式を推定する。推定手段は、通信推定手段による推定結果に基づいて、複数の実端末装置の複数の受信電力のうち誤差を抑制した複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

(構成１１)
構成１０において、推定手段は、通信推定手段によって推定された通信方式が時分割複信であるとき、複数の実端末装置の複数の受信電力のうち、他の通信装置を制御および／または管理する通信装置から送信された電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

（構成１２）
構成１０において、推定手段は、通信推定手段によって推定された通信方式が周波数分割複信であるとき、複数の実端末装置の複数の受信電力のうち、上りリンクにおける電波の複数の受信電力を用いず、下りリンクにおける電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

（構成１３）
構成１０において、推定手段は、通信推定手段によって推定された変調方式がＦＭ変調であるとき、複数の実端末装置の複数の受信電力のうち、変調信号による周波数偏移を含む帯域幅を有するフィルタを用いて検出された複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

(構成１４)
また、この発明の実施の形態によれば、プログラムは、受信手段が、移動端末または静止端末からなる実端末装置の位置情報と実端末装置における受信電力とを含む実端末情報を複数の実端末装置から受信する第１のステップと、推定手段が、複数の実端末装置の複数の実端末情報に含まれる複数の位置情報と複数の受信電力とに基づいて受信電力によって重み付けされた複数の実端末装置の重心を波源位置として推定し、その推定した波源位置に誤差があると判定したとき、少なくとも１つの仮想端末装置の位置を示す位置情報と受信電力とを追加して複数の実端末装置および少なくとも１つの仮想端末装置の位置情報および受信電力に基づいて受信電力によって重み付けされた複数の実端末装置および少なくとも１つの仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する第２のステップとをコンピュータに実行させるためのプログラムである。

（構成１５）
構成１４において、推定手段は、第２のステップにおいて、電力分布の推定結果を示す電力分布推定結果に基づいて少なくとも１つの仮想端末装置の位置情報と受信電力とを追加し、その追加した少なくとも１つの仮想端末装置の位置情報と受信電力と、複数の実端末装置の複数の位置情報と複数の受信電力とに基づいて推定処理を収束条件が満たされるまで繰り返し実行する。

(構成１６)
構成１５において、推定手段は、第２のステップにおいて、少なくとも１つの仮想端末装置をランダムに配置し、その配置した位置を示す位置情報を前記少なくとも１つの仮想端末装置の位置情報とし、電力分布推定結果に基づいて、少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を前記少なくとも１つの仮想端末装置の受信電力とする。

（構成１７）
構成１５において、推定手段は、第２のステップにおいて、電力分布推定結果によって示される等電力線上に少なくとも１つの仮想端末装置を配置し、等電力線上の電力値を少なくとも１つの仮想端末装置の受信電力とし、少なくとも１つの仮想端末装置を配置した位置を示す位置情報を少なくとも１つの仮想端末装置の位置情報とする。

（構成１８）
構成１５において、推定手段は、第２のステップにおいて、複数の実端末装置のうち、隣接する実端末装置間の距離を求め、その求めた距離が第１のしきい値以下である実端末装置が配置されている領域以外の領域に少なくとも１つの仮想端末装置を配置し、その配置した位置を示す位置情報を少なくとも１つの仮想端末装置の位置情報とし、電力分布推定結果に基づいて、少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を少なくとも１つの仮想端末装置の受信電力とする。

（構成１９）
構成１５から構成１８のいずれかにおいて、推定手段は、第２のステップにおいて、電力分布推定結果によって示される等電力線上の電力値と実端末装置における受信電力との電力差が第２のしきい値よりも大きいか否かを複数の実端末装置の全てについて判定し、電力差が第２のしきい値よりも大きい実端末装置を除外して推定処理を収束条件が満たされるまで繰り返し実行する。

（構成２０）
構成１４から構成１９のいずれかにおいて、推定手段は、第２のステップにおいて、複数の実端末装置および少なくとも１つの仮想端末装置の位置の平均が第３のしきい値以下であるとき、収束条件を満たすと判定し、推定処理を完了する。

（構成２１）
構成１４から構成１９のいずれかにおいて、推定手段は、第２のステップにおいて、推定処理の実行回数が所定回数に達したとき、収束条件を満たすと判定し、推定処理を完了する。

（構成２２）
構成１４から構成１９のいずれかにおいて、推定手段は、第２のステップにおいて、推定処理を２回以上実行したとき、推定処理をｋ−１（ｋは２以上の整数）回実行したときに推定された波源位置と、推定処理をｋ回実行したときに推定された波源位置との誤差が第４のしきい値以下であるとき、収束条件を満たすと判定し、推定処理を完了する。

（構成２３）
構成１４から構成２２のいずれかにおいて、通信推定手段が、１次利用者および２次利用者の通信方式または変調方式を推定する第３のステップを更にコンピュータに実行させ、推定手段は、第２のステップにおいて、通信推定手段による推定結果に基づいて、複数の実端末装置の複数の受信電力のうち誤差を抑制した複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

（構成２４）
構成２３において、推定手段は、第２のステップにおいて、通信推定手段によって推定された通信方式が時分割複信であるとき、複数の実端末装置の複数の受信電力のうち、他の通信装置を制御および／または管理する通信装置から送信された電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

（構成２５）
構成２３において、推定手段は、第２のステップにおいて、通信推定手段によって推定された通信方式が周波数分割複信であるとき、複数の実端末装置の複数の受信電力のうち、上りリンクにおける電波の複数の受信電力を用いず、下りリンクにおける電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

(構成２６)
構成２３において、推定手段は、第２のステップにおいて、通信推定手段によって推定された変調方式がＦＭ変調であるとき、複数の実端末装置の複数の受信電力のうち、変調信号による周波数偏移を含む帯域幅を有するフィルタを用いて検出された複数の受信電力を選択し、その選択した複数の受信電力を用いて推定処理を収束条件が満たされるまで繰り返し実行する。

（構成２７）
更に、この発明の実施の形態によれば、記録媒体は、構成１４から構成２６のいずれかに記載のプログラムを記録したコンピュータ読み取り可能な記録媒体である。

計算量を抑制して正確に波源位置を推定できる。

この発明の実施の形態１における無線通信システムを示す図である。 図１に示す波源位置推定装置の概略図である。 受信電力と波源からの距離との関係を示す図である。 重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置とする推定法における計算結果を示す図である。 重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置とする推定法の問題点を説明するための図である。 電力分布推定結果を示す図である。 仮想端末装置をランダムに配置する場合の位置情報および受信電力の決定方法を説明するための図である。 等電力線上に仮想端末装置を配置する場合の位置情報および受信電力の決定方法を説明するための図である。 端末装置２間の距離がしきい値Ｌｔｈ２以下である場合の位置情報および受信電力の決定方法を説明するための図である。 波源位置の推定誤差と判定回数との関係を示す図である。 図２に示す波源位置推定装置の動作を説明するためのフローチャートである。 図１１に示すステップＳ５の詳細な動作を説明するためのフローチャートである。 図１１に示すステップＳ５の別の詳細な動作を説明するためのフローチャートである。 図１１に示すステップＳ５の更に別の詳細な動作を説明するためのフローチャートである。 図１１のステップＳ６の詳細な動作を説明するためのフローチャートである。 図２に示す波源位置推定装置の動作を説明するための別のフローチャートである。 図２に示す波源位置推定装置の動作を説明するための更に別のフローチャートである。 実施の形態１による別の波源位置推定装置の概略図である。 受信電力の等高線を作成する方法を説明するための図である。 ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄの抽出方法を説明するための図である。 受信電力の等高線ＣＴＲの作成方法２を説明するための図である。 ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄの別の抽出方法を説明するための図である。 ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄの抽出方式と受信電力の等高線ＣＴＲの作成処理との関係を示す図である。 受信電力の等高線ＣＴＲを作成する別の方法を説明するための図である。 端末装置２における電波の受信感度以下の領域において受信電力を推定する方法を説明するための図である。 電力分布推定結果を作成する方法を説明するためのフローチャートである。 図２６のステップＳ１０１の詳細な動作を説明するためのフローチャートである。 図２６のステップＳ１０４の詳細な動作を説明するためのフローチャートである。 図１８に示す波源位置推定装置の動作を説明するためのフローチャートである。 図１８に示す波源位置推定装置の動作を説明するための別のフローチャートである。 図１８に示す波源位置推定装置の動作を説明するための更に別のフローチャートである。 波源位置の推定誤差と端末装置の個数との関係を示す図である。 実施の形態２における無線通信システムを示す図である。 周波数と電力との関係を示す図である。 図３３に示す波源位置推定装置の概略図である。 通信方式および変調方式の推定方法を説明するための図である。 自己相関の概念図である。 自己相関の別の概念図である。 データベースの概略図である。 図３３に示す波源位置推定装置の動作を説明するためのフローチャートである。 図４０に示すステップＳ１２２の詳細な動作を説明するためのフローチャートである。 図４１のステップＳ１２２７の詳細な動作を説明するためのフローチャートである。 図４１のステップＳ１２２９の詳細な動作を説明するためのフローチャートである。 図４１のステップＳ１２３１の詳細な動作を説明するためのフローチャートである。 図４０のステップＳ１２２の別の詳細な動作を説明するためのフローチャートである。 図４０のステップＳ１２４の詳細な動作を説明するためのフローチャートである。 図３３に示す波源位置推定装置の別の動作を説明するためのフローチャートである。 図３３に示す波源位置推定装置の更に別の動作を説明するためのフローチャートである。 実施の形態２による別の波源位置推定装置の概略図である。 実施の形態２における受信電力の等高線を作成する方法を説明するための図である。 実施の形態２におけるｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄの抽出方法を説明するための図である。 受信電力の等高線ＣＴＲの作成方法３を説明するための図である。 実施の形態２におけるｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄの別の抽出方法を説明するための図である。 実施の形態２におけるｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄの抽出方式と受信電力の等高線ＣＴＲの作成処理との関係を示す図である。 実施の形態２における電力分布推定結果を作成する方法を説明するためのフローチャートである。 図５５に示すステップＳ１２６の詳細な動作を説明するためのフローチャートである。 図５５に示すステップＳ１２８の詳細な動作を説明するためのフローチャートである。 図５５に示すステップＳ１２７の詳細な動作を説明するためのフローチャートである。 図５５に示すステップＳ１２９の詳細な動作を説明するためのフローチャートである。 図４９に示す波源位置推定装置の動作を説明するためのフローチャートである。 図４９に示す波源位置推定装置の動作を説明するための別のフローチャートである。 図４９に示す波源位置推定装置の動作を説明するための更に別のフローチャートである。

本発明の実施の形態について図面を参照しながら詳細に説明する。なお、図中同一または相当部分には同一符号を付してその説明は繰返さない。

［実施の形態１］
図１は、この発明の実施の形態１における無線通信システムを示す図である。図１を参照して、無線通信システム１０は、波源位置推定装置１と、複数の端末装置２とを備える。

波源位置推定装置１および複数の端末装置２は、無線通信空間に配置される。複数の端末装置２および波源Ｓは、観測エリアＲＥＧ内に配置される。観測エリアＲＥＧは、波源Ｓを中心とした円形形状を有する。

複数の端末装置２の各々は、移動端末または静止端末からなり、実際に存在する実端末装置である。そして、複数の端末装置２の各々は、例えば、ＧＰＳ（Ｇｌｏｂａｌ Ｐｏｓｉｔｉｏｎｉｎｇ Ｓｙｓｔｅｍ）によって自己の位置を検出する。

また、複数の端末装置２の各々は、波源Ｓから電波を受信し、電波を受信したときの受信電力Ｐ_ｉ（ｉ＝１，２，３，・・・）を検出する。

更に、複数の端末装置２の各々は、受信電力Ｐ_ｉを検出したときの時刻を検出する。

そして、複数の端末装置２の各々は、自己の識別情報ＩＤ_ｉと、自己の位置を示す位置情報（ｘ_ｉ，ｙ_ｉ）と、受信電力Ｐ_ｉと、受信電力Ｐ_ｉを検出したときの時刻を示す時間情報ｔ_ｉとを含む端末情報ＴＩＦ＝［ＩＤ_ｉ／（ｘ_ｉ，ｙ_ｉ）／Ｐ_ｉ／ｔ_ｉ］を生成し、その生成した端末情報ＴＩＦ＝［ＩＤ_ｉ／（ｘ_ｉ，ｙ_ｉ）／Ｐ_ｉ／ｔ_ｉ］を無線通信によって波源位置推定装置１へ送信する。

波源位置推定装置１は、複数の端末装置２の各々から端末情報ＴＩＦ＝［ＩＤ_ｉ／（ｘ_ｉ，ｙ_ｉ）／Ｐ_ｉ／ｔ_ｉ］を受信し、その受信した端末情報ＴＩＦ＝［ＩＤ_ｉ／（ｘ_ｉ，ｙ_ｉ）／Ｐ_ｉ／ｔ_ｉ］に基づいて、後述する方法によって、波源Ｓの位置である波源位置を推定する。

図２は、図１に示す波源位置推定装置１の概略図である。図２を参照して、波源位置推定装置１は、アンテナ１１と、受信手段１２と、記憶手段１３と、推定手段１４とを備える。

受信手段１２は、アンテナ１１を介して端末装置２の各々から端末情報ＴＩＦ＝［ＩＤ_ｉ／（ｘ_ｉ，ｙ_ｉ）／Ｐ_ｉ／ｔ_ｉ］を受信し、その受信した端末情報ＴＩＦ＝［ＩＤ_ｉ／（ｘ_ｉ，ｙ_ｉ）／Ｐ_ｉ／ｔ_ｉ］を記憶手段１３に記憶する。

記憶手段１３は、受信手段１２から受けた端末情報ＴＩＦ＝［ＩＤ_ｉ／（ｘ_ｉ，ｙ_ｉ）／Ｐ_ｉ／ｔ_ｉ］を記憶する。

推定手段１４は、複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］（Ｎは、端末装置２の総数を表す２以上の整数であり、ｉ＝１〜Ｎ）を記憶手段１３から読み出す。そして、推定手段１４は、複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］に含まれる複数の位置情報（ｘ_１，ｙ_１）〜（ｘ_Ｎ，ｙ_Ｎ）および複数の受信電力Ｐ_１〜Ｐ_Ｎに基づいて波源位置を推定する。

波源位置の推定方法について説明する。

推定手段１４は、次式によって、複数の端末装置２の重心（ｘ_Ｇ，ｙ_Ｇ）を算出する。

式（１）において、ｗ_ｉは、ｉ番目の端末装置２のウェイトを表し、ｉ番目の端末装置２における受信電力Ｐ_ｉからなる。

そして、推定手段１４は、式（１）によって算出した重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置として推定する。即ち、推定手段１４は、受信電力Ｐ_ｉによって重み付けされた複数の端末装置２の重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置として推定する。

図３は、受信電力と波源からの距離との関係を示す図である。図３において、縦軸は、受信電力を表し、横軸は、波源Ｓからの距離を表す。

図３を参照して、受信電力Ｐ_ｉは、波源Ｓからの距離に反比例して減衰する。その結果、受信電力Ｐ_ｉは、波源Ｓに近づくほど大きくなるので、波源Ｓの近傍の方が受信電力Ｐ_ｉによる重み付けの効果が大きい。従って、重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置とする推定法によって推定された推定値は、波源近傍に近づく。

図４は、重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置とする推定法における計算結果を示す図である。

波源Ｓは、ｘ−ｙ座標において（３００，３００）の位置に配置された。端末装置２の台数は、５０台である。図４の白丸は、端末装置２を表す。観測エリアの半径は、１０００ｍである。

式（１）を用いて計算した重心（ｘ_Ｇ，ｙ_Ｇ）の位置は、（２９５，３０６）であった。

従って、式（１）を用いて計算した重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置とすることにより、推定値は、波源Ｓの近傍に位置することが分かった。

図５は、重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置とする推定法の問題点を説明するための図である。

図５において、丸は、重心位置を表し、四角は、波源位置を表す。

図５の（ａ）を参照して、端末装置２の分布が一様である場合、重心位置は、波源位置の近傍に精度良く位置する。

一方、端末装置２の分布に偏りがある場合、重心位置は、波源位置と大きく異なる。式（１）を用いて算出された重心（ｘ_Ｇ，ｙ_Ｇ）は、端末装置２が分布している領域のほど中央部に位置するので、複数の端末装置２の分布領域が波源位置に対して一方側に偏っている場合、重心位置は、波源位置の近傍に位置しない（図５の（ｂ）参照）。

そこで、推定手段１４は、推定した波源位置と原点との距離Ｌ_Ｓ−Ｏを演算し、その演算した距離Ｌ_Ｓ−Ｏがしきい値Ｌｔｈ１よりも大きいとき、推定した波源位置に誤差があると判定し、後述する方法によって、仮想端末装置を追加し、その追加した仮想端末装置の位置を示す位置情報および受信電力と、複数の端末装置２の位置情報および受信電力とに基づいて、式（１）を用いた受信電力によって重み付けされた複数の端末装置２および仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する。

距離Ｌ_Ｓ−Ｏがしきい値Ｌｔｈ１よりも大きいか否かを判定するとき、距離Ｌ_Ｓ−Ｏをｘ軸成分Ｌ_{Ｓ−Ｏ＿ｘ}とｙ軸成分Ｌ_{Ｓ−Ｏ＿ｙ}とに分け、しきい値Ｌｔｈ１をｘ軸成分Ｌｔｈ１＿ｘとｙ軸成分Ｌｔｈ１＿ｙとに分ける。そして、Ｌ_{Ｓ−Ｏ＿ｘ}がＬｔｈ１＿ｘよりも大きいか否かを判定するとともに、Ｌ_{Ｓ−Ｏ＿ｙ}がＬｔｈ１＿ｙよりも大きいか否かを判定する。そして、Ｌ_{Ｓ−Ｏ＿ｘ}がＬｔｈ１＿ｘよりも大きく、かつ、Ｌ_{Ｓ−Ｏ＿ｙ}がＬｔｈ１＿ｙよりも大きいとき、推定した波源位置に誤差があると判定する。一方、Ｌ_{Ｓ−Ｏ＿ｘ}＞Ｌｔｈ１＿ｘおよびＬ_{Ｓ−Ｏ＿ｙ}＞Ｌｔｈ１＿ｙの少なくとも一方が満たされないとき、推定した波源位置に誤差がないと判定する。この場合、Ｌｔｈ１＿ｘおよびＬｔｈ１＿ｙの各々は、端末装置２が自己の位置をＧＰＳで検出しているとき、例えば、１０ｍに設定される。

なお、距離Ｌ_Ｓ−Ｏを演算し、その演算した距離Ｌ_Ｓ−Ｏがしきい値Ｌｔｈ１よりも大きいとき、推定した波源位置に誤差があると判定するのは、次の理由による。波源位置を推定するとき、複数の端末装置２を配置したｘ−ｙ平面の原点を波源の位置とするので、推定した波源位置と原点との距離Ｌ_Ｓ−Ｏがしきい値Ｌｔｈ１よりも大きいとき、推定した波源位置に誤差があると判定できるからである。

仮想端末装置を選択する方法について説明する。図６は、電力分布推定結果を示す図である。図６を参照して、白四角は、波源位置を示す。そして、波源位置を中心として、複数の等電力線ＣＴＲ１，ＣＴＲ２が同心円状に配置されている。そして、等電力線ＣＴＲ１は、電力値ＰＷ１を有し、等電力線ＣＴＲ２は、電力値ＰＷ２（＜ＰＷ１）を有する。従って、電力分布推定結果は、図６に示すように、複数の等電力線ＣＴＲ１，ＣＴＲ２を有し、電力分布の推定結果を示すものである。

なお、図６においては、２つの等電力線ＣＴＲ１，ＣＴＲ２が示されているが、電力分布推定結果は、３個以上の等電力線を含んでいてもよい。

推定手段１４は、次のいずれかの方法によって仮想端末装置を選択する。
（１）仮想端末装置をランダムに配置する。
（２）等電力線上に仮想端末装置を配置する。
（３）複数の端末装置２のうち、隣接する端末装置２間の距離を求め、その求めた距離がしきい値Ｌｔｈ２以下である複数の端末装置２は、偏っていると判定できるので、その偏っていると判定された複数の端末装置２が配置された領域以外の領域に仮想端末装置を配置する。なお、しきい値Ｌｔｈ２は、例えば、１００ｍである。

図７は、仮想端末装置をランダムに配置する場合の位置情報および受信電力の決定方法を説明するための図である。

図７を参照して、推定手段１４は、Ｊ（Ｊは、正の整数）個の仮想端末装置２ｖｊ（ｊ＝１〜Ｊ）をランダムに配置する。そして、推定手段１４は、各仮想端末装置２ｖｊを配置した位置を示す位置情報（ｘ_ｊｖ，ｙ_ｊｖ）を決定する。また、推定手段１４は、仮想端末装置２ｖｊから最も近い等高線を選択し、仮想端末装置２ｖｊと等高線との距離Ｌ_ｃ−ｖを演算し、その演算した距離Ｌ_ｃ−ｖに基づいて、電波の自由空間伝搬モデルを用いて仮想端末装置２ｖにおける受信電力Ｐ_ｊｖを演算する。

電波の自由空間伝搬モデルでは、受信電力Ｐ_ｉは、図３に示すように波源Ｓからの距離に反比例して減少する。

電波の自由空間伝搬モデルにおいて、波源Ｓからの距離が１０００ｍ以上である場合、受信電力Ｐ_ｉは、距離の増加に対して直線状に低下する。そして、受信電力Ｐ_ｉが１０ｄＢだけ低下するときの距離は、３．１６ｋｍである。つまり、電波の自由空間伝搬モデルにおいては、受信電力Ｐ_ｉが１０ｄＢだけ低下するときの距離は、３．１６ｋｍ／１０ｄＢｍとなる。

従って、推定手段１４Ａは、既に作成した受信電力の等高線ＣＴＲ＿ＳＹＮ（電力値ＰＷの等高線）から距離Ｌ_ｃ−ｖだけ波源から遠くなる方向に離れた位置に存在する仮想端末装置２ｖｊにおける受信電力Ｐ_ｊｖをＰＷ−１０Ｌ_ｃ−ｖ／３．１６［ｄＢｍ］によって演算する。また、推定手段１４Ａは、等高線ＣＴＲ＿ＳＹＮ（電力値がＰＷ１またはＰＷ２の等高線）から距離Ｌ_ｃ−ｖだけ波源に近づいた位置に存在する仮想端末装置２ｖｊにおける受信電力Ｐ_ｊｖをＰＷ＋１０Ｌ_ｃ−ｖ／３．１６［ｄＢｍ］によって演算する。

そして、推定手段１４は、位置情報（ｘ_ｊｖ，ｙ_ｊｖ）の決定と、受信電力Ｐ_ｊｖの算出とを全ての仮想端末装置２ｖ１〜２ｖＪについて実行する。

なお、推定手段１４は、距離Ｌ_ｃ−ｖを演算する場合、等高線に垂直になるように仮想端末装置２ｖから等高線に線を引き、その引いた線と等高線との交点の位置と位置（ｘ_ｊｖ，ｙ_ｊｖ）との距離を演算することによって距離Ｌ_ｃ−ｖを演算する。

図８は、等電力線上に仮想端末装置を配置する場合の位置情報および受信電力の決定方法を説明するための図である。

図８を参照して、推定手段１４は、Ｊ個の仮想端末装置２ｖｊを等高線ＣＴＲ１，ＣＴＲ２上に配置する。そして、推定手段１４は、各仮想端末装置２ｖｊを配置した位置を示す位置情報（ｘ_ｊｖ，ｙ_ｊｖ）を決定する。また、推定手段１４は、各仮想端末装置２ｖｊを配置した等高線（等高線ＣＴＲ１，ＣＴＲ２のいずれか）上の電力値（ＰＷ１，ＰＷ２のいずれか）を仮想端末装置２ｖｊにおける受信電力Ｐ_ｊｖとして決定する。そして、推定手段１４は、位置情報（ｘ_ｊｖ，ｙ_ｊｖ）の決定と、受信電力Ｐ_ｊｖの決定とを全ての仮想端末装置２ｖ１〜２ｖＪについて実行する。

図９は、端末装置２間の距離がしきい値Ｌｔｈ２以下である場合の位置情報および受信電力の決定方法を説明するための図である。

図９を参照して、推定手段１４は、複数の端末装置２において、隣接する端末装置２間の距離Ｌ_２−２を演算する。推定手段１４は、距離Ｌ_２−２の演算を全ての隣接する端末装置２間について実行する。そして、推定手段１４は、その演算した複数の距離Ｌ_２−２の少なくとも１つがしきい値Ｌｔｈ２以下であるとき、複数の端末装置２が偏っていると判定する。

そうすると、推定手段１４は、複数の端末装置２が配置されている領域ＲＥＧ１以外の領域にＪ個の仮想端末装置２ｖｊを配置する。この場合、推定手段１４は、仮想端末装置２ｖｊを等高線ＣＴＲ１，ＣＴＲ２上に配置してもよく、等高線ＣＴＲ１，ＣＴＲ２以外の位置に配置してもよい。

そして、推定手段１４は、各仮想端末装置２ｖｊを配置した位置を示す位置情報（ｘ_ｊｖ，ｙ_ｊｖ）を決定する。

また、推定手段１４は、各仮想端末装置２ｖｊを等高線ＣＴＲ１，ＣＴＲ２以外の位置に配置した場合、上述した距離Ｌ_ｃ−ｖを演算し、その演算した距離Ｌ_ｃ−ｖに基づいて、上述した方法によって仮想端末装置２ｖにおける受信電力Ｐ_ｊｖを演算する。一方、推定手段１４は、各仮想端末装置２ｖｊを等高線ＣＴＲ１，ＣＴＲ２上に配置した場合、各仮想端末装置２ｖｊを配置した等高線（等高線ＣＴＲ１，ＣＴＲ２のいずれか）上の電力値（ＰＷ１，ＰＷ２のいずれか）を仮想端末装置２ｖｊにおける受信電力Ｐ_ｊｖとして決定する。推定手段１４は、位置情報（ｘ_ｊｖ，ｙ_ｊｖ）の決定と、受信電力Ｐ_ｊｖの決定とを全ての仮想端末装置２ｖ１〜２ｖＪについて実行する。

推定手段１４は、仮想端末装置２ｖｊの位置情報（ｘ_ｊｖ，ｙ_ｊｖ）および受信電力Ｐ_ｊｖを決定すると、端末装置２および仮想端末装置２ｖｊの位置情報（ｘ_ｉ，ｙ_ｉ），（ｘ_ｊｖ，ｙ_ｊｖ）および受信電力Ｐ_ｉ，Ｐ_ｊｖを用いて式（１）によって重心（ｘ_Ｇ，ｙ_Ｇ）を算出する。

推定手段１４は、重心（ｘ_Ｇ，ｙ_Ｇ）を算出すると、収束条件を満たすか否かを判定する。収束条件は、次のとおりである。

（Ｃｄ１）端末装置２および仮想端末装置２ｖｊを含む全ての端末装置の位置の平均（ｘ_Ａ，ｙ_Ａ）がしきい値Δｘ_Ａ，Δｙ_Ａ以下であるとき、収束したと判定する。この場合、ｘ_Ａがしきい値Δｘ_Ａ以下であり、かつ、ｙ_Ａがしきい値Δｙ_Ａ以下である場合、収束したと判定する。しきい値Δｘ_Ａは、例えば、１０ｍであり、しきい値Δｙ_Ａは、例えば、１０ｍである。

（Ｃｄ２）判定回数が５回〜１０回のように一定回数に達したときに収束したと判定する。

（Ｃｄ３）２回目以降の判定において、ｋ−１（ｋは２以上の整数）回目の判定時において推定された波源位置と、ｋ回目の判定時において推定された波源位置との誤差ΔＧがしきい値Ｇｔｈ以下であるときに収束したと判定する。

収束条件（Ｃｄ１）が用いられる場合、位置の平均（ｘ_Ａ，ｙ_Ａ）は、次式によって算出される。

式（２）において、ＳＮ_ｒは、端末装置２を表し、ＳＮ_ｄは、仮想端末装置２ｖｊを表す。

図１０は、波源位置の推定誤差と判定回数との関係を示す図である。図１０において、縦軸は、波源位置の推定誤差を表し、横軸は、判定回数を表す。曲線ｋ１は、波源位置の推定誤差と判定回数との関係を示す。

図１０を参照して、波源位置の推定誤差ΔＧは、判定回数が増加するに従って指数関数的に減少する。そして、上記収束条件（Ｃｄ３）においては、波源位置の推定誤差ΔＧがしきい値Ｇｔｈ以下になると、収束したと判定する。なお、しきい値Ｇｔｈは、例えば、１０ｍである。

フェージングの影響を受けている端末装置２は、波源位置の推定精度を低下させる要因になる。そこで、この発明の実施の形態においては、受信電力の信頼度ＲＥＬを算出し、その算出した信頼度ＲＥＬがしきい値Ｒｔｈ以下である端末装置２を波源位置の推定に用いる端末装置２から削除する。

推定手段１４は、端末装置２が等電力線上に配置されている場合、等電力線上の電力値ＰＷ（電力値ＰＷ１，ＰＷ２のいずれか）と、各端末装置２における受信電力Ｐ_ｉとの差ΔＰ_ｉを演算する。また、推定手段１４は、端末装置２が等電力線上に配置されていない場合、等電力線上の電力値ＰＷ（電力値ＰＷ１，ＰＷ２のいずれか）に基づいて電波の自由空間伝搬モデルを用いて決定された電力値（ＰＷ−１０Ｌ_ｃ−ｖ／３．１６またはＰＷ＋１０Ｌ_ｃ−ｖ／３．１６）と、各端末装置２における受信電力Ｐ_ｉとの差ΔＰ_ｉを演算する。そして、推定手段１４は、差ΔＰ_ｉがしきい値Ｐｔｈよりも大きいか否かを判定することによって信頼度ＲＥＬ_ｉがしきい値Ｒｔｈ以下であるか否かを判定する。より具体的には、推定手段１４は、差ΔＰ_ｉがしきい値ΔＰｔｈよりも大きいとき、信頼度ＲＥＬ_ｉがしきい値Ｒｔｈ以下であると判定し、差ΔＰ_ｉがしきい値ΔＰｔｈ以下であるとき、信頼度ＲＥＬ_ｉがしきい値Ｒｔｈよりも大きいと判定する。なお、しきい値Ｐｔｈは、例えば、１０ｄＢである。

このように、この発明の実施の形態においては、推定手段１４は、複数の端末装置２の位置情報（ｘ_ｉ，ｙ_ｉ）および受信電力Ｐ_ｉに基づいて推定した波源位置に誤差がある場合、仮想端末装置２ｖｊの追加と、受信電力の信頼度ＲＥＬがしきい値Ｒｔｈ以下である端末装置２の削除とを行い、端末装置２および仮想端末装置２ｖｊの位置情報（ｘ_ｉ，ｙ_ｉ），（ｘ_ｊｖ，ｙ_ｊｖ）および受信電力Ｐ_ｉ，Ｐ_ｊｖに基づいて式（１）によって演算した重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置と推定する推定処理を収束条件が満たされるまで繰り返し行う。

図１１は、図２に示す波源位置推定装置１の動作を説明するためのフローチャートである。図１１に示すフローチャートは、端末装置２および仮想端末装置２ｖｊを含む全ての端末装置の位置の平均（ｘ_Ａ，ｙ_Ａ）がしきい値Δｘ_Ａ，Δｙ_Ａ以下であるとき、収束したと判定するフローチャートである。

図１１を参照して、波源位置推定装置１の動作が開始されると、波源位置推定装置１の受信手段１２は、アンテナ１１を介して複数の端末装置２から複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］を受信することにより、端末装置２の位置情報（ｘ_ｉ，ｙ_ｉ）と受信電力Ｐ_ｉを取得する（ステップＳ１）。そして、受信手段１２は、その受信した複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］を記憶手段１３に記憶する。

推定手段１４は、複数の位置情報（ｘ_１，ｙ_１）〜（ｘ_Ｎ，ｙ_Ｎ）および複数の受信電力Ｐ_１〜Ｐ_Ｎを記憶手段１３から読み出し、その読み出した複数の位置情報（ｘ_１，ｙ_１）〜（ｘ_Ｎ，ｙ_Ｎ）および複数の受信電力Ｐ_１〜Ｐ_Ｎを式（１）に代入して重心（ｘ_Ｇ，ｙ_Ｇ）を演算し、その演算した重心（ｘ_Ｇ，ｙ_Ｇ）を波源位置として推定する（ステップＳ２）。

そして、推定手段１４は、複数の位置情報（ｘ_１，ｙ_１）〜（ｘ_Ｎ，ｙ_Ｎ）および端末装置２の総数Ｎを式（２）に代入して位置の平均（ｘ_Ａ，ｙ_Ａ）を演算する（ステップＳ３）。

そうすると、推定手段１４は、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａであるか否かを判定する（ステップＳ４）。

ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａでないと判定されたとき、推定手段１４は、予め保持した電力分布推定結果に基づいてＪ個の仮想端末装置２ｖ１〜２ｖＪを追加する（ステップＳ５）。

引き続いて、推定手段１４は、受信電力の信頼度ＲＥＬが低い端末装置２を波源位置の推定に用いる端末装置から削除する（ステップＳ６）。

その後、一連の動作は、ステップＳ２へ移行し、ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａであると判定されるまで、ステップＳ２〜ステップＳ６が繰り返し実行される。

そして、ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａであると判定されると、波源位置を推定する動作が終了する。

なお、ステップＳ６からステップＳ２へ移行した場合、仮想端末装置２ｖｊの位置情報（ｘ_ｊｖ，ｙ_ｊｖ）および受信電力Ｐ_ｊｖと、端末装置２の位置情報（ｘ_ｉ，ｙ_ｉ）および受信電力Ｐ_ｉとを式（１）に代入して重心（ｘ_Ｇ，ｙ_Ｇ）を演算することによって波源位置が推定される。

図１２は、図１１に示すステップＳ５の詳細な動作を説明するためのフローチャートである。なお、図１２に示すフローチャートは、仮想端末装置２ｖｊをランダムに配置することによって仮想端末装置２ｖｊを追加する動作を説明するためのフローチャートである。

図１２を参照して、図１１のステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａでないと判定されたとき、推定手段１４は、Ｊ個の仮想端末装置２ｖ１〜２ｖＪをランダムに配置する（ステップＳ５１）。

そして、推定手段１４は、ｊ＝１を設定し（ステップＳ５２）、仮想端末装置２ｖｊの配置位置に基づいて仮想端末装置２ｖｊの位置情報（ｘ_ｊｖ，ｙ_ｊｖ）を決定する（ステップＳ５３）。この場合、Ｊ個の仮想端末装置２ｖ１〜２ｖＪは、ｘ−ｙ平面上に配置されているので、推定手段１４は、仮想端末装置２ｖｊが配置されたｘ座標およびｙ座標を検出することによって位置情報（ｘ_ｊｖ，ｙ_ｊｖ）を決定する。

ステップＳ５３の後、推定手段１４は、仮想端末装置２ｖｊが等電力線上に配置されているか否かを判定する（ステップＳ５４）。

ステップＳ５４において、仮想端末装置２ｖｊが等電力線上に配置されていないと判定されたとき、推定手段１４は、仮想端末装置２ｖｊと等電力線との距離Ｌ_ｃ−ｖに基づいて電波の自由空間伝搬モデルを用いて上述した方法によって仮想端末装置２ｖｊの位置における電力値を算出し、その算出した電力値を仮想端末装置２ｖｊの受信電力Ｐ_ｊｖとして決定する（ステップＳ５５）。

一方、ステップＳ５４において、仮想端末装置２ｖｊが等電力線上に配置されていると判定されたとき、推定手段１４は、等電力線上の電力値を仮想端末装置２ｖｊの受信電力Ｐ_ｊｖとして決定する（ステップＳ５６）。

そして、ステップＳ５５またはステップＳ５６の後、推定手段１４は、ｊ＝Ｊであるか否かを判定する（ステップＳ５７）。

ステップＳ５７において、ｊ＝Ｊでないと判定されたとき、推定手段１４は、ｊ＝ｊ＋１を設定する（ステップＳ５８）。その後、一連の動作は、ステップＳ５３へ移行し、ステップＳ５７において、ｊ＝Ｊであると判定されるまで、ステップＳ５３〜ステップＳ５８が繰り返し実行される。そして、ステップＳ５７において、ｊ＝Ｊであると判定されると、一連の動作は、図１１のステップＳ６へ移行する。

図１３は、図１１に示すステップＳ５の別の詳細な動作を説明するためのフローチャートである。なお、図１３に示すフローチャートは、仮想端末装置２ｖｊを等電力線上に配置することによって仮想端末装置２ｖｊを追加する動作を説明するためのフローチャートである。

図１３に示すフローチャートは、図１２に示すフローチャートのステップＳ５１をステップＳ５１Ａに変え、ステップＳ５４，Ｓ５５を削除したものであり、その他は、図１２に示すフローチャートと同じである。従って、図１３に示すフローチャートにおいて、ステップＳ５３の後に、ステップＳ５６が実行され、ステップＳ５６の後に、ステップＳ５７，Ｓ５８が順次実行される。

図１３を参照して、図１１のステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａでないと判定されたとき、推定手段１４は、Ｊ個の仮想端末装置２ｖ１〜２ｖＪを等電力線上に配置する（ステップＳ５１Ａ）。

その後、上述したステップＳ５２，Ｓ５３，Ｓ５６，Ｓ５７，Ｓ５８が順次実行される。そして、ステップＳ５７において、ｊ＝Ｊであると判定されると、一連の動作は、図１１のステップＳ６へ移行する。

図１４は、図１１に示すステップＳ５の更に別の詳細な動作を説明するためのフローチャートである。なお、図１４に示すフローチャートは、複数の端末装置２の配置位置に偏りがある場合、複数の端末装置２の配置領域以外の領域に仮想端末装置２ｖｊを配置することによって仮想端末装置２ｖｊを追加する動作を説明するためのフローチャートである。

図１４に示すフローチャートは、図１２に示すフローチャートのステップＳ５１をステップＳ６１〜ステップＳ６４に変えたものであり、その他は、図１２に示すフローチャートと同じである。

図１４を参照して、図１１のステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａでないと判定されたとき、推定手段１４は、複数の端末装置２のうち、隣接する端末装置間の距離Ｌ_２−２を全ての隣接する端末装置間について演算する（ステップＳ６１）。

そして、推定手段１４は、少なくとも１つの距離Ｌ_２−２がしきい値Ｌｔｈ２以下であるか否かを判定する（ステップＳ６２）。なお、ステップＳ６２において、少なくとも１つの距離Ｌ_２−２がしきい値Ｌｔｈ２以下であるか否かを判定するのは、端末装置２の総数Ｎは、２以上であるので、隣接する端末装置間の距離Ｌ_２−２は、１個以上であるからである。

ステップＳ６２において、少なくとも１つの距離Ｌ_２−２がしきい値Ｌｔｈ２以下であると判定されたとき、複数の端末装置２の配置位置に偏りがあるので、推定手段１４は、複数の端末装置２の配置領域以外の領域にＪ個の仮想端末装置２ｖ１〜２ｖＪを配置する（ステップＳ６３）。

一方、ステップＳ６２において、少なくとも１つの距離Ｌ_２−２がしきい値Ｌｔｈ２以下でないと判定されたとき、複数の端末装置２の配置位置に偏りがないので、推定手段１４は、複数の端末装置２の配置領域も含めてＪ個の仮想端末装置２ｖ１〜２ｖＪを配置する（ステップＳ６４）。

そして、ステップＳ６３またはステップＳ６４の後、ステップＳ５２へ移行し、上述したステップＳ５２〜ステップＳ５８が順次実行される。そして、ステップＳ５７において、ｊ＝Ｊであると判定されると、一連の動作は、図１１のステップＳ６へ移行する。

なお、図１４のステップＳ６３またはステップＳ６４においては、Ｊ個の仮想端末装置２ｖ１〜２ｖＪは、ランダムまたは等高線上に配置される。そして、Ｊ個の仮想端末装置２ｖ１〜２ｖＪがランダムまたは等高線上に配置された場合でも、ステップＳ５３を実行することによって仮想端末装置２ｖｊの位置情報（ｘ_ｊｖ，ｙ_ｊｖ）を決定でき、ステップＳ５４〜ステップＳ５６を実行することによって仮想端末装置２ｖｊの受信電力Ｐ_ｊｖを決定できる。

図１５は、図１１のステップＳ６の詳細な動作を説明するためのフローチャートである。

図１５を参照して、図１１のステップＳ５の後、推定手段１４は、ｉ＝１を設定する（ステップＳ７１）。

そして、推定手段１４は、端末装置２が等高線上に配置されているか否かを判定する（ステップＳ７２）。

ステップＳ７２において、端末装置２が等高線上に配置されていないと判定されたとき、推定手段１４は、端末装置２＿ｉに最も近い等高線ＣＴＲ＿ｎｅａｒｅｓｔ上の電力値ＰＷ１を電力分布推定結果に基づいて検出する（ステップＳ７３）。

その後、推定手段１４は、電力値ＰＷ１に基づいて電波の自由空間伝搬モデルを用いて算出した電力値と端末装置２＿ｉの受信電力Ｐ_ｉとの差ΔＰ_ｉを演算する（ステップＳ７４）。

一方、ステップＳ７２において、端末装置２が等高線上に配置されていると判定されたとき、推定手段１４は、端末装置２＿ｉが配置された等高線ＣＴＲ上の電力値ＰＷ２を電力分布推定結果に基づいて検出する（ステップＳ７５）。

その後、推定手段１４は、電力値ＰＷ２と端末装置２＿ｉの受信電力Ｐ_ｉとの差ΔＰ_ｉを演算する（ステップＳ７６）。

そして、ステップＳ７４またはステップＳ７６の後、推定手段１４は、差ΔＰ_ｉがしきい値ΔＰｔｈよりも大きいか否かを判定することによって信頼度ＲＥＬ_ｉがしきい値Ｒｔｈ以下であるか否かを判定する（ステップＳ７７）。この場合、推定手段１４は、差ΔＰ_ｉがしきい値ΔＰｔｈよりも大きいとき、信頼度ＲＥＬ_ｉがしきい値Ｒｔｈ以下であると判定し、差ΔＰ_ｉがしきい値ΔＰｔｈ以下であるとき、信頼度ＲＥＬ_ｉがしきい値Ｒｔｈよりも大きいと判定する。

ステップＳ７７において、信頼度ＲＥＬ_ｉがしきい値Ｒｔｈ以下でないと判定されたとき、推定手段１４は、端末装置２＿ｉを波源位置の推定に用いる端末装置とする（ステップＳ７８）。

一方、ステップＳ７７において、信頼度ＲＥＬ_ｉがしきい値Ｒｔｈ以下であると判定されたとき、推定手段１４は、端末装置２＿ｉを波源位置の推定に用いる端末装置から削除する（ステップＳ７９）。

そして、ステップＳ７８またはステップＳ７９の後、推定手段１４は、ｉ＝Ｎであるか否かを判定する（ステップＳ８０）。

ステップＳ８０において、ｉ＝Ｎでないと判定されたとき、推定手段１４は、ｉ＝ｉ＋１を設定する（ステップＳ８１）。その後、一連の動作は、ステップＳ７２へ移行し、ステップＳ８０において、ｉ＝Ｎであると判定されるまで、ステップＳ７２〜ステップＳ８１が繰り返し実行される。

そして、ステップＳ８０において、ｉ＝Ｎであると判定されると、一連の動作は、図１１のステップＳ２へ移行する。

このように、波源位置推定装置１０は、複数の端末装置２の位置情報（ｘ_ｉ，ｙ_ｉ）および受信電力Ｐ_ｉに基づいて推定した波源位置に誤差がある場合、Ｊ個の仮想端末装置２ｖ１〜２ｖＪを追加して複数の端末装置２およびＪ個の仮想端末装置２ｖ１〜２ｖＪの重心を波源位置として推定する推定処理を収束条件が満たれるまで繰り返し実行する。従って、仮想端末装置２ｖｊを追加せずに波源位置を推定する場合よりも波源位置を正確に推定できる。

図１６は、図２に示す波源位置推定装置１の動作を説明するための別のフローチャートである。なお、図１６に示すフローチャートは、判定回数が所定回数Ｋ（Ｋは正の整数）であるとき、収束したと判定するフローチャートである。

図１６に示すフローチャートは、図１１に示すフローチャートのステップＳ３，Ｓ４をそれぞれステップＳ３Ａ，Ｓ４Ａに変え、ステップＳ４Ｂを追加したものであり、その他は、図１１に示すフローチャートと同じである。

図１６を参照して、波源位置推定装置１の動作が開始されると、上述したステップＳ１が実行される。

そして、ステップＳ１の後、推定手段１４は、ｋ＝１を設定する（ステップＳ３Ａ）。その後、上述したステップＳ２が実行される。ステップＳ２の後、推定手段１４は、ｋ＝Ｋ（Ｋは、判定の完了を示す判定回数）であるか否かを判定する（ステップＳ４Ａ）。

ステップＳ４Ａにおいて、ｋ＝Ｋでないと判定されたとき、推定手段１４は、ｋ＝ｋ＋１を設定する（ステップＳ４Ｂ）。その後、上述したステップＳ５，Ｓ６が順次実行された後、一連の動作は、ステップＳ２へ移行し、ステップＳ４Ａにおいて、ｋ＝Ｋであると判定されるまで、ステップＳ２，Ｓ４Ａ，Ｓ４Ｂ，Ｓ５，Ｓ６が繰り返し実行される。そして、ステップＳ４Ａにおいて、ｋ＝Ｋであると判定されると、一連の動作が終了する。

図１６に示すフローチャートによれば、判定回数が所定回数Ｋに達すると波源位置の推定処理を完了するので、計算量を低減して波源位置を正確に推定できる。

図１７は、図２に示す波源位置推定装置１の動作を説明するための更に別のフローチャートである。なお、図１７に示すフローチャートは、ｋ−１（ｋは２以上の整数）回目の判定時において推定された波源位置と、ｋ回目の判定時において推定された波源位置との誤差ΔＧがしきい値Ｇｔｈ以下であるときに収束したと判定するフローチャートである。

図１７に示すフローチャートは、図１６に示すフローチャートのステップＳ４ＡをステップＳ９１〜ステップＳ９３に変えたものであり、その他は、図１６に示すフローチャートと同じである。

図１７を参照して、波源位置推定装置１の動作が開始されると、上述したステップＳ１，Ｓ３Ａ，Ｓ２が順次実行される。

そして、ステップＳ２の後、推定手段１４は、ｋが２以上であるか否かを判定する（ステップＳ９１）。

ステップＳ９１において、ｋが２以上であると判定されたとき、推定手段１４は、ｋ−１回目の判定時に推定された波源位置と、ｋ回目の判定時に推定された波源位置との誤差ΔＧを演算する（ステップＳ９２）。

その後、推定手段１４は、誤差ΔＧがしきい値Ｇｔｈ以下であるか否かを判定する（ステップＳ９３）。

ステップＳ９１において、ｋが２以上でないと判定されたとき、またはステップＳ９３において、誤差ΔＧがしきい値Ｇｔｈ以下でないと判定されたとき、推定手段１４は、ｋ＝ｋ＋１を設定する（ステップＳ４Ｂ）。そして、上述したステップＳ５，Ｓ６が順次実行された後、一連の動作は、ステップＳ２へ移行し、ステップＳ９３において、誤差ΔＧがしきい値Ｇｔｈ以下であると判定されるまで、ステップＳ２，Ｓ９１〜ステップＳ９３，Ｓ４Ｂ，Ｓ５，Ｓ６が繰り返し実行される。そして、ステップＳ９３において、誤差ΔＧがしきい値Ｇｔｈ以下であると判定されると、一連の動作が終了する。

図１７に示すフローチャートによれば、連続する２回において推定された波源位置の誤差がしきい値Ｇｔｈ以下になるまで波源位置の推定処理が繰り返し実行されるので、推定した波源位置の誤差を低減できる。

図１８は、実施の形態１による別の波源位置推定装置の概略図である。実施の形態１による波源位置推定装置は、図１８に示す波源位置推定装置１Ａであってもよい。

図１８を参照して、波源位置推定装置１Ａは、図２に示す波源位置推定装置１の推定手段１４を推定手段１４Ａに変えたものであり、その他は、波源位置推定装置１と同じである。

推定手段１４Ａは、記憶手段１３に記憶された複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］に含まれる複数の位置情報（ｘ_１，ｙ_１）〜（ｘ_Ｎ，ｙ_Ｎ）および複数の受信電力Ｐ_１〜Ｐ_Ｎに基づいて、後述する方法によって、電力分布を推定し、その推定結果である電力分布推定結果を生成する。そして、推定手段１４Ａは、その生成した電力分布推定結果を用いて上述した方法によって波源位置を推定する。

図１９は、受信電力の等高線ＣＴＲを作成する方法を説明するための図である。図１９において、白四角は、電力値がＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上である測定点を示し、黒四角は、電力値がＰＷ未満である測定点を示し、白丸は、受信電力の重心点Ｐ＿Ｇを示す。また、測定点ｓ１〜ｓ３は、フェージングの影響によって受信電力が電力値ＰＷ未満になった測定点を示す。そして、白四角および黒四角の各々は、位置情報（ｘ_ｉ，ｙ_ｉ）および受信電力Ｐ_ｉからなる。

図１９の（ａ）を参照して、推定手段１４Ａは、記憶手段１３に記録された複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］に基づいて、複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］の分布状態（図１９の（ａ）に示す分布状態）を取得する。

そして、推定手段１４Ａは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷを決定する。推定手段１４Ａは、例えば、−６０ｄＢｍ，−７０ｄＢｍ，−８０ｄＢｍ等を電力値ＰＷとして決定する。

推定手段１４Ａは、電力値ＰＷと処理範囲の半径ｒとの関係を予め保持している。例えば、推定手段１４Ａは、−６０ｄＢｍの電力値ＰＷに対して処理範囲の半径ｒ＝５００ｍ、および−７０ｄＢｍの電力値ＰＷに対して処理範囲の半径１５００ｍ等の関係を予め保持している。

推定手段１４Ａは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷを決定すると、複数の端末情報ＴＩＦ＿１〜ＴＩＦ＿Ｎ＝［ＩＤ_１／（ｘ_１，ｙ_１）／Ｐ_１／ｔ_１］〜［ＩＤ_Ｎ／（ｘ_Ｎ，ｙ_Ｎ）／Ｐ_Ｎ／ｔ_Ｎ］に含まれる複数の受信電力Ｐ_１〜Ｐ_Ｎのうち、電力値ＰＷ以上である受信電力の最大値Ｐ＿ＭＡＸを抽出する。ここで、受信電力の最大値Ｐ＿ＭＡＸを抽出するのは、受信電力の最大値Ｐ＿ＭＡＸを有する端末装置２が最も波源Ｓに近いと考えられるからである。

そして、推定手段１４Ａは、−６０ｄＢｍの電力値ＰＷを有する受信電力の等高線ＣＴＲを作成する場合、受信電力の最大値Ｐ＿ＭＡＸを抽出すると、その抽出した受信電力の最大値Ｐ＿ＭＡＸ（＝ＲＳＳＩ＿ＭＡＸ）を有する測定点を中心として半径ｒ＝５００ｍの範囲を処理範囲ＲＥＧ＿ＰＲＳとして決定する。

そうすると、推定手段１４Ａは、処理範囲ＲＥＧ＿ＰＲＳ内のｄ（ｄは２以上の整数）個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する。

図１９の（ｂ）を参照して、推定手段１４Ａは、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに含まれるｄ個の受信電力Ｐ_１〜Ｐ_ｄに基づいて、式（１）を用いて受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）（Ｐ＿Ｇ＝ＲＳＳＩ＿Ｇ）を算出する。

即ち、推定手段１４Ａは、ｄ個の受信電力Ｐ_１〜Ｐ_ｄによって重み付けされたｄ個の端末装置２の位置の平均を演算することによって受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）を算出する。

そうすると、推定手段１４Ａは、ｄ個の受信電力Ｐ_１〜Ｐ_ｄのうち、電力値ＰＷ以上の受信電力Ｐ_ｉを有し、かつ、受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）から最も遠い位置に存在する端末装置２−１と受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）との間の距離を最大半径Ｒとして求める。そして、推定手段１４Ａは、受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）を中心とし、かつ、最大半径Ｒを半径とする円形状を有する受信電力の等高線ＣＴＲを作成する。

このように、推定手段１４Ａは、電力値ＰＷ未満の受信電力（黒四角）を除外し、電力値ＰＷ以上の受信電力（白四角）のみに基づいて最大半径Ｒを決定する。黒四角によって示される測定点のうち、測定点ｓ１〜ｓ３は、フェージングの影響によって受信電力が電力値ＰＷ未満になった測定点である。

従って、推定手段１４Ａによる受信電力の等高線ＣＴＲの作成方法によれば、フェージングの影響を抑制して受信電力の等高線ＣＴＲを作成できる。

上記においては、ｄは、２以上の整数であると説明したが、これは、次の理由による。ｄ＝２である場合、２個の受信電力Ｐ_ｉのうち、１個は、受信電力の最大値Ｐ＿ＭＡＸであり、もう１個は、電力値ＰＷ以上である受信電力か電力値ＰＷ未満である受信電力かが不明である受信電力である。そして、この２個の受信電力Ｐ_ｉを用いて、受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）が式（１）によって算出される。

その結果、受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）は、受信電力の最大値Ｐ＿ＭＡＸを有する測定点からずれる。

そうすると、受信電力の最大値Ｐ＿ＭＡＸは、電力値ＰＷ以上であるので、受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）と受信電力の最大値Ｐ＿ＭＡＸを有する測定点との距離を最大半径Ｒとして求めることができ、受信電力の等高線ＣＴＲを作成することができる。従って、ｄは、２以上であればよい。

なお、電力値ＰＷ以上の受信電力Ｐ_ｉが少なくとも１個有れば、少なくとも１個の受信電力Ｐ_ｉを有する端末装置２の位置が受信電力の重心点Ｐ＿Ｇからずれているので、最大半径Ｒを決定でき、受信電力の等高線ＣＴＲを作成することができる。従って、受信電力の等高線ＣＴＲを作成するために必要な電力値ＰＷ以上の受信電力の個数は、１個以上である。

図２０は、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄの抽出方法を説明するための図である。図２０において、黒丸は、各測定点（端末情報ＴＩＦ）を示し、電力値ＰＷ以上であるか電力値ＰＷ未満であるかを識別していない。

図２０の（ａ）を参照して、推定手段１４Ａは、上述した方法によって、処理範囲ＲＥＧ＿ＰＲＳ＿１を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿１内に存在するｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する。そして、推定手段１４Ａは、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて、上述した方法によって、受信電力の等高線ＣＴＲ１を作成する。

図２０の（ｂ）を参照して、推定手段１４Ａは、その後、処理範囲ＲＥＧ＿ＰＲＳ＿１内に存在するｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを除外して、上述した方法によって、処理範囲ＲＥＧ＿ＰＲＳ＿２を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿２内に存在する新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する。そして、推定手段１４Ａは、新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて、上述した方法によって、受信電力の等高線ＣＴＲ２を作成する。

そして、推定手段１４Ａは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力がなくなるまで、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを除外しながら、上述した方法によって処理範囲ＲＥＧ＿ＰＲＳを決定し、その決定した処理範囲ＲＥＧ＿ＰＲＳ内に存在する新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出することを繰り返し実行する。

そして、推定手段１４Ａは、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出するごとに、上述した方法によって、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の等高線ＣＴＲを作成する。

この場合、推定手段１４Ａは、次の２つの方法のいずれかを用いて受信電力の等高線ＣＴＲを作成する。

（１）作成方法１
推定手段１４Ａは、既に受信電力の等高線ＣＴＲの作成に用いたｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを用いずに、抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の等高線ＣＴＲを作成する。

（２）作成方法２
推定手段１４Ａは、既に受信電力の等高線ＣＴＲの作成に用いたｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄと新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄとの両方に含まれる端末情報と、新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄとに基づいて受信電力の等高線ＣＴＲを作成する。

図２１は、受信電力の等高線ＣＴＲの作成方法２を説明するための図である。図２１を参照して、端末情報ＴＩＦ＿１は、受信電力の最大値Ｐ＿ＭＡＸ１を含み、端末情報ＴＩＦ＿２は、端末情報ＴＩＦ＿１を除いた場合における受信電力の最大値Ｐ＿ＭＡＸ２を含む。また、端末情報ＴＩＦ＿３，ＴＩＦ＿４は、電力値ＰＷ以上の受信電力を含む。更に、白丸は、受信電力の重心点Ｐ＿Ｇを示す。

推定手段１４Ａは、処理範囲ＲＥＧ＿ＰＲＳ＿３を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿３内に存在するｄ個の端末情報ＴＩＦ＿ＲＥＧ＿ＰＲＳ＿３を抽出する。

そして、推定手段１４Ａは、ｄ個の端末情報ＴＩＦ＿ＲＥＧ＿ＰＲＳ＿３に基づいて、上述した方法によって受信電力の等高線ＣＴＲ１を作成する。

その後、推定手段１４Ａは、処理範囲ＲＥＧ＿ＰＲＳ＿３内に存在するｄ個の端末情報ＴＩＦ＿ＲＥＧ＿ＰＲＳ＿３を除いて、処理範囲ＲＥＧ＿ＰＲＳ＿４を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿４内に存在するｄ個の端末情報ＴＩＦＴ＿ＲＥＧ＿ＰＲＳ＿４を抽出する。

この場合、端末情報ＴＩＦ＿３，ＴＩＦ＿４は、ｄ個の端末情報ＭＮＴ＿ＲＥＧ＿ＰＲＳ＿３およびｄ個の端末情報ＭＮＴ＿ＲＥＧ＿ＰＲＳ＿４の両方に含まれる。

そして、推定手段１４Ａは、ｄ個の端末情報ＭＮＴ＿ＲＥＧ＿ＰＲＳ＿４および端末情報ＴＩＦ＿３，ＴＩＦ＿４に基づいて受信電力の重心点Ｐ＿Ｇ１を算出し、受信電力の重心点Ｐ＿Ｇ１と端末情報ＴＩＦ＿４との距離を半径とする円形状を有する受信電力の等高線ＣＴＲ２を作成する。

端末情報ＴＩＦ＿３，ＴＩＦ＿４を含めずに受信電力の等高線ＣＴＲを作成する場合、受信電力の重心点Ｐ＿Ｇは、受信電力の重心点Ｐ＿Ｇ２になり、受信電力の最大値Ｐ＿ＭＡＸ２を含む端末情報ＴＩＦ＿２の位置からのずれが大きくなり、受信電力の等高線ＣＴＲ２よりも半径が小さい円形状を有する受信電力の等高線ＣＴＲになる。その結果、精度良く受信電力の等高線ＣＴＲを作成することが困難になる。

しかし、端末情報ＴＩＦ＿３，ＴＩＦ＿４を含めて受信電力の等高線ＣＴＲを作成することにより、受信電力の重心点Ｐ＿Ｇ１は、受信電力の最大値Ｐ＿ＭＡＸ２を含む端末情報ＴＩＦ＿２の位置の近傍に存在し、より大きい半径を有する受信電力の等高線ＣＴＲ２が作成される。その結果、精度良く受信電力の等高線ＣＴＲを作成することができる。

図２２は、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄの別の抽出方法を説明するための図である。図２２を参照して、受信電力の最大値Ｐ＿ＭＡＸ６は、受信電力の最大値Ｐ＿ＭＡＸ５よりも小さい。

推定手段１４Ａは、上述した方法によって、処理範囲ＲＥＧ＿ＰＲＳ＿５を決定し、処理範囲ＲＧ＿ＰＲＳ＿５内に存在するｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する（図２２の（ａ）参照）。

その後、推定手段１４Ａは、既に抽出した処理範囲ＲＥＧ＿ＰＲＳ＿５内のｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸ５を含む端末情報ＴＩＦのみを除外して処理範囲ＲＥＧ＿ＰＲＳ＿６を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿６内に存在するｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する（図２２の（ｂ）参照）。

このように、推定手段１４Ａは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力がなくなるまで、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸを含む端末情報ＴＩＦのみを除外しながら新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを繰り返し抽出する。

そして、推定手段１４Ａは、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出するごとに、その抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて、上述した方法によって受信電力の等高線ＣＴＲを作成する。

上述したように、推定手段１４Ａは、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸを含む端末情報ＴＩＦを少なくとも除きながらｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出することを、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力がなくなるまで繰り返し行う。

図２３は、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄの抽出方式と受信電力の等高線ＣＴＲの作成処理との関係を示す図である。

図２３を参照して、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄの抽出方式は、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力がなくなるまで、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを除外しながら新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する抽出方式１と、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力がなくなるまで、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸを含む端末情報ＴＩＦのみを除外しながら新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する抽出方式２とを含む。

また、受信電力の等高線ＣＴＲの作成処理は、抽出されたｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の重心点Ｐ＿Ｇを算出し、電力値ＰＷ以上の受信電力を有し、かつ、受信電力の重心点Ｐ＿Ｇから最も遠い位置に存在する端末装置と受信電力の重心点Ｐ＿Ｇとの距離を最大距離として求め、受信電力の重心点Ｐ＿Ｇを中心とし、かつ、最大距離を半径とする円形状を有する受信電力の等高線ＣＴＲを作成する等高線作成処理１と、既に抽出されたｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄと新たに抽出されたｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄとの両方に含まれる端末情報と、新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄとに基づいて受信電力の重心点Ｐ＿Ｇを算出し、電力値ＰＷ以上の受信電力を有し、かつ、受信電力の重心点Ｐ＿Ｇから最も遠い位置に存在する端末装置と受信電力の重心点Ｐ＿Ｇとの距離を最大距離として求め、受信電力の重心点Ｐ＿Ｇを中心とし、かつ、最大距離を半径とする円形状を有する受信電力の等高線ＣＴＲを作成する等高線作成処理２とを含む。

そして、推定手段１４Ａは、抽出方式１に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する場合、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成し、または等高線作成処理２に従って受信電力の等高線ＣＴＲを作成する。

また、推定手段１４Ａは、抽出方式２に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する場合、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する。

ここで、推定手段１４Ａが、抽出方式１に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する場合、計算量が最も少なくなり、推定手段１４Ａが、抽出方式１に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理２に従って受信電力の等高線ＣＴＲを作成する場合、計算量が２番目に少なくなり、推定手段１４Ａが、抽出方式２に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する場合、計算量が最も多くなる。

推定手段１４Ａが、抽出方式１に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する場合、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを除外しながら新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、その抽出した新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の等高線ＣＴＲを作成するので、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉの個数をＩ（Ｉは、１≦Ｉ≦Ｎを満たす整数、Ｎは、Ｎ≧ｄを満たす）個とした場合、受信電力の等高線ＣＴＲを作成する回数は、ＩＮＴ［Ｉ／ｄ］（Ｉ／ｄを整数化したもの）となる。なお、複数の端末装置２の全てが受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉを有する場合もあり得るので、Ｉは、Ｎ以下である。

推定手段１４Ａが、抽出方式１に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理２に従って受信電力の等高線ＣＴＲを作成する場合、既に抽出されたｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄと新たに抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄとの両方に含まれる端末情報と、新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄとに基づいて受信電力の等高線ＣＴＲを作成するので、受信電力の等高線ＣＴＲを作成する回数は、ＩＮＴ［Ｉ／ｄ］となるが、受信電力の重心点Ｐ＿Ｇを算出する際の端末情報（受信電力Ｐ_ｉおよび位置情報（ｘ_ｉ，ｙ_ｉ））の個数および最大距離を求める際の端末情報（受信電力Ｐ_ｉおよび位置情報（ｘ_ｉ，ｙ_ｉ））の個数が等高線作成処理１の場合よりも多くなり、計算量は、等高線作成処理１の場合よりも多くなる。

推定手段１４Ａが、抽出方式２に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する場合、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸを含む端末情報ＴＩＦのみを除外しながら新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、その抽出した新たなｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の等高線ＣＴＲを作成するので、受信電力の等高線ＣＴＲを作成する回数は、Ｉ回となる。ｄ＝２であり、かつ、Ｉ＝１である場合、推定手段１４Ａが、抽出方式２によってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出した場合の受信電力の等高線ＣＴＲの作成回数Ｉは、推定手段１４Ａが抽出方式１によってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出した場合の受信電力の等高線ＣＴＲの作成回数ＩＮＴ［Ｉ／ｄ］と等しくなる（等高線ＣＴＲの作成回数＝１回）。一方、Ｉ≧２である場合、推定手段１４Ａが抽出方式２によってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出した場合の受信電力の等高線ＣＴＲの作成回数Ｉは、推定手段１４Ａが抽出方式１によってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出した場合の受信電力の等高線ＣＴＲの作成回数ＩＮＴ［Ｉ／ｄ］よりも大きくなる。この発明の実施の形態において、推定手段１４Ａが受信電力の等高線ＣＴＲを作成する対象領域ＲＥＧは、市町村の領域または県の領域であるので、通常、Ｉ＝１である状況は、まず生じ得ないと考えられる。その結果、推定手段１４Ａが抽出方式２によってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出した場合の受信電力の等高線ＣＴＲの作成回数Ｉは、通常、最も多くなる。従って、推定手段１４Ａが、抽出方式２に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する場合、計算量が最も多くなる。

そして、推定手段１４Ａが、抽出方式１に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する方式を第１の方式とし、推定手段１４Ａが、抽出方式１に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理２に従って受信電力の等高線ＣＴＲを作成する方式を第２の方式とし、推定手段１４Ａが、抽出方式２に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出し、かつ、等高線作成処理１に従って受信電力の等高線ＣＴＲを作成する方式を第３の方式とした場合、受信電力の等高線ＣＴＲは、計算量を最少に設定する場合、第１の方式を用いて作成され、計算量が２番目に少なくなるように設定する場合、第２の方式を用いて作成され、計算量が最大の計算量になることが許容される場合、第３の方式を用いて作成される。

図２４は、受信電力の等高線ＣＴＲを作成する別の方法を説明するための図である。図２４において、白四角は、電力値がＰＷ以上である測定点を示し、黒四角は、電力値がＰＷ未満である測定点を示し、白丸は、受信電力の重心点Ｐ＿Ｇを示す。また、白四角および黒四角の各々は、位置情報（ｘ_ｉ，ｙ_ｉ）および受信電力Ｐ_ｉからなる。

推定手段１４Ａは、上述した方法によって、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の等高線ＣＴＲ３，ＣＴＲ４を作成する（２４の（ａ）参照）。

そして、推定手段１４Ａは、２つの受信電力の等高線ＣＴＲ３，ＣＴＲ４が相互に交差する場合、受信電力の等高線ＣＴＲ３，ＣＴＲ４を合成して受信電力の等高線ＣＴＲ＿ＳＹＮを作成する（図２４の（ｂ）参照）。その結果、受信電力の等高線ＣＴＲ＿ＳＹＮは、雪だるまの形状を有する。

このように、推定手段１４Ａは、波源Ｓが複数ある場合、受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）を波源群（複数の波源Ｓの集合）の仮想波源と見做して受信電力の等高線ＣＴＲを作成する。

従来の技術においては、波源ごとに受信電力の等高線を作成するため、波源の数が多くなると、計算量が多くなる。

一方、この発明の実施の形態においては、波源群ごとに受信電力の等高線ＣＴＲを作成するため、複数の波源があっても、波源位置推定装置１の処理量は、それほど増加しないため、従来の技術に比べて処理量を少なくできる。

図２５は、端末装置２における電波の受信感度以下の領域において受信電力を推定する方法を説明するための図である。

端末装置２は、受信感度以下の受信電力を取得できないため、受信感度以下の領域における受信電力の等高線ＣＴＲは、既に作成した受信電力の等高線ＣＴＲと電波の伝搬モデルとに基づいて作成する。

電波の自由空間伝搬モデルでは、受信電力Ｐ_ｉは、図３に示すように波源Ｓからの距離に反比例して減少する。

電波の自由空間伝搬モデルにおいて、波源Ｓからの距離が１０００ｍ以上である場合、受信電力Ｐ_ｉは、距離の増加に対して直線状に低下する。そして、受信電力Ｐ_ｉが１０ｄＢｍだけ低下するときの距離は、３．１６ｋｍである。つまり、電波の自由空間伝搬モデルにおいては、受信電力Ｐ_ｉが１０ｄＢｍだけ低下するときの距離は、３．１６ｋｍ／１０ｄＢｍとなる。

従って、推定手段１４Ａは、電波の受信感度以下の領域において、既に作成した受信電力の等高線ＣＴＲ＿ＳＹＮ（電力値ＰＷの等高線）から３．１６ｋｍだけ波源から遠くなる方向に離れ、かつ、受信電力の等高線ＣＴＲ＿ＳＹＮと相似形を有する受信電力の等高線ＣＴＲ＿ＬＳをＰＷ−１０[ｄＢ]の電力値を有する等高線ＣＴＲとして作成する（図２５参照）。

この場合、受信電力の等高線ＣＴＲ＿ＬＳにおける電力値は、ＰＷ−１０[ｄＢ]に限らず、電力値ＰＷから所望の電力値ΔＰＷだけ低下した電力値であればよい。従って、推定手段１４Ａは、電波の自由空間伝搬モデルにおいて、電力値ＰＷよりもΔＰＷだけ低下するときの距離を求め、その求めた距離だけ波源から遠くなる方向に離れ、かつ、受信電力の等高線ＣＴＲ＿ＳＹＮと相似形を有する受信電力の等高線ＣＴＲ＿ＬＳをＰＷ−ΔＰＷの電力値を有する等高線ＣＴＲとして作成する。

また、電波の伝搬モデルは、自由空間伝搬モデルに限らず、どのような伝搬モデルであってもよい。

従って、推定手段１４Ａは、一般的に、電波の受信感度以下の領域において、既に作成した受信電力の等高線ＣＴＲ上の電力値よりも所望の電力値だけ低い受信電力になる電力低下距離（３．１６ｋｍ等）を電波の伝搬モデルに基づいて求め、既に作成した受信電力の等高線ＣＴＲから電力低下距離だけ離れ、かつ、既に作成した受信電力の等高線ＣＴＲと相似形を有する等高線ＣＴＲを電波の受信感度以下の領域における受信電力の等高線ＣＴＲとして作成する。

このように、推定手段１４Ａは、電波の受信感度以下の領域においても受信電力の等高線ＣＴＲを作成するので、受信電力Ｐ_ｉの広い範囲において受信電力の等高線ＣＴＲを作成できる。

図２６は、電力分布推定結果を作成する方法を説明するためのフローチャートである。なお、この発明の実施の形態においては、等高線ＣＴＲを作成することは、電力分布推定結果を作成することに相当する。

図２６を参照して、電力分布推定結果の作成が開始されると、推定手段１４Ａは、記憶手段１３に記録された複数の端末情報ＴＩＦに基づいて、上述した方法によって処理範囲ＲＥＧ＿ＰＲＳ内に存在するｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する（ステップＳ１０１）。

そして、推定手段１４Ａは、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを除外し、または既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸを有する端末情報ＴＩＦを除外する（ステップＳ１０２）。

そして、推定手段１４Ａは、既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを除外した複数の端末情報、または既に抽出したｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸを有する端末情報ＴＩＦを除外した複数の端末情報に基づいて、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在しないか否かを判定する（ステップＳ１０３）。

ステップＳ１０３において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在すると判定されたとき、一連の動作は、ステップＳ１０１へ戻り、ステップＳ１０３において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在しないと判定されるまで、ステップＳ１０１〜ステップＳ１０３が繰り返し実行される。

そして、ステップＳ１０３において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在しないと判定されると、推定手段１４Ａは、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出するごとに、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて、上述した方法によって受信電力の等高線ＣＴＲを作成する（ステップＳ１０４）。これによって、電力分布推定結果の作成が終了する。

なお、推定手段１４Ａは、ステップＳ１０１において、処理範囲ＲＥＧ＿ＰＲＳ内に電力値が同じである複数の受信電力の最大値Ｐ＿ＭＡＸが存在する場合、複数の受信電力の最大値Ｐ＿ＭＡＸのうちの１つの受信電力の最大値Ｐ＿ＭＡＸを選択し、それ以外の受信電力の最大値Ｐ＿ＭＡＸを除外して処理範囲ＲＥＧ＿ＰＲＳ内のｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを選択する。

この発明の実施の形態においては、推定手段１４Ａは、抽出方式１または抽出方式２に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出するので、受信電力の最大値Ｐ＿ＭＡＸを少なくとも除外しながらｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを繰り返し抽出する。

従って、推定手段１４Ａは、ステップＳ１０１において、処理範囲ＲＥＧ＿ＰＲＳ内に電力値が同じである複数の受信電力の最大値Ｐ＿ＭＡＸ（等高線ＣＴＲの電力値ＰＷ以上の受信電力の最大値Ｐ＿ＭＡＸ）が存在する場合、複数の受信電力の最大値Ｐ＿ＭＡＸのうちの１つの受信電力の最大値Ｐ＿ＭＡＸを選択し、それ以外の受信電力の最大値Ｐ＿ＭＡＸを除外して処理範囲ＲＥＧ＿ＰＲＳ内のｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを選択することによって、複数の受信電力の最大値Ｐ＿ＭＡＸの各々を含むｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の等高線ＣＴＲを作成する。

その結果、処理範囲ＲＥＧ＿ＰＲＳ内に電力値が同じである複数の受信電力の最大値Ｐ＿ＭＡＸが存在する場合でも、受信電力の等高線ＣＴＲを正確に作成できる。

図２７は、図２６のステップＳ１０１の詳細な動作を説明するためのフローチャートである。

図２７を参照して、図２６の”スタート”の後、または図２６のステップＳ１０３において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在すると判定されたとき、推定手段１４Ａは、複数の端末情報ＴＩＦに基づいて受信電力の最大値ＰＩ＿ＭＡＸを求める（ステップＳ１０１１）。

そして、推定手段１４Ａは、受信電力の最大値Ｐ＿ＭＡＸが、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上であるか否かを判定する（ステップＳ１０１２）。なお、ステップＳ１０１１が、図２６のステップＳ１０３において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在すると判定された後に実行される場合、受信電力の最大値Ｐ＿ＭＡＸは、通常、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上になるが、この場合、ステップＳ１０１２は、より正確に受信電力の等高線ＣＴＲを作成するために実行される。

ステップＳ１０１２において、受信電力の最大値Ｐ＿ＭＡＸが、作成すべき等高線ＣＴＲの電力値ＰＷ以上でないと判定されたとき、一連の動作は、図２６の”終了”へ移行する。

一方、ステップＳ１０１２において、受信電力の最大値Ｐ＿ＭＡＸが、作成すべき等高線ＣＴＲの電力値ＰＷ以上であると判定されたとき、推定手段１４Ａは、受信電力の最大値Ｐ＿ＭＡＸを中心とした半径ｒの処理範囲ＲＥＧ＿ＰＲＳを決定する（ステップＳ１０１３）。

そして、推定手段１４Ａは、処理範囲ＲＥＧ＿ＰＲＳ内に存在するｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する（ステップＳ１０１４）。

その後、一連の動作は、図２６のステップＳ１０２へ移行する。

なお、図２７に示すフローチャートが、図２６のステップＳ１０３において作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在すると判定された後に実行される場合、推定手段１４Ａは、図２３に示す抽出方式１または抽出方式２に従ってｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄを抽出する。

図２８は、図２６のステップＳ１０４の詳細な動作を説明するためのフローチャートである。

図２８を参照して、図２６のステップＳ１０３において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上の受信電力Ｐ_ｉが存在しないと判定されたとき、推定手段１４Ａは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷが端末装置２の受信感度以上であるか否かを判定する（ステップＳ１０４１）。

ステップＳ１０４１において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷが端末装置２の受信感度以上でないと判定されたとき、推定手段１４Ａは、既に作成した受信電力の等高線ＣＴＲが有るか否かを更に判定する（ステップＳ１０４２）。

ステップＳ１０４２において、既に作成した受信電力の等高線ＣＴＲが無いと判定されたとき、推定手段１４Ａは、新たなＰＷを選択する（ステップＳ１０４３）。その後、一連の動作は、ステップＳ１０４１へ移行する。

一方、ステップＳ１０４２において、既に作成した受信電力の等高線ＣＴＲが有ると判定されたとき、推定手段１４Ａは、既に作成した受信電力の等高線ＣＴＲと電波の伝搬モデルとを用いて、上述した方法によって、電波の受信感度以下の領域において受信電力の等高線ＣＴＲを作成する（ステップＳ１０４４）。

ステップＳ１０４１において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷが端末装置２の受信感度以上であると判定されたとき、推定手段１４Ａは、ｄ個の端末情報ＴＩＦ＿１〜ＴＩＦ＿ｄに基づいて受信電力の重心点Ｐ＿Ｇを算出する（ステップＳ１０４５）。

そして、推定手段１４Ａは、電力値ＰＷ以上の受信電力を有し、かつ、受信電力の重心点Ｐ＿Ｇから最も遠い位置に存在する端末装置と受信電力の重心点Ｐ＿Ｇとの距離を最大距離として求める（ステップＳ１０４６）。

そうすると、推定手段１４Ａは、受信電力の重心点Ｐ＿Ｇを中心とし、最大距離を半径とする円形状を有する受信電力の等高線ＣＴＲを作成する（ステップＳ１０４７）。

そして、ステップＳ１０４４またはステップＳ１０４７の後、一連の動作は、図２６の“終了”へ移行する。

図２９は、図１８に示す波源位置推定装置１Ａの動作を説明するためのフローチャートである。

図２９に示すフローチャートは、図１１に示すフローチャートにステップＳ１１０を追加したものであり、その他は、図１１に示すフローチャートと同じである。

図２９を参照して、波源位置推定装置１Ａの動作が開始されると、上述したステップＳ１〜ステップＳ４が順次実行される。

そして、ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａでないと判定されたとき、推定手段１４Ａは、図２６に示すフローチャート（図２７および図２８に示すフローチャートを含む）に従って電力分布推定結果を作成する（ステップＳ１１０）。

その後、推定手段１４Ａは、電力分布推定結果を用いて、上述したステップＳ５，Ｓ６を順次実行する。そして、ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａであると判定されたとき、波源位置推定装置１Ａの動作が終了する。

図３０は、図１８に示す波源位置推定装置１Ａの動作を説明するための別のフローチャートである。

図３０に示すフローチャートは、図１６に示すフローチャートにステップＳ１１０を追加したものであり、その他は、図１６に示すフローチャートと同じである。

図３０を参照して、波源位置推定装置１Ａの動作が開始されると、上述したステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ４Ａ，Ｓ４Ｂが順次実行される。

そして、ステップＳ４Ｂの後、推定手段１４Ａは、上述したステップＳ１１０，Ｓ５，Ｓ６を順次実行し、ステップＳ４Ａにおいて、ｋ＝Ｋであると判定されると、波源位置推定装置１Ａの動作が終了する。

図３１は、図１８に示す波源位置推定装置１Ａの動作を説明するための更に別のフローチャートである。

図３１に示すフローチャートは、図１７に示すフローチャートにステップＳ１１０を追加したものであり、その他は、図１７に示すフローチャートと同じである。

図３１を参照して、波源位置推定装置１Ａの動作が開始されると、上述したステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ９１〜Ｓ９３，Ｓ４Ｂが順次実行される。

そして、ステップＳ４Ｂの後、推定手段１４Ａは、上述したステップＳ１１０，Ｓ５，Ｓ６を順次実行し、ステップＳ９３において、誤差ΔＧがしきい値Ｇｔｈ以下であると判定されると、波源位置推定装置１Ａの動作が終了する。

このように、波源位置推定装置１Ａは、波源位置の推定が収束条件を満たさないと判定された場合、収束条件を満たすと判定されるまで、電力分布を推定し、その推定した電力分布を用いて仮想端末装置２ｖｊの追加と受信電力Ｐ_ｉの信頼度ＲＥＬ_ｉが低い端末装置２の削除とを行い、波源位置を推定する処理を繰り返し行う。

図３２は、波源位置の推定誤差と端末装置の個数との関係を示す図である。図３２において、縦軸が、波源位置の推定誤差を表し、横軸は、端末装置の個数を表す。そして、黒丸および破線は、実端末装置である端末装置２の位置情報（ｘ_ｉ，ｙ_ｉ）および受信電力Ｐ_ｉのみを用いて波源位置を推定したとき（従来方式）の波源位置の推定誤差と端末装置の個数との関係を示す。また、白菱形および実線は、この発明の実施の形態による方法を用いて波源位置を推定したときの波源位置の推定誤差と端末装置の個数との関係を示す。

なお、図３２においては、端末装置２の個数を“１０”に固定して仮想端末装置２ｖｊの個数を１０〜９０の範囲で変化させたときの波源位置の推定誤差と端末装置の個数との関係を示す。従って、図３２の横軸は、端末装置２の個数と仮想端末装置２ｖｊの個数との和である。

図３２に示すように、従来方式およびこの発明の実施の形態における方式とも、波源位置の推定誤差は、端末装置の個数（総数）に対して同様に変化することが分かった。

なお、この発明の実施の形態においては、波源位置推定装置１の動作は、ソフトウェアによって実行されてもよい。

この場合、波源位置推定装置１は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）およびＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）を備える。

ＲＯＭは、図１１に示すフローチャートのステップＳ１〜Ｓ６（図１２〜図１４のいずれかに示すフローチャートと図１５に示すフローチャートとを含む）を備えるプログラムＰｒｏｇ＿Ａ、または図１６に示すフローチャートのステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ４Ａ，Ｓ４Ｂ，Ｓ５，Ｓ６（図１２〜図１４のいずれかに示すフローチャートと図１５に示すフローチャートとを含む）を備えるプログラムＰｒｏｇ＿Ｂ、または図１７に示すフローチャートのステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ９１〜Ｓ９３，Ｓ４Ｂ，Ｓ５，Ｓ６（図１２〜図１４のいずれかに示すフローチャートと図１５に示すフローチャートとを含む）を備えるプログラムＰｒｏｇ＿Ｃを格納する。ＲＡＭは、各種演算の結果を一時的に格納する。

そして、ＣＰＵは、ＲＯＭからプログラムＰｒｏｇ＿Ａ〜Ｐｒｏｇ＿Ｃのいずれかを読み出して実行し、上述した方法によって、波源位置を推定する。

また、この発明の実施の形態においては、波源位置推定装置１Ａの動作は、ソフトウェアによって実行されてもよい。

この場合、波源位置推定装置１Ａは、ＣＰＵ、ＲＯＭおよびＲＡＭを備える。

ＲＯＭは、図２９に示すフローチャートのステップＳ１〜Ｓ４，Ｓ１１０，Ｓ５，Ｓ６（図１２〜図１４のいずれかに示すフローチャートと図１５に示すフローチャートと図２６〜図２８に示すフローチャートとを含む）を備えるプログラムＰｒｏｇ＿Ｄ、または図３０に示すフローチャートのステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ４Ａ，Ｓ４Ｂ，Ｓ１１０，Ｓ５，Ｓ６（図１２〜図１４のいずれかに示すフローチャートと図１５に示すフローチャートと図２６〜図２８に示すフローチャートとを含む）を備えるプログラムＰｒｏｇ＿Ｅ、または図３１に示すフローチャートのステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ９１〜Ｓ９３，Ｓ４Ｂ，Ｓ１１０，Ｓ５，Ｓ６（図１２〜図１４のいずれかに示すフローチャートと図１５に示すフローチャートと図２６〜図２８に示すフローチャートとを含む）を備えるプログラムＰｒｏｇ＿Ｆを格納する。ＲＡＭは、各種演算の結果を一時的に格納する。

そして、ＣＰＵは、ＲＯＭからプログラムＰｒｏｇ＿Ｄ〜Ｐｒｏｇ＿Ｆのいずれかを読み出して実行し、上述した方法によって、電力分布を推定し、その推定した電力分布を用いて波源位置を推定する。

また、プログラムＰｒｏｇ＿Ａ〜Ｐｒｏｇ＿Ｆは、ＣＤ，ＤＶＤ等の記録媒体に記録されて流通してもよい。

プログラムＰｒｏｇ＿Ａ〜Ｐｒｏｇ＿Ｆのいずれかを記録した記録媒体をパーソナルコンピュータに装着すると、ＣＰＵは、記録媒体からプログラムＰｒｏｇ＿Ａ〜Ｐｒｏｇ＿Ｆのいずれかを読み出して実行し、上述した方法によって、波源位置を推定する。

従って、プログラムＰｒｏｇ＿Ａ〜Ｐｒｏｇ＿Ｆのいずれかを記録した記録媒体は、コンピュータ（ＣＰＵ）読み取り可能な記録媒体である。

なお、上述した波源位置の推定方法は、説明上の領域ＲＥＧを簡易的に設定し、その領域ＲＥＧの中心を原点とする座標系を設定し、その設定した座標系で表される位置を実際に位置情報として取得される緯度経度に変換することを前提としている。

また、この発明の実施の形態においては、信頼度ＲＥＬ_ｉがしきい値Ｒｔｈ以下である端末装置２を削除することを行わなくてもよい。仮想端末装置２ｖｊを追加すれば、複数の端末装置２が偏っていても、式（１）を用いて推定した波源位置の誤差を低減できるからである。

［実施の形態２］
図３３は、実施の形態２における無線通信システムを示す図である。図３３は、無線通信システム１０Ａは、波源位置推定装置１Ｂと、複数の端末装置２２と、複数の基地局３１，３２とを備える。

波源位置推定装置１Ｂ、複数の端末装置２２および複数の基地局３１，３２は、無線通信空間に配置される。波源Ｓおよび観測エリアＲＥＧについては、実施の形態１において説明したとおりである。

複数の端末装置２２の各々は、移動端末または静止端末からなり、実際に存在する端末装置である。そして、複数の端末装置２２の各々は、２次利用者である。基地局３１，３２の各々は、実際に存在する通信装置である。そして、波源Ｓおよび基地局３１，３２の各々は、１次利用者である。

波源位置推定装置１Ｂは、複数の端末装置２２から一定時間におけるサンプルデータＤ_ＳＰＬ＿ｉを受信する。

サンプルデータＤ_ＳＰＬ＿ｉは、端末装置（複数の端末装置２２のいずれか）の識別情報ＩＤ_ｉと、Ｊ（Ｊは、２以上の整数）個の受信情報ＩＮＦＯ_Ｒ＿１〜ＩＮＦＯ_Ｒ＿Ｊと、Ｊ個の送信情報ＩＮＦＯ_Ｔ＿１〜ＩＮＦＯ_Ｔ＿Ｊとを含む。

受信情報ＩＮＦＯ_Ｒ＿ｊ（ｊは、２≦ｊ≦Ｊを満たす整数）は、送信元の識別情報ＩＤ_ＳＲＣ＿ｊと、受信時刻ｔ_Ｒｉ＿ｊと、受信信号Ｓ_Ｒｉ＿ｊと、受信信号Ｓ_Ｒｉ＿ｊの受信電力Ｐ_ｉ＿ｊと、受信信号Ｓ_Ｒｉ＿ｊの周波数ｆ_Ｒｉ＿ｊと、識別情報ＩＤ_ｉを有する端末装置（複数の端末装置２２のいずれか）の位置情報［ｘ_ｉ，ｙ_ｉ］とを含む。

送信情報ＩＮＦＯ_Ｔ＿ｊは、送信先の識別情報ＩＤ_ＤＥＳ＿ｊと、送信時刻ｔ_Ｔｉ＿ｊと、送信信号Ｓ_Ｔｉ＿ｊと、送信信号Ｓ_Ｔｉ＿ｊの周波数ｆ_Ｔｉ＿ｊと、識別情報ＩＤ_ｉを有する端末装置（複数の端末装置２２のいずれか）の位置情報［ｘ_ｉ，ｙ_ｉ］とを含む。

従って、サンプルデータＤ_ＳＰＬ＿ｉは、Ｄ_ＳＰＬ＿ｉ＝［ＩＤ_ｉ／{ＩＤ_ＳＲＣｉ＿１／ｔ_Ｒｉ＿１／Ｓ_Ｒｉ＿１／Ｐ_ｉ＿１／［ｘ_ｉ，ｙ_ｉ］／ｆ_Ｒｉ＿１}，{ＩＤ_ＳＲＣｉ＿２／ｔ_Ｒｉ＿２／Ｓ_Ｒｉ＿２／Ｐ_ｉ＿２／［ｘ_ｉ，ｙ_ｉ］／ｆ_Ｒｉ＿２}，・・・，{ＩＤ_ＳＲＣｉ＿Ｊ／ｔ_Ｒｉ＿Ｊ／Ｓ_Ｒｉ＿Ｊ／Ｐ_ｉ＿Ｊ／［ｘ_ｉ，ｙ_ｉ］／ｆ_Ｒｉ＿Ｊ}／{ＩＤ_ＤＥＳｉ＿１／ｔ_Ｔｉ＿１／Ｓ_Ｔｉ＿１／［ｘ_ｉ，ｙ_ｉ］／ｆ_Ｔｉ＿１}，{ＩＤ_ＤＥＳｉ＿２／ｔ_Ｔｉ＿２／Ｓ_Ｔｉ＿２／［ｘ_ｉ，ｙ_ｉ］／ｆ_Ｔｉ＿２}，・・・，{ＩＤ_ＤＥＳｉ＿Ｊ／ｔ_Ｔｉ＿Ｊ／Ｓ_Ｔｉ＿Ｊ／［ｘ_ｉ，ｙ_ｉ］／ｆ_Ｔｉ＿Ｊ}］からなる。

波源位置推定装置１Ｂは、サンプルデータＤ_ＳＰＬ＿ｉに基づいて、後述する方法によって、１次利用者（波源Ｓおよび基地局３３１，３２の少なくとも１つ。以下、同じ。）と、２次利用者（複数の端末装置２２の少なくとも１つ。以下、同じ。）との通信方式または変調方式を推定する。

そして、波源位置推定装置１Ｂは、推定した通信方式または変調方式に適合するように複数の端末装置２２の複数の受信電力を選択し、その選択した複数の受信電力を用いて、実施の形態１において説明した方法によって、波源位置を推定する。

図３４は、周波数と電力との関係を示す図である。図３４を参照して、複数の端末装置２２の各々は、開始周波数ｆ_{ＳＴＡＲＴ}と終了周波数ｆ_ＥＮＤとを有する帯域幅Ｗｄで周波数スキャンを行い、帯域幅Ｗｄにおける周波数ｆと受信電力Ｐとの関係を記録する。そして、複数の端末装置２２の各々は、帯域幅Ｗｄにおける周波数ｆと受信電力Ｐとの関係に基づいて、１次利用者と２次利用者との間で共用する共用周波数を決定する。そうすると、複数の端末装置２２の各々は、その決定した共用周波数を周波数ｆ_Ｒｉ＿ｊとし、周波数ｆ_Ｒｉ＿ｊにおいて電波を受信して、上述した送信元の識別情報ＩＤ_ＳＲＣ＿ｊ、受信時刻ｔ_Ｒｉ＿ｊ、受信信号Ｓ_Ｒｉ＿ｊおよび受信信号Ｓ_Ｒｉ＿ｊの受信電力Ｐ_ｉ＿ｊを検出し、受信情報ＩＮＦＯ_Ｒ＿ｊを生成する。

また、複数の端末装置２２の各々は、共用周波数である周波数ｆ_Ｔｉ＿ｊにおいて信号を送信すると、送信先の識別情報ＩＤ_ＤＥＳ＿ｊ、送信時刻ｔ_Ｔｉ＿ｊおよび送信信号Ｓ_Ｔｉ＿ｊを検出して送信情報ＩＮＦＯ_Ｔ＿ｊを生成する。

そして、複数の端末装置２２の各々は、受信情報ＩＮＦＯ_Ｒ＿ｊと送信情報ＩＮＦＯ_Ｔ＿ｊとを含むサンプルデータＤ_ＳＰＬ＿ｉを生成し、その生成したサンプルデータＤ_ＳＰＬ＿ｉを波源位置推定装置１Ｂへ送信する。

基地局３１，３２の各々は、端末装置２２と無線通信を行う。基地局３１，３２の各々は、端末装置２２を制御および／または管理して無線通信を行うので、他の通信装置を制御および／または管理するマスター局（通信装置）に相当する。

図３５は、図３３に示す波源位置推定装置１Ｂの概略図である。図３５を参照して、波源位置推定装置１Ｂは、図２に示す波源位置推定装置１の受信手段１２、記憶手段１３および推定手段１４をそれぞれ受信手段１２Ａ、記憶手段１３Ａおよび推定手段１４Ｂに変え、通信推定手段１５を追加したものであり、その他は、波源位置推定装置１と同じである。

受信手段１２Ａは、端末装置２２からサンプルデータＤ_ＳＰＬ＿ｉを受信し、その受信したサンプルデータＤ_ＳＰＬ＿ｉを記憶手段１３Ａに格納する。

記憶手段１３Ａは、受信手段１２Ａによって格納されたサンプルデータＤ_ＳＰＬ＿ｉを記憶する。また、記憶手段１３Ａは、通信システムと、周波数帯と、通信方式と、変調方式と、システム帯域幅と、チャネル帯域幅と、許容ＳＩＮＲ（Signal to Interference Plus Noise Ratio）とを相互に対応付けたデータベースＤＢ１を予め記憶している。

通信推定手段１５は、記憶手段１３ＡからＩ（Ｉは、２以上の整数）個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿Ｉを読み出し、その読み出したＩ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿Ｉに基づいて、後述する方法によって、１次利用者および２次利用者の通信方式または変調方式を推定する。そして、通信推定手段１５は、１次利用者および２次利用者の通信方式または変調方式の推定結果を推定手段１４Ｂへ出力する。

推定手段１４Ｂは、１次利用者および２次利用者の通信方式または変調方式の推定結果を通信推定手段１５から受ける。そして、推定手段１４Ｂは、記憶手段１３ＡからＩ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿Ｉを読み出し、その読み出したＩ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿Ｉと、１次利用者および２次利用者の通信方式または変調方式の推定結果とに基づいて、複数の端末装置２２の複数の受信電力Ｐ_１＿１〜Ｐ_Ｉ＿Ｊのうち、誤差を抑制した複数の受信電力Ｐ＿ＥＲＬＷ＿ｎ（ｎは２以上の整数）を選択し、その選択した複数の受信電力Ｐ＿ＥＲＬＷ＿ｎを用いて、実施の形態１における方法によって波源位置を推定する。

［通信方式または変調方式の推定］
（基本方式）
基本方式を用いた通信方式または変調方式の推定方法について説明する。推定の対象となる通信方式は、時分割複信（ＴＤＤ：Time Division Duplex）、周波数分割複信（ＦＤＤ：Frequency Division Duplex）、セルラー方式およびテレビジョン放送である。テレビジョン放送は、ＦＰＵ（Field Pickup Unit、無線中継伝送装置）によって伝送される。従って、この明細書においては、テレビジョン放送を「ＦＰＵ」と表記することもある。また、推定の対象となる変調方式は、ＡＭ（Amplitude Modulation）変調、ＡＳＫ（Amplitude shift Keying）変調、ＦＭ(Frequency Modulation)変調、ＦＳＫ（Frequency shift Keying）変調、ＰＳＫ（Phase shift Keying）変調、ＣＤＭＡ（Code Division Multiple Access）変調およびＯＦＤＭ（Orthogonal Frequency Division Multiplexing）変調である。

図３６は、通信方式および変調方式の推定方法を説明するための図である。図３６を参照して、端末装置の識別情報（ＩＤ）と、Ｊ個の受信情報ＩＮＦＯ_Ｒ＿１〜ＩＮＦＯ_Ｒ＿ＪおよびＪ個の送信情報ＩＮＦＯ_Ｔ＿１〜ＩＮＦＯ_Ｔ＿Ｊとは、相互に対応付けられている。

通信推定手段１５は、Ｉ個のサンプルデータＤ_ＳＰＬ＿１〜サンプルデータＤ_ＳＰＬ＿Ｉを記憶手段１３Ａから読み出し、その読み出したＩ個のサンプルデータＤ_ＳＰＬ＿１〜サンプルデータＤ_ＳＰＬ＿Ｉに基づいて、図３６に示すように、Ｊ個の受信情報ＩＮＦＯ_Ｒ＿１〜ＩＮＦＯ_Ｒ＿ＪおよびＪ個の送信情報ＩＮＦＯ_Ｔ＿１〜ＩＮＦＯ_Ｔ＿ＪをそれぞれＩ個の端末装置の識別情報（ＩＤ）に対応付ける。
（ｉ）ＴＤＤの推定
そして、通信推定手段１５は、識別情報ＩＤ_１を有する端末装置について、次の項目を実行する。
（１−１）Ｊ個の送信元の識別情報ＩＤ_ＳＲＣ１＿１〜ＩＤ_ＳＲＣ１＿Ｊのうち、相互に異なり、かつ、１次利用者を表す送信元の識別情報が２個以上有るか否かを判定する。
（１−２）Ｊ個の送信先の識別情報ＩＤ_ＤＥＳ１＿１〜ＩＤ_ＤＥＳ１＿Ｊのうち、相互に異なり、かつ、１次利用者を表す送信先の識別情報が２個以上有るか否かを判定する。
（１−３）項目（１−１）の判定結果が肯定的であるときの２個以上の送信元の識別情報が、項目（１−２）の判定結果が肯定的であるときの２個以上の送信先の識別情報とそれぞれ一致するか否かを判定する。
（１−４）項目（１−３）の判定結果が肯定的であるときの相互に一致する送信元および送信先が使用する２つの周波数ｆ_Ｒｊ＿１，ｆ_Ｔｊ＿１が一致するか否かを判定する。

通信推定手段１５は、項目（１−１）〜項目（１−４）の全ての判定が肯定的であるとき、端末装置ＩＤ_１および２以上の１次利用者の通信方式がＴＤＤであると推定する。

ＴＤＤにおいては、同一チャネルで２以上の利用者が居ることが前提であり、波源も２個以上観測されることになる。一方、項目（１−１）〜項目（１−４）の全ての判定が肯定的であるとき、同一周波数で端末装置ＩＤ_１と無線通信を行っている１次利用者が２以上存在すると判定されることになるので、項目（１−１）〜項目（１−４）の全ての判定が肯定的であるとき、通信方式がＴＤＤであると推定することにしたものである。

そして、通信推定手段１５は、項目（１−１）〜項目（１−４）の少なくとも１つの判定が否定的であるとき、端末装置ＩＤ_１、および端末装置ＩＤ_１と無線通信を行っている相手方の通信方式は、ＴＤＤでないと推定する。
（ｉｉ）ＯＦＤＭ，ＡＭ変調，ＡＳＫ変調の推定
通信推定手段１５は、Ｊ個の受信電力ＲＳＳＩ_１＿１〜ＲＳＳＩ_１＿Ｊに基づいて、受信電力の変動幅ΔＲＳＳＩ_１を求める。そして、通信推定手段１５は、変動幅ΔＲＳＳＩ_１が、端末装置ＩＤ_ｉが予め保持している自信号の電力の２倍以上であるとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＯＦＤＭ変調、ＡＭ変調およびＡＳＫ変調のいずれかであると推定する。

一方、通信推定手段１５は、変動幅ΔＲＳＳＩ_１が、端末装置ＩＤ_ｉが予め保持している自信号の電力の２倍以上でないとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＯＦＤＭ変調、ＡＭ変調およびＡＳＫ変調のいずれでもないと推定する。

通信推定手段１５は、変動幅ΔＲＳＳＩ_１が、端末装置ＩＤ_ｉが予め保持している自信号の電力の２倍以上であると判定したとき、受信信号Ｓ_Ｒ１＿１〜Ｓ_Ｒ１＿Ｊの時間応答Ｓ_Ｒ１＿１（ｔ）〜Ｓ_Ｒ１＿Ｊ（ｔ）を生成し、その生成した時間応答Ｓ_Ｒ１＿１（ｔ）〜Ｓ_Ｒ１＿Ｊ（ｔ）に基づいて自己相関を演算する。そして、通信推定手段１５は、自己相関が周期的に表れる２つのピークからなるとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＯＦＤＭ変調であると推定し、自己相関が周期的に表れる２つのピークからならないとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＡＭ変調またはＡＳＫ変調であると推定する。

図３７は、自己相関の概念図である。図３７を参照して、変調方式がＯＦＤＭ変調であるとき、自己相関は、２つのピークが周期Ｔ１で現れる。これは、ＯＦＤＭの場合、サイクリックエクステンション（ガード期間）が含まれるので、相関値に周期性が現れるからである。従って、図３７に示すように、自己相関が周期的に表れる２つのピークからなるとき、変調方式をＯＦＤＭ変調であると推定することにしたものである。
（ｉｉｉ）ＣＤＭＡ，ＦＭ変調，ＦＳＫ変調の推定
通信推定手段１５は、サンプルデータの先頭の受信電力ＲＳＳＩ_１＿１の周波数ｆ_Ｒ１＿１が周波数偏移するとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＣＤＭＡ変調、ＦＭ変調およびＦＳＫ変調のいずれかであると推定する。

一方、通信推定手段１５は、先頭の受信電力ＲＳＳＩ_１＿１の周波数ｆ_Ｒ１＿１が周波数偏移しないとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＣＤＭＡ変調、ＦＭ変調およびＦＳＫ変調のいずれでもないと推定する。

通信推定手段１５は、先頭の受信電力ＲＳＳＩ_１＿１の周波数ｆ_Ｒ１＿１が周波数偏移すると推定したとき、上述した方法によって受信信号Ｓ_Ｒ１＿１〜Ｓ_Ｒ１＿Ｊの自己相関を演算する。そして、通信推定手段１５は、自己相関が周期的に表れる１つのピークからなるとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＣＤＭＡ変調であると推定し、自己相関が周期的に表れる１つのピークからならないとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＦＭ変調またはＦＳＫ変調であると判定する。

図３８は、自己相関の別の概念図である。図３８を参照して、変調方式がＣＤＭＡ変調であるとき、自己相関は、１つのピークが周期Ｔ２で現れる。これは、ＣＤＭＡ変調の場合、同一の拡散コードが使用されている間、そのコードによる相関が現れるからである。従って、図３８に示すように、自己相関が周期的に表れる１つのピークからなるとき、変調方式をＣＤＭＡ変調であると推定することにしたものである。
（ｉｖ）ＣＤＭＡ，ＰＳＫ変調の推定
通信推定手段１５は、変動幅ΔＲＳＳＩ_１が閾値ΔＲＳＳＩ＿ｔｈ以下であるとき、即ち、受信電力ＲＳＳＩ_１＿１〜ＲＳＳＩ_１＿Ｊの変動が殆ど無いとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＣＤＭＡ変調またはＰＳＫ変調であると推定する。

一方、通信推定手段１５は、変動幅ΔＲＳＳＩ_１が閾値ΔＲＳＳＩ＿ｔｈよりも大きいとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＣＤＭＡ変調およびＰＳＫ変調のいずれでもないと推定する。

通信推定手段１５は、変動幅ΔＲＳＳＩ_１が閾値ΔＲＳＳＩ＿ｔｈ以下であると判定したとき、自己相関が、周期的に表れる１つのピークからなるか否かを更に判定する。そして、通信推定手段１５は、自己相関が、周期的に表れる１つのピークからなるとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＣＤＭＡ変調であると推定し、自己相関が、周期的に表れる１つのピークからならないとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の変調方式をＰＳＫ変調であると推定する。
（ｖ）セルラー方式またはテレビジョン放送の推定
通信推定手段１５は、上述した（ｉ）〜（ｉｖ）における推定結果がＴＤＤ、ＡＭ変調、ＡＳＫ変調、ＦＭ変調、ＦＳＫ変調、ＰＳＫ変調、ＣＤＭＡ変調およびＯＦＤＭ変調のいずれでもないとき、端末装置ＩＤ_１および端末装置ＩＤ_１と無線通信を行っている相手方の通信方式をセルラー方式またはテレビジョン放送であると推定する。

なお、通信推定手段１５は、上述した（ｉ）〜（ｖ）における方法によって、端末装置ＩＤ_２〜ＩＤ_Ｉ、および端末装置ＩＤ_２〜ＩＤ_Ｉと無線通信を行っている相手方の通信方式または変調方式を推定する。

通信推定手段１５は、通信方式の推定結果がＴＤＤであるとき、通信方式または変調方式の推定結果がＴＤＤであることを示す信号Ｓ_ＴＤＤを生成して推定手段１４Ｂへ出力し、通信方式または変調方式の推定結果がＡＭ変調またはＡＳＫ変調であるとき、通信方式または変調方式の推定結果がＡＭ変調またはＡＳＫ変調であることを示す信号Ｓ_Ａを生成して推定手段１４Ｂへ出力し、通信方式または変調方式の推定結果がＦＭ変調またはＦＳＫ変調であるとき、通信方式または変調方式の推定結果がＦＭ変調またはＦＳＫ変調であることを示す信号Ｓ_Ｆを生成して推定手段１４Ｂへ出力する。

また、通信推定手段１５は、通信方式または変調方式の推定結果がＯＦＤＭ変調であるとき、通信方式または変調方式の推定結果がＯＦＤＭ変調であることを示す信号Ｓ_ＯＦを生成して推定手段１４Ｂへ出力し、通信方式または変調方式の推定結果がＣＤＭＡ変調であるとき、通信方式または変調方式の推定結果がＣＤＭＡ変調であることを示す信号Ｓ_ＣＤを生成して推定手段１４Ｂへ出力し、通信方式または変調方式の推定結果がＰＳＫ変調であるとき、通信方式または変調方式の推定結果がＰＳＫ変調であることを示す信号Ｓ_Ｐを生成して推定手段１４Ｂへ出力し、通信方式または変調方式の推定結果がセルラシステムまたはテレビジョン方法であるとき、通信方式または変調方式の推定結果がセルラシステムまたはテレビジョン方法であることを示す信号Ｓ_ＴＣを生成して推定手段１４Ｂへ出力する。

（テーブル方式）
テーブルを用いた通信方式または変調方式の推定方法について説明する。通信推定手段１５は、テーブル方式においては、データベースＤＢ１を参照して通信方式または変調方式を推定する。

図３９は、データベースＤＢ１の概略図である。図３９を参照して、データベースＤＢ１は、通信システムと、周波数帯と、通信方式と、変調方式と、システム帯域幅と、チャネル帯域幅と、許容ＳＩＮＲとを含む。通信システム、周波数帯、通信方式、変調方式、システム帯域幅、チャネル帯域幅および許容ＳＩＮＲは、相互に対応付けられる。

通信システムは、例えば、無線システムＡ、無線システムＢ、無線システムＣ、無線システムＤ、無線システムＥ、無線システムＦおよび無線システムＧからなる。無線システムＡ〜無線システムＧは、例えば、ＦＰＵ、アマチュア無線および防災無線等である。

７００ＭＨｚの周波数帯、ＦＤＤの通信方式（複信方式）、ＣＤＭＡの変調方式、１０ＭＨｚのシステム帯域幅、２０ｋＨｚのチャネル帯域幅および１０ｄＢの許容ＳＩＮＲは、無線システムＡに対応付けられる。

１．２ＧＨｚの周波数帯、ＴＤＤの通信方式（複信方式）、ＯＦＤＭの変調方式、２０ＭＨｚのシステム帯域幅、５０ｋＨｚのチャネル帯域幅および１２ｄＢの許容ＳＩＮＲは、無線システムＢに対応付けられる。

２００ＭＨｚの周波数帯、ＴＤＤの通信方式（複信方式）、ＡＳＫの変調方式、５ＭＨｚのシステム帯域幅、１００ｋＨｚのチャネル帯域幅および９ｄＢの許容ＳＩＮＲは、無線システムＣに対応付けられる。

１．７ＧＨｚの周波数帯、ＦＤＤの通信方式（複信方式）、ＡＭの変調方式、１ＭＨｚのシステム帯域幅、１０ｋＨｚのチャネル帯域幅および０ｄＢの許容ＳＩＮＲは、無線システムＤに対応付けられる。

２．３ＧＨｚの周波数帯、ＦＤＤの通信方式（複信方式）、ＦＭの変調方式、５００ｋＨｚのシステム帯域幅、１５ｋＨｚのチャネル帯域幅および１５ｄＢの許容ＳＩＮＲは、無線システムＥに対応付けられる。

３．３ＧＨｚの周波数帯、ＴＤＤの通信方式（複信方式）、ＦＳＫの変調方式、２ＭＨｚのシステム帯域幅、５ｋＨｚのチャネル帯域幅および４０ｄＢの許容ＳＩＮＲは、無線システムＦに対応付けられる。

８００ＭＨｚの周波数帯、ＴＤＤの通信方式（複信方式）、ＰＳＫの変調方式、１０ＭＨｚのシステム帯域幅、１０ｋＨｚのチャネル帯域幅および２０ｄＢの許容ＳＩＮＲは、無線システムＧに対応付けられる。

通信推定手段１５は、サンプルデータＤ_ＳＰＬ＿ｉに基づいて周波数帯およびシステム帯域幅を検出する。各サンプルデータＤ_ＳＰＬ＿ｉは、Ｊ個の受信情報ＩＮＦＯ_Ｒ＿１〜ＩＮＦＯ_Ｒ＿Ｊからなるので、Ｊ個の受信情報ＩＮＦＯ_Ｒ＿１〜ＩＮＦＯ_Ｒ＿ＪのＪ個の周波数ｆ_Ｒｉ＿１〜ｆ_Ｒｉ＿Ｊに基づいて周波数帯およびシステム帯域幅を検出することができる。

通信推定手段１５は、周波数帯およびシステム帯域幅を検出すると、その検出した周波数帯およびシステム帯域幅に対応する通信方式または変調方式を検出することによって通信方式または変調方式を推定する。例えば、通信推定手段１５は、サンプルデータＤ_ＳＰＬ＿ｉに基づいて、１．２ＧＨｚの周波数帯および２０ＭＨｚのシステム帯域幅を検出した場合、１．２ＧＨｚの周波数帯および２０ＭＨｚのシステム帯域幅に対応するＴＤＤおよびＯＦＤＭを検出することによって、通信方式がＴＤＤであり、変調方式がＯＦＤＭであると推定する。通信推定手段１５は、他の周波数帯およびシステム帯域幅が検出された場合、同様にして通信方式または変調方式を推定する。

なお、通信推定手段１５は、データベースＤＢ１を用いて通信方式または変調方式を推定する場合、通信方式および変調方式の両方を推定してもよく、通信方式および変調方式のいずれか一方だけを推定してもよい。即ち、通信推定手段１８は、通信方式および変調方式の少なくとも１つを推定すればよい。

実施の形態２においては、上述した基本方式による通信方式または変調方式の推定方法とテーブル方式による通信方式または変調方式の推定方法とを組み合わせてもよい。この場合、テーブル方式による通信方式または変調方式の推定方法を実行した後に基本方式による通信方式または変調方式の推定方法を実行してもよく、基本方式による通信方式または変調方式の推定方法を実行した後にテーブル方式による通信方式または変調方式の推定方法を実行してもよい。

［波源位置の推定］
（ＴＤＤ）
推定手段１４Ｂは、通信方式または変調方式の推定結果がＴＤＤであるとき、信号Ｓ_ＴＤＤを通信推定手段１５から受け、その受けた信号Ｓ_ＴＤＤに応じて、送信元が基地局であることを示す識別情報ＩＤ_ＳＲＣ＿ｊに対応付けられた受信電力Ｐ_ｉ＿ｊおよび位置情報［ｘ_ｉ，ｙ_ｉ］をサンプルデータＤ_ＳＰＬ＿ｉから検出する。この場合、ＴＤＤにおいては、複数の波源（＝複数の基地局）が存在するため、推定手段１４Ｂは、複数の受信電力Ｐ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］をサンプルデータＤ_ＳＰＬ＿ｉから検出する。

そして、推定手段１４Ｂは、その検出した複数の受信電力Ｐ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］を式（１）に代入して重心［ｘ_Ｇ，ｙ_Ｇ］を算出し、その算出した重心［ｘ_Ｇ，ｙ_Ｇ］を波源位置として推定する。

その後、推定手段１４Ｂは、その推定した波源位置に誤差がある場合、実施の形態１において説明した方法によって複数の仮想端末装置２ｖｊを追加し、複数の受信電力Ｐ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］と、複数の仮想端末装置２ｖｊの複数の受信電力Ｐ_ｊｖおよび複数の位置情報［ｘ_ｊｖ，ｙ_ｊｖ］とを式（１）に代入して重心［ｘ_Ｇ，ｙ_Ｇ］を算出し、その算出した重心［ｘ_Ｇ，ｙ_Ｇ］を波源位置として推定する。

このように、推定手段１４Ｂは、通信方式または変調方式の推定結果がＴＤＤであるとき、マスター局（基地局）から送信された電波の受信電力Ｐ_ｉ＿ｊに基づいて波源位置を推定する。

ＴＤＤの場合、複数の波源が存在する。また、スペクトルアナライザーによる電力測定においては、電波の送信と受信のフェーズが電力のデューティ比として現れるため、電力測定に誤差が発生する。

しかし、通信方式または変調方式の推定結果がＴＤＤであるとき、マスター局（基地局）からの受信電力Ｐ_ｉ＿ｊを用いて波源位置を推定するので、誤差を抑制した受信電力Ｐ＿ＥＲＬＷを用いて波源位置を推定できる。従って、波源位置を正確に推定できる。
（ＣＤＭＡ）
推定手段１４Ｂは、通信方式または変調方式の推定結果がＣＤＭＡ変調であるとき、信号Ｓ_ＣＤを通信推定手段１５から受け、その受けた信号Ｓ_ＣＤに応じて、送信元の識別情報ＩＤ_ＳＲＣ＿ｊを参照してフォワードリンク（＝ダウンリンク）の送信元の識別情報ＩＤ_ＳＲＣ＿ｊに対応付けられた複数の受信電力Ｐ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］をサンプルデータＤ_ＳＰＬ＿ｉから検出する。

そして、推定手段１４Ｂは、その検出した複数の受信電力ＲＳＳＩ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］を式（１）に代入して重心［ｘ_Ｇ，ｙ_Ｇ］を算出し、その算出した重心［ｘ_Ｇ，ｙ_Ｇ］を波源位置（フォワードリンクの送信元の位置）として推定する。なお、推定手段１４Ｂは、波源位置の推定においては、リバースリンク（＝アップリンク）の受信電力Ｐ_ｉ＿ｊを用いない。

その後、推定手段１４Ｂは、その推定した波源位置に誤差がある場合、実施の形態１において説明した方法によって複数の仮想端末装置２ｖｊを追加し、複数の受信電力Ｐ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］と、複数の仮想端末装置２ｖｊの複数の受信電力Ｐ_ｊｖおよび複数の位置情報［ｘ_ｊｖ，ｙ_ｊｖ］とを式（１）に代入して重心［ｘ_Ｇ，ｙ_Ｇ］を算出し、その算出した重心［ｘ_Ｇ，ｙ_Ｇ］を波源位置として推定する。

このように、推定手段１４Ｂは、通信方式または変調方式の推定結果がＣＤＭＡ変調であるとき、リバースリンク（＝アップリンク）の受信電力ＲＳＳＩ_ｉ＿ｊを用いずに、フォワードリンクの送信元からの受信電力Ｐ_ｉ＿ｊに基づいて波源を推定する。

ＣＤＭＡ変調の場合、フォワードリンクの送信電力が定電力であるので、フォワードリンクにおける受信電力Ｐ_ｉ＿ｊを用いることによって波源を正確に推定できる。
（ＦＭ変調）
推定手段１４Ｂは、通信方式または変調方式の推定結果がＦＭ変調であるとき、信号Ｓ_Ｆを通信推定手段１５から受け、その受けた信号Ｓ_Ｆに応じて、周波数偏移している複数の受信信号Ｓ_Ｒｉ＿ｊに対応する複数の受信電力Ｐ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］をサンプルデータＤ_ＳＰＬ＿ｉから検出する。なお、周波数偏移している受信信号Ｓ_Ｒｉ＿ｊに対応する受信電力Ｐ_ｉ＿ｊは、変調信号による周波数偏移の範囲を含む帯域幅を有するフィルタを用いて検出された受信電力である。

そして、推定手段１４Ｂは、その検出した複数の受信電力ＲＳＳＩ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］を式（１）に代入して重心［ｘ_Ｇ，ｙ_Ｇ］を算出し、その算出した重心［ｘ_Ｇ，ｙ_Ｇ］を波源位置として推定する。

その後、推定手段１４Ｂは、その推定した波源位置に誤差がある場合、実施の形態１において説明した方法によって複数の仮想端末装置２ｖｊを追加し、複数の受信電力Ｐ_ｉ＿ｊおよび複数の位置情報［ｘ_ｉ，ｙ_ｉ］と、複数の仮想端末装置２ｖｊの複数の受信電力Ｐ_ｊｖおよび複数の位置情報［ｘ_ｊｖ，ｙ_ｊｖ］とを式（１）に代入して重心［ｘ_Ｇ，ｙ_Ｇ］を算出し、その算出した重心［ｘ_Ｇ，ｙ_Ｇ］を波源位置として推定する。

このように、推定手段１４Ｂは、通信方式または変調方式の推定結果がＦＭ変調であるとき、変調信号による周波数偏移の範囲を含む帯域幅を有するフィルタを用いて検出された受信電力Ｐ_ｉ＿ｊを用いて波源位置を推定するので、波源位置を正確に推定できる。

受信電力を検出するときのフィルタが周波数偏移に対して狭い帯域幅を有するとき、受信電力は、変調信号に応じて変化する。

しかし、変調信号による周波数偏移の範囲を含む帯域幅を有するフィルタを用いて受信電力Ｐ_ｉ＿ｊを検出した場合、受信信力Ｐ_ｉ＿ｊは、変調信号に応じて変化しない。従って、波源位置を正確に推定できる。

なお、通信方式または変調方式の推定結果がＴＤＤ、ＣＤＭＡおよびＦＭ変調以外であるとき、推定手段１４Ｂは、実施の形態１における推定方法を用いて波源位置を推定する。

図４０は、図３３に示す波源位置推定装置１Ｂの動作を説明するためのフローチャートである。図４０に示すフローチャートは、図１１に示すフローチャートのステップＳ１をステップＳ１２１〜ステップＳ１２５に変えたものであり、その他は、図１１に示すフローチャートと同じである。

図４０を参照して、波源位置推定装置１Ｂの動作が開始されると、受信手段１２Ａは、複数の端末装置２２から一定時間におけるサンプルデータＤ_ＳＰＬ＿ｉを受信し、その受信したサンプルデータＤ_ＳＰＬ＿ｉを記憶手段１３Ａに格納する。即ち、波源位置推定装置１Ｂは、一定時間におけるサンプルデータＤ_ＳＰＬ＿ｉを取得する（ステップＳ１２１）。

そして、通信推定手段１５は、記憶手段１３ＡからサンプルデータＤ_ＳＰＬ＿ｉを読み出し、その読み出したサンプルデータＤ_ＳＰＬ＿ｉに基づいて、共用周波数を用いる１次利用者および２次利用者の通信方式または変調方式を推定し（ステップＳ１２２）、通信方式または変調方式の推定結果ＥＳＴを推定手段１４Ｂへ出力する。

その後、推定手段１４Ｂは、通信方式または変調方式の推定結果ＥＳＴを通信推定手段１５から受け、その受けた通信方式または変調方式の推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかであるか否かを判定する（ステップＳ１２３）。

ステップＳ１２３において、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかであると判定されたとき、推定手段１４Ｂは、推定結果ＥＳＴに基づいて、サンプルデータＤ_ＳＰＬ＿ｉから複数の実端末装置の複数の受信電力のうちの誤差を抑制した複数の受信電力を検出するとともに、その検出した複数の受信電力に対応する複数の位置情報を検出する（ステップＳ１２４）。

一方、ステップＳ１２３において、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれでもないと判定されたとき、推定手段１４Ｂは、受信電力を峻別することなく、サンプルデータＤ_ＳＰＬ＿ｉから複数の実端末装置の複数の受信電力および複数の位置情報を検出する（ステップＳ１２５）。

そして、ステップＳ１２４またはステップＳ１２５の後、上述したステップＳ２〜ステップＳ６が順次実行され、ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａであると判定されたとき、波源位置推定装置１Ｂにおける波源位置を推定する動作が終了する。

図４０に示すフローチャートにおいては、ステップＳ１２４からステップＳ２へ移行した場合、ステップＳ２において、推定手段１４Ｂは、誤差を抑制した複数の受信電力を用いて波源位置［ｘ_Ｇ，ｙ_Ｇ］を推定する。従って、正確に波源位置［ｘ_Ｇ，ｙ_Ｇ］を推定できる。

誤差を抑制した複数の受信電力を用いて波源位置を推定した場合でも、実端末装置の分布に偏りがあれば、その推定した波源位置［ｘ_Ｇ，ｙ_Ｇ］に誤差があることになる。

従って、ステップＳ２，Ｓ３が実行された後、ステップＳ４において、その推定した波源位置［ｘ_Ｇ，ｙ_Ｇ］に誤差があると判定されたとき（ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａでないと判定されたとき）、上述したステップＳ５，Ｓ６が順次実行された後、ステップＳ２へ移行する。

そして、ステップＳ２において、推定手段１４Ｂは、複数の実端末装置の複数の位置情報および誤差を抑制した複数の受信電力と、仮想端末装置の複数の位置情報および複数の受信電力とを用いて波源位置を推定する。従って、波源位置を精度良く推定できる。

なお、ステップＳ１２５からステップＳ２へ移行した場合、推定手段１４Ｂは、受信電力を峻別することなくサンプルデータＤ_ＳＰＬ＿ｉから検出した複数の受信電力および複数の位置情報を用いて波源位置を推定する（ステップＳ２参照）。

図４１は、図４０に示すステップＳ１２２の詳細な動作を説明するためのフローチャートである。

図４１を参照して、図４０のステップＳ１２１の後、通信推定手段１５は、サンプルデータＤ_ＳＰＬ＿ｉに基づいて受信電力の電力変動および先頭電力の周波数偏移を検出する（ステップＳ１２２１）。

そして、通信推定手段１５は、電力変動の変動幅を検出し（ステップＳ１２２２）、受信信号の自己相関を演算する（ステップＳ１２２３）。

その後、通信推定手段１５は、波源が複数であるか否かを判定する（ステップＳ１２２４）。

ステップＳ１２２４において、波源が複数であると判定されたとき、通信推定手段１５は、通信方式がＴＤＤであると推定する（ステップＳ１２２５）。

一方、ステップＳ１２２４において、波源が複数でないと判定されたとき、またはステップＳ１２２５の後、通信推定手段１５は、電力変動の変動幅が、端末装置ＩＤ_ｉにおける自信号の電力の２倍以上であるか否かを判定する（ステップＳ１２２６）。

ステップＳ１２２６において、電力変動の変動幅が、端末装置ＩＤ_ｉにおける自信号の電力の２倍以上であると判定されたとき、通信推定手段１５は、変調方式がＯＦＤＭ変調またはＡＭ変調またはＡＳＫ変調であると推定する（ステップＳ１２２７）。

一方、ステップＳ１２２６において、電力変動の変動幅が、端末装置ＩＤ_ｉにおける自信号の電力の２倍以上でないと判定されたとき、またはステップＳ１２２７の後、通信推定手段１５は、先頭電力が周波数偏移するか否かを判定する（ステップＳ１２２８）。

ステップＳ１２２８において、先頭電力が周波数偏移すると判定されたとき、通信推定手段１５は、変調方式がＣＤＭＡ変調またはＦＭ変調またはＦＳＫ変調であると推定する（ステップＳ１２２９）。

一方、ステップＳ１２２８において、先頭電力が周波数偏移しないと判定されたとき、またはステップＳ１２２９の後、通信推定手段１５は、電力変動の変動幅が閾値以下であるか否かを判定する（ステップＳ１２３０）。

ステップＳ１２３０において、電力変動の変動幅が閾値以下であると判定されたとき、通信推定手段１５は、変調方式がＣＤＭＡ変調またはＰＳＫ変調であると推定する（ステップＳ１２３１）。

一方、ステップＳ１２３０において、電力変動の変動幅が閾値以下でないと判定されたとき、通信推定手段１５は、通信方式がセルラー方式またはテレビジョン放送であると推定する（ステップＳ１２３２）。

そして、ステップＳ１２３１またはステップＳ１２３２の後、一連の動作は、図４０のステップＳ１２３へ移行する。

図４２は、図４１のステップＳ１２２７の詳細な動作を説明するためのフローチャートである。図４２を参照して、図４１のステップＳ１２２６において、電力変動の変動幅が、端末装置ＩＤ_ｉにおける自信号の電力の２倍以上であると判定されたとき、通信推定手段１５は、自己相関が周期的に現れる２つのピークからなるか否かを判定する（ステップＳ１２２７Ａ）。

ステップＳ１２２７Ａにおいて、自己相関が周期的に現れる２つのピークからなると判定されたとき、通信推定手段１５は、変調方式がＯＦＤＭ変調であると推定する（ステップＳ１２２７Ｂ）。

一方、ステップＳ１２２７Ａにおいて、自己相関が周期的に現れる２つのピークからならないと判定されたとき、通信推定手段１５は、変調方式がＡＭ変調またはＡＳＫ変調であると推定する（ステップＳ１２２７Ｃ）。そして、ステップＳ１２２７ＢまたはステップＳ１２２７Ｃの後、一連の動作は、図４１のステップＳ１２２８へ移行する。

図４３は、図４１のステップＳ１２２９の詳細な動作を説明するためのフローチャートである。図４３を参照して、図４１のステップＳ１２２８において、先頭電力が周波数偏移すると判定されたとき、通信推定手段１５は、自己相関が周期的に現れる１つのピークからなるか否かを判定する（ステップＳ１２２９Ａ）。

ステップＳ１２２９Ａにおいて、自己相関が周期的に現れる１つのピークからなると判定されたとき、通信推定手段１５は、変調方式がＣＤＭＡ変調であると推定する（ステップＳ１２２９Ｂ）。

一方、ステップＳ１２２９Ａにおいて、自己相関が周期的に現れる１つのピークからならないと判定されたとき、通信推定手段１５は、変調方式がＦＭ変調またはＦＳＫ変調であると推定する（ステップＳ１２２９Ｃ）。そして、ステップＳ１２２９ＢまたはステップＳ１２２９Ｃの後、一連の動作は、図４１のステップＳ１２３０へ移行する。

図４４は、図４１のステップＳ１２３１の詳細な動作を説明するためのフローチャートである。図４４を参照して、図４１のステップＳ１２３０において、電力変動の変動幅が閾値以下であると判定されたとき、通信推定手段１５は、自己相関が周期的に現れる１つのピークからなるか否かを判定する（ステップＳ１２３１Ａ）。

ステップＳ１２３１Ａにおいて、自己相関が周期的に現れる１つのピークからなると判定されたとき、通信推定手段１５は、変調方式がＣＤＭＡ変調であると推定する（ステップＳ１２３１Ｂ）。

一方、ステップＳ１２３１Ａにおいて、自己相関が周期的に現れる１つのピークからならないと判定されたとき、通信推定手段１５は、変調方式がＰＳＫ変調であると推定する（ステップＳ１２３１Ｃ）。そして、ステップＳ１２３１ＢまたはステップＳ１２３１Ｃの後、一連の動作は、図４０のステップＳ１２３へ移行する。

図４５は、図４０のステップＳ１２２の別の詳細な動作を説明するためのフローチャートである。なお、図４５においては、通信方式または変調方式の推定結果を「ＥＳＴ１」または「ＥＳＴ２」と表す。

図４５を参照して、図４０のステップＳ１２１の後、通信推定手段１５は、上述した図４１のステップＳ１２３１〜ステップＳ１２３２を順次実行する（ステップＳ１２３３）。そして、通信推定手段１５は、推定した通信方式または変調方式を推定結果ＥＳＴ１として保持する（ステップＳ１２３４）。

その後、通信推定手段１５は、サンプルデータＤ_ＳＰＬ＿ｉを記憶手段１３Ａから読み出し、その読み出したサンプルデータＤ_ＳＰＬ＿ｉに基づいて周波数帯およびシステム帯域幅を検出する（ステップＳ１２３５）。

引き続いて、通信推定手段１５は、データベースＤＢ１を記憶手段１３Ａから読み出し、その読み出したデータベースＤＢ１を参照して、検出した周波数帯およびシステム帯域幅に対応する通信方式または変調方式を検出し、その検出した通信方式または変調方式を推定結果ＥＳＴ２とする（ステップＳ１２３６）。

そうすると、通信推定手段１５は、推定結果ＥＳＴ１が推定結果ＥＳＴ２に一致するか否かを判定する（ステップＳ１２３７）。

ステップＳ１２３７において、推定結果ＥＳＴ１が推定結果ＥＳＴ２に一致しないと判定されたとき、一連の動作は、ステップＳ１２３３へ移行する。その後、ステップＳ１２３７において、推定結果ＥＳＴ１が推定結果ＥＳＴ２に一致すると判定されるまで、ステップＳ１２３３〜ステップＳ１２３７が繰り返し実行される。

そして、ステップＳ１２３７において、推定結果ＥＳＴ１が推定結果ＥＳＴ２に一致すると判定されたとき、通信推定手段１５は、一致した推定結果を最終的な推定結果とする（ステップＳ１２３８）。その後、一連の動作は、図４０のステップＳ１２３へ移行する。

なお、図４５に示すフローチャートにおいては、ステップＳ１２３３に代えて、ステップＳ１２３５〜ステップＳ１２３６を実行し、ステップＳ１２３５〜ステップＳ１２３６に代えて、ステップＳ１２３３を実行してもよい。この場合、ステップＳ１２３６における推定結果が推定結果ＥＳＴ１となり、ステップＳ１２３３における推定結果が推定結果ＥＳＴ２になる。

図４６は、図４０のステップＳ１２４の詳細な動作を説明するためのフローチャートである。

図４６を参照して、図４０のステップＳ１２３において、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかであると判定されたとき、推定手段１４Ｂは、通信推定手段１５から信号Ｓ_ＴＤＤを受けたか否かを判定することによって推定結果ＥＳＴがＴＤＤであるか否かを判定する（ステップＳ１２４１）。この場合、推定手段１４Ｂは、通信推定手段１５から信号Ｓ_ＴＤＤを受けているとき、推定結果ＥＳＴがＴＤＤであると判定し、通信推定手段１５から信号Ｓ_ＴＤＤを受けていないとき、推定結果ＥＳＴがＴＤＤではないと判定する。

ステップＳ１２４１において、推定結果ＥＳＴがＴＤＤであると判定されたとき、推定手段１４Ｂは、他の通信装置を制御および／または管理する通信装置から送信された複数の電波の複数の受信電力と、その複数の受信電力に対応する複数の位置情報とをサンプルデータＤ_ＳＰＬ＿ｉから検出する（ステップＳ１２４２）。

一方、ステップＳ１２４１において、推定結果ＥＳＴがＴＤＤでないと判定されたとき、推定手段１４Ｂは、通信推定手段１５から信号Ｓ_ＣＤを受けたか否かを判定することによって推定結果ＥＳＴがＣＤＭＡ変調であるか否かを判定する（ステップＳ１２４３）。この場合、推定手段１４Ｂは、通信推定手段１５から信号Ｓ_ＣＤを受けているとき、推定結果ＥＳＴがＣＤＭＡ変調であると判定し、通信推定手段１５から信号Ｓ_ＣＤを受けていないとき、推定結果ＥＳＴがＣＤＭＡではないと判定する。

ステップＳ１２４３において、推定結果ＥＳＴがＣＤＭＡ変調であると判定されたとき、推定手段１４Ｂは、フォワードリンクにおける複数の電波の複数の受信電力と、その複数の受信電力に対応する複数の位置情報とをサンプルデータＤ_ＳＰＬ＿ｉから検出する（ステップＳ１２４４）。

一方、ステップＳ１２４３において、推定結果ＥＳＴがＣＤＭＡ変調でないと判定されたとき、推定手段１４Ｂは、変調信号による周波数偏移を含む帯域幅を有するフィルタで検出された複数の電波の複数の受信電力と、その複数の受信電力に対応する複数の位置情報とをサンプルデータＤ_ＳＰＬ＿ｉから検出する（ステップＳ１２４５）。

そして、ステップＳ１２４２、ステップＳ１２４４およびステップＳ１２４５のいずれかの後、図４０のステップＳ２へ移行する。

なお、ステップＳ１２４３において、推定結果ＥＳＴがＣＤＭＡ変調でないと判定されたとき、推定結果ＥＳＴがＦＭ変調であるか否かを判定しないのは、図４１のフローチャートは、図４０のステップＳ１２３において推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかであると判定されたことを前提としており、ステップＳ１２４１において推定結果ＥＳＴがＴＤＤでないと判定されており、かつ、ステップＳ１２４３において推定結果ＥＳＴがＣＤＭＡでないと判定されているので、推定結果ＥＳＴがＦＭ変調であることが明らかであるからである。

図４７は、図３３に示す波源位置推定装置１Ｂの別の動作を説明するためのフローチャートである。図４７に示すフローチャートは、図１６に示すフローチャートのステップＳ１をステップＳ１２１〜ステップＳ１２５に変えたものであり、その他は、図１６に示すフローチャートと同じである。

図４７を参照して、波源位置推定装置１Ｂの動作が開始されると、上述したステップＳ１２１〜ステップＳ１２５が順次実行される。そして、ステップＳ１２４またはステップＳ１２５の後、上述したステップＳ３Ａ，Ｓ２，Ｓ４Ａ，Ｓ４Ｂ，Ｓ５，Ｓ６が順次実行され、ステップＳ４Ａにおいて、ｋ＝Ｋであると判定されると、波源位置推定装置１Ｂの動作が終了する。

図４８は、図３３に示す波源位置推定装置１Ｂの更に別の動作を説明するためのフローチャートである。図４８に示すフローチャートは、図１７に示すフローチャートのステップＳ１をステップＳ１２１〜ステップＳ１２５に変えたものであり、その他は、図１７に示すフローチャートと同じである。

図４８を参照して、波源位置推定装置１Ｂの動作が開始されると、上述したステップＳ１２１〜ステップＳ１２５が順次実行される。そして、ステップＳ１２４またはステップＳ１２５の後、上述したステップＳ３Ａ，Ｓ２，Ｓ９１〜Ｓ９３，Ｓ４Ｂ，Ｓ５，Ｓ６が順次実行され、ステップＳ９３において、誤差ΔＧがしきい値Ｇｔｈ以下であると判定されると、波源位置推定装置１Ｂの動作が終了する。

図４９は、実施の形態２による別の波源位置推定装置の概略図である。実施の形態２による波源位置推定装置は、図４９に示す波源位置推定装置１Ｃであってもよい。

図４９を参照して、波源位置推定装置１Ｃは、図３３に示す波源位置推定装置１Ｂの推定手段１４Ｂを推定手段１４Ｃに変えたものであり、その他は、波源位置推定装置１Ｂと同じである。

推定手段１４Ｃは、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかである場合、記憶手段１３Ａに記憶されたサンプルデータＤ_ＳＰＬ＿ｉに含まれる複数の受信電力Ｐ_１〜Ｐ_Ｎのうち誤差を抑制した複数の受信電力を検出し、その検出した複数の受信電力に対応する複数の位置情報を検出する。そして、推定手段１４Ｃは、その検出した複数の受信電力および複数の位置情報に基づいて、実施の形態１において説明した方法によって、電力分布を推定し、その推定結果である電力分布推定結果を生成する。そして、推定手段１４Ｃは、その生成した電力分布推定結果を用いて実施の形態１において説明した方法によって波源位置を推定する。推定手段１４Ｃは、その他、推定手段１４Ｂと同じ機能を果たす。

図５０は、実施の形態２における受信電力の等高線ＣＴＲを作成する方法を説明するための図である。図５０において、白三角は、電力値がＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力の測定点を示し、黒三角は、電力値がＰＷ未満であり、マスター局からの電波の受信電力の測定点を示し、菱形は、マスター局以外の通信装置からの電波の受信電力の測定点を示し、白丸は、受信電力の重心点Ｐ＿Ｇを示す。そして、白三角、黒三角および菱形の各々は、位置情報（ｘ_ｉ，ｙ_ｉ）および受信電力Ｐ_ｉからなる。

図５０の（ａ）を参照して、推定手段１４Ｃは、記憶手段１３Ａに記録された複数のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿Ｎに基づいて、複数のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿Ｎに含まれる複数の位置情報［ｘ_１，ｙ_１］〜［ｘ_Ｎ，ｙ_Ｎ］および複数の受信電力Ｐ_１〜Ｐ_Ｎの分布状態（図５０の（ａ）に示す分布状態）を取得する。

そして、推定手段１４Ｃは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷを決定する。また、推定手段１４Ｃは、推定手段１４Ａと同じように、電力値ＰＷと処理範囲の半径ｒとの関係を予め保持している。

推定手段１４Ｃは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷを決定すると、推定結果ＥＳＴがＴＤＤであるとき、複数の受信電力Ｐ_１〜Ｐ_Ｎのうち、電力値ＰＷ以上であり、かつ、マスター局からの電波の受信電力において受信電力の最大値Ｐ＿ＭＡＸを抽出する。即ち、推定手段１４Ｃは、白三角で示される測定点に含まれる受信電力Ｐのうち、最大値Ｐ＿ＭＡＸを抽出する。ここで、受信電力の最大値Ｐ＿ＭＡＸを抽出するのは、受信電力の最大値Ｐ＿ＭＡＸを有する端末装置２２が最も波源Ｓに近いと考えられるからである。

そして、推定手段１４Ｃは、−６０ｄＢｍの電力値ＰＷを有する受信電力の等高線ＣＴＲを作成する場合、受信電力の最大値Ｐ＿ＭＡＸを抽出すると、その抽出した受信電力の最大値Ｐ＿ＭＡＸを有する測定点を中心として半径ｒ＝５００ｍの範囲を処理範囲ＲＥＧ＿ＰＲＳ＿Ｍとして決定する。このように、推定手段１４Ｃは、電力値がＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力に基づいて処理範囲ＲＥＧ＿ＰＲＳ＿Ｍを決定する。

推定手段１４Ｃは、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍとして決定すると、その決定した処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ内のｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する。

図５０の（ｂ）を参照して、推定手段１４Ｃは、抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち、マスター局からの電波の受信電力を含むｓ（ｓは、１≦ｓ≦ｄを満たす整数）個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓを抽出する。そして、推定手段１４Ｃは、ｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに含まれるｓ個の受信電力Ｐ_１〜Ｐ_ｓおよびｓ個の位置情報［ｘ_１，ｙ_１］〜［ｘ_ｓ，ｙ_ｓ］基づいて、式（１）を用いて受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）を算出する。

即ち、推定手段１４Ｃは、ｓ個の受信電力Ｐ_１〜Ｐ_ｓによって重み付けされたｓ個の端末装置２２の位置の平均を演算することによって受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）を算出する。

そうすると、推定手段１４Ｃは、マスター局からの電波の受信電力のうち電力値ＰＷ以上の受信電力Ｐ_ｉを有し、かつ、受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）から最も遠い位置に存在する端末装置２２−１と受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）との間の距離を最大半径Ｒとして求める。そして、推定手段１４Ｃは、受信電力の重心点Ｐ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）を中心とし、かつ、最大半径Ｒを半径とする円形状を有する受信電力の等高線ＣＴＲを作成する。

このように、推定手段１４Ｃは、マスター局以外の通信装置からの電波の受信電力、および電力値ＰＷ未満の受信電力（黒三角）を除外し、電力値ＰＷ以上であり、かつ、マスター局からの電波の受信電力（白三角）のみに基づいて最大半径Ｒを決定する。

従って、推定手段１４Ｃによる受信電力の等高線ＣＴＲの作成方法によれば、誤差を抑制した受信電力を用いて受信電力の等高線ＣＴＲを作成できる。

なお、推定手段１４Ｃは、推定結果ＥＳＴがＣＤＭＡ変調またはＦＭ変調であるときも、図５０において説明した方法によって受信電力の等高線ＣＴＲを作成する。そして、推定結果ＥＳＴがＣＤＭＡ変調であるとき、上記の「マスター局からの電波の受信電力」を「フォワードリンクにおける電波の受信電力」と読み替えればよく、推定結果ＥＳＴがＦＭ変調であるとき、上記の「マスター局からの電波の受信電力」を「変調信号による周波数偏移を含む帯域幅を有するフィルタで検出した受信電力」と読み替えればよい。

また、推定手段１４Ｃは、通信方式または変調方式の推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調以外の場合、実施の形態１において説明した説明した方法によって受信電力の等高線ＣＴＲを作成する。

図５１は、実施の形態２におけるｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄの抽出方法を説明するための図である。図４６において、白三角は、電力値Ｐ以上、かつ、マスター局からの電波の受信電力を含む測定点を示し、黒三角は、電力値Ｐ未満、かつ、マスター局からの電波の受信電力を含む測定点を示し、菱形は、マスター局からの電波の受信電力を含む測定点を示す。なお、菱形の測定点は、電力値Ｐ以上の受信電力を含むか否かを識別していない。

図５１の（ａ）を参照して、推定手段１４Ｃは、上述した方法によって、電力値がＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力に基づいて処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ１を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ１内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する。そして、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち電力値ＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力である受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて、上述した方法によって、受信電力の等高線ＣＴＲ１を作成する。

図５１の（ｂ）を参照して、推定手段１４Ｃは、その後、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ１内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを除外して、上述した方法によって、電力値がＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力である受信電力に基づいて処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ２を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ２内に存在する新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する。そして、推定手段１４Ｃは、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち電力値ＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力である受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて、上述した方法によって、受信電力の等高線ＣＴＲ２を作成する。

そして、推定手段１４Ｃは、マスター局からの電波の受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力がなくなるまで、既に抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを除外しながら、上述した方法によって処理範囲ＲＥＧ＿ＰＲＳ＿Ｍを決定し、その決定した処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ内に存在する新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出することを繰り返し実行する。

そして、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出するごとに、上述した方法によって、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち電力値ＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力である受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて、上述した方法によって、受信電力の等高線ＣＴＲを作成する。

この場合、推定手段１４Ｃは、次の２つの方法のいずれかを用いて受信電力の等高線ＣＴＲを作成する。

（３）作成方法３
推定手段１４Ｃは、既に受信電力の等高線ＣＴＲの作成に用いたｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを用いずに、抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄに基づいて受信電力の等高線ＣＴＲを作成する。

（４）作成方法４
推定手段１４Ｃは、既に受信電力の等高線ＣＴＲの作成に用いたｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄと新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄとの両方に含まれる端末情報と、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄとに基づいて受信電力の等高線ＣＴＲを作成する。

図５２は、受信電力の等高線ＣＴＲの作成方法３を説明するための図である。図５２を参照して、サンプルデータＤ_ＳＰＬ＿１は、受信電力の最大値Ｐ＿ＭＡＸ１を含み、サンプルデータＤ_ＳＰＬ＿２は、サンプルデータＤ_ＳＰＬ＿１を除いた場合における受信電力の最大値Ｐ＿ＭＡＸ２を含む。また、サンプルデータＤ_ＳＰＬ＿３，Ｄ_ＳＰＬ＿４は、マスター局からの電波の受信電力であり、かつ、電力値ＰＷ以上である受信電力を含む。更に、白丸は、受信電力の重心点Ｐ＿Ｇを示す。

推定手段１４Ｃは、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ３を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ３内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿ＲＥＧ＿ＰＲＳ＿Ｍ３を抽出する。

そして、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿ＲＥＧ＿ＰＲＳ＿Ｍ３のうち電力値ＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力である受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて、上述した方法によって受信電力の等高線ＣＴＲ５を作成する。

その後、推定手段１４Ｃは、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ３内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿ＲＥＧ＿ＰＲＳ＿Ｍ３を除いて、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ４を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ４内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿ＲＥＧ＿ＰＲＳ＿Ｍ４を抽出する。

この場合、サンプルデータＤ_ＳＰＬ＿３，Ｄ_ＳＰＬ＿４は、ｄ個のサンプルデータＤ_ＳＰＬ＿ＲＥＧ＿ＰＲＳ＿Ｍ３およびｄ個のサンプルデータＤ_ＳＰＬ＿ＲＥＧ＿ＰＲＳ＿Ｍ４の両方に含まれる。

そして、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿ＲＥＧ＿ＰＲＳ＿Ｍ４およびサンプルデータＤ_ＳＰＬ＿３，Ｄ_ＳＰＬ＿Ｍ４のうちマスター局からの電波の受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓを抽出し、その抽出したｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて受信電力の重心点Ｐ＿Ｇ３を算出し、受信電力の重心点Ｐ＿Ｇ３とサンプルデータＤ_ＳＰＬ＿４との距離を半径とする円形状を有する受信電力の等高線ＣＴＲ６を作成する。

サンプルデータＤ_ＳＰＬ＿３，Ｄ_ＳＰＬ＿４を含めずに受信電力の等高線ＣＴＲを作成する場合、受信電力の重心点Ｐ＿Ｇは、受信電力の重心点Ｐ＿Ｇ４になり、受信電力の最大値Ｐ＿ＭＡＸ２を含むサンプルデータＤ_ＳＰＬ＿２の位置からのずれが大きくなり、受信電力の等高線ＣＴＲ６よりも半径が小さい円形状を有する受信電力の等高線ＣＴＲになる。その結果、精度良く受信電力の等高線ＣＴＲを作成することが困難になる。

しかし、サンプルデータＤ_ＳＰＬ＿３，Ｄ_ＳＰＬ＿４を含めて受信電力の等高線ＣＴＲを作成することにより、受信電力の重心点Ｐ＿Ｇ３は、受信電力の最大値Ｐ＿ＭＡＸ２を含むサンプルデータＤ_ＳＰＬ＿２の位置の近傍に存在し、より大きい半径を有する受信電力の等高線ＣＴＲ６が作成される。その結果、精度良く受信電力の等高線ＣＴＲを作成することができる。

図５３は、実施の形態２におけるｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄの別の抽出方法を説明するための図である。図５３を参照して、受信電力の最大値Ｐ＿ＭＡＸ６は、受信電力の最大値Ｐ＿ＭＡＸ５よりも小さい。

推定手段１４Ｃは、上述した方法によって、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ５を決定し、処理範囲ＲＧ＿ＰＲＳ＿Ｍ５内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する（図５３の（ａ）参照）。

その後、推定手段１４Ｃは、既に抽出した処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ５内のｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち、マスター局からの電波の受信電力であり、かつ、受信電力の最大値Ｐ＿ＭＡＸからなる受信電力を含むサンプルデータＤ_ＳＰＬのみを除外して処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ６を決定し、処理範囲ＲＥＧ＿ＰＲＳ＿Ｍ６内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する（図５３の（ｂ）参照）。

このように、推定手段１４Ｃは、マスター局からの電波の受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力がなくなるまで、既に抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち、マスター局からの電波の受信電力であり、かつ、受信電力の最大値Ｐ＿ＭＡＸからなる受信電力を含むサンプルデータＤ_ＳＰＬのみを除外しながら新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを繰り返し抽出する。

そして、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出するごとに、その抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち、電力値ＰＷ（＝受信電力の等高線ＣＴＲが有する電力値）以上であり、かつ、マスター局からの電波の受信電力である受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて、上述した方法によって受信電力の等高線ＣＴＲを作成する。

上述したように、推定手段１４Ｃは、既に抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち、マスター局からの電波の受信電力であり、かつ、受信電力の最大値Ｐ＿ＭＡＸからなる受信電力を含むサンプルデータＤ_ＳＰＬを少なくとも除きながらｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出することを、マスター局からの電波の受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力がなくなるまで繰り返し行う。

なお、推定手段１４Ｃは、通信方式または変調方式の推定結果ＥＳＴがＣＤＭＡ変調またはＦＭ変調である場合、図４６〜図４８において説明した方法によってｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出するとともに、その抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄに含まれるｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて受信電力の等高線ＣＴＲを作成する。そして、推定結果ＥＳＴがＣＤＭＡ変調であるとき、上記の「マスター局からの電波の受信電力」を「フォワードリンクにおける電波の受信電力」と読み替えればよく、推定結果ＥＳＴがＦＭ変調であるとき、上記の「マスター局からの電波の受信電力」を「変調信号による周波数偏移を含む帯域幅を有するフィルタで検出した受信電力」と読み替えればよい。

図５４は、実施の形態２におけるｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄの抽出方式と受信電力の等高線ＣＴＲの作成処理との関係を示す図である。

図５４を参照して、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄの抽出方式は、誤差を抑制した受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力がなくなるまで、既に抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを除外しながら新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する抽出方式３と、誤差を抑制した受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力がなくなるまで、既に抽出したｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち、受信電力の最大値Ｐ＿ＭＡＸを含むサンプルデータＤ_ＳＰＬのみを除外しながら新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する抽出方式４とを含む。

また、受信電力の等高線ＣＴＲの作成処理は、抽出されたｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち誤差を抑制した受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて受信電力の重心点Ｐ＿Ｇを算出し、誤差を抑制した受信電力のうち電力値ＰＷ以上の受信電力を有し、かつ、受信電力の重心点Ｐ＿Ｇから最も遠い位置に存在する端末装置と受信電力の重心点Ｐ＿Ｇとの距離を最大距離として求め、受信電力の重心点Ｐ＿Ｇを中心とし、かつ、最大距離を半径とする円形状を有する受信電力の等高線ＣＴＲを作成する等高線作成処理３と、既に抽出されたｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄと新たに抽出されたｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄとの両方に含まれるサンプルデータＤ_ＳＰＬと、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄとのうち誤差を抑制した受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて受信電力の重心点Ｐ＿Ｇを算出し、誤差を抑制した受信電力のうち電力値ＰＷ以上の受信電力を有し、かつ、受信電力の重心点Ｐ＿Ｇから最も遠い位置に存在する端末装置と受信電力の重心点Ｐ＿Ｇとの距離を最大距離として求め、受信電力の重心点Ｐ＿Ｇを中心とし、かつ、最大距離を半径とする円形状を有する受信電力の等高線ＣＴＲを作成する等高線作成処理４とを含む。

そして、推定手段１４Ｃは、抽出方式３に従ってｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する場合、等高線作成処理３に従って受信電力の等高線ＣＴＲを作成し、または等高線作成処理４に従って受信電力の等高線ＣＴＲを作成する。

また、推定手段１４Ｃは、抽出方式４に従ってｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する場合、等高線作成処理３に従って受信電力の等高線ＣＴＲを作成する。

なお、抽出方式３，４と等高線作成処理３，４との組み合わせによって等高線を作成するときの計算量についての説明は、図２３における説明と同じである。

また、推定手段１４Ｃは、通信方式または変調方式がＴＤＤ、ＣＤＭＡ変調およびＦＭ変調以外である場合、上述した図２３に示す抽出方式１，２と等高線作成処理１，２との組み合わせを用いて等高線を作成する。

更に、推定手段１４Ｃは、上述した方法によって作成した２つの等高線ＣＴＲを図２４において説明した方法によって合成し、合成した等高線ＣＴＲ＿ＳＹＮを作成する。

更に、推定手段１４Ｃは、電波の受信感度以下の領域においては、上述した方法によって作成した等高線ＣＴＲに基づいて、図２５において説明した方法によって等高線ＣＴＲを作成する。

図５５は、実施の形態２における電力分布推定結果を作成する方法を説明するためのフローチャートである。

図５５を参照して、電力分布を推定する動作が開始されると、受信手段１２Ａは、一定時間におけるサンプルデータＤ_ＳＰＬ＿ｉを複数の端末装置２２から受信し、その受信したサンプルデータＤ_ＳＰＬ＿ｉを記憶手段１３Ａに格納する。即ち、波源位置推定装置１Ｃは、一定時間におけるサンプルデータＤ_ＳＰＬ＿ｉを取得する（ステップＳ１２１）。

そして、通信推定手段１５は、記憶手段１３Ａに記憶されたサンプルデータＤ_ＳＰＬ＿ｉに基づいて、上述した方法によって、共用周波数を用いる１次利用者および２次利用者の通信方式または変調方式を推定し（ステップＳ１２２）、その推定した推定結果ＥＳＴを推定手段１４Ｃへ出力する。

推定手段１４Ｃは、推定結果ＥＳＴを通信推定手段１５から受け、その受けた推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかであるか否かを判定する（ステップＳ１２３）。

ステップＳ１２３において、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかであると判定されたとき、推定手段１４Ｃは、抽出方式３または抽出方式４によって少なくとも１組のｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを処理範囲内から抽出する（ステップＳ１２６）。

その後、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出するごとに、サンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄのうち、誤差を抑制した受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに基づいて、上述した方法によって受信電力の等高線ＣＴＲを作成する（ステップＳ１２７）。

一方、ステップＳ１２３において、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれでもないと判定されたとき、推定手段１４Ｃは、抽出方式１または抽出方式２によって少なくとも１組のｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを処理範囲内から抽出する（ステップＳ１２８）。

その後、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出するごとに、サンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄに基づいて、上述した方法によって受信電力の等高線ＣＴＲを作成する（ステップＳ１２９）。

そして、ステップＳ１２７またはステップＳ１２９の後、電力分布を推定する動作が終了する。

図５６は、図５５に示すステップＳ１２６の詳細な動作を説明するためのフローチャートである。図５６を参照して、図５５のステップＳ１２３において、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれかであると判定されたとき、推定手段１４Ｃは、複数のサンプルデータＤ_ＳＰＬ＿ｉに基づいて、誤差を抑制した受信電力の最大値を検出する（ステップＳ１２６１）。

そして、推定手段１４Ｃは、受信電力の最大値が、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上であるか否かを判定する（ステップＳ１２６２）。

ステップＳ１２６２において、受信電力の最大値が電力値ＰＷ以上であると判定されたとき、推定手段１４Ｃは、受信電力の最大値を中心した半径ｒの処理範囲を決定する（ステップＳ１２６３）。その後、推定手段１４Ｃは、処理範囲内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する（ステップＳ１２６４）。

そして、ステップＳ１２６２において、受信電力の最大値が電力値ＰＷ以上でないと判定されたとき、またはステップＳ１２６４の後、推定手段１４Ｃは、誤差を抑制した受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力が無いか否かを判定する（ステップＳ１２６５）。

ステップＳ１２６５において、誤差を抑制した受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力があると判定されたとき、一連の動作は、ステップＳ１２６１へ戻り、ステップＳ１２６５において、誤差を抑制した受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力が無いと判定されるまで、ステップＳ１２６１〜ステップＳ１２６５が繰り返し実行される。

そして、ステップＳ１２６５において、誤差を抑制した受信電力であり、かつ、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力が無いと判定されると、一連の動作は、図５５のステップＳ１２７へ移行する。

図５７は、図５５に示すステップＳ１２８の詳細な動作を説明するためのフローチャートである。図５７を参照して、図５５のステップＳ１２３において、推定結果ＥＳＴがＴＤＤ、ＣＤＭＡ変調およびＦＭ変調のいずれでもないと判定されたとき、推定手段１４Ｃは、複数のサンプルデータＤ_ＳＰＬ＿ｉに基づいて受信電力の最大値を検出する（ステップＳ１２８１）。

そして、推定手段１４Ｃは、受信電力の最大値が、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上であるか否かを判定する（ステップＳ１２８２）。

ステップＳ１２８２において、受信電力の最大値が電力値ＰＷ以上であると判定されたとき、推定手段１４Ｃは、受信電力の最大値を中心した半径ｒの処理範囲を決定する（ステップＳ１２８３）。その後、推定手段１４Ｃは、処理範囲内に存在するｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出する（ステップＳ１２８４）。

そして、ステップＳ１２８２において、受信電力の最大値が電力値ＰＷ以上でないと判定されたとき、またはステップＳ１２８４の後、推定手段１４Ｃは、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力が無いか否かを判定する（ステップＳ１２８５）。

ステップＳ１２８５において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力があると判定されたとき、一連の動作は、ステップＳ１２８１へ戻り、ステップＳ１２８５において、作成すべき受信電力の等高線ＣＴＲの電力値ＰＷ以上である受信電力が無いと判定されるまで、ステップＳ１２８１〜ステップＳ１２８５が繰り返し実行される。

図５８は、図５５に示すステップＳ１２７の詳細な動作を説明するためのフローチャートである。図５８を参照して、図５５のステップＳ１２６の後、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄに基づいて、誤差を抑制した受信電力を含むｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓを抽出する（ステップＳ１２７１）。

そして、推定手段１４Ｃは、ｓ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｓに含まれるｓ個の受信電力およびｓ個の位置情報に基づいて受信電力の重心点Ｐ＿Ｇを演算する（ステップＳ１２７２）。

その後、推定手段１４Ｃは、誤差を抑制した受信電力のうち電力値ＰＷ（作成すべき受信電力の等高線ＣＴＲの電力値）以上の受信電力を有し、かつ、受信電力の重心点Ｐ＿Ｇから最も遠い位置に存在する端末装置と受信電力の重心点Ｐ＿Ｇとの距離を最大距離として求める（ステップＳ１２７３）。

引き続いて、推定手段１４Ｃは、重心点Ｐ＿Ｇを中心とし、最大距離を半径とする円形状を有する等高線ＣＴＲを作成する（ステップＳ１２７４）。

そうすると、推定手段１４Ｃは、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出したか否かを判定する（ステップＳ１２７５）。

ステップＳ１２７５において、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出したと判定されたとき、一連の動作は、ステップＳ１２７１へ移行する。その後、ステップＳ１２７５において、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出しなかったと判定されるまで、ステップＳ１２７１〜ステップＳ１２７５が繰り返し実行される。そして、ステップＳ１２７５において、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出しなかったと判定されると、一連の動作は、図５５の“終了”へ移行する。

図５９は、図５５に示すステップＳ１２９の詳細な動作を説明するためのフローチャートである。図５９を参照して、図５５のステップＳ１２８の後、推定手段１４Ｃは、ｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄに含まれるｄ個の受信電力およびｄ個の位置情報に基づいて受信電力の重心点Ｐ＿Ｇを演算する（ステップＳ１２９１）。

そして、推定手段１４Ｃは、電力値ＰＷ（作成すべき受信電力の等高線ＣＴＲの電力値）以上の受信電力を有し、かつ、受信電力の重心点Ｐ＿Ｇから最も遠い位置に存在する端末装置と受信電力の重心点Ｐ＿Ｇとの距離を最大距離として求める（ステップＳ１２９２）。

その後、推定手段１４Ｃは、重心点Ｐ＿Ｇを中心とし、最大距離を半径とする円形状を有する等高線ＣＴＲを作成する（ステップＳ１２９３）。

そうすると、推定手段１４Ｃは、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出したか否かを判定する（ステップＳ１２９４）。

ステップＳ１２９４において、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出したと判定されたとき、一連の動作は、ステップＳ１２９１へ移行する。その後、ステップＳ１２９４において、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出しなかったと判定されるまで、ステップＳ１２９１〜ステップＳ１２９４が繰り返し実行される。そして、ステップＳ１２９４において、新たなｄ個のサンプルデータＤ_ＳＰＬ＿１〜Ｄ_ＳＰＬ＿ｄを抽出しなかったと判定されると、一連の動作は、図５５の“終了”へ移行する。

図６０は、図４９に示す波源位置推定装置１Ｃの動作を説明するためのフローチャートである。

図６０に示すフローチャートは、図２９に示すフローチャートのステップＳ１１０をステップＳ１１０Ａに変えたものであり、その他は、図２９に示すフローチャートと同じである。

図６０を参照して、波源位置推定装置１Ｃの動作が開始されると、上述したステップＳ１〜ステップＳ４が順次実行される。

そして、ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａでないと判定されたとき、推定手段１４Ｃは、図５５に示すフローチャート（図５６〜図５９に示すフローチャートを含む）に従って電力分布推定結果を作成する（ステップＳ１１０Ａ）。

その後、推定手段１４Ｃは、電力分布推定結果を用いて、上述したステップＳ５，Ｓ６を順次実行する。そして、ステップＳ４において、ｘ_Ａ＜Δｘ_Ａ、かつ、ｙ_Ａ＜Δｙ_Ａであると判定されたとき、波源位置推定装置１Ｃの動作が終了する。

図６１は、図４９に示す波源位置推定装置１Ｃの動作を説明するための別のフローチャートである。

図６１に示すフローチャートは、図３０に示すフローチャートのステップＳ１１０をステップＳ１１０Ａに変えたものであり、その他は、図３０に示すフローチャートと同じである。

図６１を参照して、波源位置推定装置１Ｃの動作が開始されると、上述したステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ４Ａ，Ｓ４Ｂが順次実行される。

そして、ステップＳ４Ｂの後、推定手段１４Ｃは、上述したステップＳ１１０Ａ，Ｓ５，Ｓ６を順次実行し、ステップＳ４Ａにおいて、ｋ＝Ｋであると判定されると、波源位置推定装置１Ｃの動作が終了する。

図６２は、図４９に示す波源位置推定装置１Ｃの動作を説明するための更に別のフローチャートである。

図６２に示すフローチャートは、図３１に示すフローチャートのステップＳ１１０をステップＳ１１０Ａに変えたものであり、その他は、図３１に示すフローチャートと同じである。

図６２を参照して、波源位置推定装置１Ｃの動作が開始されると、上述したステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ９１〜Ｓ９３，Ｓ４Ｂが順次実行される。

そして、ステップＳ４Ｂの後、推定手段１４Ｃは、上述したステップＳ１１０Ａ，Ｓ５，Ｓ６を順次実行し、ステップＳ９３において、誤差ΔＧがしきい値Ｇｔｈ以下であると判定されると、波源位置推定装置１Ｃの動作が終了する。

このように、波源位置推定装置１Ｃは、波源位置の推定が収束条件を満たさないと判定された場合、収束条件を満たすと判定されるまで、電力分布を推定し、その推定した電力分布を用いて仮想端末装置２ｖｊの追加と受信電力Ｐ_ｉの信頼度ＲＥＬ_ｉが低い端末装置２２の削除とを行い、波源位置を推定する処理を繰り返し行う。

なお、実施の形態２においては、波源位置推定装置１Ｂの動作は、ソフトウェアによって実行されてもよい。この場合、波源位置推定装置１Ｂは、ＣＰＵ、ＲＯＭおよびＲＡＭを備える。

ＲＯＭは、図４０に示すフローチャートのステップＳ１２１〜Ｓ１２５，Ｓ１〜Ｓ６（図４１に示すフローチャート、図４２に示すフローチャート、図４３に示すフローチャート、図４４に示すフローチャートおよび図４６に示すフローチャートを含む）を備えるプログラムＰｒｏｇ＿Ｇ、または図４０に示すフローチャートのステップＳ１２１〜Ｓ１２５，Ｓ１〜Ｓ６（図４５に示すフローチャートおよび図４６に示すフローチャートを含む）を備えるプログラムＰｒｏｇ＿Ｈ、または図４７に示すフローチャート（図４１に示すフローチャート、図４２に示すフローチャート、図４３に示すフローチャート、図４４に示すフローチャートおよび図４６に示すフローチャートを含む）を備えるプログラムＰｒｏｇ＿Ｉ、または図４７に示すフローチャート（図４５に示すフローチャートおよび図４６に示すフローチャートを含む）を備えるプログラムＰｒｏｇ＿Ｊ、または図４８に示すフローチャート（図４１に示すフローチャート、図４２に示すフローチャート、図４３に示すフローチャート、図４４に示すフローチャートおよび図４６に示すフローチャートを含む）を備えるプログラムＰｒｏｇ＿Ｋ、または図４８に示すフローチャート（図４５に示すフローチャートおよび図４６に示すフローチャートを含む）を備えるプログラムＰｒｏｇ＿Ｌを格納する。ＲＡＭは、各種演算の結果を一時的に格納する。

そして、ＣＰＵは、ＲＯＭからプログラムＰｒｏｇ＿Ｇ〜Ｐｒｏｇ＿Ｌのいずれかを読み出して実行し、上述した方法によって、波源位置を推定する。

また、実施の形態２においては、波源位置推定装置１Ｃの動作は、ソフトウェアによって実行されてもよい。この場合、波源位置推定装置１Ｃは、ＣＰＵ、ＲＯＭおよびＲＡＭを備える。

ＲＯＭは、図５５に示すフローチャートのステップＳ１２１〜Ｓ１２３，Ｓ１２６〜Ｓ１２９（図５６〜図５９に示すフローチャートを含む）を備えるプログラムＰｒｏｇ＿Ｍ、または図６０に示すフローチャートのステップＳ１〜Ｓ４，Ｓ１１０Ａ，Ｓ５，Ｓ６（図５０に示すフローチャート（図５６〜図５９に示すフローチャートを含む）を含む）を備えるプログラムＰｒｏｇ＿Ｎ、または図６１に示すフローチャートのステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ４Ａ，Ｓ４Ｂ，Ｓ１１０Ａ，Ｓ５，Ｓ６（図５０に示すフローチャート（図５６〜図５９に示すフローチャートを含む）を含む）を備えるプログラムＰｒｏｇ＿Ｏ、または図６２に示すフローチャートのステップＳ１，Ｓ３Ａ，Ｓ２，Ｓ９１〜Ｓ９３，Ｓ４Ｂ，Ｓ１１０Ａ，Ｓ５，Ｓ６（図５０に示すフローチャート（図５６〜図５９に示すフローチャートを含む）を含む）を備えるプログラムＰｒｏｇ＿Ｐを格納する。ＲＡＭは、各種演算の結果を一時的に格納する。

そして、ＣＰＵは、ＲＯＭからプログラムＰｒｏｇ＿Ｍを読み出して実行し、上述した方法によって、電力分布を推定する。また、ＣＰＵは、ＲＯＭからプログラムＰｒｏｇ＿Ｎ〜Ｐｒｏｇ＿Ｐのいずれかを読み出して実行し、上述した方法によって、電力分布を推定し、その推定した電力分布を用いて波源位置を推定する。

また、プログラムＰｒｏｇ＿Ｇ〜Ｐｒｏｇ＿Ｐは、ＣＤ，ＤＶＤ等の記録媒体に記録されて流通してもよい。

プログラムＰｒｏｇ＿Ｇ〜Ｐｒｏｇ＿Ｐのいずれかを記録した記録媒体をパーソナルコンピュータに装着すると、ＣＰＵは、記録媒体からプログラムＰｒｏｇ＿Ｇ〜Ｐｒｏｇ＿Ｐのいずれかを読み出して実行し、上述した方法によって、波源位置を推定する。

従って、プログラムＰｒｏｇ＿Ｇ〜Ｐｒｏｇ＿Ｐのいずれかを記録した記録媒体は、コンピュータ（ＣＰＵ）読み取り可能な記録媒体である。

実施の形態２におけるその他の説明は、実施の形態１における説明と同じである。

上述した実施の形態によれば、この発明の実施の形態による波源位置推定装置は、移動端末または静止端末からなる実端末装置の位置情報と実端末装置における受信電力とを含む実端末情報を複数の実端末装置から受信する受信手段と、複数の実端末装置の複数の実端末情報に含まれる複数の位置情報と複数の受信電力とに基づいて受信電力によって重み付けされた複数の実端末装置の重心を波源位置として推定し、その推定した波源位置に誤差があると判定したとき、少なくとも１つの仮想端末装置の位置を示す位置情報と受信電力とを追加して複数の実端末装置および少なくとも１つの仮想端末装置の位置情報および受信電力に基づいて受信電力によって重み付けされた複数の実端末装置および少なくとも１つの仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する推定手段とを備えていればよい。

また、この発明の実施の形態によるプログラムは、受信手段が、移動端末または静止端末からなる実端末装置の位置情報と実端末装置における受信電力とを含む実端末情報を複数の実端末装置から受信する第１のステップと、推定手段が、複数の実端末装置の複数の実端末情報に含まれる複数の位置情報と複数の受信電力とに基づいて受信電力によって重み付けされた複数の実端末装置の重心を波源位置として推定し、その推定した波源位置に誤差があると判定したとき、少なくとも１つの仮想端末装置の位置を示す位置情報と受信電力とを追加して複数の実端末装置および少なくとも１つの仮想端末装置の位置情報および受信電力に基づいて受信電力によって重み付けされた複数の実端末装置および少なくとも１つの仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する第２のステップとをコンピュータに実行させるものであればよい。

今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施の形態の説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。

この発明は、波源位置推定装置、コンピュータに実行させるためのプログラム、およびプログラムを記録したコンピュータ読み取り可能な記録媒体に適用される。

１，１Ａ，１Ｂ，１Ｃ 波源位置推定装置、２，２２ 端末装置、１０，１０Ａ 無線通信システム、１１ アンテナ、１２，１２Ａ 受信手段、１３，１３Ａ 記憶手段、１４，１４Ａ，１４Ｂ，１４Ｃ 推定手段、１５ 通信推定手段、３１，３２ 基地局。

Claims (27)

1. 移動端末または静止端末からなる実端末装置の位置情報と前記実端末装置における受信電力とを含む実端末情報を複数の実端末装置から受信する受信手段と、
   前記複数の実端末装置の複数の実端末情報に含まれる複数の位置情報と複数の受信電力とに基づいて前記受信電力によって重み付けされた前記複数の実端末装置の重心を波源位置として推定し、その推定した波源位置に誤差があると判定したとき、少なくとも１つの仮想端末装置の位置を示す位置情報と受信電力とを追加して前記複数の実端末装置および前記少なくとも１つの仮想端末装置の位置情報および受信電力に基づいて前記受信電力によって重み付けされた前記複数の実端末装置および前記少なくとも１つの仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する推定手段とを備える波源位置推定装置。
2. 前記推定手段は、電力分布の推定結果を示す電力分布推定結果に基づいて前記少なくとも１つの仮想端末装置の位置情報と受信電力とを追加し、その追加した前記少なくとも１つの仮想端末装置の位置情報と受信電力と、前記複数の実端末装置の複数の位置情報と複数の受信電力とに基づいて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項１に記載の波源位置推定装置。
3. 前記推定手段は、前記少なくとも１つの仮想端末装置をランダムに配置し、その配置した位置を示す位置情報を前記少なくとも１つの仮想端末装置の位置情報とし、前記電力分布推定結果に基づいて、前記少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を前記少なくとも１つの仮想端末装置の受信電力とする、請求項２に記載の波源位置推定装置。
4. 前記推定手段は、前記電力分布推定結果によって示される等電力線上に前記少なくとも１つの仮想端末装置を配置し、前記等電力線上の電力値を前記少なくとも１つの仮想端末装置の受信電力とし、前記少なくとも１つの仮想端末装置を配置した位置を示す位置情報を前記少なくとも１つの仮想端末装置の位置情報とする、請求項２に記載の波源位置推定装置。
5. 前記推定手段は、複数の実端末装置のうち、隣接する実端末装置間の距離を求め、その求めた距離が第１のしきい値以下である実端末装置が配置されている領域以外の領域に前記少なくとも１つの仮想端末装置を配置し、その配置した位置を示す位置情報を前記少なくとも１つの仮想端末装置の位置情報とし、前記電力分布推定結果に基づいて、前記少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を前記少なくとも１つの仮想端末装置の受信電力とする、請求項２に記載の波源位置推定装置。
6. 前記推定手段は、前記電力分布推定結果によって示される等電力線上の電力値と前記実端末装置における受信電力との電力差が第２のしきい値よりも大きいか否かを前記複数の実端末装置の全てについて判定し、前記電力差が前記第２のしきい値よりも大きい実端末装置を除外して前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項２から請求項５のいずれか１項に記載の波源位置推定装置。
7. 前記推定手段は、前記複数の実端末装置および前記少なくとも１つの仮想端末装置の位置の平均が第３のしきい値以下であるとき、前記収束条件を満たすと判定し、前記推定処理を完了する、請求項１から請求項６のいずれか１項に記載の波源位置推定装置。
8. 前記推定手段は、前記推定処理の実行回数が所定回数に達したとき、前記収束条件を満たすと判定し、前記推定処理を完了する、請求項１から請求項６のいずれか１項に記載の波源位置推定装置。
9. 前記推定手段は、前記推定処理を２回以上実行したとき、前記推定処理をｋ−１（ｋは２以上の整数）回実行したときに推定された波源位置と、前記推定処理をｋ回実行したときに推定された波源位置との誤差が第４のしきい値以下であるとき、前記収束条件を満たすと判定し、前記推定処理を完了する、請求項１から請求項６のいずれか１項に記載の波源位置推定装置。
10. １次利用者および２次利用者の通信方式または変調方式を推定する通信推定手段を更に備え、
    前記推定手段は、前記通信推定手段による推定結果に基づいて、前記複数の実端末装置の複数の受信電力のうち誤差を抑制した複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項１から請求項９のいずれか１項に記載の波源位置推定装置。
11. 前記推定手段は、前記通信推定手段によって推定された通信方式が時分割複信であるとき、前記複数の実端末装置の複数の受信電力のうち、他の通信装置を制御および／または管理する通信装置から送信された電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項１０に記載の波源位置推定装置。
12. 前記推定手段は、前記通信推定手段によって推定された通信方式が周波数分割複信であるとき、前記複数の実端末装置の複数の受信電力のうち、上りリンクにおける電波の複数の受信電力を用いず、下りリンクにおける電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項１０に記載の波源位置推定装置。
13. 前記推定手段は、前記通信推定手段によって推定された変調方式がＦＭ変調であるとき、前記複数の実端末装置の複数の受信電力のうち、変調信号による周波数偏移を含む帯域幅を有するフィルタを用いて検出された複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項１０に記載の波源位置推定装置。
14. 受信手段が、移動端末または静止端末からなる実端末装置の位置情報と前記実端末装置における受信電力とを含む実端末情報を複数の実端末装置から受信する第１のステップと、
    推定手段が、前記複数の実端末装置の複数の実端末情報に含まれる複数の位置情報と複数の受信電力とに基づいて前記受信電力によって重み付けされた前記複数の実端末装置の重心を波源位置として推定し、その推定した波源位置に誤差があると判定したとき、少なくとも１つの仮想端末装置の位置を示す位置情報と受信電力とを追加して前記複数の実端末装置および前記少なくとも１つの仮想端末装置の位置情報および受信電力に基づいて前記受信電力によって重み付けされた前記複数の実端末装置および前記少なくとも１つの仮想端末装置の重心を波源位置として推定する推定処理を収束条件が満たされるまで繰り返し実行する第２のステップとをコンピュータに実行させるためのプログラム。
15. 前記推定手段は、前記第２のステップにおいて、電力分布の推定結果を示す電力分布推定結果に基づいて前記少なくとも１つの仮想端末装置の位置情報と受信電力とを追加し、その追加した前記少なくとも１つの仮想端末装置の位置情報と受信電力と、前記複数の実端末装置の複数の位置情報と複数の受信電力とに基づいて前記推定処理を収束条件が満たされるまで繰り返し実行する、請求項１４に記載のコンピュータに実行させるためのプログラム。
16. 前記推定手段は、前記第２のステップにおいて、前記少なくとも１つの仮想端末装置をランダムに配置し、その配置した位置を示す位置情報を前記少なくとも１つの仮想端末装置の位置情報とし、前記電力分布推定結果に基づいて、前記少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を前記少なくとも１つの仮想端末装置の受信電力とする、請求項１５に記載のコンピュータに実行させるためのプログラム。
17. 前記推定手段は、前記第２のステップにおいて、前記電力分布推定結果によって示される等電力線上に前記少なくとも１つの仮想端末装置を配置し、前記等電力線上の電力値を前記少なくとも１つの仮想端末装置の受信電力とし、前記少なくとも１つの仮想端末装置を配置した位置を示す位置情報を前記少なくとも１つの仮想端末装置の位置情報とする、請求項１５に記載のコンピュータに実行させるためのプログラム。
18. 前記推定手段は、前記第２のステップにおいて、複数の実端末装置のうち、隣接する実端末装置間の距離を求め、その求めた距離が第１のしきい値以下である実端末装置が配置されている領域以外の領域に前記少なくとも１つの仮想端末装置を配置し、その配置した位置を示す位置情報を前記少なくとも１つの仮想端末装置の位置情報とし、前記電力分布推定結果に基づいて、前記少なくとも１つの仮想端末装置を配置した位置における受信電力を求め、その求めた受信電力を前記少なくとも１つの仮想端末装置の受信電力とする、請求項１５に記載のコンピュータに実行させるためのプログラム。
19. 前記推定手段は、前記第２のステップにおいて、前記電力分布推定結果によって示される等電力線上の電力値と前記実端末装置における受信電力との電力差が第２のしきい値よりも大きいか否かを前記複数の実端末装置の全てについて判定し、前記電力差が前記第２のしきい値よりも大きい実端末装置を除外して前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項１５から請求項１８のいずれか１項に記載のコンピュータに実行させるためのプログラム。
20. 前記推定手段は、前記第２のステップにおいて、前記複数の実端末装置および前記少なくとも１つの仮想端末装置の位置の平均が第３のしきい値以下であるとき、前記収束条件を満たすと判定し、前記推定処理を完了する、請求項１４から請求項１９のいずれか１項に記載のコンピュータに実行させるためのプログラム。
21. 前記推定手段は、前記第２のステップにおいて、前記推定処理の実行回数が所定回数に達したとき、前記収束条件を満たすと判定し、前記推定処理を完了する、請求項１４から請求項１９のいずれか１項に記載のコンピュータに実行させるためのプログラム。
22. 前記推定手段は、前記第２のステップにおいて、前記推定処理を２回以上実行したとき、前記推定処理をｊ（ｊは２以上の整数）回実行したときに推定された波源位置と、前記推定処理をｊ＋１回実行したときに推定された波源位置との誤差が第４のしきい値以下であるとき、前記収束条件を満たすと判定し、前記推定処理を完了する、請求項１４から請求項１９のいずれか１項に記載のコンピュータに実行させるためのプログラム。
23. 通信推定手段が、１次利用者および２次利用者の通信方式または変調方式を推定する第３のステップを更にコンピュータに実行させ、
    前記推定手段は、前記第２のステップにおいて、前記通信推定手段による推定結果に基づいて、前記複数の実端末装置の複数の受信電力のうち誤差を抑制した複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項１４から請求項２２のいずれか１項に記載のコンピュータに実行させるためのプログラム。
24. 前記推定手段は、前記第２のステップにおいて、前記通信推定手段によって推定された通信方式が時分割複信であるとき、前記複数の実端末装置の複数の受信電力のうち、他の通信装置を制御および／または管理する通信装置から送信された電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項２３に記載のコンピュータに実行させるためのプログラム。
25. 前記推定手段は、前記第２のステップにおいて、前記通信推定手段によって推定された通信方式が周波数分割複信であるとき、前記複数の実端末装置の複数の受信電力のうち、上りリンクにおける電波の複数の受信電力を用いず、下りリンクにおける電波の複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項２３に記載のコンピュータに実行させるためのプログラム。
26. 前記推定手段は、前記第２のステップにおいて、前記通信推定手段によって推定された変調方式がＦＭ変調であるとき、前記複数の実端末装置の複数の受信電力のうち、変調信号による周波数偏移を含む帯域幅を有するフィルタを用いて検出された複数の受信電力を選択し、その選択した複数の受信電力を用いて前記推定処理を前記収束条件が満たされるまで繰り返し実行する、請求項２３に記載のコンピュータに実行させるためのプログラム。
27. 請求項１４から請求項２６のいずれか１項に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。

Images