JP3506604B2 - 航行物体までの距離検出装置および距離検出方法 - Google Patents
航行物体までの距離検出装置および距離検出方法Info
- Publication number
- JP3506604B2 JP3506604B2 JP07817998A JP7817998A JP3506604B2 JP 3506604 B2 JP3506604 B2 JP 3506604B2 JP 07817998 A JP07817998 A JP 07817998A JP 7817998 A JP7817998 A JP 7817998A JP 3506604 B2 JP3506604 B2 JP 3506604B2
- Authority
- JP
- Japan
- Prior art keywords
- sound wave
- wave receiver
- distance
- passage
- point
- 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
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は、航行物体までの距
離検出装置および距離検出方法に関する。
離検出装置および距離検出方法に関する。
【０００２】
【従来の技術】従来、海面（または海中）を航行する船
舶までの距離を検出するものとしては、アクティブ方式
の装置がある。
舶までの距離を検出するものとしては、アクティブ方式
の装置がある。
【０００３】このアクティブ方式は、船舶が航行する海
域の海底に、船舶を検出するための音波発信器および受
信器が設けられた検出装置を配置しておき、音波発信器
から発射（送信）された音波が船舶に反射して戻ってき
た音波を受信器にて受信し、このときの音波の送受信間
隔に海中での音速を掛け算することにより、音波発信器
から船舶までの距離を検出するものであった。
域の海底に、船舶を検出するための音波発信器および受
信器が設けられた検出装置を配置しておき、音波発信器
から発射（送信）された音波が船舶に反射して戻ってき
た音波を受信器にて受信し、このときの音波の送受信間
隔に海中での音速を掛け算することにより、音波発信器
から船舶までの距離を検出するものであった。
【０００４】
【発明が解決しようとする課題】上述したアクティブ方
式の装置によると、自ら音波を発信するとともに、所定
地域の海面に向けて設置する必要があるという問題があ
った。
式の装置によると、自ら音波を発信するとともに、所定
地域の海面に向けて設置する必要があるという問題があ
った。
【０００５】そこで、本発明は、自ら音波を発信するこ
となく、かつ設置を容易に行い得る航行物体までの距離
検出装置および距離検出方法を提供することを目的とす
る。
となく、かつ設置を容易に行い得る航行物体までの距離
検出装置および距離検出方法を提供することを目的とす
る。
【０００６】
【課題を解決するための手段】上記課題を解決するため
に、本発明の航行物体までの距離検出装置は、海面また
は海中を航行する航行物体から発される航走音を海面下
にて受信する音波受信器と、この音波受信器からの音波
信号に１回および２回微分を施して音波受信器に対して
航行物体が最も接近する地点を通過する第１通過時点お
よびこの最接近位置に対して４５度方向の地点を通過す
る第２通過時点を検出する測定位置通過時点検出部と、
上記音波受信器からの音波信号を入力して所定時間毎に
音波信号の周波数スペクトル解析を行う周波数解析部
と、上記測定位置通過時点検出部からの第１および第２
通過時点並びに周波数解析部からの周波数スペクトル解
析結果をそれぞれ入力して、ドップラー効果により生じ
る周波数のずれ量を求めるとともに、このずれ量から第
２通過時点での航行物体の音波受信器方向での速度成分
を求めた後、この速度成分から航行物体の進行方向での
速度を検出する速度計算部と、上記２回微分された音波
波形の第１通過時点における値および上記速度計算部に
より求められた速度を入力して、最接近位置までの距離
を計算する距離計算部とから構成したものである。
に、本発明の航行物体までの距離検出装置は、海面また
は海中を航行する航行物体から発される航走音を海面下
にて受信する音波受信器と、この音波受信器からの音波
信号に１回および２回微分を施して音波受信器に対して
航行物体が最も接近する地点を通過する第１通過時点お
よびこの最接近位置に対して４５度方向の地点を通過す
る第２通過時点を検出する測定位置通過時点検出部と、
上記音波受信器からの音波信号を入力して所定時間毎に
音波信号の周波数スペクトル解析を行う周波数解析部
と、上記測定位置通過時点検出部からの第１および第２
通過時点並びに周波数解析部からの周波数スペクトル解
析結果をそれぞれ入力して、ドップラー効果により生じ
る周波数のずれ量を求めるとともに、このずれ量から第
２通過時点での航行物体の音波受信器方向での速度成分
を求めた後、この速度成分から航行物体の進行方向での
速度を検出する速度計算部と、上記２回微分された音波
波形の第１通過時点における値および上記速度計算部に
より求められた速度を入力して、最接近位置までの距離
を計算する距離計算部とから構成したものである。
【０００７】また、本発明の航行物体までの距離検出方
法は、海面または海中を航行する航行物体から発される
走行音を海面下に設けられた音波受信器により受信し、
この音波受信器からの音波信号に１回および２回微分を
施して、音波受信器に対して航行物体が最も接近する地
点を通過する第１通過時点およびこの最接近位置に対し
て４５度方向の地点を通過する第２通過時点を検出し、
上記音波受信器からの音波信号を入力して所定間隔毎に
周波数スペクトル解析を行い、第１通過時点と第２通過
時点との周波数スペクトル解析結果からドップラー効果
により生じる周波数のずれ量を求め、このずれ量に基づ
き第２通過時点における航行物体の音波受信器方向での
速度成分を求めた後、この速度成分から航行物体の進行
方向での速度を検出し、次に上記２回微分された音波波
形の第１通過時点における２回微分値および上記速度計
算部で求められた速度に基づき、音波受信器から最接近
位置までの距離を検出する方法である。
法は、海面または海中を航行する航行物体から発される
走行音を海面下に設けられた音波受信器により受信し、
この音波受信器からの音波信号に１回および２回微分を
施して、音波受信器に対して航行物体が最も接近する地
点を通過する第１通過時点およびこの最接近位置に対し
て４５度方向の地点を通過する第２通過時点を検出し、
上記音波受信器からの音波信号を入力して所定間隔毎に
周波数スペクトル解析を行い、第１通過時点と第２通過
時点との周波数スペクトル解析結果からドップラー効果
により生じる周波数のずれ量を求め、このずれ量に基づ
き第２通過時点における航行物体の音波受信器方向での
速度成分を求めた後、この速度成分から航行物体の進行
方向での速度を検出し、次に上記２回微分された音波波
形の第１通過時点における２回微分値および上記速度計
算部で求められた速度に基づき、音波受信器から最接近
位置までの距離を検出する方法である。
【０００８】上記の各構成によると、海面または海中を
航行する航行物体からの航走音を音波受信器にて受信す
るとともに、異なる通過地点における両周波数のずれ量
を求め、そしてこのずれ量にドップラー効果の式を適用
することにより、航行物体の速度を検出するとともに、
この速度および第１通過時点における音波信号の２回微
分値を使用して、船舶までの距離を検出するようにした
ので、速度検出のための信号を出力する必要がなく、ま
たその検出方向についても全方位的である。
航行する航行物体からの航走音を音波受信器にて受信す
るとともに、異なる通過地点における両周波数のずれ量
を求め、そしてこのずれ量にドップラー効果の式を適用
することにより、航行物体の速度を検出するとともに、
この速度および第１通過時点における音波信号の２回微
分値を使用して、船舶までの距離を検出するようにした
ので、速度検出のための信号を出力する必要がなく、ま
たその検出方向についても全方位的である。
【０００９】
【発明の実施の形態】以下、本発明の実施の形態におけ
る航行物体までの距離検出装置および距離検出方法を、
図１〜図４に基づき説明する。
る航行物体までの距離検出装置および距離検出方法を、
図１〜図４に基づき説明する。
【００１０】本実施の形態においては、航行物体とし
て、海面（水面も含む）を航行する船舶の場合について
説明する。本実施の形態における距離検出方式は、図１
に示すように、距離検出装置１の内、少なくとも音波受
信器（音響センサーの一例で、より具体的には、音圧セ
ンサーが用いられる）１１を所定海域の海底に設置して
おき、船舶２がこの海域を航行する際に、船舶２から発
される航走音の音波を受信して、航行している船舶２の
速度を検出し、そしてこの速度を所定の計算式に代入し
て、最接近位置における船舶までの距離を検出するもの
である。
て、海面（水面も含む）を航行する船舶の場合について
説明する。本実施の形態における距離検出方式は、図１
に示すように、距離検出装置１の内、少なくとも音波受
信器（音響センサーの一例で、より具体的には、音圧セ
ンサーが用いられる）１１を所定海域の海底に設置して
おき、船舶２がこの海域を航行する際に、船舶２から発
される航走音の音波を受信して、航行している船舶２の
速度を検出し、そしてこの速度を所定の計算式に代入し
て、最接近位置における船舶までの距離を検出するもの
である。
【００１１】まず、距離検出装置の構成および距離検出
の原理について説明する。図２に示すように、この距離
検出装置１は、海底（または、海中）に設置されて海面
を航行する船舶２から発される航走音などの音波を受信
する音波受信器１１と、この音波受信器１１で受信され
た音波信号（音圧信号）を増幅器１２で増幅した後入力
するとともに所定の演算を行い、船舶２の速度およびこ
の速度に基づき船舶２までの距離を検出する演算処理部
（ディジタルシグナルプロセッサー、マイクロコンピュ
ータなどから構成されている）１３と、例えば陸上に設
けられた基地局などのその検出速度を送信する送信部１
４とから構成されている。
の原理について説明する。図２に示すように、この距離
検出装置１は、海底（または、海中）に設置されて海面
を航行する船舶２から発される航走音などの音波を受信
する音波受信器１１と、この音波受信器１１で受信され
た音波信号（音圧信号）を増幅器１２で増幅した後入力
するとともに所定の演算を行い、船舶２の速度およびこ
の速度に基づき船舶２までの距離を検出する演算処理部
（ディジタルシグナルプロセッサー、マイクロコンピュ
ータなどから構成されている）１３と、例えば陸上に設
けられた基地局などのその検出速度を送信する送信部１
４とから構成されている。
【００１２】上記演算処理部１３は、増幅器１２で増幅
された音波信号を入力して、音波受信器１１に対して予
め設定された２個所の測定位置を船舶が通過した時点を
検出する測定位置通過時点検出部２１と、同じく上記増
幅器１２からの音波信号を入力して、音波信号の周波数
スペクトル解析（ＦＥＴ解析）を行う周波数解析部２２
と、上記測定位置通過時点検出部２１からの通過時点お
よび周波数解析部２２からの解析結果を入力して、ドッ
プラー効果を利用して船舶２の速度を計算する速度計算
部２３と、音波受信器１２からの音波信号を２回微分
（この構成については後述する）して得られた値の内、
第１通過時点における２回微分値を検出する微分値検出
部２４と、上記速度計算部２３で求められた速度とこの
微分値検出部２４で求められた２回微分値とに基づき所
定の計算を行い、船舶２までの距離を計算する距離計算
部２５とから構成されている。
された音波信号を入力して、音波受信器１１に対して予
め設定された２個所の測定位置を船舶が通過した時点を
検出する測定位置通過時点検出部２１と、同じく上記増
幅器１２からの音波信号を入力して、音波信号の周波数
スペクトル解析（ＦＥＴ解析）を行う周波数解析部２２
と、上記測定位置通過時点検出部２１からの通過時点お
よび周波数解析部２２からの解析結果を入力して、ドッ
プラー効果を利用して船舶２の速度を計算する速度計算
部２３と、音波受信器１２からの音波信号を２回微分
（この構成については後述する）して得られた値の内、
第１通過時点における２回微分値を検出する微分値検出
部２４と、上記速度計算部２３で求められた速度とこの
微分値検出部２４で求められた２回微分値とに基づき所
定の計算を行い、船舶２までの距離を計算する距離計算
部２５とから構成されている。
【００１３】上記測定位置通過時点検出部２１は、増幅
された音波信号から必要な帯域分を取り出すバンドパス
フィルタ（またはローパスフィルタ）などのフィルタ部
３１と、この帯域で取り出された音波信号から所定の信
号波形を取り出す検波回路部３２と、この取り出された
信号波形を微分する第１微分回路部３３と、この第１微
分回路部３３で微分された信号波形をさらに微分する第
２微分回路部３４と、上記第１微分回路部３３で得られ
た１回微分による信号波形のゼロ点を求める第１通過時
点検出部３５と、上記第２微分回路部３４で得られた２
回微分による信号波形のゼロ点を求める第２通過時点検
出部３５とから構成されている。
された音波信号から必要な帯域分を取り出すバンドパス
フィルタ（またはローパスフィルタ）などのフィルタ部
３１と、この帯域で取り出された音波信号から所定の信
号波形を取り出す検波回路部３２と、この取り出された
信号波形を微分する第１微分回路部３３と、この第１微
分回路部３３で微分された信号波形をさらに微分する第
２微分回路部３４と、上記第１微分回路部３３で得られ
た１回微分による信号波形のゼロ点を求める第１通過時
点検出部３５と、上記第２微分回路部３４で得られた２
回微分による信号波形のゼロ点を求める第２通過時点検
出部３５とから構成されている。
【００１４】ここで、上記１回微分および２回微分によ
り、船舶２の通過時点、すなわち２つの測定対象位置を
決定した理由について説明する。上述したように、船舶
２の速度をドップラー効果を用いて検出するように説明
したが、この場合、まず所定速度で航行する船舶２の音
波信号を異なる２個所で測定して得られた両周波数のず
れ量（以下、シフト量という）を求め、次にこのシフト
量をドップラー効果の式に適用して、真の（進行方向で
の）船舶２の速度を求めるのであるが、このドップラー
効果の式を適用する際に、シフトが生じていない基本周
波数を知る必要がある。このため、船舶２が音波受信器
１１に最も接近した位置（最接近位置）Ａを通過した時
点の周波数を検出することにより、基本周波数を知るこ
とができる。
り、船舶２の通過時点、すなわち２つの測定対象位置を
決定した理由について説明する。上述したように、船舶
２の速度をドップラー効果を用いて検出するように説明
したが、この場合、まず所定速度で航行する船舶２の音
波信号を異なる２個所で測定して得られた両周波数のず
れ量（以下、シフト量という）を求め、次にこのシフト
量をドップラー効果の式に適用して、真の（進行方向で
の）船舶２の速度を求めるのであるが、このドップラー
効果の式を適用する際に、シフトが生じていない基本周
波数を知る必要がある。このため、船舶２が音波受信器
１１に最も接近した位置（最接近位置）Ａを通過した時
点の周波数を検出することにより、基本周波数を知るこ
とができる。
【００１５】すなわち、この基本周波数は音波受信器１
１にて測定した音圧レベルが最も高い（強い）時点にお
ける周波数であり、この音圧レベルが最も高い位置は、
音波信号（音圧レベル信号）を１回微分したとき、その
値がゼロとなるときである。
１にて測定した音圧レベルが最も高い（強い）時点にお
ける周波数であり、この音圧レベルが最も高い位置は、
音波信号（音圧レベル信号）を１回微分したとき、その
値がゼロとなるときである。
【００１６】次に、２回微分をする理由であるが、ドッ
プラー効果により求められた船舶２の速度は、船舶２に
おける速度の内、音波受信器１１に対する速度成分であ
り、この速度成分から船舶２の進行方向での真の速度を
求める必要がある。このため、この速度成分の方向と進
行方向との交差角θを知る必要がある（図４参照）。
プラー効果により求められた船舶２の速度は、船舶２に
おける速度の内、音波受信器１１に対する速度成分であ
り、この速度成分から船舶２の進行方向での真の速度を
求める必要がある。このため、この速度成分の方向と進
行方向との交差角θを知る必要がある（図４参照）。
【００１７】この交差角θを、単に、船舶２から受信し
た音波信号から知ろうとすると、２回微分することによ
り、音圧レベルの変化点、すなわち交差角θが４５度で
ある地点を容易に知り得るからである。
た音波信号から知ろうとすると、２回微分することによ
り、音圧レベルの変化点、すなわち交差角θが４５度で
ある地点を容易に知り得るからである。
【００１８】以下、上記事項を具体的に説明する。図１
に示すように、海底に設置された速度検出装置１の音波
受信器１１の位置を、三次元座標系（Ｘ，Ｙ，Ｚ）での
原点（０，０，０）とするとともに、航行する船舶２の
三次元座標を（ｘ，ｙ，ｚ）とすると、船舶２と音波受
信器１１との相対距離（直線距離）ｒは、下記の（１）
式にて表わされる。なお、（１）式を変形して、船舶２
と音波受信器１１との最接近距離ｄで表わすと、（２）
式のようになる。
に示すように、海底に設置された速度検出装置１の音波
受信器１１の位置を、三次元座標系（Ｘ，Ｙ，Ｚ）での
原点（０，０，０）とするとともに、航行する船舶２の
三次元座標を（ｘ，ｙ，ｚ）とすると、船舶２と音波受
信器１１との相対距離（直線距離）ｒは、下記の（１）
式にて表わされる。なお、（１）式を変形して、船舶２
と音波受信器１１との最接近距離ｄで表わすと、（２）
式のようになる。
【００１９】
【数１】
【００２０】また、船舶２から発される音波の受波レベ
ルの距離減衰に関しては、下記の（３）式が成立する。 ＲＬ＝ＳＬ−ＴＬ・・・（３） （３）式中、ＲＬ：音波受信器での受波レベル（ｄＢ） ＳＬ：船舶からの音の送波レベル（ｄＢ） ＴＬ：伝搬損失（ｄＢ）で、ＴＬ＝２０logｒである。
ルの距離減衰に関しては、下記の（３）式が成立する。 ＲＬ＝ＳＬ−ＴＬ・・・（３） （３）式中、ＲＬ：音波受信器での受波レベル（ｄＢ） ＳＬ：船舶からの音の送波レベル（ｄＢ） ＴＬ：伝搬損失（ｄＢ）で、ＴＬ＝２０logｒである。
【００２１】上記（３）式に、ｘ＝ｖｔ［ｖは船舶２の
速度を示し、ｔは船舶２の現在位置（検出位置）Ｂから
水中局１１に対する最接近位置Ａまでの所要時間を示
す。］を代入して変形すると、下記（４）式のようにな
る。
速度を示し、ｔは船舶２の現在位置（検出位置）Ｂから
水中局１１に対する最接近位置Ａまでの所要時間を示
す。］を代入して変形すると、下記（４）式のようにな
る。
【００２２】
【数２】
【００２３】上記（４）式を１回微分すると、下記
（５）式となり、さらに微分して２回微分を行うと、下
記（６）式となる。
（５）式となり、さらに微分して２回微分を行うと、下
記（６）式となる。
【００２４】
【数３】
【００２５】そして、上記２回微分した（６）式の値が
ゼロとなるようにする。すなわち、分子［（ｖｔ）2−
ｄ2］の値がゼロとなる条件は、ｖｔとｄとが等しいと
きであり、このことは、船舶２の現在位置Ｂと船舶２の
最接近位置Ａとの交差角θが４５度であることを示して
いる。
ゼロとなるようにする。すなわち、分子［（ｖｔ）2−
ｄ2］の値がゼロとなる条件は、ｖｔとｄとが等しいと
きであり、このことは、船舶２の現在位置Ｂと船舶２の
最接近位置Ａとの交差角θが４５度であることを示して
いる。
【００２６】次に、上記周波数解析部２２は、周波数ス
ペクトル解析の対象となる帯域（例えば、１０００〜１
０３０Ｈｚ）だけを抽出するための周波数変換部４１
と、この帯域における音波信号をＡ／Ｄ変換するＡ／Ｄ
変換部４２と、このＡ／Ｄ変換部４２で変換された信号
の周波数スペクトル解析（ＦＥＴ解析）を行う解析部４
３と、所定時間毎の解析結果を記録する解析結果記録部
４４とから構成されている。
ペクトル解析の対象となる帯域（例えば、１０００〜１
０３０Ｈｚ）だけを抽出するための周波数変換部４１
と、この帯域における音波信号をＡ／Ｄ変換するＡ／Ｄ
変換部４２と、このＡ／Ｄ変換部４２で変換された信号
の周波数スペクトル解析（ＦＥＴ解析）を行う解析部４
３と、所定時間毎の解析結果を記録する解析結果記録部
４４とから構成されている。
【００２７】上記解析部４３では、所定時間毎に周波数
スペクトル解析が行われ、それぞれの解析結果が解析結
果記録部４４にて記録保持される。そして、速度計算部
２３では、解析結果記録部４４からの解析結果が入力さ
れるとともに、上記各通過時点検出部３５，３６で得ら
れた通過時点すなわち通過時刻が入力されて、図３に示
すように、速度の計算に必要な測定位置（Ａ，Ｂ）に応
じて、すなわち基本周波数を示している基本スペクトラ
ムｆ0および交差角θが４５度の時のシフトスペクトラ
ムｆDに基づき、両測定位置におけるスペクトラムのシ
フト量ＦDが求められる。
スペクトル解析が行われ、それぞれの解析結果が解析結
果記録部４４にて記録保持される。そして、速度計算部
２３では、解析結果記録部４４からの解析結果が入力さ
れるとともに、上記各通過時点検出部３５，３６で得ら
れた通過時点すなわち通過時刻が入力されて、図３に示
すように、速度の計算に必要な測定位置（Ａ，Ｂ）に応
じて、すなわち基本周波数を示している基本スペクトラ
ムｆ0および交差角θが４５度の時のシフトスペクトラ
ムｆDに基づき、両測定位置におけるスペクトラムのシ
フト量ＦDが求められる。
【００２８】以下、シフト量ＦDに基づき船舶２の音波
受信器１１方向での速度ｖrを求める手順を説明する。
ドップラー効果の基本式から、下記（７）式が求まる。
受信器１１方向での速度ｖrを求める手順を説明する。
ドップラー効果の基本式から、下記（７）式が求まる。
【００２９】
ｆD＝｛ｃ／（ｃ−ｖr）｝×ｆ0・・・（７）
但し、ｃは音速を表わす。ここで、シフト量ＦDは下記
（８）式にて表わされる。
（８）式にて表わされる。
【００３０】ＦD＝ｆD−ｆ0・・・（８）
上記（８）式に（７）式を代入すると、シフト量Ｆ
Dは、下記（９）式で表わされる。
Dは、下記（９）式で表わされる。
【００３１】
ＦD＝｛ｃ／（ｃ−ｖr）｝×ｆ0−ｆ0
＝｛ｖr／（ｃ−ｖr）｝×ｆ0・・・（９）
上記（９）式を変形して、ｖrを求めると、下記（１
０）式のようになる。
０）式のようになる。
【００３２】
ｖr＝ＦD×ｃ／（ＦD＋ｆ0）
＝（ＦD×ｃ）／ｆD・・・（１０）
上記（１０）式より、ｖrが求められると、下記（１
１）式により、船舶２の進行方向での速度ｖが求められ
る。
１）式により、船舶２の進行方向での速度ｖが求められ
る。
【００３３】ｖ＝ｖr／ｓｉｎθ・・・（１１）
また、上記微分値検出部２４においては、音波信号を２
回微分した波形信号の値の内、第１通過時点Ａにおける
２回微分値（値が正の場合には最大値、値が負の場合に
は最小値である）（−Ｈ）および速度計算部２３で求め
られた速度ｖに基づき、船舶３までの最接近距離（最短
距離）ｄが計算される。
回微分した波形信号の値の内、第１通過時点Ａにおける
２回微分値（値が正の場合には最大値、値が負の場合に
は最小値である）（−Ｈ）および速度計算部２３で求め
られた速度ｖに基づき、船舶３までの最接近距離（最短
距離）ｄが計算される。
【００３４】以下、この計算手順について説明する。こ
の計算においては、上記受波レベルを示す（３）式を２
回微分して得られた（６）式を使用する。
の計算においては、上記受波レベルを示す（３）式を２
回微分して得られた（６）式を使用する。
【００３５】すなわち、最接近位置Ａにおける船舶２の
位置（ｔ＝０の時）を求めるため、（６）式にｔ＝０を
代入すると、下記（１２）式が得られる。
位置（ｔ＝０の時）を求めるため、（６）式にｔ＝０を
代入すると、下記（１２）式が得られる。
【００３６】
【数４】
【００３７】ところで、ｔ＝０の時の上記（１２）式の
左辺の値は−Ｈであるため、下記（１３）式が成立す
る。（１３）式を変形して、ｄを求めると、下記（１
４）式のようになる。
左辺の値は−Ｈであるため、下記（１３）式が成立す
る。（１３）式を変形して、ｄを求めると、下記（１
４）式のようになる。
【００３８】
【数５】
【００３９】次に、上記距離検出装置により、船舶まで
の最接近距離（最短距離）の検出方法を概略的に説明す
る。まず、音波受信器１１にて海面を航行する船舶２の
航走音の音波を受信する。
の最接近距離（最短距離）の検出方法を概略的に説明す
る。まず、音波受信器１１にて海面を航行する船舶２の
航走音の音波を受信する。
【００４０】この受信された音波信号は増幅器１２で増
幅された後、測定位置通過時点検出部２１に入力され、
ここで１回微分および２回微分が施されて、音波受信器
１１からの最接近位置Ａおよび４５度位置Ｂを通過した
時刻が検出される。
幅された後、測定位置通過時点検出部２１に入力され、
ここで１回微分および２回微分が施されて、音波受信器
１１からの最接近位置Ａおよび４５度位置Ｂを通過した
時刻が検出される。
【００４１】一方、増幅器１２で増幅された音波信号
は、測定位置通過時点検出部２１と並列に、周波数解析
部２２に入力され、ここで所定時間毎の周波数スペクト
ラムが求められる。
は、測定位置通過時点検出部２１と並列に、周波数解析
部２２に入力され、ここで所定時間毎の周波数スペクト
ラムが求められる。
【００４２】そして、測定位置通過時点検出部２１にお
ける第１および第２通過時点検出部３５，３６で得られ
た最接近位置Ａおよび４５度位置Ｂの各通過時刻および
周波数解析部２２での解析結果が速度計算部２３に入力
される。
ける第１および第２通過時点検出部３５，３６で得られ
た最接近位置Ａおよび４５度位置Ｂの各通過時刻および
周波数解析部２２での解析結果が速度計算部２３に入力
される。
【００４３】ここでは、図３に示すように、各位置Ａ，
Ｂでの各通過時刻に対応する周波数スペクトラムを示す
グラフから周波数のシフト量ＦDが求められた後、上記
（１０）式に基づき、音波受信器１１方向での速度成分
ｖrが求められ、さらに（１１）式から、進行方向での
船舶２の真の速度ｖが求められる。
Ｂでの各通過時刻に対応する周波数スペクトラムを示す
グラフから周波数のシフト量ＦDが求められた後、上記
（１０）式に基づき、音波受信器１１方向での速度成分
ｖrが求められ、さらに（１１）式から、進行方向での
船舶２の真の速度ｖが求められる。
【００４４】次に、上記第２微分回路部３４で求められ
た値が微分値検出部２４に入力され、ここで第１通過時
点Ａにおける２回微分値である最小値（−Ｈ）が求めら
れる。
た値が微分値検出部２４に入力され、ここで第１通過時
点Ａにおける２回微分値である最小値（−Ｈ）が求めら
れる。
【００４５】そして、上記船舶２の速度ｖおよび最小値
（−Ｈ）が距離計算部２５に入力されて、船舶２の最接
近位置Ａまでの距離、すなわち最短距離ｄが計算され
る。なお、この求められた最短距離ｄは送信部１４を介
して陸上の基地局に送られ、例えばモニターなどに表示
される。
（−Ｈ）が距離計算部２５に入力されて、船舶２の最接
近位置Ａまでの距離、すなわち最短距離ｄが計算され
る。なお、この求められた最短距離ｄは送信部１４を介
して陸上の基地局に送られ、例えばモニターなどに表示
される。
【００４６】このように、海面を航行する船舶からの航
走音を音波受信器にて受信するとともに、異なる通過地
点における両周波数のシフト量を求め、そしてこのシフ
ト量にドップラー効果の式を適用して、船舶の速度を検
出するとともに、音波信号の２回微分値を使用して、船
舶までの最短距離を検出するようにしたので、従来のよ
うに、距離検出装置側からわざわざ検出用信号を出力す
る必要がないとともに、その検出方向については無指向
性であるため、方位角検出器を使用するものに比べて、
その設置を極めて容易に行うことができる。
走音を音波受信器にて受信するとともに、異なる通過地
点における両周波数のシフト量を求め、そしてこのシフ
ト量にドップラー効果の式を適用して、船舶の速度を検
出するとともに、音波信号の２回微分値を使用して、船
舶までの最短距離を検出するようにしたので、従来のよ
うに、距離検出装置側からわざわざ検出用信号を出力す
る必要がないとともに、その検出方向については無指向
性であるため、方位角検出器を使用するものに比べて、
その設置を極めて容易に行うことができる。
【００４７】ところで、上記実施の形態においては、海
面上を航行する船舶までの距離を検出する場合について
説明したが、海中（水中も含む）を航行する物体例えば
潜水艦までの距離の検出にも適用することができる。
面上を航行する船舶までの距離を検出する場合について
説明したが、海中（水中も含む）を航行する物体例えば
潜水艦までの距離の検出にも適用することができる。
【００４８】
【発明の効果】以上のように本発明の距離検出装置およ
び距離検出方法によると、海面または海中を航行する航
行物体からの航走音を音波受信器にて受信するととも
に、異なる通過地点における両周波数のずれ量を求め、
そしてこのずれ量にドップラー効果の式を適用して航行
物体の速度を検出するとともに、この速度および第１通
過時点における音波信号の２回微分値を使用して、船舶
までの距離を検出するようにしたので、従来のように、
距離検出装置側から検出用信号を出力する必要がないと
ともに、その検出方向については無指向性であるため、
方位角検出器を使用するものに比べて、その設置を極め
て容易に行うことができる。
び距離検出方法によると、海面または海中を航行する航
行物体からの航走音を音波受信器にて受信するととも
に、異なる通過地点における両周波数のずれ量を求め、
そしてこのずれ量にドップラー効果の式を適用して航行
物体の速度を検出するとともに、この速度および第１通
過時点における音波信号の２回微分値を使用して、船舶
までの距離を検出するようにしたので、従来のように、
距離検出装置側から検出用信号を出力する必要がないと
ともに、その検出方向については無指向性であるため、
方位角検出器を使用するものに比べて、その設置を極め
て容易に行うことができる。
【図１】本発明の実施の形態における船舶と音波受信器
との配置状態を示す模式図である。
との配置状態を示す模式図である。
【図２】同実施の形態における距離検出装置の概略構成
を示すブロック図である。
を示すブロック図である。
【図３】同距離検出装置の周波数解析部での周波数スペ
クトル解析結果を示す図である。
クトル解析結果を示す図である。
【図４】同距離検出装置の速度計算部での速度を求める
手順を説明する図である。
手順を説明する図である。
１ 距離検出装置
２ 船舶
１１ 音波受信器
１３ 演算処理部
２１ 測定位置通過時点検出部
２２ 周波数解析部
２３ 速度計算部
２４ 微分最大値検出部
２５ 距離計算部
３３ 第１微分回路部
３４ 第２微分回路部
３５ 第１通過時点検出部
３６ 第２通過時点検出部
４３ 解析部
４４ 解析結果記録部
─────────────────────────────────────────────────────
フロントページの続き
(56)参考文献 特開 平11−271446（ＪＰ，Ａ)
特開 平４−370781（ＪＰ，Ａ)
特開 平１−136081（ＪＰ，Ａ)
特開 平３−33672（ＪＰ，Ａ)
特開 平９−166664（ＪＰ，Ａ)
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01S 15/00 - 15/96
G01S 7/52 - 7/66
Claims (2)
- 【請求項１】海面または海中を航行する航行物体から発
される航走音を海面下にて受信する音波受信器と、この
音波受信器からの音波信号に１回および２回微分を施し
て音波受信器に対して航行物体が最も接近する地点を通
過する第１通過時点およびこの最接近位置に対して４５
度方向の地点を通過する第２通過時点を検出する測定位
置通過時点検出部と、上記音波受信器からの音波信号を
入力して所定時間毎に音波信号の周波数スペクトル解析
を行う周波数解析部と、上記測定位置通過時点検出部か
らの第１および第２通過時点並びに周波数解析部からの
周波数スペクトル解析結果をそれぞれ入力して、ドップ
ラー効果により生じる周波数のずれ量を求めるととも
に、このずれ量から第２通過時点での航行物体の音波受
信器方向での速度成分を求めた後、この速度成分から航
行物体の進行方向での速度を検出する速度計算部と、上
記２回微分された音波波形の第１通過時点における２回
微分値および上記速度計算部により求められた速度を入
力して、最接近位置までの距離を計算する距離計算部と
から構成したことを特徴とする航行物体までの距離検出
装置。 - 【請求項２】海面または海中を航行する航行物体から発
される走行音を海面下に設けられた音波受信器により受
信し、この音波受信器からの音波信号に１回および２回
微分を施して、音波受信器に対して航行物体が最も接近
する地点を通過する第１通過時点およびこの最接近位置
に対して４５度方向の地点を通過する第２通過時点を検
出し、上記音波受信器からの音波信号を入力して所定間
隔毎に周波数スペクトル解析を行い、第１通過時点と第
２通過時点との周波数スペクトル解析結果からドップラ
ー効果により生じる周波数のずれ量を求め、このずれ量
に基づき第２通過時点における航行物体の音波受信器方
向での速度成分を求めた後、この速度成分から航行物体
の進行方向での速度を検出し、次に上記２回微分された
音波波形の第１通過時点における２回微分値および上記
速度計算部で求められた速度に基づき、音波受信器から
最接近位置までの距離を検出することを特徴とする航行
物体までの距離検出方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP07817998A JP3506604B2 (ja) | 1998-03-26 | 1998-03-26 | 航行物体までの距離検出装置および距離検出方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP07817998A JP3506604B2 (ja) | 1998-03-26 | 1998-03-26 | 航行物体までの距離検出装置および距離検出方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11271445A JPH11271445A (ja) | 1999-10-08 |
| JP3506604B2 true JP3506604B2 (ja) | 2004-03-15 |
Family
ID=13654751
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP07817998A Expired - Fee Related JP3506604B2 (ja) | 1998-03-26 | 1998-03-26 | 航行物体までの距離検出装置および距離検出方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3506604B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4958605B2 (ja) * | 2007-04-02 | 2012-06-20 | ユニバーサル特機株式会社 | 移動体位置等推定検出方法、装置及び移動体位置等推定検出方法のプログラム |
| CN101806615B (zh) * | 2010-03-31 | 2012-12-26 | 深圳市南方泰科软件技术有限公司 | 声波测速仪 |
| JP7275472B2 (ja) * | 2018-02-05 | 2023-05-18 | 株式会社Ihi | 速度計測システム |
-
1998
- 1998-03-26 JP JP07817998A patent/JP3506604B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH11271445A (ja) | 1999-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3506604B2 (ja) | 航行物体までの距離検出装置および距離検出方法 | |
| JP3506605B2 (ja) | 航行物体の速度検出装置および速度検出方法 | |
| JP3998363B2 (ja) | 航行物体までの距離検出装置および距離検出方法 | |
| JP3069663B1 (ja) | 音源方向測定方法及び装置 | |
| CN211308915U (zh) | 一种船舶低航速水下检测系统 | |
| JP2916362B2 (ja) | 位置測定における音速補正装置及びその方法 | |
| JP2002071784A (ja) | 航走体放射雑音からの航跡標定方法及び装置 | |
| JP4211218B2 (ja) | 目標判別方法 | |
| JP3199240B2 (ja) | ドップラー補正を行う航走雑音用音響測位装置及び方法 | |
| JP2004012237A (ja) | クロスファンビームによる船舶検出方法及び船舶検出装置 | |
| JPH0850177A (ja) | 超音波式距離測定装置 | |
| JPH11352225A (ja) | 速度測定装置 | |
| JP2820311B2 (ja) | 潜水船の音響航法測定装置 | |
| JPS6050296B2 (ja) | 船舶における水温測定方式 | |
| JPH0850172A (ja) | 音波方位計測装置 | |
| JP2944481B2 (ja) | 水中アクティブ音響探知装置 | |
| JPH035552B2 (ja) | ||
| JPH02203849A (ja) | パルスドプラ計測装置 | |
| CN105807272A (zh) | 一种基于多普勒效应的电动水下航行器航速测量方法及装置 | |
| JP2865090B2 (ja) | パッシブソーナー装置 | |
| JP2976940B2 (ja) | 航走体ケーブルの形状算出方法及びその形状算出装置 | |
| JP2000275322A (ja) | 水中位置測定装置 | |
| JPS5853783A (ja) | 船舶用速度測定装置 | |
| JPH06180358A (ja) | ピンガーの位置検出方法 | |
| JPH0381682A (ja) | 移動式水中目標追跡装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20031118 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20031216 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |