JP2019197957A - 電力推定装置、それを備えた無線通信システム、コンピュータに実行させるためのプログラムおよびデータ構造 - Google Patents

Abstract

【課題】受信電力の３次元的な電力分布を推定可能な電力推定装置を提供する。【解決手段】抽出手段１３は、端末装置の位置情報と電波の受信電力とを含む複数のモニター情報から受信電力の最大値を含むｋ個のモニター情報を抽出する。処理手段１４は、ｋ個のモニター情報に基づいて、受信電力の重心点Ｇと、重心点Ｇとの最大距離とを求める処理をｍ個の高さについて実行する。演算手段１５は、ｍ個の重心点Ｇ、ｍ個の最大距離及びｍ個の高さに基づいて重心点Ｇと高さとの関係を示す第１の関係式と、最大距離と高さとの関係を示す第２の関係式とを演算する。作成手段１６は、第１及び第２の関係式に基づいて、任意の高さにおける重心点Ｇ及び最大距離を求め、その求めた重心点Ｇを中心とし、最大距離を半径とした等高線を作成する。【選択図】図２

Description

この発明は、電力推定装置、それを備えた無線通信システム、コンピュータに実行させるためのプログラムおよびデータ構造に関する。

移動通信に適した６ＧＨｚ以下の周波数帯域に対して、周波数を共用することによって、第５世代移動通信システムに必要な周波数帯域を新たに創出する研究がなされている。

その一環として、既に幅広いエリアに分布している移動端末をモニター端末とし、取得したモニター端末の位置および電波の受信電力の情報から、場所ごとの既存システムにおける電波の受信電力を求め、これらを解析することで、周波数共用の条件を決定する移動分散モニタリングに着目して研究が進められている。

移動分散モニタリングでは、モニター端末の分布にバラツキがあるため、モニター端末が存在しない場所の受信電力は、周辺の受信電力などから推定する必要がある。

受信電力を推定する方式として、モニター端末が取得したモニター情報（モニター端末の位置および電波の受信電力の情報）から受信電力の等高線を作成し、その作成した受信電力の等高線から場所ごとの受信電力を推定する方法が考えられる。

従来、受信電力の等高線を作成する方法として、非特許文献１に記載の方法が知られている。非特許文献１に記載された受信電力の等高線の作成方法は、２つのモニター端末における２つの受信電力に基づいて受信電力のベクトルを作成し、その作成したベクトルと直交する勾配（Gradient）曲線を作成して受信電力の等高線とする方法である。

また、端末装置側で、フェージング等により低下した受信電力の精度を補償する方法として非特許文献２に記載の方法が知られている。

この方法は、信号の周期性が既知である場合、実際に受信した受信信号と、フェージング等の影響を受けた場合の既知の受信信号との周期性を比較してフェージング等により低下した受信電力の精度を補償するものである。

また、この方法は、信号の形式が既知である場合、実際に受信した受信信号と、フェージング等の影響を受けた場合の既知の受信信号との信号形式が一致したことに基づいて、フェージング等により低下した受信電力の精度を補償するものである。

Yunhao Liu and Mo Li, "Iso-Map: Energy Efficient Contour Mapping in Wireless Sensor Networks," Proc. In ICDCS2007, 2007. Tecfik Yucek and Huseyin Arslan, "A Survey of Spectrum Sensing Algorithms for Cognitive Radio Applications," IEEE Communications Surveys & Tutorials, vol. 11, no.1, 2009.

しかし、非特許文献１，２に記載された受信電力の等高線の作成方法は、受信電力の２次元的な等高線であるため、受信電力の３次元的な電力分布を推定することは困難である。最近では、ドローンが様々な用途に用いられており、ドローンが発する電波が１次利用者と干渉しないようにするために受信電力の３次元的な電力分布を推定する必要がある。

そこで、この発明の実施の形態によれば、受信電力の３次元的な電力分布を推定可能な電力推定装置を提供する。

また、この発明の実施の形態によれば、受信電力の３次元的な電力分布を推定可能な電力推定装置を備える無線通信システムを提供する。

更に、この発明の実施の形態によれば、受信電力の３次元的な電力分布の推定をコンピュータに実行させるためのプログラムを提供する。

更に、この発明の実施の形態によれば、受信電力の３次元的な電力分布を推定可能なデータ構造を提供する。

（構成１）
この発明の実施の形態によれば、電力推定装置は、電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線を作成する電力推定装置であって、抽出手段と、処理手段と、演算手段と、作成手段とを備える。抽出手段は、複数の端末装置から送信され、かつ、各々が端末装置の位置を示す位置情報と端末装置における電波の受信電力とを含む複数のモニター情報から、受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する。処理手段は、抽出手段によって抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点から最も遠い位置に存在する端末装置と受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する。演算手段は、処理手段によって求められたｍ個の受信電力の重心点およびｍ個の最大距離と、ｍ個の高さとに基づいて受信電力の重心点と高さとの関係を示す第１の関係式と、受信電力の複数の電力値に対応付けられ、かつ、各々が最大距離と高さとの関係を示す複数の第２の関係式とを演算する。作成手段は、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、複数の第２の関係式のうち、作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線を作成する。

構成１によれば、任意の高さにおける受信電力の重心点は、第１の関係式に基づいて求められ、作成すべき受信電力の等高線の電力値に対応し、かつ、任意の高さにおける最大距離は、第２の関係式に基づいて求められるので、任意の高さにおいて、受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線の電力値を有する等高線が作成される。この等高線の作成を複数の高さについて行うことによって３次元的な電力分布を推定できる。

（構成２）
構成１において、作成手段は、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求めるとともに第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点を中心とし、かつ、求めた最大距離を半径とする円形形状を有する第１の等高線、求めた受信電力の重心点を中心とし、かつ、求めた最大距離を長径とする楕円形状を有する第２の等高線、および求めた受信電力の重心点を中心とし、かつ、求めた最大距離を中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを、作成すべき受信電力の等高線として作成する。

構成２によれば、各種の平面形状を有する受信電力の等高線を作成できる。

（構成３）
構成１または構成２において、電力推定装置は、推定手段を更に備える。推定手段は、作成手段によって作成された受信電力の等高線に基づいて対象領域内の所望の場所において電波の受信電力を推定する。

構成３によれば、対象領域内の所望の場所における電波の受信電力を推定できる。

（構成４）
構成３において、推定手段は、対象領域における電波環境に応じた電波伝搬モデルを用いて、対象領域内の所望の場所において電波の受信電力を推定する。

構成４によれば、建物が密集した電波環境および建物が少ない電波環境等の各種の電波環境に適して、対象領域内の所望の場所において電波の受信電力を推定できる。

（構成５）
構成１から構成４のいずれかにおいて、処理手段は、ｍ個の高さに対応して対象領域内の高さ方向をｍ個の領域に分割し、その分割した１つの領域に含まれる高さを有するｋ個のモニター情報に基づいて第１の処理をｍ個の領域の全てについて実行する。

構成５によれば、高さ方向のｍ個の領域の各々について、受信電力の重心点および最大距離を求めるので、その求めた受信電力の重心点および最大距離を用いて第１および第２の関係式を容易に算出できる。

（構成６）
構成５において、処理手段は、対象領域内の建物の高さまたは建物の平均高さよりも小さい高さ方向の長さを有するように対象領域内の高さ方向をｍ個の領域に分割する。

構成６によれば、高さ方向に精度良く受信電力の重心点および最大距離を求めることができるので、第１および第２の関係式を精度良く算出できる。

（構成７）
構成５または構成６において、ｍ個の領域は、第１の領域と２の領域とを含む。第１の領域は、建物の高さまたは建物の平均高さ以下の位置に存在し、かつ、高さ方向に第１の長さを有する。第２の領域は、建物の高さまたは建物の平均高さよりも高い位置に存在し、かつ、高さ方向に第１の長さよりも長い第２の長さを有する第２の領域とを含む。

構成７によれば、建物の高さまたは建物の平均高さ以下の位置における実際の電波環境と、建物の高さまたは建物の平均高さよりも高い位置における実際の電波環境とに応じて、第１および第２の関係式を精度良く算出できる。

（構成８）
また、この発明の実施の形態によれば、無線通信システムは、構成１から構成７のいずれかに記載の電力推定装置を備える。

構成８によれば、電力推定装置によって推定された３次元的な電波の電力分布に基づいて無線通信を行う無線通信システムを実現できる。

（構成９）
更に、この発明の実施の形態によれば、プログラムは、電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線の作成をコンピュータに実行させるためのプログラムであって、抽出手段が、複数の端末装置から送信され、かつ、各々が端末装置の位置を示す位置情報と端末装置における電波の受信電力とを含む複数のモニター情報から、受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する第１のステップと、処理手段が、第１のステップにおいて抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点から最も遠い位置に存在する端末装置と受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する第２のステップと、演算手段が、第２のステップにおいて求められたｍ個の受信電力の重心点およびｍ個の最大距離と、ｍ個の高さとに基づいて受信電力の重心点と高さとの関係を示す第１の関係式と、受信電力の複数の電力値に対応付けられ、かつ、各々が最大距離と前記高さとの関係を示す複数の第２の関係式とを演算する第３のステップと、作成手段が、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、複数の第２の関係式のうち、作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線を作成する第４のステップとをコンピュータに実行させるためのプログラムである。

構成９によれば、プログラムを実行することにより、任意の高さにおける受信電力の重心点は、第１の関係式に基づいて求められ、作成すべき受信電力の等高線の電力値に対応し、かつ、任意の高さにおける最大距離は、第２の関係式に基づいて求められるので、任意の高さにおいて、受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線の電力値を有する等高線が作成される。この等高線の作成を複数の高さについて行うことによって３次元的な電力分布を推定できる。

（構成１０）
構成９において、作成手段は、第４のステップにおいて、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求めるとともに第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点を中心とし、かつ、その求めた最大距離を半径とする円形形状を有する第１の等高線、求めた受信電力の重心点を中心とし、かつ、求めた最大距離を長径とする楕円形状を有する第２の等高線、および求めた受信電力の重心点を中心とし、かつ、求めた最大距離を中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを、作成すべき受信電力の等高線として作成する。

構成１０によれば、プログラムを実行することにより、各種の平面形状を有する受信電力の等高線を作成できる。

（構成１１）
構成９または構成１０において、プログラムは、推定手段が、第４のステップにおいて作成された受信電力の等高線に基づいて対象領域内の所望の場所において電波の受信電力を推定する第５のステップを更にコンピュータに実行させる。

構成１１によれば、プログラムを実行することにより、対象領域内の所望の場所における電波の受信電力を推定できる。

（構成１２）
構成１１において、推定手段は、第５のステップにおいて、対象領域における電波環境に応じた電波伝搬モデルを用いて、対象領域内の所望の場所において電波の受信電力を推定する。

構成１２によれば、プログラムを実行することにより、建物が密集した電波環境および建物が少ない電波環境等の各種の電波環境に適して、対象領域内の所望の場所において電波の受信電力を推定できる。

（構成１３）
構成９から構成１２のいずれかにおいて、処理手段は、第２のステップにおいて、ｍ個の高さに対応して対象領域内の高さ方向をｍ個の領域に分割し、その分割した１つの領域に含まれる高さを有するｋ個のモニター情報に基づいて第１の処理をｍ個の領域の全てについて実行する。

構成１３によれば、プログラムを実行することにより、高さ方向のｍ個の領域の各々について、受信電力の重心点および最大距離を求めるので、その求めた受信電力の重心点および最大距離を用いて第１および第２の関係式を容易に算出できる。

（構成１４）
構成１３において、処理手段は、第２のステップにおいて、対象領域内の建物の高さまたは建物の平均高さよりも小さい高さ方向の長さを有するように対象領域内の高さ方向をｍ個の領域に分割する。

構成１４によれば、プログラムを実行することにより、高さ方向に精度良く受信電力の重心点および最大距離を求めることができるので、第１および第２の関係式を精度良く算出できる。

（構成１５）
構成１３または構成１４において、ｍ個の領域は、第１の領域と第２の領域とを含む。第１の領域は、建物の高さまたは建物の平均高さ以下の位置に存在し、かつ、高さ方向に第１の長さを有する。第２の領域は、建物の高さまたは建物の平均高さよりも高い位置に存在し、かつ、高さ方向に第１の長さよりも長い第２の長さを有する。

構成１５によれば、プログラムを実行することにより、建物の高さまたは建物の平均高さ以下の位置における実際の電波環境と、建物の高さまたは建物の平均高さよりも高い位置における実際の電波環境とに応じて、第１および第２の関係式を精度良く算出できる。

（構成１６）
更に、この発明の実施の形態によれば、データ構造は、複数のモニター情報を含むデータ構造であって、複数のモニター情報の各々は、端末装置の位置を示す位置情報と、端末装置が無線通信に共用する共用周波数と、端末装置によって検出された電波の受信電力とを相互に対応付けた構造を有する。複数のモニター情報に含まれる複数の受信電力のうち、受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報は、抽出手段によって複数のモニター情報から抽出される。抽出手段によって抽出されたｋ個のモニター情報は、処理手段が受信電力によって重み付けされたｋ個の端末装置の位置の平均である受信電力の重心点を求めるのに用いられる。抽出手段によって抽出されたｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力と受信電力の等高線における電力値以上の受信電力に対応付けられた位置情報とは、処理手段が、受信電力の重心点から最も遠い位置に存在する端末装置と受信電力の重心点との距離である最大距離を求めるのに用いられる。処理手段によって求められた受信電力の重心点は、演算手段が受信電力の重心点と高さとの関係を示す第１の関係式を求めるのに用いられる。、処理手段によって求められた最大距離は、演算手段が、作成すべき受信電力の等高線の電力値に対応し、かつ、最大距離と高さとの関係を示す第２の関係式を求めるのに用いられる。演算手段によって求められた第１および第２の関係式は、作成手段が、任意の高さにおいて、任意の電力値を有する受信電力の等高線を作成するのに用いられる。

構成１６によれば、この発明の実施の形態によるデータ構造を用いることによって、受信電力の３次元的な電力分布を推定できる。

受信電力の３次元的な電力分布を推定できる。

この発明の実施の形態による無線通信システムを示す概略図である。 図１に示す電力推定装置の概略図である。 対象領域ＲＥＧを分割した分割領域の概念図である。 図３に示す１つの分割領域ＲＥＧ＿ＤＶｑを高さ方向に分割する方法を説明するための図である。 受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離Ｌ_ＭＡＸを求める方法を説明するための図である。 受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。 受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。 最大距離Ｌ_ＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。 任意の高さにおける受信電力の等高線を求める方法を説明するための図である。 等高線の概念図である。 受信電力と波源からの距離との関係を示す図である。 受信電力の等高線に基づいて任意の位置における受信電力の電力値を推定する方法を説明するための図である。 奥村−秦モデルにおける電界強度と距離との関係を示す図である。 端末装置における電波の受信感度以下の領域において受信電力を推定する方法を説明するための図である。 図１に示す電力推定装置の動作を説明するためのフローチャートである。 図１５に示すフローチャートのステップの実行状態を示す概念図である。 図１に示す電力推定装置の動作を説明するための別のフローチャートである。 図１７のステップＳ１６の詳細な動作を説明するためのフローチャートである。 図１７および図１８に示すフローチャートのステップの実行状態を示す概念図である。 図１７および図１８に示すフローチャートの実行が完了した状態を示す概念図である。 図１に示す電力推定装置の別の動作を説明するためのフローチャートである。 楕円形状を有する受信電力の等高線の作成方法を説明するための図である。 多角形の形状を有する受信電力の等高線の作成方法を説明するための図である。 この発明の実施の形態によるデータ構造の概略図である。 この発明の実施の形態による別のデータ構造の概略図である。

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

図１は、この発明の実施の形態による無線通信システムを示す概略図である。図１を参照して、この発明の実施の形態による無線通信システム１０は、電力推定装置１と、複数の端末装置２と、基地局３とを備える。

電力推定装置１は、有線ケーブル４によって基地局３に接続される。電力推定装置１は、基地局３から有線ケーブル４を介して複数の端末装置２における複数のモニター情報を受信し、その受信した複数のモニター情報を記録する。ここで、各モニター情報は、端末装置２の位置を示す位置情報、波源Ｓから送信された電波の端末装置２における受信電力、受信電力を検出したときの時刻を示す時刻情報、および波源Ｓから送信された電波の周波数を示す周波数情報を含む。そして、位置情報は、例えば、ｘ座標、ｙ座標およびｚ座標によって表される。ｚ座標は、端末装置２が配置される領域の高さ方向の座標である。

また、電力推定装置１は、図１に図示されていない基地局から基地局３へ送信された複数のモニター情報も基地局３から有線ケーブル４を介して受信し、その受信した複数のモニター情報を記録する。

即ち、電力推定装置１は、市町村の領域または県の領域（電力を推定する対象領域ＲＥＧ）に存在する端末装置２のモニター情報を基地局３から受信して記録する。

そして、電力推定装置１は、その記録した複数のモニター情報に基づいて、後述する方法によって、任意の高さにおける受信電力の等高線を作成し、その作成した受信電力の等高線に基づいて、対象領域ＲＥＧの任意の位置における電波の受信電力を推定する。

このように、電力推定装置１は、非常に広い対象領域ＲＥＧにおける３次元的な電波の受信電力を推定する。

複数の端末装置２および基地局３は、無線通信空間に配置される。複数の端末装置２は、波源Ｓの周囲に配置される。

複数の端末装置２の各々は、３ＧＨｚ以下の周波数帯域、３ＧＨｚ〜６ＧＨｚの周波数帯域、６ＧＨｚ〜３０ＧＨｚの周波数帯域および３０ＧＨｚよりも高周波数の周波数帯域のいずれかに含まれる複数の周波数のうち、複数の端末装置２が無線通信に共用する共用周波数ｆ_ｃｏｍを用いて他の端末装置２と無線通信を行う。

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

また、複数の端末装置２の各々は、共用周波数ｆ_ｃｏｍを有する電波を波源Ｓから受信し、電波を受信したときの受信電力ＲＳＳＩ_ｉ（ｉ＝１，２，３，・・・）を検出する。

更に、複数の端末装置２の各々は、タイマーを内蔵しており、受信電力ＲＳＳＩ_ｉを検出したときの時刻を検出する。

そして、複数の端末装置２の各々は、自己の位置を示す位置情報（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ）と、受信電力ＲＳＳＩ_ｉと、受信電力ＲＳＳＩ_ｉを検出したときの時刻を示す時間情報ｔ_ｉと、共用周波数ｆ_ｃｏｍとを含むモニター情報ＭＮＴ＿ｉ＝［ｔ_ｉ：（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ），ＲＳＳＩ_ｉ，ｆ_{ｃｏｍ＿ｉ}］を生成し、その生成したモニター情報ＭＮＴ＿ｉ＝［ｔ_ｉ：（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ），ＲＳＳＩ_ｉ，ｆ_{ｃｏｍ＿ｉ}］を無線通信によって基地局３へ送信する。

この場合、端末装置２は、モニター情報ＭＮＴを基地局３へ直接送信してもよく、モニター情報ＭＮＴをマルチホップによって基地局３へ送信してもよい。

複数の端末装置２の各々は、他の端末装置２と無線通信を行っているときを除いて、位置情報（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ）、受信電力ＲＳＳＩ_ｉ、時刻情報ｔ_ｉおよび共用周波数ｆ_{ｃｏｍ＿ｉ}を検出してモニター情報ＭＮＴ＿ｉ＝［ｔ_ｉ：（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ），ＲＳＳＩ_ｉ，ｆ_{ｃｏｍ＿ｉ}］を生成し、その生成したモニター情報ＭＮＴ＿ｉ＝［ｔ_ｉ：（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ），ＲＳＳＩ_ｉ，ｆ_{ｃｏｍ＿ｉ}］を基地局３へ送信する動作を定期的（例えば、１分間隔）に行う。

基地局３は、モニター情報ＭＮＴを複数の端末装置２から受信し、その受信した複数のモニター情報ＭＮＴを有線ケーブル４を介して電力推定装置１へ送信する。

また、基地局３は、図示されていない基地局から受信した複数のモニター情報ＭＮＴを有線ケーブル４を介して電力推定装置１へ送信する。

波源Ｓは、例えば、レーダ、テレビジョン放送用の無線中継伝送装置（ＦＰＵ：Field Pickup Unit）、離島へ電波を中継する中継所および端末装置等からなる。

図２は、図１に示す電力推定装置１の概略図である。図２を参照して、電力推定装置１は、受信手段１１と、記録手段１２と、抽出手段１３と、処理手段１４と、演算手段１５と、作成手段１６と、推定手段１７とを備える。

受信手段１１は、複数のモニター情報ＭＮＴを基地局３から受信し、その受信した複数のモニター情報ＭＮＴを記録手段１２に記録する。

記録手段１２は、受信手段１１から受けた複数のモニター情報ＭＮＴを記録する。より具体的には、記録手段１２は、複数のモニター情報ＭＮＴを［ｔ_１：（ｘ_１，ｙ_１，ｚ_１），ＲＳＳＩ_１，ｆ_{ｃｏｍ＿１}］，［ｔ_２：（ｘ_２，ｙ_２，ｚ_２），ＲＳＳＩ_２，ｆ_{ｃｏｍ＿１}］，［ｔ_３：（ｘ_３，ｙ_３，ｚ_３），ＲＳＳＩ_３，ｆ_{ｃｏｍ＿１}］，・・・,［ｔ_Ｔ：（ｘ_Ｍ，ｙ_Ｍ，ｚ_Ｍ），ＲＳＳＩ_Ｍ，ｆ_{ｃｏｍ＿１}］；［ｔ_１：（ｘ_１，ｙ_１，ｚ_１），ＲＳＳＩ_１，ｆ_{ｃｏｍ＿２}］，［ｔ_２：（ｘ_２，ｙ_２，ｚ_２），ＲＳＳＩ_２，ｆ_{ｃｏｍ＿２}］，［ｔ_３：（ｘ_３，ｙ_３，ｚ_３），ＲＳＳＩ_３，ｆ_{ｃｏｍ＿２}］，・・・,［ｔ_Ｔ：（ｘ_Ｍ，ｙ_Ｍ，ｚ_Ｍ），ＲＳＳＩ_Ｍ，ｆ_{ｃｏｍ＿２}］；・・・；［ｔ_１：（ｘ_１，ｙ_１，ｚ_１），ＲＳＳＩ_１，ｆ_{ｃｏｍ＿Ｆ}］，［ｔ_２：（ｘ_２，ｙ_２，ｚ_２），ＲＳＳＩ_２，ｆ_{ｃｏｍ＿Ｆ}］，［ｔ_３：（ｘ_３，ｙ_３，ｚ_３），ＲＳＳＩ_３，ｆ_{ｃｏｍ＿Ｆ}］，・・・,［ｔ_Ｔ：（ｘ_Ｍ，ｙ_Ｍ，ｚ_Ｍ），ＲＳＳＩ_Ｍ，ｆ_{ｃｏｍ＿Ｆ}］のように記録する。ここで、Ｔは、時刻情報ｔ_ｉの検出回数を示し、Ｍは、端末装置２の個数を示し、Ｆは、共用周波数ｆ_ｃｏｍの個数を示す。そして、Ｔは、２以上の整数であり、Ｍは、後述するｋよりも大きい整数であり、Ｆは、１以上の整数である。

また、記録手段１２は、移動端末であるが停止している端末装置２または静止端末である端末装置２に対しては、位置情報（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ）が同じで、時刻情報ｔ_ｉ、受信電力ＲＳＳＩ_ｉおよび周波数情報ｆ_{ｃｏｍ＿ｊ}（ｊは、１≦ｊ≦Ｆを満たす整数）が異なるようにモニター情報ＭＮＴを記録する。

更に、記録手段１２は、時刻情報ｔ_ｉの検出間隔が位置情報（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ）、受信電力ＲＳＳＩ_ｉおよび周波数情報ｆ_{ｃｏｍ＿ｊ}の検出間隔よりも長い場合、複数のモニター情報ＭＮＴを［ｔ_１：（ｘ_１，ｙ_１，ｚ_１），ＲＳＳＩ_１，ｆ_{ｃｏｍ＿１}］，［ｔ_１：（ｘ_２，ｙ_２，ｚ_２），ＲＳＳＩ_２，ｆ_{ｃｏｍ＿１}］，［ｔ_１：（ｘ_３，ｙ_３，ｚ_３），ＲＳＳＩ_３，ｆ_{ｃｏｍ＿１}］，・・・,［ｔ_１：（ｘ_Ｍ，ｙ_Ｍ，ｚ_Ｍ），ＲＳＳＩ_Ｍ，ｆ_{ｃｏｍ＿１}］；［ｔ_２：（ｘ_１，ｙ_１，ｚ_１），ＲＳＳＩ_１，ｆ_{ｃｏｍ＿２}］，［ｔ_２：（ｘ_２，ｙ_２，ｚ_２），ＲＳＳＩ_２，ｆ_{ｃｏｍ＿２}］，［ｔ_２：（ｘ_３，ｙ_３，ｚ_３），ＲＳＳＩ_３，ｆ_{ｃｏｍ＿２}］，・・・,［ｔ_２：（ｘ_Ｍ，ｙ_Ｍ，ｚ_Ｍ），ＲＳＳＩ_Ｍ，ｆ_{ｃｏｍ＿２}］；・・・；［ｔ_Ｔ：（ｘ_１，ｙ_１，ｚ_１），ＲＳＳＩ_１，ｆ_{ｃｏｍ＿Ｆ}］，［ｔ_Ｔ：（ｘ_２，ｙ_２，ｚ_２），ＲＳＳＩ_２，ｆ_{ｃｏｍ＿Ｆ}］，［ｔ_Ｔ：（ｘ_３，ｙ_３，ｚ_３），ＲＳＳＩ_３，ｆ_{ｃｏｍ＿Ｆ}］，・・・,［ｔ_Ｔ：（ｘ_Ｍ，ｙ_Ｍ，ｚ_Ｍ），ＲＳＳＩ_Ｍ，ｆ_{ｃｏｍ＿Ｆ}］のように記録する。即ち、記録手段１２は、同じ時刻情報ｔ_ｉにおいて、位置情報（ｘ_ｉ，ｙ_ｉ，ｚ_ｉ）および受信電力ＲＳＳＩ_ｉが変化するように複数のモニター情報ＭＮＴを記録する。

抽出手段１３は、後述する方法によって、記録手段１２に記録された複数のモニター情報ＭＮＴに基づいて、複数のモニター情報ＭＮＴの分布状態を作成し、その作成した複数のモニター情報ＭＮＴの分布状態に基づいて受信電力の等高線を作成する処理を行うための処理範囲を決定し、その決定した処理範囲に含まれるｋ（ｋは２以上の整数）個の端末装置２のｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを抽出する。そして、抽出手段１３は、その抽出したｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを処理手段１４へ出力する。抽出手段１３は、作成すべき受信電力の等高線ＣＴＲの電力値Ｐ以上の受信電力ＲＳＳＩ_ｉがなくなるまで、後述する方法によって、複数のモニター情報ＭＮＴからｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを抽出し、その抽出したｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを処理手段１４へ出力する処理を繰り返し行う。

処理手段１４は、抽出手段１３からｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを受け、その受けたｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに基づいて、後述する方法によって、受信電力の重心点ＲＳＳＩ＿Ｇを算出するとともに、ｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇから最も遠い位置に存在する端末装置２と受信電力の重心点ＲＳＳＩ＿Ｇとの距離を最大距離Ｌ_ＭＡＸとして求める。そして、処理手段１４は、受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離Ｌ_ＭＡＸを求める処理をｍ（ｍは、２以上の整数）個の高さｚ_１〜ｚ_ｍについて実行する。そうすると、処理手段１４は、ｍ個の受信電力の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離Ｌ_{ＭＡＸ＿１}〜Ｌ_{ＭＡＸ＿ｍ}およびｍ個の高さｚ_１〜ｚ_ｍを演算手段１５へ出力する。

演算手段１５は、ｍ個の受信電力の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離Ｌ_{ＭＡＸ＿１}〜Ｌ_{ＭＡＸ＿ｍ}およびｍ個の高さｚ_１〜ｚ_ｍを処理手段１４から受ける。そして、演算手段１５は、ｍ個の受信電力の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離Ｌ_{ＭＡＸ＿１}〜Ｌ_{ＭＡＸ＿ｍ}およびｍ個の高さｚ_１〜ｚ_ｍに基づいて、後述する方法によって、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇと高さｚとの関係を示す関係式ＥＱ１、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙ_Ｇと高さｚとの関係を示す関係式ＥＱ２および最大距離Ｌ_ＭＡＸと高さｚとの関係を示す関係式ＥＱ３を演算する。そうすると、演算手段１５は、関係式ＥＱ１，ＥＱ２，ＥＱ３を作成手段１６へ出力する。

作成手段１６は、関係式ＥＱ１，ＥＱ２，ＥＱ３を演算手段１５から受ける。そして、作成手段１６は、関係式ＥＱ１，ＥＱ２，ＥＱ３に基づいて、後述する方法によって、受信電力の３次元的な電力分布を示す等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲを推定手段１７へ出力する。

推定手段１７は、受信電力の等高線ＣＴＲを作成手段１６から受け、その受けた受信電力の等高線ＣＴＲに基づいて、後述する方法によって、対象領域ＲＥＧの任意の位置における電波の受信電力を推定する。そして、推定手段１７は、その推定した電波の受信電力に基づいて、複数の端末装置２が共用周波数ｆ_ｃｏｍを用いて無線通信を行うときの通信条件（共用周波数ｆ_ｃｏｍ、場所および送信電力等）を決定する。

図３は、対象領域ＲＥＧを分割した分割領域の概念図である。図３を参照して、対象領域ＲＥＧは、水平方向において分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑ（ｑは２以上の整数）に分割される。

分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々は、地面に平行な水平方向において長さＤを有し、地面に垂直な方向において高さＨを有する。例えば、長さＤは、ドローンの飛行領域における水平方向の長さに設定され、高さＨは、ドローンの飛行領域における垂直方向の長さに設定される。

なお、分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑは、相互に、同じ長さＤおよび高さＨを有していてもよいし、相互に、異なる長さＤおよび高さＨを有していてもよい。

図４は、図３に示す１つの分割領域ＲＥＧ＿ＤＶｑを高さ方向に分割する方法を説明するための図である。図４を参照して、抽出手段１３は、建物２０の高さＨ_Ｂよりも小さい高さ方向の長さｈを予め保持している。長さｈは、例えば、３〜１０ｍである。そして、抽出手段１３は、対象領域ＲＥＧの高さ方向を長さｈで分割して複数の領域ＳＰＣ１〜ＳＰＣｍを生成する。なお、高さＨ_Ｂは、建物の平均高さであってもよい。

抽出手段１３は、複数の領域ＳＰＣ１〜ＳＰＣｍを生成すると、モニター情報ＭＮＴに含まれる位置情報のｚ座標に基づいて複数の領域ＳＰＣ１〜ＳＰＣｍの各々に含まれるモニター情報を抽出する。この場合、抽出手段１３は、領域ＳＰＣ１に含まれるモニター情報ＭＮＴ＿ＳＰＣ１、領域ＳＰＣ２に含まれるモニター情報ＭＮＴ＿ＳＰＣ２、・・・、および領域ＳＰＣｍに含まれるモニター情報ＭＮＴ＿ＳＰＣｍを抽出する。モニター情報ＭＮＴ＿ＳＰＣ１〜ＭＮＴ＿ＳＰＣｍの各々は、複数のモニター情報からなる。

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

図５の（ａ）を参照して、電力推定装置１の抽出手段１３は、１つの領域ＳＰＣ１に含まれる複数のモニター情報ＭＮＴに基づいて、複数のモニター情報ＭＮＴの分布状態（図５の（ａ）に示す分布状態）を取得する。

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

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

抽出手段１３は、作成すべき受信電力の等高線ＣＴＲの電力値Ｐを決定すると、複数のモニター情報ＭＮＴに含まれる複数の受信電力ＲＳＳＩ_ｉのうち、電力値Ｐ以上である受信電力の最大値ＲＳＳＩ＿ＭＡＸを抽出する。ここで、受信電力の最大値ＲＳＳＩ＿ＭＡＸを抽出するのは、受信電力の最大値ＲＳＳＩ＿ＭＡＸを有する端末装置２が最も波源Ｓに近いと考えられるからである。

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

そうすると、抽出手段１３は、処理範囲ＲＥＧ＿ＰＲＳ内のｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを抽出し、その抽出したｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを処理手段１４へ出力する。

図５の（ｂ）を参照して、処理手段１４は、ｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを抽出手段１３から受け、その受けたｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに含まれるｋ個の受信電力ＲＳＳＩ_１〜ＲＳＳＩ_ｋに基づいて、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）を次式によって算出する。

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

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

次に、処理手段１４は、ｋ個の受信電力ＲＳＳＩ_ｉ（ｉ＝１〜ｋ）のうち、電力値Ｐ以上の受信電力ＲＳＳＩ_ｉを有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）から最も遠い位置に存在する端末装置２−１と受信電力の重心点ＲＳＳＩ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）との間の距離を最大距離Ｌ_ＭＡＸとして求める。

このように、処理手段１４は、電力値Ｐ未満の受信電力（黒四角）を除外し、電力値Ｐ以上の受信電力（白四角）のみに基づいて最大距離Ｌ_ＭＡＸを求める。黒四角によって示される測定点のうち、測定点ｓ１〜ｓ３は、フェージングの影響によって受信電力が電力値Ｐ未満になった測定点である。

なお、処理手段１４は、上述した方法によって、複数の電力値Ｐ（＝Ｐ_１〜Ｐ_Ｎ（Ｎは、２以上の整数））について最大距離Ｌ_{ＭＡＸ＿１}〜Ｌ_{ＭＡＸ＿Ｎ}を求める。なお、複数の最大距離Ｌ_{ＭＡＸ＿１}〜Ｌ_{ＭＡＸ＿Ｎ}は、それぞれ、作成すべき受信電力の等高線ＣＴＲの電力値Ｐ_１〜Ｐ_Ｎに対応する最大距離である。

従って、処理手段１４による受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離Ｌ_ＭＡＸを求める方法によれば、フェージングの影響を抑制して受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離Ｌ_ＭＡＸを求めることができる。

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

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

そうすると、受信電力の最大値ＲＳＳＩ＿ＭＡＸは、電力値Ｐ以上であるので、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘ_Ｇ，ｙ_Ｇ）と受信電力の最大値ＲＳＳＩ＿ＭＡＸを有する測定点との距離を最大距離Ｌ_ＭＡＸとして求めることができる。従って、ｋは、２以上であればよい。

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

抽出手段１３は、領域ＳＰＣ２〜ＳＰＣｍの各々においても、上述した方法によってｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを抽出し、処理手段１４は、抽出手段１３によって抽出されたｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに基づいて、上述した方法によって受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離Ｌ_ＭＡＸ（＝Ｌ_{ＭＡＸ＿１}〜Ｌ_{ＭＡＸ＿Ｎ}）を求める。

抽出手段１３および処理手段１４は、上述した処理を複数の領域ＳＰＣ１〜ＳＰＣｍの全てについて実行し、複数の領域ＳＰＣ１〜ＳＰＣｍの各々について受信電力の重心点ＲＳＳＩ＿Ｇおよび複数の最大距離Ｌ_{ＭＡＸ＿１}〜Ｌ_{ＭＡＸ＿Ｎ}を求める。その結果、処理手段１４は、領域ＳＰＣ１において、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ１および複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ１}〜Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ１}を求め、領域ＳＰＣ２において、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ２および複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ２}〜Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ２}を求め、・・・、領域ＳＰＣｍにおいて、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣｍおよび複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣｍ}〜Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣｍ}を求める。

そうすると、処理手段１４は、領域ＳＰＣ１、領域ＳＰＣ１におけるｋ個のモニター情報ＭＮＴ_{１＿ＳＰＣ１}〜ＭＮＴ_{ｋ＿ＳＰＣ１}、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ１および複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ１}〜Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ１}を相互に対応付けた情報ＩＮＦ１＝［ＳＰＣ１／モニター情報ＭＮＴ_{１＿ＳＰＣ１}〜ＭＮＴ_{ｋ＿ＳＰＣ１}／ＲＳＳＩ＿Ｇ＿ＳＰＣ１／Ｌ_{ＭＡＸ＿１＿ＳＰＣ１}〜Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ１}］を生成する。また、処理手段１４は、同様にして、領域ＳＰＣ２に対応する情報ＩＮＦ２＝［ＳＰＣ２／モニター情報ＭＮＴ_{１＿ＳＰＣ２}〜ＭＮＴ_{ｋ＿ＳＰＣ２}／ＲＳＳＩ＿Ｇ＿ＳＰＣ２／Ｌ_{ＭＡＸ＿１＿ＳＰＣ２}〜Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ２}］、・・・、領域ＳＰＣｍに対応する情報ＩＮＦｍ＝［ＳＰＣｍ／モニター情報ＭＮＴ_{１＿ＳＰＣｍ}〜ＭＮＴ_{ｋ＿ＳＰＣｍ}／ＲＳＳＩ＿Ｇ＿ＳＰＣｍ／Ｌ_{ＭＡＸ＿１＿ＳＰＣｍ}〜Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣｍ}］を生成する。そして、処理手段１４は、情報ＩＮＦ１〜ＩＮＦ４を演算手段１５へ出力する。

図６は、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。図６においては、高さ方向の領域ＳＰＣが４個（ｍ＝４）である場合について、重心点ＲＳＳＩ＿Ｇのｘ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を説明する。

図６において、縦軸は、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標を表し、横軸は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚを表す。各領域ＳＰＣ１〜ＳＰＣ４の高さは、例えば、各領域ＳＰＣ１〜ＳＰＣ４において抽出されたｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに含まれるｋ個の位置情報のｋ個のｚ座標の平均値として求められる。また、各領域ＳＰＣ１〜ＳＰＣ４の高さは、各領域ＳＰＣ１〜ＳＰＣ４において抽出されたｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに含まれるｋ個の位置情報のｋ個のｚ座標をｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに含まれるｋ個の受信電力ＲＳＳＩ_１〜ＲＳＳＩ_ｋによって重み付けした重み付け平均（ｚ_Ｇ＝（ｗ_１ｚ_１＋ｗ_２ｚ_２＋・・・＋ｗ_ｋｚ_ｋ）／（ｗ_１＋ｗ_２＋・・・＋ｗ_ｋ）、ｗ_ｉ＝ＲＳＳＩ_ｉ）によって求められる。

演算手段１５は、情報ＩＮＦ１〜ＩＮＦ４を処理手段１４から受ける。演算手段１５は、情報ＩＮＦ１＝［ＳＰＣ１／モニター情報ＭＮＴ_{１＿ＳＰＣ１}〜ＭＮＴ_{ｋ＿ＳＰＣ１}／ＲＳＳＩ＿Ｇ＿ＳＰＣ１／Ｌ_{ＭＡＸ＿１＿ＳＰＣ１}〜Ｌ_{ＭＡＸ＿ｎ＿ＳＰＣ１}］からｋ個のモニター情報ＭＮＴ_{１＿ＳＰＣ１}〜ＭＮＴ_{ｋ＿ＳＰＣ１}を検出し、その検出したｋ個のモニター情報ＭＮＴ_{１＿ＳＰＣ１}〜ＭＮＴ_{ｋ＿ＳＰＣ１}に含まれるｋ個の位置情報のｋ個のｚ座標ｚ_１〜ｚ_ｋ（またはｋ個のｚ座標ｚ_１〜ｚ_ｋおよびｋ個の受信電力ＲＳＳＩ_１〜ＲＳＳＩ_ｋ）を検出する。そして、演算手段１５は、ｋ個のｚ座標ｚ_１〜ｚ_ｋ（またはｋ個のｚ座標ｚ_１〜ｚ_ｋおよびｋ個の受信電力ＲＳＳＩ_１〜ＲＳＳＩ_ｋ）に基づいて、上述した方法によって領域ＳＰＣ１の高さｚ_ＳＰＣ１を演算する。同様にして、演算手段１５は、情報ＩＮＦ２〜ＩＮＦ４に基づいて、それぞれ、領域ＳＰＣ２〜ＳＰＣ４の高さｚ_ＳＰＣ２〜ｚ_ＳＰＣ４を演算する。

そうすると、演算手段１５は、情報ＩＮＦ１〜ＩＮＦ４からそれぞれ受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ１〜ＲＳＳＩ＿Ｇ＿ＳＰＣ４のｘ座標ｘ_{Ｇ＿ＳＰＣ１}〜ｘ_{Ｇ＿ＳＰＣ４}およびｙ座標ｙ_{Ｇ＿ＳＰＣ１}〜ｙ_{Ｇ＿ＳＰＣ４}を検出する。

そして、演算手段１５は、ｘ座標ｘ_{Ｇ＿ＳＰＣ１}〜ｘ_{Ｇ＿ＳＰＣ４}および高さｚ_ＳＰＣ１〜ｚ_ＳＰＣ４に基づいて、点（ｘ_{Ｇ＿ＳＰＣ１}，ｚ_ＳＰＣ１），点（ｘ_{Ｇ＿ＳＰＣ２}，ｚ_ＳＰＣ２），点（ｘ_{Ｇ＿ＳＰＣ３}，ｚ_ＳＰＣ３），点（ｘ_{Ｇ＿ＳＰＣ４}，ｚ_ＳＰＣ４）をプロットする（図６の黒丸参照）。

その後、演算手段１５は、プロットした点に基づいて受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_{Ｇ＿ＳＰＣ１}〜ｘ_{Ｇ＿ＳＰＣ４}と各領域ＳＰＣ１〜ＳＰＣ４の高さｚ_ＳＰＣ１〜ｚ_ＳＰＣ４との関係を示す関係式ＥＱ１（曲線ｋ１参照）を求める。従って、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ（ｚ）は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚの関数である。

図７は、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。

図７において、縦軸は、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標を表し、横軸は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚを表す。各領域ＳＰＣ１〜ＳＰＣ４の高さは、図６において説明した方法によって求められる。

演算手段１５は、ｙ座標ｙ_{Ｇ＿ＳＰＣ１}〜ｙ_{Ｇ＿ＳＰＣ４}および高さｚ_ＳＰＣ１〜ｚ_ＳＰＣ４に基づいて、点（ｙ_{Ｇ＿ＳＰＣ１}，ｚ_ＳＰＣ１），点（ｙ_{Ｇ＿ＳＰＣ２}，ｚ_ＳＰＣ２），点（ｙ_{Ｇ＿ＳＰＣ３}，ｚ_ＳＰＣ３），点（ｙ_{Ｇ＿ＳＰＣ４}，ｚ_ＳＰＣ４）をプロットする（図７の黒丸参照）。

その後、演算手段１５は、プロットした点に基づいて受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙ_{Ｇ＿ＳＰＣ１}〜ｙ_{Ｇ＿ＳＰＣ４}と各領域ＳＰＣ１〜ＳＰＣ４の高さｚ_ＳＰＣ１〜ｚ_ＳＰＣ４との関係を示す関係式ＥＱ２（曲線ｋ２参照）を求める。従って、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙ（ｚ）は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚの関数である。

図８は、最大距離Ｌ_ＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。図８において、縦軸は、最大距離Ｌ_ＭＡＸを表し、横軸は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚを表す。各領域ＳＰＣ１〜ＳＰＣ４の高さは、図６において説明した方法によって求められる。

演算手段１５は、情報ＩＮＦ１〜ＩＮＦ４から、それぞれ、複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ１}〜Ｌ_{ＭＡＸ＿ｎ＿ＳＰＣ１}、複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ２}〜Ｌ_{ＭＡＸ＿ｎ＿ＳＰＣ２}、複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ３}〜Ｌ_{ＭＡＸ＿ｎ＿ＳＰＣ３}、および複数の最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ４}〜Ｌ_{ＭＡＸ＿ｎ＿ＳＰＣ４}を検出する。

そして、演算手段１５は、同じ電力値Ｐ_１［ｄＢｍ］に対応する最大距離Ｌ_{ＭＡＸ＿１＿ＳＰＣ１}，Ｌ_{ＭＡＸ＿１＿ＳＰＣ２}，Ｌ_{ＭＡＸ＿１＿ＳＰＣ３}，Ｌ_{ＭＡＸ＿１＿ＳＰＣ４}を検出する。

そうすると、演算手段１５は、点（Ｌ_{ＭＡＸ＿１＿ＳＰＣ１}，ｚ_ＳＰＣ１）点（Ｌ_{ＭＡＸ＿１＿ＳＰＣ２}，ｚ_ＳＰＣ２），点（Ｌ_{ＭＡＸ＿１＿ＳＰＣ３}，ｚ_ＳＰＣ３），点（Ｌ_{ＭＡＸ＿１＿ＳＰＣ４}，ｚ_ＳＰＣ４）をプロットする（図８の黒丸参照）。

その後、演算手段１５は、プロットした点に基づいて最大距離Ｌ_ＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さｚとの関係を示す関係式ＥＱ３−１（直線ｋ３参照）を求める。

演算手段１５は、同様にして、点（Ｌ_{ＭＡＸ＿２＿ＳＰＣ１}，ｚ_ＳＰＣ１）点（Ｌ_{ＭＡＸ＿２＿ＳＰＣ２}，ｚ_ＳＰＣ２），点（Ｌ_{ＭＡＸ＿２＿ＳＰＣ３}，ｚ_ＳＰＣ３），点（Ｌ_{ＭＡＸ＿２＿ＳＰＣ４}，ｚ_ＳＰＣ４）をプロットして、電力値Ｐ_２について、最大距離Ｌ_ＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す関係式ＥＱ３−２を求める。

以下、同様にして、演算手段１５は、点（Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ１}，ｚ_ＳＰＣ１）点（Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ２}，ｚ_ＳＰＣ２），点（Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ３}，ｚ_ＳＰＣ３），点（Ｌ_{ＭＡＸ＿Ｎ＿ＳＰＣ４}，ｚ_ＳＰＣ４）をプロットして、電力値Ｐ_Ｎについて、最大距離Ｌ_ＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す関係式ＥＱ３−Ｎを求める。

従って、図８に示す縦軸の最大距離Ｌ_ＭＡＸ（Ｐ，ｚ）は、作成すべき等高線ＣＴＲの電力値Ｐと各領域ＳＰＣ１〜ＳＰＣ４の高さ（ｚ座標）との関数である。つまり、演算手段１５は、作成すべき等高線ＣＴＲの電力値Ｐ_１〜Ｐ_Ｎの個数Ｎに等しい個数Ｎの最大距離Ｌ_ＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す関係式ＥＱ３＿１〜ＥＱ３＿Ｎを求める。なお、関係式ＥＱ３＿１〜ＥＱ３＿Ｎは、最大距離Ｌ_ＭＡＸと高さｚとの関係を示す関係式ＥＱ３を構成する。

そして、演算手段１５は、関係式ＥＱ１、関係式ＥＱ２および関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎ）を作成手段１６へ出力する。

図９は、任意の高さにおける受信電力の等高線を求める方法を説明するための図である。

図９を参照して、任意の高さｚにおける受信電力の重心点ＲＳＳＩ＿Ｇを決定する関係式は、関係式ＥＱ１（＝ｘ（ｚ））と関係式ＥＱ２（＝ｙ（ｚ））とからなる。また、最大距離Ｌ_ＭＡＸを決定する関係式は、関係式ＥＱ３（＝Ｌ_ＭＡＸ（Ｐ，ｚ））からなる。

関係式ＥＱ３（＝Ｌ_ＭＡＸ（Ｐ，ｚ））は、作成すべき等高線ＣＴＲの電力値Ｐ_１〜Ｐ_Ｎにそれぞれ対応付けられた関係式ＥＱ３＿１（＝Ｌ_{ＭＡＸ＿１}（ｚ））〜ＥＱ３＿Ｎ（＝Ｌ_{ＭＡＸ＿Ｎ}（ｚ））からなる。

作成手段１６は、関係式ＥＱ１、関係式ＥＱ２および関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎ）を演算手段１５から受ける。また、作成手段１６は、作成すべき等高線ＣＴＲの高さｚおよび電力値Ｐ（Ｐ_１〜Ｐ_Ｎの少なくとも１つ）を外部から受ける。

作成手段１６は、作成すべき等高線ＣＴＲの高さｚ_ａと、作成すべき等高線ＣＴＲの電力値Ｐ_ｎとを受けると、関係式ＥＱ１（＝ｘ（ｚ））に基づいて高さｚ_ａに対応する受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇ（ｚ_ａ）を求めるとともに関係式ＥＱ２（＝ｙ（ｚ））に基づいて高さｚ_ａに対応する受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙ_Ｇ（ｚ_ａ）を求める。

次に、作成手段１６は、関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿ｎ）に基づいて電力値Ｐ_ｎに対応する関係式ＥＱ３＿ｎ（Ｌ_{ＭＡＸ＿ｎ}（ｚ））を検出し、その検出した関係式ＥＱ３＿ｎ（Ｌ_{ＭＡＸ＿ｎ}（ｚ））に基づいて高さｚ_ａに対応する最大距離Ｌ_{ＭＡＸ＿ｎ}（ｚ_ａ）を求める。

そうすると、作成手段１６は、受信電力の重心点（ｘ_Ｇ（ｚ_ａ），ｙ_Ｇ（ｚ_ａ））を中心とし、最大距離Ｌ_{ＭＡＸ＿ｎ}（ｚ_ａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲ１を作成する。

また、作成手段１６は、作成すべき等高線ＣＴＲの高さｚ_ａと、作成すべき等高線ＣＴＲの電力値Ｐ_ｎ，Ｐ_ｎ＋２とを受けると、関係式ＥＱ１，ＥＱ２に基づいて、上述した方法によって、それぞれ、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇ（ｚ_ａ）およびｙ座標ｙ_Ｇ（ｚ_ａ）を求める。

次に、作成手段１６は、電力値Ｐ_ｎに基づいて、上述した方法によって最大距離Ｌ_{ＭＡＸ＿ｎ}（ｚ_ａ）を求める。そして、作成手段１６は、関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿ｎ）に基づいて電力値Ｐ_ｎ＋２に対応する関係式ＥＱ３＿ｎ＋２（Ｌ_{ＭＡＸ＿ｎ＋２}（ｚ））を検出し、その検出した関係式ＥＱ３＿ｎ＋２（Ｌ_{ＭＡＸ＿ｎ＋２}（ｚ））に基づいて高さｚ_ａに対応する最大距離Ｌ_{ＭＡＸ＿ｎ＋２}（ｚ_ａ）を求める。

そうすると、作成手段１６は、受信電力の重心点（ｘ_Ｇ（ｚ_ａ），ｙ_Ｇ（ｚ_ａ））を中心とし、最大距離Ｌ_{ＭＡＸ＿ｎ}（ｚ_ａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲ１を作成するとともに受信電力の重心点（ｘ_Ｇ（ｚ_ａ），ｙ_Ｇ（ｚ_ａ））を中心とし、最大距離Ｌ_{ＭＡＸ＿ｎ＋２}（ｚ_ａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲ２を作成する。

作成手段１６は、作成すべき等高線ＣＴＲの高さｚ_ａと、作成すべき等高線ＣＴＲの電力値Ｐ_１〜Ｐ_ｎとを受けた場合、上述した方法によって、高さｚ_ａにおいて電力値Ｐ_１〜Ｐ_ｎをそれぞれ有する受信電力のｎ個の等高線ＣＴＲ１〜ＣＴＲｎを作成する。

また、作成手段１６は、作成すべき等高線ＣＴＲの高さとして高さｚ_ａ以外の高さｚ_ｊと、作成すべき等高線ＣＴＲの電力値Ｐ_ｎ（Ｐ_１〜Ｐ_Ｎの少なくとも１つ）とを受けた場合、上述した方法によって、高さｚ_ｊにおいて電力値Ｐ_ｎ（Ｐ_１〜Ｐ_Ｎの少なくとも１つ）を有する受信電力の少なくとも１つの等高線を作成する。

図１０は、等高線の概念図である。図１０を参照して、等高線ＣＴＲ１は、高さｚ_ｊにおいて受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、最大距離Ｌ_{ＭＡＸ＿１}を半径とする電力値Ｐ_１［ｄＢｍ］を有する等高線であり、等高線ＣＴＲ２は、高さｚ_ｊにおいて受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、最大距離Ｌ_{ＭＡＸ＿２}を半径とする電力値Ｐ_２［ｄＢｍ］を有する等高線であり、等高線ＣＴＲ３は、高さｚ_ｊにおいて受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、最大距離Ｌ_{ＭＡＸ＿３}を半径とする電力値Ｐ_３［ｄＢｍ］を有する等高線である。

図１０においては、高さｚ_ｊにおいて、３個の等高線ＣＴＲ１〜ＣＴＲ３が図示されているが、高さｚ_ｊにおいて少なくとも１つの等高線ＣＴＲが作成される。同様にして、高さｚ_１〜ｚ_３においても、それぞれ、少なくとも１つの等高線が作成される。

このように、この発明の実施の形態においては、任意の高さにおいて少なくとも１つの等高線ＣＴＲが作成される。

作成手段１６は、上述した方法によって、任意の高さｚにおける受信電力の等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲを推定手段１７へ出力する。

図１１は、受信電力と波源からの距離との関係を示す図である。図１１において、縦軸は、受信電力を表し、横軸は、波源Ｓ（受信電力の重心点ＲＳＳＩ＿Ｇ）からの距離を表す。

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

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

図１２は、受信電力の等高線ＣＴＲに基づいて任意の位置における受信電力の電力値を推定する方法を説明するための図である。

図１２を参照して、等高線ＣＴＲは、作成手段１６によって作成された等高線である。受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）から等高線ＣＴＲよりも３．１６ｋｍだけ遠い位置においては、受信電力の電力値は、等高線ＣＴＲの電力値をＰとすると、Ｐ−１０［ｄＢｍ］となる。

また、等高線ＣＴＲよりも３．１６ｋｍだけ受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）に近い位置においては、受信電力の電力値は、Ｐ＋１０［ｄＢｍ］となる。

推定手段１７は、受信電力の等高線ＣＴＲ（図１０参照）を作成手段１６から受ける。そして、推定手段１７は、受信電力の等高線ＣＴＲ（図１０参照）に基づいて、受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）から等高線ＣＴＲよりも３．１６ｋｍだけ遠い位置における受信電力の電力値をＰ−１０［ｄＢｍ］と推定し、その推定した電力値Ｐ−１０［ｄＢｍ］を有する受信電力の等高線ＣＴＲ１−１を作成する。

また、推定手段１７は、受信電力の等高線ＣＴＲ（図１０参照）に基づいて、等高線ＣＴＲよりも３．１６ｋｍだけ受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）に近い位置における受信電力の電力値をＰ＋１０［ｄＢｍ］と推定し、その推定した電力値Ｐ＋１０［ｄＢｍ］を有する受信電力の等高線ＣＴＲ１−２を作成する。

受信電力の等高線ＣＴＲ１−１，ＣＴＲ１−２は、等高線ＣＴＲと相似な形状を有する。

なお、推定手段１７は、単位距離当たりの受信電力の電力値の変化割合ＲＳＳＩ＿ｒａｔｅ（＝１０／３．１６＝３．１６［ｄＢｍ／ｋｍ］）を求め、変化割合ＲＳＳＩ＿ｒａｔｅを用いて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定してもよい。より具体的には、推定手段１７は、受信電力の重心点ＲＳＳＩ＿Ｇから遠ざかるように受信電力の等高線ＣＴＲから距離Ｌ_１だけ離れた位置における受信電力の電力値をＰ−ＲＳＳＩ＿ｒａｔｅ×Ｌ_１によって推定する。また、推定手段１７は、受信電力の重心点ＲＳＳＩ＿Ｇに近づくように受信電力の等高線ＣＴＲから距離Ｌ_１だけ離れた位置における受信電力の電力値をＰ＋ＲＳＳＩ＿ｒａｔｅ×Ｌ_１［ｄＢｍ］によって推定する。これによって、推定手段１７は、ある高さｚにおいて、水平方向の任意の位置における受信電力の電力値を推定できる。

推定手段１７は、電波の自由空間伝搬モデルに代えて奥村−秦モデルを用いて水平方向の任意の位置における受信電力の電力値を推定してもよい。

推定手段１７は、奥村−秦モデルを用いて受信電力の電力値を推定する場合、奥村−秦モデルにおける複数の電波伝搬モデルを地形データと対応付けて保持しており、標高地図または航空写真を外部から受ける。そして、推定手段１７は、標高地図または航空写真に基づいて抽出した地形データに最も近い電波伝搬モデルを複数の電波伝搬モデルから選択し、その選択した電波伝搬モデルを用いて受信電力の電力値を推定する。

図１３は、奥村−秦モデルにおける電界強度と距離との関係を示す図である。図１３において、縦軸は、電界強度を表し、横軸は、距離を表す。また、曲線ｋ４は、開放地における電界強度と距離との関係を示し、曲線ｋ５は、郊外における電界強度と距離との関係を示し、曲線ｋ６は、中小都市における電界強度と距離との関係を示し、曲線ｋ７は、大都市における電界強度と距離との関係を示す。

図１３に示す曲線ｋ４〜ｋ７は、次式によって表される。

式（２）において、Ａ，Ｂは、曲線ｋ４〜ｋ７に共通であり、ａ（ｈ_ｍ）およびＣは、エリア毎に異なる。また、Ｌ_ＣＨは、電界強度であり、ｄは、距離であり、ｈ_ｍは、２次利用者（端末装置２）のアンテナ高さである。

式（２）のＡ，Ｂは、次式によって表される。

式（３）において、ｆは、伝搬する電波の周波数であり、ｈ_ｂは、１次利用者（波源Ｓ）のアンテナ高さである。

電波の伝搬環境が開放地である場合（曲線ｋ４）、ａ（ｈ_ｍ）およびＣは、次式によって表される。

電波の伝搬環境が郊外である場合（曲線ｋ５）、ａ（ｈ_ｍ）およびＣは、次式によって表される。

電波の伝搬環境が中小都市である場合（曲線ｋ６）、ａ（ｈ_ｍ）およびＣは、次式によって表される。

電波の伝搬環境が大都市である場合（曲線ｋ７）、ａ（ｈ_ｍ）およびＣは、次式によって表される。

式（７）に示すように、大都市においては、周波数によってａ（ｈ_ｍ）が異なる。

推定手段１７は、式（２）〜式（７）を保持しており、地形データが開放地に最も近いとき、式（２），（３），（４）によって曲線ｋ４に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。

また、推定手段１７は、地形データが郊外に最も近いとき、式（２），（３），（５）によって曲線ｋ５に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。

更に、推定手段１７は、地形データが中小都市に最も近いとき、式（２），（３），（６）によって曲線ｋ６に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。

更に、推定手段１７は、地形データが大都市に最も近いとき、式（２），（３），（７）によって曲線ｋ７に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。

図１３に示す曲線ｋ４〜ｋ７においては、電界強度が距離に対して直線的に減少する領域があるので、電波の自由空間伝搬モデルにおいて説明したように、推定手段１７は、変化割合ＲＳＳＩ＿ｒａｔｅを算出し、その算出した変化割合ＲＳＳＩ＿ｒａｔｅを用いて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定してもよい。

このように、推定手段１７は、電波の自由空間伝搬モデルおよび奥村−秦モデル等の電波の伝搬モデルを用いて受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。

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

端末装置２は、受信感度以下の受信電力を取得できないため、受信感度以下の領域における受信電力の等高線ＣＴＲは、既に作成した受信電力の等高線ＣＴＲと電波の伝搬モデルとに基づいて作成する。ここで、電波の伝搬モデルとしては、上述した電波の自由空間伝搬モデルおよび奥村−秦モデル等が用いられる。

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

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

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

また、作成手段１６は、同様にして、奥村−秦モデルを用いて電力値ＰよりもΔＰだけ低下するときの距離を求め、その求めた距離だけ受信電力の重心点ＲＳＳＩ＿Ｇから遠くなる方向に離れ、かつ、受信電力の等高線ＣＴＲと相似形を有する受信電力の等高線ＣＴＲ＿ＬＳをＰ−ΔＰの電力値を有する等高線ＣＴＲとして作成する。

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

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

図１５は、図１に示す電力推定装置１の動作を説明するためのフローチャートである。図１５を参照して、電力推定装置１の動作が開始されると、ｎ＝１が設定され（ステップＳ１）、ｓ＝１が設定される（ステップＳ２）。

そして、抽出手段１３は、高さ方向の領域ＳＰＣ（ｓ）に含まれる複数のモニター情報を記録手段１２に記録されたモニター情報から抽出する（ステップＳ３）。

その後、抽出手段１３は、作成すべき受信電力の等高線ＣＴＲの電力値Ｐ_ｎに応じて処理範囲ＲＥＧ＿ＰＲＳを決定し（ステップＳ４）、その決定した処理範囲ＲＥＧ＿ＰＲＳ内に存在するｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを抽出する（ステップＳ５）。そして、抽出手段１３は、その抽出したｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを処理手段１４へ出力する（ステップＳ６）。

処理手段１４は、ｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋを抽出手段１３から受け、その受けたｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに基づいて受信電力の重心点ＲＳＳＩ＿Ｇ（ｓ）および最大距離Ｌ_ＭＡＸ（ｓ，Ｐ_ｎ）を算出する（ステップＳ７）。

引き続いて、処理手段１４は、ｋ個のモニター情報ＭＮＴ_１〜ＭＮＴ_ｋに基づいて領域ＳＰＣ（ｓ）の高さ位置ｚ_ｓ＝（ｚ_１＋・・・＋ｚ_ｋ）／ｋを算出する（ステップＳ８）。なお、処理手段１４は、ｚ_ｓ＝（ｗ_１ｚ_１＋ｗ_２ｚ_２＋・・・＋ｗ_ｋｚ_ｋ）／（ｗ_１＋ｗ_２＋・・・＋ｗ_ｋ）によって高さ位置ｚ_ｓを算出してもよい。ここで、ｗ_ｉ＝ＲＳＳＩ_ｉ（ｉ＝１〜ｋ）である。

その後、処理手段１４は、ｓ＝ｍ（高さ方向の領域ＳＰＣの最大個数）であるか否かを判定する（ステップＳ９）。

ステップＳ９において、ｓ＝ｍでないと判定されたとき、処理手段１４は、ｓ＝ｓ＋１を設定する（ステップＳ１０）。その後、一連の動作は、ステップＳ３へ戻り、ステップＳ９において、ｓ＝ｍであると判定されるまで、ステップＳ３〜ステップＳ１０が繰り返し実行される。

そして、ステップＳ９において、ｓ＝ｍであると判定されると、処理手段１４は、ｍ個の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離Ｌ_{ＭＡＸ＿１}（Ｐ_ｎ）〜Ｌ_{ＭＡＸ＿ｍ}（Ｐ_ｎ）およびｍ個の高さｚ_１〜ｚ_ｍを演算手段１５へ出力する。

演算手段１５は、ｍ個の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離Ｌ_{ＭＡＸ＿１}（Ｐ_ｎ）〜Ｌ_{ＭＡＸ＿ｍ}（Ｐ_ｎ）およびｍ個の高さｚ_１〜ｚ_ｍを処理手段１４から受け、その受けたｍ個の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離Ｌ_{ＭＡＸ＿１}（Ｐ_ｎ）〜Ｌ_{ＭＡＸ＿ｍ}（Ｐ_ｎ）およびｍ個の高さｚ_１〜ｚ_ｍに基づいて、上述した方法によって、重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇと高さｚとの関係式ＥＱ１と、重心点ＲＳＳＩ＿Ｇのｙ座標ｙ_Ｇと高さｚとの関係式ＥＱ２と、最大距離Ｌ_ＭＡＸと高さｚとの関係式ＥＱ３（等高線ＣＴＲの電力値Ｐ_ｎに対応した関係式ＥＱ３−ｎ）とを求める（ステップＳ１１）。そして、演算手段１５は、関係式ＥＱ１〜ＥＱ３を作成手段１６へ出力する。

作成手段１６は、関係式ＥＱ１〜ＥＱ３を演算手段１５から受ける。そして、作成手段１６は、任意の高さｚ_ａを外部から受けると、関係式ＥＱ１，ＥＱ２に基づいて、それぞれ、任意の高さｚ_ａにおける重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇ（ｚ_ａ）およびｙ座標ｙ_Ｇ（ｚ_ａ）を求め、関係式ＥＱ３（等高線ＣＴＲの電力値Ｐ_ｎに対応した関係式ＥＱ３−ｎ）に基づいて任意の高さｚ_ａにおける最大距離Ｌ_{ＭＡＸ＿ｎ}（ｚ_ａ）を求める（ステップＳ１２）。

そうすると、作成手段１６は、点（ｘ_Ｇ（ｚ_ａ），ｙ_Ｇ（ｚ_ａ））を中心とし、最大距離Ｌ_{ＭＡＸ＿ｎ}（ｚ_ａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲを作成する（ステップＳ１３）。

そして、作成手段１６は、ｎ＝Ｎ（Ｎは、作成すべき受信電力の等高線ＣＴＲの電力値Ｐ_ｎの最大個数）であるか否かを判定する（ステップＳ１４）。

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

そして、ステップＳ１４において、ｎ＝Ｎであると判定されると、一連の動作は、終了する。なお、図１５に示すフローチャートは、上述した分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々において実行される。

図１６は、図１５に示すフローチャートのステップの実行状態を示す概念図である。図１６を参照して、図１５に示すフローチャートのステップＳ３〜ステップＳ１０がｍ回実行されると、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐ_ｎにおいて、高さ方向のｍ個の領域ＳＰＣ（１）〜ＳＰＣ（ｍ）の全てについて、高さ位置ｚ、受信電力の重心点ＲＳＳＩ＿Ｇ（（ｘ_Ｇ（ｚ_１），ｙ_Ｇ（ｚ_１），ｘ_Ｇ（ｚ_２），ｙ_Ｇ（ｚ_２），・・・，ｘ_Ｇ（ｚ_ｍ），ｙ_Ｇ（ｚ_ｍ））および最大距離Ｌ_ＭＡＸ（ｚ_１），Ｌ_ＭＡＸ（ｚ_２），・・・，Ｌ_ＭＡＸ（ｚ_ｍ）が算出される（図１６の（ａ）参照）。

また、ステップＳ１１が実行されると、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐ_ｎにおいて、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇとの関係式ＥＱ１（＝ｘ_Ｇ（ｚ））、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙ_Ｇとの関係式ＥＱ２（＝ｙ_Ｇ（ｚ））、および高さｚおよび電力値Ｐ_ｎと最大距離Ｌ_ＭＡＸとの関係式ＥＱ３（＝Ｌ_{ＭＡＸ＿ｎ}（Ｐ_ｎ，ｚ））が求められる（図１６の（ｂ）参照）。ここで、Ｌ_{ＭＡＸ＿ｎ}（Ｐ_ｎ，ｚ）は、等高線ＣＴＲの１つの電力値Ｐ_ｎに対応し、かつ、最大距離Ｌ_ＭＡＸと高さｚとの関係を示す関係式である。

