JP4798531B2 - 受信局、それを用いた信号送受信方式 - Google Patents
受信局、それを用いた信号送受信方式 Download PDFInfo
- Publication number
- JP4798531B2 JP4798531B2 JP2005327242A JP2005327242A JP4798531B2 JP 4798531 B2 JP4798531 B2 JP 4798531B2 JP 2005327242 A JP2005327242 A JP 2005327242A JP 2005327242 A JP2005327242 A JP 2005327242A JP 4798531 B2 JP4798531 B2 JP 4798531B2
- Authority
- JP
- Japan
- Prior art keywords
- radio wave
- beam pattern
- pattern
- receiving station
- station
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Description
本発明は、水中からの微弱電波を受信したときに、水中から空中へ伝搬するときの影響や障害物の影響を修正する受信局、および前記受信局を用い、アンテナの回転角度に応じた受信電界強度の変化を測定し、発信源の位置を特定する信号送受信方式に関する。
（水中からの電波を受信する位置を最適化する方法）
従来、野生生物、特に魚類の自動追尾を行う方法は、魚類に電波発信機を装着し、行動を追跡するシステムを調査地内に設置し、
（１）目的とするサービスエリアから発信される電波を効率よく受信できる地点を絞り込み、
（２）受信局設置候補地点におけるカバーエリアを特定する手順からなる。
上記手順（１）（２）は、電波発信源と受信局設置候補地点間の電波伝搬解析に基づいて行われる手順である。従来、上記（１）（２）の電波伝搬解析方法は、大地伝搬式や奥村式といった比較的伝搬過程を単純化した電波伝搬モデルにより電波伝搬解析を行っていた。
従来、野生生物、特に魚類の自動追尾を行う方法は、魚類に電波発信機を装着し、行動を追跡するシステムを調査地内に設置し、
（１）目的とするサービスエリアから発信される電波を効率よく受信できる地点を絞り込み、
（２）受信局設置候補地点におけるカバーエリアを特定する手順からなる。
上記手順（１）（２）は、電波発信源と受信局設置候補地点間の電波伝搬解析に基づいて行われる手順である。従来、上記（１）（２）の電波伝搬解析方法は、大地伝搬式や奥村式といった比較的伝搬過程を単純化した電波伝搬モデルにより電波伝搬解析を行っていた。
（水中から発信される電波到来角を高精度で推定する方法）
また、指向性アンテナを回転させ、間欠発信を行う電波発信機からの電波到来角を推定する方式（例えば、特許文献１参照）では、間欠発信波の発信間隔の間は信号が無くなり、受信特性曲線の頂点を検出することは困難になる。このため、高精度の電波到来角の推定ができなかった。
特開平７−６３５５９号公報
また、指向性アンテナを回転させ、間欠発信を行う電波発信機からの電波到来角を推定する方式（例えば、特許文献１参照）では、間欠発信波の発信間隔の間は信号が無くなり、受信特性曲線の頂点を検出することは困難になる。このため、高精度の電波到来角の推定ができなかった。
電波発信機を用いて魚類の行動を自動追尾する場合、自然環境下、特に水中から空中へ電波が伝搬する電波伝搬過程は複雑な状態となる。
例えば、空中の電波伝搬特性の場合、電波発信源から電波受信源へ発信される電波は、直達波が最も強いが反射屈折波もかなりの強さで伝搬する。伝搬路は複数パス存在するが、発信源近くでの散乱が少なく比較的単純な伝搬特性を示す。
これに対し、水中の電波伝搬特性の場合、直達波が反射屈折波と比較して強く伝搬する。このため、受信波のサイドローブが大きく減少する。また、水中から出る時点で激しく散乱し、水中からの発信のため、電波自体が弱く、散乱し複雑な伝搬パスを形成する。反射・屈折等の影響も受け易い。
例えば、空中の電波伝搬特性の場合、電波発信源から電波受信源へ発信される電波は、直達波が最も強いが反射屈折波もかなりの強さで伝搬する。伝搬路は複数パス存在するが、発信源近くでの散乱が少なく比較的単純な伝搬特性を示す。
これに対し、水中の電波伝搬特性の場合、直達波が反射屈折波と比較して強く伝搬する。このため、受信波のサイドローブが大きく減少する。また、水中から出る時点で激しく散乱し、水中からの発信のため、電波自体が弱く、散乱し複雑な伝搬パスを形成する。反射・屈折等の影響も受け易い。
また、指向性アンテナを回転させ電波到来角を推定する方法において、遮蔽物があった場合や水中からの電波発信等の特殊な電波伝搬を阻害する要因があった場合、ビームパターンの変形や一部情報の欠損が生じる。
本発明は、前記問題点に鑑み、信号欠損部分の状態を判断し、電波到来角度の推定精度を高めると共に、電波伝搬路を構成する水中の過程を補正するための新しい電波伝搬モデルを組み込んだ受信局、それを用いた信号送受信方式を提供することを目的とする。
本発明は、前記問題点に鑑み、信号欠損部分の状態を判断し、電波到来角度の推定精度を高めると共に、電波伝搬路を構成する水中の過程を補正するための新しい電波伝搬モデルを組み込んだ受信局、それを用いた信号送受信方式を提供することを目的とする。
本発明の受信局は、
予め、遮蔽物のない条件での受信ビームパターンをモデルビームパターンとして取り込んでおく手段１、
次に、実測データを取り込み、所定受信期間内の全受信データの内から最大値を検索する手段２、
検索した最大値を含む放物線の推定ビームパターンをモデルビームパターンと対比して特定する手段３、
前記手段３における推定ビームパターンを成す放物線の末端に近い受信データから、近い受信データ間の階差をとり、局所トレンドを推定する手段４、
前記手段４における局所トレンドの推定処理の結果、前記推定ビームパターンの正常範囲内の値となる受信データを候補に設定する手段５、
前記手段５で抽出した候補の内で、前記推定ビームパターンに近い受信データを推定母集団の大勢を決めるサンプルと仮定し、そのサンプルの中から関数を生成する手段６、
前記手段６で生成した関数を最も可能性があるビームパターンの関数とし、そのビームパターンの関数の頂点を電波到来角として抽出する手段７、
を有することを特徴とする。
予め、遮蔽物のない条件での受信ビームパターンをモデルビームパターンとして取り込んでおく手段１、
次に、実測データを取り込み、所定受信期間内の全受信データの内から最大値を検索する手段２、
検索した最大値を含む放物線の推定ビームパターンをモデルビームパターンと対比して特定する手段３、
前記手段３における推定ビームパターンを成す放物線の末端に近い受信データから、近い受信データ間の階差をとり、局所トレンドを推定する手段４、
前記手段４における局所トレンドの推定処理の結果、前記推定ビームパターンの正常範囲内の値となる受信データを候補に設定する手段５、
前記手段５で抽出した候補の内で、前記推定ビームパターンに近い受信データを推定母集団の大勢を決めるサンプルと仮定し、そのサンプルの中から関数を生成する手段６、
前記手段６で生成した関数を最も可能性があるビームパターンの関数とし、そのビームパターンの関数の頂点を電波到来角として抽出する手段７、
を有することを特徴とする。
また、上記受信局と制御局を備えた信号送受信方式であって、
前記受信局各局により電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を受信し、受信信号を前記制御局へ送信し、
前記制御局は前記受信信号を取り込み、効率よく受信できる地点を演算し、
受信局設置候補地点におけるカバーエリアを特定し、
前記カバーエリア同士の重なり合いを判断（閾値：例えばカバーエリアに対するサービスエリアの比率等の閾値により）し、目的とするサービスエリアを確保するための受信局の最適配置を求める。
前記受信局各局により電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を受信し、受信信号を前記制御局へ送信し、
前記制御局は前記受信信号を取り込み、効率よく受信できる地点を演算し、
受信局設置候補地点におけるカバーエリアを特定し、
前記カバーエリア同士の重なり合いを判断（閾値：例えばカバーエリアに対するサービスエリアの比率等の閾値により）し、目的とするサービスエリアを確保するための受信局の最適配置を求める。
また、信号送受信方式は、上記受信局を用い、魚類に装着した前記電波発信機の発信電波を受信し、前記受信局と前記制御局とにより電波発信機の位置を追跡する信号送受信方式であって、
前記受信局各局は、前記信号送受信方式により前記電波伝搬モデルを用いて電波到来角度θを求めて前記制御局へ送信し、
前記制御局は、前記受信局から前記電波到来角度θを取り込み、三角計測により、発信源の位置を抽出する。
前記受信局各局は、前記信号送受信方式により前記電波伝搬モデルを用いて電波到来角度θを求めて前記制御局へ送信し、
前記制御局は、前記受信局から前記電波到来角度θを取り込み、三角計測により、発信源の位置を抽出する。
（２）上記（１）記載の受信局を用い、魚類に装着した電波発信機の発信電波を受信し、前記受信局と前記制御局とにより前記電波発信機の位置を追跡する信号送受信方式において、
前記受信局各局は、前記信号送受信方式により前記電波伝搬モデルを用いて電波到来角度θを求めて前記制御局へ送信し、前記制御局は、前記受信局から前記電波到来角度θを取り込み、三角計測により、発信源の位置を抽出することを特徴とする。
前記受信局各局は、前記信号送受信方式により前記電波伝搬モデルを用いて電波到来角度θを求めて前記制御局へ送信し、前記制御局は、前記受信局から前記電波到来角度θを取り込み、三角計測により、発信源の位置を抽出することを特徴とする。
（３）上記（２）記載の信号送受信方式において、
前記受信局は、受信信号のビームパターンをビーム伝送経路に遮蔽物がある場合のパターン、前記ビーム伝送経路により歪みが生じた場合のパターン、およびその他の場合のパターンのいずれかに判定し、前記各パターンに応じて前記各パターンに補正を加えることにより、修正したビームパターンを求め、該修正したビームパターンから前記電波到来角度θを求めることを特徴とする。
前記受信局は、受信信号のビームパターンをビーム伝送経路に遮蔽物がある場合のパターン、前記ビーム伝送経路により歪みが生じた場合のパターン、およびその他の場合のパターンのいずれかに判定し、前記各パターンに応じて前記各パターンに補正を加えることにより、修正したビームパターンを求め、該修正したビームパターンから前記電波到来角度θを求めることを特徴とする。
（４）上記（２）又は３記載の受信局と制御局を備えた信号送受信方式において、
シミュレーションにより、前記受信局各局により電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を受信し、受信信号を前記制御局へ送信し、前記制御局は前記受信信号を取り込み、効率よく受信できる地点を演算し、受信局設置候補地点におけるカバーエリアを特定し、前記カバーエリア同士の重なり合いを判断し、目的とするサービスエリアを確保するための前記受信局の最適配置位置を求めることを特徴とする。
シミュレーションにより、前記受信局各局により電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を受信し、受信信号を前記制御局へ送信し、前記制御局は前記受信信号を取り込み、効率よく受信できる地点を演算し、受信局設置候補地点におけるカバーエリアを特定し、前記カバーエリア同士の重なり合いを判断し、目的とするサービスエリアを確保するための前記受信局の最適配置位置を求めることを特徴とする。
本発明の信号送受信方式は、水中から空気中への電波伝搬特性の影響、および障害物による受信信号の影響を修正した電波伝搬モデルを用いるので、正確な受信信号を得ることができる。
また、本発明の電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を効率よく受信できる地点を絞り込み、該地点に受信局を配置する方法を採用することにより、発信電波を効率よく受信することができる。
また、上記修正した電波伝搬モデルを有する信号送受信方式を用いて電波到来角度を推定する方式にしたので、正確な電波到来角度を推定することができる。
また、本発明の電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を効率よく受信できる地点を絞り込み、該地点に受信局を配置する方法を採用することにより、発信電波を効率よく受信することができる。
また、上記修正した電波伝搬モデルを有する信号送受信方式を用いて電波到来角度を推定する方式にしたので、正確な電波到来角度を推定することができる。
本発明の実施の形態を以下図に基づいて詳細に説明する。
本発明の信号送受信方式は、基本的に、以下の手順（ステップ）をとる。
図１は本発明の信号送受信方式の動作態様を示すフローチャートである。
フローチャートは、大きく分けて、２つの手順をとる。まず、受信局設置位置の検討の手順と、その次のシステムの運用の手順である。
受信局設置位置の検討の手順の内容は、
（手順１）電波伝搬シミュレーションにより受信局設置予定位置を含むシステムサービスエリアの抽出を行う（ステップＳ１）：
調査対象地域を選定し、システムサービスエリアを選定する。
（手順２）設置予定位置を含むシステムサービスエリアと測定手法との整合性の検証を行う（ステップＳ２）：
サービスエリア内の任意の設置予定位置から電波伝搬シミュレーションを行い、受信局候補エリアにおける受信エリア解析を行う。
（手順３）検証結果を反映してシステムの設置・校正を行う（ステップＳ３）：
受信エリア解析の結果から受信局をサービスエリア内のどの位置に設置すべきか決定する。
本発明の信号送受信方式は、基本的に、以下の手順（ステップ）をとる。
図１は本発明の信号送受信方式の動作態様を示すフローチャートである。
フローチャートは、大きく分けて、２つの手順をとる。まず、受信局設置位置の検討の手順と、その次のシステムの運用の手順である。
受信局設置位置の検討の手順の内容は、
（手順１）電波伝搬シミュレーションにより受信局設置予定位置を含むシステムサービスエリアの抽出を行う（ステップＳ１）：
調査対象地域を選定し、システムサービスエリアを選定する。
（手順２）設置予定位置を含むシステムサービスエリアと測定手法との整合性の検証を行う（ステップＳ２）：
サービスエリア内の任意の設置予定位置から電波伝搬シミュレーションを行い、受信局候補エリアにおける受信エリア解析を行う。
（手順３）検証結果を反映してシステムの設置・校正を行う（ステップＳ３）：
受信エリア解析の結果から受信局をサービスエリア内のどの位置に設置すべきか決定する。
このように、受信局設置位置の検討の手順を実行して、受信局の最適な設置場所を見つける。受信局の設置位置の決定が終わった段階で次のシステムの運用の手順を実行し、受信を始め、電波発信機の位置を所定のアルゴリズムに基づくプログラムにより求める。すなわち、
（手順４）三角測量による電波発信機の位置を特定する（ステップＳ４）：
（手順５）電波到来角度の推定（演算）を高精度に行う（ステップ）Ｓ５：
（手順６）三角測量手法により電波発信機の位置を特定する（ステップＳ６）：
（手順４）三角測量による電波発信機の位置を特定する（ステップＳ４）：
（手順５）電波到来角度の推定（演算）を高精度に行う（ステップ）Ｓ５：
（手順６）三角測量手法により電波発信機の位置を特定する（ステップＳ６）：
以上の手順は、すべてシミュレーションにより実行する場合を説明する。これとは別に、実際の現場の予め候補地点に実際の受信機を配置し、移動体モデルを動かしてサンプリングデータを取得することにより実行することもできる。
詳しくは、以下前記電波伝搬モデルを用いた位置特定原理に基づいて説明する。
(1) 予め、地理的条件による伝搬損失を考慮した特定の電波伝搬モデルを作成する。
受信局に特定の電波伝搬モデルを取り込み、
（２）仮想の電波発信機を想定調査地に設置し、前記特定の電波伝搬モデルを用いて調査候補地内の任意の位置への電波伝搬状態をシミュレーションする。
（３）伝搬路の状態をＧＩＳを用いて解析（見通し解析、伝搬路内の地形要素による修正を加えた解析等）して修正条件を抽出し、前記修正条件により前記特定の電波伝搬モデルを修正する。
（４）土木研究所が算定した水中での電波減衰量データに基づき、想定調査地内の水中での電波減衰量や遮蔽物による電波減衰量、水中から空中への電波減衰量等を推定し、前記特定の電波伝搬モデルを修正する。
（５）修正の済んだ電波伝搬モデルを用いて電波発信機を調査地内の複数の調査地点に設置し、調査地内のサンプリング位置への電波伝搬状態を再度シミュレーションする。
（６）調査地内のサンプリング位置の電波受信特性に基づき、該調査地内で最も電波受信確率が高いエリアＡを抽出する。
（７）前記エリアＡが、ＡＴＳの位置特定方法を実行した場合の適切な位置となるか否か判断する。
（８）上記（７）のステップにおいて、エリアＡが、ＡＴＳの位置特定方法を実行した場合の適切な位置と判断されたとき、エリアＡを設置場所の候補とする。
（９）制御器の指令によりエリアＡの設置場所の特定が済んだ受信局によりシミュレーションして発信電波を受信し、到来角度θを計測し、制御局へ送信する。
（１０）制御局は、各受信局からの到来角度θを基に、三点計測（三角測量）演算を行い、発信位置を特定する。
(1) 予め、地理的条件による伝搬損失を考慮した特定の電波伝搬モデルを作成する。
受信局に特定の電波伝搬モデルを取り込み、
（２）仮想の電波発信機を想定調査地に設置し、前記特定の電波伝搬モデルを用いて調査候補地内の任意の位置への電波伝搬状態をシミュレーションする。
（３）伝搬路の状態をＧＩＳを用いて解析（見通し解析、伝搬路内の地形要素による修正を加えた解析等）して修正条件を抽出し、前記修正条件により前記特定の電波伝搬モデルを修正する。
（４）土木研究所が算定した水中での電波減衰量データに基づき、想定調査地内の水中での電波減衰量や遮蔽物による電波減衰量、水中から空中への電波減衰量等を推定し、前記特定の電波伝搬モデルを修正する。
（５）修正の済んだ電波伝搬モデルを用いて電波発信機を調査地内の複数の調査地点に設置し、調査地内のサンプリング位置への電波伝搬状態を再度シミュレーションする。
（６）調査地内のサンプリング位置の電波受信特性に基づき、該調査地内で最も電波受信確率が高いエリアＡを抽出する。
（７）前記エリアＡが、ＡＴＳの位置特定方法を実行した場合の適切な位置となるか否か判断する。
（８）上記（７）のステップにおいて、エリアＡが、ＡＴＳの位置特定方法を実行した場合の適切な位置と判断されたとき、エリアＡを設置場所の候補とする。
（９）制御器の指令によりエリアＡの設置場所の特定が済んだ受信局によりシミュレーションして発信電波を受信し、到来角度θを計測し、制御局へ送信する。
（１０）制御局は、各受信局からの到来角度θを基に、三点計測（三角測量）演算を行い、発信位置を特定する。
本発明の電波伝搬モデルを用いた発信源位置特定方式は、以下の手順を有する。
魚類に電波発信機を装着し、受信局および制御局により、行動を追跡する方法および方式であって、
（１）電波伝搬路の状況、即ち、伝搬特性に影響を与える状況、例えば水中から空中への伝搬、障害物に遮られる伝搬等の状況、に応じた独自の電波伝搬モデルを作成する。
（２）前記電波伝搬モデルを用いて、目的とするサービスエリアから発信される電波を効率よく受信できる地点を絞り込む。
（３）受信局設置候補地点におけるカバーエリアを特定する。
（４）カバーエリア同士の重なり合いを判断し、目的とするサービスエリアを確保するための受信局の最適配置を行う。
（５）各受信局において、前記電波伝搬モデルを用いて、発信電波の到来角度θを抽出し、制御局へ送信する。
（６）制御局は、各受信局からの前記到来角度θを基に、三角計測（三角測量）により発信位置を特定する。
魚類に電波発信機を装着し、受信局および制御局により、行動を追跡する方法および方式であって、
（１）電波伝搬路の状況、即ち、伝搬特性に影響を与える状況、例えば水中から空中への伝搬、障害物に遮られる伝搬等の状況、に応じた独自の電波伝搬モデルを作成する。
（２）前記電波伝搬モデルを用いて、目的とするサービスエリアから発信される電波を効率よく受信できる地点を絞り込む。
（３）受信局設置候補地点におけるカバーエリアを特定する。
（４）カバーエリア同士の重なり合いを判断し、目的とするサービスエリアを確保するための受信局の最適配置を行う。
（５）各受信局において、前記電波伝搬モデルを用いて、発信電波の到来角度θを抽出し、制御局へ送信する。
（６）制御局は、各受信局からの前記到来角度θを基に、三角計測（三角測量）により発信位置を特定する。
上記発信源の位置特定のためには、電波伝搬モデルを必要とする。その必要性について説明する。
図２は遮蔽物に遮られる場合の電波伝搬特性の説明図である。
発信機と受信局との間に遮蔽物が無い場合の受信特性を検討する。
水中から発信された連続発信波を、アンテナを回転させながら受信すると、図２（ａ）に示すアンテナ回転角に対する受信電界強度特性のビームパターンとなり、指向性のよい特性になる。
一方、水中からの発信電波が間欠発信波の場合、受信信号パターンは図２（ｂ）のようになる。
図２（ｂ）の離散値信号形状は、図２（ａ）に示す連続発信信号の変化の特徴を有する。
図２は遮蔽物に遮られる場合の電波伝搬特性の説明図である。
発信機と受信局との間に遮蔽物が無い場合の受信特性を検討する。
水中から発信された連続発信波を、アンテナを回転させながら受信すると、図２（ａ）に示すアンテナ回転角に対する受信電界強度特性のビームパターンとなり、指向性のよい特性になる。
一方、水中からの発信電波が間欠発信波の場合、受信信号パターンは図２（ｂ）のようになる。
図２（ｂ）の離散値信号形状は、図２（ａ）に示す連続発信信号の変化の特徴を有する。
次に、発信機と受信局との間の、発信機周辺に大きな遮蔽物（植生、地形）が有る場合の受信特性を検討する。
図２ｃは、水中から発信された連続発信波が途中の遮蔽物を介して受信された場合の、受信信号をサンプリング表示したアンテナ回転角に応じた受信電界強度特性を示す。図中枠で囲った領域の特性は、先の図２（ａ）の特性を保持することができず、崩れた形になっている。
図２（ｄ）は、水中から発信された間欠発信波が途中の遮蔽物を介して受信された場合のアンテナ回転角に応じた受信電界強度特性を示す。図中枠で囲った領域は、受信データのバラツキが生じている。間欠発信波を用いる場合、連続発信波を用いる場合に比べて受信データのバラツキが少ない傾向を示す。
図２ｃは、水中から発信された連続発信波が途中の遮蔽物を介して受信された場合の、受信信号をサンプリング表示したアンテナ回転角に応じた受信電界強度特性を示す。図中枠で囲った領域の特性は、先の図２（ａ）の特性を保持することができず、崩れた形になっている。
図２（ｄ）は、水中から発信された間欠発信波が途中の遮蔽物を介して受信された場合のアンテナ回転角に応じた受信電界強度特性を示す。図中枠で囲った領域は、受信データのバラツキが生じている。間欠発信波を用いる場合、連続発信波を用いる場合に比べて受信データのバラツキが少ない傾向を示す。
以上述べたように、遮蔽物によって受信電界強度特性にバラツキが発生する現象が起こる。この現象に共通することは、
（１）電波の遮蔽、瞬時値変動によりビームパターンの形が微妙に変動する、
（２）受信電界強度の最大値を含む受信電圧強度領域のサンプルデータは、ビームパターンの原型に近いか、若しくは、ビームパターンの原型に近い関数に乗る。
以上の検討結果に着目すると、以下のようにまとめることができる。
（１）受信データが離散値データで、且つ、回転角が変化するにつれて受信電界強度が指向性を持って、換言すると、概略、先鋭な放物線で変化するので、隣接離散値データを順次前記先鋭な放物線に従って設定することにより、適切な受信データを選択して受信精度を向上することができる。
（２）電波伝搬過程の影響を受けたデータを除外すれば、数値補間の精度が向上する。
（１）電波の遮蔽、瞬時値変動によりビームパターンの形が微妙に変動する、
（２）受信電界強度の最大値を含む受信電圧強度領域のサンプルデータは、ビームパターンの原型に近いか、若しくは、ビームパターンの原型に近い関数に乗る。
以上の検討結果に着目すると、以下のようにまとめることができる。
（１）受信データが離散値データで、且つ、回転角が変化するにつれて受信電界強度が指向性を持って、換言すると、概略、先鋭な放物線で変化するので、隣接離散値データを順次前記先鋭な放物線に従って設定することにより、適切な受信データを選択して受信精度を向上することができる。
（２）電波伝搬過程の影響を受けたデータを除外すれば、数値補間の精度が向上する。
以上の検討結果を踏まえると、水中から発信される電波の到来角度を高精度で推定する際の基本的な考え方は、以下のようになる。
図３はビームパターンの説明図である。図３（ａ）はパターンＡの説明図、図３（ｂ）はメインローブとサイドローブを有する基本波形のビームパターンの説明図、図３（ｃ）はパターンＢの説明図である。図中、実線は実測ビームパターン、点線は電波伝搬による撹乱がない場合のビームパターンを表す。
以下に述べるビームパターンの判定手段およびビームパターンの修正手段は、ビーム伝搬経路に遮蔽物がある場合と無い場合（理想的な場合）、またはビーム伝搬経路上の問題でビームパターンがゆがむ場合とその問題が無くビームパターンがゆがまない場合との比較を行い、必要な修正データを採取し、予めテーブルに記憶してある。
図３はビームパターンの説明図である。図３（ａ）はパターンＡの説明図、図３（ｂ）はメインローブとサイドローブを有する基本波形のビームパターンの説明図、図３（ｃ）はパターンＢの説明図である。図中、実線は実測ビームパターン、点線は電波伝搬による撹乱がない場合のビームパターンを表す。
以下に述べるビームパターンの判定手段およびビームパターンの修正手段は、ビーム伝搬経路に遮蔽物がある場合と無い場合（理想的な場合）、またはビーム伝搬経路上の問題でビームパターンがゆがむ場合とその問題が無くビームパターンがゆがまない場合との比較を行い、必要な修正データを採取し、予めテーブルに記憶してある。
電波到来角度を推定する場合問題となる受信パターンが２種類（Ａ、Ｂ）ある。
パターンＡ：受信前のビーム経路に遮蔽物がある場合のビームパターン、
パターンＢ：電波伝搬上の問題でビームパターンがゆがむビームパターン。
パターンＡ：受信前のビーム経路に遮蔽物がある場合のビームパターン、
パターンＢ：電波伝搬上の問題でビームパターンがゆがむビームパターン。
パターンＡは、図３（ａ）に示されているように、撹乱がない場合点線の特性Ａ１となるが、遮蔽物があるため、実線のように２つのビームパターンになる。そのうちの受信電界強度が大きい方の波形の最大値の回転角θ１が実際の現地計測で得られるビームパターンとなる。なお、実線の小さいビームパターンの最大値の回転角はθ２とする。一方、電波伝搬による攪乱がない場合の点線のビームパターンの受信電界強度が最大の回転角はθＴとなる。パターンＡの場合、θ１とθＴとが一致しないため、実測値のθ１をアルゴリズムを修正してθＴに一致させる。
パターンＢは、図３（ｃ）に示されているように、実線の実際の現地計測で得られるビームパターンの最大値の回転角はθＣとなり、点線の電波伝搬による攪乱がない場合のビームパターンの最大値の回転角θＴとが一致しないため、実測値のθＣをアルゴリズムを修正してθＴに一致させる。点線のビームパターンの最大値の回転角はθＴで、実線のビームパターンの最大値（Ｖ１Ｃ）の回転角はθＣとする。
θＣとＶ１Ｃを結んだ線と点線のビームパターンとが交差する点を通る等受信電界強度直線と実線のビームパターンとの交わる点（Ｖ１Ａ）の回転角θ１を求め、回転角の軸線上で、θ２−θＣ＝θＣ−θ１となるようにθ２を設け、回転角θ２の直線と実線のビームパターンの交わる点の受信電界強度Ｖ２Ａを求める。
θＣとＶ１Ｃを結んだ線と点線のビームパターンとが交差する点を通る等受信電界強度直線と実線のビームパターンとの交わる点（Ｖ１Ａ）の回転角θ１を求め、回転角の軸線上で、θ２−θＣ＝θＣ−θ１となるようにθ２を設け、回転角θ２の直線と実線のビームパターンの交わる点の受信電界強度Ｖ２Ａを求める。
このアルゴリズム修正のためには、下記（１）（２）を行う。
（１）ビームパターンの分類による処理アルゴリズムの選択、
図３（ｂ）に示すように、真の電波到来方向のビームパターンは、サイドローブのビームパターンよりも大きく、物理的に決まる。図中、Ｖ１：受信電界強度の最大値、Ｖ２：受信電界強度最大値から２番目の値、Ｗ：メインローブ幅を表す。また、テレメトリシステムのアンテナ回転速度によりその基本形は影響を受け、測定条件により基本的なビームパターンの放物線に近い関数が決定される。
（１）ビームパターンの分類による処理アルゴリズムの選択、
図３（ｂ）に示すように、真の電波到来方向のビームパターンは、サイドローブのビームパターンよりも大きく、物理的に決まる。図中、Ｖ１：受信電界強度の最大値、Ｖ２：受信電界強度最大値から２番目の値、Ｗ：メインローブ幅を表す。また、テレメトリシステムのアンテナ回転速度によりその基本形は影響を受け、測定条件により基本的なビームパターンの放物線に近い関数が決定される。
（ビームパターン判断のアルゴリズム）
（ａ）ビームパターンの判断（パターンＡ）：
図３（ａ）および図３（ｂ）において、θｄ＝θ１−θ２とすると、
条件１：θｄ＜Ｗ
且つ
条件２：Ｖｄａ＜Ｖｐ
であれば、パターンＡといえる。
（ａ）ビームパターンの判断（パターンＡ）：
図３（ａ）および図３（ｂ）において、θｄ＝θ１−θ２とすると、
条件１：θｄ＜Ｗ
且つ
条件２：Ｖｄａ＜Ｖｐ
であれば、パターンＡといえる。
（ｂ）ビームパターンの判断（パターンＢ）：
図Ｃ（ｃ）において、
条件１：
但し、数１式中のＣ１は現地の状況に応じて変化する値で、例えば１．２をとる。
又は
但し、数２式中のＣ２は現地の状況に応じて変化する値で、例えば１．２をとる。
且つ、
図Ｃ（ｃ）において、
条件１：
又は
且つ、
条件２：
但し、数３式中のＣ３は現地の状況に応じて変化する値で、例えば１．２をとる。
又は
但し、数４式中のＣ４は現地の状況に応じて変化する値で、例えば１．２をとる。
であれば、パターンBといえる。
上記（ａ）または（ｂ）の条件を満たさない場合には、図３（ｂ）のビームパターンとする。
又は
であれば、パターンBといえる。
上記（ａ）または（ｂ）の条件を満たさない場合には、図３（ｂ）のビームパターンとする。
（２）実際の処理アルゴリズムのアレンジ（例えば、パターンＡに対する補足、等）を必要とする。→すなわち、簡単な数値処理から典型的なエラーを生じさせるビームパターンか否かを判断し、アルゴリズムを選択する。パターンＡ、またはパターンＢに該当する場合には、それぞれに対応するアルゴリズムを選択・判断し、必要に応じて補正処理等を行う。
（ビームパターン補正アルゴリズム）
（ａ）パターンＡ：
図３（ａ）に基づいて説明する。
（１）実線の実測特性の各頂点（最大値）を基準としてノイズレベルの回転角θＡ、θＢを求める。ノイズレベルは、例えば最大値の２０パーセントの受信電界強度とする。
（２）θＴ＝（θＡ＋θＢ）／２を求める。
（３）θＡ＜θ＜θＴの条件を満たし、Δ（Ｖ／θ）＞０になる回転角θ部分を検索・選択する。
（４）θＴ＜θ＜θＢの条件を満たし、Δ（Ｖ／θ）＜０になる部分を検索・選択する。
（５）（３）と（４）を組み合わせてビームパターンを合成する（（３）と（４）により、実測の特性曲線Ｐ１の左側領域Ａ１の傾向および実測の特性曲線Ｐ２の右側領域Ａ２の傾向を含み左右対称な特性曲線を求める）。
（６）合成したビームパターンの最大値の回転角を電波到来角として求める。
（ビームパターン補正アルゴリズム）
（ａ）パターンＡ：
図３（ａ）に基づいて説明する。
（１）実線の実測特性の各頂点（最大値）を基準としてノイズレベルの回転角θＡ、θＢを求める。ノイズレベルは、例えば最大値の２０パーセントの受信電界強度とする。
（２）θＴ＝（θＡ＋θＢ）／２を求める。
（３）θＡ＜θ＜θＴの条件を満たし、Δ（Ｖ／θ）＞０になる回転角θ部分を検索・選択する。
（４）θＴ＜θ＜θＢの条件を満たし、Δ（Ｖ／θ）＜０になる部分を検索・選択する。
（５）（３）と（４）を組み合わせてビームパターンを合成する（（３）と（４）により、実測の特性曲線Ｐ１の左側領域Ａ１の傾向および実測の特性曲線Ｐ２の右側領域Ａ２の傾向を含み左右対称な特性曲線を求める）。
（６）合成したビームパターンの最大値の回転角を電波到来角として求める。
（ｂ）パターンＢ：
図３（ｃ）に基づいて説明する。
（１）探索限界値の算出を行う。算出式は（計測最大値−ノイズ平均）／２とする。ノイズ平均は頻度分布に基づいて求める。
（２）探索限界値以上のパルスデータを抽出する。
（３）抽出したパルスデータを用い、以下のビームパターン推定アルゴリズムに基づきビームパターンの最大値を求め、その回転角を電波到来角として求める（θAとθBを端点として（θAとθBの間隔）、最大値をＶ１Ｃとする左右対称な関数式をテーブルから求める）。
図３（ｃ）に基づいて説明する。
（１）探索限界値の算出を行う。算出式は（計測最大値−ノイズ平均）／２とする。ノイズ平均は頻度分布に基づいて求める。
（２）探索限界値以上のパルスデータを抽出する。
（３）抽出したパルスデータを用い、以下のビームパターン推定アルゴリズムに基づきビームパターンの最大値を求め、その回転角を電波到来角として求める（θAとθBを端点として（θAとθBの間隔）、最大値をＶ１Ｃとする左右対称な関数式をテーブルから求める）。
以上の検討結果を踏まえると、水中から発信される電波の到来角度を高精度で推定する方法のアルゴリズムは以下のようになる。
受信局において、
（１）予め、遮蔽物のない条件での受信ビームパターンをモデルビームパターンとして取り込んでおく。
（２）次に、実測データを取り込み、所定受信期間内の全受信データの内から最大値を検索する。
（３）検索した最大値を含む放物線の推定ビームパターンをモデルビームパターンと対比して特定する。
（４）（３）の推定ビームパターンを成す放物線の末端に近い受信データから、近い受信データ間のｎ階差（ｎ次の移動平均、但しｎは任意の自然数）をとり、局所トレンドを推定する。
（５）（４）の局所トレンドの推定処理の結果、推定ビームパターンの正常範囲内の値となる受信データを候補に設定する。
（６）（５）で抽出した候補の内で、推定ビームパターンに近い受信データ（例えば、最大値から±２０°〜４０°の中のサンプル）を推定母集団の大勢を決めるサンプルと仮定し、そのサンプルの中から関数を生成する。
（７）（６）で生成した関数を最も可能性があるビームパターンの関数とし、そのビームパターンの関数の頂点を電波到来角として抽出する。
受信局において、
（１）予め、遮蔽物のない条件での受信ビームパターンをモデルビームパターンとして取り込んでおく。
（２）次に、実測データを取り込み、所定受信期間内の全受信データの内から最大値を検索する。
（３）検索した最大値を含む放物線の推定ビームパターンをモデルビームパターンと対比して特定する。
（４）（３）の推定ビームパターンを成す放物線の末端に近い受信データから、近い受信データ間のｎ階差（ｎ次の移動平均、但しｎは任意の自然数）をとり、局所トレンドを推定する。
（５）（４）の局所トレンドの推定処理の結果、推定ビームパターンの正常範囲内の値となる受信データを候補に設定する。
（６）（５）で抽出した候補の内で、推定ビームパターンに近い受信データ（例えば、最大値から±２０°〜４０°の中のサンプル）を推定母集団の大勢を決めるサンプルと仮定し、そのサンプルの中から関数を生成する。
（７）（６）で生成した関数を最も可能性があるビームパターンの関数とし、そのビームパターンの関数の頂点を電波到来角として抽出する。
以上述べた水中から空中へ発信される電波の電波伝搬特性を表す電波伝搬モデルを作成する方法、およびそれを用いた電波到来角度推定方法をそれぞれ方式として構成した実施例を次に説明する。
図４は、本発明のＡＴＳ（Ａｄｖａｎｃｅｄ ｔｅｌｅｍｅｔｒｙ ｓｙｓｔｅｍ）の説明図である。
ＡＴＳ１は、発信器を装着した河川内の魚類の位置を追跡し、魚類生体を観測するシステムである。
河川の流れに沿って、上流から下流に向かって、受信局２（クライアント）Ｎｏ．１〜Ｎｏ．３を適宜間隔で配置する。また、受信局２から適当に離れた位置に制御局３（ホスト）を設け、受信局２Ｎｏ．１〜Ｎｏ．３からの送信信号を受信するようになっている。
図４のシステムは、以下の手順で動作を実行する。
（１）複数の受信局２（実施例の場合、Ｎｏ．１〜Ｎｏ．３）を調査地内に電波受信上、効率的に配置する。
（２）制御局３から各受信局２へ測定命令を無線ＬＡＮを通じて送信し、各受信局２は指向性アンテナ（図５参照）を回転させて電波が最も強い方向を探査する。
（３）電波が最も強い方向を電波発信機４がある角度（電波到来角度）として計測する。
（４）各受信局２は、電波到来角度計測後、制御局３へ電波到来角度を返信する。
（５）制御局３は、三角測量の原理により電波発信機４の位置（魚類の位置）を特定する。
ＡＴＳ１は、発信器を装着した河川内の魚類の位置を追跡し、魚類生体を観測するシステムである。
河川の流れに沿って、上流から下流に向かって、受信局２（クライアント）Ｎｏ．１〜Ｎｏ．３を適宜間隔で配置する。また、受信局２から適当に離れた位置に制御局３（ホスト）を設け、受信局２Ｎｏ．１〜Ｎｏ．３からの送信信号を受信するようになっている。
図４のシステムは、以下の手順で動作を実行する。
（１）複数の受信局２（実施例の場合、Ｎｏ．１〜Ｎｏ．３）を調査地内に電波受信上、効率的に配置する。
（２）制御局３から各受信局２へ測定命令を無線ＬＡＮを通じて送信し、各受信局２は指向性アンテナ（図５参照）を回転させて電波が最も強い方向を探査する。
（３）電波が最も強い方向を電波発信機４がある角度（電波到来角度）として計測する。
（４）各受信局２は、電波到来角度計測後、制御局３へ電波到来角度を返信する。
（５）制御局３は、三角測量の原理により電波発信機４の位置（魚類の位置）を特定する。
図５は、本発明の受信局および制御局の構成図である。
受信局２は、指向性の強いアンテナ、例えば八木アンテナ５をローテータ６で回転駆動し、八木アンテナ５の受信波を受信機７およびＡＤボード８を介して受信電界強度Ｖとして解析用演算装置、例えば解析用パーソナルコンピュータ９に取り込むと共に、ローテータ６の回転を制御するリレー制御ボード１０から出力されるローテータ６の回転角θの値を解析用演算装置、例えば解析用パーソナルコンピュータ９に取り込む。また、受信局２は制御局３との間で無線ＬＡＮを構成する通信手段（図示省略）を備える。
解析用パーソナルコンピュータ９は、以下の処理を行う。
（１）受信電界強度Ｖと回転角度θの抽出を行う。
（２）抽出した受信電界強度Ｖと回転角度θの散布図（離散ビームパターン）を作成する。
（３）離散ビームパターンの数値補間を行い、連続ビームパターンになおし、このパターンの関数式を求める。
（４）前記関数式を解析し、電波到来角度θの算出を行う。
一方、制御局３は、演算装置、例えばホストパーソナルコンピュータ１１と、受信局２との間で無線ＬＡＮを構成する通信手段（図示省略）を備える。
ホストパーソナルコンピュータ１１は、以下の処理を行う。
（１）受信局２の制御、計測スタート制御、アンテナ回転角度制御を行う。
（２）受信局分の電波到来角度θの受信処理を行う。
（３）受信局分の電波到来角度θから三角測量の原理で発信源の位置を特定する。
以下、受信局２の動作は解析用演算装置９が行い、制御局３の動作は演算装置１１が行う。
受信局２は、指向性の強いアンテナ、例えば八木アンテナ５をローテータ６で回転駆動し、八木アンテナ５の受信波を受信機７およびＡＤボード８を介して受信電界強度Ｖとして解析用演算装置、例えば解析用パーソナルコンピュータ９に取り込むと共に、ローテータ６の回転を制御するリレー制御ボード１０から出力されるローテータ６の回転角θの値を解析用演算装置、例えば解析用パーソナルコンピュータ９に取り込む。また、受信局２は制御局３との間で無線ＬＡＮを構成する通信手段（図示省略）を備える。
解析用パーソナルコンピュータ９は、以下の処理を行う。
（１）受信電界強度Ｖと回転角度θの抽出を行う。
（２）抽出した受信電界強度Ｖと回転角度θの散布図（離散ビームパターン）を作成する。
（３）離散ビームパターンの数値補間を行い、連続ビームパターンになおし、このパターンの関数式を求める。
（４）前記関数式を解析し、電波到来角度θの算出を行う。
一方、制御局３は、演算装置、例えばホストパーソナルコンピュータ１１と、受信局２との間で無線ＬＡＮを構成する通信手段（図示省略）を備える。
ホストパーソナルコンピュータ１１は、以下の処理を行う。
（１）受信局２の制御、計測スタート制御、アンテナ回転角度制御を行う。
（２）受信局分の電波到来角度θの受信処理を行う。
（３）受信局分の電波到来角度θから三角測量の原理で発信源の位置を特定する。
以下、受信局２の動作は解析用演算装置９が行い、制御局３の動作は演算装置１１が行う。
前記受信局２の位置決めは以下のようにして決められる。
図６は本発明の魚類行動追跡システムにおけるシステム設置フロー説明図である。
図６（ａ）は調査地域を選定する説明図、図６（ｂ）は仮想上でシステムサービスエリアの選定を説明する図、図６（ｃ）は仮想上のサービスエリアにおける電波伝搬シミュレーションの説明図、図６（ｄ）は仮想上の受信局設置候補エリアの抽出を説明する説明図、図６（ｅ）は受信局候補エリアからの受信エリア解析の説明図、図６（ｆ）は最終的に受信局を設置する説明図である。
魚類行動追跡システムにおけるシステム設置フローは、
（１）調査対象地域を選定する（ＳＴＥＰ（ステップ）１：図６（ａ））。
（２）システムサービスエリアを選定する（仮想（シミュレーション）上でシステムサービスエリア（システムが調査できる範囲）の選定を行う（ＳＴＥＰ１：図６（ｂ））。
（３）サービスエリアから電波伝搬シミュレーションを行う（仮想上で電波発信機をサービスエリア中におき、伝搬可能性を検証する）（ＳＴＥＰ１：図６（ｃ））。
（４）サービスエリアから電波伝搬シミュレーションを行う（電波が安定して受信できる受信局設置候補エリアを選定する（ＳＴＥＰ１：図６（ｄ）））。
（５）受信局設置候補エリアからの受信エリア解析を行う（受信局設置候補エリアからサービスエリアを高確率で見通し、三角測量上適切な箇所の判断を行う。その際、水中からの電波伝搬条件と位置特定測量の関係から受信局の位置を判断する）（ＳＴＥＰ（ステップ）２：図６（ｅ））。
（６）最終的に受信局１、２の箇所に受信局を設置する（結果に基づいて最終的に受信局位置を決定する（ＳＴＥＰ（ステップ）３：図６（ｆ）））。
図６は本発明の魚類行動追跡システムにおけるシステム設置フロー説明図である。
図６（ａ）は調査地域を選定する説明図、図６（ｂ）は仮想上でシステムサービスエリアの選定を説明する図、図６（ｃ）は仮想上のサービスエリアにおける電波伝搬シミュレーションの説明図、図６（ｄ）は仮想上の受信局設置候補エリアの抽出を説明する説明図、図６（ｅ）は受信局候補エリアからの受信エリア解析の説明図、図６（ｆ）は最終的に受信局を設置する説明図である。
魚類行動追跡システムにおけるシステム設置フローは、
（１）調査対象地域を選定する（ＳＴＥＰ（ステップ）１：図６（ａ））。
（２）システムサービスエリアを選定する（仮想（シミュレーション）上でシステムサービスエリア（システムが調査できる範囲）の選定を行う（ＳＴＥＰ１：図６（ｂ））。
（３）サービスエリアから電波伝搬シミュレーションを行う（仮想上で電波発信機をサービスエリア中におき、伝搬可能性を検証する）（ＳＴＥＰ１：図６（ｃ））。
（４）サービスエリアから電波伝搬シミュレーションを行う（電波が安定して受信できる受信局設置候補エリアを選定する（ＳＴＥＰ１：図６（ｄ）））。
（５）受信局設置候補エリアからの受信エリア解析を行う（受信局設置候補エリアからサービスエリアを高確率で見通し、三角測量上適切な箇所の判断を行う。その際、水中からの電波伝搬条件と位置特定測量の関係から受信局の位置を判断する）（ＳＴＥＰ（ステップ）２：図６（ｅ））。
（６）最終的に受信局１、２の箇所に受信局を設置する（結果に基づいて最終的に受信局位置を決定する（ＳＴＥＰ（ステップ）３：図６（ｆ）））。
このシステム設置フローにおける受信局設置候補エリアからの受信エリア解析の際には、水中からの電波伝搬条件と位置特定測量の関係から次の様にして解析する。
（１）水中の発信機４から空中の受信局２への電波伝搬特性は、以下のようにして求める。
まず、本発明の電波伝搬特性を求めるにつき、参考となる従来例を２つ検討する。
１つ目である、電波伝搬モデルの基本形Ａは、以下の既応研究による平面大地反射モデルとなる。
伝搬損失 ＝２０ｌｏｇ（Ａｔ） 伝搬経路の状況により、係数を変更するものとする。
ただし、
ｄ１：直接波の伝搬距離［ｍ］
ｄ２：反射波の伝搬距離［ｍ］ （Δｄ＝ｄｒｉ―ｄＤ）
λ ：波長［ｍ］
Γ ：反射係数（０≦Γ≦１）
２つ目である、森林帯通過に伴う伝搬損の基本形Ｂは、以下のようになる
「無線通信の電波伝搬（進士昌明編著 社団法人電子情報通信学会 Ｐ．４９）」によれば、送受信アンテナの一方が奥行４００ｍ以下の樹林のごく近くにあって、伝搬路の大部分が樹林である場合、樹林による損失Ｌ（ｄＢ）は次式で与えられるとしている。
（２００ < ｆ < ９５，０００ＭＨｚ）
ここで、
ΔＬ ：樹林による損失［ｄＢ］
ｆ ：周波数［ＭＨｚ］
ｄ ：樹林の深さ［ｍ］
（１）水中の発信機４から空中の受信局２への電波伝搬特性は、以下のようにして求める。
まず、本発明の電波伝搬特性を求めるにつき、参考となる従来例を２つ検討する。
１つ目である、電波伝搬モデルの基本形Ａは、以下の既応研究による平面大地反射モデルとなる。
伝搬損失 ＝２０ｌｏｇ（Ａｔ） 伝搬経路の状況により、係数を変更するものとする。
ただし、
ｄ１：直接波の伝搬距離［ｍ］
ｄ２：反射波の伝搬距離［ｍ］ （Δｄ＝ｄｒｉ―ｄＤ）
λ ：波長［ｍ］
Γ ：反射係数（０≦Γ≦１）
２つ目である、森林帯通過に伴う伝搬損の基本形Ｂは、以下のようになる
「無線通信の電波伝搬（進士昌明編著 社団法人電子情報通信学会 Ｐ．４９）」によれば、送受信アンテナの一方が奥行４００ｍ以下の樹林のごく近くにあって、伝搬路の大部分が樹林である場合、樹林による損失Ｌ（ｄＢ）は次式で与えられるとしている。
ここで、
ΔＬ ：樹林による損失［ｄＢ］
ｆ ：周波数［ＭＨｚ］
ｄ ：樹林の深さ［ｍ］
上記基本形Ａおよび基本形Ｂにおいて、直接波の伝搬距離ｄ１［ｍ］、反射波の伝搬距離ｄ２［ｍ］や樹林の深さｄ［ｍ］は、ＧＩＳデータで提供できるので、上記基本形Ａおよび基本形ＢはこれらＧＩＳデータにより修正することが可能となる。
そこで、本発明の水中の発信局から空中の受信局への電波伝搬特性は、これら従来例をベースにして特定の電波伝搬モデルを作成する。これにより、従来例と同様に、ＧＩＳデータにより電波伝搬特性を修正することが可能となる。
また、本発明の水中の発信局から空中の受信局への電波伝搬特性は、直接波の伝搬距離ｄ１［ｍ］、反射波の伝搬距離ｄ２［ｍ］や樹林の深さｄ［ｍ］等のＧＩＳデータや、伝搬媒質の相違による反射係数の変化等の電波伝搬特性におけるＧＩＳデータや反射係数等の変化に応じて求められることになるため、受信局の位置を変えると、電波伝搬特性に影響がでることになる。
そこで、本発明の水中の発信局から空中の受信局への電波伝搬特性は、これら従来例をベースにして特定の電波伝搬モデルを作成する。これにより、従来例と同様に、ＧＩＳデータにより電波伝搬特性を修正することが可能となる。
また、本発明の水中の発信局から空中の受信局への電波伝搬特性は、直接波の伝搬距離ｄ１［ｍ］、反射波の伝搬距離ｄ２［ｍ］や樹林の深さｄ［ｍ］等のＧＩＳデータや、伝搬媒質の相違による反射係数の変化等の電波伝搬特性におけるＧＩＳデータや反射係数等の変化に応じて求められることになるため、受信局の位置を変えると、電波伝搬特性に影響がでることになる。
隣接受信局のカバーエリアは、相互に重複する領域がサービスエリアとして存在し、このサービスエリア内の発信源を隣接受信局からの三角測量により特定する。
具体的には、受信局と制御局を備えた信号送受信方式において、
前記受信局各局により前記電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を受信し、受信信号を制御局へ送信し、
制御局は前記受信信号を取り込み、効率よく受信できる地点を演算し、
受信局設置候補地点におけるカバーエリアを特定し、
カバーエリア同士の重なり合いを判断し、目的とするサービスエリアを確保するための受信局の最適配置を求める。
具体的には、受信局と制御局を備えた信号送受信方式において、
前記受信局各局により前記電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を受信し、受信信号を制御局へ送信し、
制御局は前記受信信号を取り込み、効率よく受信できる地点を演算し、
受信局設置候補地点におけるカバーエリアを特定し、
カバーエリア同士の重なり合いを判断し、目的とするサービスエリアを確保するための受信局の最適配置を求める。
水中から発信される電波の到来角度を高精度で推定する方法は以下の手順で行う。
（１）予め、遮蔽物のない条件での受信ビームパターンをモデルビームパターンとして取り込んでおく（ステップ１）。
モデルビームパターン：水中からの電波は、電波発信源の方向（電波発信源と受信点を直線で結んだ方向）が他の伝搬路（サイドローブ）と比較して極端に強くサンプリング点内の最大値の近傍に求める関数の最大値のある場合がほとんどである。このような特徴を有するモデルビームパターンを、設定条件を替えて複数取り込んでおく。
（１）予め、遮蔽物のない条件での受信ビームパターンをモデルビームパターンとして取り込んでおく（ステップ１）。
モデルビームパターン：水中からの電波は、電波発信源の方向（電波発信源と受信点を直線で結んだ方向）が他の伝搬路（サイドローブ）と比較して極端に強くサンプリング点内の最大値の近傍に求める関数の最大値のある場合がほとんどである。このような特徴を有するモデルビームパターンを、設定条件を替えて複数取り込んでおく。
（２）次に、実測し、所定受信期間内の全受信データの内から最大値を検索する。
最大値：予め、信号とノイズを定義する方法を作成しておく。例えば、１回の計測サンプリングにおいて得られた受信データの最大値を「１」、受信データの最小値を「０」として各サンプリングデータを標準化する。標準化したデータの内、最小値から約２０％をノイズレベルとする。サンプリングデータの最大値から標準化する方法で定義したノイズレベルまで最小値を検索する方法で放物線の末端を検索する。
（３）検索した最大値を含む放物線の推定ビームパターンをモデルビームパターンと対比して特定する。
最大値：予め、信号とノイズを定義する方法を作成しておく。例えば、１回の計測サンプリングにおいて得られた受信データの最大値を「１」、受信データの最小値を「０」として各サンプリングデータを標準化する。標準化したデータの内、最小値から約２０％をノイズレベルとする。サンプリングデータの最大値から標準化する方法で定義したノイズレベルまで最小値を検索する方法で放物線の末端を検索する。
（３）検索した最大値を含む放物線の推定ビームパターンをモデルビームパターンと対比して特定する。
（４）（３）の推定ビームパターンを成す放物線の末端に近い受信データから、近い受信データ間のｍ階差（ｍ次の移動平均、但しｍは任意の自然数）をとり、局所トレンドを推定する。
末端： 上記（２）で説明したように、サンプリングデータの最大値から標準化する方法で定義したノイズレベルまで最小値を検索する方法で放物線の末端を検索する。
移動平均：ノイズの幅よりは、指向性アンテナ（実施例では、八木アンテナ）のアンテナ回転速度と電波発信機の発信間隔の関係から移動平均の次数を決める。この関係に着目することで離散化の度合いを評価することができる。例えば、アンテナの回転速度に対して電波発信機の発信間隔が短い場合、離散化される率は少なくなる。従って、移動平均の次数を上げなくてもスムーズ（連続的）なビームパターンを復元することが可能になり、結果として電波到来角度の推定精度は向上する。
局所トレンド：局所トレンドは階差で近似される関数で、階差を算出した地点から階差の傾きで線を延長（延線）して評価する。階差は１つの算出値となる（例えば、１次直線の傾き）。階差がプラスのとき、関数は増加傾向になる。関数がマイナスのとき、関数は減少傾向になる。この階差は、曲線変化トレンドを抽出するのに優れた指標で、この値が一致する場合には曲線変化は同一の関数形を持つ可能性が高くなる。
末端： 上記（２）で説明したように、サンプリングデータの最大値から標準化する方法で定義したノイズレベルまで最小値を検索する方法で放物線の末端を検索する。
移動平均：ノイズの幅よりは、指向性アンテナ（実施例では、八木アンテナ）のアンテナ回転速度と電波発信機の発信間隔の関係から移動平均の次数を決める。この関係に着目することで離散化の度合いを評価することができる。例えば、アンテナの回転速度に対して電波発信機の発信間隔が短い場合、離散化される率は少なくなる。従って、移動平均の次数を上げなくてもスムーズ（連続的）なビームパターンを復元することが可能になり、結果として電波到来角度の推定精度は向上する。
局所トレンド：局所トレンドは階差で近似される関数で、階差を算出した地点から階差の傾きで線を延長（延線）して評価する。階差は１つの算出値となる（例えば、１次直線の傾き）。階差がプラスのとき、関数は増加傾向になる。関数がマイナスのとき、関数は減少傾向になる。この階差は、曲線変化トレンドを抽出するのに優れた指標で、この値が一致する場合には曲線変化は同一の関数形を持つ可能性が高くなる。
（５）（４）の局所トレンドの推定処理の結果、推定ビームパターンの正常範囲内の値となる受信データを候補に設定する。
正常範囲内：正常範囲内又は正常範囲外を容易に判断するために、局所トレンドに幅を付けている。階差の算出においてできた延線から外れるサンプルを特定し延線とサンプルの距離を算出し、距離のデータサンプルの大きさから一定のマージンを設ける。
正常範囲内：正常範囲内又は正常範囲外を容易に判断するために、局所トレンドに幅を付けている。階差の算出においてできた延線から外れるサンプルを特定し延線とサンプルの距離を算出し、距離のデータサンプルの大きさから一定のマージンを設ける。
（６）（５）で抽出した候補の内で、推定ビームパターンに近い受信データを推定母集団の大勢を決めるサンプルと仮定し、そのサンプルから関数を生成する。
延線との間隔で定義する。延線と一致することは、受信電界強度が最大付近のビームパターンを形成する潜在的な関数の変形であることを意味する。この関数は、現在、２次関数（Ｖ＝ａθ２＋ｂθ＋ｃ）で近似する。
延線との間隔で定義する。延線と一致することは、受信電界強度が最大付近のビームパターンを形成する潜在的な関数の変形であることを意味する。この関数は、現在、２次関数（Ｖ＝ａθ２＋ｂθ＋ｃ）で近似する。
（７）（６）で生成した関数を最も可能性があるビームパターンの関数とし、そのビームパターンの関数の頂点を電波到来角θとして抽出する。
最も可能性があるビームパターンの関数：現在の段階では、２次関数（ｙ＝ａｘ２＋ｂｘ＋ｃ）で近似する。
最も可能性があるビームパターンの関数：現在の段階では、２次関数（ｙ＝ａｘ２＋ｂｘ＋ｃ）で近似する。
以上説明した各局の制御動作は、最終的に所期の動作を行う限りにおいて、手順の入れ替えや代替手順で実行することができる。
１ ＡＴＳ
２ 受信局
３ 制御局
４ 発信機
５ 八木アンテナ
６ ローテータ
７ 受信機
８ ＡＤボード
９ 解析用演算装置（解析用パーソナルコンピュータ）
１０ リレー制御ボード
１１ 演算装置（ホストパーソナルコンピュータ）
２ 受信局
３ 制御局
４ 発信機
５ 八木アンテナ
６ ローテータ
７ 受信機
８ ＡＤボード
９ 解析用演算装置（解析用パーソナルコンピュータ）
１０ リレー制御ボード
１１ 演算装置（ホストパーソナルコンピュータ）
Claims (4)
- 予め、遮蔽物のない条件での受信ビームパターンをモデルビームパターンとして取り込んでおく手段１、
次に、実測データを取り込み、所定受信期間内の全受信データの内から最大値を検索する手段２、
検索した最大値を含む放物線の推定ビームパターンをモデルビームパターンと対比して特定する手段３、
前記手段３における推定ビームパターンを成す放物線の末端に近い受信データから、近い受信データ間の階差をとり、局所トレンドを推定する手段４、
前記手段４における局所トレンドの推定処理の結果、前記推定ビームパターンの正常範囲内の値となる受信データを候補に設定する手段５、
前記手段５で抽出した候補の内で、前記推定ビームパターンに近い受信データを推定母集団の大勢を決めるサンプルと仮定し、そのサンプルの中から関数を生成する手段６、
前記手段６で生成した関数を最も可能性があるビームパターンの関数とし、そのビームパターンの関数の頂点を電波到来角として抽出する手段７、
を有することを特徴とする受信局。 - 請求項１記載の受信局を用い、魚類に装着した電波発信機の発信電波を受信し、前記受信局と制御局とにより前記電波発信機の位置を追跡する信号送受信方式であって、
前記受信局各局は、前記信号送受信方式により電波伝搬モデルを用いて電波到来角度θを求めて前記制御局へ送信し、前記制御局は、前記受信局から前記電波到来角度θを取り込み、三角計測により、発信源の位置を抽出することを特徴とする信号送受信方式。 - 前記受信局は、受信信号のビームパターンをビーム伝送経路に遮蔽物がある場合のパターン、前記ビーム伝送経路により歪みが生じた場合のパターン、およびその他の場合のパターンのいずれかに判定し、前記各パターンに応じて前記各パターンに補正を加えることにより、修正したビームパターンを求め、該修正したビームパターンから前記電波到来角度θを求めることを特徴とする請求項２記載の信号送受信方式。
- 請求項２又は３記載の受信局と制御局を備えた信号送受信方式であって、
シミュレーションにより、前記受信局各局により前記電波伝搬モデルを用いて目的とするサービスエリアから発信される電波を受信し、受信信号を前記制御局へ送信し、前記制御局は前記受信信号を取り込み、効率よく受信できる地点を演算し、受信局設置候補地点におけるカバーエリアを特定し、前記カバーエリア同士の重なり合いを判断し、目的とするサービスエリアを確保するための前記受信局の最適配置位置を求めることを特徴とする信号送受信方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005327242A JP4798531B2 (ja) | 2005-11-11 | 2005-11-11 | 受信局、それを用いた信号送受信方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005327242A JP4798531B2 (ja) | 2005-11-11 | 2005-11-11 | 受信局、それを用いた信号送受信方式 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2007132838A JP2007132838A (ja) | 2007-05-31 |
| JP4798531B2 true JP4798531B2 (ja) | 2011-10-19 |
Family
ID=38154621
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005327242A Expired - Fee Related JP4798531B2 (ja) | 2005-11-11 | 2005-11-11 | 受信局、それを用いた信号送受信方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4798531B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5555347B1 (ja) * | 2013-03-13 | 2014-07-23 | 中国電力株式会社 | 位置検知システム及び位置検知方法 |
| JP7459457B2 (ja) * | 2019-04-11 | 2024-04-02 | 日本電気株式会社 | 位置推定装置、移動体装置、空間データ提示装置、位置推定方法及びプログラム |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS58180967A (ja) * | 1982-04-16 | 1983-10-22 | Toshiba Corp | 方位測定装置 |
| JP3718342B2 (ja) * | 1998-05-12 | 2005-11-24 | 株式会社アドバンテスト | 広域電波監視方法及び装置 |
| JP4214477B2 (ja) * | 2004-03-26 | 2009-01-28 | 独立行政法人土木研究所 | 信号復元方法、そのためのプログラム、記録媒体、信号復元装置、およびそれらを備えた受信局、受信システム |
-
2005
- 2005-11-11 JP JP2005327242A patent/JP4798531B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2007132838A (ja) | 2007-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Merino-Martínez et al. | A review of acoustic imaging methods using phased microphone arrays: Part of the “Aircraft Noise Generation and Assessment” Special Issue | |
| JP4207081B2 (ja) | 電波伝搬特性推定システム及びその方法並びにプログラム | |
| RU2465616C2 (ru) | Способ и устройство для трилатерации с использованием прогнозирования линий связи в пределах прямой видимости и фильтрации трасс в пределах прямой видимости до проведения измерений | |
| US8380150B2 (en) | Received electric field intensity estimating device and received electric field intensity estimating program | |
| JP4608515B2 (ja) | 受信レベル推定システム | |
| CN102692627B (zh) | 用于处理多普勒测量信号的仪器和方法 | |
| Kidera et al. | A robust and fast imaging algorithm with an envelope of circles for UWB pulse radars | |
| US11073498B2 (en) | Detection system, detection device, and detection method | |
| KR101879627B1 (ko) | 가변 송신 출력 및 지향성 안테나를 이용한 근거리 무선통신 기반의 실시간 측위 시스템 및 방법 | |
| JP4798531B2 (ja) | 受信局、それを用いた信号送受信方式 | |
| CN109715071A (zh) | 用于检测介入工具的装置和方法 | |
| JP2003318811A (ja) | 受信電界強度推定計算装置及び方法並びにプログラム及び記録媒体 | |
| JP2009290333A (ja) | 無線ネットワーク評価方法、無線ネットワーク評価システム、無線ネットワーク監視装置、及び無線ネットワーク評価プログラム、プログラムを記録した記録媒体 | |
| KR101868594B1 (ko) | 레이더 빔을 이용한 지형 가시선 분석 장치 및 장거리 레이더 시스템 | |
| JP5664869B2 (ja) | 測定装置、測定システム、測定方法、及びプログラム | |
| KR101308101B1 (ko) | 저주파/고주파 복합형 지하 탐사 레이더 시스템 및 그 제어 방법 | |
| WO2018142629A1 (ja) | 物体検知装置、物体検知方法、及びコンピュータ読み取り可能な記録媒体 | |
| KR101668531B1 (ko) | 전술항법장치에 대한 분석 시뮬레이션 장치 및 방법 | |
| JP2016153728A (ja) | レーダ信号処理装置及びレーダ装置 | |
| Caceres et al. | WLAN-based real time vehicle locating system | |
| JP7279801B2 (ja) | 電波発信源位置推定システム及び電波発信源位置推定方法 | |
| Kao et al. | Feature extraction from a broadband sonar sensor for mapping structured environments efficiently | |
| JP5991599B2 (ja) | 目標探知装置 | |
| JP2007163271A (ja) | 地中レーダ画像処理法 | |
| JP4376733B2 (ja) | レーダ装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20081022 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100811 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100824 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101023 |
|
| 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: 20110712 |
|
| 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: 20110723 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140812 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |