(57)【要約】【課題】衛星からの干渉電力をより容易にかつより正確に推定する。【解決手段】干渉電力推定装置は、球体の上空を周回する与干渉局及び球体に設置された被干渉局の位置の情報に基づいて、被干渉局から与干渉局に至る局間ベクトルを算出する。干渉電力推定装置は、被干渉局の位置から与干渉局の位置への見通しがあると判定した場合に、被干渉局のアンテナの方向と被干渉局の位置の情報とに基づいて、被干渉局の位置からのアンテナの方向を表すアンテナ方向ベクトルを算出する。干渉電力推定装置は、アンテナ方向ベクトルと局間ベクトルとの角度の差に基づいて算出されるアンテナ指向性減衰量、及び、被干渉局と与干渉局との間の距離に基づいて算出される伝搬損失の総和を算出する。干渉電力推定装置は、与干渉局がとり得る複数の位置の情報それぞれについて算出した総和の最小値に基づいて、与干渉局からの干渉電力を算出する。【選択図】図３
2019054346.tif 000002
【特許請求の範囲】
【請求項１】
球体の上空を周回する与干渉局の位置の情報と、前記球体に設置された被干渉局の位置の情報とに基づいて、前記被干渉局から前記与干渉局に至るベクトルである局間ベクトルを算出する局間ベクトル算出ステップと、
前記被干渉局の位置から前記与干渉局の位置への見通しがあるか否かを判定する見通し判定ステップと、
前記見通し判定ステップにおいて見通しありと判定した場合に、前記被干渉局のアンテナの方向と、前記被干渉局の位置の情報とに基づいて、前記被干渉局の位置からの前記アンテナの方向を表すベクトルであるアンテナ方向ベクトルを算出するアンテナ方向ベクトル算出ステップと、
前記アンテナ方向ベクトル算出ステップにおいて算出された前記アンテナ方向ベクトルと前記局間ベクトルとの角度の差に基づいて算出されるアンテナ指向性減衰量、及び、前記被干渉局と前記与干渉局との間の距離に基づいて算出される伝搬損失の総和を算出する総和算出ステップと、
前記局間ベクトル算出ステップ、前記見通し判定ステップ、前記アンテナ方向ベクトル算出ステップ及び前記総和算出ステップを、前記与干渉局がとり得る複数の位置の情報それぞれについて行う繰り返しステップと、
前記与干渉局がとり得る前記複数の位置の情報それぞれについて算出された前記総和の最小値を探索する探索ステップと、
前記探索ステップにおける探索により得られた前記総和の最小値に基づいて、前記与干渉局からの干渉電力を算出する干渉電力算出ステップと、
を有する干渉電力推定方法。
【請求項２】
前記与干渉局が存在し得る緯度、経度及び標高と、前記被干渉局の緯度、経度及び標高との少なくとも一部に基づいて、前記干渉電力の算出対象となる前記与干渉局の位置の情報を複数取得する位置取得ステップをさらに備え、
前記繰り返しステップにおいては、前記局間ベクトル算出ステップ、前記見通し判定ステップ、前記アンテナ方向ベクトル算出ステップ及び前記総和算出ステップを、前記位置取得ステップにおいて取得された複数の前記位置の情報それぞれについて行う、
請求項１に記載の干渉電力推定方法。
【請求項３】
前記干渉電力算出ステップにおいては、前記総和の最小値と、前記与干渉局及び前記被干渉局のアンテナ利得及び給電の損失とに基づいて前記与干渉局からの干渉電力を算出する、
請求項１又は請求項２に記載の干渉電力推定方法。
【請求項４】
前記見通し判定ステップにおいては、前記被干渉局から前記与干渉局へ向かう直線と前記球体を表す球の方程式とから求められる判別式の値に基づいて見通しがあるか否かを判定する、
請求項１から請求項３のいずれか一項に記載の干渉電力推定方法。
【請求項５】
球体の上空を周回する与干渉局の位置の情報と、前記球体に設置された被干渉局の位置の情報とに基づいて、前記被干渉局から前記与干渉局に至るベクトルである局間ベクトルを算出する処理を、前記与干渉局がとり得る複数の位置の情報それぞれについて行う局間ベクトル算出部と、
前記被干渉局の位置から前記与干渉局の位置への見通しがあるか否かを判定する処理を、前記与干渉局がとり得る前記複数の位置の情報それぞれについて行う見通し判定部と、
前記見通し判定部が見通しありと判定した場合に、前記被干渉局のアンテナの方向と、前記被干渉局の位置の情報とに基づいて、前記被干渉局の位置からの前記アンテナの方向を表すベクトルであるアンテナ方向ベクトルを算出する処理を、前記与干渉局がとり得る前記複数の位置の情報それぞれについて行うアンテナ方向ベクトル算出部と、
前記アンテナ方向ベクトル算出部により算出された前記アンテナ方向ベクトルと前記局間ベクトルとの角度の差に基づいて算出されるアンテナ指向性減衰量、及び、前記被干渉局と前記与干渉局との間の距離に基づいて算出される伝搬損失の総和を算出する処理を、前記与干渉局がとり得る前記複数の位置の情報それぞれについて行う総和算出部と、
前記与干渉局がとり得る前記複数の位置の情報それぞれについて前記総和算出部が算出した総和の最小値を探索する探索部と、
前記探索部の探索により得られた前記総和の最小値に基づいて、前記与干渉局からの干渉電力を算出する干渉電力算出部と、
を備える干渉電力推定装置。
【請求項６】
コンピュータに、
請求項１から請求項４のいずれか一項に記載の干渉電力推定方法を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、干渉電力推定方法、干渉電力推定装置及びプログラムに関する。
【背景技術】
【０００２】
近年、月面や小惑星の探査、気象や温室効果ガスの観測等、様々な目的での衛星打ち上げが活発化している。そのため、衛星通信に利用されている数ＧＨｚから数十ＧＨｚといったマイクロ波帯において、他の通信システムと共用されている周波数では、事業者間での干渉調整が増加してきている。一般的に、干渉調整では、被干渉局のシステムにおける干渉電力の許容値と与干渉局から到来する干渉電力との大小関係を比較し、共用可否を判断する。衛星が正常な状態の場合については、例えば第４世代移動通信システムに対する固定衛星業務システムからの干渉調整では、無線通信規則／電波法施工規則中の地表面電力束密度（ＰＦＤ）制限によって共用可否が検討されてきた（例えば、非特許文献１参照）。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】“平成24年度  情報通信審議会  情報通信技術分科会  携帯電話等高度化委員会報告（案）”，総務省，2012年，p.86-89
【発明の概要】
【発明が解決しようとする課題】
【０００４】
月面や小惑星の探査に用いられる衛星は地球から遠く離れた目的地から地球へ向かって情報を伝えるため、地球近傍よりも利得の高いアンテナおよび高い送信電力で電波を放射する必要がある。このような衛星では、異常状態になった場合に地球近傍にもかかわらず利得の高いアンテナおよび高い送信電力で電波が放射される可能性があるため、異常状態の検討が必要となる。簡易に見積ろうとすると、被干渉局に指向性アンテナを使用している場合は、指向性減衰量が最も小さくなるアンテナ正対方向に対して干渉電力を検討することが考えられる。しかし、干渉電力は指向性減衰量だけでなく、２局間の距離に基づく伝搬損失も考慮する必要があるので、以下のように干渉電力を小さく見積る可能性がある。
【０００５】
図１５は、衛星が与干渉局Ａｉ、別の無線通信を行う地上の固定局が被干渉局Ｂである場合に、衛星の位置（経度）に対する伝搬損失及び指向性減衰の変化を示す図である。同図のグラフでは、与干渉局Ａｉとなる衛星軌道の条件を緯度４２°で計算している。ここで、図１５での条件として、与干渉局Ａｉの標高（高度）を僅か１０ｋｍとした。このような条件は、後述する課題が顕著に数値として示せる例を挙げるために設定した。この高度では、通常の衛星として地球を周回できない軌道である。衛星打上げ途中の高度、あるいは地球周回軌道へ投入できず大気圏へ突入した場合などに当る。またあるいは、この後で述べる本件で提案する干渉計算方法を衛星だけでなく、航空機の航空路等へ拡張させて適用するケースにもなる。
【０００６】
再び、表示グラフの内容に関する説明に戻ると、図１５での符号Ｌ９１は指向性減衰を、符号Ｌ９２は伝搬損失を、符号Ｌ９３は伝搬損失と指向性減衰の和を示している。同図によれば、指向性減衰（符号Ｌ９１）が最も低い場合に、衛星がある位置は経度１３０°となることが分かる。しかし、伝搬損失（符号Ｌ９２）も考慮して両者を足し合わせると（符号Ｌ９３）、経度１４０°での与干渉局Ａｉが被干渉局Ｂへ最も影響を与えている。つまり、先に挙げた単に指向性減衰のみが最小となる位置で見積もる場合は、伝搬損失を含めて精確に見積もる場合と比較して、干渉電力を小さく見積もる可能性がある。本条件においてはその差は８ｄＢ程度となり、大きな差が生じる。これに対して、精密に見積ろうとすると、衛星の実軌跡上で、かつ、全てのポイントに対して検討を行う必要があり、膨大な計算量となる。
【０００７】
上記事情に鑑み、本発明は、衛星からの干渉電力をより容易にかつより正確に推定することができる干渉電力推定方法、干渉電力推定装置及びプログラムを提供することを目的としている。
【課題を解決するための手段】
【０００８】
本発明の一態様は、球体の上空を周回する与干渉局の位置の情報と、前記球体に設置された被干渉局の位置の情報とに基づいて、前記被干渉局から前記与干渉局に至るベクトルである局間ベクトルを算出する局間ベクトル算出ステップと、前記被干渉局の位置から前記与干渉局の位置への見通しがあるか否かを判定する見通し判定ステップと、前記見通し判定ステップにおいて見通しありと判定した場合に、前記被干渉局のアンテナの方向と、前記被干渉局の位置の情報とに基づいて、前記被干渉局の位置からの前記アンテナの方向を表すベクトルであるアンテナ方向ベクトルを算出するアンテナ方向ベクトル算出ステップと、前記アンテナ方向ベクトル算出ステップにおいて算出された前記アンテナ方向ベクトルと前記局間ベクトルとの角度の差に基づいて算出されるアンテナ指向性減衰量、及び、前記被干渉局と前記与干渉局との間の距離に基づいて算出される伝搬損失の総和を算出する総和算出ステップと、前記局間ベクトル算出ステップ、前記見通し判定ステップ、前記アンテナ方向ベクトル算出ステップ及び前記総和算出ステップを、前記与干渉局がとり得る複数の位置の情報それぞれについて行う繰り返しステップと、前記与干渉局がとり得る前記複数の位置の情報それぞれについて算出された前記総和の最小値を探索する探索ステップと、前記探索ステップにおける探索により得られた前記総和の最小値に基づいて、前記与干渉局からの干渉電力を算出する干渉電力算出ステップと、を有する干渉電力推定方法である。
【０００９】
また、本発明の一態様は、上述の干渉電力推定方法であって、前記与干渉局が存在し得る緯度、経度及び標高と、前記被干渉局の緯度、経度及び標高との少なくとも一部に基づいて、前記干渉電力の算出対象となる前記与干渉局の位置の情報を複数取得する位置取得ステップをさらに備え、前記繰り返しステップにおいては、前記局間ベクトル算出ステップ、前記見通し判定ステップ、前記アンテナ方向ベクトル算出ステップ及び前記総和算出ステップを、前記位置取得ステップにおいて取得された複数の前記位置の情報それぞれについて行う。
【００１０】
また、本発明の一態様は、上述の干渉電力推定方法であって、前記干渉電力算出ステップにおいては、前記総和の最小値と、前記与干渉局及び前記被干渉局のアンテナ利得及び給電の損失とに基づいて前記与干渉局からの干渉電力を算出する。
【００１１】
また、本発明の一態様は、上述の干渉電力推定方法であって、前記見通し判定ステップにおいては、前記被干渉局から前記与干渉局へ向かう直線と前記球体を表す球の方程式とから求められる判別式の値に基づいて見通しがあるか否かを判定する。
【００１２】
本発明の一態様は、球体の上空を周回する与干渉局の位置の情報と、前記球体に設置された被干渉局の位置の情報とに基づいて、前記被干渉局から前記与干渉局に至るベクトルである局間ベクトルを算出する処理を、前記与干渉局がとり得る複数の位置の情報それぞれについて行う局間ベクトル算出部と、前記被干渉局の位置から前記与干渉局の位置への見通しがあるか否かを判定する処理を、前記与干渉局がとり得る前記複数の位置の情報それぞれについて行う見通し判定部と、前記見通し判定部が見通しありと判定した場合に、前記被干渉局のアンテナの方向と、前記被干渉局の位置の情報とに基づいて、前記被干渉局の位置からの前記アンテナの方向を表すベクトルであるアンテナ方向ベクトルを算出する処理を、前記与干渉局がとり得る前記複数の位置の情報それぞれについて行うアンテナ方向ベクトル算出部と、前記アンテナ方向ベクトル算出部により算出された前記アンテナ方向ベクトルと前記局間ベクトルとの角度の差に基づいて算出されるアンテナ指向性減衰量、及び、前記被干渉局と前記与干渉局との間の距離に基づいて算出される伝搬損失の総和を算出する処理を、前記与干渉局がとり得る前記複数の位置の情報それぞれについて行う総和算出部と、前記与干渉局がとり得る前記複数の位置の情報それぞれについて前記総和算出部が算出した総和の最小値を探索する探索部と、前記探索部の探索により得られた前記総和の最小値に基づいて、前記与干渉局からの干渉電力を算出する干渉電力算出部と、を備える干渉電力推定装置である。
【００１３】
本発明の一態様は、コンピュータに、上述の干渉電力推定方法を実行させるプログラムである。
【発明の効果】
【００１４】
本発明により、衛星からの干渉電力をより容易にかつより正確に推定することが可能となる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施形態による探査衛星の航行軌道と、干渉推定のための簡易な形状の衛星軌道とを示す図である。
【図２】同実施形態による干渉電力推定装置の機能ブロック図である。
【図３】同実施形態による干渉電力推定装置が実行する干渉電力推定の簡略フローを示す図である。
【図４】同実施形態による干渉電力推定装置が実行する干渉電力推定の詳細フローを示す図である。
【図５】同実施形態による干渉電力推定装置が実行する干渉電力推定の詳細フローを示す図である。
【図６】同実施形態による干渉電力推定装置が実行する干渉電力推定の詳細フローを示す図である。
【図７】同実施形態による被干渉局および与干渉局の位置関係の例を示す図である。
【図８】同実施形態による地球の中心Ｏと主要なベクトルを示す図である。
【図９】同実施形態による判別式の値と見通しの有無との関係を示す図である。
【図１０】同実施形態による俯瞰図で各種の位置とベクトルを示す図である。
【図１１】同実施形態による点Ｔの算出方法をまとめた図である。
【図１２】同実施形態による被干渉局から与干渉局へ向かうベクトルと被干渉局のアンテナ方向を示すベクトルの例を示す図である。
【図１３】同実施形態による伝搬損失および指向性減衰量の算出例を示す図である。
【図１４】同実施形態による伝搬損失および指向性減衰量の算出例を示す図である。
【図１５】簡易見積による干渉電力の推定誤差を示す図である。
【発明を実施するための形態】
【００１６】
以下、図面を参照しながら本発明の実施形態を詳細に説明する。以下では、ベクトルを表すため上部に→を付した文字列（例：Ｘ）を、「（Ｘ）→」のように記載する。
【００１７】
図１は、探査衛星の航行軌道と、干渉推定のための簡易な形状の衛星軌道とを示す図である。同図には、探査衛星が地球からロケットで打上げられ、探査ターゲットの月や惑星を周回する軌道となるまでを示している。まず、衛星打上げ後１〜２週間での地球周回軌道がある。この地球周回軌道は、衛星位置Ａ１，Ａ２，Ａ３，Ａ４，…，Ａｋを含む。この地球周回軌道の次に、衛星は、月などの探査ターゲットへ向けた航行軌道（衛星位置Ａｌ）となり、最後には月や惑星などの探査ターゲットを周回する軌道（衛星位置Ａｍ）になる。最初の地球周回軌道は、地軸に対して斜めの断面で楕円の場合（例えば、衛星位置Ａｋ）もあるなど、実際は不規則である。ここでは、簡単に干渉量を推定するため、衛星の軌道を以下のようにモデル化する。すなわち、緯度θａｉ、標高（高度）ｈａｉとし、地球の半径Ｒと併せて地球の中心から距離Ｒ＋ｈａｉで周回する簡易な円状モデル（同図における衛星位置Ａ１，Ａ３）として検討する。なお、楕円状モデルとして検討してもよい。
【００１８】
図２は、干渉電力推定装置１の機能ブロック図である。干渉電力推定装置１は、入力部２、繰り返し処理部３、局間ベクトル算出部４、見通し判定部５、アンテナ方向ベクトル算出部６、総和算出部７、探索部８及び干渉電力算出部９を備える。
【００１９】
入力部２は、円状モデル上の与干渉局Ａｉ（ｉ＝１〜ｎ）及び被干渉局Ｂに関する情報の入力を受ける。本実施形態では、与干渉局Ａｉは衛星局であり、被干渉局Ｂは地球上の局である。入力される情報は、与干渉局Ａｉ及び被干渉局Ｂが存在しうる位置、アンテナ利得及び給電等の損失、与干渉局Ａｉの送信電力（又は送信電力密度）、被干渉局Ｂのアンテナ方向α及びアンテナ指向性パターン等である。
【００２０】
繰り返し処理部３は、繰り返し算出の指示を出力し、円状モデル上のｎ個の異なる位置の与干渉局Ａｉ（ｉ＝１〜ｎ）それぞれについて、局間ベクトル算出部４、見通し判定部５、アンテナ方向ベクトル算出部６及び総和算出部７により伝搬損失および指向性減衰量の総和を算出する処理を行うよう制御する。与干渉局Ａｉの位置は、入力部２により入力された情報から得られる緯度、経度及び標高の組合せで表される。繰り返し処理部３は、全ての与干渉局Ａｉについて総和の算出を終了すると、繰り返し算出の終了を探索部８に通知する。
【００２１】
局間ベクトル算出部４は、局間ベクトルであるベクトル（Ｃｉ）→を算出する。局間ベクトルは、与干渉局Ａｉから被干渉局Ｂへ向かうベクトルである。見通し判定部５は、被干渉局Ｂの位置から与干渉局Ａｉの位置への見通しの有無を判定する。アンテナ方向ベクトル算出部６は、見通し判定部５が見通しありと判定した場合、被干渉局Ｂのアンテナ方向を表すベクトル（Ｔ１）→を算出する。総和算出部７は、与干渉局Ａｉの伝搬損失と指向性減衰量の総和Ｌｉを算出する。伝搬損失は、与干渉局Ａｉ及び被干渉局Ｂの距離｜（Ｃｉ）→｜に基づいて算出される。指向性減衰量は、与干渉局Ａｉから被干渉局Ｂへ向かうベクトル（Ｃｉ）→及び被干渉局Ｂのアンテナ方向（Ｔ１）→がなす角γと、被干渉局Ｂの指向性パターンとに基づいて算出される。
また、見通し判定部５が、もし被干渉局Ｂの位置から与干渉局Ａｉの位置への見通しはないと判定した場合、見通しなしの判定が繰り返し処理部３へ返される。この見通しなしの判定となる場合には、アンテナ方向ベクトル算出部６と総和算出部７はそれぞれの算出の処理を実行しない。
【００２２】
探索部８は、繰り返し処理部３から繰り返し算出の終了の通知を受け、総和算出部７が算出した与干渉局Ａ１〜Ａｎそれぞれの伝搬損失と指向性減衰量の総和Ｌ１〜Ｌｎのうち、最小の総和Ｌｍｉｎを探索する。干渉電力算出部９は、探索部８が探索した最小の総和Ｌｍｉｎと、入力部２により入力された与干渉局Ａｉの送信電力（又は送信電力密度）、与干渉局Ａｉ及び被干渉局Ｂのアンテナ利得及び給電等の損失等に基づいて干渉電力を算出する。
【００２３】
図３は、干渉電力推定装置１が実行する干渉電力推定の簡略フローを示す図である。まず、入力部２は、与干渉局Ａｉ及び被干渉局Ｂに関する情報を入力する（ステップＳ１０５）。繰り返し処理部３は、入力部２により入力された与干渉局Ａｉ、被干渉局Ｂそれぞれが存在し得る位置の情報に基づいて、与干渉局Ａｉの円状の軌道における位置を表す経度、経度及び標高の組合せをｎ個生成する。繰り返し処理部３は、それらの組のそれぞれについて、伝搬損失および指向性減衰量の総和を算出するため、ｉを１からｎまで１ずつ増加させて、ステップＳ１１５〜Ｓ１３０までの処理を繰り返すよう制御する（ステップＳ１１０）。
【００２４】
局間ベクトル算出部４は、与干渉局Ａｉと被干渉局Ｂの位置に基づいて被干渉局Ｂから与干渉局Ａｉへ向かうベクトル（Ｃｉ）→を算出する（ステップＳ１２０）。見通し判定部５は、ベクトル（Ｃｉ）→に基づいて、被干渉局Ｂから与干渉局Ａｉへの見通しがあるか否かを判定する（ステップＳ１２０）。見通し判定部５が見通しありと判定した場合、アンテナ方向ベクトル算出部６は、被干渉局Ｂのアンテナ方向を示すベクトル（Ｔ１）→を求める（ステップＳ１２５）。詳細な算出方法については後述する。
【００２５】
総和算出部７は、与干渉局Ａｉ及び被干渉局Ｂの距離｜（Ｃｉ）→｜に基づいて伝搬損失を算出する。さらに、総和算出部７は、与干渉局Ａｉから被干渉局Ｂへ向かうベクトル（Ｃｉ）→及び被干渉局Ｂのアンテナ方向（Ｔ１）→がなす角γと、被干渉局Ｂの指向性パターンとに基づいて指向性減衰量を算出する。総和算出部７は、伝搬損失と指向性減衰量の総和Ｌｉを算出する（ステップＳ１３０）。
【００２６】
繰り返し処理部３は、ｉが繰り返し回数ｎに達したと判定すると、繰り返し算出の終了を探索部８に通知する。探索部８は、全ての与干渉局Ａｉそれぞれについて算出した伝搬損失及び指向性減衰量の総和Ｌ１〜Ｌｎのうち最小値のＬｍｉｎを探索する（ステップＳ１３５）。干渉電力算出部９は、ステップＳ１１０で入力された各値と、ステップＳ１３５において探索部８が求めた伝搬損失及び指向性減衰量の総和の最小値Ｌｍｉｎとを用いて、被干渉局Ｂが受信する干渉電力Ｐｒを算出する（ステップＳ１４０）。
【００２７】
図４〜図６は、図３に示す干渉電力推定処理の詳細フローを示す図である。図３のステップＳ１０５（両局の情報入力）は図４のステップＳ２０５〜ステップＳ２１０に、図３のステップＳ１１０（繰り返し）は図４のステップＳ２１５に、図３のステップＳ１１５（ベクトル（Ｃｉ）→算出）は図４のステップＳ２２０に相当する。また、図３のステップＳ１２０（見通し判定）は図４のステップＳ２２５〜図５のステップＳ３０５に、図３のステップＳ１２５（ベクトル（Ｔ１）→算出）は図５のステップＳ３１０〜図６のステップＳ４２０に、図３のステップＳ１３０（総和Ｌｉ算出）は図６のステップＳ４２５〜ステップＳ４３５に相当する。そして、図３のステップＳ１３５（最小値Ｌｍｉｎ探索）は図６のステップＳ４４０に、図３のステップＳ１４０（干渉電力Ｐｒ算出）は図６のステップＳ４４５に相当する。
【００２８】
まず、図４において、衛星局を与干渉局Ａｉ（ｉは１以上の整数）、地球上の局を被干渉局Ｂと設定する（ステップＳ２０５）。入力部２は、与干渉局Ａｉ及び被干渉局Ｂの情報を入力する（ステップＳ２１０）。入力される情報は、与干渉局Ａｉのアンテナ利得Ｇｔ、給電等の損失Ｌｔ及び送信電力（又は送信電力密度）、与干渉局Ａｉが存在し得る緯度θａｉ、経度φａｉ及び標高ｈａｉ、被干渉局Ｂのアンテナ利得Ｇｒ及び給電等の損失Ｌｒ、被干渉局Ｂが存在し得る緯度θｂ、経度φｂ及び標高ｈｂ、被干渉局Ｂのアンテナ方向α及びアンテナ指向性パターン等である。アンテナ方向αは、北を０°とし、時計回りの方向の角度が大きいほど、値が増加する。
【００２９】
図７は、与干渉局Ａｉおよび被干渉局Ｂの位置関係の例を示す図である。同図では、与干渉局Ａｉの緯度θａｉ＝−３０〜３０°、経度φａｉ＝０〜３６０°、標高ｈａｉ＝５５０ｋｍであり、被干渉局Ｂの緯度θｂ＝４２°、経度φｂ＝１４１°、標高ｈｂ＝８９１ｍ、アンテナ方向α＝２６４°（北を０°として、時計回りとした角度）である。
【００３０】
図４において、繰り返し処理部３は、与干渉局Ａｉ（ｉ＝１〜ｎ）の位置を表す緯度、経度及び標高のｎ個の組全てについて、伝搬損失および指向性減衰量の総和を算出するため、ｉを１からｎまで１ずつ増加させて、図４のステップＳ２２０〜図６のステップＳ４３５までの処理を繰り返す（ステップＳ２１５）。図４のステップＳ２２０〜図６のステップＳ４３５までの処理を総和算出処理と記載する。
【００３１】
総和算出処理は、与干渉局Ａｉのパラメータである衛星の軌道に当る緯度θａｉ、経度φａｉ及び標高（高度）ｈａｉについて行われる。例えば、具体例として、与干渉局Ａｉの緯度θａｉが３０°、標高（高度）ｈａｉ≧５５０ｋｍであり、被干渉局Ｂの位置が北海道（緯度θｂ＝４２°、経度φｂ＝１４１°）であることが入力された場合、繰り返し処理部３は、次の様にする。繰り返し処理部３は、入力された被干渉局Ｂの位置を考慮して、与干渉局Ａｉの標高（高度）ｈａｉを取り得る最低標高である５５０ｋｍ、緯度θａｉを３０°とし、経度φａｉについては、８０°（≒被干渉局Ｂの経度φｂ−６０°）〜１８０°（≒被干渉局Ｂの経度φｂ＋４０°）の間を１°間隔で変えた値とする。これにより、経度φａｉを変えながら、総和算出処理により１００箇所の各位置の与干渉局Ａｉについて伝搬損失および指向性減衰量の総和が計算される。この場合、繰り返し処理部３は、ｉを１から１００まで１ずつ変化させる。この例では、与干渉局Ａｉの経度のみを変えているが、入力の条件によっては与干渉局Ａｉの緯度や与干渉局Ａｉの標高（高度）を変えて干渉量の推定を行う。
【００３２】
なお、上述した具体例の両局ＡｉとＢに関する数値（緯度θａｉ＝３０°、標高ｈａｉ≧５５０ｋｍ、緯度θｂ＝４２°、経度φｂ＝１４１°）は、先の図７および後述する計算例（図１３、図１４）で同じ条件を挙げている。また、後述する計算例では、グラフ表示した計算結果の変化状況を考慮し、８０°〜１８０°の間を１°間隔とした経度φａｉの１００箇所（ｉ＝１，２，…，１００）の与干渉局Ａｉに対し、伝搬損失および指向性減衰量の総和の計算を繰り返している。
【００３３】
図４において、局間ベクトル算出部４は、被干渉局Ｂから与干渉局Ａｉへ向かうベクトル（Ｃｉ）→を算出する（ステップＳ２２０）。局間ベクトル算出部４は、地球半径Ｒを用いて与干渉局Ａｉ、被干渉局Ｂそれぞれの緯度経度情報を直交座標である点Ａｉ（Ｘａｉ，Ｙａｉ，Ｚａｉ）および点Ｂ（Ｘｂ，Ｙｂ，Ｚｂ）に変換し、点Ｂから点Ａｉに向かうベクトル（Ｃｉ）→を以下の式（１）のように算出する。
【００３４】
【数１】
2019054346.tif 000003
【００３５】
図８は、地球の中心Ｏと主要なベクトルを示す図である。同図は、地球の中心Ｏと標高ｈｂにある被干渉局Ｂおよび与干渉局Ａｉを通る断面を示している。また、同図には、地球の中心Ｏから被干渉局Ｂに向かうベクトル（Ｂ）→、被干渉局Ｂから与干渉局Ａｉに向かうベクトル（Ｃｉ）→に加え、ベクトル（Ｂ）→の単位ベクトル（ａｂ）→も示している。
【００３６】
見通し判定部５は、被干渉局Ｂから与干渉局Ａｉへ向かう直線と球の方程式から求まる判別式により、被干渉局Ｂから与干渉局Ａｉへの見通しがあるか否かを判定する。見通し判定部５が見通しありと判定した場合、干渉電力推定装置１は、与干渉局Ａｉの総和算出処理を続行し、見通しなしと判定した場合は、次の与干渉局Ａｉ＋１の総和算出処理へ移る。
【００３７】
具体的には、まず、図４において、見通し判定部５は、被干渉局Ｂから与干渉局Ａｉを通る直線上の点を、以下の式（２）に示すベクトル（Ｂ）→とする。ただし、ｔは実数である。
【００３８】
【数２】
2019054346.tif 000004
【００３９】
式（２）に示すこの点が、半径Ｒ（地球半径）の球上に存在すると仮定して球の方程式に代入すると、式（３）となる。
【００４０】
【数３】
2019054346.tif 000005
【００４１】
見通し判定部５は、式（３）をｔについて整理し、以下の式（４）の２次方程式を得る（ステップＳ２２５）。
【００４２】
【数４】
2019054346.tif 000006
【００４３】
図５において、見通し判定部５は、式（４）の方程式から、式（５）に示す判別式Ｄ／４を求める。
【００４４】
【数５】
2019054346.tif 000007
【００４５】
見通し判定部５は、式（５）の値を算出し、算出した判別式Ｄ／４の値に基づいて、被干渉局Ｂから与干渉局Ａｉへの見通しの有無を判定する（ステップＳ３０５）。
【００４６】
図９は、式（５）に示す判別式の値と、被干渉局Ｂから与干渉局Ａｉへの見通しの有無との関係を示す図である。見通し判定部５は、判別式Ｄ／４の値が０より大きい場合、見通しなしと判断する。見通しなしの場合、干渉電力推定装置１は、現在の与干渉局Ａｉについては、図５のステップＳ３１０以降の総和算出処理を行わない。一方、見通し判定部５は、判別式Ｄ／４の値が０以下であれば見通しありと判断する。見通しありの場合、干渉電力推定装置１は、図５のステップＳ３１０以降の総和算出処理を行う。
【００４７】
図１０は、俯瞰図で各種の位置とベクトルを示す図である。同図では、この後で述べる接平面の外周上の点Ｔ、ベクトル（Ｖ）→、被干渉局Ｂのアンテナ方向αを示すベクトル（Ｔ１）→、アンテナ方向αと後述する方向βｍｉｎとの差δを示している。
【００４８】
アンテナ方向ベクトル算出部６は、図５のステップＳ３１０〜図６のステップＳ４２０の処理により、被干渉局Ｂのアンテナ方向を示すベクトル（Ｔ１）→を求める。具体的には、まず半径Ｒ＋Ｚｂの球について点Ｂを通る式（６）の接平面を考える。
【００４９】
【数６】
2019054346.tif 000008
【００５０】
次に、アンテナ方向ベクトル算出部６は、β=０°(北)、９０°(東)、１８０°(南)、２７０°(西)のうち、αとの差の絶対値(|β−α|)が最小となるβを選択してβｍｉｎとするとともに、その差をδ＝βｍｉｎ−αとする（ステップＳ３１０）。ここからは被干渉局Ｂの位置により処理が異なる。
【００５１】
先ず、アンテナ方向ベクトル算出部６は、被干渉局Ｂの位置が北半球であるか南半球であるかを判断する（ステップＳ３１５）。アンテナ方向ベクトル算出部６は、被干渉局Ｂの緯度θｂ＝０〜＋９０°であり、被干渉局Ｂが北半球に位置すると判定した場合、ステップＳ３１０で求めたβｍｉｎの値を判定する（ステップＳ３２０）。
【００５２】
アンテナ方向ベクトル算出部６は、ステップＳ３２０においてβｍｉｎが０°(北)であると判定した場合、Ｘ＝Ｘｔ＝０、Ｙ＝Ｙｔ＝０として、接線の式により点Ｔ（０，０，Ｚｔ）を求める（ステップＳ３２５）。
【００５３】
アンテナ方向ベクトル算出部６は、ステップＳ３２０においてβｍｉｎが９０°(東)であると判定した場合、被干渉局Ｂの経度φｂを判定する（ステップＳ３３０）。アンテナ方向ベクトル算出部６は、ステップＳ３３０において被干渉局Ｂの経度φｂが東経０〜９０°又は西経９０〜１８０°の範囲であると判定した場合、Ｘ＝Ｘｔ＝０、Ｚ＝Ｚｔ＝（Ｒ＋ｈｂ）ｓｉｎθｂとして、接線の式により点Ｔ（０，Ｙｔ，Ｚｔ）を求める（ステップＳ３３５）。
【００５４】
アンテナ方向ベクトル算出部６は、ステップＳ３３０において被干渉局Ｂの経度φｂが東経９０〜１８０°又は西経０〜９０°の範囲であると判定した場合、Ｙ＝Ｙｔ＝０，Ｚ＝Ｚｔ＝（Ｒ＋ｈｂ）ｓｉｎθｂとして、接線の式により点Ｔ（Ｘｔ，０，Ｚｔ）を求める（ステップＳ３４０）。
【００５５】
アンテナ方向ベクトル算出部６は、ステップＳ３２０においてβｍｉｎが１８０°(南)であると判定した場合、Ｙ＝Ｙｔ＝（Ｙｂ／Ｘｂ）Ｘｔ、Ｚｔ＝０として、接線の式により点Ｔ（Ｘｔ，Ｙｔ，０）を求める（ステップＳ３４５）。
【００５６】
アンテナ方向ベクトル算出部６は、ステップＳ３２０においてβｍｉｎが２７０°(西)であると判定した場合、被干渉局Ｂの経度φｂを判定する（ステップＳ３５０）。アンテナ方向ベクトル算出部６は、ステップＳ３５０において被干渉局Ｂの経度φｂが東経９０〜１８０°又は西経０〜９０°の範囲であると判定した場合、Ｘ＝Ｘｔ＝０、Ｚ＝Ｚｔ＝（Ｒ＋ｈｂ）ｓｉｎθｂとして、接線の式により点Ｔ（０，Ｙｔ，Ｚｔ）を求める（ステップＳ３３５）。
【００５７】
アンテナ方向ベクトル算出部６は、ステップＳ３５０において被干渉局Ｂの経度φｂが東経０〜９０°又は西経９０〜１８０°の範囲であると判定した場合、Ｙ＝Ｙｔ＝０，Ｚ＝Ｚｔ＝（Ｒ＋ｈｂ）ｓｉｎθｂとして、接線の式により点Ｔ（Ｘｔ，０，Ｚｔ）を求める（ステップＳ３４０）。
【００５８】
アンテナ方向ベクトル算出部６は、ステップＳ３１５において被干渉局Ｂの緯度θｂ＝０〜−９０°であり、被干渉局Ｂが南半球に位置すると判定した場合、ステップＳ３１０で求めたβｍｉｎの値を判定する（ステップＳ３５５）。
【００５９】
アンテナ方向ベクトル算出部６は、ステップＳ３５５においてβｍｉｎが１８０°(南)であると判定した場合、Ｘ＝Ｘｔ＝０、Ｙ＝Ｙｔ＝０として、接線の式により点Ｔ（０，０，Ｚｔ）を求める（ステップＳ３２５）。
【００６０】
アンテナ方向ベクトル算出部６は、ステップＳ３５５においてβｍｉｎが９０°(東)であると判定した場合、上記のステップＳ３３０以降の処理を行う。
【００６１】
アンテナ方向ベクトル算出部６は、ステップＳ３５５においてβｍｉｎが０°(北)であると判定した場合、Ｙ＝Ｙｔ＝（Ｙｂ／Ｘｂ）Ｘｔ＝０、Ｚｔ＝０として、接線の式により点Ｔ（Ｘｔ，Ｙｔ，Ｚｔ）を求める（ステップＳ３４５）。
【００６２】
アンテナ方向ベクトル算出部６は、ステップＳ３５５においてβｍｉｎが２７０°(西)であると判定した場合、上記のステップＳ３５０以降の処理を行う。
【００６３】
図１１は、上記のステップＳ３１５〜ステップＳ３５５における点Ｔの座標（Ｘｔ，Ｙｔ，Ｚｔ）の算出方法をまとめた図である。図５のステップＳ３１５〜ステップＳ３５５において点Ｔが算出されると、アンテナ方向ベクトル算出部６は、点Ｂから点Ｔに向かうベクトル（Ｔ０）→＝（Ｘｔ−Ｘｂ，Ｙｔ−Ｙｂ，Ｚｔ−Ｚｂ）を求める（ステップＳ３６０）。
【００６４】
さらに、図６において、アンテナ方向ベクトル算出部６は、位置ベクトル（Ｂ）→の単位ベクトル（ａｂ）→を、以下の式（７）により求める（ステップＳ４０５）。
【００６５】
【数７】
2019054346.tif 000009
【００６６】
アンテナ方向ベクトル算出部６は、単位ベクトル（ａｂ）→および（Ｔ０）→に垂直なベクトル（Ｖ）→を、以下の式（８）により求める（ステップＳ４１０）。
【００６７】
【数８】
2019054346.tif 000010
【００６８】
アンテナ方向ベクトル算出部６は、単位ベクトル（ａｂ）→を軸として、位置ベクトル（Ｔ）→を反時計回りにδだけ回転させたベクトル（Ｔ’）→を、以下の式（９）により導出する（ステップＳ４１５）。
【００６９】
【数９】
2019054346.tif 000011
【００７０】
最後に、アンテナ方向ベクトル算出部６は、点Ｂから点Ｔ’に向かうベクトル（Ｔ１）→を、以下の式（１０）により算出する（ステップＳ４２０）。
【００７１】
【数１０】
2019054346.tif 000012
【００７２】
このような衛星からの干渉推定では、被干渉局Ｂの位置の直交座標を点Ｂとしたときに、半径Ｒ＋ｈｂの球上での点Ｂの接平面を考える。仮に、ベクトル（Ｔ１）→を直接的に求めようとした場合、接平面内でα方向に当る直線を考え、その直線上の他の１点を求める必要がある。例えば、直線をＺｔ＝０とすると、Ｘｔ，Ｙｔの両方を求めるため、さらに別の式を考えて連立方程式を解かなければならなくなる。これに対して、βｍｉｎから間接的にベクトル（Ｔ１）→を求める本実施形態の干渉推定の場合は、接平面内で南北および東西を結ぶ直線を考える。その直線上の他の１点を求めるには、本実施形態の図５に示すアルゴリズムの通り処理が容易である。従って、このような理由から本干渉推定ではβｍｉｎを求めてベクトルを回転させる方法を用いている。
【００７３】
次に、総和算出部７は、伝搬損失及び指向性減衰量の総和Ｌｉを算出する。
図１２は、被干渉局Ｂから与干渉局Ａｉへ向かうベクトル（Ｃｉ）→と、被干渉局Ｂのアンテナ方向を示すベクトル（Ｔ１）→の例を示す図である。総和算出部７は、局間ベクトル算出部４が求めたベクトル（Ｃｉ）→の絶対値を距離として伝搬損失を算出する。さらに総和算出部７は、ベクトル（Ｃｉ）→と、アンテナ方向ベクトル算出部６が求めた被干渉局Ｂのアンテナ方向を示すベクトル（Ｔ１）→のなす角γに基づいて指向性減衰量を算出する。総和算出部７は、これらの総和であるＬｉを算出する。
【００７４】
具体的には、図６において、総和算出部７は、ベクトル（Ｃｉ）→とベクトル（Ｔ１）→のなす角γを、以下の式（１１）により算出する。
【００７５】
【数１１】
2019054346.tif 000013
【００７６】
総和算出部７は、ステップＳ２０５において入力されたアンテナ指向性パターンを参照して、角γであるときのアンテナ指向性減衰量を導出する（ステップＳ４２５）。総和算出部７は、被干渉局Ｂと与干渉局Ａｉの距離｜（Ｃｉ）→｜に基づいて伝搬損失を算出する（ステップＳ４３０）。伝搬損失は、例えば自由空間伝搬損失であれば、以下の式（１２）により算出される。
【００７７】
【数１２】
2019054346.tif 000014
【００７８】
総和算出部７は、ステップＳ４３０において求めた伝搬損失と、ステップＳ４２５において求めたアンテナ指向性減衰量との総和Ｌｉを算出する（ステップＳ４３５）。
【００７９】
繰り返し処理部３は、ｉが繰り返し回数ｎに達し、図４のステップＳ２２０〜図６のステップＳ４３５までの総和算出処理をｎ回行ったと判定し、繰り返し処理の終了を探索部８に通知する。探索部８は、全ての伝搬損失及び指向性減衰量の総和Ｌ１〜Ｌｎのうち最小値のＬｍｉｎを探索する（ステップＳ４４０）。
【００８０】
図１３および図１４は、干渉電力推定装置１による伝搬損失および指向性減衰量の算出例を示す図である。図１３（ａ）及び図１４（ａ）は、経度に対する損失の変化を示すグラフであり、図１３（ｂ）及び図１４（ｂ）は、与干渉局Ａｉの緯度θａｉに対する伝搬損失及び指向性減衰量の総和Ｌｉを示す図である。図１３（ａ）は、経度φａｉを、図１３（ｂ）のグラフにおいて総和Ｌｉが最小となる３０°として計算している。一方、図１４（ａ）は、経度φａｉを、図１４（ｂ）のグラフにおいて総和Ｌｉが最小となる２２°として計算している。符号Ｌ１１、Ｌ２１は指向性減衰量を、符号Ｌ１２、Ｌ２２は伝搬損失を、符号Ｌ１３、Ｌ２３は伝搬損失と指向性減衰量の総和Ｌｉを示す。
【００８１】
図６において、干渉電力算出部９は、先に図４のステップＳ２０５、Ｓ２１０（図３のステップＳ１０５に相当）で入力された各値および探索部８が求めた（ステップＳ４４０）伝搬損失及び指向性減衰量の総和の最小値Ｌｍｉｎを用いて、被干渉局Ｂが受信する干渉電力Ｐｒを算出する（ステップＳ４４５）。すなわち、干渉電力Ｐｒは、与干渉局Ａｉの送信電力（もしくは送信電力密度)、アンテナ利得Ｇｔ、給電線等の損失Ｌｔと、被干渉局Ｂのアンテナ利得Ｇｒ、給電線等の損失とＴｒ、算出した伝搬損失と指向性減衰量の総和の最小値Ｌｍｉｎとを、各数値が「ｄＢ」単位である場合は、足し合わせることで求められる。また、この総和の最小値Ｌｍｉｎを導出する過程では、被干渉局Ｂのアンテナ方向およびアンテナ指向性パターンが既に考慮されている。従って、その算出した干渉電力Ｐｒを許容干渉電力と比較し、許容干渉電力より干渉電力Ｐｒが低い場合は、衛星からの干渉電力が許容できると判断できる。
【００８２】
なお、本実施形態では、与干渉局Ａｉのアンテナの向き及び指向性は無視している。これは以下の理由による。与干渉局Ａｉである衛星から被干渉局Ｂへ影響する干渉となる無線信号は、異常状況で自律的に送信される。この衛星の異常状況とは、地球側のコントロールセンタから衛星を制御できてなかったり、衛星の正確な軌道（位置）を見失ったりする状況である。従って、この状況下では、衛星である与干渉局Ａｉのアンテナの向きも不確定であり、干渉の影響量は最悪を想定し、最大のアンテナ利得である正面方向を向いた時を考える。このような理由から、本実施形態の干渉推定において、被干渉局Ｂのアンテナの向き及び指向性は考慮しているが、与干渉局Ａｉについては無視している。
【００８３】
上述した実施形態によれば、衛星の軌道を最低標高と円状もしくは楕円状等の簡易な形状でモデル化する。干渉電力推定装置１は、与干渉局である衛星および被干渉局の緯度経度および標高（これらの１つ或いは複数のパラメータを範囲として与えてもよい）、アンテナ利得、給電線等の損失、与干渉局の送信電力（もしくは送信電力密度）、被干渉局のアンテナ方向およびアンテナ指向性パターン、許容干渉電力等の情報を元に、伝搬損失および指向性減衰量の総和の最小値を求める。伝搬損失は、２局間の距離から算出され、指向性減衰量は、伝搬損失およびアンテナ正面方向からのズレから算出される。干渉電力推定装置１は、伝搬損失および指向性減衰量の総和の最小値を求めることにより、衛星からの干渉電力を容易かつ正確に見積もることができる。
【００８４】
なお、上述した実施形態として、衛星の軌道を最低高度（標高）で円状の簡易なモデルを示したが、円状ではなく楕円状等の数式として与えられれば、同様に干渉電力を見積もることができる。また、衛星以外にも最低高度がより低い軌道（航路）も数式として与えることで航空機にも対応が可能である。
また、上述した実施形態で説明したように、見通しを自動的に判定して、見通しがない場合の干渉算出を省略した。この省略した算出処理のリソース（算出時間や演算能力）を有効な干渉算出に振向けられるので、例えば省略された分を活用して円状の軌道の分割数を１°毎に干渉計算する条件を倍の０．５°毎にする条件でより細かくし、さらに正確な干渉量の見積もりへ繋げられる。
【００８５】
従来の干渉推定は，被干渉局Ｂのアンテナ方向に近い軌道に与干渉局Ａｉがある時を対象に干渉量を算出していた。一方、本実施形態では、与干渉局Ａｉの軌道の位置を変えて干渉量を計算し、減衰が最小になる干渉電力を導く。このような算出は通常困難であるが、本実施形態では、与干渉局Ａｉの位置を変えながら、被干渉局Ｂからみた与干渉局Ａｉの位置を表すベクトル(Ｃｉ）→及び被干渉局Ｂのアンテナ方向（Ｔ１）→を用いて繰り返し干渉電力を算出している。
【００８６】
上述した実施形態によれば、干渉電力推定装置は、局間ベクトル算出ステップと、見通し判定ステップと、アンテナ方向ベクトル算出ステップと、総和算出ステップと、繰り返しステップと、探索ステップと、干渉電力算出ステップとを実行する。例えば、局間ベクトル算出ステップは局間ベクトル算出部４により実行され、見通し判定ステップは見通し判定部５により実行され、アンテナ方向ベクトル算出ステップはアンテナ方向ベクトル算出部６により実行され、総和算出ステップは総和算出部７により実行され、繰り返しステップは繰り返し処理部３により実行され、探索ステップは探索部８により実行され、干渉電力算出ステップは干渉電力算出部９により実行される。
【００８７】
局間ベクトル算出ステップでは、球体の上空を周回する与干渉局の位置の情報と、球体に設置された被干渉局の位置の情報とに基づいて、被干渉局から与干渉局に至るベクトルである局間ベクトルを算出する。例えば、球体は、地球などの天体であり、局間ベクトルは、局間ベクトル（Ｃｉ）→である。アンテナ方向ベクトル算出ステップでは、被干渉局から与干渉局へ向かう直線と球の方程式とから求められる判別式の値に基づいて、被干渉局の位置から与干渉局の位置への見通しがあるか否かを判定する。見通しありと判定された場合、アンテナ方向ベクトル算出ステップでは、被干渉局のアンテナの方向と、被干渉局の位置の情報とに基づいて、被干渉局の位置からのアンテナの方向を表すベクトルであるアンテナ方向ベクトルを算出する。アンテナ方向ベクトルは、例えば、ベクトル（Ｔ１）→である。総和算出ステップでは、アンテナ方向ベクトルと局間ベクトルとの角度の差に基づいて算出されるアンテナ指向性減衰量、及び、被干渉局と与干渉局との間の距離に基づいて算出される伝搬損失の総和を算出する。繰り返しステップでは、局間ベクトル算出ステップ、見通し判定ステップ、アンテナ方向ベクトル算出ステップ及び総和算出ステップを、与干渉局がとり得る複数の位置の情報それぞれについて行う。探索ステップでは、与干渉局がとり得る複数の位置の情報それぞれについて算出された総和の最小値を探索する。干渉電力算出ステップでは、総和の最小値に基づいて、与干渉局からの干渉電力を算出する。
【００８８】
干渉電力推定装置は、与干渉局が存在し得る緯度、経度及び標高と、被干渉局の緯度、経度及び標高との少なくとも一部に基づいて、干渉電力の算出対象となる与干渉局の位置の情報を複数取得する位置取得ステップをさらに実行してもよい。位置取得ステップは、例えば、繰り返し処理部３により実行される。繰り返しステップにおいては、局間ベクトル算出ステップ、見通し判定ステップ、アンテナ方向ベクトル算出ステップ及び総和算出ステップを、位置取得ステップにおいて取得された複数の位置の情報それぞれについて行う。
【００８９】
上述した実施形態における干渉電力推定装置１の機能をコンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【００９０】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【００９１】
１…干渉電力推定装置，  ２…入力部，  ３…繰り返し処理部，  ４…局間ベクトル算出部，  ５…見通し判定部，  ６…アンテナ方向ベクトル算出部，  ７…総和算出部，  ８…探索部，  ９…干渉電力算出部
【図１】
2019054346.tif 000015
【図２】
2019054346.tif 000016
【図３】
2019054346.tif 000017
【図４】
2019054346.tif 000018
【図５】
2019054346.tif 000019
【図６】
2019054346.tif 000020
【図７】
2019054346.tif 000021
【図８】
2019054346.tif 000022
【図９】
2019054346.tif 000023
【図１０】
2019054346.tif 000024
【図１１】
2019054346.tif 000025
【図１２】
2019054346.tif 000026
【図１３】
2019054346.tif 000027
【図１４】
2019054346.tif 000028
【図１５】
2019054346.tif 000029