(57)【要約】    【課題】      複数の計算機を用いた時ｎｉ各計算機間の計算量の偏りを少なくし、高速かつ高精度に電波の伝搬解析を行うこと。    【解決手段】  送信点１０１を中心として離散角度Δθの間隔で擬似光線が放射する方向を求め、送信点１０１を中心として離散角度ΔθａのＮ個の方位に分割して放射グループ１からＮまでを作成する。擬似光線１−１を端末計算機１０４−１に、擬似光線１−２を端末計算機１０４−２に割り当てていき、Ｋ本の擬似光線２０１を端末計算機１０４−１から端末計算機１０４−ｍの順に繰り返し割り当てる。送信点１０１からの擬似光線２０１をＭ台の端末計算機に割り当てる際に、割り当てた端末計算機において以前に割り当てた擬似光線２０１の計算が終了していない場合、計算が終了している次の端末計算機に擬似光線２０１を割り当て、Ｎ個の放射グループ１からＮに含まれる全ての擬似光線２０１の電波伝搬解析を実行する。  【選択図】      図１
2005331391.tif 000002
【特許請求の範囲】
【請求項１】
送信点より複数の擬似光線を放射し、その擬似光線が受信点へ到達するまでの軌跡を追跡することにより、所定の電波伝搬モデルにおける送受信点間の電波伝搬特性を幾何光学的に算出する電波伝搬特性解析方法であって、
前記送信点から擬似光線を放射するときに、放射する方位を複数の方向に分割する分割手順と、
分割された方向に含まれる擬似光線を、一本ごとに順に複数の計算機に個別に割り振る割振り手順と、
前記複数の計算機が、それぞれ、自己に割り振られた擬似光線の軌跡を追跡してその擬似光線の伝搬特性を計算する計算手順と、
を含むことを特徴とする電波伝搬特性解析方法。
【請求項２】
有限差分時間領域法を用いてアンテナ複素指向性を導出する手段を備え、
前記送信点を中心とした所定の離散的角度の方位に対する前記アンテナ複素指向性の振幅があらかじめ設定した閾値以下の場合は、
前記割振り手順において、分割された方位に含まれる擬似光線を計算機に割り振らず、
前記計算手順において、前記擬似光線の伝搬特性の計算を行わない
ことを特徴とする請求項１に記載の電波伝搬特性解析方法。
【請求項３】
送信点より複数の擬似光線を放射し、その擬似光線が受信点へ到達するまでの軌跡を追跡することにより、所定の電波伝搬モデルにおける送受信点間の電波伝搬特性を幾何光学的に算出する電波伝搬特性解析装置であって、
前記送信点から擬似光線を放射するときに放射方位を複数の方向に分割し、分割された方向に含まれる擬似光線を一本ごとに分配する主計算手段と、
前記主計算手段から自己に割り振られた擬似光線を取得し、その擬似光線の軌跡を追跡することによって電波伝搬特性を計算する複数の端末計算手段とを備え、
前記主計算手段が、前記複数の端末計算手段の計算結果に基づいて、所定の電波伝搬モデルにおける送受信点間の電波伝搬特性を幾何光学的に算出することを特徴とする電波伝搬特性解析装置。
【請求項４】
さらに、有限差分時間領域法を用いてアンテナ複素指向性を導出する導出手段を備え、
前記主計算手段が、前記複数の端末計算手段の計算結果と前記導出手段によって導出されたアンテナ複素指向性とに基づいて、前記電波伝搬特性を幾何光学的に算出することを特徴とする請求項３に記載の電波伝搬特性解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、送信点から送信された電波の伝搬特性を解析し、その解析結果を出力することを可能にする電波伝搬特性解析方法及びその方法を実現する電波伝搬特性解析装置に関し、特に、無線局から送信された複数の光線を探索する幾何光学的手法を用いて電波の伝搬特性を解析する電波伝搬特性解析方法と電波伝搬特性解析装置に関する。
【背景技術】
【０００２】
近年広く利用されている携帯電話機は、さまざまな使用状態や電波伝搬環境に適応して指向性や偏波を制御できるようになっている。このような制御に対応できる携帯電話機用のアンテナの開発を行うには、人体への影響やアンテナの微細構造を考慮して計算されたアンテナの全立体角複素指向性と多重波伝搬パラメータとを統合して、実際の多重波伝搬環境における実効的なアンテナ性能を評価できる系（つまり、手段）が必要となる。
【０００３】
このようなアンテナ性能の評価系としては従来、ランダムフィールド法などの実験的手法が用いられているが、この方法では大規模な実験装置や多大な労力が必要となるなどの課題がある。そこで、このような課題に対応するためには計算機によるアンテナ特性評価方法が広く利用されている。このようなアンテナ特性評価方法は、一般的には、無線局から送信された電磁波の電界強度分布を求めるための手段として、電磁界解析手法あるいは幾何光学的計算手法を用いたものがある。例えば、電磁波の電界強度分布を求める手段として、下記の非特許文献１に記載されているＦＤＴＤ（Finite Difference Time Domain：有限差分時間領域）法を用いた室内伝搬解析法や、複数の擬似光線の軌跡を探索する幾何光学的計算手法のレイラウンチング法などがある。レイラウンチング法においては、下記の特許文献１に示されているような、送信点から放射される全ての擬似光線を一台の計算機で計算する基本的な計算手段や、下記の特許文献２に示されているような、複数の計算機を用いて送信点から放射される擬似光線の追跡を行う方法などが知られている。
【特許文献１】特開平８−１９４０２９号公報
【特許文献２】特開２００２−２３２３４８号公報
【非特許文献１】IEICE TRANS. CMMUN.,VOL.E86-B, NO.6 JUNE 2003
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、ＦＤＴＤ法を用いた電波伝搬解析方法では、一般的に、電波の周波数が高い、又は電界強度分布の解析領域が大きいと演算量が多くなり、屋外伝搬モデルなどの大規模解析領域において電波の伝搬解析を行うことが困難となる。また、大規模解析領域を解析する手法として幾何光学的手法を用いた電波伝搬解析手法があるが、この解析手法を一台の計算機で処理すると多大な計算時間を要するなどの問題がある。なお、幾何光学的手法を用いた電波伝搬解析手法を複数の計算機で行う解析手法もあるが、この解析手法では、電波の送信点を中心として放射方向ごとに含まれる擬似光線を個々の計算機に割り振って電波の伝搬解析を行うため、各計算機間に計算量の偏りが生じ計算効率が低下するなどの問題がある。
【０００５】
本発明は上記の事情に鑑みてなされたものであり、複数の計算機を用いても各計算機間の計算量の偏りを少なくすると共に高速かつ高精度に電波の伝搬解析ができるような電波伝搬特性解析方法及び電波伝搬特性解析装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の電波伝搬特性解析方法は、送信点より複数の擬似光線を放射し、その擬似光線が受信点へ到達するまでの軌跡を追跡することにより、所定の電波伝搬モデルにおける送受信点間の電波伝搬特性を幾何光学的に算出する電波伝搬特性解析方法であって、送信点から擬似光線を放射するときに放射する方位を複数の方向に分割する分割手順と、分割された方向に含まれる擬似光線を一本ごとに順に複数の計算機に個別に割り振る割振り手順と、複数の計算機が、それぞれ自己に割り振られた擬似光線の軌跡を追跡してその擬似光線の伝搬特性を計算する計算手順とを含むことを特徴とする。
【０００７】
また、本発明の電波伝搬特性解析方法は、有限差分時間領域法を用いてアンテナ複素指向性を導出する手段を備え、送信点を中心とした所定の離散的角度の方位に対するアンテナ複素指向性の振幅があらかじめ設定した閾値以下の場合は、割振り手順においては分割された方位に含まれる擬似光線を計算機に割り振らず、計算手順においては擬似光線の伝搬特性の計算を行わないことを特徴とする。
【０００８】
また、本発明の電波伝搬特性解析装置は、送信点より複数の擬似光線を放射し、その擬似光線が受信点へ到達するまでの軌跡を追跡することにより、所定の電波伝搬モデルにおける送受信点間の電波伝搬特性を幾何光学的に算出する電波伝搬特性解析装置であって、送信点から擬似光線を放射するときに放射方位を複数の方向に分割し、分割された方向に含まれる擬似光線を一本ごとに分配する主計算手段と、主計算手段から自己に割り振られた擬似光線を取得し、その擬似光線の軌跡を追跡することによって電波伝搬特性を計算する複数の端末計算手段とを備え、主計算手段が複数の端末計算手段の計算結果に基づいて、所定の電波伝搬モデルにおける送受信点間の電波伝搬特性を幾何光学的に算出するように構成されている。
【０００９】
さらに、本発明の電波伝搬特性解析装置は、有限差分時間領域法を用いてアンテナ複素指向性を導出する導出手段を備え、主計算手段が複数の端末計算手段の計算結果と導出手段によって導出されたアンテナ複素指向性とに基づいて、電波伝搬特性を幾何光学的に算出するように構成されている。
【発明の効果】
【００１０】
本発明によれば、送信点から擬似光線を放射するときに放射する方位を複数の方向に分割し、分割された方向に含まれる擬似光線を一本ずつ順に複数の計算機に割り振り、それぞれの計算機によって擬似光線の軌跡を追跡するようにしている。これによって、ある特定方向の擬似光線の計算量が多くなるような電波伝搬環境においても、各計算機の間で計算量の偏りがなくなり、かつ計算時間を短縮化することができるので、幾何光学的な電波伝搬の解析計算を短時間で行うことができる。したがって、本発明によって高速処理が可能な幾何光学的な電波伝搬特性解析装置を構築することができる。
【００１１】
また、本発明によれば、有限時間領域法を用いてアンテナ複素指向性を導出し、送信点を中心としてある離散的角度方位に対するアンテナ複素指向性の振幅が所定の閾値より小さい場合は、分割された方位に含まれる擬似光線の計算処理を計算機に割り振らないようにしたので、電波伝搬特性の解析に影響を与えないような擬似光線の計算を省略することができる。したがって、電波伝搬特性解析の処理時間をさらに高速化することができる。
【発明を実施するための最良の形態】
【００１２】
本発明にかかる電波伝搬特性解析方法及び電波伝搬特性解析装置の骨子は、無線局から送信された複数の光線を探索する幾何光学的手法を用い、送信点より複数の擬似光線を放射して、それらの擬似光線が受信点へ到達するまでの軌跡を複数の計算機によって追跡することにより、所定の電波伝搬モデルにおける送受信点間の電波の伝搬特性を算出するものである。これにより、各計算機間における計算量の偏りがなくなり、かつ電波伝搬特性解析の計算時間を短縮することができる。
【００１３】
以下、図面を参照しながら、本発明における電波伝搬特性解析方法及び電波伝搬特性解析装置の一実施の形態を詳細に説明する。なお、以下の説明では、各図面において共通する構成要素は同一符号を付し、かつ重複する説明は省略する。
【００１４】
図１は、本発明の一実施の形態にかかる幾何光学的手法による電波伝搬特性解析方法を示す概念図である。つまり、図１は、平面上において送信点から放射される擬似光線を端末計算機に取り込む状態を示している。また、図２は、図１の立体的概念図であり、送信点を中心として放射方向を分割し、その放射方向に含まれる擬似光線を端末計算機に割り振る状態を示す概念図である。つまり、図２は、送信点１０１を中心としてＮ個の方位に分割して放射グループ１からＮまでを作成し、Ｎ個に分割された各方位に含まれる擬似光線２０１をＭ個の端末計算機１０４−１、１０４−２、…、１０４−ｍに割り振った幾何光学的手法による電波伝搬特性解析方法を示している。
【００１５】
図３は、本発明の一実施の形態にかかる電波伝搬特性解析装置の主計算機と端末計算機の関係を示すブロック図である。つまり、図３は、幾何光学的手法によって実現される電波伝搬特性解析装置に使用される計算機システムの構成を示しており、主計算機３０１とＭ個の端末計算機１０４−１、１０４−２、…、１０４−ｍは、ネットワーク３０２を介してデータの送受信ができるように構成されている。図４は、図３に示す各計算機の構成を示すシステム構成図である。図４に示すように、各計算機（つまり、主計算機３０１及び各端末計算機１０４−１、１０４−２、…、１０４−ｍ）は、電波伝搬特性の演算処理を行うパーソナルコンピュータ（以下、パソコンと云う）などの制御装置４０１と、制御装置４０１の情報を取得して表示を行うＣＲＴなどの表示装置４０２と、制御装置４０１の情報を取得して印刷を行うプリンタなどの印刷装置４０３と、データの入力を行うキーボードやマウスなどの入力装置４０４と、電波伝搬特性の解析結果を記憶するハードディスクなどの記憶装置４０５と、各計算機間で通信を行う通信装置４０６とを備えた構成となっている。また、図４の計算機システムは、脱着可能なＣＤやフロッピーディスク（登録商標）などの記憶媒体４０７を使用することができる。さらに、記憶装置４０５及び記憶媒体４０７は、幾何光学的に電波の伝搬解析を実行するためのプログラムやパラメータなどを記憶することができる。
【００１６】
次に、図１から図４を用いて、本発明における幾何光学的な電波伝搬特性解析方法を概略的に説明する。なお、図１は、図２の説明を簡便にするために平面上に示した図であるので、以下の説明では図１を主体的に用いて、各端末計算機１０４−１、１０４−２、…、１０４−ｍが送信点１０１からの擬似光線２０１を取り込む概念を説明する。
【００１７】
まず、送信点１０１を中心として、離散角度Δθの間隔で各方向へ放射する擬似光線２０１を放射する方向を求める。さらに、送信点１０１を中心として、離散角度ΔθａでＮ個の方位に分割した放射グループ１、２、…、Ｎを作成する。そして、放射グループ１に含まれる擬似光線１−１を端末計算機１０４−１に割り当て、擬似光線１−２を端末計算機１０４−２に割り当て、擬似光線１−Ｋを端末計算機１０４−ｍに割り当てるというように、Ｋ本の擬似光線２０１を端末計算機１０４−１から順に端末計算機１０４−ｍまで割り当てる。なお、擬似光線を総称する場合は擬似光線２０１と符号を付し、個々の擬似光線を示す場合は、擬似光線１−１、擬似光線１−２、…、擬似光線１−Ｋなどというように個別の符号を付すことにする。
【００１８】
このようにして、放射グループ１に含まれるＫ本の擬似光線２０１をＭ個の端末計算機１０４−１、１０４−２、…、１０４−ｍに割り当てたら、放射グループ２に含まれるＫ本の擬似光線２０１をＭ個の端末計算機１０４−１、１０４−２、…、１０４−ｍに順に割り当てる。このような手法によって、放射グループＮまでに含まれるすべての擬似光線２０１をＭ台の端末計算機１０４−１、１０４−２、…、１０４−ｍに順に割り当てる。なお、送信点１０１から放射される擬似光線２０１を割り当てる際に、割り当てた端末計算機において以前に割り当てた擬似光線２０１の計算が終了していない場合は、計算が終了している次の端末計算機に対して先に擬似光線２０１を割り当てる。
【００１９】
次に、フローチャートを用いて、本発明における幾何光学的な電波伝搬特性解析方法をさらに詳細に説明する。図５は、本発明における幾何光学的な電波伝搬特性解析方法を実現するための主計算機の前処理を示すフローチャートである。つまり、図５は主計算機３０１が送信点１０１から送信される全擬似光線２０１をＭ個の端末計算機１０４−１、１０４−２、…、１０４−ｍに割り振るときの前処理の流れを示すフローチャートである。なお、以下の説明では、送信点から送信される全光線を計算する処理を前処理と定義し、受信点に到達する光線を探索して受信強度を求める処理を後処理と定義することにする。
【００２０】
図６は、本発明における幾何光学的な電波伝搬特性解析方法を実現するための各端末計算機の前処理を示すフローチャートである。つまり、図６は、Ｍ台の端末計算機１０４−１、１０４−２、…、１０４−ｍの各々が、自己に割り当てられた擬似光線について伝搬解析を行うときの前処理の流れを示すフローチャートである。図７は、図６に示す各端末計算機の前処理によって記憶装置に記憶されるデータ構成を示す構成図である。図７に示すように、各端末計算機１０４−１、１０４−２、…、１０４−ｍの前処理によって、記憶装置４０５には、全光線の方程式７０１、全光線の進行角度７０２、全光線の偏波方程式７０３、全光線の反射面番号７０４、及び全光線の壁面番号７０５が記憶される。
【００２１】
図８は、本発明における幾何光学的な電波伝搬特性解析方法を実現するための主計算機の後処理を示すフローチャートである。つまり、図８は主計算機３０１が各受信点に到達した擬似光線２０１のデータを各端末計算機１０４−１、１０４−２、…、１０４−ｍから取得して受信強度を求める後処理の流れを示すフローチャートである。図９は、本発明における幾何光学的な電波伝搬特性解析方法を実現するための端末計算機の後処理を示すフローチャートである。つまり、図９は、各端末計算機１０４−１、１０４−２、…、１０４−ｍが、各放射グループ１、２、…、Ｎの受信点に到達する擬似光線２０１を探索する後処理の流れを示すフローチャートである。
【００２２】
以下、図５から図１０を用いて、本発明における一実施の形態の幾何光学的な電波伝搬特性解析の計算手順について説明する。なお、必要に応じて図１から図４の各図面を参照しながらフローチャートの流れを説明する。
【００２３】
最初に、図５に示す主計算機が行う電波伝搬特性解析の前処理の流れを説明する。まず、パラメータ入力の手順において、オペレータは、光線の最大反射回数Ｒmax、閾値、周波数、壁面材質、放射擬似光線数Ｊ、送信位置、放射グループ数Ｎ、及びアンテナ複素指向性などの各種のパラメータを入力装置４０４から入力する（ステップＳＴ１）。この処理により、壁面方程式導出の手順において、制御装置４０１は、入力されたパラメータに対応する壁面の形状に基づいて壁面方程式を導き出し、この壁面方程式を記憶装置４０５に記憶させる（ステップＳＴ２）。
【００２４】
さらに、放射擬似光線方程式導出の手順において、送信点１０１を中心として先のパラメータ入力手順のステップＳＴ１で入力した放射擬似光線数Ｊによって求めたある一定の離散角度Δθごとに、全方位へ放射する擬似光線２０１の放射擬似光線方程式を求める（ステップＳＴ３）。ここで、擬似光線番号ｊが入力した放射擬似光線数Ｊより小さければ、その擬似光線番号ｊを１つ加算して（つまり、擬似光線番号をｊ＋１として）、擬似光線番号ｊが放射擬似光線数Ｊに到達するまで放射擬似光線方程式の導き出しを続ける（ステップＳＴ４）。
【００２５】
次に、送信点１０１を中心として放射される擬似光線２０１の放射方向をある一定の離散角度Δθａで区切り、放射方向ごとにＮ個の放射グループにグループ分けする（ステップＳＴ５）。そして、各放射グループの擬似光線２０１をＭ個の端末計算機１０４−１、１０４−２…１０４−ｍに順に割り振る。つまり、各放射グループに含まれる擬似光線２０１について、図１に示す放射グループ１に含まれる擬似光線１−１から１−ＫまでのＫ本の擬似光線２０１を、端末計算機１から端末計算機Ｍまでの処理が行われていない各端末計算機に順に割り当てる（ステップＳＴ６）。
【００２６】
このようにして、放射グループ１の擬似光線の処理が終了したら、次に、放射グループ２に含まれる擬似光線の処理を端末計算機に割り当てるというように、放射グループＮになるまで擬似光線番号ｊを１ずつ加算してゆき（つまり、擬似光線番号をｊ＋１としてゆき）、放射擬似光線数Ｊに到達するまで同様の処理（すなわち、各擬似光線をそれぞれの端末計算機へ割り当てる処理）を繰り返す（ステップＳＴ７）。このとき、各端末計算機１０４−１、１０４−２…１０４−ｍに擬似光線２０１を割り当てる際、送信点１０１から放射される擬似光線２０１の放射擬似光線方程式や、ステップＳＴ１のパラメータ入力手順で入力された各種のパラメータなどのデータを各端末計算機１０４−１、１０４−２、…、１０４−ｍへ送信する処理を併せて行う。
【００２７】
次に、図６に示す各端末計算機１０４−１、１０４−２、…、１０４−ｍの前処理の流れを説明する。図６に示す各端末計算機１０４−１、１０４−２、…、１０４−ｍの前処理では、反射光線方程式導出の手順において、ある特定の擬似光線２０１の反射壁面を決定して反射光線方程式を導き出す（ステップＳＴ１１）。そして、最大反射回数Ｒになるまで反射回数ｒを１ずつ加算してゆき（つまり、反射回数をｒ＋１と加算してゆき）、反射光線方程式の導出処理を繰り返す（ステップＳＴ１２）。
【００２８】
図６に示すような各端末計算機１０４−１、１０４−２、…、１０４−ｍの前処理により、図５のステップＳＴ１におけるパラメータ入力の手順によって入力されたデータと、図７に示すような全光線の方程式７０１、全光線の進行角度７０２、全光線の偏波方程式７０３、全光線の反射面番号７０４、及び全光線の壁面番号７０５とが、各端末計算機１０４−１、１０４−２、…、１０４−ｍの記憶装置４０５に記憶される。このようにして、各端末計算機１０４−１、１０４−２、…、１０４−ｍにおいて擬似光線の反射光線方程式の導出処理（ステップＳＴ１１）が終了するごとに、各端末計算機１０４−１、１０４−２、…、１０４−ｍは、ネットワーク３０５を介して主計算機３０１に対して処理が終了したことを示す情報を伝える。
【００２９】
次に、図８に示す主計算機における後処理の流れを説明する。主計算機３０１における後処理のパラメータ入力の手順において、オペレータは、アンテナの複素指向性、受信位置、及び周波数を入力装置４０４より入力する（ステップＳＴ２１）。そして、通信装置４０６が、パラメータ入力によって得られた各種のデータをＭ個の端末計算機１０４−１、１０４−２、…、１０４−ｍに送信する。主計算機３０１は、Ｍ個の端末計算機１０４−１、１０４−２、…、１０４−ｍから擬似光線２０１のヒットレイ情報を取得する（ステップＳＴ２２）。
【００３０】
このようにして、主計算機３０１が、Ｍ個の端末計算機１０４−１、１０４−２、…、１０４−ｍから受信位置に到達した擬似光線２０１を取得すると、位相合成の手順において、受信位置に到達した全ての擬似光線２０１を位相合成する（ステップＳＴ２３）。そして、受信位置番号ｐが受信位置数Ｐになるまで受信位置番号ｐに１を加算し（つまり、受信位置番号をｐ＋１とし）、前述と同様の処理（すなわち、位相合成の処理）を繰り返す（ステップＳＴ２４）。そして、全ての受信位置数Ｐについて位相合成が完了すると、結果出力の手順において、各々の受信地点における受信強度を記憶装置４０５に記憶させる（ステップＳＴ２５）。
【００３１】
次に、図９に示す端末計算機における後処理の流れを説明する。まず、制御装置４０１は、受信球導出の手順において、記憶装置４０５から読み込んだ光線番号ａの受信地点までの距離ｄを導き出し、送信角度間隔をφとしたとき、後処理における受信球の半径ｒを次の式（１）によって求める（ステップＳＴ３１）。
ｒ＝２１/２ｄ・sin（φ/２）                                        （１）
【００３２】
次に、制御装置４０１は、反射光線方程式と受信点を中心とする半径ｒの球との交点を求めることにより到達光線の探索を行い、電波が受信点に到達したか否かを判定する（ステップＳＴ３２）。ここで、電波が受信点に到達した場合は（ステップＳＴ３２でＹＥＳの場合）、受信点に到達した光線に光線到来角度に対応した複素指向性を乗算する（ステップＳＴ３３）。一方、光線番号ａの光線が受信点に到達していないときは（ステップＳＴ３２でＮＯの場合）、光線番号ａに１を加算してａ＋１にして（ステップＳＴ３４）、ステップＳＴ３１の受信球導出の手順に戻り、次の光線番号（ａ＋１）の受信球導出の処理を行う。
【００３３】
ステップＳＴ３３の到達擬似光線に複素指向性を乗算する手順では、受信位置番号ｐごとに受信点に到達した擬似光線に対して擬似光線到来角度における複素指向性を乗算し、到来した擬似光線の電界強度と受信位置を記憶装置４０５に記憶させる。そして、受信位置番号ｐが受信位置数Ｐになるまで受信位置番号ｐに１を加算し（つまり、受信位置番号をｐ＋１とし）、前述と同様の処理（すなわち、到来光線に複素指向性を乗算する処理）を繰り返し、それぞれ到来した擬似光線の電界強度と受信位置を記憶装置４０５に記憶させる。（ステップＳＴ３５）。
【００３４】
上述したような本発明の実施の形態における幾何光学的な電波伝搬解析方法によれば、ある特定の方向の電波の計算量が多くなるような電波伝搬環境においても、各計算機の間で計算量の偏りがなく、電波伝搬の解析時間を従来の方法の８０％程度の時間に短縮して計算することができる。
【００３５】
なお、上記の実施に形態の説明においては、図５のフローチャートにおけるステップＳＴ６の各放射グループの擬似光線をＭ個の端末計算機に順に割り振る手順では、放射方向の複素指向性の振幅の違いによらず、全ての放射グループに含まれる擬似光線をＭ個の端末計算機に順に割り振っていた。しかし、このような処理に限定されることなく、ある放射グループの放射方向の複素指向性の振幅が小さい場合は擬似光線を端末計算機に割り振らず、その放射グループに含まれる擬似光線の処理を行わないようにすることによって、電波伝搬解析の計算時間を３０％程度減少させることができる。
【００３６】
また、図５のフローチャートにおけるステップＳＴ１や図８のフローチャートにおけるステップＳＴ２１のパラメータ入力の手順において、アンテナ複素指向性にＦＤＴＤ法を使用することにより、実際の多重波環境における人体への影響を含んだアンテナの性能評価を比較的高速に行うことが可能となる。さらに、上記の実施の形態の説明においては、送信点から擬似光線を放射することによって電波伝搬解析を行っているが、これに限らず、受信点から擬似光線を放射して、送信点に到達する擬似光線を探索する場合でも前述と同様に電波伝搬解析を実施することが可能である。
【産業上の利用可能性】
【００３７】
以上説明したように、本発明にかかる幾何光学的な電波伝搬解析方法は、送信点を中心としてある一定角度ごとに放射グループをグループ分けし、その放射グループに含まれる擬似光線を複数の端末計算機に割り当てている。このようにして、各端末計算機が割り当てられた擬似光線を個別に処理することにより、各端末計算機の間で計算量の偏りが生じることなく短時間で電波伝搬解析を行うことが可能となる。従って、本発明は複数の端末計算機を用いて高速に電波伝搬解析を行うような電波伝搬解析システムに利用することができる。
【図面の簡単な説明】
【００３８】
【図１】本発明の一実施の形態にかかる幾何光学的手法による電波伝搬特性解析方法を示す概念図
【図２】図１の立体的概念図であり、送信点を中心として放射方向を分割し、その放射方向に含まれる擬似光線を端末計算機に割り振る状態を示す概念図
【図３】本発明の一実施の形態にかかる電波伝搬特性解析装置の主計算機と端末計算機の関係を示すブロック図
【図４】図３に示す各計算機の構成を示すシステム構成図
【図５】本発明における幾何光学的な電波伝搬特性解析方法を実現するための主計算機の前処理を示すフローチャート
【図６】本発明における幾何光学的な電波伝搬特性解析方法を実現するための各端末計算機の前処理を示すフローチャート
【図７】図６に示す各端末計算機の前処理によって記憶装置に記憶されるデータ構成を示す構成図
【図８】本発明における幾何光学的な電波伝搬特性解析方法を実現するための主計算機の後処理を示すフローチャート
【図９】本発明における幾何光学的な電波伝搬特性解析方法を実現するための端末計算機の後処理を示すフローチャート
【符号の説明】
【００３９】
１〜Ｎ  放射グループ
１０１  送信点
１０４−１，１０４−２…１０４−ｍ  端末計算機
２０１  擬似光線
３０１  主計算機
３０２  ネットワーク
４０１  制御装置
４０２  表示装置
４０３  印刷装置
４０４  入力装置
４０５  記憶装置
４０６  通信装置
４０７  記憶媒体
７０１  全光線の方程式
７０２  全光線の進行角度
７０３  全光線の偏波方程式
７０４  全光線の反射面番号
７０５  全光線の壁面番号
【図１】
2005331391.tif 000003
【図２】
2005331391.tif 000004
【図３】
2005331391.tif 000005
【図４】
2005331391.tif 000006
【図５】
2005331391.tif 000007
【図６】
2005331391.tif 000008
【図７】
2005331391.tif 000009
【図８】
2005331391.tif 000010
【図９】
2005331391.tif 000011