更に、ステップＳ１２が実行されると、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐ_ｎにおいて、任意の高さｚ_ａにおける受信電力の重心点ＲＳＳＩ＿Ｇの座標（（ｘ_Ｇ（ｚ_ａ），ｙ_Ｇ（ｚ_ａ））および任意の高さｚ_ａにおける最大距離Ｌ_{ＭＡＸ＿ｎ}（Ｐ_ｎ，ｚ_ａ）が求められる（図１６の（ｃ）参照）。この場合、重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇ（ｚ_ａ）は、関係式ＥＱ１＝ｘ_Ｇ（ｚ）にｚ＝ｚ_ａを代入して求められ、重心点ＲＳＳＩ＿Ｇのｙ座標ｙ_Ｇ（ｚ_ａ）は、関係式ＥＱ２＝ｙ_Ｇ（ｚ）にｚ＝ｚ_ａを代入して求められ、最大距離Ｌ_ＭＡＸｎ（Ｐ_ｎ，ｚ_ａ）は、関係式ＥＱ３（等高線ＣＴＲの電力値Ｐ_ｎに対応した関係式ＥＱ３−ｎ）＝Ｌ_{ＭＡＸ＿ｎ}（Ｐ_ｎ，ｚ）にｚ＝ｚ_ａを代入して求められる。

更に、ステップＳ１３が実行されると、高さｚ_ａにおいて１つの電力値Ｐ_ｎを有する等高線ＣＴＲ＿ｎが作成される。等高線ＣＴＲ＿ｎは、受信電力の重心点（（ｘ_Ｇ（ｚ_ａ），ｙ_Ｇ（ｚ_ａ））を中心とし、最大距離Ｌ_{ＭＡＸ＿ｎ}（Ｐ_ｎ，ｚ_ａ）を半径とする円形形状を有する（図１６の（ｄ）参照）。

更に、ステップＳ２〜ステップＳ１５がＮ回実行されると、高さｚ_ａにおいて、それぞれ、電力値Ｐ_１〜Ｐ_Ｎを有するＮ個の等高線ＣＴＲ＿１〜ＣＴＲ＿Ｎが作成される（図１６の（ｅ）参照）。等高線ＣＴＲ＿１〜ＣＴＲ＿Ｎは、受信電力の重心点（（ｘ_Ｇ（ｚ_ａ），ｙ_Ｇ（ｚ_ａ））を中心とし、それぞれ、最大距離Ｌ_{ＭＡＸ＿１}（Ｐ_１，ｚ_ａ）〜Ｌ_{ＭＡＸ＿Ｎ}（Ｐ_Ｎ，ｚ_ａ）を半径とする円形形状を有する。そして、最大距離Ｌ_{ＭＡＸ＿１}（Ｐ_１，ｚ_ａ）〜Ｌ_{ＭＡＸ＿Ｎ}（Ｐ_Ｎ，ｚ_ａ）は、それぞれ、電力値Ｐ_１〜Ｐ_ｎに対応する関係式ＥＱ３＿１＝Ｌ_{ＭＡＸ＿１}（Ｐ_１，ｚ）〜ＥＱ３＿Ｎ＝Ｌ_{ＭＡＸ＿Ｎ}（Ｐ_Ｎ，ｚ）を用いて求められる。

従って、図１５に示すフローチャートは、高さｚ_ａにおいて、１つの電力値Ｐ_ｎ毎に受信電力の等高線ＣＴＲを作成する動作を示すものである。

図１７は、図１に示す電力推定装置１の動作を説明するための別のフローチャートである。

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

図１７を参照して、電力推定装置１の動作が開始されると、上述したステップＳ１〜Ｓ１１，Ｓ１４，Ｓ１５が順次実行される。

そして、ステップＳ１４において、ｎ＝Ｎであると判定されると、作成手段１６は、関係式ＥＱ１〜ＥＱ３に基づいて、任意の高さｚ_ａにおける任意の電力値Ｐ_ｎを有する受信電力の複数の等高線ＣＴＲを作成する（ステップＳ１６）。これによって、電力推定装置１の動作が終了する。

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

図１８を参照して、図１７に示すステップＳ１４において、ｎ＝Ｎであると判定されると、作成手段１６は、高さｚ_ａと電力値Ｐ_ｎとのＶ（Ｖは，２≦Ｖ≦Ｎを満たす整数）個の組（ｚ_ａ１，Ｐ_１），（ｚ_ａ２，Ｐ_２），・・・，（ｚ_ａＶ，Ｐ_Ｖ）を外部から受ける（ステップＳ１６１）。

そして、作成手段１６は、ｖ＝１を設定する（ステップＳ１６２）。その後、作成手段１６は、関係式ＥＱ１，ＥＱ２に基づいて、それぞれ高さｚ_ａｒにおける受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇ（ｚ_ａｒ）およびｙ座標ｙ_Ｇ（ｚ_ａｒ）を求め、関係式ＥＱ３（等高線ＣＴＲの電力値Ｐ_ｖに対応する関係式ＥＱ３−ｖ）に基づいて、電力値Ｐ_ｖおよび高さｚ_ａｖにおける最大距離Ｌ_{ＭＡＸ＿ｖ}（Ｐ_ｖ，ｚ_ａｖ）を求める（ステップＳ１６３）。

そうすると、作成手段１６は、受信電力の重心点（ｘ_Ｇ（ｚ_ａｒ），ｙ_Ｇ（ｚ_ａｒ））を中心とし、最大距離Ｌ_{ＭＡＸ＿ｖ}（Ｐ_ｖ，ｚ_ａｖ）を半径とする受信電力の等高線ＣＴＲ＿ｖを作成する（ステップＳ１６４）。

そして、作成手段１６は、ｒ＝Ｒであるか否かを判定する（ステップＳ１６５）。ステップＳ１６５において、ｖ＝Ｖでないと判定されたとき、作成手段１６は、ｖ＝ｖ＋１を設定する（ステップＳ１６６）。その後、一連の動作は、ステップＳ１６３に移行し、ステップＳ１６５において、ｖ＝Ｖであると判定されるまで、ステップＳ１６３〜ステップＳ１６６が繰り返し実行される。

そして、ステップＳ１６５において、ｖ＝Ｖであると判定されると、一連の動作は、図１７の“終了”へ移行する。なお、図１７に示すフローチャート（図１８に示すフローチャートを含む）は、上述した分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々において実行される。

図１９は、図１７および図１８に示すフローチャートのステップの実行状態を示す概念図である。図２０は、図１７および図１８に示すフローチャートの実行が完了した状態を示す概念図である。

図１９を参照して、図１７に示すフローチャートのステップＳ３〜ステップＳ１０がｍ回実行されると、図１６において説明したとおり、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐ_ｎにおいて、高さ方向のｍ個の領域ＳＰＣ（１）〜ＳＰＣ（ｍ）の全てについて、高さ位置ｚ、受信電力の重心点ＲＳＳＩ＿Ｇ（（ｘ_Ｇ（ｚ_１），ｙ_Ｇ（ｚ_１），ｘ_Ｇ（ｚ_２），ｙ_Ｇ（ｚ_２），・・・，ｘ_Ｇ（ｚ_ｍ），ｙ_Ｇ（ｚ_ｍ））および最大距離Ｌ_ＭＡＸ（ｚ_１），Ｌ_ＭＡＸ（ｚ_２），・・・，Ｌ_ＭＡＸ（ｚ_ｍ）が算出される（図１９の（ａ）参照）。

また、図１７に示すフローチャートのステップＳ２〜Ｓ１１，Ｓ１４，Ｓ１５がＮ回実行されると、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ_Ｇとの関係式ＥＱ１（＝ｘ_Ｇ（ｚ））、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙ_Ｇとの関係式ＥＱ２（＝ｙ_Ｇ（ｚ））、および電力値Ｐ_１〜Ｐ_Ｎに対応して、それぞれ、高さｚおよび電力値Ｐ_ｎと最大距離Ｌ_ＭＡＸとの関係式ＥＱ３＿１（＝Ｌ_{ＭＡＸ＿１}（Ｐ_１，ｚ））〜ＥＱ３＿Ｎ（＝Ｌ_{ＭＡＸ＿Ｎ}（Ｐ_Ｎ，ｚ））が求められる（図１９の（ｂ）参照）。ここで、Ｌ_{ＭＡＸ＿１}（Ｐ_１，ｚ）〜Ｌ_{ＭＡＸ＿Ｎ}（Ｐ_Ｎ，ｚ）は、それぞれ、等高線ＣＴＲの電力値Ｐ_１〜Ｐ_Ｎに対応し、かつ、最大距離Ｌ_ＭＡＸと高さｚとの関係を示す関係式である。

更に、図１８に示すステップＳ１６３が実行されると、関係式ＥＱ１，ＥＱ２に基づいて、高さｚ_ａｖにおける受信電力の重心点の座標（ｘ_Ｇ（ｚ_ａｖ），ｙ_Ｇ（ｚ_ａｖ））が求められ、電力値Ｐ_ｖに対応する関係式ＥＱ３＿ｖ（ＥＱ３＿１〜ＥＱ３＿Ｎのいずれか）に基づいて、高さｚ_ａｖにおける最大距離Ｌ_{ＭＡＸ＿ｖ}（Ｐ_ｖ，ｚ_ａｖ）が求められる（図１９の（ｃ）参照）。

更に、図１８に示すステップＳ１６４が実行されると、高さｚ_ａｖにおいて１つの電力値Ｐ_ｖを有する等高線ＣＴＲ＿ｖが作成される。等高線ＣＴＲ＿ｖは、受信電力の重心点（（ｘ_Ｇ（ｚ_ａｖ），ｙ_Ｇ（ｚ_ａｖ））を中心とし、最大距離Ｌ_{ＭＡＸ＿ｖ}（Ｐ_ｖ，ｚ_ａｖ）を半径とする円形形状を有する（図１９の（ｄ）参照）。

更に、図１７および図１８に示すフローチャートの実行が完了した状態においては、高さｚ_ａ１，ｚ_ａ２，ｚ_ａ３，・・・，ｚ_ａＶにおいて、それぞれ、電力値Ｐ_１，Ｐ_２，Ｐ_３，・・・，Ｐ_Ｖを有する受信電力の等高線ＣＴＲ＿１，ＣＴＲ＿２，ＣＴＲ＿３，・・・，ＣＴＲ＿Ｖが作成される。等高線ＣＴＲ＿１は、受信電力の重心点（（ｘ_Ｇ（ｚ_ａ１），ｙ_Ｇ（ｚ_ａ１））を中心とし、最大距離Ｌ_{ＭＡＸ＿１}（Ｐ_１，ｚ_ａ１）を半径とする円形形状を有する。等高線ＣＴＲ＿２は、受信電力の重心点（（ｘ_Ｇ（ｚ_ａ２），ｙ_Ｇ（ｚ_ａ２））を中心とし、最大距離Ｌ_{ＭＡＸ＿２}（Ｐ_２，ｚ_ａ２）を半径とする円形形状を有する。等高線ＣＴＲ＿３は、受信電力の重心点（（ｘ_Ｇ（ｚ_ａ３），ｙ_Ｇ（ｚ_ａ３））を中心とし、最大距離Ｌ_{ＭＡＸ＿３}（Ｐ_３，ｚ_ａ３）を半径とする円形形状を有する。以下、同様にして、等高線ＣＴＲ＿Ｖは、受信電力の重心点（（ｘ_Ｇ（ｚ_ａＶ），ｙ_Ｇ（ｚ_ａＶ））を中心とし、最大距離Ｌ_{ＭＡＸ＿Ｖ}（Ｐ_Ｖ，ｚ_ａＶ）を半径とする円形形状を有する（図２０参照）。

受信電力の等高線ＣＴＲ＿１〜ＣＴＲ＿Ｖは、異なる高さｚ_ａ１〜ｚ_ａＶにおいて作成され、それぞれ、異なる半径（最大距離Ｌ_{ＭＡＸ＿１}（Ｐ_１，ｚ_ａ１）〜Ｌ_{ＭＡＸ＿Ｖ}（Ｐ_Ｖ，ｚ_ａＶ））および異なる電力値Ｐ_１〜Ｐ_Ｖを有する。

従って、図１７および図１８に示すフローチャートは、任意の高さにおいて任意の電力値を有する受信電力の複数の等高線ＣＴＲを作成する動作を示すものである。

図２１は、図１に示す電力推定装置１の別の動作を説明するためのフローチャートである。

図２１を参照して、電力推定装置１の動作が開始されると、作成手段１６は、図１５に示すフローチャート、または図１７に示すフローチャート（図１８に示すフローチャートを含む）を実行する（ステップＳ２１）。そして、作成手段１６は、作成した受信電力の等高線ＣＴＲを推定手段１７へ出力する。

推定手段１７は、受信電力の等高線ＣＴＲを作成手段１６から受け、その受けた受信電力の等高線ＣＴＲに基づいて受信電力を推定する（ステップＳ２２）。この場合、推定手段１７は、上述した方法によって、電波の伝搬モデルを用いて受信電力を推定する。そして、ステップＳ２２の後、電力推定装置１の動作が終了する。

なお、電力推定装置１の動作は、ソフトウェアによって実行されてもよい。この場合、電力推定装置１は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）およびＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）を備える。

ＲＯＭは、図１５に示すフローチャートからなるプログラムＰｒｏｇ＿Ａ、図１７に示すフローチャート（図１８に示すフローチャートを含む）からなるプログラムＰｒｏｇ＿Ｂ、および図２１に示すフローチャートからなるプログラムＰｒｏｇ＿Ｃのいずれかを格納する。

そして、ＣＰＵは、プログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿ＣのいずれかをＲＯＭから読み出し、その読み出したプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかを実行して受信電力の等高線ＣＴＲの作成、または受信電力の等高線ＣＴＲの作成および受信電力の推定を行う。

この場合、ＲＡＭは、端末装置２から受信したモニター情報ＭＮＴおよび関係式ＥＱ１〜ＥＱ３等の記憶に用いられるとともに各種の演算に用いられる。

