JP2002214267A - 電波伝搬特性推定方法及びレイの空間分解能制御方法並びにシステム及びプログラム - Google Patents
電波伝搬特性推定方法及びレイの空間分解能制御方法並びにシステム及びプログラムInfo
- Publication number
- JP2002214267A JP2002214267A JP2001008368A JP2001008368A JP2002214267A JP 2002214267 A JP2002214267 A JP 2002214267A JP 2001008368 A JP2001008368 A JP 2001008368A JP 2001008368 A JP2001008368 A JP 2001008368A JP 2002214267 A JP2002214267 A JP 2002214267A
- Authority
- JP
- Japan
- Prior art keywords
- ray
- regular hexagonal
- pyramid
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 75
- 230000004304 visual acuity Effects 0.000 title 1
- 238000001514 detection method Methods 0.000 claims description 6
- 230000005684 electric field Effects 0.000 claims description 6
- 230000000903 blocking effect Effects 0.000 claims 1
- 230000004888 barrier function Effects 0.000 abstract description 22
- 238000007796 conventional method Methods 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 7
- 230000005855 radiation Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 230000000644 propagated effect Effects 0.000 abstract description 2
- 238000004904 shortening Methods 0.000 abstract description 2
- 230000000593 degrading effect Effects 0.000 abstract 4
- 230000005540 biological transmission Effects 0.000 description 33
- 238000010586 diagram Methods 0.000 description 19
- 238000012545 processing Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 10
- 230000006866 deterioration Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000005484 gravity Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000007619 statistical method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 230000006872 improvement Effects 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
- 230000009467 reduction Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003936 working memory Effects 0.000 description 1
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)
- Radar Systems Or Details Thereof (AREA)
- Mobile Radio Communication Systems (AREA)
- Image Generation (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
(57)【要約】
【課題】 ラウンチング法において伝搬推定の対象環境
が与えられた場合に、推定精度の劣化を抑えながら所要
推定時間を短縮化する。 【解決手段】 与えられた観測領域内をレイが伝搬する
場合、その経路上に存在する障壁と衝突する際に、所定
の条件が満たされるまで当該レイ近傍の空間を複数の部
分空間に分割する。推定精度の劣化を招くと考えられる
レイ近傍の空間は、反射の際に、劣化の影響が少なくな
るまで分割し、細分化された各部分空間を個別に取扱う
ことが可能となり、電波伝搬の推定精度を従来方法より
も向上させることができる。一方、推定精度の劣化の影
響が小さいと考えられる場合は、空間分割を行わないた
め、送信アンテナから粗くレイを放射すれば、従来方法
よりもレイの総数が低減し、伝搬推定に要する全体の計
算時間を短縮することができる。
が与えられた場合に、推定精度の劣化を抑えながら所要
推定時間を短縮化する。 【解決手段】 与えられた観測領域内をレイが伝搬する
場合、その経路上に存在する障壁と衝突する際に、所定
の条件が満たされるまで当該レイ近傍の空間を複数の部
分空間に分割する。推定精度の劣化を招くと考えられる
レイ近傍の空間は、反射の際に、劣化の影響が少なくな
るまで分割し、細分化された各部分空間を個別に取扱う
ことが可能となり、電波伝搬の推定精度を従来方法より
も向上させることができる。一方、推定精度の劣化の影
響が小さいと考えられる場合は、空間分割を行わないた
め、送信アンテナから粗くレイを放射すれば、従来方法
よりもレイの総数が低減し、伝搬推定に要する全体の計
算時間を短縮することができる。
Description
【０００１】
【発明の属する技術分野】本発明は電波伝搬特性推定方
法及びレイの空間分解能制御方法並びにシステム及びプ
ログラムに関し、特に幾何光学的手法による電波伝搬特
性予測における電波伝搬特性推定方式及びレイの空間分
解能制御方式に関するものである。
法及びレイの空間分解能制御方法並びにシステム及びプ
ログラムに関し、特に幾何光学的手法による電波伝搬特
性予測における電波伝搬特性推定方式及びレイの空間分
解能制御方式に関するものである。
【０００２】
【従来の技術】無線通信システムにおける基地局や親機
等の配置を支援するために電波伝搬シミュレータが用い
られる。電波伝搬シミュレータによって任意の受信点に
おける受信電力や遅延拡がりを評価できる。この情報を
もとに効果的な送信局の設置場所を決定し、配置が必要
な基地局数を削減する等の効率化が達成される。
等の配置を支援するために電波伝搬シミュレータが用い
られる。電波伝搬シミュレータによって任意の受信点に
おける受信電力や遅延拡がりを評価できる。この情報を
もとに効果的な送信局の設置場所を決定し、配置が必要
な基地局数を削減する等の効率化が達成される。
【０００３】電波伝搬シミュレーションは統計的手法と
決定論的手法とに大別できる。統計的手法は、距離や周
波数などを引数とする伝搬損失推定式を用いる手法であ
り、推定式中のパラメータは、伝搬損失の実測定で得ら
れた多数のデータをもとに多変量解析等の統計的手法に
より決定される。一方、決定論的手法は、アンテナから
放射される電波を多数の電波線（レイ）の集まりとして
表現する手法であり、各レイは幾何光学的に反射・透過
を繰り返しながら伝搬するものとする。観測点における
伝搬損および遅延量は、観測点に到達する各レイの電力
と伝搬時間を合成することにより求めることができる。
決定論的手法とに大別できる。統計的手法は、距離や周
波数などを引数とする伝搬損失推定式を用いる手法であ
り、推定式中のパラメータは、伝搬損失の実測定で得ら
れた多数のデータをもとに多変量解析等の統計的手法に
より決定される。一方、決定論的手法は、アンテナから
放射される電波を多数の電波線（レイ）の集まりとして
表現する手法であり、各レイは幾何光学的に反射・透過
を繰り返しながら伝搬するものとする。観測点における
伝搬損および遅延量は、観測点に到達する各レイの電力
と伝搬時間を合成することにより求めることができる。
【０００４】決定論的手法は、伝搬経路を追跡する方法
の違いから、更にイメージング法とラウンチング法とに
大別できる。イメージング法は送受信点間を結ぶレイの
反射・透過経路を反射面に対する鏡映点を求めて決定す
る手法である。反射・透過経路は、送受信点と反射・透
過障壁の位置が決まれば一意に求まるために、イメージ
ング法は厳密なレイの伝搬経路を探索する手法である。
の違いから、更にイメージング法とラウンチング法とに
大別できる。イメージング法は送受信点間を結ぶレイの
反射・透過経路を反射面に対する鏡映点を求めて決定す
る手法である。反射・透過経路は、送受信点と反射・透
過障壁の位置が決まれば一意に求まるために、イメージ
ング法は厳密なレイの伝搬経路を探索する手法である。
【０００５】一方、ラウンチング法は、受信点の位置に
関わり無く、アンテナから離散的な角度間隔でレイを放
射し、反射・透過等を繰返しながら受信点近傍を通過す
るレイを当該受信点に到達したレイとみなす手法であ
る。このラウンチング法は、送受信点間を結ぶレイの伝
搬経路の解をイメージング法のように厳密に求めるので
はなく、近似的に求めるために、伝搬経路探索に要する
時間を短縮できるという特徴がある。
関わり無く、アンテナから離散的な角度間隔でレイを放
射し、反射・透過等を繰返しながら受信点近傍を通過す
るレイを当該受信点に到達したレイとみなす手法であ
る。このラウンチング法は、送受信点間を結ぶレイの伝
搬経路の解をイメージング法のように厳密に求めるので
はなく、近似的に求めるために、伝搬経路探索に要する
時間を短縮できるという特徴がある。
【０００６】ラウンチング法では、まず、送信アンテナ
の周囲に閉領域を定め、その閉領域を部分空間に分割す
る。次に、各部分空間に対して１本のレイを割当て、部
分空間と割当てたレイを同一視する。その後、送信アン
テナから放射される有限な本数のレイの伝搬経路を追跡
し、その結果を基にアンテナの周囲全空間の電波伝搬状
況を推定する。
の周囲に閉領域を定め、その閉領域を部分空間に分割す
る。次に、各部分空間に対して１本のレイを割当て、部
分空間と割当てたレイを同一視する。その後、送信アン
テナから放射される有限な本数のレイの伝搬経路を追跡
し、その結果を基にアンテナの周囲全空間の電波伝搬状
況を推定する。
【０００７】アンテナ周囲の閉領域を部分空間に分割す
る方法の一例が、文献「Scott Y. Seidel,et al ：“Si
te-Specific Propagation Prediction for Wireless In
-Building Personal Communication System Design”,
IEEE TRANSACTIONS VEHICULAR TECHNOLOGY，VOL.43，N
O.4 ,NOVEMBER 1994 ，pp.879-891」に説明されてい
る。
る方法の一例が、文献「Scott Y. Seidel,et al ：“Si
te-Specific Propagation Prediction for Wireless In
-Building Personal Communication System Design”,
IEEE TRANSACTIONS VEHICULAR TECHNOLOGY，VOL.43，N
O.4 ,NOVEMBER 1994 ，pp.879-891」に説明されてい
る。
【０００８】この文献では、先ず、図１２のように、送
信アンテナ３０１の周囲に、正２０面体の３次元閉領域
を設ける。次に、正２０面体を構成する面、すなわち図
１３のように、頂点４０６，４０７，４０８からなる正
三角形の平面を取出した後、点４０６，４１０，４１１
を用いて各辺を等間隔に分割する。頂点４０６，４０
７，４０８からなる正三角形の各辺に平行で、かつ分割
点を通る線分を引くことにより、もとの正三角形と相似
な三角形が新しく内部に作られる。以上の処理を図１２
の正２０面体をなす全ての面に対して行い、正２０面体
の重心と新しく生成された各正三角形の頂点とを結ぶ方
向へ、重心からの距離が等しくなるように各正三角形の
頂点を移動させると、例えば図１４のようになる。
信アンテナ３０１の周囲に、正２０面体の３次元閉領域
を設ける。次に、正２０面体を構成する面、すなわち図
１３のように、頂点４０６，４０７，４０８からなる正
三角形の平面を取出した後、点４０６，４１０，４１１
を用いて各辺を等間隔に分割する。頂点４０６，４０
７，４０８からなる正三角形の各辺に平行で、かつ分割
点を通る線分を引くことにより、もとの正三角形と相似
な三角形が新しく内部に作られる。以上の処理を図１２
の正２０面体をなす全ての面に対して行い、正２０面体
の重心と新しく生成された各正三角形の頂点とを結ぶ方
向へ、重心からの距離が等しくなるように各正三角形の
頂点を移動させると、例えば図１４のようになる。
【０００９】図１４は図１２の正２０面体の各面をなす
正三角形の１辺を２等分した場合の図である。もとの正
２０面体の重心に位置する送信アンテナ５０１から放射
されるレイは、送信アンテナ５０１と図１４の多面体の
各頂点とを結ぶ各方向に放射される。図１４には、一例
として頂点５０２を通るレイ５０４を示してある。この
とき、図１３の多角形と送信アンテナ５０１から定まる
部分空間５０３はレイ５０４と同一視される。
正三角形の１辺を２等分した場合の図である。もとの正
２０面体の重心に位置する送信アンテナ５０１から放射
されるレイは、送信アンテナ５０１と図１４の多面体の
各頂点とを結ぶ各方向に放射される。図１４には、一例
として頂点５０２を通るレイ５０４を示してある。この
とき、図１３の多角形と送信アンテナ５０１から定まる
部分空間５０３はレイ５０４と同一視される。
【００１０】図１５はレイ６０１と同一視される部分空
間６０５を示す図である。点６０２において、レイ６０
１に垂直な断面６０３の面積を、今後、点６０２におけ
る空間分解能と定義する。点６０２における空間分解能
は送信アンテナ６０４と点６０２との間の長さ、すなわ
ち当該レイの伝搬距離が長くなるに従い増大する。この
ように広い空間を１本のレイと同一視するようになるた
め、ラウンチング法による電波伝搬推定精度は、伝搬距
離が増大するに従い低下することになる。
間６０５を示す図である。点６０２において、レイ６０
１に垂直な断面６０３の面積を、今後、点６０２におけ
る空間分解能と定義する。点６０２における空間分解能
は送信アンテナ６０４と点６０２との間の長さ、すなわ
ち当該レイの伝搬距離が長くなるに従い増大する。この
ように広い空間を１本のレイと同一視するようになるた
め、ラウンチング法による電波伝搬推定精度は、伝搬距
離が増大するに従い低下することになる。
【００１１】そこで、空間分解能を伝搬距離に関わら
ず、常に一定の値以下に維持することを目的として、所
定の距離を伝搬するごとにレイを分割する手法が提案さ
れた。この手法の一例が、文献「Steven Fortune：“EF
FICIENT ALGORITHMS FOR PREDICTION OF INDOOR RADIO
PROPAGATION ”,in Proceedings of the 48th IEEE Veh
icular Technology Conference,May,1998,pp.572-576」
に記されている。図１６はこの文献の手法を説明する図
である。
ず、常に一定の値以下に維持することを目的として、所
定の距離を伝搬するごとにレイを分割する手法が提案さ
れた。この手法の一例が、文献「Steven Fortune：“EF
FICIENT ALGORITHMS FOR PREDICTION OF INDOOR RADIO
PROPAGATION ”,in Proceedings of the 48th IEEE Veh
icular Technology Conference,May,1998,pp.572-576」
に記されている。図１６はこの文献の手法を説明する図
である。
【００１２】図１６の三角錐７１２は、割当てられたレ
イ７０２の送信アンテナ７０１からの伝搬距離が所定の
値に達すると、三角錐７０８から三角錐７１１に細分さ
れる。細分された三角錐には、それぞれレイ７０４から
レイ７０７が割当てられ、以降、同様な処理が繰返され
る。この結果、空間分解能は伝搬距離に関わらず常に一
定値以下に保たれることになるのである。
イ７０２の送信アンテナ７０１からの伝搬距離が所定の
値に達すると、三角錐７０８から三角錐７１１に細分さ
れる。細分された三角錐には、それぞれレイ７０４から
レイ７０７が割当てられ、以降、同様な処理が繰返され
る。この結果、空間分解能は伝搬距離に関わらず常に一
定値以下に保たれることになるのである。
【００１３】図１７は観測領域０１７と、送信点０１５
と、受信点０１６、２つの物体００１，００２とが与え
られた場合のラウンチング法の動作について説明する図
である。図１７では、簡単のため２次元平面に限定して
動作の説明をするが、実際には３次元空間内で動作させ
ることは勿論である。
と、受信点０１６、２つの物体００１，００２とが与え
られた場合のラウンチング法の動作について説明する図
である。図１７では、簡単のため２次元平面に限定して
動作の説明をするが、実際には３次元空間内で動作させ
ることは勿論である。
【００１４】離散的な間隔で伝搬経路００３の方向へ放
射されたレイは、点０１８で物体００１へ入射し、その
結果、反射レイ００５と透過レイ００４とが生成され
る。反射レイ００５は、さらに点０１９で物体００２へ
入射し、再び反射レイ００６と透過レイ００７とが生成
される。生成された反射レイ００６は受信点０１６の近
傍を通過するため、レイ００６を受信点に到来したレイ
とみなし、経路００３，００５，００６からなる経路を
送信点０１５と受信点０１６とを結ぶ伝搬経路の一つと
みなす。
射されたレイは、点０１８で物体００１へ入射し、その
結果、反射レイ００５と透過レイ００４とが生成され
る。反射レイ００５は、さらに点０１９で物体００２へ
入射し、再び反射レイ００６と透過レイ００７とが生成
される。生成された反射レイ００６は受信点０１６の近
傍を通過するため、レイ００６を受信点に到来したレイ
とみなし、経路００３，００５，００６からなる経路を
送信点０１５と受信点０１６とを結ぶ伝搬経路の一つと
みなす。
【００１５】伝搬経路００３，００５，００６から求ま
る受信電力ならびに到来遅延時間は図１８のように記録
される。図１８の横軸１０３はレイが送信点０１５から
受信点０１６まで経路００３，００５，００６を経由し
て到達するのに要した遅延時間を示し、縦軸１０２は当
該経路を経て到達したレイが有する電力強度を表す。透
過レイ００４ならびに００７についても、伝搬経路００
３，００５，００６と同様の反射・透過レイ探索を繰り
返し、受信点０１６の近傍をレイが通過した場合には、
伝搬経路００６の場合と同様に到来レイとして扱い、以
上の処理を探索終了条件を満たすまで続ける。
る受信電力ならびに到来遅延時間は図１８のように記録
される。図１８の横軸１０３はレイが送信点０１５から
受信点０１６まで経路００３，００５，００６を経由し
て到達するのに要した遅延時間を示し、縦軸１０２は当
該経路を経て到達したレイが有する電力強度を表す。透
過レイ００４ならびに００７についても、伝搬経路００
３，００５，００６と同様の反射・透過レイ探索を繰り
返し、受信点０１６の近傍をレイが通過した場合には、
伝搬経路００６の場合と同様に到来レイとして扱い、以
上の処理を探索終了条件を満たすまで続ける。
【００１６】探索終了条件は、反射・透過点での受信電
界強度が所定値を下回った場合、または反射・透過の総
数が所定の回数に達した場合などとする。送信点０１５
から伝搬経路００３の方向へ放射されたレイの反射・透
過経路探索が終了した後、方向００８から０１４のよう
に他の放射方向へ向けてレイを放射した場合についても
同様に反射・透過経路探索を行う。事前に定めた全ての
方向へレイを放射し、以上の伝搬経路探索を行うことに
より、受信点０１６における遅延プロファイルである図
１９が得られる。図１９の横軸２０３は送信点０１５か
ら受信点０１６へ到達するのに要した時間を示し、縦軸
２０２は当該経路を通過したレイが有する電力強度を表
す。
界強度が所定値を下回った場合、または反射・透過の総
数が所定の回数に達した場合などとする。送信点０１５
から伝搬経路００３の方向へ放射されたレイの反射・透
過経路探索が終了した後、方向００８から０１４のよう
に他の放射方向へ向けてレイを放射した場合についても
同様に反射・透過経路探索を行う。事前に定めた全ての
方向へレイを放射し、以上の伝搬経路探索を行うことに
より、受信点０１６における遅延プロファイルである図
１９が得られる。図１９の横軸２０３は送信点０１５か
ら受信点０１６へ到達するのに要した時間を示し、縦軸
２０２は当該経路を通過したレイが有する電力強度を表
す。
【００１７】受信電力は図１９に示した全ての経路の電
力強度の和で与えられ、歪の程度を示す遅延拡がりは、
各遅延時間において受信電力に対する当該電力強度の電
力比を当該遅延時間の出現確率として、遅延時間の標準
偏差により与えられる。
力強度の和で与えられ、歪の程度を示す遅延拡がりは、
各遅延時間において受信電力に対する当該電力強度の電
力比を当該遅延時間の出現確率として、遅延時間の標準
偏差により与えられる。
【００１８】
【発明が解決しようとする課題】上述したラウンチング
法では、レイと近傍の部分空間が同一視されるために、
レイが障壁と反射する場合は、レイ近傍の部分空間でも
同様な反射が生じているとみなされる。そのため、推定
精度の観点からは、反射点における空間分解能を当該障
壁の面積に比べて大きくなり過ぎないように設定するこ
とが望ましい。しかし、推定精度の向上を目的として全
体の空間分解能の上限値を各障壁との反射の際に要求さ
れる空間分解能の最小値に設定すると、電波伝搬推定に
要する計算時間が増大するという問題点がある。その理
由は、空間分解能を小さくすると放射するレイの総数が
増大し、経路探索に要する全体の時間が増大するからで
ある。
法では、レイと近傍の部分空間が同一視されるために、
レイが障壁と反射する場合は、レイ近傍の部分空間でも
同様な反射が生じているとみなされる。そのため、推定
精度の観点からは、反射点における空間分解能を当該障
壁の面積に比べて大きくなり過ぎないように設定するこ
とが望ましい。しかし、推定精度の向上を目的として全
体の空間分解能の上限値を各障壁との反射の際に要求さ
れる空間分解能の最小値に設定すると、電波伝搬推定に
要する計算時間が増大するという問題点がある。その理
由は、空間分解能を小さくすると放射するレイの総数が
増大し、経路探索に要する全体の時間が増大するからで
ある。
【００１９】一方、推定に要する計算時間を抑えるため
に空間分解能の上限値を比較的大きく設定すると、全体
の推定精度に劣化が生じるという問題点がある。その理
由は、所要空間分解能が小さな障壁でレイが反射する場
合には、実際には反射が生じていないレイ近傍の領域で
もレイと同様に反射が生じているとみなすからである。
に空間分解能の上限値を比較的大きく設定すると、全体
の推定精度に劣化が生じるという問題点がある。その理
由は、所要空間分解能が小さな障壁でレイが反射する場
合には、実際には反射が生じていないレイ近傍の領域で
もレイと同様に反射が生じているとみなすからである。
【００２０】本発明は、ラウンチング法において伝搬推
定の対象環境が与えられた場合に、推定精度の劣化を抑
えながら所要推定時間を短縮化する手法を提供すること
を目的とする。
定の対象環境が与えられた場合に、推定精度の劣化を抑
えながら所要推定時間を短縮化する手法を提供すること
を目的とする。
【００２１】
【課題を解決するための手段】本発明によれば、３次元
空間上に規定される観測領域内に複数の物体が与えら
れ、前記観測領域内に電波線（レイ）が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物の検出を行って電
波伝搬特性の推定をなす電波伝搬特性推定方法であっ
て、前記レイが前記遮蔽物と衝突する場面において、前
記レイを包含する閉空間を複数の部分空間に分割し、こ
れ等各部分空間により包含されるレイを新たに生成する
ことを特徴とする電波伝搬特性推定方法が得られる。
空間上に規定される観測領域内に複数の物体が与えら
れ、前記観測領域内に電波線（レイ）が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物の検出を行って電
波伝搬特性の推定をなす電波伝搬特性推定方法であっ
て、前記レイが前記遮蔽物と衝突する場面において、前
記レイを包含する閉空間を複数の部分空間に分割し、こ
れ等各部分空間により包含されるレイを新たに生成する
ことを特徴とする電波伝搬特性推定方法が得られる。
【００２２】また、本発明によれば、３次元空間上に規
定される観測領域内に複数の物体が与えられ、前記観測
領域内に電波線（レイ）が与えられると共に、前記レイ
の近傍に前記レイと同一視される空間が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突
した際に、予め定められた基準に基づき前記レイの空間
分解能の調節を行う電波伝搬特性推定におけるレイの空
間分解能制御方法であって、前記遮蔽物と前記レイとの
衝突点を検出する第一のステップと、前記レイと同一視
される前記空間を分割するか否かの判定を、所定判定基
準に従って行う第二のステップと、この第二のステップ
で前記レイと同一視される前記空間を分割すると判定さ
れた場合に、前記レイと同一視される前記空間を複数の
部分空間に分割する第三のステップと、この第三のステ
ップにおいて分割された前記部分空間の各々に対して同
一視するレイを割当てる第四のステップとを含むことを
特徴とするレイの空間分解能制御方法が得られる。
定される観測領域内に複数の物体が与えられ、前記観測
領域内に電波線（レイ）が与えられると共に、前記レイ
の近傍に前記レイと同一視される空間が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突
した際に、予め定められた基準に基づき前記レイの空間
分解能の調節を行う電波伝搬特性推定におけるレイの空
間分解能制御方法であって、前記遮蔽物と前記レイとの
衝突点を検出する第一のステップと、前記レイと同一視
される前記空間を分割するか否かの判定を、所定判定基
準に従って行う第二のステップと、この第二のステップ
で前記レイと同一視される前記空間を分割すると判定さ
れた場合に、前記レイと同一視される前記空間を複数の
部分空間に分割する第三のステップと、この第三のステ
ップにおいて分割された前記部分空間の各々に対して同
一視するレイを割当てる第四のステップとを含むことを
特徴とするレイの空間分解能制御方法が得られる。
【００２３】そして、前記第三のステップにおける分割
数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番
目（ｉは１〜Ｍ）の前記部分空間に対して前記第四のス
テップにおいてレイを割当てた後、割当てた各レイに対
して再び前記第一から第四のステップを適用し、前記第
二のステップにおいて当該レイと同一視される前記空間
を分割しないと判定されるまで前記第一から第四のステ
ップを、順次繰返すことを特徴とする。
数はＭ個（Ｍは２以上の整数）に設定されており、ｉ番
目（ｉは１〜Ｍ）の前記部分空間に対して前記第四のス
テップにおいてレイを割当てた後、割当てた各レイに対
して再び前記第一から第四のステップを適用し、前記第
二のステップにおいて当該レイと同一視される前記空間
を分割しないと判定されるまで前記第一から第四のステ
ップを、順次繰返すことを特徴とする。
【００２４】また、前記第二のステップにおける前記判
定基準は、前記遮蔽物との衝突点における前記レイの電
界強度が所定の値以上であるか否かを調べる第一の判定
基準と、前記遮蔽物の衝突面を含む平面による前記レイ
近傍の前記空間の断面に対する前記遮蔽物の衝突面の面
積比が所定値以下であるか否かを調べる第二の判定基準
とからなることを特徴とし、更には、前記判定基準は、
前記断面の面積が所定の値以上であり、かつ前記レイと
前記衝突面の法線との成す角（入射角）が所定値以上で
あるか否かを調べる第三の判定基準を含むことを特徴と
する。
定基準は、前記遮蔽物との衝突点における前記レイの電
界強度が所定の値以上であるか否かを調べる第一の判定
基準と、前記遮蔽物の衝突面を含む平面による前記レイ
近傍の前記空間の断面に対する前記遮蔽物の衝突面の面
積比が所定値以下であるか否かを調べる第二の判定基準
とからなることを特徴とし、更には、前記判定基準は、
前記断面の面積が所定の値以上であり、かつ前記レイと
前記衝突面の法線との成す角（入射角）が所定値以上で
あるか否かを調べる第三の判定基準を含むことを特徴と
する。
【００２５】また、円錐に内接する錐体のうち、前記円
錐との体積誤差が小さい錐体を設け、その中心線を前記
レイの経路とすることにより前記レイ近傍の前記空間を
前記錐体によって形成し、前記第三のステップにおい
て、前記レイ近傍の空間を形成する前記錐体に相似な複
数の錐体を用いて前記レイ近傍の空間を複数の部分空間
に分割する際に、前記レイ近傍の空間を形成する前記錐
体と、この錐体を分割する前記錐体に相似な前記複数の
錐体の和との体積誤差が小さくなるように、前記レイの
近傍に前記複数の錐体を配置することを特徴とする。
錐との体積誤差が小さい錐体を設け、その中心線を前記
レイの経路とすることにより前記レイ近傍の前記空間を
前記錐体によって形成し、前記第三のステップにおい
て、前記レイ近傍の空間を形成する前記錐体に相似な複
数の錐体を用いて前記レイ近傍の空間を複数の部分空間
に分割する際に、前記レイ近傍の空間を形成する前記錐
体と、この錐体を分割する前記錐体に相似な前記複数の
錐体の和との体積誤差が小さくなるように、前記レイの
近傍に前記複数の錐体を配置することを特徴とする。
【００２６】更に、前記第三のステップにおいて、円錐
に内接する前記錐体として三角錐を用い、前記三角錐に
相似な複数の同一三角錐を前記三角錐内部に隙間無く配
置できるように前記相似三角錐の大きさを設定すること
を特徴とし、また前記第三のステップにおいて、円錐に
内接する前記錐体として正六角錐を用い、前記正六角錐
に相似な正六角錐を前記正六角錐の中心に配置し、その
前記相似正六角錐の各側面に前記相似正六角錐と同一な
６個の正六角錐を各々側面が接するように配置し、全７
個の前記相似正六角錐の和の体積と、もとの前記正六角
錐の体積誤差が最小となるように前記相似六角錐の大き
さを設定し配置することを特徴とする。
に内接する前記錐体として三角錐を用い、前記三角錐に
相似な複数の同一三角錐を前記三角錐内部に隙間無く配
置できるように前記相似三角錐の大きさを設定すること
を特徴とし、また前記第三のステップにおいて、円錐に
内接する前記錐体として正六角錐を用い、前記正六角錐
に相似な正六角錐を前記正六角錐の中心に配置し、その
前記相似正六角錐の各側面に前記相似正六角錐と同一な
６個の正六角錐を各々側面が接するように配置し、全７
個の前記相似正六角錐の和の体積と、もとの前記正六角
錐の体積誤差が最小となるように前記相似六角錐の大き
さを設定し配置することを特徴とする。
【００２７】更にはまた、前記第三のステップにおい
て、円錐に内接する前記錐体として正六角錐を用い、前
記正六角錐に相似な正六角錐を前記正六角錐の中心に配
置し、その前記相似正六角錐の各側面に前記相似正六角
錐と同一な６個の正六角錐を各々側面が接するように配
置し、全７個の前記相似正六角錐の和からなる錐体の最
外郭側面がもとの前記正六角錐の側面に内接するように
前記相似六角錐の大きさを設定し配置することを特徴と
する。また、前記第三のステップにおいて、円錐に内接
する前記錐体として正六角錐を用い、前記正六角錐に相
似な正六角錐を前記正六角錐の中心に配置し、その前記
相似正六角錐の各側面に前記相似正六角錐と同一な６個
の正六角錐を各々側面が接するように配置し、全７個の
前記相似正六角錐の和からなる錐体の最外郭側面に、も
との前記正六角錐が内接するように前記相似六角錐の大
きさを設定し配置することを特徴とする。
て、円錐に内接する前記錐体として正六角錐を用い、前
記正六角錐に相似な正六角錐を前記正六角錐の中心に配
置し、その前記相似正六角錐の各側面に前記相似正六角
錐と同一な６個の正六角錐を各々側面が接するように配
置し、全７個の前記相似正六角錐の和からなる錐体の最
外郭側面がもとの前記正六角錐の側面に内接するように
前記相似六角錐の大きさを設定し配置することを特徴と
する。また、前記第三のステップにおいて、円錐に内接
する前記錐体として正六角錐を用い、前記正六角錐に相
似な正六角錐を前記正六角錐の中心に配置し、その前記
相似正六角錐の各側面に前記相似正六角錐と同一な６個
の正六角錐を各々側面が接するように配置し、全７個の
前記相似正六角錐の和からなる錐体の最外郭側面に、も
との前記正六角錐が内接するように前記相似六角錐の大
きさを設定し配置することを特徴とする。
【００２８】本発明によれば、３次元空間上に規定され
る観測領域内に複数の物体が与えられ、前記観測領域内
に電波線（レイ）が与えられ、前記物体のうち前記レイ
を遮蔽する遮蔽物の検出を行って電波伝搬特性の推定を
なす電波伝搬特性推定システムであって、前記レイが前
記遮蔽物と衝突する場面において、前記レイを包含する
閉空間を複数の部分空間に分割し、これ等各部分空間に
より包含されるレイを新たに生成することを特徴とする
電波伝搬特性推定システムが得られる。
る観測領域内に複数の物体が与えられ、前記観測領域内
に電波線（レイ）が与えられ、前記物体のうち前記レイ
を遮蔽する遮蔽物の検出を行って電波伝搬特性の推定を
なす電波伝搬特性推定システムであって、前記レイが前
記遮蔽物と衝突する場面において、前記レイを包含する
閉空間を複数の部分空間に分割し、これ等各部分空間に
より包含されるレイを新たに生成することを特徴とする
電波伝搬特性推定システムが得られる。
【００２９】また、本発明によれば、３次元空間上に規
定される観測領域内に複数の物体が与えられ、前記観測
領域内に電波線（レイ）が与えられると共に、前記レイ
の近傍に前記レイと同一視される空間が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突
した際に、予め定められた基準に基づき前記レイの空間
分解能の調節を行う電波伝搬特性推定におけるレイの空
間分解能制御システムであって、前記遮蔽物と前記レイ
との衝突点を検出する衝突点検出手段と、前記レイと同
一視される前記空間を分割するか否かの判定を、所定判
定基準に従って行う判定手段と、この判定手段で前記レ
イと同一視される前記空間を分割すると判定された場合
に、前記レイと同一視される前記空間を複数の部分空間
に分割する分割手段と、この分割手段において分割され
た前記部分空間の各々に対して同一視するレイを割当て
る割当手段とを含むことを特徴とするレイの空間分解能
制御システムが得られる。
定される観測領域内に複数の物体が与えられ、前記観測
領域内に電波線（レイ）が与えられると共に、前記レイ
の近傍に前記レイと同一視される空間が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突
した際に、予め定められた基準に基づき前記レイの空間
分解能の調節を行う電波伝搬特性推定におけるレイの空
間分解能制御システムであって、前記遮蔽物と前記レイ
との衝突点を検出する衝突点検出手段と、前記レイと同
一視される前記空間を分割するか否かの判定を、所定判
定基準に従って行う判定手段と、この判定手段で前記レ
イと同一視される前記空間を分割すると判定された場合
に、前記レイと同一視される前記空間を複数の部分空間
に分割する分割手段と、この分割手段において分割され
た前記部分空間の各々に対して同一視するレイを割当て
る割当手段とを含むことを特徴とするレイの空間分解能
制御システムが得られる。
【００３０】そして、前記分割手段における分割数はＭ
個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉ
は１〜Ｍ）の前記部分空間に対して前記割当手段におい
てレイを割当てた後、割当てた各レイに対して再び前記
衝突点検出手段、前記判定手段、前記分割手段及び前記
割当手段を実行制御し、前記判定手段において当該レイ
と同一視される前記空間を分割しないと判定されるま
で、前記衝突点検出手段、前記判定手段、前記分割手段
及び前記割当手段を、順次繰返し実行制御する制御手段
とを含むことを特徴とする。
個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉ
は１〜Ｍ）の前記部分空間に対して前記割当手段におい
てレイを割当てた後、割当てた各レイに対して再び前記
衝突点検出手段、前記判定手段、前記分割手段及び前記
割当手段を実行制御し、前記判定手段において当該レイ
と同一視される前記空間を分割しないと判定されるま
で、前記衝突点検出手段、前記判定手段、前記分割手段
及び前記割当手段を、順次繰返し実行制御する制御手段
とを含むことを特徴とする。
【００３１】また、前記判定手段における前記判定基準
は、前記遮蔽物との衝突点における前記レイの電界強度
が所定の値以上であるか否かを調べる第一の判定基準
と、前記遮蔽物の衝突面を含む平面による前記レイ近傍
の前記空間の断面に対する前記遮蔽物の衝突面の面積比
が所定値以下であるか否かを調べる第二の判定基準とか
らなることを特徴とし、更に、前記判定基準は、前記断
面の面積が所定の値以上であり、かつ前記レイと前記衝
突面の法線との成す角（入射角）が所定値以上であるか
否かを調べる第三の判定基準を含むことを特徴とする。
は、前記遮蔽物との衝突点における前記レイの電界強度
が所定の値以上であるか否かを調べる第一の判定基準
と、前記遮蔽物の衝突面を含む平面による前記レイ近傍
の前記空間の断面に対する前記遮蔽物の衝突面の面積比
が所定値以下であるか否かを調べる第二の判定基準とか
らなることを特徴とし、更に、前記判定基準は、前記断
面の面積が所定の値以上であり、かつ前記レイと前記衝
突面の法線との成す角（入射角）が所定値以上であるか
否かを調べる第三の判定基準を含むことを特徴とする。
【００３２】また、円錐に内接する錐体のうち、前記円
錐との体積誤差が小さい錐体を設け、その中心線を前記
レイの経路とすることにより前記レイ近傍の前記空間を
前記錐体によって形成し、前記分割手段において、前記
レイ近傍の空間を形成する前記錐体に相似な複数の錐体
を用いて前記レイ近傍の空間を複数の部分空間に分割す
る際に、前記レイ近傍の空間を形成する前記錐体と、こ
の錐体を分割する前記錐体に相似な前記複数の錐体の和
との体積誤差が小さくなるように、前記レイの近傍に前
記複数の錐体を配置することを特徴とする。
錐との体積誤差が小さい錐体を設け、その中心線を前記
レイの経路とすることにより前記レイ近傍の前記空間を
前記錐体によって形成し、前記分割手段において、前記
レイ近傍の空間を形成する前記錐体に相似な複数の錐体
を用いて前記レイ近傍の空間を複数の部分空間に分割す
る際に、前記レイ近傍の空間を形成する前記錐体と、こ
の錐体を分割する前記錐体に相似な前記複数の錐体の和
との体積誤差が小さくなるように、前記レイの近傍に前
記複数の錐体を配置することを特徴とする。
【００３３】更に、前記分割手段において、円錐に内接
する前記錐体として三角錐を用い、前記三角錐に相似な
複数の同一三角錐を前記三角錐内部に隙間無く配置でき
るように前記相似三角錐の大きさを設定することを特徴
とし、また、前記分割手段において、円錐に内接する前
記錐体として正六角錐を用い、前記正六角錐に相似な正
六角錐を前記正六角錐の中心に配置し、その前記相似正
六角錐の各側面に前記相似正六角錐と同一な６個の正六
角錐を各々側面が接するように配置し、全７個の前記相
似正六角錐の和の体積と、もとの前記正六角錐の体積誤
差が最小となるように前記相似六角錐の大きさを設定し
配置することを特徴とする。
する前記錐体として三角錐を用い、前記三角錐に相似な
複数の同一三角錐を前記三角錐内部に隙間無く配置でき
るように前記相似三角錐の大きさを設定することを特徴
とし、また、前記分割手段において、円錐に内接する前
記錐体として正六角錐を用い、前記正六角錐に相似な正
六角錐を前記正六角錐の中心に配置し、その前記相似正
六角錐の各側面に前記相似正六角錐と同一な６個の正六
角錐を各々側面が接するように配置し、全７個の前記相
似正六角錐の和の体積と、もとの前記正六角錐の体積誤
差が最小となるように前記相似六角錐の大きさを設定し
配置することを特徴とする。
【００３４】また、前記分割手段において、円錐に内接
する前記錐体として正六角錐を用い、前記正六角錐に相
似な正六角錐を前記正六角錐の中心に配置し、その前記
相似正六角錐の各側面に前記相似正六角錐と同一な６個
の正六角錐を各々側面が接するように配置し、全７個の
前記相似正六角錐の和からなる錐体の最外郭側面がもと
の前記正六角錐の側面に内接するように前記相似六角錐
の大きさを設定し配置することを特徴とする。更に、前
記分割手段において、円錐に内接する前記錐体として正
六角錐を用い、前記正六角錐に相似な正六角錐を前記正
六角錐の中心に配置し、その前記相似正六角錐の各側面
に前記相似正六角錐と同一な６個の正六角錐を各々側面
が接するように配置し、全７個の前記相似正六角錐の和
からなる錐体の最外郭側面に、もとの前記正六角錐が内
接するように前記相似六角錐の大きさを設定し配置する
ことを特徴とする。
する前記錐体として正六角錐を用い、前記正六角錐に相
似な正六角錐を前記正六角錐の中心に配置し、その前記
相似正六角錐の各側面に前記相似正六角錐と同一な６個
の正六角錐を各々側面が接するように配置し、全７個の
前記相似正六角錐の和からなる錐体の最外郭側面がもと
の前記正六角錐の側面に内接するように前記相似六角錐
の大きさを設定し配置することを特徴とする。更に、前
記分割手段において、円錐に内接する前記錐体として正
六角錐を用い、前記正六角錐に相似な正六角錐を前記正
六角錐の中心に配置し、その前記相似正六角錐の各側面
に前記相似正六角錐と同一な６個の正六角錐を各々側面
が接するように配置し、全７個の前記相似正六角錐の和
からなる錐体の最外郭側面に、もとの前記正六角錐が内
接するように前記相似六角錐の大きさを設定し配置する
ことを特徴とする。
【００３５】本発明によれば、３次元空間上に規定され
る観測領域内に複数の物体が与えられ、前記観測領域内
に電波線（レイ）が与えられ、前記物体のうち前記レイ
を遮蔽する遮蔽物の検出を行って電波伝搬特性の推定を
なす電波伝搬特性推定方法をコンピュータに実行させる
ためのプログラムであって、前記レイが前記遮蔽物と衝
突する場面において、前記レイを包含する閉空間を複数
の部分空間に分割し、これ等各部分空間により包含され
るレイを新たに生成することを特徴とするプログラムが
得られる。
る観測領域内に複数の物体が与えられ、前記観測領域内
に電波線（レイ）が与えられ、前記物体のうち前記レイ
を遮蔽する遮蔽物の検出を行って電波伝搬特性の推定を
なす電波伝搬特性推定方法をコンピュータに実行させる
ためのプログラムであって、前記レイが前記遮蔽物と衝
突する場面において、前記レイを包含する閉空間を複数
の部分空間に分割し、これ等各部分空間により包含され
るレイを新たに生成することを特徴とするプログラムが
得られる。
【００３６】また、本発明によれば、３次元空間上に規
定される観測領域内に複数の物体が与えられ、前記観測
領域内に電波線（レイ）が与えられると共に、前記レイ
の近傍に前記レイと同一視される空間が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突
した際に、予め定められた基準に基づき前記レイの空間
分解能の調節を行う電波伝搬特性推定におけるレイの空
間分解能制御方法をコンピュータに実行させるためのプ
ログラムであって、前記遮蔽物と前記レイとの衝突点を
検出する第一のステップと、前記レイと同一視される前
記空間を分割するか否かの判定を、所定判定基準に従っ
て行う第二のステップと、この第二のステップで前記レ
イと同一視される前記空間を分割すると判定された場合
に、前記レイと同一視される前記空間を複数の部分空間
に分割する第三のステップと、この第三のステップにお
いて分割された前記部分空間の各々に対して同一視する
レイを割当てる第四のステップと含むことを特徴とする
プログラムが得られる。
定される観測領域内に複数の物体が与えられ、前記観測
領域内に電波線（レイ）が与えられると共に、前記レイ
の近傍に前記レイと同一視される空間が与えられ、前記
物体のうち前記レイを遮蔽する遮蔽物と前記レイが衝突
した際に、予め定められた基準に基づき前記レイの空間
分解能の調節を行う電波伝搬特性推定におけるレイの空
間分解能制御方法をコンピュータに実行させるためのプ
ログラムであって、前記遮蔽物と前記レイとの衝突点を
検出する第一のステップと、前記レイと同一視される前
記空間を分割するか否かの判定を、所定判定基準に従っ
て行う第二のステップと、この第二のステップで前記レ
イと同一視される前記空間を分割すると判定された場合
に、前記レイと同一視される前記空間を複数の部分空間
に分割する第三のステップと、この第三のステップにお
いて分割された前記部分空間の各々に対して同一視する
レイを割当てる第四のステップと含むことを特徴とする
プログラムが得られる。
【００３７】本発明の作用を述べる。本発明では、与え
られた観測領域内をレイが伝搬する場合、その経路上に
存在する障壁とレイが衝突する際に、所定の条件が満た
されるまで当該レイ近傍の部分空間を複数の部分空間に
分割する。従来の方法では、反射とは無関係に伝搬距離
に応じて一律に部分空間を分割していたが、本発明の方
法によれば、一律に扱った場合に推定精度の劣化を招く
と考えられる部分空間は、劣化の影響が少なくなるまで
分割し、細分化された各部分空間を個別に取り扱うこと
が可能となる。このような手法を採用することにより、
電波伝搬の推定精度を従来方法よりも向上させることが
可能となる。一方、推定精度の劣化の影響が小さいと考
えられる場合は、当該部分空間を分割しないため、空間
分解能の初期値を比較的大きめに設定しておけば、従来
の方法において推定精度の向上を狙って一律に空間分解
能を小さく設定した場合よりもレイの総数が低減し、伝
搬推定に要する全体の計算時間を短縮させることができ
る。
られた観測領域内をレイが伝搬する場合、その経路上に
存在する障壁とレイが衝突する際に、所定の条件が満た
されるまで当該レイ近傍の部分空間を複数の部分空間に
分割する。従来の方法では、反射とは無関係に伝搬距離
に応じて一律に部分空間を分割していたが、本発明の方
法によれば、一律に扱った場合に推定精度の劣化を招く
と考えられる部分空間は、劣化の影響が少なくなるまで
分割し、細分化された各部分空間を個別に取り扱うこと
が可能となる。このような手法を採用することにより、
電波伝搬の推定精度を従来方法よりも向上させることが
可能となる。一方、推定精度の劣化の影響が小さいと考
えられる場合は、当該部分空間を分割しないため、空間
分解能の初期値を比較的大きめに設定しておけば、従来
の方法において推定精度の向上を狙って一律に空間分解
能を小さく設定した場合よりもレイの総数が低減し、伝
搬推定に要する全体の計算時間を短縮させることができ
る。
【００３８】これらの特徴により、本発明では、与えら
れた推定領域に応じて空間分解能を適応制御することが
可能なことから、従来の方法に比べて推定精度の向上と
計算時間の短縮を両立させることが可能である。
れた推定領域に応じて空間分解能を適応制御することが
可能なことから、従来の方法に比べて推定精度の向上と
計算時間の短縮を両立させることが可能である。
【００３９】
【発明の実施の形態】以下に、図面を参照しながら本発
明の実施例を詳細に説明する。図１は本発明の実施例の
動作を示す流れ図であり、２つのルーチンＡ，Ｂ（８１
８，８１９）よりなるものとする。図１では、３次元空
間上に観測領域が与えられ、当該観測領域内に壁や什器
等の物体が与えられ、当該観測領域内の点Ｐから当該伝
搬方向へ放射されるレイが与えられることを前提とし
（ステップ８０２）、反射点における反射・透過レイを
生成するまでの一連の流れを示す。
明の実施例を詳細に説明する。図１は本発明の実施例の
動作を示す流れ図であり、２つのルーチンＡ，Ｂ（８１
８，８１９）よりなるものとする。図１では、３次元空
間上に観測領域が与えられ、当該観測領域内に壁や什器
等の物体が与えられ、当該観測領域内の点Ｐから当該伝
搬方向へ放射されるレイが与えられることを前提とし
（ステップ８０２）、反射点における反射・透過レイを
生成するまでの一連の流れを示す。
【００４０】ここで、点Ｐとは、ラウンチング法におけ
る反射・透過点もしくは送信アンテナの位置を指し、具
体的には従来技術の説明において参照した図１７におけ
る点０１８、点０１９、送信アンテナ０１５等を指すも
のとする。また、当該伝搬方向とは、点Ｐが送信アンテ
ナである場合は、送信アンテナの周囲に設けた離散的な
放射方向を指し、具体的には従来の技術で説明したよう
に図１７の方向００３、００８〜０１４等、もしくは図
１４の送信アンテナ５０１と多面体の各頂点を結ぶ方
向、例えばレイ５０４の方向等を指すものとする。一
方、点Ｐが反射・透過点である場合は、反射レイもしく
は透過レイの方向を指し、具体的には図１７における伝
搬経路００３〜００７等を指すものとする。
る反射・透過点もしくは送信アンテナの位置を指し、具
体的には従来技術の説明において参照した図１７におけ
る点０１８、点０１９、送信アンテナ０１５等を指すも
のとする。また、当該伝搬方向とは、点Ｐが送信アンテ
ナである場合は、送信アンテナの周囲に設けた離散的な
放射方向を指し、具体的には従来の技術で説明したよう
に図１７の方向００３、００８〜０１４等、もしくは図
１４の送信アンテナ５０１と多面体の各頂点を結ぶ方
向、例えばレイ５０４の方向等を指すものとする。一
方、点Ｐが反射・透過点である場合は、反射レイもしく
は透過レイの方向を指し、具体的には図１７における伝
搬経路００３〜００７等を指すものとする。
【００４１】まず、点Ｐから当該伝搬方向へ放射したレ
イが観測領域内の物体と衝突する点、すなわち反射点を
探索し検出する（ステップ８０３，８０４）。次に、反
射点に入射する当該レイと同一視されている部分空間
を、より小さな部分空間に細分化するか否かの判定を事
前に定めた判定基準（後述）に基づいて行う（ステップ
８０５）。
イが観測領域内の物体と衝突する点、すなわち反射点を
探索し検出する（ステップ８０３，８０４）。次に、反
射点に入射する当該レイと同一視されている部分空間
を、より小さな部分空間に細分化するか否かの判定を事
前に定めた判定基準（後述）に基づいて行う（ステップ
８０５）。
【００４２】当該部分空間をより小さな部分空間に細分
する場合は、事前に定めた方法（後述）に基づいて当該
部分空間をＭ個（Ｍは整数）の部分空間に分割し（ステ
ップ８０６）、ルーチンＢ（８１９）へ移行する（ステ
ップ８０７）。当該部分空間の細分化を行わない場合に
は、当該反射点における反射レイと透過レイを生成し終
了する（ステップ８０８，８０９）。また、ステップ８
０４で反射点が検出されなければ終了となる。
する場合は、事前に定めた方法（後述）に基づいて当該
部分空間をＭ個（Ｍは整数）の部分空間に分割し（ステ
ップ８０６）、ルーチンＢ（８１９）へ移行する（ステ
ップ８０７）。当該部分空間の細分化を行わない場合に
は、当該反射点における反射レイと透過レイを生成し終
了する（ステップ８０８，８０９）。また、ステップ８
０４で反射点が検出されなければ終了となる。
【００４３】ルーチンＢ（８１９）では、まず内部変数
ｉ（ｉ＝１〜Ｍ）を定義し、初期値としてｉ＝１を与え
る（ステップ８１１）。次に、ルーチンＡ（８１８）の
ステップ８０６でＭ個に分割された部分空間のうちｉ番
目の部分空間を選択する（ステップ８１２）。そして、
選択したｉ番目の部分空間にレイを割当て（ステップ８
２０）、そのレイの伝搬方向を事前に定めた方法（後
述）に基づいて算出する（ステップ８１３）。
ｉ（ｉ＝１〜Ｍ）を定義し、初期値としてｉ＝１を与え
る（ステップ８１１）。次に、ルーチンＡ（８１８）の
ステップ８０６でＭ個に分割された部分空間のうちｉ番
目の部分空間を選択する（ステップ８１２）。そして、
選択したｉ番目の部分空間にレイを割当て（ステップ８
２０）、そのレイの伝搬方向を事前に定めた方法（後
述）に基づいて算出する（ステップ８１３）。
【００４４】ここで、ステップ８０６で細分化された部
分空間にステップ８２０で割当てられたレイが、ステッ
プ８１３で求めた当該伝搬方向へ起点Ｐから放射される
ことを前提とし、再度ルーチンＡ（８１８）へ移行す
る。ルーチンＡ（８１８）では、細分化された部分空間
をさらに細分化するかどうかの判定を行い、細分化する
場合は、ステップ８０５で再度細分化処理を行い、ルー
チンＢ（８１９）へ移行する。細分化を行わない場合に
は、反射・透過レイを生成して（ステップ８０８）終了
する（ステップ８０９）。
分空間にステップ８２０で割当てられたレイが、ステッ
プ８１３で求めた当該伝搬方向へ起点Ｐから放射される
ことを前提とし、再度ルーチンＡ（８１８）へ移行す
る。ルーチンＡ（８１８）では、細分化された部分空間
をさらに細分化するかどうかの判定を行い、細分化する
場合は、ステップ８０５で再度細分化処理を行い、ルー
チンＢ（８１９）へ移行する。細分化を行わない場合に
は、反射・透過レイを生成して（ステップ８０８）終了
する（ステップ８０９）。
【００４５】ステップ８１４でルーチン８１８へ移行
し、ステップ８０９でルーチン８１８を終了した場合
は、その直後にステップ８１５へ移行する。ステップ８
１５では、内部変数ｉの値がＭより小さいか否かの判定
を行い、Ｍよりも小さい場合は、ｉに１を加算して（ス
テップ８１７）他の部分空間を選択する（ステップ８１
２）。ｉの値がＭ以上である場合は、ルーチン８１９を
終了する（ステップ８１６）。
し、ステップ８０９でルーチン８１８を終了した場合
は、その直後にステップ８１５へ移行する。ステップ８
１５では、内部変数ｉの値がＭより小さいか否かの判定
を行い、Ｍよりも小さい場合は、ｉに１を加算して（ス
テップ８１７）他の部分空間を選択する（ステップ８１
２）。ｉの値がＭ以上である場合は、ルーチン８１９を
終了する（ステップ８１６）。
【００４６】ステップ８０７でルーチン８１９へ移行
し、ステップ８１６でルーチン８１９を終了した場合
は、その直後にステップ８０９へ移行し、入れ子となっ
ているルーチン８１８の一つを終了する。以上の処理
を、ステップ８０６で部分空間を細分化する必要がなく
なるまで実行する。
し、ステップ８１６でルーチン８１９を終了した場合
は、その直後にステップ８０９へ移行し、入れ子となっ
ているルーチン８１８の一つを終了する。以上の処理
を、ステップ８０６で部分空間を細分化する必要がなく
なるまで実行する。
【００４７】図１に示した発明の実施例の具体的動作例
を、図２〜図５を用いて説明する。図２は部分空間と障
壁の交差面の生成過程を説明する図である。点９０１か
ら放射されたレイ９０２は、点９０３において障壁９０
６と衝突する。線９０８は障壁９０６に対する法線であ
り、角度９０７は障壁９０６に対するレイ９０２の入射
角を指す。このとき、レイ９０２近傍の部分空間９０４
と障壁９０６とによって、交差面９０５が生じる。図２
はアンテナ周囲の空間を三角錐によって分割した場合を
示す。
を、図２〜図５を用いて説明する。図２は部分空間と障
壁の交差面の生成過程を説明する図である。点９０１か
ら放射されたレイ９０２は、点９０３において障壁９０
６と衝突する。線９０８は障壁９０６に対する法線であ
り、角度９０７は障壁９０６に対するレイ９０２の入射
角を指す。このとき、レイ９０２近傍の部分空間９０４
と障壁９０６とによって、交差面９０５が生じる。図２
はアンテナ周囲の空間を三角錐によって分割した場合を
示す。
【００４８】図３は本発明における部分空間の細分化を
説明する図であり、図２においてθ＝０の場合を示して
いる。図３１（ａ）は、レイ９０２が点１２１３へ入射
した結果、障壁１２１２と部分空間の断面１２１１とに
よって生じる交差面１２１２（斜線部分）を示してい
る。点９０３において、部分空間９０４を細分するか否
かの判定基準（図１のフローのステップ８０５の判定基
準）の一例として、以下の項目が挙げられる。
説明する図であり、図２においてθ＝０の場合を示して
いる。図３１（ａ）は、レイ９０２が点１２１３へ入射
した結果、障壁１２１２と部分空間の断面１２１１とに
よって生じる交差面１２１２（斜線部分）を示してい
る。点９０３において、部分空間９０４を細分するか否
かの判定基準（図１のフローのステップ８０５の判定基
準）の一例として、以下の項目が挙げられる。
【００４９】（１）点１２１３における当該レイの電界
強度が所定値以上である： （２）部分空間断面１２１１に対する交差面１２１２の
面積比が所定値以下である： （３）部分空間断面１２１１の面積が所定の値以上であ
り、かつレイ９０２の入射角９０７が所定値以上であ
る：基本的には、上記の３項目の基準のうち（１）及び
（２）を同時に満たす場にを部分空間９０４の分割（細
分）を行うものとするが、これ等２つの基準に加えて、
（３）の基準を更に付加して、３つの基準を同時に満た
す場合に部分空間の分割（細分）を行うものとしても良
い。
強度が所定値以上である： （２）部分空間断面１２１１に対する交差面１２１２の
面積比が所定値以下である： （３）部分空間断面１２１１の面積が所定の値以上であ
り、かつレイ９０２の入射角９０７が所定値以上であ
る：基本的には、上記の３項目の基準のうち（１）及び
（２）を同時に満たす場にを部分空間９０４の分割（細
分）を行うものとするが、これ等２つの基準に加えて、
（３）の基準を更に付加して、３つの基準を同時に満た
す場合に部分空間の分割（細分）を行うものとしても良
い。
【００５０】細分した部分空間に割当てられるレイの伝
搬方向の決定方法（図１のフローのステップ８１３の方
向算出法）としては、例えば以下の方法が挙げられる。
まず、点９０３と点９０１とを結ぶ線分上の任意の点に
おいてレイ９０２に対し垂直な面で部分空間を切断す
る。その結果生じる、図４のような部分空間の断面１３
０９をＭ個（この例では、Ｍ＝４）の相似な領域１３０
５〜１３０８に分割し、点１３０１〜１３０４のように
各領域の重心点を求める。同様な重心点を部分空間の別
の切断面についても求める。対向する重心点同士を結ぶ
と、図５のように所望の伝搬方向ベクトル１９０１が得
られる。
搬方向の決定方法（図１のフローのステップ８１３の方
向算出法）としては、例えば以下の方法が挙げられる。
まず、点９０３と点９０１とを結ぶ線分上の任意の点に
おいてレイ９０２に対し垂直な面で部分空間を切断す
る。その結果生じる、図４のような部分空間の断面１３
０９をＭ個（この例では、Ｍ＝４）の相似な領域１３０
５〜１３０８に分割し、点１３０１〜１３０４のように
各領域の重心点を求める。同様な重心点を部分空間の別
の切断面についても求める。対向する重心点同士を結ぶ
と、図５のように所望の伝搬方向ベクトル１９０１が得
られる。
【００５１】なお、図３〜図５はＭ＝４の場合を示した
ものであるが、Ｍが他の値の場合も同様である。
ものであるが、Ｍが他の値の場合も同様である。
【００５２】このように部分空間９０４を細分すると、
反射点９０３における部分空間の断面は図３（ｂ）の様
になる。図９のフロー及び上記の分割判定基準に基づく
と、例えば、部分空間１２０２はこれ以上分割されず
に、部分空間１２０２に割当てられたレイは障壁１２１
２で反射することなく、点１２０６を通過して他の障壁
に衝突するまで進行を続ける。この部分空間１２０２が
分割されないのは、図１のフローにおけるステップ８０
４で反射点が検出されないからであり、詳述すると、こ
の部分空間１２０２に割当てられたレイは、点１２０６
を通過してしまい、障壁１２１２で反射しないからであ
る。
反射点９０３における部分空間の断面は図３（ｂ）の様
になる。図９のフロー及び上記の分割判定基準に基づく
と、例えば、部分空間１２０２はこれ以上分割されず
に、部分空間１２０２に割当てられたレイは障壁１２１
２で反射することなく、点１２０６を通過して他の障壁
に衝突するまで進行を続ける。この部分空間１２０２が
分割されないのは、図１のフローにおけるステップ８０
４で反射点が検出されないからであり、詳述すると、こ
の部分空間１２０２に割当てられたレイは、点１２０６
を通過してしまい、障壁１２１２で反射しないからであ
る。
【００５３】一方、部分空間１２０３は上記基準を満た
しているとして、図３（ｃ）のように取出された後、さ
らに細分され、以後同様に必要なだけ細分化が行われ
る。
しているとして、図３（ｃ）のように取出された後、さ
らに細分され、以後同様に必要なだけ細分化が行われ
る。
【００５４】次に、本実施例の効果について説明する。
本実施例では、当該部分空間と障壁との交差面と反射点
における当該部分空間の断面の面積比を分割判定の基準
とし、実際には反射が生じていない領域を分割して別扱
いにすることができる（上記判定基準の（２））。ま
た、図６に示すとおり、電波の反射率は一般に入射角が
大きくなるにつれて急増することが多いが、当該レイの
入射角度を判定基準とすることにより、本来は異なる角
度で障壁に入射し反射率が大きく異なる領域を分割して
別扱いにすることができる（上記判定基準の（３））。
本実施例では、当該部分空間と障壁との交差面と反射点
における当該部分空間の断面の面積比を分割判定の基準
とし、実際には反射が生じていない領域を分割して別扱
いにすることができる（上記判定基準の（２））。ま
た、図６に示すとおり、電波の反射率は一般に入射角が
大きくなるにつれて急増することが多いが、当該レイの
入射角度を判定基準とすることにより、本来は異なる角
度で障壁に入射し反射率が大きく異なる領域を分割して
別扱いにすることができる（上記判定基準の（３））。
【００５５】これらの作用によって、ラウンチング法よ
る推定精度を向上させることができることになるのであ
る。一方、本実施の形態では細分した後の部分空間はも
との部分空間と相似であるため、以降の細分化処理にお
いても同様な処理を繰り返し適用することができ、処理
を簡易化することが可能である。さらに、分割前の部分
領域は分割後の部分空間の和で与えられるために、アン
テナ周囲の全空間を隙間無く推定することが可能であ
る。
る推定精度を向上させることができることになるのであ
る。一方、本実施の形態では細分した後の部分空間はも
との部分空間と相似であるため、以降の細分化処理にお
いても同様な処理を繰り返し適用することができ、処理
を簡易化することが可能である。さらに、分割前の部分
領域は分割後の部分空間の和で与えられるために、アン
テナ周囲の全空間を隙間無く推定することが可能であ
る。
【００５６】以上述べた実施例の処理動作から、この処
理動作を実現するための機能ブロックが、図７に示す様
に得られることになる。図７を参照すると、衝突検出部
１と、部分空間分割判定部２と、部分空間分割部３と、
部分空間へのレイ割当部４と、制御部（ＣＰＵ）５と、
ＲＡＭ（ランダムアクセスメモリ）６と、ＲＯＭ（リー
ドオンリメモリ）７とを含んで構成される。
理動作を実現するための機能ブロックが、図７に示す様
に得られることになる。図７を参照すると、衝突検出部
１と、部分空間分割判定部２と、部分空間分割部３と、
部分空間へのレイ割当部４と、制御部（ＣＰＵ）５と、
ＲＡＭ（ランダムアクセスメモリ）６と、ＲＯＭ（リー
ドオンリメモリ）７とを含んで構成される。
【００５７】衝突検出部１はレイと遮蔽物（障壁）との
衝突を検出する機能を有しており、例えば、図３の例で
は、障壁１２１２に対してレイが衝突する関係にあり、
これは障壁１２１２の面を表す３次元空間上に規定され
る式と、同じくレイを表す３次元空間上に規定される式
とを連立させて解く手法を採用することで、衝突の有無
が判定可能である。
衝突を検出する機能を有しており、例えば、図３の例で
は、障壁１２１２に対してレイが衝突する関係にあり、
これは障壁１２１２の面を表す３次元空間上に規定され
る式と、同じくレイを表す３次元空間上に規定される式
とを連立させて解く手法を採用することで、衝突の有無
が判定可能である。
【００５８】部分空間判定部２はレイと同一視される空
間を分割するか否かの判定を、上述した判定基準に従っ
て判定する機能を有しており、部分空間分割部３は上述
した予め定められた規定に従って当該空間を複数の部分
空間に分割する機能を有する。部分空間へのレイ割当部
４は分割された部分空間の各々に対して同一視するレイ
を割当てる機能を有する。
間を分割するか否かの判定を、上述した判定基準に従っ
て判定する機能を有しており、部分空間分割部３は上述
した予め定められた規定に従って当該空間を複数の部分
空間に分割する機能を有する。部分空間へのレイ割当部
４は分割された部分空間の各々に対して同一視するレイ
を割当てる機能を有する。
【００５９】制御部５はＣＰＵ（コンピュータの中央演
算処理部に相当）からなり、ＲＯＭ７の記録媒体に予め
格納されている動作制御プログラムを読込んで、このプ
ログラムの制御手順に従って、作業用メモリであるＲＡ
Ｍ６を使用しつつ各部１〜４を制御し、図１に示した動
作処理を実行制御するものである。すなわち、そのプロ
グラムがＣＰＵの動作を制御するものであり、プログラ
ム制御されたそのＣＰＵが、当該プログラムにより指令
されて、上記の各処理を行うようになっている。
算処理部に相当）からなり、ＲＯＭ７の記録媒体に予め
格納されている動作制御プログラムを読込んで、このプ
ログラムの制御手順に従って、作業用メモリであるＲＡ
Ｍ６を使用しつつ各部１〜４を制御し、図１に示した動
作処理を実行制御するものである。すなわち、そのプロ
グラムがＣＰＵの動作を制御するものであり、プログラ
ム制御されたそのＣＰＵが、当該プログラムにより指令
されて、上記の各処理を行うようになっている。
【００６０】次に、本発明の第２の実施例について、図
面を参照して説明する。第１の実施例においては、レイ
に垂直な平面で部分空間を切断したときに、図４に示す
形態となるようにアンテナ周囲の空間を錐体で分割して
いた。これに対して、第２の実施例では、部分空間の断
面が図８となるようにアンテナ周囲の空間を分割してい
る点で異なる。図８では、正６角形１４０８と相似な正
６角形１４０９〜１４１５で部分空間を細分化するＭ＝
７の場合を示している。正６角形１４０９〜１４１４は
それぞれ辺の中点１４１６〜１４２１で正６角形１４０
８の頂点と接しており、分割後に各部分空間に割当てら
れるレイは正６角形１４０９〜１４１５の重心１４０１
〜１４０７を通る。
面を参照して説明する。第１の実施例においては、レイ
に垂直な平面で部分空間を切断したときに、図４に示す
形態となるようにアンテナ周囲の空間を錐体で分割して
いた。これに対して、第２の実施例では、部分空間の断
面が図８となるようにアンテナ周囲の空間を分割してい
る点で異なる。図８では、正６角形１４０８と相似な正
６角形１４０９〜１４１５で部分空間を細分化するＭ＝
７の場合を示している。正６角形１４０９〜１４１４は
それぞれ辺の中点１４１６〜１４２１で正６角形１４０
８の頂点と接しており、分割後に各部分空間に割当てら
れるレイは正６角形１４０９〜１４１５の重心１４０１
〜１４０７を通る。
【００６１】この第２の実施例で用いた正６角形のよう
に比較的円形に近い図形を用いた場合、受信アンテナ近
傍をレイが通過した場合にそのレイを取込むか否かの判
定を簡易化でき処理速度を向上させることができる。そ
の理由は、一般にこの判定が、当該レイ近傍の部分空間
の内部に受信アンテナが含まれるか否かを調べることに
よって行われることが多いためである。そこで、当該部
分空間を図９のように正６角形１６０２の外接円１６０
１を当該レイに垂直な切断面とする部分空間で近似する
ことによって、この判定は受信アンテナと当該レイの間
の距離を調べることに帰着でき、判定処理の簡易化が達
成できる。
に比較的円形に近い図形を用いた場合、受信アンテナ近
傍をレイが通過した場合にそのレイを取込むか否かの判
定を簡易化でき処理速度を向上させることができる。そ
の理由は、一般にこの判定が、当該レイ近傍の部分空間
の内部に受信アンテナが含まれるか否かを調べることに
よって行われることが多いためである。そこで、当該部
分空間を図９のように正６角形１６０２の外接円１６０
１を当該レイに垂直な切断面とする部分空間で近似する
ことによって、この判定は受信アンテナと当該レイの間
の距離を調べることに帰着でき、判定処理の簡易化が達
成できる。
【００６２】次に、本発明の第３の実施例について図面
を参照して説明する。当該部分空間の断面が図８となる
ようにアンテナ周囲の空間を分割した場合、分割前の断
面１４０８と分割後の部分空間の和領域を比較すると、
図８で斜線で示した部分のような過不足領域が生じる。
そこで、分割後の部分空間が分割前の部分空間の領域外
にはみ出さないような分割形態を図１０に示す。
を参照して説明する。当該部分空間の断面が図８となる
ようにアンテナ周囲の空間を分割した場合、分割前の断
面１４０８と分割後の部分空間の和領域を比較すると、
図８で斜線で示した部分のような過不足領域が生じる。
そこで、分割後の部分空間が分割前の部分空間の領域外
にはみ出さないような分割形態を図１０に示す。
【００６３】一方、分割後の部分空間と分割前の部分空
間との間に隙間が生じないように設定した分割形態を図
１１に示す。図１０では、頂点１８０１〜１８１２が正
６角形１８１３に内接するように、分割後の正６角形の
大きさを設定している。一方、図１１では、分割前の正
６角形の各辺の中点１５０１〜１５０６において、分割
後の正６角形と交差するように、分割後の正６角形の大
きさを設定している。
間との間に隙間が生じないように設定した分割形態を図
１１に示す。図１０では、頂点１８０１〜１８１２が正
６角形１８１３に内接するように、分割後の正６角形の
大きさを設定している。一方、図１１では、分割前の正
６角形の各辺の中点１５０１〜１５０６において、分割
後の正６角形と交差するように、分割後の正６角形の大
きさを設定している。
【００６４】特に、図１１の分割形態を用いた場合は、
第１の実施例と同様に、アンテナ周囲の空間の電波伝搬
状況を隙間なく推定することができる。
第１の実施例と同様に、アンテナ周囲の空間の電波伝搬
状況を隙間なく推定することができる。
【００６５】部分空間分割の例として上述の様に種々示
しているが、要は、円錐に内接する錐体のうち、この円
錐との体積誤差が小さい錐体を設け、その中心線をレイ
の経路とすることによりレイ近傍の空間を当該錐体によ
って形成し、レイ近傍の空間を形成するこの錐体に相似
な複数の錐体を用いてレイ近傍の空間を複数の部分空間
に分割する際に、レイ近傍の空間を形成する錐体と、こ
の錐体を分割する錐体に相似な複数の錐体の和との体積
誤差が小さくなるように、レイの近傍に複数の錐体を配
置すれば良い。
しているが、要は、円錐に内接する錐体のうち、この円
錐との体積誤差が小さい錐体を設け、その中心線をレイ
の経路とすることによりレイ近傍の空間を当該錐体によ
って形成し、レイ近傍の空間を形成するこの錐体に相似
な複数の錐体を用いてレイ近傍の空間を複数の部分空間
に分割する際に、レイ近傍の空間を形成する錐体と、こ
の錐体を分割する錐体に相似な複数の錐体の和との体積
誤差が小さくなるように、レイの近傍に複数の錐体を配
置すれば良い。
【００６６】
【発明の効果】以上述べた様に、本発明によれば、従来
の方法に比べてラウンチング法におけるレイの密度を柔
軟に制御し、推定精度の向上と計算時間の短縮を効率よ
く両立できるという効果がある。その理由は、従来の方
法では伝搬距離に応じて一律にレイの密度を制御してい
たのに対し、本発明では、反射・透過の際にレイの密度
が過小で推定精度の劣化を招くと考えられる場合に、劣
化が小さいと判断されるまでレイの密度を増大させるた
めである。これにより、推定精度を損なわない範囲で不
必要なレイの生成による計算時間の増大を抑えることが
できる。
の方法に比べてラウンチング法におけるレイの密度を柔
軟に制御し、推定精度の向上と計算時間の短縮を効率よ
く両立できるという効果がある。その理由は、従来の方
法では伝搬距離に応じて一律にレイの密度を制御してい
たのに対し、本発明では、反射・透過の際にレイの密度
が過小で推定精度の劣化を招くと考えられる場合に、劣
化が小さいと判断されるまでレイの密度を増大させるた
めである。これにより、推定精度を損なわない範囲で不
必要なレイの生成による計算時間の増大を抑えることが
できる。
【図１】本発明の実施例の動作を示すフロー図である。
【図２】部分空間と障壁の交差面の生成過程を説明する
図である。
図である。
【図３】本発明における部分空間の細分化を説明する図
である。
である。
【図４】本発明の第１の実施例における部分空間の分割
例を示す断面図である。
例を示す断面図である。
【図５】分割された部分空間に割当てられたレイの伝搬
方向を求める過程を説明する図である。
方向を求める過程を説明する図である。
【図６】反射率の入射角特性の例を示す図である。
【図７】本発明の機能ブロックの例を示す図である。
【図８】本発明の第２の実施例における部分空間の分割
例を示す断面図である。
例を示す断面図である。
【図９】部分空間の円錐近似を説明する図である。
【図１０】本発明の第３の実施例における部分空間の分
割の一例を示す断面図である。
割の一例を示す断面図である。
【図１１】本発明の第３の実施例における部分空間の分
割の他の例を示す断面図である。
割の他の例を示す断面図である。
【図１２】送信アンテナ周囲の閉領域の例を示す図であ
る。
る。
【図１３】送信アンテナ周囲の閉領域の分割過程を説明
する図である。
する図である。
【図１４】離散的間隔によるレイの放射を説明する図で
ある。
ある。
【図１５】レイと部分空間の例を示す図である。
【図１６】従来の手法による部分空間の分割を説明する
図である。
図である。
【図１７】ラウンチング法の説明図である。
【図１８】一本のレイに対する観測点におけるパスプロ
ファイルの例を示す図である。
ファイルの例を示す図である。
【図１９】ラウンチング法によって得られる観測点にお
けるパスプロファイルの例を示す図である。
けるパスプロファイルの例を示す図である。
１ 衝突検出部 ２ 部分空間分割判定部 ３ 部分空間分割部 ４ 部分空間へのレイ割当部 ５ 制御部 ６ ＲＡＭ ７ ＲＯＭ ００１，００２ 物体（遮蔽物） ００３〜０１４ レイ伝搬方向 ０１５ 送信アンテナ位置 ０１６ 受信点 ０１８，０１９，９０１，９０３ 反射・透過点 ９０２，１９０１ レイ ９０４ 部分空間 ９０５ 交差面 ９０６ 障壁
Claims (23)
- 【請求項１】 ３次元空間上に規定される観測領域内に
複数の物体が与えられ、前記観測領域内に電波線（レ
イ）が与えられ、前記物体のうち前記レイを遮蔽する遮
蔽物の検出を行って電波伝搬特性の推定をなす電波伝搬
特性推定方法であって、 前記レイが前記遮蔽物と衝突する場面において、前記レ
イを包含する閉空間を複数の部分空間に分割し、これ等
各部分空間により包含されるレイを新たに生成すること
を特徴とする電波伝搬特性推定方法。 - 【請求項２】 ３次元空間上に規定される観測領域内に
複数の物体が与えられ、前記観測領域内に電波線（レ
イ）が与えられると共に、前記レイの近傍に前記レイと
同一視される空間が与えられ、前記物体のうち前記レイ
を遮蔽する遮蔽物と前記レイが衝突した際に、予め定め
られた基準に基づき前記レイの空間分解能の調節を行う
電波伝搬特性推定におけるレイの空間分解能制御方法で
あって、 前記遮蔽物と前記レイとの衝突点を検出する第一のステ
ップと、 前記レイと同一視される前記空間を分割するか否かの判
定を、所定判定基準に従って行う第二のステップと、 この第二のステップで前記レイと同一視される前記空間
を分割すると判定された場合に、前記レイと同一視され
る前記空間を複数の部分空間に分割する第三のステップ
と、 この第三のステップにおいて分割された前記部分空間の
各々に対して同一視するレイを割当てる第四のステップ
と、を含むことを特徴とするレイの空間分解能制御方
法。 - 【請求項３】 前記第三のステップにおける分割数はＭ
個（Ｍは２以上の整数）に設定されており、ｉ番目（ｉ
は１〜Ｍ）の前記部分空間に対して前第四のステップに
おいてレイを割当てた後、割当てた各レイに対して再び
前記第一から第四のステップを適用し、前記第二のステ
ップにおいて当該レイと同一視される前記空間を分割し
ないと判定されるまで前記第一から第四のステップを、
順次繰返すことを特徴とする請求項２記載のレイの空間
分解能制御方法。 - 【請求項４】 前記第二のステップにおける前記判定基
準は、前記遮蔽物との衝突点における前記レイの電界強
度が所定の値以上であるか否かを調べる第一の判定基準
と、前記遮蔽物の衝突面を含む平面による前記レイ近傍
の前記空間の断面に対する前記遮蔽物の衝突面の面積比
が所定値以下であるか否かを調べる第二の判定基準とか
らなることを特徴とする請求項２または３記載のレイの
空間分解能制御方法。 - 【請求項５】 前記第二のステップにおける前記判定基
準は、前記断面の面積が所定の値以上であり、かつ前記
レイと前記衝突面の法線との成す角（入射角）が所定値
以上であるか否かを調べる第三の判定基準を更に含むこ
とを特徴とする請求項４記載のレイの空間分解能制御方
法。 - 【請求項６】 円錐に内接する錐体のうち、前記円錐と
の体積誤差が小さい錐体を設け、その中心線を前記レイ
の経路とすることにより前記レイ近傍の前記空間を前記
錐体によって形成し、 前記第三のステップにおいて、前記レイ近傍の空間を形
成する前記錐体に相似な複数の錐体を用いて前記レイ近
傍の空間を複数の部分空間に分割する際に、前記レイ近
傍の空間を形成する前記錐体と、この錐体を分割する前
記錐体に相似な前記複数の錐体の和との体積誤差が小さ
くなるように、前記レイの近傍に前記複数の錐体を配置
することを特徴とする請求項２〜５いずれか記載のレイ
の空間分解能制御方法。 - 【請求項７】 前記第三のステップにおいて、円錐に内
接する前記錐体として三角錐を用い、前記三角錐に相似
な複数の同一三角錐を前記三角錐内部に隙間無く配置で
きるように前記相似三角錐の大きさを設定することを特
徴とする請求項２〜５いずれか記載のレイの空間分解能
制御方法。 - 【請求項８】 前記第三のステップにおいて、円錐に内
接する前記錐体として正六角錐を用い、前記正六角錐に
相似な正六角錐を前記正六角錐の中心に配置し、その前
記相似正六角錐の各側面に前記相似正六角錐と同一な６
個の正六角錐を各々側面が接するように配置し、全７個
の前記相似正六角錐の和の体積と、もとの前記正六角錐
の体積誤差が最小となるように前記相似六角錐の大きさ
を設定し配置することを特徴とする請求項２〜５いずれ
か記載のレイの空間分解能制御方法。 - 【請求項９】 前記第三のステップにおいて、円錐に内
接する前記錐体として正六角錐を用い、前記正六角錐に
相似な正六角錐を前記正六角錐の中心に配置し、その前
記相似正六角錐の各側面に前記相似正六角錐と同一な６
個の正六角錐を各々側面が接するように配置し、全７個
の前記相似正六角錐の和からなる錐体の最外郭側面がも
との前記正六角錐の側面に内接するように前記相似六角
錐の大きさを設定し配置することを特徴とする請求項２
〜５いずれか記載のレイの空間分解能制御方法。 - 【請求項１０】 前記第三のステップにおいて、円錐に
内接する前記錐体として正六角錐を用い、前記正六角錐
に相似な正六角錐を前記正六角錐の中心に配置し、その
前記相似正六角錐の各側面に前記相似正六角錐と同一な
６個の正六角錐を各々側面が接するように配置し、全７
個の前記相似正六角錐の和からなる錐体の最外郭側面
に、もとの前記正六角錐が内接するように前記相似六角
錐の大きさを設定し配置することを特徴とする請求項２
〜５いずれか記載のレイの空間分解能制御方法。 - 【請求項１１】 ３次元空間上に規定される観測領域内
に複数の物体が与えられ、前記観測領域内に電波線（レ
イ）が与えられ、前記物体のうち前記レイを遮蔽する遮
蔽物の検出を行って電波伝搬特性の推定をなす電波伝搬
特性推定システムであって、 前記レイが前記遮蔽物と衝突する場面において、前記レ
イを包含する閉空間を複数の部分空間に分割し、これ等
各部分空間により包含されるレイを新たに生成すること
を特徴とする電波伝搬特性推定システム。 - 【請求項１２】 ３次元空間上に規定される観測領域内
に複数の物体が与えられ、前記観測領域内に電波線（レ
イ）が与えられると共に、前記レイの近傍に前記レイと
同一視される空間が与えられ、前記物体のうち前記レイ
を遮蔽する遮蔽物と前記レイが衝突した際に、予め定め
られた基準に基づき前記レイの空間分解能の調節を行う
電波伝搬特性推定におけるレイの空間分解能制御システ
ムであって、 前記遮蔽物と前記レイとの衝突点を検出する衝突点検出
手段と、 前記レイと同一視される前記空間を分割するか否かの判
定を、所定判定基準に従って行う判定手段と、 この判定手段で前記レイと同一視される前記空間を分割
すると判定された場合に、前記レイと同一視される前記
空間を複数の部分空間に分割する分割手段と、 この分割手段において分割された前記部分空間の各々に
対して同一視するレイを割当てる割当手段と、を含むこ
とを特徴とするレイの空間分解能制御システム。 - 【請求項１３】 前記分割手段における分割数はＭ個
（Ｍは２以上の整数）に設定されており、ｉ番目（ｉは
１〜Ｍ）の前記部分空間に対して前記割当手段において
レイを割当てた後、割当てた各レイに対して再び前記衝
突点検出手段、前記判定手段、前記分割手段及び前記割
当手段を実行制御し、前記判定手段において当該レイと
同一視される前記空間を分割しないと判定されるまで、
前記衝突点検出手段、前記判定手段、前記分割手段及び
前記割当手段を、順次繰返し実行制御する制御手段とを
含むことを特徴とする請求項１２記載のレイの空間分解
能制御システム。 - 【請求項１４】 前記判定手段における前記判定基準
は、前記遮蔽物との衝突点における前記レイの電界強度
が所定の値以上であるか否かを調べる第一の判定基準
と、前記遮蔽物の衝突面を含む平面による前記レイ近傍
の前記空間の断面に対する前記遮蔽物の衝突面の面積比
が所定値以下であるか否かを調べる第二の判定基準とか
らなることを特徴とする請求項１２または１３記載のレ
イの空間分解能制御システム。 - 【請求項１５】 前記判定手段における前記判定基準
は、前記断面の面積が所定の値以上であり、かつ前記レ
イと前記衝突面の法線との成す角（入射角）が所定値以
上であるか否かを調べる第三の判定基準を更に含むこと
を特徴とする請求項１４記載のレイの空間分解能制御シ
ステム。 - 【請求項１６】 円錐に内接する錐体のうち、前記円錐
との体積誤差が小さい錐体を設け、その中心線を前記レ
イの経路とすることにより前記レイ近傍の前記空間を前
記錐体によって形成し、 前記分割手段において、前記レイ近傍の空間を形成する
前記錐体に相似な複数の錐体を用いて前記レイ近傍の空
間を複数の部分空間に分割する際に、前記レイ近傍の空
間を形成する前記錐体と、この錐体を分割する前記錐体
に相似な前記複数の錐体の和との体積誤差が小さくなる
ように、前記レイの近傍に前記複数の錐体を配置するこ
とを特徴とする請求項１２〜１５いずれか記載のレイの
空間分解能制御システム。 - 【請求項１７】 前記分割手段において、円錐に内接す
る前記錐体として三角錐を用い、前記三角錐に相似な複
数の同一三角錐を前記三角錐内部に隙間無く配置できる
ように前記相似三角錐の大きさを設定することを特徴と
する請求項１２〜１５いずれか記載のレイの空間分解能
制御システム。 - 【請求項１８】 前記分割手段において、円錐に内接す
る前記錐体として正六角錐を用い、前記正六角錐に相似
な正六角錐を前記正六角錐の中心に配置し、その前記相
似正六角錐の各側面に前記相似正六角錐と同一な６個の
正六角錐を各々側面が接するように配置し、全７個の前
記相似正六角錐の和の体積と、もとの前記正六角錐の体
積誤差が最小となるように前記相似六角錐の大きさを設
定し配置することを特徴とする請求項１２〜１５いずれ
か記載のレイの空間分解能制御システム。 - 【請求項１９】 前記分割手段において、円錐に内接す
る前記錐体として正六角錐を用い、前記正六角錐に相似
な正六角錐を前記正六角錐の中心に配置し、その前記相
似正六角錐の各側面に前記相似正六角錐と同一な６個の
正六角錐を各々側面が接するように配置し、全７個の前
記相似正六角錐の和からなる錐体の最外郭側面がもとの
前記正六角錐の側面に内接するように前記相似六角錐の
大きさを設定し配置することを特徴とする請求項１２〜
１５いずれか記載のレイの空間分解能制御システム。 - 【請求項２０】 前記分割手段において、円錐に内接す
る前記錐体として正六角錐を用い、前記正六角錐に相似
な正六角錐を前記正六角錐の中心に配置し、その前記相
似正六角錐の各側面に前記相似正六角錐と同一な６個の
正六角錐を各々側面が接するように配置し、全７個の前
記相似正六角錐の和からなる錐体の最外郭側面に、もと
の前記正六角錐が内接するように前記相似六角錐の大き
さを設定し配置することを特徴とする請求項１２〜１５
いずれか記載のレイの空間分解能制御システム。 - 【請求項２１】 ３次元空間上に規定される観測領域内
に複数の物体が与えられ、前記観測領域内に電波線（レ
イ）が与えられ、前記物体のうち前記レイを遮蔽する遮
蔽物の検出を行って電波伝搬特性の推定をなす電波伝搬
特性推定方法をコンピュータに実行させるためのプログ
ラムであって、 前記レイが前記遮蔽物と衝突する場面において、前記レ
イを包含する閉空間を複数の部分空間に分割し、これ等
各部分空間により包含されるレイを新たに生成すること
を特徴とするプログラム。 - 【請求項２２】 ３次元空間上に規定される観測領域内
に複数の物体が与えられ、前記観測領域内に電波線（レ
イ）が与えられると共に、前記レイの近傍に前記レイと
同一視される空間が与えられ、前記物体のうち前記レイ
を遮蔽する遮蔽物と前記レイが衝突した際に、予め定め
られた基準に基づき前記レイの空間分解能の調節を行う
電波伝搬特性推定におけるレイの空間分解能制御方法を
コンピュータに実行させるためのプログラムであって、 前記遮蔽物と前記レイとの衝突点を検出する第一のステ
ップと、 前記レイと同一視される前記空間を分割するか否かの判
定を、所定判定基準に従って行う第二のステップと、 この第二のステップで前記レイと同一視される前記空間
を分割すると判定された場合に、前記レイと同一視され
る前記空間を複数の部分空間に分割する第三のステップ
と、 この第三のステップにおいて分割された前記部分空間の
各々に対して同一視するレイを割当てる第四のステップ
と、を含むことを特徴とするプログラム。 - 【請求項２３】 前記第三のステップにおける分割数は
Ｍ個（Ｍは２以上の整数）に設定されており、ｉ番目
（ｉは１〜Ｍ）の前記部分空間に対して前記第四のステ
ップにおいてレイを割当てた後、割当てた各レイに対し
て再び前記第一から第四のステップを適用し、前記第二
のステップにおいて当該レイと同一視される前記空間を
分割しないと判定されるまで前記第一から第四のステッ
プを、順次繰返すことを特徴とする請求項２２記載のプ
ログラム。
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 true JP2002214267A (ja) | 2002-07-31 |
| JP4273665B2 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006203526A (ja) * | 2005-01-20 | 2006-08-03 | Canon Inc | 情報処理方法および装置 |
| JP2006287685A (ja) * | 2005-04-01 | 2006-10-19 | Hitachi Ltd | 電波伝搬の推定プログラム、電波伝搬の推定方法、この方法を実行する装置 |
| JPWO2005088868A1 (ja) * | 2004-03-17 | 2008-01-31 | 日本電気株式会社 | 電波伝搬特性推定システム及びその方法並びにプログラム |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPWO2005088868A1 (ja) * | 2004-03-17 | 2008-01-31 | 日本電気株式会社 | 電波伝搬特性推定システム及びその方法並びにプログラム |
| US7634265B2 (en) | 2004-03-17 | 2009-12-15 | Nec Corporation | Radio wave propagation characteristic estimation system, and its method and program |
| JP2006203526A (ja) * | 2005-01-20 | 2006-08-03 | Canon Inc | 情報処理方法および装置 |
| JP4481840B2 (ja) * | 2005-01-20 | 2010-06-16 | キヤノン株式会社 | 情報処理方法および装置 |
| JP2006287685A (ja) * | 2005-04-01 | 2006-10-19 | Hitachi Ltd | 電波伝搬の推定プログラム、電波伝搬の推定方法、この方法を実行する装置 |
| JP4530898B2 (ja) * | 2005-04-01 | 2010-08-25 | 株式会社日立製作所 | 電波伝搬の推定プログラム、電波伝搬の推定方法、この方法を実行する装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4273665B2 (ja) | 2009-06-03 |
| US20020094809A1 (en) | 2002-07-18 |
| US6876851B2 (en) | 2005-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4819303B2 (ja) | 移動通信システムにおける基地局設置設計方法及び基地局設置設計装置並びにプログラム | |
| JP2002214267A (ja) | 電波伝搬特性推定方法及びレイの空間分解能制御方法並びにシステム及びプログラム | |
| KR100288560B1 (ko) | 도심지 협곡 모델에서 편파 효과를 고려한 전파의 전파특성 예측 방법 | |
| Wölfle et al. | A fast and enhanced ray optical propagation model for indoor and urban scenarios, based on an intelligent preprocessing of the database | |
| JPH10503915A (ja) | 無線通信システム計画のための方法 | |
| US8285219B2 (en) | Apparatus and method for searching for wireless station locations | |
| JPH07288495A (ja) | 電磁波伝播を迅速に予測するための技法 | |
| Xiong et al. | Customizing indoor wireless coverage via 3D-fabricated reflectors | |
| Hussain | Efficient ray-tracing algorithms for radio wave propagation in urban environments | |
| US7933558B2 (en) | Radio wave arrival state estimation system, radio wave arrival state estimation method, and program | |
| JP3823917B2 (ja) | レイの受信判定方法並びにシステム及びそれを用いた電波伝搬特性推定方法 | |
| US7379710B2 (en) | Radio-wave propagation characteristic forecasting system and its method, and program | |
| Kim et al. | Accelerated three dimensional ray tracing techniques using ray frustums for wireless propagation models | |
| CN109861775B (zh) | 一种传播路径搜索方法以及装置 | |
| Hussain et al. | An efficient ray tracing method for propagation prediction along a mobile route in urban environments | |
| US8358305B2 (en) | Method for 3D preprocessing for ray tracing | |
| JP3654197B2 (ja) | 電波伝搬特性予測システム及びその方法並びにプログラム | |
| Hussain et al. | An image visibility based pre-processing method for fast ray tracing in urban environments | |
| KR100329518B1 (ko) | 무선망 설계시스템에서의 전파해석방법 | |
| Kreuzgruber et al. | Prediction of indoor radio propagation with the ray splitting model including edge diffraction and rough surfaces | |
| JP2002107397A (ja) | 電波伝搬特性予測における遮蔽物の検出方法及びシステム並びに記録媒体 | |
| JP7221333B2 (ja) | 電波反射部材の設置位置を決定する情報処理装置、方法及びプログラム | |
| Wu et al. | Efficient Allocation of Intelligent Surfaces in Programmable Wireless Environments | |
| Schweins et al. | Optimization and Evaluation of a 3-D Ray Tracing Channel Predictor Individually for Each Propagation Effect | |
| Ben-Moshe | Geometric heuristics for rural radio maps approximation |
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 |