(57)【要約】【課題】受信電力の３次元的な電力分布を推定可能な電力推定装置を提供する。【解決手段】抽出手段１３は、端末装置の位置情報と電波の受信電力とを含む複数のモニター情報から受信電力の最大値を含むｋ個のモニター情報を抽出する。処理手段１４は、ｋ個のモニター情報に基づいて、受信電力の重心点Ｇと、重心点Ｇとの最大距離とを求める処理をｍ個の高さについて実行する。演算手段１５は、ｍ個の重心点Ｇ、ｍ個の最大距離及びｍ個の高さに基づいて重心点Ｇと高さとの関係を示す第１の関係式と、最大距離と高さとの関係を示す第２の関係式とを演算する。作成手段１６は、第１及び第２の関係式に基づいて、任意の高さにおける重心点Ｇ及び最大距離を求め、その求めた重心点Ｇを中心とし、最大距離を半径とした等高線を作成する。【選択図】図２
2019197957.tif 000002
【特許請求の範囲】
【請求項１】
電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線を作成する電力推定装置であって、
前記複数の端末装置から送信され、かつ、各々が前記端末装置の位置を示す位置情報と前記端末装置における前記電波の受信電力とを含む複数のモニター情報から、前記受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する抽出手段と、
前記抽出手段によって抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、前記ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、前記受信電力の重心点から最も遠い位置に存在する端末装置と前記受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する処理手段と、
前記処理手段によって求められたｍ個の受信電力の重心点およびｍ個の最大距離と、前記ｍ個の高さとに基づいて前記受信電力の重心点と前記高さとの関係を示す第１の関係式と、前記受信電力の複数の電力値に対応付けられ、かつ、各々が前記最大距離と前記高さとの関係を示す複数の第２の関係式とを演算する演算手段と、
前記第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、前記複数の第２の関係式のうち、前記作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて前記任意の高さにおける前記最大距離を求め、その求めた受信電力の重心点および最大距離を用いて前記作成すべき受信電力の等高線を作成する作成手段とを備える電力推定装置。
【請求項２】
前記作成手段は、前記第１の関係式に基づいて任意の高さにおける受信電力の重心点を求めるとともに前記第２の関係式に基づいて前記任意の高さにおける最大距離を求め、その求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を半径とする円形形状を有する第１の等高線、前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を長径とする楕円形状を有する第２の等高線、および前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを前記作成すべき受信電力の等高線として作成する、請求項１に記載の電力推定装置。
【請求項３】
前記作成手段によって作成された受信電力の等高線に基づいて前記対象領域内の所望の場所において電波の受信電力を推定する推定手段を更に備える、請求項１または請求項２に記載の電力推定装置。
【請求項４】
前記推定手段は、前記対象領域における電波環境に応じた電波伝搬モデルを用いて、前記対象領域内の所望の場所において電波の受信電力を推定する、請求項３に記載の電力推定装置。
【請求項５】
前記処理手段は、前記ｍ個の高さに対応して前記対象領域内の高さ方向をｍ個の領域に分割し、その分割した１つの領域に含まれる高さを有する前記ｋ個のモニター情報に基づいて前記第１の処理を前記ｍ個の領域の全てについて実行する、請求項１から請求項４のいずれか１項に記載の電力推定装置。
【請求項６】
前記処理手段は、前記対象領域内の建物の高さまたは建物の平均高さよりも小さい高さ方向の長さを有するように前記対象領域内の高さ方向をｍ個の領域に分割する、請求項５に記載の電力推定装置。
【請求項７】
前記ｍ個の領域は、
建物の高さまたは建物の平均高さ以下の位置に存在し、かつ、高さ方向に第１の長さを有する第１の領域と、
前記建物の高さまたは建物の平均高さよりも高い位置に存在し、かつ、高さ方向に前記第１の長さよりも長い第２の長さを有する第２の領域とを含む、請求項５または請求項６に記載の電力推定装置。
【請求項８】
請求項１から請求項７のいずれか１項に記載の電力推定装置を備えた無線通信システム。
【請求項９】
電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線の作成をコンピュータに実行させるためのプログラムであって、
抽出手段が、前記複数の端末装置から送信され、かつ、各々が前記端末装置の位置を示す位置情報と前記端末装置における前記電波の受信電力とを含む複数のモニター情報から、前記受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する第１のステップと、
処理手段が、前記第１のステップにおいて抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、前記ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、前記受信電力の重心点から最も遠い位置に存在する端末装置と前記受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する第２のステップと、
演算手段が、前記第２のステップにおいて求められたｍ個の受信電力の重心点およびｍ個の最大距離と、前記ｍ個の高さとに基づいて前記受信電力の重心点と前記高さとの関係を示す第１の関係式と、前記受信電力の複数の電力値に対応付けられ、かつ、各々が前記最大距離と前記高さとの関係を示す複数の第２の関係式とを演算する第３のステップと、
作成手段が、前記第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、前記複数の第２の関係式のうち、前記作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて前記任意の高さにおける前記最大距離を求め、その求めた受信電力の重心点および最大距離を用いて前記作成すべき受信電力の等高線を作成する第４のステップとをコンピュータに実行させるためのプログラム。
【請求項１０】
前記作成手段は、前記第４のステップにおいて、前記第１の関係式に基づいて前記任意の高さにおける受信電力の重心点を求めるとともに前記第２の関係式に基づいて前記任意の高さにおける最大距離を求め、前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を半径とする円形形状を有する第１の等高線、前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を長径とする楕円形状を有する第２の等高線、および前記求めた受信電力の重心点を中心とし、かつ、前記求めた最大距離を中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを前記作成すべき受信電力の等高線として作成する、請求項９に記載のコンピュータに実行させるためのプログラム。
【請求項１１】
推定手段が、前記第４のステップにおいて作成された受信電力の等高線に基づいて前記対象領域内の所望の場所において電波の受信電力を推定する第５のステップを更にコンピュータに実行させる、請求項９または請求項１０に記載のコンピュータに実行させるためのプログラム。
【請求項１２】
前記推定手段は、前記第５のステップにおいて、前記対象領域における電波環境に応じた電波伝搬モデルを用いて、前記対象領域内の所望の場所において電波の受信電力を推定する、請求項１１に記載のコンピュータに実行させるためのプログラム。
【請求項１３】
前記処理手段は、前記第２のステップにおいて、前記ｍ個の高さに対応して前記対象領域内の高さ方向をｍ個の領域に分割し、その分割した１つの領域に含まれる高さを有する前記ｋ個のモニター情報に基づいて前記第１の処理を前記ｍ個の領域の全てについて実行する、請求項９から請求項１２のいずれか１項に記載のコンピュータに実行させるためのプログラム。
【請求項１４】
前記処理手段は、前記第２のステップにおいて、前記対象領域内の建物の高さまたは建物の平均高さよりも小さい高さ方向の長さを有するように前記対象領域内の高さ方向をｍ個の領域に分割する、請求項１３に記載のコンピュータに実行させるためのプログラム。
【請求項１５】
前記ｍ個の領域は、
建物の高さまたは建物の平均高さ以下の位置に存在し、かつ、高さ方向に第１の長さを有する第１の領域と、
前記建物の高さまたは建物の平均高さよりも高い位置に存在し、かつ、高さ方向に前記第１の長さよりも長い第２の長さを有する第２の領域とを含む、請求項１３または請求項１４に記載のコンピュータに実行させるためのプログラム。
【請求項１６】
複数のモニター情報を含むデータ構造であって、
前記複数のモニター情報の各々は、端末装置の位置を示す位置情報と、前記端末装置が無線通信に共用する共用周波数と、前記端末装置によって検出された電波の受信電力とを相互に対応付けた構造を有し、
前記複数のモニター情報に含まれる複数の受信電力のうち、前記受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報は、抽出手段によって前記複数のモニター情報から抽出され、
前記抽出手段によって抽出されたｋ個のモニター情報は、処理手段が受信電力によって重み付けされたｋ個の端末装置の位置の平均である受信電力の重心点を求めるのに用いられ、
前記抽出手段によって抽出されたｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力と前記受信電力の等高線における電力値以上の受信電力に対応付けられた位置情報とは、前記処理手段が、前記受信電力の重心点から最も遠い位置に存在する端末装置と前記受信電力の重心点との距離である最大距離を求めるのに用いられ、
前記処理手段によって求められた受信電力の重心点は、演算手段が前記受信電力の重心点と高さとの関係を示す第１の関係式を求めるのに用いられ、
前記処理手段によって求められた最大距離は、演算手段が、作成すべき受信電力の等高線の電力値に対応し、かつ、前記最大距離と前記高さとの関係を示す第２の関係式を求めるのに用いられ、
前記演算手段によって求められた第１および第２の関係式は、作成手段が、任意の高さにおいて、任意の電力値を有する受信電力の等高線を作成するのに用いられる、データ構造。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、電力推定装置、それを備えた無線通信システム、コンピュータに実行させるためのプログラムおよびデータ構造に関する。
【背景技術】
【０００２】
移動通信に適した６ＧＨｚ以下の周波数帯域に対して、周波数を共用することによって、第５世代移動通信システムに必要な周波数帯域を新たに創出する研究がなされている。
【０００３】
その一環として、既に幅広いエリアに分布している移動端末をモニター端末とし、取得したモニター端末の位置および電波の受信電力の情報から、場所ごとの既存システムにおける電波の受信電力を求め、これらを解析することで、周波数共用の条件を決定する移動分散モニタリングに着目して研究が進められている。
【０００４】
移動分散モニタリングでは、モニター端末の分布にバラツキがあるため、モニター端末が存在しない場所の受信電力は、周辺の受信電力などから推定する必要がある。
【０００５】
受信電力を推定する方式として、モニター端末が取得したモニター情報（モニター端末の位置および電波の受信電力の情報）から受信電力の等高線を作成し、その作成した受信電力の等高線から場所ごとの受信電力を推定する方法が考えられる。
【０００６】
従来、受信電力の等高線を作成する方法として、非特許文献１に記載の方法が知られている。非特許文献１に記載された受信電力の等高線の作成方法は、２つのモニター端末における２つの受信電力に基づいて受信電力のベクトルを作成し、その作成したベクトルと直交する勾配（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）、離島へ電波を中継する中継所および端末装置等からなる。
【００６８】
図２は、図１に示す電力推定装置１の概略図である。図２を参照して、電力推定装置１は、受信手段１１と、記録手段１２と、抽出手段１３と、処理手段１４と、演算手段１５と、作成手段１６と、推定手段１７とを備える。
【００６９】
受信手段１１は、複数のモニター情報ＭＮＴを基地局３から受信し、その受信した複数のモニター情報ＭＮＴを記録手段１２に記録する。
【００７０】
記録手段１２は、受信手段１１から受けた複数のモニター情報ＭＮＴを記録する。より具体的には、記録手段１２は、複数のモニター情報ＭＮＴを［ｔ１：（ｘ１，ｙ１，ｚ１），ＲＳＳＩ１，ｆｃｏｍ＿１］，［ｔ２：（ｘ２，ｙ２，ｚ２），ＲＳＳＩ２，ｆｃｏｍ＿１］，［ｔ３：（ｘ３，ｙ３，ｚ３），ＲＳＳＩ３，ｆｃｏｍ＿１］，・・・,［ｔＴ：（ｘＭ，ｙＭ，ｚＭ），ＲＳＳＩＭ，ｆｃｏｍ＿１］；［ｔ１：（ｘ１，ｙ１，ｚ１），ＲＳＳＩ１，ｆｃｏｍ＿２］，［ｔ２：（ｘ２，ｙ２，ｚ２），ＲＳＳＩ２，ｆｃｏｍ＿２］，［ｔ３：（ｘ３，ｙ３，ｚ３），ＲＳＳＩ３，ｆｃｏｍ＿２］，・・・,［ｔＴ：（ｘＭ，ｙＭ，ｚＭ），ＲＳＳＩＭ，ｆｃｏｍ＿２］；・・・；［ｔ１：（ｘ１，ｙ１，ｚ１），ＲＳＳＩ１，ｆｃｏｍ＿Ｆ］，［ｔ２：（ｘ２，ｙ２，ｚ２），ＲＳＳＩ２，ｆｃｏｍ＿Ｆ］，［ｔ３：（ｘ３，ｙ３，ｚ３），ＲＳＳＩ３，ｆｃｏｍ＿Ｆ］，・・・,［ｔＴ：（ｘＭ，ｙＭ，ｚＭ），ＲＳＳＩＭ，ｆｃｏｍ＿Ｆ］のように記録する。ここで、Ｔは、時刻情報ｔｉの検出回数を示し、Ｍは、端末装置２の個数を示し、Ｆは、共用周波数ｆｃｏｍの個数を示す。そして、Ｔは、２以上の整数であり、Ｍは、後述するｋよりも大きい整数であり、Ｆは、１以上の整数である。
【００７１】
また、記録手段１２は、移動端末であるが停止している端末装置２または静止端末である端末装置２に対しては、位置情報（ｘｉ，ｙｉ，ｚｉ）が同じで、時刻情報ｔｉ、受信電力ＲＳＳＩｉおよび周波数情報ｆｃｏｍ＿ｊ（ｊは、１≦ｊ≦Ｆを満たす整数）が異なるようにモニター情報ＭＮＴを記録する。
【００７２】
更に、記録手段１２は、時刻情報ｔｉの検出間隔が位置情報（ｘｉ，ｙｉ，ｚｉ）、受信電力ＲＳＳＩｉおよび周波数情報ｆｃｏｍ＿ｊの検出間隔よりも長い場合、複数のモニター情報ＭＮＴを［ｔ１：（ｘ１，ｙ１，ｚ１），ＲＳＳＩ１，ｆｃｏｍ＿１］，［ｔ１：（ｘ２，ｙ２，ｚ２），ＲＳＳＩ２，ｆｃｏｍ＿１］，［ｔ１：（ｘ３，ｙ３，ｚ３），ＲＳＳＩ３，ｆｃｏｍ＿１］，・・・,［ｔ１：（ｘＭ，ｙＭ，ｚＭ），ＲＳＳＩＭ，ｆｃｏｍ＿１］；［ｔ２：（ｘ１，ｙ１，ｚ１），ＲＳＳＩ１，ｆｃｏｍ＿２］，［ｔ２：（ｘ２，ｙ２，ｚ２），ＲＳＳＩ２，ｆｃｏｍ＿２］，［ｔ２：（ｘ３，ｙ３，ｚ３），ＲＳＳＩ３，ｆｃｏｍ＿２］，・・・,［ｔ２：（ｘＭ，ｙＭ，ｚＭ），ＲＳＳＩＭ，ｆｃｏｍ＿２］；・・・；［ｔＴ：（ｘ１，ｙ１，ｚ１），ＲＳＳＩ１，ｆｃｏｍ＿Ｆ］，［ｔＴ：（ｘ２，ｙ２，ｚ２），ＲＳＳＩ２，ｆｃｏｍ＿Ｆ］，［ｔＴ：（ｘ３，ｙ３，ｚ３），ＲＳＳＩ３，ｆｃｏｍ＿Ｆ］，・・・,［ｔＴ：（ｘＭ，ｙＭ，ｚＭ），ＲＳＳＩＭ，ｆｃｏｍ＿Ｆ］のように記録する。即ち、記録手段１２は、同じ時刻情報ｔｉにおいて、位置情報（ｘｉ，ｙｉ，ｚｉ）および受信電力ＲＳＳＩｉが変化するように複数のモニター情報ＭＮＴを記録する。
【００７３】
抽出手段１３は、後述する方法によって、記録手段１２に記録された複数のモニター情報ＭＮＴに基づいて、複数のモニター情報ＭＮＴの分布状態を作成し、その作成した複数のモニター情報ＭＮＴの分布状態に基づいて受信電力の等高線を作成する処理を行うための処理範囲を決定し、その決定した処理範囲に含まれるｋ（ｋは２以上の整数）個の端末装置２のｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを抽出する。そして、抽出手段１３は、その抽出したｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを処理手段１４へ出力する。抽出手段１３は、作成すべき受信電力の等高線ＣＴＲの電力値Ｐ以上の受信電力ＲＳＳＩｉがなくなるまで、後述する方法によって、複数のモニター情報ＭＮＴからｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを抽出し、その抽出したｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを処理手段１４へ出力する処理を繰り返し行う。
【００７４】
処理手段１４は、抽出手段１３からｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを受け、その受けたｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに基づいて、後述する方法によって、受信電力の重心点ＲＳＳＩ＿Ｇを算出するとともに、ｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇから最も遠い位置に存在する端末装置２と受信電力の重心点ＲＳＳＩ＿Ｇとの距離を最大距離ＬＭＡＸとして求める。そして、処理手段１４は、受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離ＬＭＡＸを求める処理をｍ（ｍは、２以上の整数）個の高さｚ１〜ｚｍについて実行する。そうすると、処理手段１４は、ｍ個の受信電力の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離ＬＭＡＸ＿１〜ＬＭＡＸ＿ｍおよびｍ個の高さｚ１〜ｚｍを演算手段１５へ出力する。
【００７５】
演算手段１５は、ｍ個の受信電力の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離ＬＭＡＸ＿１〜ＬＭＡＸ＿ｍおよびｍ個の高さｚ１〜ｚｍを処理手段１４から受ける。そして、演算手段１５は、ｍ個の受信電力の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離ＬＭＡＸ＿１〜ＬＭＡＸ＿ｍおよびｍ個の高さｚ１〜ｚｍに基づいて、後述する方法によって、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧと高さｚとの関係を示す関係式ＥＱ１、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙＧと高さｚとの関係を示す関係式ＥＱ２および最大距離ＬＭＡＸと高さｚとの関係を示す関係式ＥＱ３を演算する。そうすると、演算手段１５は、関係式ＥＱ１，ＥＱ２，ＥＱ３を作成手段１６へ出力する。
【００７６】
作成手段１６は、関係式ＥＱ１，ＥＱ２，ＥＱ３を演算手段１５から受ける。そして、作成手段１６は、関係式ＥＱ１，ＥＱ２，ＥＱ３に基づいて、後述する方法によって、受信電力の３次元的な電力分布を示す等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲを推定手段１７へ出力する。
【００７７】
推定手段１７は、受信電力の等高線ＣＴＲを作成手段１６から受け、その受けた受信電力の等高線ＣＴＲに基づいて、後述する方法によって、対象領域ＲＥＧの任意の位置における電波の受信電力を推定する。そして、推定手段１７は、その推定した電波の受信電力に基づいて、複数の端末装置２が共用周波数ｆｃｏｍを用いて無線通信を行うときの通信条件（共用周波数ｆｃｏｍ、場所および送信電力等）を決定する。
【００７８】
図３は、対象領域ＲＥＧを分割した分割領域の概念図である。図３を参照して、対象領域ＲＥＧは、水平方向において分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑ（ｑは２以上の整数）に分割される。
【００７９】
分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々は、地面に平行な水平方向において長さＤを有し、地面に垂直な方向において高さＨを有する。例えば、長さＤは、ドローンの飛行領域における水平方向の長さに設定され、高さＨは、ドローンの飛行領域における垂直方向の長さに設定される。
【００８０】
なお、分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑは、相互に、同じ長さＤおよび高さＨを有していてもよいし、相互に、異なる長さＤおよび高さＨを有していてもよい。
【００８１】
図４は、図３に示す１つの分割領域ＲＥＧ＿ＤＶｑを高さ方向に分割する方法を説明するための図である。図４を参照して、抽出手段１３は、建物２０の高さＨＢよりも小さい高さ方向の長さｈを予め保持している。長さｈは、例えば、３〜１０ｍである。そして、抽出手段１３は、対象領域ＲＥＧの高さ方向を長さｈで分割して複数の領域ＳＰＣ１〜ＳＰＣｍを生成する。なお、高さＨＢは、建物の平均高さであってもよい。
【００８２】
抽出手段１３は、複数の領域ＳＰＣ１〜ＳＰＣｍを生成すると、モニター情報ＭＮＴに含まれる位置情報のｚ座標に基づいて複数の領域ＳＰＣ１〜ＳＰＣｍの各々に含まれるモニター情報を抽出する。この場合、抽出手段１３は、領域ＳＰＣ１に含まれるモニター情報ＭＮＴ＿ＳＰＣ１、領域ＳＰＣ２に含まれるモニター情報ＭＮＴ＿ＳＰＣ２、・・・、および領域ＳＰＣｍに含まれるモニター情報ＭＮＴ＿ＳＰＣｍを抽出する。モニター情報ＭＮＴ＿ＳＰＣ１〜ＭＮＴ＿ＳＰＣｍの各々は、複数のモニター情報からなる。
【００８３】
図５は、受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離ＬＭＡＸを求める方法を説明するための図である。図５において、白四角は、電力値がＰ（＝受信電力の等高線ＣＴＲが有する電力値）以上である測定点を示し、黒四角は、電力値がＰ未満である測定点を示し、白丸は、受信電力の重心点ＲＳＳＩ＿Ｇを示す。また、測定点ｓ１〜ｓ３は、フェージングの影響によって受信電力が電力値Ｐ未満になった測定点を示す。そして、白四角および黒四角の各々は、位置情報（ｘｉ，ｙｉ，ｚｉ）および受信電力ＲＳＳＩｉからなる。
【００８４】
図５の（ａ）を参照して、電力推定装置１の抽出手段１３は、１つの領域ＳＰＣ１に含まれる複数のモニター情報ＭＮＴに基づいて、複数のモニター情報ＭＮＴの分布状態（図５の（ａ）に示す分布状態）を取得する。
【００８５】
そして、抽出手段１３は、作成すべき受信電力の等高線ＣＴＲの電力値Ｐを決定する。抽出手段１３は、例えば、−６０ｄＢｍ，−７０ｄＢｍ，−８０ｄＢｍ等を電力値Ｐとして決定する。
【００８６】
抽出手段１３は、電力値Ｐと処理範囲の半径ｒとの関係を予め保持している。例えば、抽出手段１３は、−６０ｄＢｍの電力値Ｐに対して処理範囲の半径ｒ＝５００ｍ、および−７０ｄＢｍの電力値Ｐに対して処理範囲の半径１５００ｍ等の関係を予め保持している。
【００８７】
抽出手段１３は、作成すべき受信電力の等高線ＣＴＲの電力値Ｐを決定すると、複数のモニター情報ＭＮＴに含まれる複数の受信電力ＲＳＳＩｉのうち、電力値Ｐ以上である受信電力の最大値ＲＳＳＩ＿ＭＡＸを抽出する。ここで、受信電力の最大値ＲＳＳＩ＿ＭＡＸを抽出するのは、受信電力の最大値ＲＳＳＩ＿ＭＡＸを有する端末装置２が最も波源Ｓに近いと考えられるからである。
【００８８】
そして、抽出手段１３は、−６０ｄＢｍの電力値Ｐを有する受信電力の等高線ＣＴＲを作成する場合、受信電力の最大値ＲＳＳＩ＿ＭＡＸを抽出すると、その抽出した受信電力の最大値ＲＳＳＩ＿ＭＡＸを有する測定点を中心として半径ｒ＝５００ｍの範囲を処理範囲ＲＥＧ＿ＰＲＳとして決定する。
【００８９】
そうすると、抽出手段１３は、処理範囲ＲＥＧ＿ＰＲＳ内のｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを抽出し、その抽出したｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを処理手段１４へ出力する。
【００９０】
図５の（ｂ）を参照して、処理手段１４は、ｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを抽出手段１３から受け、その受けたｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに含まれるｋ個の受信電力ＲＳＳＩ１〜ＲＳＳＩｋに基づいて、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘＧ，ｙＧ）を次式によって算出する。
【００９１】
【数１】
2019197957.tif 000003
【００９２】
式（１）において、ｗｉは、ｉ番目の端末装置２のウェイトを表し、ｉ番目の端末装置２における受信電力ＲＳＳＩｉからなる。
【００９３】
即ち、処理手段１４は、ｋ個の受信電力ＲＳＳＩｉ（ｉ＝１〜ｋ）によって重み付けされたｋ個の端末装置２の位置の平均を演算することによって受信電力の重心点ＲＳＳＩ＿Ｇ（ｘＧ，ｙＧ）を算出する。
【００９４】
次に、処理手段１４は、ｋ個の受信電力ＲＳＳＩｉ（ｉ＝１〜ｋ）のうち、電力値Ｐ以上の受信電力ＲＳＳＩｉを有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘＧ，ｙＧ）から最も遠い位置に存在する端末装置２−１と受信電力の重心点ＲＳＳＩ＿Ｇ（ｘＧ，ｙＧ）との間の距離を最大距離ＬＭＡＸとして求める。
【００９５】
このように、処理手段１４は、電力値Ｐ未満の受信電力（黒四角）を除外し、電力値Ｐ以上の受信電力（白四角）のみに基づいて最大距離ＬＭＡＸを求める。黒四角によって示される測定点のうち、測定点ｓ１〜ｓ３は、フェージングの影響によって受信電力が電力値Ｐ未満になった測定点である。
【００９６】
なお、処理手段１４は、上述した方法によって、複数の電力値Ｐ（＝Ｐ１〜ＰＮ（Ｎは、２以上の整数））について最大距離ＬＭＡＸ＿１〜ＬＭＡＸ＿Ｎを求める。なお、複数の最大距離ＬＭＡＸ＿１〜ＬＭＡＸ＿Ｎは、それぞれ、作成すべき受信電力の等高線ＣＴＲの電力値Ｐ１〜ＰＮに対応する最大距離である。
【００９７】
従って、処理手段１４による受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離ＬＭＡＸを求める方法によれば、フェージングの影響を抑制して受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離ＬＭＡＸを求めることができる。
【００９８】
上記においては、ｋは、２以上の整数であると説明したが、これは、次の理由による。ｋ＝２である場合、２個の受信電力ＲＳＳＩｉのうち、１個は、受信電力の最大値ＲＳＳＩ＿ＭＡＸであり、もう１個は、電力値Ｐ以上である受信電力か電力値Ｐ未満である受信電力かが不明である受信電力である。そして、この２個の受信電力ＲＳＳＩｉを用いて、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘＧ，ｙＧ）が式（１）によって算出される。
【００９９】
その結果、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘＧ，ｙＧ）は、受信電力の最大値ＲＳＳＩ＿ＭＡＸを有する測定点からずれる。
【０１００】
そうすると、受信電力の最大値ＲＳＳＩ＿ＭＡＸは、電力値Ｐ以上であるので、受信電力の重心点ＲＳＳＩ＿Ｇ（ｘＧ，ｙＧ）と受信電力の最大値ＲＳＳＩ＿ＭＡＸを有する測定点との距離を最大距離ＬＭＡＸとして求めることができる。従って、ｋは、２以上であればよい。
【０１０１】
なお、電力値Ｐ以上の受信電力ＲＳＳＩｉが少なくとも１個有れば、少なくとも１個の受信電力ＲＳＳＩｉを有する端末装置２の位置が受信電力の重心点ＲＳＳＩ＿Ｇからずれているので、最大距離ＬＭＡＸを決定できる。従って、受信電力の等高線ＣＴＲを作成するために必要な電力値Ｐ以上の受信電力の個数は、１個以上である。
【０１０２】
抽出手段１３は、領域ＳＰＣ２〜ＳＰＣｍの各々においても、上述した方法によってｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを抽出し、処理手段１４は、抽出手段１３によって抽出されたｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに基づいて、上述した方法によって受信電力の重心点ＲＳＳＩ＿Ｇおよび最大距離ＬＭＡＸ（＝ＬＭＡＸ＿１〜ＬＭＡＸ＿Ｎ）を求める。
【０１０３】
抽出手段１３および処理手段１４は、上述した処理を複数の領域ＳＰＣ１〜ＳＰＣｍの全てについて実行し、複数の領域ＳＰＣ１〜ＳＰＣｍの各々について受信電力の重心点ＲＳＳＩ＿Ｇおよび複数の最大距離ＬＭＡＸ＿１〜ＬＭＡＸ＿Ｎを求める。その結果、処理手段１４は、領域ＳＰＣ１において、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ１および複数の最大距離ＬＭＡＸ＿１＿ＳＰＣ１〜ＬＭＡＸ＿Ｎ＿ＳＰＣ１を求め、領域ＳＰＣ２において、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ２および複数の最大距離ＬＭＡＸ＿１＿ＳＰＣ２〜ＬＭＡＸ＿Ｎ＿ＳＰＣ２を求め、・・・、領域ＳＰＣｍにおいて、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣｍおよび複数の最大距離ＬＭＡＸ＿１＿ＳＰＣｍ〜ＬＭＡＸ＿Ｎ＿ＳＰＣｍを求める。
【０１０４】
そうすると、処理手段１４は、領域ＳＰＣ１、領域ＳＰＣ１におけるｋ個のモニター情報ＭＮＴ１＿ＳＰＣ１〜ＭＮＴｋ＿ＳＰＣ１、受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ１および複数の最大距離ＬＭＡＸ＿１＿ＳＰＣ１〜ＬＭＡＸ＿Ｎ＿ＳＰＣ１を相互に対応付けた情報ＩＮＦ１＝［ＳＰＣ１／モニター情報ＭＮＴ１＿ＳＰＣ１〜ＭＮＴｋ＿ＳＰＣ１／ＲＳＳＩ＿Ｇ＿ＳＰＣ１／ＬＭＡＸ＿１＿ＳＰＣ１〜ＬＭＡＸ＿Ｎ＿ＳＰＣ１］を生成する。また、処理手段１４は、同様にして、領域ＳＰＣ２に対応する情報ＩＮＦ２＝［ＳＰＣ２／モニター情報ＭＮＴ１＿ＳＰＣ２〜ＭＮＴｋ＿ＳＰＣ２／ＲＳＳＩ＿Ｇ＿ＳＰＣ２／ＬＭＡＸ＿１＿ＳＰＣ２〜ＬＭＡＸ＿Ｎ＿ＳＰＣ２］、・・・、領域ＳＰＣｍに対応する情報ＩＮＦｍ＝［ＳＰＣｍ／モニター情報ＭＮＴ１＿ＳＰＣｍ〜ＭＮＴｋ＿ＳＰＣｍ／ＲＳＳＩ＿Ｇ＿ＳＰＣｍ／ＬＭＡＸ＿１＿ＳＰＣｍ〜ＬＭＡＸ＿Ｎ＿ＳＰＣｍ］を生成する。そして、処理手段１４は、情報ＩＮＦ１〜ＩＮＦ４を演算手段１５へ出力する。
【０１０５】
図６は、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。図６においては、高さ方向の領域ＳＰＣが４個（ｍ＝４）である場合について、重心点ＲＳＳＩ＿Ｇのｘ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を説明する。
【０１０６】
図６において、縦軸は、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標を表し、横軸は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚを表す。各領域ＳＰＣ１〜ＳＰＣ４の高さは、例えば、各領域ＳＰＣ１〜ＳＰＣ４において抽出されたｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに含まれるｋ個の位置情報のｋ個のｚ座標の平均値として求められる。また、各領域ＳＰＣ１〜ＳＰＣ４の高さは、各領域ＳＰＣ１〜ＳＰＣ４において抽出されたｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに含まれるｋ個の位置情報のｋ個のｚ座標をｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに含まれるｋ個の受信電力ＲＳＳＩ１〜ＲＳＳＩｋによって重み付けした重み付け平均（ｚＧ＝（ｗ１ｚ１＋ｗ２ｚ２＋・・・＋ｗｋｚｋ）／（ｗ１＋ｗ２＋・・・＋ｗｋ）、ｗｉ＝ＲＳＳＩｉ）によって求められる。
【０１０７】
演算手段１５は、情報ＩＮＦ１〜ＩＮＦ４を処理手段１４から受ける。演算手段１５は、情報ＩＮＦ１＝［ＳＰＣ１／モニター情報ＭＮＴ１＿ＳＰＣ１〜ＭＮＴｋ＿ＳＰＣ１／ＲＳＳＩ＿Ｇ＿ＳＰＣ１／ＬＭＡＸ＿１＿ＳＰＣ１〜ＬＭＡＸ＿ｎ＿ＳＰＣ１］からｋ個のモニター情報ＭＮＴ１＿ＳＰＣ１〜ＭＮＴｋ＿ＳＰＣ１を検出し、その検出したｋ個のモニター情報ＭＮＴ１＿ＳＰＣ１〜ＭＮＴｋ＿ＳＰＣ１に含まれるｋ個の位置情報のｋ個のｚ座標ｚ１〜ｚｋ（またはｋ個のｚ座標ｚ１〜ｚｋおよびｋ個の受信電力ＲＳＳＩ１〜ＲＳＳＩｋ）を検出する。そして、演算手段１５は、ｋ個のｚ座標ｚ１〜ｚｋ（またはｋ個のｚ座標ｚ１〜ｚｋおよびｋ個の受信電力ＲＳＳＩ１〜ＲＳＳＩｋ）に基づいて、上述した方法によって領域ＳＰＣ１の高さｚＳＰＣ１を演算する。同様にして、演算手段１５は、情報ＩＮＦ２〜ＩＮＦ４に基づいて、それぞれ、領域ＳＰＣ２〜ＳＰＣ４の高さｚＳＰＣ２〜ｚＳＰＣ４を演算する。
【０１０８】
そうすると、演算手段１５は、情報ＩＮＦ１〜ＩＮＦ４からそれぞれ受信電力の重心点ＲＳＳＩ＿Ｇ＿ＳＰＣ１〜ＲＳＳＩ＿Ｇ＿ＳＰＣ４のｘ座標ｘＧ＿ＳＰＣ１〜ｘＧ＿ＳＰＣ４およびｙ座標ｙＧ＿ＳＰＣ１〜ｙＧ＿ＳＰＣ４を検出する。
【０１０９】
そして、演算手段１５は、ｘ座標ｘＧ＿ＳＰＣ１〜ｘＧ＿ＳＰＣ４および高さｚＳＰＣ１〜ｚＳＰＣ４に基づいて、点（ｘＧ＿ＳＰＣ１，ｚＳＰＣ１），点（ｘＧ＿ＳＰＣ２，ｚＳＰＣ２），点（ｘＧ＿ＳＰＣ３，ｚＳＰＣ３），点（ｘＧ＿ＳＰＣ４，ｚＳＰＣ４）をプロットする（図６の黒丸参照）。
【０１１０】
その後、演算手段１５は、プロットした点に基づいて受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧ＿ＳＰＣ１〜ｘＧ＿ＳＰＣ４と各領域ＳＰＣ１〜ＳＰＣ４の高さｚＳＰＣ１〜ｚＳＰＣ４との関係を示す関係式ＥＱ１（曲線ｋ１参照）を求める。従って、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘ（ｚ）は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚの関数である。
【０１１１】
図７は、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標と各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。
【０１１２】
図７において、縦軸は、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標を表し、横軸は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚを表す。各領域ＳＰＣ１〜ＳＰＣ４の高さは、図６において説明した方法によって求められる。
【０１１３】
演算手段１５は、ｙ座標ｙＧ＿ＳＰＣ１〜ｙＧ＿ＳＰＣ４および高さｚＳＰＣ１〜ｚＳＰＣ４に基づいて、点（ｙＧ＿ＳＰＣ１，ｚＳＰＣ１），点（ｙＧ＿ＳＰＣ２，ｚＳＰＣ２），点（ｙＧ＿ＳＰＣ３，ｚＳＰＣ３），点（ｙＧ＿ＳＰＣ４，ｚＳＰＣ４）をプロットする（図７の黒丸参照）。
【０１１４】
その後、演算手段１５は、プロットした点に基づいて受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙＧ＿ＳＰＣ１〜ｙＧ＿ＳＰＣ４と各領域ＳＰＣ１〜ＳＰＣ４の高さｚＳＰＣ１〜ｚＳＰＣ４との関係を示す関係式ＥＱ２（曲線ｋ２参照）を求める。従って、受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙ（ｚ）は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚの関数である。
【０１１５】
図８は、最大距離ＬＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す図である。図８において、縦軸は、最大距離ＬＭＡＸを表し、横軸は、各領域ＳＰＣ１〜ＳＰＣ４の高さｚを表す。各領域ＳＰＣ１〜ＳＰＣ４の高さは、図６において説明した方法によって求められる。
【０１１６】
演算手段１５は、情報ＩＮＦ１〜ＩＮＦ４から、それぞれ、複数の最大距離ＬＭＡＸ＿１＿ＳＰＣ１〜ＬＭＡＸ＿ｎ＿ＳＰＣ１、複数の最大距離ＬＭＡＸ＿１＿ＳＰＣ２〜ＬＭＡＸ＿ｎ＿ＳＰＣ２、複数の最大距離ＬＭＡＸ＿１＿ＳＰＣ３〜ＬＭＡＸ＿ｎ＿ＳＰＣ３、および複数の最大距離ＬＭＡＸ＿１＿ＳＰＣ４〜ＬＭＡＸ＿ｎ＿ＳＰＣ４を検出する。
【０１１７】
そして、演算手段１５は、同じ電力値Ｐ１［ｄＢｍ］に対応する最大距離ＬＭＡＸ＿１＿ＳＰＣ１，ＬＭＡＸ＿１＿ＳＰＣ２，ＬＭＡＸ＿１＿ＳＰＣ３，ＬＭＡＸ＿１＿ＳＰＣ４を検出する。
【０１１８】
そうすると、演算手段１５は、点（ＬＭＡＸ＿１＿ＳＰＣ１，ｚＳＰＣ１）点（ＬＭＡＸ＿１＿ＳＰＣ２，ｚＳＰＣ２），点（ＬＭＡＸ＿１＿ＳＰＣ３，ｚＳＰＣ３），点（ＬＭＡＸ＿１＿ＳＰＣ４，ｚＳＰＣ４）をプロットする（図８の黒丸参照）。
【０１１９】
その後、演算手段１５は、プロットした点に基づいて最大距離ＬＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さｚとの関係を示す関係式ＥＱ３−１（直線ｋ３参照）を求める。
【０１２０】
演算手段１５は、同様にして、点（ＬＭＡＸ＿２＿ＳＰＣ１，ｚＳＰＣ１）点（ＬＭＡＸ＿２＿ＳＰＣ２，ｚＳＰＣ２），点（ＬＭＡＸ＿２＿ＳＰＣ３，ｚＳＰＣ３），点（ＬＭＡＸ＿２＿ＳＰＣ４，ｚＳＰＣ４）をプロットして、電力値Ｐ２について、最大距離ＬＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す関係式ＥＱ３−２を求める。
【０１２１】
以下、同様にして、演算手段１５は、点（ＬＭＡＸ＿Ｎ＿ＳＰＣ１，ｚＳＰＣ１）点（ＬＭＡＸ＿Ｎ＿ＳＰＣ２，ｚＳＰＣ２），点（ＬＭＡＸ＿Ｎ＿ＳＰＣ３，ｚＳＰＣ３），点（ＬＭＡＸ＿Ｎ＿ＳＰＣ４，ｚＳＰＣ４）をプロットして、電力値ＰＮについて、最大距離ＬＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す関係式ＥＱ３−Ｎを求める。
【０１２２】
従って、図８に示す縦軸の最大距離ＬＭＡＸ（Ｐ，ｚ）は、作成すべき等高線ＣＴＲの電力値Ｐと各領域ＳＰＣ１〜ＳＰＣ４の高さ（ｚ座標）との関数である。つまり、演算手段１５は、作成すべき等高線ＣＴＲの電力値Ｐ１〜ＰＮの個数Ｎに等しい個数Ｎの最大距離ＬＭＡＸと各領域ＳＰＣ１〜ＳＰＣ４の高さとの関係を示す関係式ＥＱ３＿１〜ＥＱ３＿Ｎを求める。なお、関係式ＥＱ３＿１〜ＥＱ３＿Ｎは、最大距離ＬＭＡＸと高さｚとの関係を示す関係式ＥＱ３を構成する。
【０１２３】
そして、演算手段１５は、関係式ＥＱ１、関係式ＥＱ２および関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎ）を作成手段１６へ出力する。
【０１２４】
図９は、任意の高さにおける受信電力の等高線を求める方法を説明するための図である。
【０１２５】
図９を参照して、任意の高さｚにおける受信電力の重心点ＲＳＳＩ＿Ｇを決定する関係式は、関係式ＥＱ１（＝ｘ（ｚ））と関係式ＥＱ２（＝ｙ（ｚ））とからなる。また、最大距離ＬＭＡＸを決定する関係式は、関係式ＥＱ３（＝ＬＭＡＸ（Ｐ，ｚ））からなる。
【０１２６】
関係式ＥＱ３（＝ＬＭＡＸ（Ｐ，ｚ））は、作成すべき等高線ＣＴＲの電力値Ｐ１〜ＰＮにそれぞれ対応付けられた関係式ＥＱ３＿１（＝ＬＭＡＸ＿１（ｚ））〜ＥＱ３＿Ｎ（＝ＬＭＡＸ＿Ｎ（ｚ））からなる。
【０１２７】
作成手段１６は、関係式ＥＱ１、関係式ＥＱ２および関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎ）を演算手段１５から受ける。また、作成手段１６は、作成すべき等高線ＣＴＲの高さｚおよび電力値Ｐ（Ｐ１〜ＰＮの少なくとも１つ）を外部から受ける。
【０１２８】
作成手段１６は、作成すべき等高線ＣＴＲの高さｚａと、作成すべき等高線ＣＴＲの電力値Ｐｎとを受けると、関係式ＥＱ１（＝ｘ（ｚ））に基づいて高さｚａに対応する受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧ（ｚａ）を求めるとともに関係式ＥＱ２（＝ｙ（ｚ））に基づいて高さｚａに対応する受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙＧ（ｚａ）を求める。
【０１２９】
次に、作成手段１６は、関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿ｎ）に基づいて電力値Ｐｎに対応する関係式ＥＱ３＿ｎ（ＬＭＡＸ＿ｎ（ｚ））を検出し、その検出した関係式ＥＱ３＿ｎ（ＬＭＡＸ＿ｎ（ｚ））に基づいて高さｚａに対応する最大距離ＬＭＡＸ＿ｎ（ｚａ）を求める。
【０１３０】
そうすると、作成手段１６は、受信電力の重心点（ｘＧ（ｚａ），ｙＧ（ｚａ））を中心とし、最大距離ＬＭＡＸ＿ｎ（ｚａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲ１を作成する。
【０１３１】
また、作成手段１６は、作成すべき等高線ＣＴＲの高さｚａと、作成すべき等高線ＣＴＲの電力値Ｐｎ，Ｐｎ＋２とを受けると、関係式ＥＱ１，ＥＱ２に基づいて、上述した方法によって、それぞれ、受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧ（ｚａ）およびｙ座標ｙＧ（ｚａ）を求める。
【０１３２】
次に、作成手段１６は、電力値Ｐｎに基づいて、上述した方法によって最大距離ＬＭＡＸ＿ｎ（ｚａ）を求める。そして、作成手段１６は、関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿ｎ）に基づいて電力値Ｐｎ＋２に対応する関係式ＥＱ３＿ｎ＋２（ＬＭＡＸ＿ｎ＋２（ｚ））を検出し、その検出した関係式ＥＱ３＿ｎ＋２（ＬＭＡＸ＿ｎ＋２（ｚ））に基づいて高さｚａに対応する最大距離ＬＭＡＸ＿ｎ＋２（ｚａ）を求める。
【０１３３】
そうすると、作成手段１６は、受信電力の重心点（ｘＧ（ｚａ），ｙＧ（ｚａ））を中心とし、最大距離ＬＭＡＸ＿ｎ（ｚａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲ１を作成するとともに受信電力の重心点（ｘＧ（ｚａ），ｙＧ（ｚａ））を中心とし、最大距離ＬＭＡＸ＿ｎ＋２（ｚａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲ２を作成する。
【０１３４】
作成手段１６は、作成すべき等高線ＣＴＲの高さｚａと、作成すべき等高線ＣＴＲの電力値Ｐ１〜Ｐｎとを受けた場合、上述した方法によって、高さｚａにおいて電力値Ｐ１〜Ｐｎをそれぞれ有する受信電力のｎ個の等高線ＣＴＲ１〜ＣＴＲｎを作成する。
【０１３５】
また、作成手段１６は、作成すべき等高線ＣＴＲの高さとして高さｚａ以外の高さｚｊと、作成すべき等高線ＣＴＲの電力値Ｐｎ（Ｐ１〜ＰＮの少なくとも１つ）とを受けた場合、上述した方法によって、高さｚｊにおいて電力値Ｐｎ（Ｐ１〜ＰＮの少なくとも１つ）を有する受信電力の少なくとも１つの等高線を作成する。
【０１３６】
図１０は、等高線の概念図である。図１０を参照して、等高線ＣＴＲ１は、高さｚｊにおいて受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、最大距離ＬＭＡＸ＿１を半径とする電力値Ｐ１［ｄＢｍ］を有する等高線であり、等高線ＣＴＲ２は、高さｚｊにおいて受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、最大距離ＬＭＡＸ＿２を半径とする電力値Ｐ２［ｄＢｍ］を有する等高線であり、等高線ＣＴＲ３は、高さｚｊにおいて受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、最大距離ＬＭＡＸ＿３を半径とする電力値Ｐ３［ｄＢｍ］を有する等高線である。
【０１３７】
図１０においては、高さｚｊにおいて、３個の等高線ＣＴＲ１〜ＣＴＲ３が図示されているが、高さｚｊにおいて少なくとも１つの等高線ＣＴＲが作成される。同様にして、高さｚ１〜ｚ３においても、それぞれ、少なくとも１つの等高線が作成される。
【０１３８】
このように、この発明の実施の形態においては、任意の高さにおいて少なくとも１つの等高線ＣＴＲが作成される。
【０１３９】
作成手段１６は、上述した方法によって、任意の高さｚにおける受信電力の等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲを推定手段１７へ出力する。
【０１４０】
図１１は、受信電力と波源からの距離との関係を示す図である。図１１において、縦軸は、受信電力を表し、横軸は、波源Ｓ（受信電力の重心点ＲＳＳＩ＿Ｇ）からの距離を表す。
【０１４１】
電波の自由空間伝搬モデルでは、受信電力ＲＳＳＩｉは、図１１に示すように波源Ｓからの距離に反比例して減少する。
【０１４２】
電波の自由空間伝搬モデルにおいて、波源Ｓからの距離が１０００ｍ以上である場合、受信電力ＲＳＳＩｉは、距離の増加に対して直線状に低下する。そして、受信電力ＲＳＳＩｉが１０ｄＢｍだけ低下するときの距離は、３．１６ｋｍである。つまり、電波の自由空間伝搬モデルにおいては、受信電力ＲＳＳＩｉが１０ｄＢｍだけ低下するときの距離は、３．１６ｋｍ／１０ｄＢｍとなる。
【０１４３】
図１２は、受信電力の等高線ＣＴＲに基づいて任意の位置における受信電力の電力値を推定する方法を説明するための図である。
【０１４４】
図１２を参照して、等高線ＣＴＲは、作成手段１６によって作成された等高線である。受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）から等高線ＣＴＲよりも３．１６ｋｍだけ遠い位置においては、受信電力の電力値は、等高線ＣＴＲの電力値をＰとすると、Ｐ−１０［ｄＢｍ］となる。
【０１４５】
また、等高線ＣＴＲよりも３．１６ｋｍだけ受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）に近い位置においては、受信電力の電力値は、Ｐ＋１０［ｄＢｍ］となる。
【０１４６】
推定手段１７は、受信電力の等高線ＣＴＲ（図１０参照）を作成手段１６から受ける。そして、推定手段１７は、受信電力の等高線ＣＴＲ（図１０参照）に基づいて、受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）から等高線ＣＴＲよりも３．１６ｋｍだけ遠い位置における受信電力の電力値をＰ−１０［ｄＢｍ］と推定し、その推定した電力値Ｐ−１０［ｄＢｍ］を有する受信電力の等高線ＣＴＲ１−１を作成する。
【０１４７】
また、推定手段１７は、受信電力の等高線ＣＴＲ（図１０参照）に基づいて、等高線ＣＴＲよりも３．１６ｋｍだけ受信電力の重心点ＲＳＳＩ＿Ｇ（＝波源Ｓの位置）に近い位置における受信電力の電力値をＰ＋１０［ｄＢｍ］と推定し、その推定した電力値Ｐ＋１０［ｄＢｍ］を有する受信電力の等高線ＣＴＲ１−２を作成する。
【０１４８】
受信電力の等高線ＣＴＲ１−１，ＣＴＲ１−２は、等高線ＣＴＲと相似な形状を有する。
【０１４９】
なお、推定手段１７は、単位距離当たりの受信電力の電力値の変化割合ＲＳＳＩ＿ｒａｔｅ（＝１０／３．１６＝３．１６［ｄＢｍ／ｋｍ］）を求め、変化割合ＲＳＳＩ＿ｒａｔｅを用いて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定してもよい。より具体的には、推定手段１７は、受信電力の重心点ＲＳＳＩ＿Ｇから遠ざかるように受信電力の等高線ＣＴＲから距離Ｌ１だけ離れた位置における受信電力の電力値をＰ−ＲＳＳＩ＿ｒａｔｅ×Ｌ１によって推定する。また、推定手段１７は、受信電力の重心点ＲＳＳＩ＿Ｇに近づくように受信電力の等高線ＣＴＲから距離Ｌ１だけ離れた位置における受信電力の電力値をＰ＋ＲＳＳＩ＿ｒａｔｅ×Ｌ１［ｄＢｍ］によって推定する。これによって、推定手段１７は、ある高さｚにおいて、水平方向の任意の位置における受信電力の電力値を推定できる。
【０１５０】
推定手段１７は、電波の自由空間伝搬モデルに代えて奥村−秦モデルを用いて水平方向の任意の位置における受信電力の電力値を推定してもよい。
【０１５１】
推定手段１７は、奥村−秦モデルを用いて受信電力の電力値を推定する場合、奥村−秦モデルにおける複数の電波伝搬モデルを地形データと対応付けて保持しており、標高地図または航空写真を外部から受ける。そして、推定手段１７は、標高地図または航空写真に基づいて抽出した地形データに最も近い電波伝搬モデルを複数の電波伝搬モデルから選択し、その選択した電波伝搬モデルを用いて受信電力の電力値を推定する。
【０１５２】
図１３は、奥村−秦モデルにおける電界強度と距離との関係を示す図である。図１３において、縦軸は、電界強度を表し、横軸は、距離を表す。また、曲線ｋ４は、開放地における電界強度と距離との関係を示し、曲線ｋ５は、郊外における電界強度と距離との関係を示し、曲線ｋ６は、中小都市における電界強度と距離との関係を示し、曲線ｋ７は、大都市における電界強度と距離との関係を示す。
【０１５３】
図１３に示す曲線ｋ４〜ｋ７は、次式によって表される。
【０１５４】
【数２】
2019197957.tif 000004
【０１５５】
式（２）において、Ａ，Ｂは、曲線ｋ４〜ｋ７に共通であり、ａ（ｈｍ）およびＣは、エリア毎に異なる。また、ＬＣＨは、電界強度であり、ｄは、距離であり、ｈｍは、２次利用者（端末装置２）のアンテナ高さである。
【０１５６】
式（２）のＡ，Ｂは、次式によって表される。
【０１５７】
【数３】
2019197957.tif 000005
【０１５８】
式（３）において、ｆは、伝搬する電波の周波数であり、ｈｂは、１次利用者（波源Ｓ）のアンテナ高さである。
【０１５９】
電波の伝搬環境が開放地である場合（曲線ｋ４）、ａ（ｈｍ）およびＣは、次式によって表される。
【０１６０】
【数４】
2019197957.tif 000006
【０１６１】
電波の伝搬環境が郊外である場合（曲線ｋ５）、ａ（ｈｍ）およびＣは、次式によって表される。
【０１６２】
【数５】
2019197957.tif 000007
【０１６３】
電波の伝搬環境が中小都市である場合（曲線ｋ６）、ａ（ｈｍ）およびＣは、次式によって表される。
【０１６４】
【数６】
2019197957.tif 000008
【０１６５】
電波の伝搬環境が大都市である場合（曲線ｋ７）、ａ（ｈｍ）およびＣは、次式によって表される。
【０１６６】
【数７】
2019197957.tif 000009
【０１６７】
式（７）に示すように、大都市においては、周波数によってａ（ｈｍ）が異なる。
【０１６８】
推定手段１７は、式（２）〜式（７）を保持しており、地形データが開放地に最も近いとき、式（２），（３），（４）によって曲線ｋ４に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。
【０１６９】
また、推定手段１７は、地形データが郊外に最も近いとき、式（２），（３），（５）によって曲線ｋ５に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。
【０１７０】
更に、推定手段１７は、地形データが中小都市に最も近いとき、式（２），（３），（６）によって曲線ｋ６に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。
【０１７１】
更に、推定手段１７は、地形データが大都市に最も近いとき、式（２），（３），（７）によって曲線ｋ７に示す電波伝搬特性を演算し、その演算した電波伝搬特性に基づいて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。
【０１７２】
図１３に示す曲線ｋ４〜ｋ７においては、電界強度が距離に対して直線的に減少する領域があるので、電波の自由空間伝搬モデルにおいて説明したように、推定手段１７は、変化割合ＲＳＳＩ＿ｒａｔｅを算出し、その算出した変化割合ＲＳＳＩ＿ｒａｔｅを用いて、受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定してもよい。
【０１７３】
このように、推定手段１７は、電波の自由空間伝搬モデルおよび奥村−秦モデル等の電波の伝搬モデルを用いて受信電力の等高線ＣＴＲから任意の距離だけ離れた位置における受信電力の電力値を推定する。
【０１７４】
図１４は、端末装置２における電波の受信感度以下の領域において受信電力を推定する方法を説明するための図である。
【０１７５】
端末装置２は、受信感度以下の受信電力を取得できないため、受信感度以下の領域における受信電力の等高線ＣＴＲは、既に作成した受信電力の等高線ＣＴＲと電波の伝搬モデルとに基づいて作成する。ここで、電波の伝搬モデルとしては、上述した電波の自由空間伝搬モデルおよび奥村−秦モデル等が用いられる。
【０１７６】
上述したように、電波の自由空間伝搬モデルにおいては、波源Ｓからの距離が１０００ｍ以上である場合、受信電力ＲＳＳＩｉは、距離の増加に対して直線状に低下する。そして、受信電力ＲＳＳＩｉが１０ｄＢｍだけ低下するときの距離は、３．１６ｋｍ／１０ｄＢｍとなる。
【０１７７】
従って、作成手段１６は、電波の受信感度以下の領域において、既に作成した受信電力の等高線ＣＴＲ（電力値Ｐの等高線）から３．１６ｋｍだけ受信電力の重心点ＲＳＳＩ＿Ｇから遠くなる方向に離れ、かつ、受信電力の等高線ＣＴＲと相似形を有する受信電力の等高線ＣＴＲ＿ＬＳをＰ−１０の電力値を有する等高線ＣＴＲとして作成する（図１４参照）。
【０１７８】
この場合、受信電力の等高線ＣＴＲ＿ＬＳにおける電力値は、Ｐ−１０に限らず、電力値Ｐから所望の電力値ΔＰだけ低下した電力値であればよい。従って、作成手段１６は、電波の自由空間伝搬モデルにおいて、電力値ＰよりもΔＰだけ低下するときの距離を求め、その求めた距離だけ受信電力の重心点ＲＳＳＩ＿Ｇから遠くなる方向に離れ、かつ、受信電力の等高線ＣＴＲと相似形を有する受信電力の等高線ＣＴＲ＿ＬＳをＰ−ΔＰの電力値を有する等高線ＣＴＲとして作成する。
【０１７９】
また、作成手段１６は、同様にして、奥村−秦モデルを用いて電力値ＰよりもΔＰだけ低下するときの距離を求め、その求めた距離だけ受信電力の重心点ＲＳＳＩ＿Ｇから遠くなる方向に離れ、かつ、受信電力の等高線ＣＴＲと相似形を有する受信電力の等高線ＣＴＲ＿ＬＳをＰ−ΔＰの電力値を有する等高線ＣＴＲとして作成する。
【０１８０】
従って、作成手段１６は、一般的に、電波の受信感度以下の領域において、既に作成した受信電力の等高線ＣＴＲ上の電力値よりも所望の電力値だけ低い受信電力になる電力低下距離（３．１６ｋｍ等）を電波の伝搬モデルに基づいて求め、既に作成した受信電力の等高線ＣＴＲから電力低下距離だけ離れ、かつ、既に作成した受信電力の等高線ＣＴＲと相似形を有する等高線ＣＴＲを電波の受信感度以下の領域における受信電力の等高線ＣＴＲとして作成する。
【０１８１】
このように、作成手段１６は、電波の受信感度以下の領域においても受信電力の等高線ＣＴＲを作成するので、受信電力ＲＳＳＩｉの広い範囲において受信電力の等高線ＣＴＲを作成できる。
【０１８２】
図１５は、図１に示す電力推定装置１の動作を説明するためのフローチャートである。図１５を参照して、電力推定装置１の動作が開始されると、ｎ＝１が設定され（ステップＳ１）、ｓ＝１が設定される（ステップＳ２）。
【０１８３】
そして、抽出手段１３は、高さ方向の領域ＳＰＣ（ｓ）に含まれる複数のモニター情報を記録手段１２に記録されたモニター情報から抽出する（ステップＳ３）。
【０１８４】
その後、抽出手段１３は、作成すべき受信電力の等高線ＣＴＲの電力値Ｐｎに応じて処理範囲ＲＥＧ＿ＰＲＳを決定し（ステップＳ４）、その決定した処理範囲ＲＥＧ＿ＰＲＳ内に存在するｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを抽出する（ステップＳ５）。そして、抽出手段１３は、その抽出したｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを処理手段１４へ出力する（ステップＳ６）。
【０１８５】
処理手段１４は、ｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋを抽出手段１３から受け、その受けたｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに基づいて受信電力の重心点ＲＳＳＩ＿Ｇ（ｓ）および最大距離ＬＭＡＸ（ｓ，Ｐｎ）を算出する（ステップＳ７）。
【０１８６】
引き続いて、処理手段１４は、ｋ個のモニター情報ＭＮＴ１〜ＭＮＴｋに基づいて領域ＳＰＣ（ｓ）の高さ位置ｚｓ＝（ｚ１＋・・・＋ｚｋ）／ｋを算出する（ステップＳ８）。なお、処理手段１４は、ｚｓ＝（ｗ１ｚ１＋ｗ２ｚ２＋・・・＋ｗｋｚｋ）／（ｗ１＋ｗ２＋・・・＋ｗｋ）によって高さ位置ｚｓを算出してもよい。ここで、ｗｉ＝ＲＳＳＩｉ（ｉ＝１〜ｋ）である。
【０１８７】
その後、処理手段１４は、ｓ＝ｍ（高さ方向の領域ＳＰＣの最大個数）であるか否かを判定する（ステップＳ９）。
【０１８８】
ステップＳ９において、ｓ＝ｍでないと判定されたとき、処理手段１４は、ｓ＝ｓ＋１を設定する（ステップＳ１０）。その後、一連の動作は、ステップＳ３へ戻り、ステップＳ９において、ｓ＝ｍであると判定されるまで、ステップＳ３〜ステップＳ１０が繰り返し実行される。
【０１８９】
そして、ステップＳ９において、ｓ＝ｍであると判定されると、処理手段１４は、ｍ個の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離ＬＭＡＸ＿１（Ｐｎ）〜ＬＭＡＸ＿ｍ（Ｐｎ）およびｍ個の高さｚ１〜ｚｍを演算手段１５へ出力する。
【０１９０】
演算手段１５は、ｍ個の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離ＬＭＡＸ＿１（Ｐｎ）〜ＬＭＡＸ＿ｍ（Ｐｎ）およびｍ個の高さｚ１〜ｚｍを処理手段１４から受け、その受けたｍ個の重心点ＲＳＳＩ＿Ｇ＿１〜ＲＳＳＩ＿Ｇ＿ｍ、ｍ個の最大距離ＬＭＡＸ＿１（Ｐｎ）〜ＬＭＡＸ＿ｍ（Ｐｎ）およびｍ個の高さｚ１〜ｚｍに基づいて、上述した方法によって、重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧと高さｚとの関係式ＥＱ１と、重心点ＲＳＳＩ＿Ｇのｙ座標ｙＧと高さｚとの関係式ＥＱ２と、最大距離ＬＭＡＸと高さｚとの関係式ＥＱ３（等高線ＣＴＲの電力値Ｐｎに対応した関係式ＥＱ３−ｎ）とを求める（ステップＳ１１）。そして、演算手段１５は、関係式ＥＱ１〜ＥＱ３を作成手段１６へ出力する。
【０１９１】
作成手段１６は、関係式ＥＱ１〜ＥＱ３を演算手段１５から受ける。そして、作成手段１６は、任意の高さｚａを外部から受けると、関係式ＥＱ１，ＥＱ２に基づいて、それぞれ、任意の高さｚａにおける重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧ（ｚａ）およびｙ座標ｙＧ（ｚａ）を求め、関係式ＥＱ３（等高線ＣＴＲの電力値Ｐｎに対応した関係式ＥＱ３−ｎ）に基づいて任意の高さｚａにおける最大距離ＬＭＡＸ＿ｎ（ｚａ）を求める（ステップＳ１２）。
【０１９２】
そうすると、作成手段１６は、点（ｘＧ（ｚａ），ｙＧ（ｚａ））を中心とし、最大距離ＬＭＡＸ＿ｎ（ｚａ）を半径とする円形形状を有する受信電力の等高線ＣＴＲを作成する（ステップＳ１３）。
【０１９３】
そして、作成手段１６は、ｎ＝Ｎ（Ｎは、作成すべき受信電力の等高線ＣＴＲの電力値Ｐｎの最大個数）であるか否かを判定する（ステップＳ１４）。
【０１９４】
ステップＳ１４において、ｎ＝Ｎでないと判定されたとき、作成手段１６は、ｎ＝ｎ＋１を設定する（ステップＳ１５）。その後、一連の動作は、ステップＳ２へ移行し、ステップＳ１４において、ｎ＝Ｎであると判定されるまで、ステップＳ２〜ステップＳ１５が繰り返し実行される。
【０１９５】
そして、ステップＳ１４において、ｎ＝Ｎであると判定されると、一連の動作は、終了する。なお、図１５に示すフローチャートは、上述した分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々において実行される。
【０１９６】
図１６は、図１５に示すフローチャートのステップの実行状態を示す概念図である。図１６を参照して、図１５に示すフローチャートのステップＳ３〜ステップＳ１０がｍ回実行されると、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐｎにおいて、高さ方向のｍ個の領域ＳＰＣ（１）〜ＳＰＣ（ｍ）の全てについて、高さ位置ｚ、受信電力の重心点ＲＳＳＩ＿Ｇ（（ｘＧ（ｚ１），ｙＧ（ｚ１），ｘＧ（ｚ２），ｙＧ（ｚ２），・・・，ｘＧ（ｚｍ），ｙＧ（ｚｍ））および最大距離ＬＭＡＸ（ｚ１），ＬＭＡＸ（ｚ２），・・・，ＬＭＡＸ（ｚｍ）が算出される（図１６の（ａ）参照）。
【０１９７】
また、ステップＳ１１が実行されると、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐｎにおいて、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧとの関係式ＥＱ１（＝ｘＧ（ｚ））、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙＧとの関係式ＥＱ２（＝ｙＧ（ｚ））、および高さｚおよび電力値Ｐｎと最大距離ＬＭＡＸとの関係式ＥＱ３（＝ＬＭＡＸ＿ｎ（Ｐｎ，ｚ））が求められる（図１６の（ｂ）参照）。ここで、ＬＭＡＸ＿ｎ（Ｐｎ，ｚ）は、等高線ＣＴＲの１つの電力値Ｐｎに対応し、かつ、最大距離ＬＭＡＸと高さｚとの関係を示す関係式である。
【０１９８】
更に、ステップＳ１２が実行されると、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐｎにおいて、任意の高さｚａにおける受信電力の重心点ＲＳＳＩ＿Ｇの座標（（ｘＧ（ｚａ），ｙＧ（ｚａ））および任意の高さｚａにおける最大距離ＬＭＡＸ＿ｎ（Ｐｎ，ｚａ）が求められる（図１６の（ｃ）参照）。この場合、重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧ（ｚａ）は、関係式ＥＱ１＝ｘＧ（ｚ）にｚ＝ｚａを代入して求められ、重心点ＲＳＳＩ＿Ｇのｙ座標ｙＧ（ｚａ）は、関係式ＥＱ２＝ｙＧ（ｚ）にｚ＝ｚａを代入して求められ、最大距離ＬＭＡＸｎ（Ｐｎ，ｚａ）は、関係式ＥＱ３（等高線ＣＴＲの電力値Ｐｎに対応した関係式ＥＱ３−ｎ）＝ＬＭＡＸ＿ｎ（Ｐｎ，ｚ）にｚ＝ｚａを代入して求められる。
【０１９９】
更に、ステップＳ１３が実行されると、高さｚａにおいて１つの電力値Ｐｎを有する等高線ＣＴＲ＿ｎが作成される。等高線ＣＴＲ＿ｎは、受信電力の重心点（（ｘＧ（ｚａ），ｙＧ（ｚａ））を中心とし、最大距離ＬＭＡＸ＿ｎ（Ｐｎ，ｚａ）を半径とする円形形状を有する（図１６の（ｄ）参照）。
【０２００】
更に、ステップＳ２〜ステップＳ１５がＮ回実行されると、高さｚａにおいて、それぞれ、電力値Ｐ１〜ＰＮを有するＮ個の等高線ＣＴＲ＿１〜ＣＴＲ＿Ｎが作成される（図１６の（ｅ）参照）。等高線ＣＴＲ＿１〜ＣＴＲ＿Ｎは、受信電力の重心点（（ｘＧ（ｚａ），ｙＧ（ｚａ））を中心とし、それぞれ、最大距離ＬＭＡＸ＿１（Ｐ１，ｚａ）〜ＬＭＡＸ＿Ｎ（ＰＮ，ｚａ）を半径とする円形形状を有する。そして、最大距離ＬＭＡＸ＿１（Ｐ１，ｚａ）〜ＬＭＡＸ＿Ｎ（ＰＮ，ｚａ）は、それぞれ、電力値Ｐ１〜Ｐｎに対応する関係式ＥＱ３＿１＝ＬＭＡＸ＿１（Ｐ１，ｚ）〜ＥＱ３＿Ｎ＝ＬＭＡＸ＿Ｎ（ＰＮ，ｚ）を用いて求められる。
【０２０１】
従って、図１５に示すフローチャートは、高さｚａにおいて、１つの電力値Ｐｎ毎に受信電力の等高線ＣＴＲを作成する動作を示すものである。
【０２０２】
図１７は、図１に示す電力推定装置１の動作を説明するための別のフローチャートである。
【０２０３】
図１７に示すフローチャートは、図１５に示すフローチャートのステップＳ１２，Ｓ１３をステップＳ１６に変えたものであり、その他は、図１５に示すフローチャートと同じである。
【０２０４】
図１７を参照して、電力推定装置１の動作が開始されると、上述したステップＳ１〜Ｓ１１，Ｓ１４，Ｓ１５が順次実行される。
【０２０５】
そして、ステップＳ１４において、ｎ＝Ｎであると判定されると、作成手段１６は、関係式ＥＱ１〜ＥＱ３に基づいて、任意の高さｚａにおける任意の電力値Ｐｎを有する受信電力の複数の等高線ＣＴＲを作成する（ステップＳ１６）。これによって、電力推定装置１の動作が終了する。
【０２０６】
図１８は、図１７のステップＳ１６の詳細な動作を説明するためのフローチャートである。
【０２０７】
図１８を参照して、図１７に示すステップＳ１４において、ｎ＝Ｎであると判定されると、作成手段１６は、高さｚａと電力値ＰｎとのＶ（Ｖは，２≦Ｖ≦Ｎを満たす整数）個の組（ｚａ１，Ｐ１），（ｚａ２，Ｐ２），・・・，（ｚａＶ，ＰＶ）を外部から受ける（ステップＳ１６１）。
【０２０８】
そして、作成手段１６は、ｖ＝１を設定する（ステップＳ１６２）。その後、作成手段１６は、関係式ＥＱ１，ＥＱ２に基づいて、それぞれ高さｚａｒにおける受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧ（ｚａｒ）およびｙ座標ｙＧ（ｚａｒ）を求め、関係式ＥＱ３（等高線ＣＴＲの電力値Ｐｖに対応する関係式ＥＱ３−ｖ）に基づいて、電力値Ｐｖおよび高さｚａｖにおける最大距離ＬＭＡＸ＿ｖ（Ｐｖ，ｚａｖ）を求める（ステップＳ１６３）。
【０２０９】
そうすると、作成手段１６は、受信電力の重心点（ｘＧ（ｚａｒ），ｙＧ（ｚａｒ））を中心とし、最大距離ＬＭＡＸ＿ｖ（Ｐｖ，ｚａｖ）を半径とする受信電力の等高線ＣＴＲ＿ｖを作成する（ステップＳ１６４）。
【０２１０】
そして、作成手段１６は、ｒ＝Ｒであるか否かを判定する（ステップＳ１６５）。ステップＳ１６５において、ｖ＝Ｖでないと判定されたとき、作成手段１６は、ｖ＝ｖ＋１を設定する（ステップＳ１６６）。その後、一連の動作は、ステップＳ１６３に移行し、ステップＳ１６５において、ｖ＝Ｖであると判定されるまで、ステップＳ１６３〜ステップＳ１６６が繰り返し実行される。
【０２１１】
そして、ステップＳ１６５において、ｖ＝Ｖであると判定されると、一連の動作は、図１７の“終了”へ移行する。なお、図１７に示すフローチャート（図１８に示すフローチャートを含む）は、上述した分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々において実行される。
【０２１２】
図１９は、図１７および図１８に示すフローチャートのステップの実行状態を示す概念図である。図２０は、図１７および図１８に示すフローチャートの実行が完了した状態を示す概念図である。
【０２１３】
図１９を参照して、図１７に示すフローチャートのステップＳ３〜ステップＳ１０がｍ回実行されると、図１６において説明したとおり、作成すべき受信電力の等高線ＣＴＲの１つの電力値Ｐｎにおいて、高さ方向のｍ個の領域ＳＰＣ（１）〜ＳＰＣ（ｍ）の全てについて、高さ位置ｚ、受信電力の重心点ＲＳＳＩ＿Ｇ（（ｘＧ（ｚ１），ｙＧ（ｚ１），ｘＧ（ｚ２），ｙＧ（ｚ２），・・・，ｘＧ（ｚｍ），ｙＧ（ｚｍ））および最大距離ＬＭＡＸ（ｚ１），ＬＭＡＸ（ｚ２），・・・，ＬＭＡＸ（ｚｍ）が算出される（図１９の（ａ）参照）。
【０２１４】
また、図１７に示すフローチャートのステップＳ２〜Ｓ１１，Ｓ１４，Ｓ１５がＮ回実行されると、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｘ座標ｘＧとの関係式ＥＱ１（＝ｘＧ（ｚ））、高さｚと受信電力の重心点ＲＳＳＩ＿Ｇのｙ座標ｙＧとの関係式ＥＱ２（＝ｙＧ（ｚ））、および電力値Ｐ１〜ＰＮに対応して、それぞれ、高さｚおよび電力値Ｐｎと最大距離ＬＭＡＸとの関係式ＥＱ３＿１（＝ＬＭＡＸ＿１（Ｐ１，ｚ））〜ＥＱ３＿Ｎ（＝ＬＭＡＸ＿Ｎ（ＰＮ，ｚ））が求められる（図１９の（ｂ）参照）。ここで、ＬＭＡＸ＿１（Ｐ１，ｚ）〜ＬＭＡＸ＿Ｎ（ＰＮ，ｚ）は、それぞれ、等高線ＣＴＲの電力値Ｐ１〜ＰＮに対応し、かつ、最大距離ＬＭＡＸと高さｚとの関係を示す関係式である。
【０２１５】
更に、図１８に示すステップＳ１６３が実行されると、関係式ＥＱ１，ＥＱ２に基づいて、高さｚａｖにおける受信電力の重心点の座標（ｘＧ（ｚａｖ），ｙＧ（ｚａｖ））が求められ、電力値Ｐｖに対応する関係式ＥＱ３＿ｖ（ＥＱ３＿１〜ＥＱ３＿Ｎのいずれか）に基づいて、高さｚａｖにおける最大距離ＬＭＡＸ＿ｖ（Ｐｖ，ｚａｖ）が求められる（図１９の（ｃ）参照）。
【０２１６】
更に、図１８に示すステップＳ１６４が実行されると、高さｚａｖにおいて１つの電力値Ｐｖを有する等高線ＣＴＲ＿ｖが作成される。等高線ＣＴＲ＿ｖは、受信電力の重心点（（ｘＧ（ｚａｖ），ｙＧ（ｚａｖ））を中心とし、最大距離ＬＭＡＸ＿ｖ（Ｐｖ，ｚａｖ）を半径とする円形形状を有する（図１９の（ｄ）参照）。
【０２１７】
更に、図１７および図１８に示すフローチャートの実行が完了した状態においては、高さｚａ１，ｚａ２，ｚａ３，・・・，ｚａＶにおいて、それぞれ、電力値Ｐ１，Ｐ２，Ｐ３，・・・，ＰＶを有する受信電力の等高線ＣＴＲ＿１，ＣＴＲ＿２，ＣＴＲ＿３，・・・，ＣＴＲ＿Ｖが作成される。等高線ＣＴＲ＿１は、受信電力の重心点（（ｘＧ（ｚａ１），ｙＧ（ｚａ１））を中心とし、最大距離ＬＭＡＸ＿１（Ｐ１，ｚａ１）を半径とする円形形状を有する。等高線ＣＴＲ＿２は、受信電力の重心点（（ｘＧ（ｚａ２），ｙＧ（ｚａ２））を中心とし、最大距離ＬＭＡＸ＿２（Ｐ２，ｚａ２）を半径とする円形形状を有する。等高線ＣＴＲ＿３は、受信電力の重心点（（ｘＧ（ｚａ３），ｙＧ（ｚａ３））を中心とし、最大距離ＬＭＡＸ＿３（Ｐ３，ｚａ３）を半径とする円形形状を有する。以下、同様にして、等高線ＣＴＲ＿Ｖは、受信電力の重心点（（ｘＧ（ｚａＶ），ｙＧ（ｚａＶ））を中心とし、最大距離ＬＭＡＸ＿Ｖ（ＰＶ，ｚａＶ）を半径とする円形形状を有する（図２０参照）。
【０２１８】
受信電力の等高線ＣＴＲ＿１〜ＣＴＲ＿Ｖは、異なる高さｚａ１〜ｚａＶにおいて作成され、それぞれ、異なる半径（最大距離ＬＭＡＸ＿１（Ｐ１，ｚａ１）〜ＬＭＡＸ＿Ｖ（ＰＶ，ｚａＶ））および異なる電力値Ｐ１〜ＰＶを有する。
【０２１９】
従って、図１７および図１８に示すフローチャートは、任意の高さにおいて任意の電力値を有する受信電力の複数の等高線ＣＴＲを作成する動作を示すものである。
【０２２０】
図２１は、図１に示す電力推定装置１の別の動作を説明するためのフローチャートである。
【０２２１】
図２１を参照して、電力推定装置１の動作が開始されると、作成手段１６は、図１５に示すフローチャート、または図１７に示すフローチャート（図１８に示すフローチャートを含む）を実行する（ステップＳ２１）。そして、作成手段１６は、作成した受信電力の等高線ＣＴＲを推定手段１７へ出力する。
【０２２２】
推定手段１７は、受信電力の等高線ＣＴＲを作成手段１６から受け、その受けた受信電力の等高線ＣＴＲに基づいて受信電力を推定する（ステップＳ２２）。この場合、推定手段１７は、上述した方法によって、電波の伝搬モデルを用いて受信電力を推定する。そして、ステップＳ２２の後、電力推定装置１の動作が終了する。
【０２２３】
なお、電力推定装置１の動作は、ソフトウェアによって実行されてもよい。この場合、電力推定装置１は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）およびＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）を備える。
【０２２４】
ＲＯＭは、図１５に示すフローチャートからなるプログラムＰｒｏｇ＿Ａ、図１７に示すフローチャート（図１８に示すフローチャートを含む）からなるプログラムＰｒｏｇ＿Ｂ、および図２１に示すフローチャートからなるプログラムＰｒｏｇ＿Ｃのいずれかを格納する。
【０２２５】
そして、ＣＰＵは、プログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿ＣのいずれかをＲＯＭから読み出し、その読み出したプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかを実行して受信電力の等高線ＣＴＲの作成、または受信電力の等高線ＣＴＲの作成および受信電力の推定を行う。
【０２２６】
この場合、ＲＡＭは、端末装置２から受信したモニター情報ＭＮＴおよび関係式ＥＱ１〜ＥＱ３等の記憶に用いられるとともに各種の演算に用いられる。
【０２２７】
また、プログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかは、ＣＤ，ＤＶＤ等の記録媒体に記録されて流通してもよい。この場合、ＣＰＵは、装着された記録媒体からプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかを読み出して実行し、受信電力の等高線ＣＴＲの作成、または受信電力の等高線ＣＴＲの作成および受信電力の推定を行う。従って、プログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれかを記録したＣＤ，ＤＶＤ等の記録媒体は、コンピュータ（ＣＰＵ）が読み取り可能な記録媒体である。
【０２２８】
上記においては、受信電力の等高線ＣＴＲは、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離ＬＭＡＸを半径とする円形形状を有すると説明したが、この発明の実施の形態においては、これに限らず、受信電力の等高線ＣＴＲは、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離ＬＭＡＸを長径とする楕円形状を有していてもよい。
【０２２９】
図２２は、楕円形状を有する受信電力の等高線ＣＴＲの作成方法を説明するための図である。
【０２３０】
図２２において、白四角は、電力値Ｐ以上の受信電力を有する測定点を示し、黒四角は、電力値Ｐ未満の受信電力を有する測定点を示し、白丸は、受信電力の重心点ＲＳＳＩ＿Ｇを示す。
【０２３１】
図２２を参照して、作成手段１６は、上述した方法によって、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離ＬＭＡＸを半径とする円形形状を有する受信電力の等高線ＣＴＲ５を作成する。
【０２３２】
そして、作成手段１６は、受信電力の等高線ＣＴＲ５に基づいて、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離ＬＭＡＸを長径とするとともに距離Ｌ’ＭＡＸを短径とする楕円形状を有する受信電力の等高線ＣＴＲ６を作成する。ここで、距離Ｌ’ＭＡＸは、電力値Ｐ以上の受信電力を有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇからの距離が最大距離ＬＭＡＸの次に長い測定点２−２と受信電力の重心点ＲＳＳＩ＿Ｇとの距離である。
【０２３３】
また、この発明の実施の形態においては、受信電力の等高線ＣＴＲは、多角形の形状を有していてもよい。
【０２３４】
図２３は、多角形の形状を有する受信電力の等高線ＣＴＲの作成方法を説明するための図である。
【０２３５】
図２３において、白四角は、電力値Ｐ以上の受信電力を有する測定点を示し、黒四角は、電力値Ｐ未満の受信電力を有する測定点を示し、白丸は、受信電力の重心点ＲＳＳＩ＿Ｇを示す。
【０２３６】
図２３を参照して、作成手段１６は、上述した方法によって、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離ＬＭＡＸを半径とする円形形状を有する受信電力の等高線ＣＴＲ５を作成する。
【０２３７】
そして、作成手段１６は、円形形状を有する受信信号の等高線ＣＴＲ５に内接する正八角形３０を形成する（図２３の（ａ）参照）。
【０２３８】
その後、作成手段１６は、正八角形の各頂点が、電力値Ｐ以上の受信電力を有し、かつ、受信電力の重心点ＲＳＳＩ＿Ｇからできるだけ遠い位置に存在する測定点２−１〜２−８に位置するように正八角形３０を変形し、八角形の形状を有する受信電力の等高線ＣＴＲ７を作成する（図２３の（ｂ）参照）。
【０２３９】
作成手段１６は、図２３において説明した作成方法と同じ作成方法によって、八角形以外の多角形の形状を有する受信電力の等高線ＣＴＲを作成する。
【０２４０】
このように、作成手段１６は、円形形状を有する受信電力の等高線ＣＴＲに限らず、楕円形状を有する受信電力の等高線ＣＴＲまたは多角形の形状を有する受信電力の等高線ＣＴＲを作成する。
【０２４１】
従って、図１５に示すフローチャートのステップＳ１３または図１８に示すフローチャートのステップＳ１６４において、作成手段１６は、円形形状を有する受信電力の等高線ＣＴＲに代えて、楕円形状を有する受信電力の等高線ＣＴＲまたは多角形の形状を有する受信電力の等高線ＣＴＲを作成してもよく、一般的には、作成手段１６は、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、受信電力の重心点ＲＳＳＩ＿Ｇからの距離が最も長い最大距離ＬＭＡＸを半径とする円形形状を有する第１の等高線、受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離ＬＭＡＸを長径とする楕円形状を有する第２の等高線、および受信電力の重心点ＲＳＳＩ＿Ｇを中心とし、かつ、最大距離ＬＭＡＸを中心から１つの頂点までの距離とする多角形の形状を有する第３の等高線のいずれかを任意の高さｚにおいて作成すればよい。
【０２４２】
モニター情報ＭＮＴは、共用周波数ｆｃｏｍ＿ｉを含むため、電力推定装置１は、上述した方法によって、各共用周波数ｆｃｏｍ＿ｉを用いて電波を送信した場合における受信電力の等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲに基づいて受信電力を推定する。
【０２４３】
また、電力推定装置１は、例えば、端末装置２が１５：００に無線通信を行う場合、１４：００〜１５：００までの間に複数の端末装置２から送信された複数のモニター情報ＭＮＴに基づいて、上述した方法によって、任意の高さｚにおける受信電力の等高線ＣＴＲを作成し（端末装置２における電波の受信感度以下の領域における受信電力の等高線ＣＴＲの作成も含む）、その作成した受信電力の等高線ＣＴＲに基づいて、各場所における受信電力を推定し、その推定結果に基づいて、共用周波数ｆｃｏｍ＿ｉを用いて無線通信を行うときの通信条件（共用周波数ｆｃｏｍ、場所および送信電力等）を決定する。
【０２４４】
そして、各端末装置２は、電力推定装置１によって決定された通信条件を用いて無線通信を行う。
【０２４５】
その結果、電力推定装置１によって推定された受信電力は、フェージングを抑制して推定された受信電力であるため、送信電力は、フェージングが抑制された電波環境において信号を送信先へ送信できるように決定される。従って、スループットを高くして信号を送信先へ送信できる。
【０２４６】
また、上記においては、分割領域ＲＥＧ＿ＤＶ１〜ＲＥＧ＿ＤＶｑの各々において、高さ方向の領域ＳＰＣ１〜ＳＰＣｍの各々は、高さ方向の同じ長さｈを有すると説明したが、この発明の実施の形態においては、これに限らず、建物２０の高さＨＢ（または平均高さ）以下の位置に存在する領域ＳＰＣ１，ＳＰＣ２においては、高さ方向に長さｈ１を有し、建物２０の高さＨＢ（または平均高さ）よりも高い位置に存在する領域ＳＰＣｍにおいては、高さ方向に長さｈ１よりも長い長さｈ２を有するようにしてもよい。
【０２４７】
建物２０の高さＨＢ（または平均高さ）よりも高い領域においては、電波が伝搬するときの障害物が少ないので、高さ方向に長さを短くしなくても、電波の受信電力の電力分布を正確に推定することができる。
【０２４８】
一方、建物２０の高さＨＢ（または平均高さ）よりも低い領域においては、電波が伝搬するときの障害物が多いので、高さ方向に長さを短くした方が電波の受信電力の電力分布を正確に推定することができる。
【０２４９】
図２４は、この発明の実施の形態によるデータ構造の概略図である。図２４を参照して、データ構造Ｄ＿ＳＴＲ１は、端末装置ＩＤ、時間情報、端末装置２の位置情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤおよび受信電力が相互に対応付けられた構造を有する。端末装置ＩＤ、時間情報、端末装置２の位置情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤおよび受信電力は、モニター情報ＭＮＴを構成する。
【０２５０】
端末装置ＩＤは、端末装置２の識別情報である。時間情報は、受信電力が検出された時刻を示す。そして、時間情報は、年／月／日／時間／分／秒（ＹＹＹＹ／ＭＭ／ＤＤ／ＨＨ／ＭＭ／ＳＳ）の形式によって表される。端末装置２の位置情報は、緯度と、経度と、高さとを含む。緯度および経度は、ＧＰＳによって検出され、高さは、高度計によって検出される。そして、高さは、端末装置２が静止または移動している地面からの高さを示す。
【０２５１】
周波数帯域ＩＤは、複数の端末装置２が共用する周波数帯域の識別情報である。サブキャリアＮｏ．は、サブキャリアの番号を示す。波源ＩＤは、１次利用者の電波を発する波源の識別情報である。受信電力は、各端末装置２が受信した電波の受信電力である。
【０２５２】
端末装置ＩＤは、抽出手段１３が、対応する受信電力が同じ端末装置２で検出されたか否かを判定するのに用いられる。時間情報は、抽出手段１３が、対応する受信電力が同じ時刻に検出されたか否かを判定するのに用いられる。端末装置２の位置情報は、抽出手段１３が、対応する受信電力が同じ位置で検出されたか否か、または対応する受信電力が所定の領域ＳＰＣ内で検出された否かを判定するのに用いられる。また、位置情報は、処理手段１４が、複数の端末装置２間の距離、端末装置２と波源との距離および受信電力の重心点から最も遠い位置に存在する端末装置を検出するのに用いられる。
【０２５３】
周波数帯域ＩＤは、抽出手段１３が、同じ周波数帯域の受信電力を抽出するために用いられる。サブキャリアＮｏ．は、抽出手段１３が、対応する受信電力を有する電波が伝送されたサブキャリアを判定するのに用いられる。波源ＩＤは、抽出手段１３が、同じ波源から送信された電波の受信電力を抽出するのに用いられる。
【０２５４】
受信電力は、処理手段１４が、複数の受信電力のうちの最大の受信電力を検出するのに用いられる。また、受信電力は、作成手段１６が、作成すべき電力の等高線の電力値以上の受信電力を検出するのに用いられる。
【０２５５】
複数のモニター情報ＭＮＴに含まれる複数の受信電力のうち、受信電力の最大値を含むｋ個のモニター情報は、抽出手段１３によって複数のモニター情報から抽出されるのに用いられる。
【０２５６】
抽出手段１３によって抽出されたｋ個のモニター情報は、処理手段１４が受信電力によって重み付けされたｋ個の端末装置２の位置の平均である受信電力の重心点（ｘＧ，ｙＧ）を求めるのに用いられる。
【０２５７】
抽出手段１３によって抽出されたｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力と、受信電力の等高線における電力値以上の受信電力に対応付けられた位置情報とは、処理手段１４が、受信電力の重心点（ｘＧ，ｙＧ）から最も遠い位置に存在する端末装置と受信電力の重心点（ｘＧ，ｙＧ）との距離である最大距離ＬＭＡＸを求めるのに用いられる。
【０２５８】
処理手段１４によって求められた受信電力の重心点（ｘＧ，ｙＧ）は、演算手段１５が受信電力の重心点（ｘＧ，ｙＧ）と高さｚとの関係を示す第１の関係式（関係式ＥＱ１，ＥＱ２）を求めるのに用いられる。
【０２５９】
処理手段１４によって求められた最大距離ＬＭＡＸは、演算手段１５が、作成すべき受信電力の等高線の電力値に対応し、かつ、最大距離ＬＭＡＸと高さｚとの関係を示す第２の関係式（関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎを含む））を求めるのに用いられる。
【０２６０】
演算手段１５によって求められた第１の関係式（関係式ＥＱ１，ＥＱ２）および第２の関係式（関係式ＥＱ３（ＥＱ３＿１〜ＥＱ３＿Ｎを含む））は、作成手段１６が、任意の高さにおいて、任意の電力値を有する受信電力の等高線を作成するのに用いられる。
【０２６１】
第１の等高線は、受信電力の重心点を中心とし、かつ、最大距離を半径とする円形形状を有し、第２の等高線は、受信電力の重心点を中心とし、かつ、最大距離を長径とする楕円形状を有し、第３の等高線は、受信電力の重心点を中心とし、かつ、最大距離を中心から１つの頂点までの距離とする多角形の形状を有する。
【０２６２】
好ましくは、抽出手段１３によって抽出されたｋ個のモニター情報に含まれるｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値よりも小さい受信電力と、受信電力の等高線における電力値よりも小さい受信電力に対応付けられた位置情報とは、推定手段１７が、受信電力の等高線に基づいて、対象領域ＳＰＣ内の所望の場所において電波の受信電力を推定するのに用いられる。
【０２６３】
好ましくは、受信電力の等高線の電力値を有する受信電力は、作成手段１６が、端末装置２における電波の受信感度以下の領域において、既に作成した受信電力の等高線上の電力値よりも所望の電力値だけ低い受信電力になる電力低下距離を電波の伝搬モデルに基づいて求め、既に作成した受信電力の等高線から電力低下距離だけ離れ、かつ、既に作成した受信電力の等高線と相似形を有する等高線を電波の受信感度以下の領域における受信電力の等高線として作成するのに用いられる。
【０２６４】
データ構造Ｄ＿ＳＴＲ１は、電力推定装置１（またはプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれか）が上述した方法によって受信電力の等高線を作成するのに用いられる。そして、データ構造Ｄ＿ＳＴＲ１は、端末装置ＩＤ、時間情報、端末装置２の位置情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤおよび受信電力の全てを含む必要はなく、端末装置２の位置情報、周波数帯域ＩＤおよび受信電力を少なくとも含み、端末装置２の位置情報、周波数帯域ＩＤおよび受信電力が相互に対応付けられた構造を有していればよい。データ構造Ｄ＿ＳＴＲ１が端末装置２の位置情報、周波数帯域ＩＤおよび受信電力を少なくとも含んでいれば、電力推定装置１（またはプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれか）が、上述した方法によって、任意の高さにおける受信電力の等高線を作成することができるからである。この場合、端末装置２の位置情報、周波数帯域ＩＤおよび受信電力は、モニター情報を構成する。
【０２６５】
図２５は、この発明の実施の形態による別のデータ構造の概略図である。この発明の実施の形態によるデータ構造は、図２５に示すデータ構造Ｄ＿ＳＴＲ２であってもよい。
【０２６６】
図２５を参照して、データ構造Ｄ＿ＳＴＲ２は、推定１次利用者ＩＤ、時間情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤ、等高線の中心位置、等高線の電力値および等高線の半径が相互に対応付けられた構造を有する。周波数帯域ＩＤ、サブキャリアＮｏ．および波源ＩＤについては、図２４において説明したとおりである。
【０２６７】
時間情報は、受信電力の等高線が作成された時刻を示す。そして、時間情報は、年／月／日／時間／分／秒（ＹＹＹＹ／ＭＭ／ＤＤ／ＨＨ／ＭＭ／ＳＳ）の形式によって表される。推定１次利用者ＩＤは、推定された１次利用者の識別情報である。等高線の中心位置は、緯度と、経度と、高さとを含む。緯度および経度は、ＧＰＳによって検出され、高さは、高度計によって検出される。そして、高さは、等高線の中心位置の地面からの高さを示す。等高線の電力値は、作成された等高線上の電力値である。等高線の半径は、作成された等高線が円形状である場合の円の半径であり、受信電力の重心点（ｘＧ，ｙＧ）から等高線までの距離である。
【０２６８】
時間情報は、作成手段１６が受信電力の等高線を作成した時刻を特定するのに用いられる。周波数帯域ＩＤは、作成手段１６が、同じ周波数帯域の受信電力に基づいて受信電力の等高線を作成したことを示すのに用いられる。サブキャリアＮｏ．は、作成手段１６が受信電力の等高線を作成するときの元になった受信電力を有する電波が伝送されたサブキャリアを判定するのに用いられる。波源ＩＤは、作成手段１６が作成した受信電力の等高線の波源を特定するのに用いられる。
【０２６９】
周波数帯域ＩＤは、抽出手段１３が、ｋ個の端末装置が共用する共用周波数を識別するのに用いられる。等高線の中心位置を示す位置情報は、ｋ個の端末装置から受信したｋ個のモニター情報に含まれるｋ個の受信電力とｋ個の端末装置のｋ個の位置情報とに基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を処理手段１４が演算した結果である受信電力の重心点（ｘＧ，ｙＧ）からなる。等高線の半径は、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点（ｘＧ，ｙＧ）から最も遠い位置に存在する端末装置と受信電力の重心点（ｘＧ，ｙＧ）との距離を処理手段１４が演算した結果である最大距離ＬＭＡＸからなる。等高線の電力値は、周波数帯域ＩＤ（共用周波数）を有する電波の電力値であり、作成手段１６によって検出された電力値である。
【０２７０】
好ましくは、共用周波数が既に割り当てられている１次利用者の波源の識別情報は、作成手段１６によって、共用周波数、等高線の中心位置を示す位置情報、等高線の電力値および等高線の半径に更に対応付けられた構造を有する。
【０２７１】
データ構造Ｄ＿ＳＴＲ２は、推定１次利用者ＩＤ、時間情報、周波数帯域ＩＤ、サブキャリアＮｏ．、波源ＩＤ、等高線の中心位置、等高線の電力値および等高線の半径の全てを含む必要はなく、周波数帯域ＩＤ、等高線の中心位置を示す位置情報、等高線の電力値および等高線の半径を少なくとも含み、周波数帯域ＩＤ、等高線の中心位置を示す位置情報、等高線の電力値および等高線が相互に対応付けられた構造を有していればよい。データ構造Ｄ＿ＳＴＲ２が周波数帯域ＩＤ、等高線の中心位置を示す位置情報、等高線の電力値および等高線を少なくとも含んでいれば、作成手段１６によって作成された受信電力の等高線を表すことができるからである。
【０２７２】
データ構造Ｄ＿ＳＴＲ２は、電力推定装置１（またはプログラムＰｒｏｇ＿Ａ，Ｐｒｏｇ＿Ｂ，Ｐｒｏｇ＿Ｃのいずれか）が上述した方法によって作成した受信電力の等高線を表すのに用いられる。データ構造Ｄ＿ＳＴＲ２は、コンピュータに入力され、ＣＰＵによって記憶手段に記憶される。そして、データ構造Ｄ＿ＳＴＲ２は、ＣＰＵによって表示手段に表示され、または紙媒体等に印刷される。また、データ構造Ｄ＿ＳＴＲ２は、画像ソフトを実行したＣＰＵによって受信電力の等高線を示す画像を表示手段に表示するのに用いられてもよい。更に、データ構造Ｄ＿ＳＴＲ２は、受信電力の等高線から所望の距離だけ離れた位置における電波の電力値をＣＰＵが演算するのに用いられてもよい。
【０２７３】
上記においては、電力推定装置１は、複数の端末装置２において検出された複数のモニター情報ＭＮＴに基づいて任意の高さにおける受信電力の等高線ＣＴＲを作成し、その作成した受信電力の等高線ＣＴＲに基づいて各場所における受信電力を推定すると説明したが、この発明の実施の形態においては、これに限らず、電力推定装置１は、上述した方法によって、複数のモニター情報ＭＮＴに基づいて任意の高さにおける受信電力の等高線ＣＴＲを作成するだけでもよい。異なる電力値Ｐを有する受信電力の複数の等高線ＣＴＲを作成すれば、その作成された受信電力の複数の等高線ＣＴＲに基づいて、最も大きい電力値Ｐを有する受信電力の等高線ＣＴＲよりも波源Ｓ側の領域ＲＥＧ１、および電力値Ｐが異なる２つの受信電力の等高線ＣＴＲの間の領域ＲＥＧ２における受信電力が分かるので、これらの領域ＲＥＧ１，ＲＥＧ２における受信電力を推定できるからである。
【０２７４】
従って、この発明の実施の形態による電力推定装置は、電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線を作成する電力推定装置であって、複数の端末装置から送信され、かつ、各々が端末装置の位置を示す位置情報と端末装置における電波の受信電力とを含む複数のモニター情報から、受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する抽出手段と、抽出手段によって抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点から最も遠い位置に存在する端末装置と受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する処理手段と、処理手段によって求められたｍ個の受信電力の重心点およびｍ個の最大距離と、ｍ個の高さとに基づいて受信電力の重心点と高さとの関係を示す第１の関係式と、受信電力の複数の電力値に対応付けられ、かつ、各々が最大距離と前記高さとの関係を示す複数の第２の関係式とを演算する演算手段と、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、複数の第２の関係式のうち、作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線を作成する作成手段とを備えていればよい。
【０２７５】
また、この発明の実施の形態による無線通信システムは、任意の高さにおいて受信電力の等高線ＣＴＲを作成する電力推定装置を備えていればよい。
【０２７６】
更に、この発明の実施の形態によるプログラムは、電波の受信電力を推定する対象領域において無線通信に用いられる複数の周波数のうち、複数の端末装置が無線通信に共用する共用周波数を有する電波の受信電力の等高線の作成をコンピュータに実行させるためのプログラムであって、抽出手段が、複数の端末装置から送信され、かつ、各々が端末装置の位置を示す位置情報と端末装置における電波の受信電力とを含む複数のモニター情報から、受信電力の最大値を含むｋ（ｋは、２以上の整数）個のモニター情報を抽出する第１のステップと、処理手段が、第１のステップにおいて抽出されたｋ個のモニター情報に基づいて受信電力によって重み付けされたｋ個の端末装置の位置の平均を受信電力の重心点として求めるとともに、ｋ個の受信電力のうち、作成すべき受信電力の等高線における電力値以上の受信電力を有し、かつ、受信電力の重心点から最も遠い位置に存在する端末装置と受信電力の重心点との距離を最大距離として求める第１の処理をｍ（ｍは、２以上の整数）個の高さについて実行する第２のステップと、演算手段が、第２のステップにおいて求められたｍ個の受信電力の重心点およびｍ個の最大距離と、ｍ個の高さとに基づいて受信電力の重心点と高さとの関係を示す第１の関係式と、受信電力の複数の電力値に対応付けられ、かつ、各々が最大距離と高さとの関係を示す複数の第２の関係式とを演算する第３のステップと、作成手段が、第１の関係式に基づいて任意の高さにおける受信電力の重心点を求め、複数の第２の関係式のうち、作成すべき受信電力の等高線における電力値に対応する第２の関係式に基づいて任意の高さにおける最大距離を求め、その求めた受信電力の重心点および最大距離を用いて、作成すべき受信電力の等高線を作成する第４のステップとをコンピュータに実行させればよい。
【０２７７】
今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施の形態の説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
【産業上の利用可能性】
【０２７８】
この発明は、電力推定装置、それを備えた無線通信システム、コンピュータに実行させるためのプログラムおよびデータ構造に適用される。
【符号の説明】
【０２７９】
１  電力推定装置、２  端末装置、３  基地局、４  有線ケーブル、１０  無線通信システム、１１  受信手段、１２  記録手段、１３  抽出手段、１４  処理手段、１５  演算手段、１６  作成手段、１７  推定手段。
【図１】
2019197957.tif 000010
【図２】
2019197957.tif 000011
【図３】
2019197957.tif 000012
【図４】
2019197957.tif 000013
【図５】
2019197957.tif 000014
【図６】
2019197957.tif 000015
【図７】
2019197957.tif 000016
【図８】
2019197957.tif 000017
【図９】
2019197957.tif 000018
【図１０】
2019197957.tif 000019
【図１１】
2019197957.tif 000020
【図１２】
2019197957.tif 000021
【図１３】
2019197957.tif 000022
【図１４】
2019197957.tif 000023
【図１５】
2019197957.tif 000024
【図１６】
2019197957.tif 000025
【図１７】
2019197957.tif 000026
【図１８】
2019197957.tif 000027
【図１９】
2019197957.tif 000028
【図２０】
2019197957.tif 000029
【図２１】
2019197957.tif 000030
【図２２】
2019197957.tif 000031
【図２３】
2019197957.tif 000032
【図２４】
2019197957.tif 000033
【図２５】
2019197957.tif 000034