また、プログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかは、ＣＤ，ＤＶＤ等の記録媒体に記録されて流通してもよい。この場合、ＣＰＵは、装着された記録媒体からプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかを読み出して実行し、受信電力の等高線ＣＴＲの作成、または受信電力の等高線ＣＴＲの作成および受信電力の推定を行う。従って、プログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかを記録したＣＤ，ＤＶＤ等の記録媒体は、コンピュータ（ＣＰＵ）が読み取り可能な記録媒体である。

上記においては、受信電力の等高線ＣＴＲは、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離Ｌ_ＭＡＸを半径とする円形形状を有すると説明したが、この発明の実施の形態においては、これに限らず、受信電力の等高線ＣＴＲは、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離Ｌ_ＭＡＸを長径とする楕円形状を有していてもよい。

図２２は、楕円形状を有する受信電力の等高線ＣＴＲの作成方法を説明するための図である。

図２２において、白四角は、電力値Ｐ以上の受信電力を有する測定点を示し、黒四角は、電力値Ｐ未満の受信電力を有する測定点を示し、白丸は、受信電力の重心点ＲＳＳＩ＿Ｇを示す。

図２２を参照して、作成手段１６は、上述した方法によって、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離Ｌ_ＭＡＸを半径とする円形形状を有する受信電力の等高線ＣＴＲ５を作成する。

そして、作成手段１６は、受信電力の等高線ＣＴＲ５に基づいて、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離Ｌ_ＭＡＸを長径とするとともに距離Ｌ’_ＭＡＸを短径とする楕円形状を有する受信電力の等高線ＣＴＲ６を作成する。ここで、距離Ｌ’_ＭＡＸは、電力値Ｐ以上の受信電力を有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇからの距離が最大距離Ｌ_ＭＡＸの次に長い測定点２−２と受信電力の重心点ＲＳＳＩ＿Ｇとの距離である。

また、この発明の実施の形態においては、受信電力の等高線ＣＴＲは、多角形の形状を有していてもよい。

図２３は、多角形の形状を有する受信電力の等高線ＣＴＲの作成方法を説明するための図である。

図２３において、白四角は、電力値Ｐ以上の受信電力を有する測定点を示し、黒四角は、電力値Ｐ未満の受信電力を有する測定点を示し、白丸は、受信電力の重心点ＲＳＳＩ＿Ｇを示す。

図２３を参照して、作成手段１６は、上述した方法によって、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離Ｌ_ＭＡＸを半径とする円形形状を有する受信電力の等高線ＣＴＲ５を作成する。

そして、作成手段１６は、円形形状を有する受信信号の等高線ＣＴＲ５に内接する正八角形３０を形成する（図２３の（ａ）参照）。

その後、作成手段１６は、正八角形の各頂点が、電力値Ｐ以上の受信電力を有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇからできるだけ遠い位置に存在する測定点２−１〜２−８に位置するように正八角形３０を変形し、八角形の形状を有する受信電力の等高線ＣＴＲ７を作成する（図２３の（ｂ）参照）。

作成手段１６は、図２３において説明した作成方法と同じ作成方法によって、八角形以外の多角形の形状を有する受信電力の等高線ＣＴＲを作成する。

このように、作成手段１６は、円形形状を有する受信電力の等高線ＣＴＲに限らず、楕円形状を有する受信電力の等高線ＣＴＲまたは多角形の形状を有する受信電力の等高線ＣＴＲを作成する。

従って、図１５に示すフローチャートのステップＳ１３または図１８に示すフローチャートのステップＳ１６４において、作成手段１６は、円形形状を有する受信電力の等高線ＣＴＲに代えて、楕円形状を有する受信電力の等高線ＣＴＲまたは多角形の形状を有する受信電力の等高線ＣＴＲを作成してもよく、一般的には、作成手段１６は、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、受信電力の重心点ＲＳＳＩ＿Ｇからの距離が最も長い最大距離Ｌ_ＭＡＸを半径とする円形形状を有する第１の等高線、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離Ｌ_ＭＡＸを長径とする楕円形状を有する第２の等高線、および受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離Ｌ_ＭＡＸを中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを任意の高さｚにおいて作成すればよい。

モニター情報ＭＮＴは、共用周波数ｆ_{ｃｏｍ＿ｉ}を含むため、電力推定装置１は、上述した方法によって、各共用周波数ｆ_{ｃｏｍ＿ｉ}を用いて電波を送信した場合における受信電力の等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲに基づいて受信電力を推定する。

また、電力推定装置１は、例えば、端末装置２が１５：００に無線通信を行う場合、１４：００〜１５：００までの間に複数の端末装置２から送信された複数のモニター情報ＭＮＴに基づいて、上述した方法によって、任意の高さｚにおける受信電力の等高線ＣＴＲを作成し（端末装置２における電波の受信感度以下の領域における受信電力の等高線ＣＴＲの作成も含む）、その作成した受信電力の等高線ＣＴＲに基づいて、各場所における受信電力を推定し、その推定結果に基づいて、共用周波数ｆ_{ｃｏｍ＿ｉ}を用いて無線通信を行うときの通信条件（共用周波数ｆ_ｃｏｍ、場所および送信電力等）を決定する。

そして、各端末装置２は、電力推定装置１によって決定された通信条件を用いて無線通信を行う。

その結果、電力推定装置１によって推定された受信電力は、フェージングを抑制して推定された受信電力であるため、送信電力は、フェージングが抑制された電波環境において信号を送信先へ送信できるように決定される。従って、スループットを高くして信号を送信先へ送信できる。

また、上記においては、分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々において、高さ方向の領域ＳＰＣ１〜ＳＰＣｍの各々は、高さ方向の同じ長さｈを有すると説明したが、この発明の実施の形態においては、これに限らず、建物２０の高さＨ_Ｂ（または平均高さ）以下の位置に存在する領域ＳＰＣ１，ＳＰＣ２においては、高さ方向に長さｈ１を有し、建物２０の高さＨ_Ｂ（または平均高さ）よりも高い位置に存在する領域ＳＰＣｍにおいては、高さ方向に長さｈ１よりも長い長さｈ２を有するようにしてもよい。

建物２０の高さＨ_Ｂ（または平均高さ）よりも高い領域においては、電波が伝搬するときの障害物が少ないので、高さ方向に長さを短くしなくても、電波の受信電力の電力分布を正確に推定することができる。

一方、建物２０の高さＨ_Ｂ（または平均高さ）よりも低い領域においては、電波が伝搬するときの障害物が多いので、高さ方向に長さを短くした方が電波の受信電力の電力分布を正確に推定することができる。

図２４は、この発明の実施の形態によるデータ構造の概略図である。図２４を参照して、データ構造Ｄ＿ＳＴＲ１は、端末装置ＩＤ、時間情報、端末装置２の位置情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤおよび受信電力が相互に対応付けられた構造を有する。端末装置ＩＤ、時間情報、端末装置２の位置情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤおよび受信電力は、モニター情報ＭＮＴを構成する。

端末装置ＩＤは、端末装置２の識別情報である。時間情報は、受信電力が検出された時刻を示す。そして、時間情報は、年／月／日／時間／分／秒（ＹＹＹＹ／ＭＭ／ＤＤ／ＨＨ／ＭＭ／ＳＳ）の形式によって表される。端末装置２の位置情報は、緯度と、経度と、高さとを含む。緯度および経度は、ＧＰＳによって検出され、高さは、高度計によって検出される。そして、高さは、端末装置２が静止または移動している地面からの高さを示す。

周波数帯域ＩＤは、複数の端末装置２が共用する周波数帯域の識別情報である。サブキャリアＮｏ．は、サブキャリアの番号を示す。波源ＩＤは、１次利用者の電波を発する波源の識別情報である。受信電力は、各端末装置２が受信した電波の受信電力である。

端末装置ＩＤは、抽出手段１３が、対応する受信電力が同じ端末装置２で検出されたか否かを判定するのに用いられる。時間情報は、抽出手段１３が、対応する受信電力が同じ時刻に検出されたか否かを判定するのに用いられる。端末装置２の位置情報は、抽出手段１３が、対応する受信電力が同じ位置で検出されたか否か、または対応する受信電力が所定の領域ＳＰＣ内で検出された否かを判定するのに用いられる。また、位置情報は、処理手段１４が、複数の端末装置２間の距離、端末装置２と波源との距離および受信電力の重心点から最も遠い位置に存在する端末装置を検出するのに用いられる。

周波数帯域ＩＤは、抽出手段１３が、同じ周波数帯域の受信電力を抽出するために用いられる。サブキャリアＮｏ．は、抽出手段１３が、対応する受信電力を有する電波が伝送されたサブキャリアを判定するのに用いられる。波源ＩＤは、抽出手段１３が、同じ波源から送信された電波の受信電力を抽出するのに用いられる。

受信電力は、処理手段１４が、複数の受信電力のうちの最大の受信電力を検出するのに用いられる。また、受信電力は、作成手段１６が、作成すべき電力の等高線の電力値以上の受信電力を検出するのに用いられる。

複数のモニター情報ＭＮＴに含まれる複数の受信電力のうち、受信電力の最大値を含むｋ個のモニター情報は、抽出手段１３によって複数のモニター情報から抽出されるのに用いられる。

抽出手段１３によって抽出されたｋ個のモニター情報は、処理手段１４が受信電力によって重み付けされたｋ個の端末装置２の位置の平均である受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）を求めるのに用いられる。

抽出手段１３によって抽出されたｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力と、受信電力の等高線における電力値以上の受信電力に対応付けられた位置情報とは、処理手段１４が、受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）から最も遠い位置に存在する端末装置と受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）との距離である最大距離Ｌ_ＭＡＸを求めるのに用いられる。

処理手段１４によって求められた受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）は、演算手段１５が受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）と高さｚとの関係を示す第１の関係式（関係式ＥＱ１，ＥＱ２）を求めるのに用いられる。

処理手段１４によって求められた最大距離Ｌ_ＭＡＸは、演算手段１５が、作成すべき受信電力の等高線の電力値に対応し、かつ、最大距離Ｌ_ＭＡＸと高さｚとの関係を示す第２の関係式（関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎを含む））を求めるのに用いられる。

演算手段１５によって求められた第１の関係式（関係式ＥＱ１，ＥＱ２）および第２の関係式（関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎを含む））は、作成手段１６が、任意の高さにおいて、任意の電力値を有する受信電力の等高線を作成するのに用いられる。

第１の等高線は、受信電力の重心点を中心とし、かつ、最大距離を半径とする円形形状を有し、第２の等高線は、受信電力の重心点を中心とし、かつ、最大距離を長径とする楕円形状を有し、第３の等高線は、受信電力の重心点を中心とし、かつ、最大距離を中心から１つの頂点までの距離とする多角形の形状を有する。

好ましくは、抽出手段１３によって抽出されたｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値よりも小さい受信電力と、受信電力の等高線における電力値よりも小さい受信電力に対応付けられた位置情報とは、推定手段１７が、受信電力の等高線に基づいて、対象領域ＳＰＣ内の所望の場所において電波の受信電力を推定するのに用いられる。

好ましくは、受信電力の等高線の電力値を有する受信電力は、作成手段１６が、端末装置２における電波の受信感度以下の領域において、既に作成した受信電力の等高線上の電力値よりも所望の電力値だけ低い受信電力になる電力低下距離を電波の伝搬モデルに基づいて求め、既に作成した受信電力の等高線から電力低下距離だけ離れ、かつ、既に作成した受信電力の等高線と相似形を有する等高線を電波の受信感度以下の領域における受信電力の等高線として作成するのに用いられる。

データ構造Ｄ＿ＳＴＲ１は、電力推定装置１（またはプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれか）が上述した方法によって受信電力の等高線を作成するのに用いられる。そして、データ構造Ｄ＿ＳＴＲ１は、端末装置ＩＤ、時間情報、端末装置２の位置情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤおよび受信電力の全てを含む必要はなく、端末装置２の位置情報、周波数帯域ＩＤおよび受信電力を少なくとも含み、端末装置２の位置情報、周波数帯域ＩＤおよび受信電力が相互に対応付けられた構造を有していればよい。データ構造Ｄ＿ＳＴＲ１が端末装置２の位置情報、周波数帯域ＩＤおよび受信電力を少なくとも含んでいれば、電力推定装置１（またはプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれか）が、上述した方法によって、任意の高さにおける受信電力の等高線を作成することができるからである。この場合、端末装置２の位置情報、周波数帯域ＩＤおよび受信電力は、モニター情報を構成する。

図２５は、この発明の実施の形態による別のデータ構造の概略図である。この発明の実施の形態によるデータ構造は、図２５に示すデータ構造Ｄ＿ＳＴＲ２であってもよい。

図２５を参照して、データ構造Ｄ＿ＳＴＲ２は、推定１次利用者ＩＤ、時間情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤ、等高線の中心位置、等高線の電力値および等高線の半径が相互に対応付けられた構造を有する。周波数帯域ＩＤ、サブキャリアＮｏ．および波源ＩＤについては、図２４において説明したとおりである。

時間情報は、受信電力の等高線が作成された時刻を示す。そして、時間情報は、年／月／日／時間／分／秒（ＹＹＹＹ／ＭＭ／ＤＤ／ＨＨ／ＭＭ／ＳＳ）の形式によって表される。推定１次利用者ＩＤは、推定された１次利用者の識別情報である。等高線の中心位置は、緯度と、経度と、高さとを含む。緯度および経度は、ＧＰＳによって検出され、高さは、高度計によって検出される。そして、高さは、等高線の中心位置の地面からの高さを示す。等高線の電力値は、作成された等高線上の電力値である。等高線の半径は、作成された等高線が円形状である場合の円の半径であり、受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）から等高線までの距離である。

