JP4273665B2 - レイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラム - Google Patents
レイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラム Download PDFInfo
- Publication number
- JP4273665B2 JP4273665B2 JP2001008368A JP2001008368A JP4273665B2 JP 4273665 B2 JP4273665 B2 JP 4273665B2 JP 2001008368 A JP2001008368 A JP 2001008368A JP 2001008368 A JP2001008368 A JP 2001008368A JP 4273665 B2 JP4273665 B2 JP 4273665B2
- Authority
- JP
- Japan
- Prior art keywords
- ray
- regular hexagonal
- cone
- similar
- hexagonal pyramid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 75
- 230000008569 process Effects 0.000 claims description 22
- 230000005684 electric field Effects 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 description 36
- 230000004888 barrier function Effects 0.000 description 22
- 238000010586 diagram Methods 0.000 description 20
- 238000012545 processing Methods 0.000 description 17
- 238000004364 calculation method Methods 0.000 description 7
- 238000007796 conventional method Methods 0.000 description 7
- 230000006866 deterioration Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 238000000491 multivariate analysis Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003936 working memory Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
- H04W16/18—Network planning tools
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Radar Systems Or Details Thereof (AREA)
- Image Generation (AREA)
- Mobile Radio Communication Systems (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明はレイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラムに関し、特に幾何光学的手法による電波伝搬特性予測におけるレイの空間分解能制御方式に関するものである。
【０００２】
【従来の技術】
無線通信システムにおける基地局や親機等の配置を支援するために電波伝搬シミュレータが用いられる。電波伝搬シミュレータによって任意の受信点における受信電力や遅延拡がりを評価できる。この情報をもとに効果的な送信局の設置場所を決定し、配置が必要な基地局数を削減する等の効率化が達成される。
【０００３】
電波伝搬シミュレーションは統計的手法と決定論的手法とに大別できる。統計的手法は、距離や周波数などを引数とする伝搬損失推定式を用いる手法であり、推定式中のパラメータは、伝搬損失の実測定で得られた多数のデータをもとに多変量解析等の統計的手法により決定される。一方、決定論的手法は、アンテナから放射される電波を多数の電波線（レイ）の集まりとして表現する手法であり、各レイは幾何光学的に反射・透過を繰り返しながら伝搬するものとする。観測点における伝搬損および遅延量は、観測点に到達する各レイの電力と伝搬時間を合成することにより求めることができる。
【０００４】
決定論的手法は、伝搬経路を追跡する方法の違いから、更にイメージング法とラウンチング法とに大別できる。イメージング法は送受信点間を結ぶレイの反射・透過経路を反射面に対する鏡映点を求めて決定する手法である。反射・透過経路は、送受信点と反射・透過障壁の位置が決まれば一意に求まるために、イメージング法は厳密なレイの伝搬経路を探索する手法である。
【０００５】
一方、ラウンチング法は、受信点の位置に関わり無く、アンテナから離散的な角度間隔でレイを放射し、反射・透過等を繰返しながら受信点近傍を通過するレイを当該受信点に到達したレイとみなす手法である。このラウンチング法は、送受信点間を結ぶレイの伝搬経路の解をイメージング法のように厳密に求めるのではなく、近似的に求めるために、伝搬経路探索に要する時間を短縮できるという特徴がある。
【０００６】
ラウンチング法では、まず、送信アンテナの周囲に閉領域を定め、その閉領域を部分空間に分割する。次に、各部分空間に対して１本のレイを割当て、部分空間と割当てたレイを同一視する。その後、送信アンテナから放射される有限な本数のレイの伝搬経路を追跡し、その結果を基にアンテナの周囲全空間の電波伝搬状況を推定する。
【０００７】
アンテナ周囲の閉領域を部分空間に分割する方法の一例が、文献「Scott Y. Seidel,et al ：“Site-Specific Propagation Prediction for Wireless In-Building Personal Communication System Design”, IEEE TRANSACTIONS VEHICULAR TECHNOLOGY，VOL.43，NO.4 ,NOVEMBER 1994 ，pp.879-891」に説明されている。
【０００８】
この文献では、先ず、図１２のように、送信アンテナ３０１の周囲に、正２０面体の３次元閉領域を設ける。次に、正２０面体を構成する面、すなわち図１３のように、頂点４０６，４０７，４０８からなる正三角形の平面を取出した後、点４０６，４１０，４１１を用いて各辺を等間隔に分割する。頂点４０６，４０７，４０８からなる正三角形の各辺に平行で、かつ分割点を通る線分を引くことにより、もとの正三角形と相似な三角形が新しく内部に作られる。以上の処理を図１２の正２０面体をなす全ての面に対して行い、正２０面体の重心と新しく生成された各正三角形の頂点とを結ぶ方向へ、重心からの距離が等しくなるように各正三角形の頂点を移動させると、例えば図１４のようになる。
【０００９】
図１４は図１２の正２０面体の各面をなす正三角形の１辺を２等分した場合の図である。もとの正２０面体の重心に位置する送信アンテナ５０１から放射されるレイは、送信アンテナ５０１と図１４の多面体の各頂点とを結ぶ各方向に放射される。図１４には、一例として頂点５０２を通るレイ５０４を示してある。このとき、図１３の多角形と送信アンテナ５０１から定まる部分空間５０３はレイ５０４と同一視される。
【００１０】
図１５はレイ６０１と同一視される部分空間６０５を示す図である。点６０２において、レイ６０１に垂直な断面６０３の面積を、今後、点６０２における空間分解能と定義する。点６０２における空間分解能は送信アンテナ６０４と点６０２との間の長さ、すなわち当該レイの伝搬距離が長くなるに従い増大する。このように広い空間を１本のレイと同一視するようになるため、ラウンチング法による電波伝搬推定精度は、伝搬距離が増大するに従い低下することになる。
【００１１】
そこで、空間分解能を伝搬距離に関わらず、常に一定の値以下に維持することを目的として、所定の距離を伝搬するごとにレイを分割する手法が提案された。この手法の一例が、文献「Steven Fortune：“EFFICIENT ALGORITHMS FOR PREDICTION OF INDOOR RADIO PROPAGATION ”,in Proceedings of the 48th IEEE Vehicular Technology Conference,May,1998,pp.572-576」に記されている。図１６はこの文献の手法を説明する図である。
【００１２】
図１６の三角錐７１２は、割当てられたレイ７０２の送信アンテナ７０１からの伝搬距離が所定の値に達すると、三角錐７０８から三角錐７１１に細分される。細分された三角錐には、それぞれレイ７０４からレイ７０７が割当てられ、以降、同様な処理が繰返される。この結果、空間分解能は伝搬距離に関わらず常に一定値以下に保たれることになるのである。
【００１３】
図１７は観測領域０１７と、送信点０１５と、受信点０１６、２つの物体００１，００２とが与えられた場合のラウンチング法の動作について説明する図である。図１７では、簡単のため２次元平面に限定して動作の説明をするが、実際には３次元空間内で動作させることは勿論である。
【００１４】
離散的な間隔で伝搬経路００３の方向へ放射されたレイは、点０１８で物体００１へ入射し、その結果、反射レイ００５と透過レイ００４とが生成される。反射レイ００５は、さらに点０１９で物体００２へ入射し、再び反射レイ００６と透過レイ００７とが生成される。生成された反射レイ００６は受信点０１６の近傍を通過するため、レイ００６を受信点に到来したレイとみなし、経路００３，００５，００６からなる経路を送信点０１５と受信点０１６とを結ぶ伝搬経路の一つとみなす。
【００１５】
伝搬経路００３，００５，００６から求まる受信電力ならびに到来遅延時間は図１８のように記録される。図１８の横軸１０３はレイが送信点０１５から受信点０１６まで経路００３，００５，００６を経由して到達するのに要した遅延時間を示し、縦軸１０２は当該経路を経て到達したレイが有する電力強度を表す。透過レイ００４ならびに００７についても、伝搬経路００３，００５，００６と同様の反射・透過レイ探索を繰り返し、受信点０１６の近傍をレイが通過した場合には、伝搬経路００６の場合と同様に到来レイとして扱い、以上の処理を探索終了条件を満たすまで続ける。
【００１６】
探索終了条件は、反射・透過点での受信電界強度が所定値を下回った場合、または反射・透過の総数が所定の回数に達した場合などとする。送信点０１５から伝搬経路００３の方向へ放射されたレイの反射・透過経路探索が終了した後、方向００８から０１４のように他の放射方向へ向けてレイを放射した場合についても同様に反射・透過経路探索を行う。事前に定めた全ての方向へレイを放射し、以上の伝搬経路探索を行うことにより、受信点０１６における遅延プロファイルである図１９が得られる。図１９の横軸２０３は送信点０１５から受信点０１６へ到達するのに要した時間を示し、縦軸２０２は当該経路を通過したレイが有する電力強度を表す。
【００１７】
受信電力は図１９に示した全ての経路の電力強度の和で与えられ、歪の程度を示す遅延拡がりは、各遅延時間において受信電力に対する当該電力強度の電力比を当該遅延時間の出現確率として、遅延時間の標準偏差により与えられる。
【００１８】
【発明が解決しようとする課題】
上述したラウンチング法では、レイと近傍の部分空間が同一視されるために、レイが障壁と反射する場合は、レイ近傍の部分空間でも同様な反射が生じているとみなされる。そのため、推定精度の観点からは、反射点における空間分解能を当該障壁の面積に比べて大きくなり過ぎないように設定することが望ましい。しかし、推定精度の向上を目的として全体の空間分解能の上限値を各障壁との反射の際に要求される空間分解能の最小値に設定すると、電波伝搬推定に要する計算時間が増大するという問題点がある。その理由は、空間分解能を小さくすると放射するレイの総数が増大し、経路探索に要する全体の時間が増大するからである。
【００１９】
一方、推定に要する計算時間を抑えるために空間分解能の上限値を比較的大きく設定すると、全体の推定精度に劣化が生じるという問題点がある。その理由は、所要空間分解能が小さな障壁でレイが反射する場合には、実際には反射が生じていないレイ近傍の領域でもレイと同様に反射が生じているとみなすからである。
【００２０】
本発明は、ラウンチング法において伝搬推定の対象環境が与えられた場合に、推定精度の劣化を抑えながら所要推定時間を短縮化する手法を提供することを目的とする。
【００２１】
【課題を解決するための手段】
本発明によれば、３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御方法であって、前記遮蔽物と前記レイとの衝突点を検出する第一のステップと、前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う第二のステップと、この第二のステップで前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する第三のステップと、この第三のステップにおいて分割された前記部分空間の各々に対して同一視するレイを割当てる第四のステップとを含むことを特徴とするレイの空間分解能制御方法が得られる。
【００２３】
そして、前記第三のステップにおける分割数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは１〜Ｍ）の前記部分空間に対して前記第四のステップにおいてレイを割当てた後、割当てた各レイに対して再び前記第一から第四のステップを適用し、前記第二のステップにおいて当該レイと同一視される前記空間を分割しないと判定されるまで前記第一から第四のステップを、順次繰返すことを特徴とする。
【００２４】
また、前記第二のステップにおける前記判定基準は、前記遮蔽物との衝突点における前記レイの電界強度が所定の値以上であるか否かを調べる第一の判定基準と、前記遮蔽物の衝突面を含む平面による前記レイ近傍の前記空間の断面に対する前記遮蔽物の衝突面の面積比が所定値以下であるか否かを調べる第二の判定基準とからなることを特徴とし、更には、前記判定基準は、前記断面の面積が所定の値以上であり、かつ前記レイと前記衝突面の法線との成す角（入射角）が所定値以上であるか否かを調べる第三の判定基準を含むことを特徴とする。
【００２５】
また、錐体の中心線を前記レイの経路とすることにより前記レイ近傍の前記空間を前記錐体によって形成し、前記第三のステップにおいて、前記レイ近傍の空間を形成する前記錐体に相似な複数の錐体を用いて前記レイ近傍の空間を複数の部分空間に分割する際に、前記レイ近傍の空間を形成する前記錐体の体積と、この錐体を分割する前記錐体に相似な前記複数の錐体の和の体積との体積誤差が最小となるように、前記レイの近傍に前記複数の錐体を配置することを特徴とする。
【００２６】
更に、前記第三のステップにおいて、円錐に内接する前記錐体として三角錐を用い、前記三角錐に相似な複数の同一三角錐を前記三角錐内部に隙間無く配置できるように前記相似三角錐の大きさを設定することを特徴とし、また前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和の体積と、もとの前記正六角錐の体積誤差が最小となるように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００２７】
更にはまた、前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面がもとの前記正六角錐の側面に内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。また、前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面に、もとの前記正六角錐が内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００２９】
本発明によれば、３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御システムであって、前記遮蔽物と前記レイとの衝突点を検出する衝突点検出手段と、前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う判定手段と、この判定手段で前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する分割手段と、この分割手段において分割された前記部分空間の各々に対して同一視するレイを割当てる割当手段とを含むことを特徴とするレイの空間分解能制御システムが得られる。
【００３０】
そして、前記分割手段における分割数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは１〜Ｍ）の前記部分空間に対して前記割当手段においてレイを割当てた後、割当てた各レイに対して再び前記衝突点検出手段、前記判定手段、前記分割手段及び前記割当手段を実行制御し、前記判定手段において当該レイと同一視される前記空間を分割しないと判定されるまで、前記衝突点検出手段、前記判定手段、前記分割手段及び前記割当手段を、順次繰返し実行制御する制御手段とを含むことを特徴とする。
【００３１】
また、前記判定手段における前記判定基準は、前記遮蔽物との衝突点における前記レイの電界強度が所定の値以上であるか否かを調べる第一の判定基準と、前記遮蔽物の衝突面を含む平面による前記レイ近傍の前記空間の断面に対する前記遮蔽物の衝突面の面積比が所定値以下であるか否かを調べる第二の判定基準とからなることを特徴とし、更に、前記判定基準は、前記断面の面積が所定の値以上であり、かつ前記レイと前記衝突面の法線との成す角（入射角）が所定値以上であるか否かを調べる第三の判定基準を含むことを特徴とする。
【００３２】
また、錐体の中心線を前記レイの経路とすることにより前記レイ近傍の前記空間を前記錐体によって形成し、前記分割手段において、前記レイ近傍の空間を形成する前記錐体に相似な複数の錐体を用いて前記レイ近傍の空間を複数の部分空間に分割する際に、前記レイ近傍の空間を形成する前記錐体の体積と、この錐体を分割する前記錐体に相似な前記複数の錐体の和の体積との体積誤差が最小となるように、前記レイの近傍に前記複数の錐体を配置することを特徴とする。
【００３３】
更に、前記分割手段において、円錐に内接する前記錐体として三角錐を用い、前記三角錐に相似な複数の同一三角錐を前記三角錐内部に隙間無く配置できるように前記相似三角錐の大きさを設定することを特徴とし、また、前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和の体積と、もとの前記正六角錐の体積誤差が最小となるように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００３４】
また、前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面がもとの前記正六角錐の側面に内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。更に、前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面に、もとの前記正六角錐が内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００３６】
本発明によれば、３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御方法をコンピュータに実行させるためのプログラムであって、前記遮蔽物と前記レイとの衝突点を検出する第一の処理と、前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う第二の処理と、この第二の処理で前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する第三の処理と、この第三の処理において分割された前記部分空間の各々に対して同一視するレイを割当てる第四の処理とを含むことを特徴とするプログラムが得られる。
【００３７】
本発明の作用を述べる。本発明では、与えられた観測領域内をレイが伝搬する場合、その経路上に存在する障壁とレイが衝突する際に、所定の条件が満たされるまで当該レイ近傍の部分空間を複数の部分空間に分割する。従来の方法では、反射とは無関係に伝搬距離に応じて一律に部分空間を分割していたが、本発明の方法によれば、一律に扱った場合に推定精度の劣化を招くと考えられる部分空間は、劣化の影響が少なくなるまで分割し、細分化された各部分空間を個別に取り扱うことが可能となる。このような手法を採用することにより、電波伝搬の推定精度を従来方法よりも向上させることが可能となる。一方、推定精度の劣化の影響が小さいと考えられる場合は、当該部分空間を分割しないため、空間分解能の初期値を比較的大きめに設定しておけば、従来の方法において推定精度の向上を狙って一律に空間分解能を小さく設定した場合よりもレイの総数が低減し、伝搬推定に要する全体の計算時間を短縮させることができる。
【００３８】
これらの特徴により、本発明では、与えられた推定領域に応じて空間分解能を適応制御することが可能なことから、従来の方法に比べて推定精度の向上と計算時間の短縮を両立させることが可能である。
【００３９】
【発明の実施の形態】
以下に、図面を参照しながら本発明の実施例を詳細に説明する。図１は本発明の実施例の動作を示す流れ図であり、２つのルーチンＡ，Ｂ（８１８，８１９）よりなるものとする。図１では、３次元空間上に観測領域が与えられ、当該観測領域内に壁や什器等の物体が与えられ、当該観測領域内の点Ｐから当該伝搬方向へ放射されるレイが与えられることを前提とし（ステップ８０２）、反射点における反射・透過レイを生成するまでの一連の流れを示す。
【００４０】
ここで、点Ｐとは、ラウンチング法における反射・透過点もしくは送信アンテナの位置を指し、具体的には従来技術の説明において参照した図１７における点０１８、点０１９、送信アンテナ０１５等を指すものとする。また、当該伝搬方向とは、点Ｐが送信アンテナである場合は、送信アンテナの周囲に設けた離散的な放射方向を指し、具体的には従来の技術で説明したように図１７の方向００３、００８〜０１４等、もしくは図１４の送信アンテナ５０１と多面体の各頂点を結ぶ方向、例えばレイ５０４の方向等を指すものとする。一方、点Ｐが反射・透過点である場合は、反射レイもしくは透過レイの方向を指し、具体的には図１７における伝搬経路００３〜００７等を指すものとする。
【００４１】
まず、点Ｐから当該伝搬方向へ放射したレイが観測領域内の物体と衝突する点、すなわち反射点を探索し検出する（ステップ８０３，８０４）。次に、反射点に入射する当該レイと同一視されている部分空間を、より小さな部分空間に細分化するか否かの判定を事前に定めた判定基準（後述）に基づいて行う（ステップ８０５）。
【００４２】
当該部分空間をより小さな部分空間に細分する場合は、事前に定めた方法（後述）に基づいて当該部分空間をＭ個（Ｍは整数）の部分空間に分割し（ステップ８０６）、ルーチンＢ（８１９）へ移行する（ステップ８０７）。当該部分空間の細分化を行わない場合には、当該反射点における反射レイと透過レイを生成し終了する（ステップ８０８，８０９）。また、ステップ８０４で反射点が検出されなければ終了となる。
【００４３】
ルーチンＢ（８１９）では、まず内部変数ｉ（ｉ＝１〜Ｍ）を定義し、初期値としてｉ＝１を与える（ステップ８１１）。次に、ルーチンＡ（８１８）のステップ８０６でＭ個に分割された部分空間のうちｉ番目の部分空間を選択する（ステップ８１２）。そして、選択したｉ番目の部分空間にレイを割当て（ステップ８２０）、そのレイの伝搬方向を事前に定めた方法（後述）に基づいて算出する（ステップ８１３）。
【００４４】
ここで、ステップ８０６で細分化された部分空間にステップ８２０で割当てられたレイが、ステップ８１３で求めた当該伝搬方向へ起点Ｐから放射されることを前提とし、再度ルーチンＡ（８１８）へ移行する。ルーチンＡ（８１８）では、細分化された部分空間をさらに細分化するかどうかの判定を行い、細分化する場合は、ステップ８０５で再度細分化処理を行い、ルーチンＢ（８１９）へ移行する。細分化を行わない場合には、反射・透過レイを生成して（ステップ８０８）終了する（ステップ８０９）。
【００４５】
ステップ８１４でルーチン８１８へ移行し、ステップ８０９でルーチン８１８を終了した場合は、その直後にステップ８１５へ移行する。ステップ８１５では、内部変数ｉの値がＭより小さいか否かの判定を行い、Ｍよりも小さい場合は、ｉに１を加算して（ステップ８１７）他の部分空間を選択する（ステップ８１２）。ｉの値がＭ以上である場合は、ルーチン８１９を終了する（ステップ８１６）。
【００４６】
ステップ８０７でルーチン８１９へ移行し、ステップ８１６でルーチン８１９を終了した場合は、その直後にステップ８０９へ移行し、入れ子となっているルーチン８１８の一つを終了する。以上の処理を、ステップ８０６で部分空間を細分化する必要がなくなるまで実行する。
【００４７】
図１に示した発明の実施例の具体的動作例を、図２〜図５を用いて説明する。図２は部分空間と障壁の交差面の生成過程を説明する図である。点９０１から放射されたレイ９０２は、点９０３において障壁９０６と衝突する。線９０８は障壁９０６に対する法線であり、角度９０７は障壁９０６に対するレイ９０２の入射角を指す。このとき、レイ９０２近傍の部分空間９０４と障壁９０６とによって、交差面９０５が生じる。図２はアンテナ周囲の空間を三角錐によって分割した場合を示す。
【００４８】
図３は本発明における部分空間の細分化を説明する図であり、図２においてθ＝０の場合を示している。図３１（ａ）は、レイ９０２が点１２１３へ入射した結果、障壁１２１２と部分空間の断面１２１１とによって生じる交差面１２１２（斜線部分）を示している。点９０３において、部分空間９０４を細分するか否かの判定基準（図１のフローのステップ８０５の判定基準）の一例として、以下の項目が挙げられる。
【００４９】
（１）点１２１３における当該レイの電界強度が所定値以上である：
（２）部分空間断面１２１１に対する交差面１２１２の面積比が所定値以下である：
（３）部分空間断面１２１１の面積が所定の値以上であり、かつレイ９０２の入射角９０７が所定値以上である：
基本的には、上記の３項目の基準のうち（１）及び（２）を同時に満たす場にを部分空間９０４の分割（細分）を行うものとするが、これ等２つの基準に加えて、（３）の基準を更に付加して、３つの基準を同時に満たす場合に部分空間の分割（細分）を行うものとしても良い。
【００５０】
細分した部分空間に割当てられるレイの伝搬方向の決定方法（図１のフローのステップ８１３の方向算出法）としては、例えば以下の方法が挙げられる。まず、点９０３と点９０１とを結ぶ線分上の任意の点においてレイ９０２に対し垂直な面で部分空間を切断する。その結果生じる、図４のような部分空間の断面１３０９をＭ個（この例では、Ｍ＝４）の相似な領域１３０５〜１３０８に分割し、点１３０１〜１３０４のように各領域の重心点を求める。同様な重心点を部分空間の別の切断面についても求める。対向する重心点同士を結ぶと、図５のように所望の伝搬方向ベクトル１９０１が得られる。
【００５１】
なお、図３〜図５はＭ＝４の場合を示したものであるが、Ｍが他の値の場合も同様である。
【００５２】
このように部分空間９０４を細分すると、反射点９０３における部分空間の断面は図３（ｂ）の様になる。図９のフロー及び上記の分割判定基準に基づくと、例えば、部分空間１２０２はこれ以上分割されずに、部分空間１２０２に割当てられたレイは障壁１２１２で反射することなく、点１２０６を通過して他の障壁に衝突するまで進行を続ける。この部分空間１２０２が分割されないのは、図１のフローにおけるステップ８０４で反射点が検出されないからであり、詳述すると、この部分空間１２０２に割当てられたレイは、点１２０６を通過してしまい、障壁１２１２で反射しないからである。
【００５３】
一方、部分空間１２０３は上記基準を満たしているとして、図３（ｃ）のように取出された後、さらに細分され、以後同様に必要なだけ細分化が行われる。
【００５４】
次に、本実施例の効果について説明する。本実施例では、当該部分空間と障壁との交差面と反射点における当該部分空間の断面の面積比を分割判定の基準とし、実際には反射が生じていない領域を分割して別扱いにすることができる（上記判定基準の（２））。また、図６に示すとおり、電波の反射率は一般に入射角が大きくなるにつれて急増することが多いが、当該レイの入射角度を判定基準とすることにより、本来は異なる角度で障壁に入射し反射率が大きく異なる領域を分割して別扱いにすることができる（上記判定基準の（３））。
【００５５】
これらの作用によって、ラウンチング法よる推定精度を向上させることができることになるのである。一方、本実施の形態では細分した後の部分空間はもとの部分空間と相似であるため、以降の細分化処理においても同様な処理を繰り返し適用することができ、処理を簡易化することが可能である。さらに、分割前の部分領域は分割後の部分空間の和で与えられるために、アンテナ周囲の全空間を隙間無く推定することが可能である。
【００５６】
以上述べた実施例の処理動作から、この処理動作を実現するための機能ブロックが、図７に示す様に得られることになる。図７を参照すると、衝突検出部１と、部分空間分割判定部２と、部分空間分割部３と、部分空間へのレイ割当部４と、制御部（ＣＰＵ）５と、ＲＡＭ（ランダムアクセスメモリ）６と、ＲＯＭ（リードオンリメモリ）７とを含んで構成される。
【００５７】
衝突検出部１はレイと遮蔽物（障壁）との衝突を検出する機能を有しており、例えば、図３の例では、障壁１２１２に対してレイが衝突する関係にあり、これは障壁１２１２の面を表す３次元空間上に規定される式と、同じくレイを表す３次元空間上に規定される式とを連立させて解く手法を採用することで、衝突の有無が判定可能である。
【００５８】
部分空間判定部２はレイと同一視される空間を分割するか否かの判定を、上述した判定基準に従って判定する機能を有しており、部分空間分割部３は上述した予め定められた規定に従って当該空間を複数の部分空間に分割する機能を有する。部分空間へのレイ割当部４は分割された部分空間の各々に対して同一視するレイを割当てる機能を有する。
【００５９】
制御部５はＣＰＵ（コンピュータの中央演算処理部に相当）からなり、ＲＯＭ７の記録媒体に予め格納されている動作制御プログラムを読込んで、このプログラムの制御手順に従って、作業用メモリであるＲＡＭ６を使用しつつ各部１〜４を制御し、図１に示した動作処理を実行制御するものである。すなわち、そのプログラムがＣＰＵの動作を制御するものであり、プログラム制御されたそのＣＰＵが、当該プログラムにより指令されて、上記の各処理を行うようになっている。
【００６０】
次に、本発明の第２の実施例について、図面を参照して説明する。第１の実施例においては、レイに垂直な平面で部分空間を切断したときに、図４に示す形態となるようにアンテナ周囲の空間を錐体で分割していた。これに対して、第２の実施例では、部分空間の断面が図８となるようにアンテナ周囲の空間を分割している点で異なる。図８では、正６角形１４０８と相似な正６角形１４０９〜１４１５で部分空間を細分化するＭ＝７の場合を示している。正６角形１４０９〜１４１４はそれぞれ辺の中点１４１６〜１４２１で正６角形１４０８の頂点と接しており、分割後に各部分空間に割当てられるレイは正６角形１４０９〜１４１５の重心１４０１〜１４０７を通る。
【００６１】
この第２の実施例で用いた正６角形のように比較的円形に近い図形を用いた場合、受信アンテナ近傍をレイが通過した場合にそのレイを取込むか否かの判定を簡易化でき処理速度を向上させることができる。その理由は、一般にこの判定が、当該レイ近傍の部分空間の内部に受信アンテナが含まれるか否かを調べることによって行われることが多いためである。そこで、当該部分空間を図９のように正６角形１６０２の外接円１６０１を当該レイに垂直な切断面とする部分空間で近似することによって、この判定は受信アンテナと当該レイの間の距離を調べることに帰着でき、判定処理の簡易化が達成できる。
【００６２】
次に、本発明の第３の実施例について図面を参照して説明する。当該部分空間の断面が図８となるようにアンテナ周囲の空間を分割した場合、分割前の断面１４０８と分割後の部分空間の和領域を比較すると、図８で斜線で示した部分のような過不足領域が生じる。そこで、分割後の部分空間が分割前の部分空間の領域外にはみ出さないような分割形態を図１０に示す。
【００６３】
一方、分割後の部分空間と分割前の部分空間との間に隙間が生じないように設定した分割形態を図１１に示す。図１０では、頂点１８０１〜１８１２が正６角形１８１３に内接するように、分割後の正６角形の大きさを設定している。一方、図１１では、分割前の正６角形の各辺の中点１５０１〜１５０６において、分割後の正６角形と交差するように、分割後の正６角形の大きさを設定している。
【００６４】
特に、図１１の分割形態を用いた場合は、第１の実施例と同様に、アンテナ周囲の空間の電波伝搬状況を隙間なく推定することができる。
【００６５】
部分空間分割の例として上述の様に種々示しているが、要は、円錐に内接する錐体のうち、この円錐との体積誤差が小さい錐体を設け、その中心線をレイの経路とすることによりレイ近傍の空間を当該錐体によって形成し、レイ近傍の空間を形成するこの錐体に相似な複数の錐体を用いてレイ近傍の空間を複数の部分空間に分割する際に、レイ近傍の空間を形成する錐体と、この錐体を分割する錐体に相似な複数の錐体の和との体積誤差が小さくなるように、レイの近傍に複数の錐体を配置すれば良い。
【００６６】
【発明の効果】
以上述べた様に、本発明によれば、従来の方法に比べてラウンチング法におけるレイの密度を柔軟に制御し、推定精度の向上と計算時間の短縮を効率よく両立できるという効果がある。その理由は、従来の方法では伝搬距離に応じて一律にレイの密度を制御していたのに対し、本発明では、反射・透過の際にレイの密度が過小で推定精度の劣化を招くと考えられる場合に、劣化が小さいと判断されるまでレイの密度を増大させるためである。これにより、推定精度を損なわない範囲で不必要なレイの生成による計算時間の増大を抑えることができる。
【図面の簡単な説明】
【図１】本発明の実施例の動作を示すフロー図である。
【図２】部分空間と障壁の交差面の生成過程を説明する図である。
【図３】本発明における部分空間の細分化を説明する図である。
【図４】本発明の第１の実施例における部分空間の分割例を示す断面図である。
【図５】分割された部分空間に割当てられたレイの伝搬方向を求める過程を説明する図である。
【図６】反射率の入射角特性の例を示す図である。
【図７】本発明の機能ブロックの例を示す図である。
【図８】本発明の第２の実施例における部分空間の分割例を示す断面図である。
【図９】部分空間の円錐近似を説明する図である。
【図１０】本発明の第３の実施例における部分空間の分割の一例を示す断面図である。
【図１１】本発明の第３の実施例における部分空間の分割の他の例を示す断面図である。
【図１２】送信アンテナ周囲の閉領域の例を示す図である。
【図１３】送信アンテナ周囲の閉領域の分割過程を説明する図である。
【図１４】離散的間隔によるレイの放射を説明する図である。
【図１５】レイと部分空間の例を示す図である。
【図１６】従来の手法による部分空間の分割を説明する図である。
【図１７】ラウンチング法の説明図である。
【図１８】一本のレイに対する観測点におけるパスプロファイルの例を示す図である。
【図１９】ラウンチング法によって得られる観測点におけるパスプロファイルの例を示す図である。
【符号の説明】
１ 衝突検出部
２ 部分空間分割判定部
３ 部分空間分割部
４ 部分空間へのレイ割当部
５ 制御部
６ ＲＡＭ
７ ＲＯＭ
００１，００２ 物体（遮蔽物）
００３〜０１４ レイ伝搬方向
０１５ 送信アンテナ位置
０１６ 受信点
０１８，０１９，９０１，９０３ 反射・透過点
９０２，１９０１ レイ
９０４ 部分空間
９０５ 交差面
９０６ 障壁
【発明の属する技術分野】
本発明はレイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラムに関し、特に幾何光学的手法による電波伝搬特性予測におけるレイの空間分解能制御方式に関するものである。
【０００２】
【従来の技術】
無線通信システムにおける基地局や親機等の配置を支援するために電波伝搬シミュレータが用いられる。電波伝搬シミュレータによって任意の受信点における受信電力や遅延拡がりを評価できる。この情報をもとに効果的な送信局の設置場所を決定し、配置が必要な基地局数を削減する等の効率化が達成される。
【０００３】
電波伝搬シミュレーションは統計的手法と決定論的手法とに大別できる。統計的手法は、距離や周波数などを引数とする伝搬損失推定式を用いる手法であり、推定式中のパラメータは、伝搬損失の実測定で得られた多数のデータをもとに多変量解析等の統計的手法により決定される。一方、決定論的手法は、アンテナから放射される電波を多数の電波線（レイ）の集まりとして表現する手法であり、各レイは幾何光学的に反射・透過を繰り返しながら伝搬するものとする。観測点における伝搬損および遅延量は、観測点に到達する各レイの電力と伝搬時間を合成することにより求めることができる。
【０００４】
決定論的手法は、伝搬経路を追跡する方法の違いから、更にイメージング法とラウンチング法とに大別できる。イメージング法は送受信点間を結ぶレイの反射・透過経路を反射面に対する鏡映点を求めて決定する手法である。反射・透過経路は、送受信点と反射・透過障壁の位置が決まれば一意に求まるために、イメージング法は厳密なレイの伝搬経路を探索する手法である。
【０００５】
一方、ラウンチング法は、受信点の位置に関わり無く、アンテナから離散的な角度間隔でレイを放射し、反射・透過等を繰返しながら受信点近傍を通過するレイを当該受信点に到達したレイとみなす手法である。このラウンチング法は、送受信点間を結ぶレイの伝搬経路の解をイメージング法のように厳密に求めるのではなく、近似的に求めるために、伝搬経路探索に要する時間を短縮できるという特徴がある。
【０００６】
ラウンチング法では、まず、送信アンテナの周囲に閉領域を定め、その閉領域を部分空間に分割する。次に、各部分空間に対して１本のレイを割当て、部分空間と割当てたレイを同一視する。その後、送信アンテナから放射される有限な本数のレイの伝搬経路を追跡し、その結果を基にアンテナの周囲全空間の電波伝搬状況を推定する。
【０００７】
アンテナ周囲の閉領域を部分空間に分割する方法の一例が、文献「Scott Y. Seidel,et al ：“Site-Specific Propagation Prediction for Wireless In-Building Personal Communication System Design”, IEEE TRANSACTIONS VEHICULAR TECHNOLOGY，VOL.43，NO.4 ,NOVEMBER 1994 ，pp.879-891」に説明されている。
【０００８】
この文献では、先ず、図１２のように、送信アンテナ３０１の周囲に、正２０面体の３次元閉領域を設ける。次に、正２０面体を構成する面、すなわち図１３のように、頂点４０６，４０７，４０８からなる正三角形の平面を取出した後、点４０６，４１０，４１１を用いて各辺を等間隔に分割する。頂点４０６，４０７，４０８からなる正三角形の各辺に平行で、かつ分割点を通る線分を引くことにより、もとの正三角形と相似な三角形が新しく内部に作られる。以上の処理を図１２の正２０面体をなす全ての面に対して行い、正２０面体の重心と新しく生成された各正三角形の頂点とを結ぶ方向へ、重心からの距離が等しくなるように各正三角形の頂点を移動させると、例えば図１４のようになる。
【０００９】
図１４は図１２の正２０面体の各面をなす正三角形の１辺を２等分した場合の図である。もとの正２０面体の重心に位置する送信アンテナ５０１から放射されるレイは、送信アンテナ５０１と図１４の多面体の各頂点とを結ぶ各方向に放射される。図１４には、一例として頂点５０２を通るレイ５０４を示してある。このとき、図１３の多角形と送信アンテナ５０１から定まる部分空間５０３はレイ５０４と同一視される。
【００１０】
図１５はレイ６０１と同一視される部分空間６０５を示す図である。点６０２において、レイ６０１に垂直な断面６０３の面積を、今後、点６０２における空間分解能と定義する。点６０２における空間分解能は送信アンテナ６０４と点６０２との間の長さ、すなわち当該レイの伝搬距離が長くなるに従い増大する。このように広い空間を１本のレイと同一視するようになるため、ラウンチング法による電波伝搬推定精度は、伝搬距離が増大するに従い低下することになる。
【００１１】
そこで、空間分解能を伝搬距離に関わらず、常に一定の値以下に維持することを目的として、所定の距離を伝搬するごとにレイを分割する手法が提案された。この手法の一例が、文献「Steven Fortune：“EFFICIENT ALGORITHMS FOR PREDICTION OF INDOOR RADIO PROPAGATION ”,in Proceedings of the 48th IEEE Vehicular Technology Conference,May,1998,pp.572-576」に記されている。図１６はこの文献の手法を説明する図である。
【００１２】
図１６の三角錐７１２は、割当てられたレイ７０２の送信アンテナ７０１からの伝搬距離が所定の値に達すると、三角錐７０８から三角錐７１１に細分される。細分された三角錐には、それぞれレイ７０４からレイ７０７が割当てられ、以降、同様な処理が繰返される。この結果、空間分解能は伝搬距離に関わらず常に一定値以下に保たれることになるのである。
【００１３】
図１７は観測領域０１７と、送信点０１５と、受信点０１６、２つの物体００１，００２とが与えられた場合のラウンチング法の動作について説明する図である。図１７では、簡単のため２次元平面に限定して動作の説明をするが、実際には３次元空間内で動作させることは勿論である。
【００１４】
離散的な間隔で伝搬経路００３の方向へ放射されたレイは、点０１８で物体００１へ入射し、その結果、反射レイ００５と透過レイ００４とが生成される。反射レイ００５は、さらに点０１９で物体００２へ入射し、再び反射レイ００６と透過レイ００７とが生成される。生成された反射レイ００６は受信点０１６の近傍を通過するため、レイ００６を受信点に到来したレイとみなし、経路００３，００５，００６からなる経路を送信点０１５と受信点０１６とを結ぶ伝搬経路の一つとみなす。
【００１５】
伝搬経路００３，００５，００６から求まる受信電力ならびに到来遅延時間は図１８のように記録される。図１８の横軸１０３はレイが送信点０１５から受信点０１６まで経路００３，００５，００６を経由して到達するのに要した遅延時間を示し、縦軸１０２は当該経路を経て到達したレイが有する電力強度を表す。透過レイ００４ならびに００７についても、伝搬経路００３，００５，００６と同様の反射・透過レイ探索を繰り返し、受信点０１６の近傍をレイが通過した場合には、伝搬経路００６の場合と同様に到来レイとして扱い、以上の処理を探索終了条件を満たすまで続ける。
【００１６】
探索終了条件は、反射・透過点での受信電界強度が所定値を下回った場合、または反射・透過の総数が所定の回数に達した場合などとする。送信点０１５から伝搬経路００３の方向へ放射されたレイの反射・透過経路探索が終了した後、方向００８から０１４のように他の放射方向へ向けてレイを放射した場合についても同様に反射・透過経路探索を行う。事前に定めた全ての方向へレイを放射し、以上の伝搬経路探索を行うことにより、受信点０１６における遅延プロファイルである図１９が得られる。図１９の横軸２０３は送信点０１５から受信点０１６へ到達するのに要した時間を示し、縦軸２０２は当該経路を通過したレイが有する電力強度を表す。
【００１７】
受信電力は図１９に示した全ての経路の電力強度の和で与えられ、歪の程度を示す遅延拡がりは、各遅延時間において受信電力に対する当該電力強度の電力比を当該遅延時間の出現確率として、遅延時間の標準偏差により与えられる。
【００１８】
【発明が解決しようとする課題】
上述したラウンチング法では、レイと近傍の部分空間が同一視されるために、レイが障壁と反射する場合は、レイ近傍の部分空間でも同様な反射が生じているとみなされる。そのため、推定精度の観点からは、反射点における空間分解能を当該障壁の面積に比べて大きくなり過ぎないように設定することが望ましい。しかし、推定精度の向上を目的として全体の空間分解能の上限値を各障壁との反射の際に要求される空間分解能の最小値に設定すると、電波伝搬推定に要する計算時間が増大するという問題点がある。その理由は、空間分解能を小さくすると放射するレイの総数が増大し、経路探索に要する全体の時間が増大するからである。
【００１９】
一方、推定に要する計算時間を抑えるために空間分解能の上限値を比較的大きく設定すると、全体の推定精度に劣化が生じるという問題点がある。その理由は、所要空間分解能が小さな障壁でレイが反射する場合には、実際には反射が生じていないレイ近傍の領域でもレイと同様に反射が生じているとみなすからである。
【００２０】
本発明は、ラウンチング法において伝搬推定の対象環境が与えられた場合に、推定精度の劣化を抑えながら所要推定時間を短縮化する手法を提供することを目的とする。
【００２１】
【課題を解決するための手段】
本発明によれば、３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御方法であって、前記遮蔽物と前記レイとの衝突点を検出する第一のステップと、前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う第二のステップと、この第二のステップで前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する第三のステップと、この第三のステップにおいて分割された前記部分空間の各々に対して同一視するレイを割当てる第四のステップとを含むことを特徴とするレイの空間分解能制御方法が得られる。
【００２３】
そして、前記第三のステップにおける分割数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは１〜Ｍ）の前記部分空間に対して前記第四のステップにおいてレイを割当てた後、割当てた各レイに対して再び前記第一から第四のステップを適用し、前記第二のステップにおいて当該レイと同一視される前記空間を分割しないと判定されるまで前記第一から第四のステップを、順次繰返すことを特徴とする。
【００２４】
また、前記第二のステップにおける前記判定基準は、前記遮蔽物との衝突点における前記レイの電界強度が所定の値以上であるか否かを調べる第一の判定基準と、前記遮蔽物の衝突面を含む平面による前記レイ近傍の前記空間の断面に対する前記遮蔽物の衝突面の面積比が所定値以下であるか否かを調べる第二の判定基準とからなることを特徴とし、更には、前記判定基準は、前記断面の面積が所定の値以上であり、かつ前記レイと前記衝突面の法線との成す角（入射角）が所定値以上であるか否かを調べる第三の判定基準を含むことを特徴とする。
【００２５】
また、錐体の中心線を前記レイの経路とすることにより前記レイ近傍の前記空間を前記錐体によって形成し、前記第三のステップにおいて、前記レイ近傍の空間を形成する前記錐体に相似な複数の錐体を用いて前記レイ近傍の空間を複数の部分空間に分割する際に、前記レイ近傍の空間を形成する前記錐体の体積と、この錐体を分割する前記錐体に相似な前記複数の錐体の和の体積との体積誤差が最小となるように、前記レイの近傍に前記複数の錐体を配置することを特徴とする。
【００２６】
更に、前記第三のステップにおいて、円錐に内接する前記錐体として三角錐を用い、前記三角錐に相似な複数の同一三角錐を前記三角錐内部に隙間無く配置できるように前記相似三角錐の大きさを設定することを特徴とし、また前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和の体積と、もとの前記正六角錐の体積誤差が最小となるように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００２７】
更にはまた、前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面がもとの前記正六角錐の側面に内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。また、前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面に、もとの前記正六角錐が内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００２９】
本発明によれば、３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御システムであって、前記遮蔽物と前記レイとの衝突点を検出する衝突点検出手段と、前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う判定手段と、この判定手段で前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する分割手段と、この分割手段において分割された前記部分空間の各々に対して同一視するレイを割当てる割当手段とを含むことを特徴とするレイの空間分解能制御システムが得られる。
【００３０】
そして、前記分割手段における分割数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは１〜Ｍ）の前記部分空間に対して前記割当手段においてレイを割当てた後、割当てた各レイに対して再び前記衝突点検出手段、前記判定手段、前記分割手段及び前記割当手段を実行制御し、前記判定手段において当該レイと同一視される前記空間を分割しないと判定されるまで、前記衝突点検出手段、前記判定手段、前記分割手段及び前記割当手段を、順次繰返し実行制御する制御手段とを含むことを特徴とする。
【００３１】
また、前記判定手段における前記判定基準は、前記遮蔽物との衝突点における前記レイの電界強度が所定の値以上であるか否かを調べる第一の判定基準と、前記遮蔽物の衝突面を含む平面による前記レイ近傍の前記空間の断面に対する前記遮蔽物の衝突面の面積比が所定値以下であるか否かを調べる第二の判定基準とからなることを特徴とし、更に、前記判定基準は、前記断面の面積が所定の値以上であり、かつ前記レイと前記衝突面の法線との成す角（入射角）が所定値以上であるか否かを調べる第三の判定基準を含むことを特徴とする。
【００３２】
また、錐体の中心線を前記レイの経路とすることにより前記レイ近傍の前記空間を前記錐体によって形成し、前記分割手段において、前記レイ近傍の空間を形成する前記錐体に相似な複数の錐体を用いて前記レイ近傍の空間を複数の部分空間に分割する際に、前記レイ近傍の空間を形成する前記錐体の体積と、この錐体を分割する前記錐体に相似な前記複数の錐体の和の体積との体積誤差が最小となるように、前記レイの近傍に前記複数の錐体を配置することを特徴とする。
【００３３】
更に、前記分割手段において、円錐に内接する前記錐体として三角錐を用い、前記三角錐に相似な複数の同一三角錐を前記三角錐内部に隙間無く配置できるように前記相似三角錐の大きさを設定することを特徴とし、また、前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和の体積と、もとの前記正六角錐の体積誤差が最小となるように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００３４】
また、前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面がもとの前記正六角錐の側面に内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。更に、前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面に、もとの前記正六角錐が内接するように前記相似六角錐の大きさを設定し配置することを特徴とする。
【００３６】
本発明によれば、３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御方法をコンピュータに実行させるためのプログラムであって、前記遮蔽物と前記レイとの衝突点を検出する第一の処理と、前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う第二の処理と、この第二の処理で前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する第三の処理と、この第三の処理において分割された前記部分空間の各々に対して同一視するレイを割当てる第四の処理とを含むことを特徴とするプログラムが得られる。
【００３７】
本発明の作用を述べる。本発明では、与えられた観測領域内をレイが伝搬する場合、その経路上に存在する障壁とレイが衝突する際に、所定の条件が満たされるまで当該レイ近傍の部分空間を複数の部分空間に分割する。従来の方法では、反射とは無関係に伝搬距離に応じて一律に部分空間を分割していたが、本発明の方法によれば、一律に扱った場合に推定精度の劣化を招くと考えられる部分空間は、劣化の影響が少なくなるまで分割し、細分化された各部分空間を個別に取り扱うことが可能となる。このような手法を採用することにより、電波伝搬の推定精度を従来方法よりも向上させることが可能となる。一方、推定精度の劣化の影響が小さいと考えられる場合は、当該部分空間を分割しないため、空間分解能の初期値を比較的大きめに設定しておけば、従来の方法において推定精度の向上を狙って一律に空間分解能を小さく設定した場合よりもレイの総数が低減し、伝搬推定に要する全体の計算時間を短縮させることができる。
【００３８】
これらの特徴により、本発明では、与えられた推定領域に応じて空間分解能を適応制御することが可能なことから、従来の方法に比べて推定精度の向上と計算時間の短縮を両立させることが可能である。
【００３９】
【発明の実施の形態】
以下に、図面を参照しながら本発明の実施例を詳細に説明する。図１は本発明の実施例の動作を示す流れ図であり、２つのルーチンＡ，Ｂ（８１８，８１９）よりなるものとする。図１では、３次元空間上に観測領域が与えられ、当該観測領域内に壁や什器等の物体が与えられ、当該観測領域内の点Ｐから当該伝搬方向へ放射されるレイが与えられることを前提とし（ステップ８０２）、反射点における反射・透過レイを生成するまでの一連の流れを示す。
【００４０】
ここで、点Ｐとは、ラウンチング法における反射・透過点もしくは送信アンテナの位置を指し、具体的には従来技術の説明において参照した図１７における点０１８、点０１９、送信アンテナ０１５等を指すものとする。また、当該伝搬方向とは、点Ｐが送信アンテナである場合は、送信アンテナの周囲に設けた離散的な放射方向を指し、具体的には従来の技術で説明したように図１７の方向００３、００８〜０１４等、もしくは図１４の送信アンテナ５０１と多面体の各頂点を結ぶ方向、例えばレイ５０４の方向等を指すものとする。一方、点Ｐが反射・透過点である場合は、反射レイもしくは透過レイの方向を指し、具体的には図１７における伝搬経路００３〜００７等を指すものとする。
【００４１】
まず、点Ｐから当該伝搬方向へ放射したレイが観測領域内の物体と衝突する点、すなわち反射点を探索し検出する（ステップ８０３，８０４）。次に、反射点に入射する当該レイと同一視されている部分空間を、より小さな部分空間に細分化するか否かの判定を事前に定めた判定基準（後述）に基づいて行う（ステップ８０５）。
【００４２】
当該部分空間をより小さな部分空間に細分する場合は、事前に定めた方法（後述）に基づいて当該部分空間をＭ個（Ｍは整数）の部分空間に分割し（ステップ８０６）、ルーチンＢ（８１９）へ移行する（ステップ８０７）。当該部分空間の細分化を行わない場合には、当該反射点における反射レイと透過レイを生成し終了する（ステップ８０８，８０９）。また、ステップ８０４で反射点が検出されなければ終了となる。
【００４３】
ルーチンＢ（８１９）では、まず内部変数ｉ（ｉ＝１〜Ｍ）を定義し、初期値としてｉ＝１を与える（ステップ８１１）。次に、ルーチンＡ（８１８）のステップ８０６でＭ個に分割された部分空間のうちｉ番目の部分空間を選択する（ステップ８１２）。そして、選択したｉ番目の部分空間にレイを割当て（ステップ８２０）、そのレイの伝搬方向を事前に定めた方法（後述）に基づいて算出する（ステップ８１３）。
【００４４】
ここで、ステップ８０６で細分化された部分空間にステップ８２０で割当てられたレイが、ステップ８１３で求めた当該伝搬方向へ起点Ｐから放射されることを前提とし、再度ルーチンＡ（８１８）へ移行する。ルーチンＡ（８１８）では、細分化された部分空間をさらに細分化するかどうかの判定を行い、細分化する場合は、ステップ８０５で再度細分化処理を行い、ルーチンＢ（８１９）へ移行する。細分化を行わない場合には、反射・透過レイを生成して（ステップ８０８）終了する（ステップ８０９）。
【００４５】
ステップ８１４でルーチン８１８へ移行し、ステップ８０９でルーチン８１８を終了した場合は、その直後にステップ８１５へ移行する。ステップ８１５では、内部変数ｉの値がＭより小さいか否かの判定を行い、Ｍよりも小さい場合は、ｉに１を加算して（ステップ８１７）他の部分空間を選択する（ステップ８１２）。ｉの値がＭ以上である場合は、ルーチン８１９を終了する（ステップ８１６）。
【００４６】
ステップ８０７でルーチン８１９へ移行し、ステップ８１６でルーチン８１９を終了した場合は、その直後にステップ８０９へ移行し、入れ子となっているルーチン８１８の一つを終了する。以上の処理を、ステップ８０６で部分空間を細分化する必要がなくなるまで実行する。
【００４７】
図１に示した発明の実施例の具体的動作例を、図２〜図５を用いて説明する。図２は部分空間と障壁の交差面の生成過程を説明する図である。点９０１から放射されたレイ９０２は、点９０３において障壁９０６と衝突する。線９０８は障壁９０６に対する法線であり、角度９０７は障壁９０６に対するレイ９０２の入射角を指す。このとき、レイ９０２近傍の部分空間９０４と障壁９０６とによって、交差面９０５が生じる。図２はアンテナ周囲の空間を三角錐によって分割した場合を示す。
【００４８】
図３は本発明における部分空間の細分化を説明する図であり、図２においてθ＝０の場合を示している。図３１（ａ）は、レイ９０２が点１２１３へ入射した結果、障壁１２１２と部分空間の断面１２１１とによって生じる交差面１２１２（斜線部分）を示している。点９０３において、部分空間９０４を細分するか否かの判定基準（図１のフローのステップ８０５の判定基準）の一例として、以下の項目が挙げられる。
【００４９】
（１）点１２１３における当該レイの電界強度が所定値以上である：
（２）部分空間断面１２１１に対する交差面１２１２の面積比が所定値以下である：
（３）部分空間断面１２１１の面積が所定の値以上であり、かつレイ９０２の入射角９０７が所定値以上である：
基本的には、上記の３項目の基準のうち（１）及び（２）を同時に満たす場にを部分空間９０４の分割（細分）を行うものとするが、これ等２つの基準に加えて、（３）の基準を更に付加して、３つの基準を同時に満たす場合に部分空間の分割（細分）を行うものとしても良い。
【００５０】
細分した部分空間に割当てられるレイの伝搬方向の決定方法（図１のフローのステップ８１３の方向算出法）としては、例えば以下の方法が挙げられる。まず、点９０３と点９０１とを結ぶ線分上の任意の点においてレイ９０２に対し垂直な面で部分空間を切断する。その結果生じる、図４のような部分空間の断面１３０９をＭ個（この例では、Ｍ＝４）の相似な領域１３０５〜１３０８に分割し、点１３０１〜１３０４のように各領域の重心点を求める。同様な重心点を部分空間の別の切断面についても求める。対向する重心点同士を結ぶと、図５のように所望の伝搬方向ベクトル１９０１が得られる。
【００５１】
なお、図３〜図５はＭ＝４の場合を示したものであるが、Ｍが他の値の場合も同様である。
【００５２】
このように部分空間９０４を細分すると、反射点９０３における部分空間の断面は図３（ｂ）の様になる。図９のフロー及び上記の分割判定基準に基づくと、例えば、部分空間１２０２はこれ以上分割されずに、部分空間１２０２に割当てられたレイは障壁１２１２で反射することなく、点１２０６を通過して他の障壁に衝突するまで進行を続ける。この部分空間１２０２が分割されないのは、図１のフローにおけるステップ８０４で反射点が検出されないからであり、詳述すると、この部分空間１２０２に割当てられたレイは、点１２０６を通過してしまい、障壁１２１２で反射しないからである。
【００５３】
一方、部分空間１２０３は上記基準を満たしているとして、図３（ｃ）のように取出された後、さらに細分され、以後同様に必要なだけ細分化が行われる。
【００５４】
次に、本実施例の効果について説明する。本実施例では、当該部分空間と障壁との交差面と反射点における当該部分空間の断面の面積比を分割判定の基準とし、実際には反射が生じていない領域を分割して別扱いにすることができる（上記判定基準の（２））。また、図６に示すとおり、電波の反射率は一般に入射角が大きくなるにつれて急増することが多いが、当該レイの入射角度を判定基準とすることにより、本来は異なる角度で障壁に入射し反射率が大きく異なる領域を分割して別扱いにすることができる（上記判定基準の（３））。
【００５５】
これらの作用によって、ラウンチング法よる推定精度を向上させることができることになるのである。一方、本実施の形態では細分した後の部分空間はもとの部分空間と相似であるため、以降の細分化処理においても同様な処理を繰り返し適用することができ、処理を簡易化することが可能である。さらに、分割前の部分領域は分割後の部分空間の和で与えられるために、アンテナ周囲の全空間を隙間無く推定することが可能である。
【００５６】
以上述べた実施例の処理動作から、この処理動作を実現するための機能ブロックが、図７に示す様に得られることになる。図７を参照すると、衝突検出部１と、部分空間分割判定部２と、部分空間分割部３と、部分空間へのレイ割当部４と、制御部（ＣＰＵ）５と、ＲＡＭ（ランダムアクセスメモリ）６と、ＲＯＭ（リードオンリメモリ）７とを含んで構成される。
【００５７】
衝突検出部１はレイと遮蔽物（障壁）との衝突を検出する機能を有しており、例えば、図３の例では、障壁１２１２に対してレイが衝突する関係にあり、これは障壁１２１２の面を表す３次元空間上に規定される式と、同じくレイを表す３次元空間上に規定される式とを連立させて解く手法を採用することで、衝突の有無が判定可能である。
【００５８】
部分空間判定部２はレイと同一視される空間を分割するか否かの判定を、上述した判定基準に従って判定する機能を有しており、部分空間分割部３は上述した予め定められた規定に従って当該空間を複数の部分空間に分割する機能を有する。部分空間へのレイ割当部４は分割された部分空間の各々に対して同一視するレイを割当てる機能を有する。
【００５９】
制御部５はＣＰＵ（コンピュータの中央演算処理部に相当）からなり、ＲＯＭ７の記録媒体に予め格納されている動作制御プログラムを読込んで、このプログラムの制御手順に従って、作業用メモリであるＲＡＭ６を使用しつつ各部１〜４を制御し、図１に示した動作処理を実行制御するものである。すなわち、そのプログラムがＣＰＵの動作を制御するものであり、プログラム制御されたそのＣＰＵが、当該プログラムにより指令されて、上記の各処理を行うようになっている。
【００６０】
次に、本発明の第２の実施例について、図面を参照して説明する。第１の実施例においては、レイに垂直な平面で部分空間を切断したときに、図４に示す形態となるようにアンテナ周囲の空間を錐体で分割していた。これに対して、第２の実施例では、部分空間の断面が図８となるようにアンテナ周囲の空間を分割している点で異なる。図８では、正６角形１４０８と相似な正６角形１４０９〜１４１５で部分空間を細分化するＭ＝７の場合を示している。正６角形１４０９〜１４１４はそれぞれ辺の中点１４１６〜１４２１で正６角形１４０８の頂点と接しており、分割後に各部分空間に割当てられるレイは正６角形１４０９〜１４１５の重心１４０１〜１４０７を通る。
【００６１】
この第２の実施例で用いた正６角形のように比較的円形に近い図形を用いた場合、受信アンテナ近傍をレイが通過した場合にそのレイを取込むか否かの判定を簡易化でき処理速度を向上させることができる。その理由は、一般にこの判定が、当該レイ近傍の部分空間の内部に受信アンテナが含まれるか否かを調べることによって行われることが多いためである。そこで、当該部分空間を図９のように正６角形１６０２の外接円１６０１を当該レイに垂直な切断面とする部分空間で近似することによって、この判定は受信アンテナと当該レイの間の距離を調べることに帰着でき、判定処理の簡易化が達成できる。
【００６２】
次に、本発明の第３の実施例について図面を参照して説明する。当該部分空間の断面が図８となるようにアンテナ周囲の空間を分割した場合、分割前の断面１４０８と分割後の部分空間の和領域を比較すると、図８で斜線で示した部分のような過不足領域が生じる。そこで、分割後の部分空間が分割前の部分空間の領域外にはみ出さないような分割形態を図１０に示す。
【００６３】
一方、分割後の部分空間と分割前の部分空間との間に隙間が生じないように設定した分割形態を図１１に示す。図１０では、頂点１８０１〜１８１２が正６角形１８１３に内接するように、分割後の正６角形の大きさを設定している。一方、図１１では、分割前の正６角形の各辺の中点１５０１〜１５０６において、分割後の正６角形と交差するように、分割後の正６角形の大きさを設定している。
【００６４】
特に、図１１の分割形態を用いた場合は、第１の実施例と同様に、アンテナ周囲の空間の電波伝搬状況を隙間なく推定することができる。
【００６５】
部分空間分割の例として上述の様に種々示しているが、要は、円錐に内接する錐体のうち、この円錐との体積誤差が小さい錐体を設け、その中心線をレイの経路とすることによりレイ近傍の空間を当該錐体によって形成し、レイ近傍の空間を形成するこの錐体に相似な複数の錐体を用いてレイ近傍の空間を複数の部分空間に分割する際に、レイ近傍の空間を形成する錐体と、この錐体を分割する錐体に相似な複数の錐体の和との体積誤差が小さくなるように、レイの近傍に複数の錐体を配置すれば良い。
【００６６】
【発明の効果】
以上述べた様に、本発明によれば、従来の方法に比べてラウンチング法におけるレイの密度を柔軟に制御し、推定精度の向上と計算時間の短縮を効率よく両立できるという効果がある。その理由は、従来の方法では伝搬距離に応じて一律にレイの密度を制御していたのに対し、本発明では、反射・透過の際にレイの密度が過小で推定精度の劣化を招くと考えられる場合に、劣化が小さいと判断されるまでレイの密度を増大させるためである。これにより、推定精度を損なわない範囲で不必要なレイの生成による計算時間の増大を抑えることができる。
【図面の簡単な説明】
【図１】本発明の実施例の動作を示すフロー図である。
【図２】部分空間と障壁の交差面の生成過程を説明する図である。
【図３】本発明における部分空間の細分化を説明する図である。
【図４】本発明の第１の実施例における部分空間の分割例を示す断面図である。
【図５】分割された部分空間に割当てられたレイの伝搬方向を求める過程を説明する図である。
【図６】反射率の入射角特性の例を示す図である。
【図７】本発明の機能ブロックの例を示す図である。
【図８】本発明の第２の実施例における部分空間の分割例を示す断面図である。
【図９】部分空間の円錐近似を説明する図である。
【図１０】本発明の第３の実施例における部分空間の分割の一例を示す断面図である。
【図１１】本発明の第３の実施例における部分空間の分割の他の例を示す断面図である。
【図１２】送信アンテナ周囲の閉領域の例を示す図である。
【図１３】送信アンテナ周囲の閉領域の分割過程を説明する図である。
【図１４】離散的間隔によるレイの放射を説明する図である。
【図１５】レイと部分空間の例を示す図である。
【図１６】従来の手法による部分空間の分割を説明する図である。
【図１７】ラウンチング法の説明図である。
【図１８】一本のレイに対する観測点におけるパスプロファイルの例を示す図である。
【図１９】ラウンチング法によって得られる観測点におけるパスプロファイルの例を示す図である。
【符号の説明】
１ 衝突検出部
２ 部分空間分割判定部
３ 部分空間分割部
４ 部分空間へのレイ割当部
５ 制御部
６ ＲＡＭ
７ ＲＯＭ
００１，００２ 物体（遮蔽物）
００３〜０１４ レイ伝搬方向
０１５ 送信アンテナ位置
０１６ 受信点
０１８，０１９，９０１，９０３ 反射・透過点
９０２，１９０１ レイ
９０４ 部分空間
９０５ 交差面
９０６ 障壁
Claims (20)
- ３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御方法であって、 前記遮蔽物と前記レイとの衝突点を検出する第一のステップと、
前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う第二のステップと、
この第二のステップで前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する第三のステップと、
この第三のステップにおいて分割された前記部分空間の各々に対して同一視するレイを割当てる第四のステップと、
を含むことを特徴とするレイの空間分解能制御方法。 - 前記第三のステップにおける分割数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは１〜Ｍ）の前記部分空間に対して前第四のステップにおいてレイを割当てた後、割当てた各レイに対して再び前記第一から第四のステップを適用し、前記第二のステップにおいて当該レイと同一視される前記空間を分割しないと判定されるまで前記第一から第四のステップを、順次繰返すことを特徴とする請求項１記載のレイの空間分解能制御方法。
- 前記第二のステップにおける前記判定基準は、前記遮蔽物との衝突点における前記レイの電界強度が所定の値以上であるか否かを調べる第一の判定基準と、前記遮蔽物の衝突面を含む平面による前記レイ近傍の前記空間の断面に対する前記遮蔽物の衝突面の面積比が所定値以下であるか否かを調べる第二の判定基準とからなることを特徴とする請求項１または２記載のレイの空間分解能制御方法。
- 前記第二のステップにおける前記判定基準は、前記断面の面積が所定の値以上であり、かつ前記レイと前記衝突面の法線との成す角（入射角）が所定値以上であるか否かを調べる第三の判定基準を更に含むことを特徴とする請求項３記載のレイの空間分解能制御方法。
- 錐体の中心線を前記レイの経路とすることにより前記レイ近傍の前記空間を前記錐体によって形成し、
前記第三のステップにおいて、前記レイ近傍の空間を形成する前記錐体に相似な複数の錐体を用いて前記レイ近傍の空間を複数の部分空間に分割する際に、前記レイ近傍の空間を形成する前記錐体の体積と、この錐体を分割する前記錐体に相似な前記複数の錐体の和の体積との体積誤差が最小となるように、前記レイの近傍に前記複数の錐体を配置することを特徴とする請求項１〜４いずれか記載のレイの空間分解能制御方法。 - 前記第三のステップにおいて、円錐に内接する前記錐体として三角錐を用い、前記三角錐に相似な複数の同一三角錐を前記三角錐内部に隙間無く配置できるように前記相似三角錐の大きさを設定することを特徴とする請求項１〜４いずれか記載のレイの空間分解能制御方法。
- 前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和の体積と、もとの前記正六角錐の体積誤差が最小となるように前記相似六角錐の大きさを設定し配置することを特徴とする請求項１〜４いずれか記載のレイの空間分解能制御方法。
- 前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その前記相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面がもとの前記正六角錐の側面に内接するように前記相似六角錐の大きさを設定し配置することを特徴とする請求項１〜４いずれか記載のレイの空間分解能制御方法。
- 前記第三のステップにおいて、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面に、もとの前記正六角錐が内接するように前記相似六角錐の大きさを設定し配置することを特徴とする請求項１〜４いずれか記載のレイの空間分解能制御方法。
- ３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御システムであって、
前記遮蔽物と前記レイとの衝突点を検出する衝突点検出手段と、
前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う判定手段と、
この判定手段で前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する分割手段と、 この分割手段において分割された前記部分空間の各々に対して同一視するレイを割当てる割当手段と、
を含むことを特徴とするレイの空間分解能制御システム。 - 前記分割手段における分割数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは１〜Ｍ）の前記部分空間に対して前記割当手段においてレイを割当てた後、割当てた各レイに対して再び前記衝突点検出手段、前記判定手段、前記分割手段及び前記割当手段を実行制御し、前記判定手段において当該レイと同一視される前記空間を分割しないと判定されるまで、前記衝突点検出手段、前記判定手段、前記分割手段及び前記割当手段を、順次繰返し実行制御する制御手段とを含むことを特徴とする請求項１０記載のレイの空間分解能制御システム。
- 前記判定手段における前記判定基準は、前記遮蔽物との衝突点における前記レイの電界強度が所定の値以上であるか否かを調べる第一の判定基準と、前記遮蔽物の衝突面を含む平面による前記レイ近傍の前記空間の断面に対する前記遮蔽物の衝突面の面積比が所定値以下であるか否かを調べる第二の判定基準とからなることを特徴とする請求項１０または１１記載のレイの空間分解能制御システム。
- 前記判定手段における前記判定基準は、前記断面の面積が所定の値以上であり、かつ前記レイと前記衝突面の法線との成す角（入射角）が所定値以上であるか否かを調べる第三の判定基準を更に含むことを特徴とする請求項１２記載のレイの空間分解能制御システム。
- 錐体の中心線を前記レイの経路とすることにより前記レイ近傍の前記空間を前記錐体によって形成し、
前記分割手段において、前記レイ近傍の空間を形成する前記錐体に相似な複数の錐体を用いて前記レイ近傍の空間を複数の部分空間に分割する際に、前記レイ近傍の空間を形成する前記錐体の体積と、この錐体を分割する前記錐体に相似な前記複数の錐体の和の体積との体積誤差が最小となるように、前記レイの近傍に前記複数の錐体を配置することを特徴とする請求項１０〜１３いずれか記載のレイの空間分解能制御システム。 - 前記分割手段において、円錐に内接する前記錐体として三角錐を用い、前記三角錐に相似な複数の同一三角錐を前記三角錐内部に隙間無く配置できるように前記相似三角錐の大きさを設定することを特徴とする請求項１０〜１３いずれか記載のレイの空間分解能制御システム。
- 前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和の体積と、もとの前記正六角錐の体積誤差が最小となるように前記相似六角錐の大きさを設定し配置することを特徴とする請求項１０〜１３いずれか記載のレイの空間分解能制御システム。
- 前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面がもとの前記正六角錐の側面に内接するように前記相似六角錐の大きさを設定し配置することを特徴とする請求項１０〜１３いずれか記載のレイの空間分解能制御システム。
- 前記分割手段において、円錐に内接する前記錐体として正六角錐を用い、前記正六角錐に相似な正六角錐を前記正六角錐の中心に配置し、その相似正六角錐の各側面に前記相似正六角錐と同一な６個の正六角錐を各々側面が接するように配置し、全７個の前記相似正六角錐の和からなる錐体の最外郭側面に、もとの前記正六角錐が内接するように前記相似六角錐の大きさを設定し配置することを特徴とする請求項１０〜１３いずれか記載のレイの空間分解能制御システム。
- ３次元空間上に規定される観測領域内に複数の物体が与えられ、前記観測領域内に電波線（レイ）が与えられると共に、前記レイの近傍に前記レイと同一視される空間が与えられ、前記物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突した際に、予め定められた基準に基づき前記レイの空間分解能の調節を行う電波伝搬特性推定におけるレイの空間分解能制御方法をコンピュータに実行させるためのプログラムであって、
前記遮蔽物と前記レイとの衝突点を検出する第一の処理と、
前記レイと同一視される前記空間を分割するか否かの判定を、所定判定基準に従って行う第二の処理と、
この第二の処理で前記レイと同一視される前記空間を分割すると判定された場合に、前記レイと同一視される前記空間を複数の部分空間に分割する第三の処理と、
この第三の処理において分割された前記部分空間の各々に対して同一視するレイを割当てる第四の処理と、
を含むことを特徴とするプログラム。 - 前記第三の処理における分割数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは１〜Ｍ）の前記部分空間に対して前記第四のステップにおいてレイを割当てた後、割当てた各レイに対して再び前記第一から第四の処理を適用し、前記第二の処理において当該レイと同一視される前記空間を分割しないと判定されるまで前記第一から第四の処理を、順次繰返すことを特徴とする請求項１９記載のプログラム。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001008368A JP4273665B2 (ja) | 2001-01-17 | 2001-01-17 | レイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラム |
| US10/041,575 US6876851B2 (en) | 2001-01-17 | 2002-01-10 | Method and system for radio wave propagation characteristics estimation and ray spatial resolution control |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001008368A JP4273665B2 (ja) | 2001-01-17 | 2001-01-17 | レイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002214267A JP2002214267A (ja) | 2002-07-31 |
| JP4273665B2 true JP4273665B2 (ja) | 2009-06-03 |
Family
ID=18875999
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001008368A Expired - Fee Related JP4273665B2 (ja) | 2001-01-17 | 2001-01-17 | レイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラム |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US6876851B2 (ja) |
| JP (1) | JP4273665B2 (ja) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7634265B2 (en) | 2004-03-17 | 2009-12-15 | Nec Corporation | Radio wave propagation characteristic estimation system, and its method and program |
| JP4481840B2 (ja) * | 2005-01-20 | 2010-06-16 | キヤノン株式会社 | 情報処理方法および装置 |
| JP4530898B2 (ja) * | 2005-04-01 | 2010-08-25 | 株式会社日立製作所 | 電波伝搬の推定プログラム、電波伝搬の推定方法、この方法を実行する装置 |
| US7962102B2 (en) | 2006-11-30 | 2011-06-14 | Motorola Mobility, Inc. | Method and system for adaptive ray launching |
| JP5120786B2 (ja) * | 2007-02-16 | 2013-01-16 | 日本電気株式会社 | 電波伝搬特性推定システム及びその方法並びにプログラム |
| US8332196B2 (en) * | 2007-11-30 | 2012-12-11 | Motorola Mobility Llc | Method and apparatus for enhancing the accuracy and speed of a ray launching simulation tool |
| US20090167756A1 (en) * | 2007-12-31 | 2009-07-02 | Motorola, Inc. | Method and apparatus for computation of wireless signal diffraction in a three-dimensional space |
| JP5509666B2 (ja) * | 2008-05-08 | 2014-06-04 | 日本電気株式会社 | 電波伝搬特性推測支援システム、電波伝搬特性推測支援方法及び電波伝搬特性推測支援装置 |
| US8208916B2 (en) * | 2009-12-16 | 2012-06-26 | Electronics And Telecommunications Research Institute | Apparatus and method of predicting radio wave environment |
| JP2019114906A (ja) * | 2017-12-22 | 2019-07-11 | ルネサスエレクトロニクス株式会社 | 半導体装置及び半導体システム |
| US11398072B1 (en) * | 2019-12-16 | 2022-07-26 | Siemens Healthcare Gmbh | Method of obtaining a set of values for a respective set of parameters for use in a physically based path tracing process and a method of rendering using a physically based path tracing process |
| CN114928420B (zh) * | 2022-05-23 | 2023-08-29 | 南京捷希科技有限公司 | 一种基于正向射线追踪的无线信道多径仿真精度优化方法 |
| CN115102650A (zh) * | 2022-06-20 | 2022-09-23 | 南京邮电大学 | 一种新型电波传播追踪方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5450615A (en) * | 1993-12-22 | 1995-09-12 | At&T Corp. | Prediction of indoor electromagnetic wave propagation for wireless indoor systems |
| US5623429A (en) * | 1994-04-06 | 1997-04-22 | Lucent Technologies Inc. | Techniques for expeditiously predicting electromagnetic wave propagation |
| US5949988A (en) * | 1996-07-16 | 1999-09-07 | Lucent Technologies Inc. | Prediction system for RF power distribution |
| US6119009A (en) * | 1997-09-18 | 2000-09-12 | Lucent Technologies, Inc. | Method and apparatus for modeling the propagation of wireless signals in buildings |
-
2001
- 2001-01-17 JP JP2001008368A patent/JP4273665B2/ja not_active Expired - Fee Related
-
2002
- 2002-01-10 US US10/041,575 patent/US6876851B2/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| US6876851B2 (en) | 2005-04-05 |
| US20020094809A1 (en) | 2002-07-18 |
| JP2002214267A (ja) | 2002-07-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4273665B2 (ja) | レイの空間分解能制御方法及びレイの空間分解能制御システム並びにプログラム | |
| KR100288560B1 (ko) | 도심지 협곡 모델에서 편파 효과를 고려한 전파의 전파특성 예측 방법 | |
| US7962102B2 (en) | Method and system for adaptive ray launching | |
| JPH10503915A (ja) | 無線通信システム計画のための方法 | |
| CN102546039B (zh) | 一种无线电磁波传播预测方法及装置 | |
| JPH07288495A (ja) | 電磁波伝播を迅速に予測するための技法 | |
| EP2252099B1 (en) | Apparatus and method for searching for wireless station locations | |
| JP3823917B2 (ja) | レイの受信判定方法並びにシステム及びそれを用いた電波伝搬特性推定方法 | |
| Mohtashami et al. | Efficient shooting and bouncing ray tracing using decomposition of wavefronts | |
| Kim et al. | Accelerated three dimensional ray tracing techniques using ray frustums for wireless propagation models | |
| US8358305B2 (en) | Method for 3D preprocessing for ray tracing | |
| CN109861775B (zh) | 一种传播路径搜索方法以及装置 | |
| JP5301521B2 (ja) | 電波推定方法、及び電波推定プログラム、並びに電波推定装置 | |
| Bhatia et al. | Efficient ray-tracing channel emulation in industrial environments: An analysis of propagation model impact | |
| JP3654197B2 (ja) | 電波伝搬特性予測システム及びその方法並びにプログラム | |
| Hussain et al. | An image visibility based pre-processing method for fast ray tracing in urban environments | |
| KR100329518B1 (ko) | 무선망 설계시스템에서의 전파해석방법 | |
| CN112199865B (zh) | 基于sbr的多层密闭舱体射线跟踪方法、存储介质及终端 | |
| Kreuzgruber et al. | Prediction of indoor radio propagation with the ray splitting model including edge diffraction and rough surfaces | |
| JP2002107397A (ja) | 電波伝搬特性予測における遮蔽物の検出方法及びシステム並びに記録媒体 | |
| TW202133633A (zh) | 細胞佈建之方法與相關的電腦程式產品 | |
| KR100205957B1 (ko) | 전기 영상 트리를 이용한 전파 경로 추적 방법 | |
| KR19990057146A (ko) | 다중 수신점 동시 전파해석에 의한 광선추적 기법의 속도 개선방법 | |
| KR100614408B1 (ko) | 광선 구조체를 이용한 점대점 전파해석 방법 | |
| EP2249594B1 (en) | Estimation of the shadowing component in a wireless communication system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061013 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081021 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20081222 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090210 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090223 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120313 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120313 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130313 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130313 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140313 Year of fee payment: 5 |
|
| LAPS | Cancellation because of no payment of annual fees |