JP2002064440A - 電波監視方法及び装置 - Google Patents
電波監視方法及び装置Info
- Publication number
- JP2002064440A JP2002064440A JP2000250028A JP2000250028A JP2002064440A JP 2002064440 A JP2002064440 A JP 2002064440A JP 2000250028 A JP2000250028 A JP 2000250028A JP 2000250028 A JP2000250028 A JP 2000250028A JP 2002064440 A JP2002064440 A JP 2002064440A
- Authority
- JP
- Japan
- Prior art keywords
- radio wave
- point
- emitted
- electric field
- radio
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
(57)【要約】
【課題】 電波伝搬シミュレーションによるデータベー
スの作成時間を短縮し、電波監視の効率を向上する電波
監視方法及び装置を提供する。 【解決手段】 センサ局１２ａにおいて電波の到来方向
及び強度のパターンを観測し、センサ局１２ａにおける
観測結果を、センサ局１２ａから電波を発射したときに
他の複数の地点において観測される電波の到来方向及び
強度のパターンのシミュレーション結果と比較し、複数
の地点のうちセンサ局１２ａにおける前記観測結果と最
も相関の強い電波の到来方向及び強度のパターンを与え
る地点を電波発射源と特定する。
スの作成時間を短縮し、電波監視の効率を向上する電波
監視方法及び装置を提供する。 【解決手段】 センサ局１２ａにおいて電波の到来方向
及び強度のパターンを観測し、センサ局１２ａにおける
観測結果を、センサ局１２ａから電波を発射したときに
他の複数の地点において観測される電波の到来方向及び
強度のパターンのシミュレーション結果と比較し、複数
の地点のうちセンサ局１２ａにおける前記観測結果と最
も相関の強い電波の到来方向及び強度のパターンを与え
る地点を電波発射源と特定する。
Description
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は、一定の地域内で不
法無線局などの電波発射源の位置を特定し、電波利用環
境を監視しうる電波監視方法及び装置に関する。
法無線局などの電波発射源の位置を特定し、電波利用環
境を監視しうる電波監視方法及び装置に関する。
【０００２】
【従来の技術】従来、ある地域内にある不法無線局など
の電波発射源の位置を特定する際には、複数の地点にセ
ンサ局（監視局）を配置して、各センサ局において八木
・宇田アンテナやゴニオメータなどを用いてその電波発
射源からの電波の到来方向を観測し、それぞれのセンサ
局での到来方向を地図上にプロットしてその交点の位置
からその電波発射源の位置を推定していた。
の電波発射源の位置を特定する際には、複数の地点にセ
ンサ局（監視局）を配置して、各センサ局において八木
・宇田アンテナやゴニオメータなどを用いてその電波発
射源からの電波の到来方向を観測し、それぞれのセンサ
局での到来方向を地図上にプロットしてその交点の位置
からその電波発射源の位置を推定していた。
【０００３】上述した従来の方法は、監視対象の電波発
射源からの電波がセンサ局まで単一の経路で直進して伝
搬することを仮定している。しかしながら、電波発射源
からの電波は、実際には、地形やビルなどの地物による
回折や反射を経て、センサ局に伝搬してくることがあ
る。また、電波発射源からの電波は、地形や地物の影響
で複数の伝搬経路（マルチパス）を通って、センサ局に
到来することがある。従来の方法によって監視対象の電
波発射源の位置を推定した場合、電波が回折や反射を経
て伝搬することや複数の経路で伝搬することにより、電
波発射源の位置の推定精度が悪くなるという問題点を生
じる。
射源からの電波がセンサ局まで単一の経路で直進して伝
搬することを仮定している。しかしながら、電波発射源
からの電波は、実際には、地形やビルなどの地物による
回折や反射を経て、センサ局に伝搬してくることがあ
る。また、電波発射源からの電波は、地形や地物の影響
で複数の伝搬経路（マルチパス）を通って、センサ局に
到来することがある。従来の方法によって監視対象の電
波発射源の位置を推定した場合、電波が回折や反射を経
て伝搬することや複数の経路で伝搬することにより、電
波発射源の位置の推定精度が悪くなるという問題点を生
じる。
【０００４】そこで、本願発明者は、上記従来の方法に
おける問題点を解決する広域電波監視方法及び装置をす
でに提案している（例えば、特開平１１−３２６４８２
号公報を参照）。
おける問題点を解決する広域電波監視方法及び装置をす
でに提案している（例えば、特開平１１−３２６４８２
号公報を参照）。
【０００５】特開平１１−３２６４８２号公報に記載の
広域電波監視方法では、一定の地域内に１以上のセンサ
局を配置し、地域内の複数の地点の各々について、当該
地点に電波発射源があるとしたときのセンサ局での電波
の到来方向を、地図情報を利用する計算機シミュレーシ
ョンによって算出する。一方、センサ局において、伝搬
経路ごとに分離して電波発射源からの電波の到来方向を
観測する。センサ局によって観測した電波の到来方向と
計算機シミュレーションの結果とを比較し、計算機シミ
ュレーションの結果の中から、観測した電波の到来方向
と最も類似性が高い到来方向を示すものを探索して対応
する地点を候補位置として判別する。更に、候補位置の
周辺に電波発射源があるとして電波発射源の位置を調整
しながら計算機シミュレーションを実行して電波発射源
の位置を決定する。
広域電波監視方法では、一定の地域内に１以上のセンサ
局を配置し、地域内の複数の地点の各々について、当該
地点に電波発射源があるとしたときのセンサ局での電波
の到来方向を、地図情報を利用する計算機シミュレーシ
ョンによって算出する。一方、センサ局において、伝搬
経路ごとに分離して電波発射源からの電波の到来方向を
観測する。センサ局によって観測した電波の到来方向と
計算機シミュレーションの結果とを比較し、計算機シミ
ュレーションの結果の中から、観測した電波の到来方向
と最も類似性が高い到来方向を示すものを探索して対応
する地点を候補位置として判別する。更に、候補位置の
周辺に電波発射源があるとして電波発射源の位置を調整
しながら計算機シミュレーションを実行して電波発射源
の位置を決定する。
【０００６】上記の方法によれば、回折や反射が存在す
るような場合であっても、電波発射源の位置を正確に決
定することが可能となる。
るような場合であっても、電波発射源の位置を正確に決
定することが可能となる。
【０００７】
【発明が解決しようとする課題】しかしながら、特開平
１１−３２６４８２号公報に記載の広域電波監視方法及
び装置では、地域内の複数の地点の各々について、当該
地点に電波発射源があるとしたときのセンサ局での電波
の到来方向を計算機シミュレーションによって算出する
ため、不特定な電波発射源に対して、電波伝搬の計算機
シミュレーションによるデータベースを複数組用意する
必要があった。更に、これらのデータベースを作成する
のに時間を要してしまっていた。
１１−３２６４８２号公報に記載の広域電波監視方法及
び装置では、地域内の複数の地点の各々について、当該
地点に電波発射源があるとしたときのセンサ局での電波
の到来方向を計算機シミュレーションによって算出する
ため、不特定な電波発射源に対して、電波伝搬の計算機
シミュレーションによるデータベースを複数組用意する
必要があった。更に、これらのデータベースを作成する
のに時間を要してしまっていた。
【０００８】本発明の目的は、電波伝搬シミュレーショ
ンによるデータベースの作成時間を短縮し、電波監視の
効率を向上する電波監視方法及び装置を提供することに
ある。
ンによるデータベースの作成時間を短縮し、電波監視の
効率を向上する電波監視方法及び装置を提供することに
ある。
【０００９】
【課題を解決するための手段】上記目的は、一の地点に
おいて電波の到来方向及び強度のパターンを観測し、前
記一の地点における観測結果を、前記一の地点から電波
を発射したときに他の複数の地点において観測される電
波の到来方向及び強度のパターンのシミュレーション結
果と比較し、前記複数の地点のうち前記一の地点におけ
る前記観測結果と最も相関の強い電波の到来方向及び強
度のパターンを与える地点を電波発射源と特定すること
を特徴とする電波監視方法により達成される。
おいて電波の到来方向及び強度のパターンを観測し、前
記一の地点における観測結果を、前記一の地点から電波
を発射したときに他の複数の地点において観測される電
波の到来方向及び強度のパターンのシミュレーション結
果と比較し、前記複数の地点のうち前記一の地点におけ
る前記観測結果と最も相関の強い電波の到来方向及び強
度のパターンを与える地点を電波発射源と特定すること
を特徴とする電波監視方法により達成される。
【００１０】また、上記の電波監視方法において、前記
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、観測地
域を２次元的に複数の領域に分割し、分割した各領域に
おいて観測される電界強度を計算するようにしてもよ
い。
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、観測地
域を２次元的に複数の領域に分割し、分割した各領域に
おいて観測される電界強度を計算するようにしてもよ
い。
【００１１】また、上記の電波監視方法において、前記
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、観測地
域を３次元的に複数の空間に分割し、分割した各空間に
おいて観測される電界強度を計算するようにしてもよ
い。
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、観測地
域を３次元的に複数の空間に分割し、分割した各空間に
おいて観測される電界強度を計算するようにしてもよ
い。
【００１２】また、上記の電波監視方法において、前記
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、前記一
の地点から発射される電波の発射方向を変えて前記各領
域又は前記各空間において観測される電界強度をそれぞ
れ計算し、電波の発射方向毎に電界強度分布を求めるよ
うにしてもよい。
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、前記一
の地点から発射される電波の発射方向を変えて前記各領
域又は前記各空間において観測される電界強度をそれぞ
れ計算し、電波の発射方向毎に電界強度分布を求めるよ
うにしてもよい。
【００１３】また、上記の電波監視方法において、前記
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、前記観
測地域の地形や地物を考慮して電波の到来経路を計算す
るようにしてもよい。
一の地点から電波を発射したときの電波の到来方向及び
強度のパターンをシミュレーションする際には、前記観
測地域の地形や地物を考慮して電波の到来経路を計算す
るようにしてもよい。
【００１４】また、上記の電波監視方法において、特定
された前記電波発射源の位置と前記シミュレーション結
果に基づいて、前記一の地点から前記電波発射源までの
電波の伝搬経路をトレースするようにしてもよい。
された前記電波発射源の位置と前記シミュレーション結
果に基づいて、前記一の地点から前記電波発射源までの
電波の伝搬経路をトレースするようにしてもよい。
【００１５】また、上記の電波監視方法において、前記
伝搬経路のトレース結果に基づき、前記電波発射源のア
ンテナ指向性を推定するようにしてもよい。
伝搬経路のトレース結果に基づき、前記電波発射源のア
ンテナ指向性を推定するようにしてもよい。
【００１６】また、上記の電波監視方法において、推定
された前記電波発射源のアンテナ指向性に基づき、前記
電波発射源から発せられた電波の電界強度分布を求める
ようにしてもよい。
された前記電波発射源のアンテナ指向性に基づき、前記
電波発射源から発せられた電波の電界強度分布を求める
ようにしてもよい。
【００１７】また、上記目的は、一の地点に設けられ、
電波の到来方向及び強度のパターンを観測する電波観測
手段と、前記一の地点から電波を発射したときに他の複
数の地点において観測される電波の到来方向及び強度の
パターンのシミュレーション結果を記憶する記憶手段
と、前記電波観測手段における観測結果と、前記記憶手
段に記憶されているシミュレーション結果とを比較し、
前記複数の地点のうち前記電波観測手段における前記観
測結果と最も相関の強い電波の到来方向及び強度のパタ
ーンを与える地点を電波発射源と特定する電波発射源特
定手段とを有することを特徴とする電波監視装置により
達成される。
電波の到来方向及び強度のパターンを観測する電波観測
手段と、前記一の地点から電波を発射したときに他の複
数の地点において観測される電波の到来方向及び強度の
パターンのシミュレーション結果を記憶する記憶手段
と、前記電波観測手段における観測結果と、前記記憶手
段に記憶されているシミュレーション結果とを比較し、
前記複数の地点のうち前記電波観測手段における前記観
測結果と最も相関の強い電波の到来方向及び強度のパタ
ーンを与える地点を電波発射源と特定する電波発射源特
定手段とを有することを特徴とする電波監視装置により
達成される。
【００１８】また、上記の電波監視装置において、前記
記憶手段は、観測地域を２次元的に複数の領域に分割
し、前記一の地点から電波を発射したときに、分割した
各領域において観測される電界強度を計算するようにし
てもよい。
記憶手段は、観測地域を２次元的に複数の領域に分割
し、前記一の地点から電波を発射したときに、分割した
各領域において観測される電界強度を計算するようにし
てもよい。
【００１９】また、上記の電波監視装置において、前記
記憶手段は、観測地域を３次元的に複数の空間に分割
し、前記一の地点から電波を発射したときに、分割した
各空間において観測される電界強度を計算するようにし
てもよい。
記憶手段は、観測地域を３次元的に複数の空間に分割
し、前記一の地点から電波を発射したときに、分割した
各空間において観測される電界強度を計算するようにし
てもよい。
【００２０】また、上記の電波監視装置において、前記
記憶手段は、前記一の地点から発射される電波の発射方
向を変えたときに前記各領域又は前記各空間において観
測される電界強度を、電波の発射方向毎に記憶するよう
にしてもよい。
記憶手段は、前記一の地点から発射される電波の発射方
向を変えたときに前記各領域又は前記各空間において観
測される電界強度を、電波の発射方向毎に記憶するよう
にしてもよい。
【００２１】また、上記の電波監視装置において、前記
電波発射源特定手段により特定された前記電波発射源の
位置と前記シミュレーション結果に基づいて前記一の地
点から前記電波発射源までの電波の伝搬経路をトレース
する伝搬経路トレース手段を更に有するようにしてもよ
い。
電波発射源特定手段により特定された前記電波発射源の
位置と前記シミュレーション結果に基づいて前記一の地
点から前記電波発射源までの電波の伝搬経路をトレース
する伝搬経路トレース手段を更に有するようにしてもよ
い。
【００２２】また、上記の電波監視装置において、前記
伝播経路トレース手段による前記伝搬経路のトレース結
果に基づいて前記電波発射源のアンテナ指向性を推定す
るアンテナ指向性推定手段を更に有するようにしてもよ
い。
伝播経路トレース手段による前記伝搬経路のトレース結
果に基づいて前記電波発射源のアンテナ指向性を推定す
るアンテナ指向性推定手段を更に有するようにしてもよ
い。
【００２３】また、上記の電波監視装置において、前記
アンテナ指向性推定手段により推定された前記電波発射
源のアンテナ指向性に基づいて前記電波発射源から発せ
られた電波の電界強度分布を求める電界強度分布算出手
段を更に有するようにしてもよい。
アンテナ指向性推定手段により推定された前記電波発射
源のアンテナ指向性に基づいて前記電波発射源から発せ
られた電波の電界強度分布を求める電界強度分布算出手
段を更に有するようにしてもよい。
【００２４】
【発明の実施の形態】本発明の一実施形態による電波監
視方法及び装置について図１乃至図６を用いて説明す
る。図１は、本実施形態による電波監視方法を示す概略
図、図２は、電波監視装置の構成を示すブロック図、図
３は、電波伝搬シミュレーションモデルを示す概略図、
図４は、回折係数を計算するための座標系を示す概略
図、図５は、電波発射源の送信アンテナの指向性の推定
を説明する図、図６は、電波監視方法の手順を示すフロ
ーチャートである。
視方法及び装置について図１乃至図６を用いて説明す
る。図１は、本実施形態による電波監視方法を示す概略
図、図２は、電波監視装置の構成を示すブロック図、図
３は、電波伝搬シミュレーションモデルを示す概略図、
図４は、回折係数を計算するための座標系を示す概略
図、図５は、電波発射源の送信アンテナの指向性の推定
を説明する図、図６は、電波監視方法の手順を示すフロ
ーチャートである。
【００２５】（電波監視装置）まず、本実施形態による
電波監視装置の概要について図１を用いて説明する。
電波監視装置の概要について図１を用いて説明する。
【００２６】電波発射源１０は特定すべき電波を発射す
る発射源であり、周波数ｆ１の電波を発射するものとす
る。このような電波発射源１０から発射された電波の電
波監視を行うために、監視区域内には、発射された電波
の電波ホログラム観測を行うセンサ局１２ａ、１２ｂ、
１２ｃが配置されている。センサ局１２ａ、１２ｂ、１
２ｃは、回線１４を介して、電波伝搬シミュレーション
を行い、観測結果と比比較し電波発射源の位置１０を特
定するセンタ局１６に接続されている。
る発射源であり、周波数ｆ１の電波を発射するものとす
る。このような電波発射源１０から発射された電波の電
波監視を行うために、監視区域内には、発射された電波
の電波ホログラム観測を行うセンサ局１２ａ、１２ｂ、
１２ｃが配置されている。センサ局１２ａ、１２ｂ、１
２ｃは、回線１４を介して、電波伝搬シミュレーション
を行い、観測結果と比比較し電波発射源の位置１０を特
定するセンタ局１６に接続されている。
【００２７】各センサ局１２ａ、１２ｂ、１２ｃには、
電波ホログラム観測を行って電波再生像を出力する電波
ホログラム観測再生装置１８が設けられている。電波ホ
ログラム観測再生装置１８は、例えば、固定アンテナと
走査観測面内を走査する走査アンテナとを備えている。
所定の観測周波数において固定アンテナでの受信信号と
走査アンテナでの受信信号を相関させることによって走
査観測面の各点での相関値（２次元複素インタフェログ
ラム）を求め、この２次元複素インタフェログラムを再
生することによって、電波再生像を得る。
電波ホログラム観測を行って電波再生像を出力する電波
ホログラム観測再生装置１８が設けられている。電波ホ
ログラム観測再生装置１８は、例えば、固定アンテナと
走査観測面内を走査する走査アンテナとを備えている。
所定の観測周波数において固定アンテナでの受信信号と
走査アンテナでの受信信号を相関させることによって走
査観測面の各点での相関値（２次元複素インタフェログ
ラム）を求め、この２次元複素インタフェログラムを再
生することによって、電波再生像を得る。
【００２８】各センサ局１２ａ、１２ｂ、１２ｃには、
電波発射源１０からの電波が伝搬するが、山岳やその他
の地物の影響により、電波発射源１０からセンサ局へ直
接伝搬する成分（直接波）の他、回折して伝搬する成分
（回折波）や反射して伝搬する成分（反射波）が存在す
る。そこで、本実施形態では、各センサ局１２ａ、１２
ｂ、１２ｃにおいて、伝搬経路ごとの成分を分離してそ
の到来方向や強度を求めるために、電波ホログラム観測
を行い、電波像を再生することとしている。電波も波動
の一種であるから、光のホログラムの場合と同様にホロ
グラム観測を行うことができ、ホログラムを再生するこ
とによって電波再生像が得られ、この電波再生像から波
源分布や強度などを調べることができる。なお、伝搬経
路ごとの成分を分離して観測できるのであれば、電波ホ
ログラム観測によらずに従来の方向探知技術を用いても
よい。
電波発射源１０からの電波が伝搬するが、山岳やその他
の地物の影響により、電波発射源１０からセンサ局へ直
接伝搬する成分（直接波）の他、回折して伝搬する成分
（回折波）や反射して伝搬する成分（反射波）が存在す
る。そこで、本実施形態では、各センサ局１２ａ、１２
ｂ、１２ｃにおいて、伝搬経路ごとの成分を分離してそ
の到来方向や強度を求めるために、電波ホログラム観測
を行い、電波像を再生することとしている。電波も波動
の一種であるから、光のホログラムの場合と同様にホロ
グラム観測を行うことができ、ホログラムを再生するこ
とによって電波再生像が得られ、この電波再生像から波
源分布や強度などを調べることができる。なお、伝搬経
路ごとの成分を分離して観測できるのであれば、電波ホ
ログラム観測によらずに従来の方向探知技術を用いても
よい。
【００２９】センタ局１６には、電波伝搬を計算機シミ
ュレーションするためのシミュレーション装置２０と、
シミュレーション装置２０によって予めシミュレーショ
ンを行った結果を格納しておくデータベース２６と、各
センサ局１２ａ、１２ｂ、１２ｃで観測された電波再生
像とシミュレーションによる電界強度分布とを比較する
比較部２８とが設けられている。また、比較部２７によ
る比較結果から電波発射源１０の位置を特定するととも
に方向探知トレース図を作成、出力する方向探知トレー
ス図作成部３０と、位置が特定された電波発射源１０か
らの電波の電界強度分布を電波勢力地図として作成、出
力する電波勢力地図作成部３２が設けられている。
ュレーションするためのシミュレーション装置２０と、
シミュレーション装置２０によって予めシミュレーショ
ンを行った結果を格納しておくデータベース２６と、各
センサ局１２ａ、１２ｂ、１２ｃで観測された電波再生
像とシミュレーションによる電界強度分布とを比較する
比較部２８とが設けられている。また、比較部２７によ
る比較結果から電波発射源１０の位置を特定するととも
に方向探知トレース図を作成、出力する方向探知トレー
ス図作成部３０と、位置が特定された電波発射源１０か
らの電波の電界強度分布を電波勢力地図として作成、出
力する電波勢力地図作成部３２が設けられている。
【００３０】シミュレーション装置２０は、地図、地
形、地物等に関する情報（地図情報）を格納した地図情
報格納部２２と、地図情報を参照し、地図上のセンサ局
１２から仮想的な電波発射を行い電波伝搬のシミュレー
ションを行う電波伝搬シミュレータ２４とを備えてい
る。電波伝搬シミュレータ２４は、例えば、レイトレー
ス法（Ray Tracing）（例えば、IEEE Network Magazin
e, pp.27-30, November, 1991 を参照）やモーメント法
（Moment Method）（例えば、R. F. Harrington, "Fiel
d Computation by Moment Methods", IEEE Press, 1993
を参照）などを用い、地図データ上のセンサ局１２か
ら仮想的に電波を発射したときに生じる電界強度分布を
計算機シミュレーションする。
形、地物等に関する情報（地図情報）を格納した地図情
報格納部２２と、地図情報を参照し、地図上のセンサ局
１２から仮想的な電波発射を行い電波伝搬のシミュレー
ションを行う電波伝搬シミュレータ２４とを備えてい
る。電波伝搬シミュレータ２４は、例えば、レイトレー
ス法（Ray Tracing）（例えば、IEEE Network Magazin
e, pp.27-30, November, 1991 を参照）やモーメント法
（Moment Method）（例えば、R. F. Harrington, "Fiel
d Computation by Moment Methods", IEEE Press, 1993
を参照）などを用い、地図データ上のセンサ局１２か
ら仮想的に電波を発射したときに生じる電界強度分布を
計算機シミュレーションする。
【００３１】データベース２６は、シミュレーション装
置２０によって計算された電界強度分布の行列を格納す
る。データベース２６に格納された電界強度分布の行列
は、電波発射源１０の位置の特定の際に、比較部２８に
よって参照される。
置２０によって計算された電界強度分布の行列を格納す
る。データベース２６に格納された電界強度分布の行列
は、電波発射源１０の位置の特定の際に、比較部２８に
よって参照される。
【００３２】比較部２８は、シミュレーション装置２０
によって計算されデータベース２６に格納されている、
センサ局１２ａから仮想的な電波を発射したときの電界
強度分布の計算機シミュレーションの結果と、実際のセ
ンサ局１２ａでの電波ホログラム観測の結果とを比較す
る。この比較結果より、電波発射源１０の位置を特定す
る。
によって計算されデータベース２６に格納されている、
センサ局１２ａから仮想的な電波を発射したときの電界
強度分布の計算機シミュレーションの結果と、実際のセ
ンサ局１２ａでの電波ホログラム観測の結果とを比較す
る。この比較結果より、電波発射源１０の位置を特定す
る。
【００３３】方向探知トレース図作成部３０は、比較部
２８によって特定された電波発射源１０位置から各セン
サ局１２ａ、１２ｂ、１２ｃまでの電波伝搬路のトレー
スを行って方向探知トレース図３４を作成、出力し、さ
らに、トレース結果をもとに、電波発射源１０でのアン
テナ指向性を推定する。
２８によって特定された電波発射源１０位置から各セン
サ局１２ａ、１２ｂ、１２ｃまでの電波伝搬路のトレー
スを行って方向探知トレース図３４を作成、出力し、さ
らに、トレース結果をもとに、電波発射源１０でのアン
テナ指向性を推定する。
【００３４】電波勢力地図作成部３２は、比較部２８に
よって特定された電波発射源１０の位置と、方向探知ト
レース図作成部３０による電波発射源１０のアンテナ指
向性の推定結果とに基づき、シミュレーション装置２０
に電波伝搬の計算機シミュレーションを実行させ、地図
上でのその電波発射源１０からの電波の電界強度分布を
計算させる。この計算結果を、その電波発射源１０の電
波勢力地図３６として表示、印刷する。
よって特定された電波発射源１０の位置と、方向探知ト
レース図作成部３０による電波発射源１０のアンテナ指
向性の推定結果とに基づき、シミュレーション装置２０
に電波伝搬の計算機シミュレーションを実行させ、地図
上でのその電波発射源１０からの電波の電界強度分布を
計算させる。この計算結果を、その電波発射源１０の電
波勢力地図３６として表示、印刷する。
【００３５】本実施形態による電波監視装置は、上記の
構成をもとに、地形データや地物データに基づき、セン
サ局１２ａ側から仮想的な電波発射をしたときの電波発
射方向に対する各エリアの電界強度分布を、予め計算機
シミュレーションによって算出する。そして、実際にセ
ンサ局１２ａで観測された電波再生像から得られる電波
の到来方向及びそれらの強度と、計算機シミュレーショ
ンによる電界強度分布とを比較することにより、電波発
射源１０の位置を特定することに特徴がある。
構成をもとに、地形データや地物データに基づき、セン
サ局１２ａ側から仮想的な電波発射をしたときの電波発
射方向に対する各エリアの電界強度分布を、予め計算機
シミュレーションによって算出する。そして、実際にセ
ンサ局１２ａで観測された電波再生像から得られる電波
の到来方向及びそれらの強度と、計算機シミュレーショ
ンによる電界強度分布とを比較することにより、電波発
射源１０の位置を特定することに特徴がある。
【００３６】（電波監視方法）次に、本実施形態による
電波監視方法について図１乃至図６を用いて詳述する電
波監視に先立ち、地図データ上のセンサ局１２ａから仮
想的な電波発射を行ったときの２次元モデルの電界強度
分布のシミュレーションを行う。以下に、電界強度分布
の計算機シミュレーションについて図３及び図４を用い
て詳述する。以下に述べる計算機シミュレーションはレ
イトレース法を用いたものである。レイトレース法と
は、電波源から無数のレイ（光線）が発射されているも
のとして、そのレイを一定角度ごとに逐次奇跡の追跡を
行う手法である。従って、そのレイが地面や建物にぶつ
かった場合にはその点を反射点として反射方向へのレイ
を更に追跡し、建物のエッジ部分にぶつかった場合はそ
の点を回折点としてそこから再び複数のレイを発生させ
る方法である。
電波監視方法について図１乃至図６を用いて詳述する電
波監視に先立ち、地図データ上のセンサ局１２ａから仮
想的な電波発射を行ったときの２次元モデルの電界強度
分布のシミュレーションを行う。以下に、電界強度分布
の計算機シミュレーションについて図３及び図４を用い
て詳述する。以下に述べる計算機シミュレーションはレ
イトレース法を用いたものである。レイトレース法と
は、電波源から無数のレイ（光線）が発射されているも
のとして、そのレイを一定角度ごとに逐次奇跡の追跡を
行う手法である。従って、そのレイが地面や建物にぶつ
かった場合にはその点を反射点として反射方向へのレイ
を更に追跡し、建物のエッジ部分にぶつかった場合はそ
の点を回折点としてそこから再び複数のレイを発生させ
る方法である。
【００３７】まず、電波伝搬シミュレータ２４によっ
て、地図情報格納部２２から、シミュレーションを行う
地域の地図データを参照し、図３に示すように、シミュ
レーションの対象とする地域を格子状にＭ×Ｎ個のエリ
アに分割する。分割した各エリアにはエリア番号１、
２、…、ＭＮ−１、ＭＮを付与する。地図データには、
地形や地物に関する情報がマッピングされている。図３
では、建物１３ａ、建物１３ｂがある場合を考える。
て、地図情報格納部２２から、シミュレーションを行う
地域の地図データを参照し、図３に示すように、シミュ
レーションの対象とする地域を格子状にＭ×Ｎ個のエリ
アに分割する。分割した各エリアにはエリア番号１、
２、…、ＭＮ−１、ＭＮを付与する。地図データには、
地形や地物に関する情報がマッピングされている。図３
では、建物１３ａ、建物１３ｂがある場合を考える。
【００３８】続いて、図３に示すように、地図データ上
のセンサ局１２ａからＬk方向に一本のレイ（ Ray ）を
発射し、このレイの伝搬経路と各エリアを通過したとき
の電界強度を計算する。計算されたレイの電界強度を、
行番号に発射したレイのレイ番号Ｌkを有し、列番号に
エリア番号を有する行列の要素として記録する。この電
界強度の分布を記録する行列を次の表１に示す。
のセンサ局１２ａからＬk方向に一本のレイ（ Ray ）を
発射し、このレイの伝搬経路と各エリアを通過したとき
の電界強度を計算する。計算されたレイの電界強度を、
行番号に発射したレイのレイ番号Ｌkを有し、列番号に
エリア番号を有する行列の要素として記録する。この電
界強度の分布を記録する行列を次の表１に示す。
【００３９】
【表１】
【００４０】このとき、センサ局１２ａを起点として、
レイの移動距離に比例した太さをレイにもたせるものと
し、太さを有するレイが各エリアの一部にかかった場合
も、そのレイがそのエリアを通過したものとみなす。ま
た、一つのエリアに複数のレイが通過した場合は、その
うちの最大強度を与えるレイの電界強度を表１の行列に
記録するものとする。
レイの移動距離に比例した太さをレイにもたせるものと
し、太さを有するレイが各エリアの一部にかかった場合
も、そのレイがそのエリアを通過したものとみなす。ま
た、一つのエリアに複数のレイが通過した場合は、その
うちの最大強度を与えるレイの電界強度を表１の行列に
記録するものとする。
【００４１】上記のシミュレーションの手順を、図３に
示すように、センサ局１２ａを起点としてＬk方向に発
射されたレイが、建物１３ａのｐ点で反射し、建物１３
ｂのｑ点に到達する場合を用いて具体的に説明する。こ
こで、センサ局とｐ点との間を結ぶ線分を線分ａとし、
ｐ点とｑ点との間を結ぶ線分を線分ｂとする。
示すように、センサ局１２ａを起点としてＬk方向に発
射されたレイが、建物１３ａのｐ点で反射し、建物１３
ｂのｑ点に到達する場合を用いて具体的に説明する。こ
こで、センサ局とｐ点との間を結ぶ線分を線分ａとし、
ｐ点とｑ点との間を結ぶ線分を線分ｂとする。
【００４２】まず、レイがセンサ局１２ａから発射さ
れ、レイが建物１３ａのｐ点に到達する間の線分ａが通
過したエリアについて、対応する行列の要素（Ｌk，エ
リア番号）に次の数１により計算される電界強度を記録
する。
れ、レイが建物１３ａのｐ点に到達する間の線分ａが通
過したエリアについて、対応する行列の要素（Ｌk，エ
リア番号）に次の数１により計算される電界強度を記録
する。
【００４３】
【数１】
【００４４】ただし、Ｅ0はセンサ局１２ａ発射時の電
波の電界強度、ｒはセンサ局１２ａから各エリアの中心
までの距離、ｋは電波の波数である。
波の電界強度、ｒはセンサ局１２ａから各エリアの中心
までの距離、ｋは電波の波数である。
【００４５】続いて、レイが建物１３ａのｐ点で反射
し、反射波が建物１３ｂのｑ点に到達する間の線分ｂが
通過したエリアについて、対応する行列の要素（Ｌk，
エリア番号）に次の数２により計算される電界強度を記
録する。
し、反射波が建物１３ｂのｑ点に到達する間の線分ｂが
通過したエリアについて、対応する行列の要素（Ｌk，
エリア番号）に次の数２により計算される電界強度を記
録する。
【００４６】
【数２】
【００４７】ただし、ａは線分ａの長さであり、ｒ′は
ｐ点を基点とした各エリアの中心までの距離である。
ｐ点を基点とした各エリアの中心までの距離である。
【００４８】Ｒは反射係数であり、以下のようにして求
められる。
められる。
【００４９】まず、スネルの法則より媒質１から媒質２
に平面波が入射するときの複素屈折率ｎは次の数３で表
される。
に平面波が入射するときの複素屈折率ｎは次の数３で表
される。
【００５０】
【数３】
【００５１】ここで、ε1、μ1、σ1は、それぞれ媒質
１の誘電率、透磁率、導電率である。ε2、μ2、σ
2は、それぞれ媒質２の誘電率、透磁率、導電率であ
る。求めようとする反射係数Ｒは、複素屈折率ｎを用い
て、入射電界が入射面内にあるとき、すなわちＴＭ（Tr
ansverse Magnetic）入射の場合は次の数４で表され
る。
１の誘電率、透磁率、導電率である。ε2、μ2、σ
2は、それぞれ媒質２の誘電率、透磁率、導電率であ
る。求めようとする反射係数Ｒは、複素屈折率ｎを用い
て、入射電界が入射面内にあるとき、すなわちＴＭ（Tr
ansverse Magnetic）入射の場合は次の数４で表され
る。
【００５２】
【数４】
【００５３】一方、入射電界が入射面に垂直であると
き、すなわちＴＥ（Transverse Electric）入射の場合
は次の数５で表される。
き、すなわちＴＥ（Transverse Electric）入射の場合
は次の数５で表される。
【００５４】
【数５】
【００５５】続いて、建物１３ａのｐ点で反射したレイ
が建物１３ｂのｑ点に到達し、ｑ点で回折したときの回
折波の線分ｃ1、ｃ2、ｃ3、…を求める。ここで、レイ
が建物１３ｂの角に入射した場合、すなわちｑ点が建物
１３ｂの角に位置する場合には、その角から複数のレイ
を面放射するものとする。
が建物１３ｂのｑ点に到達し、ｑ点で回折したときの回
折波の線分ｃ1、ｃ2、ｃ3、…を求める。ここで、レイ
が建物１３ｂの角に入射した場合、すなわちｑ点が建物
１３ｂの角に位置する場合には、その角から複数のレイ
を面放射するものとする。
【００５６】線分ｃ1、ｃ2、ｃ3、…が通過したエリア
について、対応する行列の要素（Ｌk，エリア番号）に
次の数６を記録する。
について、対応する行列の要素（Ｌk，エリア番号）に
次の数６を記録する。
【００５７】
【数６】
【００５８】ただし、ｂは線分ｂの長さであり、ｒ″は
ｑ点を基点とした各エリアの中心までの距離、Ｄは回折
係数であり、図４に示す座標系を用いて次の数７から計
算される。
ｑ点を基点とした各エリアの中心までの距離、Ｄは回折
係数であり、図４に示す座標系を用いて次の数７から計
算される。
【００５９】
【数７】
【００６０】ただし、β0、β0′φ、φ′、ρd、ｓ
dは、それぞれ図４に示すような位置関係を表す物理的
なパラメータであり、図４でＥ0は入射波、Ｅdは回折波
である。また、Ｄs及びＤhは、ＵＴＤ法（R. G. Kouyou
mjian, P. H. Pathak, Proceedingof the IEEE, no.11,
pp.1448-1461, December, 1974 を参照）によりフレネ
ル積分関数Ｆを用いて次の数８のように表される。
dは、それぞれ図４に示すような位置関係を表す物理的
なパラメータであり、図４でＥ0は入射波、Ｅdは回折波
である。また、Ｄs及びＤhは、ＵＴＤ法（R. G. Kouyou
mjian, P. H. Pathak, Proceedingof the IEEE, no.11,
pp.1448-1461, December, 1974 を参照）によりフレネ
ル積分関数Ｆを用いて次の数８のように表される。
【００６１】
【数８】
【００６２】ここで、Ｌ、ａ+（ｘ）、ａ-（ｘ）はそれ
ぞれ次の数９のように表される。
ぞれ次の数９のように表される。
【００６３】
【数９】
【００６４】ただし、Ｎ+、Ｎ-は次の数１０で表され
る。
る。
【００６５】
【数１０】
【００６６】上述の操作を、センサ局１２ａから発射す
るレイの方向Ｌkを変化しながら繰り返し、表１に示す
Ｌkに対する電界強度分布の行列を完成する。得られた
電界強度分布の行列をデータベース２６に記録してお
く。こうして、電波伝搬のシミュレーションを終了す
る。
るレイの方向Ｌkを変化しながら繰り返し、表１に示す
Ｌkに対する電界強度分布の行列を完成する。得られた
電界強度分布の行列をデータベース２６に記録してお
く。こうして、電波伝搬のシミュレーションを終了す
る。
【００６７】図１に示すように複数のセンサ局１２ａ、
１２ｂ、１２ｃを配置した場合は、センサ局１２ａ以外
についても上述のシミュレーションを行い、複数センサ
局による同時監視を行えるようにしてもよい。センサ局
を複数用いることによって不感地域をなくすることがで
きるほか、測定データを多くすることで電波発射源の位
置の推定精度を向上させることができる。
１２ｂ、１２ｃを配置した場合は、センサ局１２ａ以外
についても上述のシミュレーションを行い、複数センサ
局による同時監視を行えるようにしてもよい。センサ局
を複数用いることによって不感地域をなくすることがで
きるほか、測定データを多くすることで電波発射源の位
置の推定精度を向上させることができる。
【００６８】以上のように、予め電波伝搬のシミュレー
ションを行ったうえで、電波監視を開始する。電波監視
は、図６に示すフローチャートの手順に従って行われ
る。
ションを行ったうえで、電波監視を開始する。電波監視
は、図６に示すフローチャートの手順に従って行われ
る。
【００６９】まず、センサ局１２ａにより電波ホログラ
ム観測を行い、電波再生像を取得する（ステップＳ１
１）。電波ホログラム観測により、電波の到来方向と強
度とを関連づけて測定することができる。
ム観測を行い、電波再生像を取得する（ステップＳ１
１）。電波ホログラム観測により、電波の到来方向と強
度とを関連づけて測定することができる。
【００７０】続いて、比較部２８によって、センサ局１
２ａで観測された電波再生像から得られる電波の到来方
向及び強度パターンと、データベース２６に格納されて
いる電界強度分布の行列の列ベクトルそれぞれのパター
ンとの比較を行う（ステップＳ１２）。
２ａで観測された電波再生像から得られる電波の到来方
向及び強度パターンと、データベース２６に格納されて
いる電界強度分布の行列の列ベクトルそれぞれのパター
ンとの比較を行う（ステップＳ１２）。
【００７１】電波伝搬では、一般に、送信と受信につい
て電波伝搬経路の可逆性が成立し、また、伝搬減衰量の
可逆性、すなわち送信と受信を入れ換えても伝搬減衰量
は同じであるということが成立する。
て電波伝搬経路の可逆性が成立し、また、伝搬減衰量の
可逆性、すなわち送信と受信を入れ換えても伝搬減衰量
は同じであるということが成立する。
【００７２】したがって、電界強度分布の行列の列ベク
トルのうち、最も一致度の高い列ベクトルのエリア番号
を有するエリアを電波発射源１０の位置とすることがで
きる。
トルのうち、最も一致度の高い列ベクトルのエリア番号
を有するエリアを電波発射源１０の位置とすることがで
きる。
【００７３】このようにして電波発射源１０を特定する
方法を採用することにより、電波監視を行う地域内の複
数の地点の各々について、当該地点に電波発射源がある
としたときのセンサ局で観測される電波再生像を、地図
情報を利用する計算機シミュレーションによって算出す
る必要がなく、電波監視の効率を向上することができ
る。
方法を採用することにより、電波監視を行う地域内の複
数の地点の各々について、当該地点に電波発射源がある
としたときのセンサ局で観測される電波再生像を、地図
情報を利用する計算機シミュレーションによって算出す
る必要がなく、電波監視の効率を向上することができ
る。
【００７４】上述の一致度の比較では、通常のパターン
マッチングの手法を用いることができる。例えば、電波
再生像における電波の最大レベルが０ｄＢであるとすれ
ば、電波レベルが−４０ｄＢ以上の部分のパターンを比
較し、図形としての相関係数を求めるといった手法を用
いる。電波発射源１０のアンテナ指向性が分からないの
で、電波の到来方向に主に着目して一致性の高いものを
探索するようにする。
マッチングの手法を用いることができる。例えば、電波
再生像における電波の最大レベルが０ｄＢであるとすれ
ば、電波レベルが−４０ｄＢ以上の部分のパターンを比
較し、図形としての相関係数を求めるといった手法を用
いる。電波発射源１０のアンテナ指向性が分からないの
で、電波の到来方向に主に着目して一致性の高いものを
探索するようにする。
【００７５】センサ局１２ａだけでなく、センサ局１２
ｂ、センサ局１２ｃを用いて同時監視を行う場合は、セ
ンサ局１２ａの場合と同様に、予め仮想的な電波発射を
行ったときの電界強度分布のシミュレーションの結果を
データベース２６に格納しておき、実際に電波ホログラ
ム観測によって得られた結果との比較を行う。
ｂ、センサ局１２ｃを用いて同時監視を行う場合は、セ
ンサ局１２ａの場合と同様に、予め仮想的な電波発射を
行ったときの電界強度分布のシミュレーションの結果を
データベース２６に格納しておき、実際に電波ホログラ
ム観測によって得られた結果との比較を行う。
【００７６】次いで、電波発射源１０の位置を特定した
後、方向探知トレース図作成部３０により、比較部２８
によって特定された電波発射源１０位置から各センサ局
１２ａ、１２ｂ、１２ｃまでの電波伝搬路のトレースを
行って方向探知トレース図３４を作成、出力する（ステ
ップＳ１３）。
後、方向探知トレース図作成部３０により、比較部２８
によって特定された電波発射源１０位置から各センサ局
１２ａ、１２ｂ、１２ｃまでの電波伝搬路のトレースを
行って方向探知トレース図３４を作成、出力する（ステ
ップＳ１３）。
【００７７】なお、地図データ上のセンサ局１２ａから
仮想的な電波発射を行ったときの電界強度分布のシミュ
レーションの際に、得られた電界強度分布の行列に対応
してレイの経路情報をデータベース２６に記録しておく
ことにより、方向探知トレース図作成部３０を用いるこ
となく、電波発射源１０の位置の特定と同時に、電波伝
搬経路のトレース結果を得ることができる。
仮想的な電波発射を行ったときの電界強度分布のシミュ
レーションの際に、得られた電界強度分布の行列に対応
してレイの経路情報をデータベース２６に記録しておく
ことにより、方向探知トレース図作成部３０を用いるこ
となく、電波発射源１０の位置の特定と同時に、電波伝
搬経路のトレース結果を得ることができる。
【００７８】さらに、上述のトレース結果をもとに、電
波発射源１０でのアンテナ指向性を以下のようにして推
定する（ステップＳ１４）。
波発射源１０でのアンテナ指向性を以下のようにして推
定する（ステップＳ１４）。
【００７９】まず、各センサ局１２ａ、１２ｂ、１２ｃ
によって電波ホログラム観測を行い、電波再生像を取得
する。各センサ局１２ａ、１２ｂ、１２ｃで得られた電
波再生像から、伝搬経路ごとの電波の到来方向とその強
度を抽出する。
によって電波ホログラム観測を行い、電波再生像を取得
する。各センサ局１２ａ、１２ｂ、１２ｃで得られた電
波再生像から、伝搬経路ごとの電波の到来方向とその強
度を抽出する。
【００８０】図５（ａ）〜（ｃ）は、それぞれセンサ局
１２ａ、１２ｂ、１２ｃでの観測による電波の振幅を示
している。太線の矢印の向きがその伝搬経路の電波の到
来方向を示し、その矢印の長さが電波の振幅を示してい
る。図示した例では、伝搬経路ごとの振幅が、ａ1、
ａ2、ｂ1、ｂ2、ｃ1、ｃ2となっている。一方、シミュ
レータ装置２３によって、すでに求めた電波発射源１０
の位置に無指向性の送信アンテナがあってそこから電波
が発射されているものとして、電波伝搬シミュレーショ
ンを行い、伝搬経路ごとに、各センサ局１２ａ、１２
ｂ、１２ｃでの電波の振幅を求める。
１２ａ、１２ｂ、１２ｃでの観測による電波の振幅を示
している。太線の矢印の向きがその伝搬経路の電波の到
来方向を示し、その矢印の長さが電波の振幅を示してい
る。図示した例では、伝搬経路ごとの振幅が、ａ1、
ａ2、ｂ1、ｂ2、ｃ1、ｃ2となっている。一方、シミュ
レータ装置２３によって、すでに求めた電波発射源１０
の位置に無指向性の送信アンテナがあってそこから電波
が発射されているものとして、電波伝搬シミュレーショ
ンを行い、伝搬経路ごとに、各センサ局１２ａ、１２
ｂ、１２ｃでの電波の振幅を求める。
【００８１】図５（ｄ）〜（ｆ）は、各センサ局１２
ａ、１２ｂ、１２ｃでのシミュレーションによる電波の
振幅を示しており、伝搬経路ごとの振幅は、ａ1’、
ａ2’、ｂ 1’、ｂ2’、ｃ1’、ｃ2’となっている。そ
して、伝搬経路Ｐ1、Ｐ2、Ｐ3、Ｐ4、Ｐ5、Ｐ6ごとに観
測による振幅をシミュレーションによる振幅で除算する
ことによって、電波発射源の位置での各伝搬経路の方向
への、送信アンテナの指向性を求めることができる。
ａ、１２ｂ、１２ｃでのシミュレーションによる電波の
振幅を示しており、伝搬経路ごとの振幅は、ａ1’、
ａ2’、ｂ 1’、ｂ2’、ｃ1’、ｃ2’となっている。そ
して、伝搬経路Ｐ1、Ｐ2、Ｐ3、Ｐ4、Ｐ5、Ｐ6ごとに観
測による振幅をシミュレーションによる振幅で除算する
ことによって、電波発射源の位置での各伝搬経路の方向
への、送信アンテナの指向性を求めることができる。
【００８２】図５（ｇ）において黒丸は各伝搬経路の方
向への指向性を示している。そして、想定されるアンテ
ナの種類などを考慮して補間計算を行えば、任意の方向
に対する指向特性を算出することができる。図５（ｇ）
の破線はこのようにして算出した指向特性を示してい
る。
向への指向性を示している。そして、想定されるアンテ
ナの種類などを考慮して補間計算を行えば、任意の方向
に対する指向特性を算出することができる。図５（ｇ）
の破線はこのようにして算出した指向特性を示してい
る。
【００８３】また、電波勢力地図作成部３２によって、
方向探知トレース図作成部３０による電波発射源１０の
アンテナ指向性の推定結果とともに、シミュレーション
装置２０に電波伝搬の計算機シミュレーションを実行さ
せ、地図上でのその電波発射源１０からの電波の電界強
度分布を計算させる。この計算結果を、その電波発射源
１０の電波勢力地図３６として表示、印刷する（ステッ
プＳ１５）。
方向探知トレース図作成部３０による電波発射源１０の
アンテナ指向性の推定結果とともに、シミュレーション
装置２０に電波伝搬の計算機シミュレーションを実行さ
せ、地図上でのその電波発射源１０からの電波の電界強
度分布を計算させる。この計算結果を、その電波発射源
１０の電波勢力地図３６として表示、印刷する（ステッ
プＳ１５）。
【００８４】こうして特定地域内の電波発射源１０の位
置を特定し、同時に、方向探知トレース図の作成、電波
発射源１０の送信アンテナの指向性の推定、電波勢力地
図の作成を行い、電波監視を終了する。
置を特定し、同時に、方向探知トレース図の作成、電波
発射源１０の送信アンテナの指向性の推定、電波勢力地
図の作成を行い、電波監視を終了する。
【００８５】このように、本実施形態によれば、地域内
の複数の地点の各々について、当該地点に電波発射源が
あるとしたときのセンサ局で観測される電波再生像を計
算機シミュレーションによって算出する必要がないの
で、電波伝搬シミュレーションによるデータベースの作
成時間を短縮し、電波監視の効率を向上することができ
る。
の複数の地点の各々について、当該地点に電波発射源が
あるとしたときのセンサ局で観測される電波再生像を計
算機シミュレーションによって算出する必要がないの
で、電波伝搬シミュレーションによるデータベースの作
成時間を短縮し、電波監視の効率を向上することができ
る。
【００８６】次に、本発明の他の実施形態による電波監
視方法及び装置について図７を用いて説明する。図７
は、３次元の電波伝搬シミュレーションモデルを示す概
略図である。
視方法及び装置について図７を用いて説明する。図７
は、３次元の電波伝搬シミュレーションモデルを示す概
略図である。
【００８７】本実施形態による電波監視方法及び装置
は、地図データ上のセンサ局１２ａから仮想的な電波発
射を行ったときの電界強度分布のシミュレーションを３
次元モデルを用いて行うことに特徴がある。電界強度分
布のシミュレーションを３次元モデルを用いて行う点を
除いては、上記実施形態と同様である。以下に本実施形
態における電界強度分布のシミュレーションについて説
明する。
は、地図データ上のセンサ局１２ａから仮想的な電波発
射を行ったときの電界強度分布のシミュレーションを３
次元モデルを用いて行うことに特徴がある。電界強度分
布のシミュレーションを３次元モデルを用いて行う点を
除いては、上記実施形態と同様である。以下に本実施形
態における電界強度分布のシミュレーションについて説
明する。
【００８８】まず、図７に示すように、監視区域をｘ×
ｙ×ｚ個の単位空間に分割し、各単位空間に空間番号
１、２、３、…、ｘｙｚ−１、ｘｙｚを付与する。
ｙ×ｚ個の単位空間に分割し、各単位空間に空間番号
１、２、３、…、ｘｙｚ−１、ｘｙｚを付与する。
【００８９】続いて、地図データ上のセンサ局１２ａか
らＬk方向に一本のレイを発射し、このレイの伝搬経路
と各単位空間を通過したときの電界強度を、上記実施形
態における２次元モデルでの場合と同様に計算する。計
算されたレイの電界強度を、行番号に発射したレイのレ
イ番号Ｌkを有し、列番号に単位空間番号を有する行列
の要素として記録する。
らＬk方向に一本のレイを発射し、このレイの伝搬経路
と各単位空間を通過したときの電界強度を、上記実施形
態における２次元モデルでの場合と同様に計算する。計
算されたレイの電界強度を、行番号に発射したレイのレ
イ番号Ｌkを有し、列番号に単位空間番号を有する行列
の要素として記録する。
【００９０】このとき、センサ局１２ａを起点として、
レイの移動距離に比例した太さをレイにもたせるものと
し、太さを有するレイが各エリアの一部にかかった場合
も、そのレイがそのエリアを通過したものとみなす。ま
た、一つのエリアに複数のレイが通過した場合は、その
うちの最大強度を与えるレイの電界強度を上述の行列に
記録する。
レイの移動距離に比例した太さをレイにもたせるものと
し、太さを有するレイが各エリアの一部にかかった場合
も、そのレイがそのエリアを通過したものとみなす。ま
た、一つのエリアに複数のレイが通過した場合は、その
うちの最大強度を与えるレイの電界強度を上述の行列に
記録する。
【００９１】センサ局１２ａを起点としてＬk方向に発
射されたレイが、建物１３ａのｐ点で反射する場合や、
建物１３ｂのｑ点で回折する場合についても、上記実施
形態と同様にレイの電界強度を計算する。この際、建物
１３ａ、１３ｂの３次元情報を考慮して反射や回折の方
向を決定する。
射されたレイが、建物１３ａのｐ点で反射する場合や、
建物１３ｂのｑ点で回折する場合についても、上記実施
形態と同様にレイの電界強度を計算する。この際、建物
１３ａ、１３ｂの３次元情報を考慮して反射や回折の方
向を決定する。
【００９２】上述の操作を、センサ局１２ａから発射す
るレイの方向Ｌkを３次元的に変化しながら繰り返し、
２次元モデルでの場合と同様に、Ｌkに対する電界強度
分布の行列を完成する。得られた電界強度分布の行列を
データベース２６に記録しておく。こうして、電波伝搬
のシミュレーションを終了する。
るレイの方向Ｌkを３次元的に変化しながら繰り返し、
２次元モデルでの場合と同様に、Ｌkに対する電界強度
分布の行列を完成する。得られた電界強度分布の行列を
データベース２６に記録しておく。こうして、電波伝搬
のシミュレーションを終了する。
【００９３】上述の電波伝搬のシミュレーションの結果
を用いて、以後、図６に示すフローチャートに従って、
上記実施形態と同様に電波監視を行うことができる。
を用いて、以後、図６に示すフローチャートに従って、
上記実施形態と同様に電波監視を行うことができる。
【００９４】このように、本実施形態によれば、地域内
の複数の地点の各々について、当該地点に電波発射源が
あるとしたときのセンサ局で観測される電波再生像を計
算機シミュレーションによって算出する必要がないの
で、電波伝搬シミュレーションによるデータベースの作
成時間を短縮し、電波監視の効率を向上することができ
る。
の複数の地点の各々について、当該地点に電波発射源が
あるとしたときのセンサ局で観測される電波再生像を計
算機シミュレーションによって算出する必要がないの
で、電波伝搬シミュレーションによるデータベースの作
成時間を短縮し、電波監視の効率を向上することができ
る。
【００９５】
【発明の効果】以上の通り、本発明によれば、一の地点
における観測結果を、一の地点から電波を発射したとき
に他の複数の地点において観測される電波の到来方向及
び強度のパターンのシミュレーション結果と比較し、複
数の地点のうち一の地点における観測結果と最も相関の
強い電波の到来方向及び強度のパターンを与える地点を
電波発射源と特定するので、電波伝搬シミュレーション
によるデータベースの作成時間を短縮し、電波監視の効
率を向上することができる。
における観測結果を、一の地点から電波を発射したとき
に他の複数の地点において観測される電波の到来方向及
び強度のパターンのシミュレーション結果と比較し、複
数の地点のうち一の地点における観測結果と最も相関の
強い電波の到来方向及び強度のパターンを与える地点を
電波発射源と特定するので、電波伝搬シミュレーション
によるデータベースの作成時間を短縮し、電波監視の効
率を向上することができる。
【図１】本発明の一実施形態による電波監視方法を示す
概略図である。
概略図である。
【図２】本発明の一実施形態による電波監視装置の構成
を示すブロック図である。
を示すブロック図である。
【図３】電波伝搬シミュレーションモデルを示す概略図
である。
である。
【図４】回折係数を求めるための座標系を示す概略図で
ある。
ある。
【図５】電波発射源の送信アンテナの指向性の推定を説
明する図である。
明する図である。
【図６】本発明の一実施形態による電波監視方法の手順
を示すフローチャートである。
を示すフローチャートである。
【図７】３次元の電波伝搬シミュレーションモデルを示
す概略図である。
す概略図である。
１０…電波発射源 １２ａ、１２ｂ、１２ｃ…センサ局 １３ａ、１３ｂ…建物 １４…回線 １６…センタ局 １８…電波ホログラム観測再生装置 ２０…シミュレーション装置 ２２…地図情報格納部 ２４…電波伝搬シミュレータ ２６…データベース ２８…比較部 ３０…方向探知トレース図作成部 ３２…電波勢力地図作成部 ３４…方向探知トレース図 ３６…電波勢力地図 Ｓ１１、Ｓ１２、Ｓ１３、Ｓ１４、Ｓ１５…ステップ
Claims (15)
- 【請求項１】 一の地点において電波の到来方向及び強
度のパターンを観測し、 前記一の地点における観測結果を、前記一の地点から電
波を発射したときに他の複数の地点において観測される
電波の到来方向及び強度のパターンのシミュレーション
結果と比較し、 前記複数の地点のうち前記一の地点における前記観測結
果と最も相関の強い電波の到来方向及び強度のパターン
を与える地点を電波発射源と特定することを特徴とする
電波監視方法。 - 【請求項２】 請求項１記載の電波監視方法において、 前記一の地点から電波を発射したときの電波の到来方向
及び強度のパターンをシミュレーションする際には、観
測地域を２次元的に複数の領域に分割し、分割した各領
域において観測される電界強度を計算することを特徴と
する電波監視方法。 - 【請求項３】 請求項１記載の電波監視方法において、 前記一の地点から電波を発射したときの電波の到来方向
及び強度のパターンをシミュレーションする際には、観
測地域を３次元的に複数の空間に分割し、分割した各空
間において観測される電界強度を計算することを特徴と
する電波監視方法。 - 【請求項４】 請求項２又は３記載の電波監視方法にお
いて、 前記一の地点から電波を発射したときの電波の到来方向
及び強度のパターンをシミュレーションする際には、前
記一の地点から発射される電波の発射方向を変えて前記
各領域又は前記各空間において観測される電界強度をそ
れぞれ計算し、電波の発射方向毎に電界強度分布を求め
ることを特徴とする電波監視方法。 - 【請求項５】 請求項２乃至４のいずれか１項に記載の
電波監視方法において、 前記一の地点から電波を発射したときの電波の到来方向
及び強度のパターンをシミュレーションする際には、前
記観測地域の地形や地物を考慮して電波の到来経路を計
算することを特徴とする電波監視方法。 - 【請求項６】 請求項１乃至５のいずれか１項に記載の
電波監視方法において、 特定された前記電波発射源の位置と前記シミュレーショ
ン結果に基づいて、前記一の地点から前記電波発射源ま
での電波の伝搬経路をトレースすることを特徴とする電
波監視方法。 - 【請求項７】 請求項６記載の電波監視方法において、 前記伝搬経路のトレース結果に基づき、前記電波発射源
のアンテナ指向性を推定することを特徴とする電波監視
方法。 - 【請求項８】 請求項７記載の電波監視方法において、 推定された前記電波発射源のアンテナ指向性に基づき、
前記電波発射源から発せられた電波の電界強度分布を求
めることを特徴とする電波監視方法。 - 【請求項９】 一の地点に設けられ、電波の到来方向及
び強度のパターンを観測する電波観測手段と、 前記一の地点から電波を発射したときに他の複数の地点
において観測される電波の到来方向及び強度のパターン
のシミュレーション結果を記憶する記憶手段と、 前記電波観測手段における観測結果と、前記記憶手段に
記憶されているシミュレーション結果とを比較し、前記
複数の地点のうち前記電波観測手段における前記観測結
果と最も相関の強い電波の到来方向及び強度のパターン
を与える地点を電波発射源と特定する電波発射源特定手
段とを有することを特徴とする電波監視装置。 - 【請求項１０】 請求項９記載の電波監視装置におい
て、 前記記憶手段は、観測地域を２次元的に複数の領域に分
割し、前記一の地点から電波を発射したときに、分割し
た各領域において観測される電界強度を計算することを
特徴とする電波監視装置。 - 【請求項１１】 請求項９記載の電波監視装置におい
て、 前記記憶手段は、観測地域を３次元的に複数の空間に分
割し、前記一の地点から電波を発射したときに、分割し
た各空間において観測される電界強度を計算することを
特徴とする電波監視装置。 - 【請求項１２】 請求項１０又は１１記載の電波監視装
置において、 前記記憶手段は、前記一の地点から発射される電波の発
射方向を変えたときに前記各領域又は前記各空間におい
て観測される電界強度を、電波の発射方向毎に記憶する
ことを特徴とする電波監視装置。 - 【請求項１３】 請求項９乃至１２のいずれか１項に記
載の電波監視装置において、 前記電波発射源特定手段により特定された前記電波発射
源の位置と前記シミュレーション結果に基づいて前記一
の地点から前記電波発射源までの電波の伝搬経路をトレ
ースする伝搬経路トレース手段を更に有することを特徴
とする電波監視装置。 - 【請求項１４】 請求項１３記載の電波監視装置におい
て、 前記伝播経路トレース手段による前記伝搬経路のトレー
ス結果に基づいて前記電波発射源のアンテナ指向性を推
定するアンテナ指向性推定手段を更に有することを特徴
とする電波監視装置。 - 【請求項１５】 請求項１４記載の電波監視装置におい
て、 前記アンテナ指向性推定手段により推定された前記電波
発射源のアンテナ指向性に基づいて前記電波発射源から
発せられた電波の電界強度分布を求める電界強度分布算
出手段を更に有することを特徴とする電波監視装置。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000250028A JP3648140B2 (ja) | 2000-08-21 | 2000-08-21 | 電波監視方法及び装置 |
| US09/924,235 US7187324B2 (en) | 2000-08-21 | 2001-08-08 | Radiowave monitoring method and apparatus |
| DE10140813A DE10140813A1 (de) | 2000-08-21 | 2001-08-20 | Radiowellenüberwachungsverfahren und -Anordnung |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000250028A JP3648140B2 (ja) | 2000-08-21 | 2000-08-21 | 電波監視方法及び装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002064440A true JP2002064440A (ja) | 2002-02-28 |
| JP3648140B2 JP3648140B2 (ja) | 2005-05-18 |
Family
ID=18739676
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000250028A Expired - Fee Related JP3648140B2 (ja) | 2000-08-21 | 2000-08-21 | 電波監視方法及び装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7187324B2 (ja) |
| JP (1) | JP3648140B2 (ja) |
| DE (1) | DE10140813A1 (ja) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004531736A (ja) * | 2001-06-25 | 2004-10-14 | ハリス コーポレイション | 受動的な反射器又は屈折器を代理受信機として使用するとともに、データベースへの問合せを利用した送信機の位置を測定するためのシステム及び方法 |
| JP2005207847A (ja) * | 2004-01-22 | 2005-08-04 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波ホログラフィ法とｍｕｓｉｃ法を用い電波源の特定に用いられる電波監視装置 |
| JP2006515727A (ja) * | 2002-10-22 | 2006-06-01 | クゥアルコム・インコーポレイテッド | 無線通信システムにおいて予測電力を用いて送信機を識別する方法および装置 |
| JP2006519389A (ja) * | 2003-03-04 | 2006-08-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | オブジェクトロケーション |
| WO2007058301A1 (ja) * | 2005-11-21 | 2007-05-24 | Nec Corporation | 検知対象の識別・位置推定システムと、その方法及びプログラム |
| JP2009111952A (ja) * | 2007-11-01 | 2009-05-21 | Clarion Co Ltd | 放送信号受信機試験装置 |
| JP2010525300A (ja) * | 2007-03-13 | 2010-07-22 | アルカテル−ルーセント ユーエスエー インコーポレーテッド | 逆方向レイ追跡法を用いて送信機を位置特定する方法 |
| JP2012118024A (ja) * | 2010-12-03 | 2012-06-21 | Brother Ind Ltd | 電波推定方法、及び電波推定プログラム、並びに電波推定装置 |
| JP2012244590A (ja) * | 2011-05-24 | 2012-12-10 | Nippon Telegr & Teleph Corp <Ntt> | 到来方向推定装置、到来方向推定方法、及び到来方向推定プログラム |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102006014546B4 (de) * | 2006-03-21 | 2010-08-19 | Universität Tübingen | Verfahren und Vorrichtung zum sensorbasierten Überwachen einer Umgebung |
| WO2008011712A1 (en) * | 2006-07-28 | 2008-01-31 | Michael Tin Yau Chan | Wide-area wireless network topology |
| DE102008002892B3 (de) * | 2008-06-18 | 2009-12-31 | Forschungszentrum Karlsruhe Gmbh | Verfahren zur Bestimmung von Eigenschaften des Erdreichs mit Hilfe von freien Radiowellen |
| CN117580081B (zh) * | 2024-01-15 | 2024-03-19 | 成都大公博创信息技术有限公司 | 一种重大活动保障指挥控制系统 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3256085B2 (ja) * | 1994-06-21 | 2002-02-12 | 株式会社日立製作所 | 電波受信強度シミュレーション方法 |
| GB2311697B (en) * | 1996-03-22 | 1999-07-28 | Matsushita Electric Ind Co Ltd | Wireless communication system and method and system for detection of position of radio mobile station |
| US6249680B1 (en) * | 1997-01-08 | 2001-06-19 | U.S. Wireless Corporation | Radio transmitter location finding in CDMA wireless communication systems |
| US6026304A (en) * | 1997-01-08 | 2000-02-15 | U.S. Wireless Corporation | Radio transmitter location finding for wireless communication network services and management |
| JP3718342B2 (ja) | 1998-05-12 | 2005-11-24 | 株式会社アドバンテスト | 広域電波監視方法及び装置 |
| US6496701B1 (en) * | 1998-08-25 | 2002-12-17 | Lucent Technologies Inc. | Pattern-recognition-based geolocation |
| US6393294B1 (en) * | 1998-09-22 | 2002-05-21 | Polaris Wireless, Inc. | Location determination using RF fingerprinting |
| FI107366B (fi) * | 1998-10-29 | 2001-07-13 | Nokia Networks Oy | Menetelmä ja laitteisto verkkosuunnittelun toteuttamiseksi |
| US6487417B1 (en) * | 1999-11-24 | 2002-11-26 | Verizon Laboratories Inc. | Method and system for characterizing propagation of radiofrequency signal |
-
2000
- 2000-08-21 JP JP2000250028A patent/JP3648140B2/ja not_active Expired - Fee Related
-
2001
- 2001-08-08 US US09/924,235 patent/US7187324B2/en not_active Expired - Fee Related
- 2001-08-20 DE DE10140813A patent/DE10140813A1/de not_active Withdrawn
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004531736A (ja) * | 2001-06-25 | 2004-10-14 | ハリス コーポレイション | 受動的な反射器又は屈折器を代理受信機として使用するとともに、データベースへの問合せを利用した送信機の位置を測定するためのシステム及び方法 |
| US8073464B2 (en) | 2002-10-22 | 2011-12-06 | Qualcomm Incorporated | Method and apparatus for identifying transmitters in a wireless communication system using power predictions |
| JP2006515727A (ja) * | 2002-10-22 | 2006-06-01 | クゥアルコム・インコーポレイテッド | 無線通信システムにおいて予測電力を用いて送信機を識別する方法および装置 |
| JP2006519389A (ja) * | 2003-03-04 | 2006-08-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | オブジェクトロケーション |
| JP4667368B2 (ja) * | 2003-03-04 | 2011-04-13 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 無線局の位置決め方法、コンピュータプログラム、無線局及び無線システム |
| JP4544869B2 (ja) * | 2004-01-22 | 2010-09-15 | 総務大臣 | 電波ホログラフィ法とｍｕｓｉｃ法を用い電波源の特定に用いられる電波監視装置 |
| JP2005207847A (ja) * | 2004-01-22 | 2005-08-04 | Ministry Of Public Management Home Affairs Posts & Telecommunications | 電波ホログラフィ法とｍｕｓｉｃ法を用い電波源の特定に用いられる電波監視装置 |
| WO2007058301A1 (ja) * | 2005-11-21 | 2007-05-24 | Nec Corporation | 検知対象の識別・位置推定システムと、その方法及びプログラム |
| JP5170389B2 (ja) * | 2005-11-21 | 2013-03-27 | 日本電気株式会社 | 検知対象の識別・位置推定システムと、その方法及びプログラム |
| US8624707B2 (en) | 2005-11-21 | 2014-01-07 | Nec Corporation | Detection target identifying/position estimating system, its method, and program |
| JP2010525300A (ja) * | 2007-03-13 | 2010-07-22 | アルカテル−ルーセント ユーエスエー インコーポレーテッド | 逆方向レイ追跡法を用いて送信機を位置特定する方法 |
| JP2009111952A (ja) * | 2007-11-01 | 2009-05-21 | Clarion Co Ltd | 放送信号受信機試験装置 |
| JP2012118024A (ja) * | 2010-12-03 | 2012-06-21 | Brother Ind Ltd | 電波推定方法、及び電波推定プログラム、並びに電波推定装置 |
| JP2012244590A (ja) * | 2011-05-24 | 2012-12-10 | Nippon Telegr & Teleph Corp <Ntt> | 到来方向推定装置、到来方向推定方法、及び到来方向推定プログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| DE10140813A1 (de) | 2002-04-25 |
| US20020115411A1 (en) | 2002-08-22 |
| US7187324B2 (en) | 2007-03-06 |
| JP3648140B2 (ja) | 2005-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3648140B2 (ja) | 電波監視方法及び装置 | |
| Schillebeeckx et al. | Biomimetic sonar: Binaural 3D localization using artificial bat pinnae | |
| CN105388478A (zh) | 用于检测声学和光学信息的方法和装置、以及对应的计算机程序和对应的计算机可读存储介质 | |
| CN113064147B (zh) | 一种低信噪比下新型匹配场被动定位方法 | |
| KR20070034205A (ko) | 3차원 광선 추적법을 이용한 전파 모델 생성 장치, 방법 및이를 구현한 컴퓨터로 실행 가능한 기록매체 | |
| JP3718342B2 (ja) | 広域電波監視方法及び装置 | |
| CN114003981A (zh) | 一种基于时空一体数字地球的电磁频谱可视化分析方法 | |
| CN112098938B (zh) | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 | |
| CN111107626A (zh) | 一种基于时间反演的doa定位方法 | |
| Norrdine et al. | Radar wave based positioning of spatially distributed MIMO radar antenna systems for near-field nondestructive testing | |
| MX2013015016A (es) | Indicador de repetibilidad 4d basado en iluminacion por disparo par adquision sismica. | |
| White et al. | Localisation of sperm whales using bottom-mounted sensors | |
| JP5997007B2 (ja) | 音源位置推定装置 | |
| JP7116949B2 (ja) | 非接触形状・誘電率測定装置、非接触形状・誘電率測定方法及びプログラム | |
| US10375501B2 (en) | Method and device for quickly determining location-dependent pulse responses in signal transmission from or into a spatial volume | |
| US6188356B1 (en) | Radio monitoring method and apparatus | |
| Åström et al. | Extension of time-difference-of-arrival self calibration solutions using robust multilateration | |
| CN113534161A (zh) | 一种用于远距离定位水下声源的波束镜像聚焦方法 | |
| Dubrovinskaya et al. | Anchorless underwater acoustic localization | |
| CN115840192A (zh) | 一种基于空间估计谱置信度估计的室内定位方法 | |
| KR100927962B1 (ko) | 경로 전처리를 이용한 광선 추적 장치 및 방법 | |
| KR102030526B1 (ko) | 연속 압축센싱 기반 도래각 추정 방법 및 장치 | |
| JP2002107440A (ja) | 方向検出方法及び方向検出装置 | |
| JP7241993B2 (ja) | 電波源位置推定装置、電波源位置推定システム、及び電波源位置推定方法 | |
| Esrafilian et al. | Global Scale Self-Supervised Channel Charting with Sensor Fusion |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040318 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20041228 |
|
| 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: 20050208 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050210 |
|
| 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: 20090218 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090218 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100218 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110218 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120218 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120218 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130218 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140218 Year of fee payment: 9 |
|
| LAPS | Cancellation because of no payment of annual fees |