時間情報は、作成手段１６が受信電力の等高線を作成した時刻を特定するのに用いられる。周波数帯域ＩＤは、作成手段１６が、同じ周波数帯域の受信電力に基づいて受信電力の等高線を作成したことを示すのに用いられる。サブキャリアＮｏ．は、作成手段１６が受信電力の等高線を作成するときの元になった受信電力を有する電波が伝送されたサブキャリアを判定するのに用いられる。波源ＩＤは、作成手段１６が作成した受信電力の等高線の波源を特定するのに用いられる。

周波数帯域ＩＤは、抽出手段１３が、ｋ個の端末装置が共用する共用周波数を識別するのに用いられる。等高線の中心位置を示す位置情報は、ｋ個の端末装置から受信したｋ個のモニター情報に含まれるｋ個の受信電力とｋ個の端末装置のｋ個の位置情報とに基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を処理手段１４が演算した結果である受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）からなる。等高線の半径は、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）から最も遠い位置に存在する端末装置と受信電力の重心点（ｘ_Ｇ，ｙ_Ｇ）との距離を処理手段１４が演算した結果である最大距離Ｌ_ＭＡＸからなる。等高線の電力値は、周波数帯域ＩＤ（共用周波数）を有する電波の電力値であり、作成手段１６によって検出された電力値である。

好ましくは、共用周波数が既に割り当てられている１次利用者の波源の識別情報は、作成手段１６によって、共用周波数、等高線の中心位置を示す位置情報、等高線の電力値および等高線の半径に更に対応付けられた構造を有する。

データ構造Ｄ＿ＳＴＲ２は、推定１次利用者ＩＤ、時間情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤ、等高線の中心位置、等高線の電力値および等高線の半径の全てを含む必要はなく、周波数帯域ＩＤ、等高線の中心位置を示す位置情報、等高線の電力値および等高線の半径を少なくとも含み、周波数帯域ＩＤ、等高線の中心位置を示す位置情報、等高線の電力値および等高線が相互に対応付けられた構造を有していればよい。データ構造Ｄ＿ＳＴＲ２が周波数帯域ＩＤ、等高線の中心位置を示す位置情報、等高線の電力値および等高線を少なくとも含んでいれば、作成手段１６によって作成された受信電力の等高線を表すことができるからである。

データ構造Ｄ＿ＳＴＲ２は、電力推定装置１（またはプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれか）が上述した方法によって作成した受信電力の等高線を表すのに用いられる。データ構造Ｄ＿ＳＴＲ２は、コンピュータに入力され、ＣＰＵによって記憶手段に記憶される。そして、データ構造Ｄ＿ＳＴＲ２は、ＣＰＵによって表示手段に表示され、または紙媒体等に印刷される。また、データ構造Ｄ＿ＳＴＲ２は、画像ソフトを実行したＣＰＵによって受信電力の等高線を示す画像を表示手段に表示するのに用いられてもよい。更に、データ構造Ｄ＿ＳＴＲ２は、受信電力の等高線から所望の距離だけ離れた位置における電波の電力値をＣＰＵが演算するのに用いられてもよい。

上記においては、電力推定装置１は、複数の端末装置２において検出された複数のモニター情報ＭＮＴに基づいて任意の高さにおける受信電力の等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲに基づいて各場所における受信電力を推定すると説明したが、この発明の実施の形態においては、これに限らず、電力推定装置１は、上述した方法によって、複数のモニター情報ＭＮＴに基づいて任意の高さにおける受信電力の等高線ＣＴＲを作成するだけでもよい。異なる電力値Ｐを有する受信電力の複数の等高線ＣＴＲを作成すれば、その作成された受信電力の複数の等高線ＣＴＲに基づいて、最も大きい電力値Ｐを有する受信電力の等高線ＣＴＲよりも波源Ｓ側の領域ＲＥＧ１、および電力値Ｐが異なる２つの受信電力の等高線ＣＴＲの間の領域ＲＥＧ２における受信電力が分かるので、これらの領域ＲＥＧ１，ＲＥＧ２における受信電力を推定できるからである。

従って、この発明の実施の形態による電力推定装置は、電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線を作成する電力推定装置であって、複数の端末装置から送信され、かつ、各々が端末装置の位置を示す位置情報と端末装置における電波の受信電力とを含む複数のモニター情報から、受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する抽出手段と、抽出手段によって抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点から最も遠い位置に存在する端末装置と受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する処理手段と、処理手段によって求められたｍ個の受信電力の重心点およびｍ個の最大距離と、ｍ個の高さとに基づいて受信電力の重心点と高さとの関係を示す第１の関係式と、受信電力の複数の電力値に対応付けられ、かつ、各々が最大距離と前記高さとの関係を示す複数の第２の関係式とを演算する演算手段と、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、複数の第２の関係式のうち、作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線を作成する作成手段とを備えていればよい。

また、この発明の実施の形態による無線通信システムは、任意の高さにおいて受信電力の等高線ＣＴＲを作成する電力推定装置を備えていればよい。

更に、この発明の実施の形態によるプログラムは、電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線の作成をコンピュータに実行させるためのプログラムであって、抽出手段が、複数の端末装置から送信され、かつ、各々が端末装置の位置を示す位置情報と端末装置における電波の受信電力とを含む複数のモニター情報から、受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する第１のステップと、処理手段が、第１のステップにおいて抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点から最も遠い位置に存在する端末装置と受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する第２のステップと、演算手段が、第２のステップにおいて求められたｍ個の受信電力の重心点およびｍ個の最大距離と、ｍ個の高さとに基づいて受信電力の重心点と高さとの関係を示す第１の関係式と、受信電力の複数の電力値に対応付けられ、かつ、各々が最大距離と高さとの関係を示す複数の第２の関係式とを演算する第３のステップと、作成手段が、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、複数の第２の関係式のうち、作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線を作成する第４のステップとをコンピュータに実行させればよい。

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

この発明は、電力推定装置、それを備えた無線通信システム、コンピュータに実行させるためのプログラムおよびデータ構造に適用される。

１ 電力推定装置、２ 端末装置、３ 基地局、４ 有線ケーブル、１０ 無線通信システム、１１ 受信手段、１２ 記録手段、１３ 抽出手段、１４ 処理手段、１５ 演算手段、１６ 作成手段、１７ 推定手段。

Claims (16)

1. 電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線を作成する電力推定装置であって、
   前記複数の端末装置から送信され、かつ、各々が前記端末装置の位置を示す位置情報と前記端末装置における前記電波の受信電力とを含む複数のモニター情報から、前記受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する抽出手段と、
   前記抽出手段によって抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、前記ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、前記受信電力の重心点から最も遠い位置に存在する端末装置と前記受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する処理手段と、
   前記処理手段によって求められたｍ個の受信電力の重心点およびｍ個の最大距離と、前記ｍ個の高さとに基づいて前記受信電力の重心点と前記高さとの関係を示す第１の関係式と、前記受信電力の複数の電力値に対応付けられ、かつ、各々が前記最大距離と前記高さとの関係を示す複数の第２の関係式とを演算する演算手段と、
   前記第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、前記複数の第２の関係式のうち、前記作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて前記任意の高さにおける前記最大距離を求め、その求めた受信電力の重心点および最大距離を用いて前記作成すべき受信電力の等高線を作成する作成手段とを備える電力推定装置。
2. 前記作成手段は、前記第１の関係式に基づいて任意の高さにおける受信電力の重心点を求めるとともに前記第２の関係式に基づいて前記任意の高さにおける最大距離を求め、その求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を半径とする円形形状を有する第１の等高線、前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を長径とする楕円形状を有する第２の等高線、および前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを前記作成すべき受信電力の等高線として作成する、請求項１に記載の電力推定装置。
3. 前記作成手段によって作成された受信電力の等高線に基づいて前記対象領域内の所望の場所において電波の受信電力を推定する推定手段を更に備える、請求項１または請求項２に記載の電力推定装置。
4. 前記推定手段は、前記対象領域における電波環境に応じた電波伝搬モデルを用いて、前記対象領域内の所望の場所において電波の受信電力を推定する、請求項３に記載の電力推定装置。
5. 前記処理手段は、前記ｍ個の高さに対応して前記対象領域内の高さ方向をｍ個の領域に分割し、その分割した１つの領域に含まれる高さを有する前記ｋ個のモニター情報に基づいて前記第１の処理を前記ｍ個の領域の全てについて実行する、請求項１から請求項４のいずれか１項に記載の電力推定装置。
6. 前記処理手段は、前記対象領域内の建物の高さまたは建物の平均高さよりも小さい高さ方向の長さを有するように前記対象領域内の高さ方向をｍ個の領域に分割する、請求項５に記載の電力推定装置。
7. 前記ｍ個の領域は、
   建物の高さまたは建物の平均高さ以下の位置に存在し、かつ、高さ方向に第１の長さを有する第１の領域と、
   前記建物の高さまたは建物の平均高さよりも高い位置に存在し、かつ、高さ方向に前記第１の長さよりも長い第２の長さを有する第２の領域とを含む、請求項５または請求項６に記載の電力推定装置。
8. 請求項１から請求項７のいずれか１項に記載の電力推定装置を備えた無線通信システム。
9. 電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線の作成をコンピュータに実行させるためのプログラムであって、
   抽出手段が、前記複数の端末装置から送信され、かつ、各々が前記端末装置の位置を示す位置情報と前記端末装置における前記電波の受信電力とを含む複数のモニター情報から、前記受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する第１のステップと、
   処理手段が、前記第１のステップにおいて抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、前記ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、前記受信電力の重心点から最も遠い位置に存在する端末装置と前記受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する第２のステップと、
   演算手段が、前記第２のステップにおいて求められたｍ個の受信電力の重心点およびｍ個の最大距離と、前記ｍ個の高さとに基づいて前記受信電力の重心点と前記高さとの関係を示す第１の関係式と、前記受信電力の複数の電力値に対応付けられ、かつ、各々が前記最大距離と前記高さとの関係を示す複数の第２の関係式とを演算する第３のステップと、
   作成手段が、前記第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、前記複数の第２の関係式のうち、前記作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて前記任意の高さにおける前記最大距離を求め、その求めた受信電力の重心点および最大距離を用いて前記作成すべき受信電力の等高線を作成する第４のステップとをコンピュータに実行させるためのプログラム。
10. 前記作成手段は、前記第４のステップにおいて、前記第１の関係式に基づいて前記任意の高さにおける受信電力の重心点を求めるとともに前記第２の関係式に基づいて前記任意の高さにおける最大距離を求め、前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を半径とする円形形状を有する第１の等高線、前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を長径とする楕円形状を有する第２の等高線、および前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを前記作成すべき受信電力の等高線として作成する、請求項９に記載のコンピュータに実行させるためのプログラム。
11. 推定手段が、前記第４のステップにおいて作成された受信電力の等高線に基づいて前記対象領域内の所望の場所において電波の受信電力を推定する第５のステップを更にコンピュータに実行させる、請求項９または請求項１０に記載のコンピュータに実行させるためのプログラム。
12. 前記推定手段は、前記第５のステップにおいて、前記対象領域における電波環境に応じた電波伝搬モデルを用いて、前記対象領域内の所望の場所において電波の受信電力を推定する、請求項１１に記載のコンピュータに実行させるためのプログラム。
13. 前記処理手段は、前記第２のステップにおいて、前記ｍ個の高さに対応して前記対象領域内の高さ方向をｍ個の領域に分割し、その分割した１つの領域に含まれる高さを有する前記ｋ個のモニター情報に基づいて前記第１の処理を前記ｍ個の領域の全てについて実行する、請求項９から請求項１２のいずれか１項に記載のコンピュータに実行させるためのプログラム。
14. 前記処理手段は、前記第２のステップにおいて、前記対象領域内の建物の高さまたは建物の平均高さよりも小さい高さ方向の長さを有するように前記対象領域内の高さ方向をｍ個の領域に分割する、請求項１３に記載のコンピュータに実行させるためのプログラム。
15. 前記ｍ個の領域は、
    建物の高さまたは建物の平均高さ以下の位置に存在し、かつ、高さ方向に第１の長さを有する第１の領域と、
    前記建物の高さまたは建物の平均高さよりも高い位置に存在し、かつ、高さ方向に前記第１の長さよりも長い第２の長さを有する第２の領域とを含む、請求項１３または請求項１４に記載のコンピュータに実行させるためのプログラム。
16. 複数のモニター情報を含むデータ構造であって、
    前記複数のモニター情報の各々は、端末装置の位置を示す位置情報と、前記端末装置が無線通信に共用する共用周波数と、前記端末装置によって検出された電波の受信電力とを相互に対応付けた構造を有し、
    前記複数のモニター情報に含まれる複数の受信電力のうち、前記受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報は、抽出手段によって前記複数のモニター情報から抽出され、
    前記抽出手段によって抽出されたｋ個のモニター情報は、処理手段が受信電力によって重み付けされたｋ個の端末装置の位置の平均である受信電力の重心点を求めるのに用いられ、
    前記抽出手段によって抽出されたｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力と前記受信電力の等高線における電力値以上の受信電力に対応付けられた位置情報とは、前記処理手段が、前記受信電力の重心点から最も遠い位置に存在する端末装置と前記受信電力の重心点との距離である最大距離を求めるのに用いられ、
    前記処理手段によって求められた受信電力の重心点は、演算手段が前記受信電力の重心点と高さとの関係を示す第１の関係式を求めるのに用いられ、
    前記処理手段によって求められた最大距離は、演算手段が、作成すべき受信電力の等高線の電力値に対応し、かつ、前記最大距離と前記高さとの関係を示す第２の関係式を求めるのに用いられ、
    前記演算手段によって求められた第１および第２の関係式は、作成手段が、任意の高さにおいて、任意の電力値を有する受信電力の等高線を作成するのに用いられる、データ構造。

Images