JP2006194627A - 音源位置推定方法及び装置並びにソーナー - Google Patents
音源位置推定方法及び装置並びにソーナー Download PDFInfo
- Publication number
- JP2006194627A JP2006194627A JP2005004243A JP2005004243A JP2006194627A JP 2006194627 A JP2006194627 A JP 2006194627A JP 2005004243 A JP2005004243 A JP 2005004243A JP 2005004243 A JP2005004243 A JP 2005004243A JP 2006194627 A JP2006194627 A JP 2006194627A
- Authority
- JP
- Japan
- Prior art keywords
- sound source
- value
- sound
- arrival time
- virtual
- 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.)
- Pending
Links
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
【課題】 現実に沿わない仮定、条件等に起因する音源位置の推定誤差をなくし、実際の条件等に基づいて音源の位置を推定できる方法等を得る。
【解決手段】 推定海域において複数の仮想音源を配置したときの、海域の環境に応じた、仮想音源から受波地点までの複数の経路による音波の各到来時刻の理論値に基づく値を各仮想音源に対して到来時間差Table作成器3が算出し、観測した音波に基づいて、マルチパス波到来時間差測定器2が音源からの音波に関する値を算出し、音源からの音波に関する値と理論値に基づく値からコスト関数算出器4が算出したコストにより、音源位置推定器5が音源からの音波による到来時刻と一致又は最も近い到来時刻の仮想音源を決定し、仮想音源の配置位置を音源の位置と推定する。あるいは、決定した仮想音源及びその近傍の仮想音源の配置位置とこれらのコストから音源位置を推定する。
【選択図】 図1
【解決手段】 推定海域において複数の仮想音源を配置したときの、海域の環境に応じた、仮想音源から受波地点までの複数の経路による音波の各到来時刻の理論値に基づく値を各仮想音源に対して到来時間差Table作成器3が算出し、観測した音波に基づいて、マルチパス波到来時間差測定器2が音源からの音波に関する値を算出し、音源からの音波に関する値と理論値に基づく値からコスト関数算出器4が算出したコストにより、音源位置推定器5が音源からの音波による到来時刻と一致又は最も近い到来時刻の仮想音源を決定し、仮想音源の配置位置を音源の位置と推定する。あるいは、決定した仮想音源及びその近傍の仮想音源の配置位置とこれらのコストから音源位置を推定する。
【選択図】 図1
Description
本発明は、例えば、水中でパルスを放射する音源の位置を推定する音源位置推定方法及び装置、逆探ソーナー等を含むソーナーに関するものである。
例えば、海中で発せられる音の音源を推定するための方法がある。この方法を行うために逆探ソーナーが用いられる。逆探ソーナーとは、音源から発せられる音(例えばアクティブソーナーが送信する音。以下、探信音という)を信号として受信し、音源の位置を推定するためのパッシブソーナーのことである(非特許文献1参照)。ここで逆探ソーナーは、探信音を受信し、サンプリング等を行って、時系列にディジタルデータ化する音響センサシステムとそのデジタルデータを処理し、音源位置推定を行う音源位置推定装置を備えているものとする。
海中において、音源からの探信音(以下、パルスと呼ぶ)の波(音波)は、例えば海面、海底等での反射等により、直接波以外にも複数の伝搬経路(マルチパスと呼ぶ)を伝搬して、時間差を持って受波センサに到来する。例えば、逆探ソーナーにおいて音源を位置推定する場合、それぞれの経路を伝搬して到来したパルス(マルチパス波と呼ぶ)の到来時間差を観測し、下記の2つの仮定及び1つの条件の下に、音源位置を解析的に推定するという方法がある。
[仮定]
(1)音波は直進する
(2)最初に到来するパルスを直接波とし、2番目に到来するパルスを海底1回反射波とする
[条件]
(1)音源深度及び受波点の深度は既知である
海洋音響用語辞典、海洋音響学会編、p.49、(1999) 海洋音響−基礎と応用−、海洋音響学会編、p.121−122、(1984)
[仮定]
(1)音波は直進する
(2)最初に到来するパルスを直接波とし、2番目に到来するパルスを海底1回反射波とする
[条件]
(1)音源深度及び受波点の深度は既知である
海洋音響用語辞典、海洋音響学会編、p.49、(1999) 海洋音響−基礎と応用−、海洋音響学会編、p.121−122、(1984)
しかしながら、海水温、水圧等の影響で、海中内において音速は空間的に変化するため、音波はスネルの法則に従って屈折しながら伝搬し、受波センサに到来する。つまり、海中では音波は直進しない。さらに、音波が屈折しながら伝搬することにより、直接波が受波点(受波センサ)に到達しない場合がある(つまり、受波センサに最初に到来するパルスが直接波でない可能性がある)。これらの問題があるにも関わらず、従来は、上記の仮定(1)、(2)及び条件(1)の下、パルスの到来時間差を経路長に換算し、音源位置を推定していた。
この場合、実際には音波は直進しないので、仮定(1)、(2)が成り立たず、音源の位置推定に誤差が生じる。また、直接波が到達しない場合、最初に到来した海底M回反射波を直接波とし、次の反射波を海底1回反射波として音源位置を推定してしまうため、音源位置の推定に生じる誤差はさらに大きくなる。さらに、音源深度が未知の場合には、条件(1)を満たさないため、音源位置の推定が困難となってしまう。
そこで、現実に沿わない仮定、条件等に起因する音源位置の推定誤差をなくし、より実際の条件等に基づいて音源の位置を推定できる方法、装置等の実現が望まれていた。
本発明に係る音源位置推定方法は、音源の位置を推定する海域において複数の仮想音源を仮想的に配置したときの、海域の環境に応じた、仮想音源から受波地点までの複数の伝搬経路による音波の各到来時刻の理論値に基づく値を各仮想音源に対して算出する工程と、受波地点において観測した音波に基づいて、推定対象の音源からの音波に関する値を算出する工程と、推定対象の音源からの音波に関する値と理論値に基づく値から、推定対象の音源からの音波による到来時刻と一致又は最も近い到来時刻の仮想音源を決定し、仮想音源の配置位置を推定対象の音源の位置と推定する工程とを有するものである。
本発明においては、音源の位置を推定する海域において複数の仮想音源を仮想的に設定、配置をして、海域の環境に応じた複数の音波の伝搬経路、その経路に基づく受波地点に音波が到来する時刻の理論値に基づく値を算出する。一方、観測によって得られる推定対象の音源からの音波に関する値(例えば、各経路の到来時刻、あらかじめ定めたレプリカの音波との相関関係に基づく値、時系列での音波の大きさ等)を算出する。そして、推定対象の音源からの音波に関する値と理論値に基づく値から、推定対象の音源からの音波による到来時刻と一致又は最も近い到来時刻の仮想音源を決定する。この仮想音源の配置位置を推定対象の音源の位置と推定することにより、あらかじめ設定した仮定、条件の下での音源位置推定ではなく、より現実に沿った音源位置推定を行って、精度を高める。
本発明においては、音源の位置を推定する海域において複数の仮想音源を仮想的に設定、配置をして、海域の環境に応じた複数の音波の伝搬経路、その経路に基づく受波地点に音波が到来する時刻の理論値に基づく値を算出する。一方、観測によって得られる推定対象の音源からの音波に関する値(例えば、各経路の到来時刻、あらかじめ定めたレプリカの音波との相関関係に基づく値、時系列での音波の大きさ等)を算出する。そして、推定対象の音源からの音波に関する値と理論値に基づく値から、推定対象の音源からの音波による到来時刻と一致又は最も近い到来時刻の仮想音源を決定する。この仮想音源の配置位置を推定対象の音源の位置と推定することにより、あらかじめ設定した仮定、条件の下での音源位置推定ではなく、より現実に沿った音源位置推定を行って、精度を高める。
また、本発明に係る音源位置推定装置は、音源の位置を推定する海域において複数の仮想音源を仮想的に配置したときの、海域の環境に応じた、仮想音源から受波地点までの複数の伝搬経路による音波の到来時刻の理論値に基づいて、最初の音波が到来する時刻と他の伝搬経路による音波が到来する時刻との差の理論値を各仮想音源に対して算出し、データとして記憶する到来時間差理論値算出手段と、受波地点において観測した音波のデータに基づいて、推定対象の音源からの最初の音波が到来した時刻と他の経路により到来した時刻との差の観測値を算出する到来時間差観測値算出手段と、推定対象の音源による差の観測値と各仮想音源による差の理論値とにより算出した一致度をコストとして算出するコスト関数算出手段と、コスト関数算出手段が算出したコストに基づいて最も高い一致度が算出された仮想音源の配置位置を推定対象の音源の位置と推定する音源位置推定手段とを備えたものである。
本発明においては、到来時間差理論値算出手段が、音源の位置を推定する海域において複数の仮想音源を仮想的に配置したときの、海域の環境に応じた、仮想音源から受波地点までの複数の伝搬経路による音波の到来時刻の理論値に基づいて、最初の音波が到来する時刻と他の伝搬経路による音波が到来する時刻との差の理論値を各仮想音源に対して算出し、データとして記憶する。一方、到来時間差観測値算出手段が、推定対象の音源からの最初の音波が到来した時刻と他の経路により到来した時刻との差の観測値を算出する。そして、コスト関数算出手段が、差の観測値と各仮想音源による差の理論値とに基づいて、各仮想音源と推定対象の音源との一致度をコストとして算出し、音源位置推定手段がその一致度が最も高い仮想音源の配置位置を推定対象の音源の位置と推定することで、より現実に沿った音源位置推定を行い、精度を高める。
本発明においては、到来時間差理論値算出手段が、音源の位置を推定する海域において複数の仮想音源を仮想的に配置したときの、海域の環境に応じた、仮想音源から受波地点までの複数の伝搬経路による音波の到来時刻の理論値に基づいて、最初の音波が到来する時刻と他の伝搬経路による音波が到来する時刻との差の理論値を各仮想音源に対して算出し、データとして記憶する。一方、到来時間差観測値算出手段が、推定対象の音源からの最初の音波が到来した時刻と他の経路により到来した時刻との差の観測値を算出する。そして、コスト関数算出手段が、差の観測値と各仮想音源による差の理論値とに基づいて、各仮想音源と推定対象の音源との一致度をコストとして算出し、音源位置推定手段がその一致度が最も高い仮想音源の配置位置を推定対象の音源の位置と推定することで、より現実に沿った音源位置推定を行い、精度を高める。
本発明によれば、音源の位置を推定する海域において複数の仮想音源を仮想的に配置して、環境に応じた仮想音源による音波の到来時刻の理論値に基づく値を算出しておき、推定対象の音源からの音波に基づいて得られる値と理論値に基づく値との関係から、複数の仮想音源のうち、推定対象の音源と一致又は最も近い仮想音源を判断し、その仮想音源の配置位置を推定対象の音源の位置と推定するようにしたので、現実的でない仮定にとらわれず、また、音源深度が既知であるという条件を満たさない環境下でも、より現実の環境条件に近く、精度の高い音源位置の推定を行うことができる。
また、本発明によれば、到来時間差理論値算出手段が、音源の位置を推定する海域において配置した各仮想音源に対して差の理論値を算出し、一方で、到来時間差観測値算出手段が差の観測値を算出し、差の理論値及び差の観測値に基づいてコスト関数算出手段が一致度をコストとして算出し、音源位置推定手段がコストに基づいて、一致度が最も高い仮想音源の配置位置を推定対象の音源の位置と推定するようにしたので、より現実に沿った精度の高い音源位置推定を行うことができる装置を得ることができる。
実施の形態1.
図1は本発明の第1の実施の形態に係る音源位置推定装置のブロック図である。本実施の形態及び以下の実施の形態において、音源位置推定装置は音響センサ(図示せず)等と共に、ソーナー(例えば逆探ソーナー。アクティブソーナー又はパッシブソーナーは問わない)に備えられている。図1において、入力端子1には音響センサからディジタルデータを含む信号が入力される。マルチパス波到来時間差測定器2は、ディジタルデータに基づいてマルチパス波の到来時間差観測値をデータとして算出する。到来時間差Table作成器3は設定した各仮想音源に対する到来時間差Table値(理論値)を算出し、到来時間差Table(テーブル形式のデータ)として記憶する。
図1は本発明の第1の実施の形態に係る音源位置推定装置のブロック図である。本実施の形態及び以下の実施の形態において、音源位置推定装置は音響センサ(図示せず)等と共に、ソーナー(例えば逆探ソーナー。アクティブソーナー又はパッシブソーナーは問わない)に備えられている。図1において、入力端子1には音響センサからディジタルデータを含む信号が入力される。マルチパス波到来時間差測定器2は、ディジタルデータに基づいてマルチパス波の到来時間差観測値をデータとして算出する。到来時間差Table作成器3は設定した各仮想音源に対する到来時間差Table値(理論値)を算出し、到来時間差Table(テーブル形式のデータ)として記憶する。
また、コスト関数算出器4は、到来時間差観測値と仮想音源に基づく到来時間差Table値とに基づいて、それらの一致度をコストとして表すためのコスト関数の値をデータとして算出する。音源位置推定器5はコスト関数の値に基づいて音源位置を推定する。そして、出力端子6は推定した音源位置のデータを含む信号を出力するための端子である。
図2はコスト関数算出器4の詳細構成ブロック図である。図2において、入力端子9には到来時間差観測値のデータを含む信号が入力される。一方、入力端子10には到来時間差Table値のデータを含む信号が入力される。Table値−観測値差算出器11は到来時間差観測値と到来時間差Table値との差を算出する。自乗和算出器12は到来時間差観測値と到来時間差Table値との差に基づいてコスト関数の値をデータとして算出する。そして、接続端子13はコスト関数の値のデータを含む信号を音源位置推定器5に送信するための端子である。
ここで、ディジタルデータの処理を行うマルチパス波到来時間差測定器2、到来時間差Table作成器3、コスト関数算出器4(Table値−観測値差算出器11、自乗和算出器12)及び音源位置推定器5は、それぞれコンピュータ等のデータ処理手段(装置)で構成されているものとする。それぞれの手段(又は装置内のある記憶手段)には、後述する処理手順のプログラムが記憶されており、各手段はそのプログラムに基づいて処理を実行する。ここでは、複数の手段がそれぞれの処理を行っているが、また、1つのデータ処理手段が各手段のプログラムを実行して処理を行うようにしてもよい。これは後述する実施の形態における各算出器等においても同様である。
本実施の形態及び以下の実施の形態は、マルチパスにより各経路を経由して到来したパルスの時刻を理論値として算出した値(又はその値に基づく値)を有する複数の仮想音源を音源の位置を推定する海域に配置し、その値と実際に観測された値とに基づいて、推定対象の音源と一致又は最も近い仮想音源を判断することで、その仮想音源の配置位置を推定対象の音源と推定するようにしたものである。
次に本発明の実施の形態1に係る音源位置推定装置の動作について図1、2のブロック図を用い、図3〜7に沿って説明する。まず音源位置推定装置が位置推定処理を実際に行う前の到来時間差Tableの作成手順について説明する。到来時間差Tableの作成に際し、あらかじめ電気伝導度水温水深計CTD(Conductivity-Temperature-Depth profiler )やXBT(Expendable-Bathy-Thermometer)等により、音源位置を推定対象の海域における水温及び塩分濃度を測定し、その測定結果に基づいて音速プロファイル(深度と音速との関係)を算出し、これをデータ化しておく。音速プロファイルの作成については、既存の方法を用いることができる。ここでは、到来時間差Table作成器3は外部の装置によって作成された音速プロファイルのデータを用いるが、到来時間差Table作成器3が測定結果に基づいて音速プロファイルのデータを作成するようにしてもよい。
図3は仮想音源の設定を表す図である。ここでは、図3に示すように、推定対象の海域を仮想的にメッシュ(格子)状に区切り、各交点には仮想音源Sij(i=1,2,…,p、j=1,2,…,q)が存在するものと想定する。本実施の形態では、メッシュ幅(仮想音源の間隔)は特に限定しない。所望する音源位置の推定精度に応じて幅を設定し、仮想的に仮想音源を配置することができる。また、例えば、幅が狭くなるほど、データ量が多くなり、演算量も多くなるので、例えばリアルタイムでの処理を行う必要があれば、幅を広くして演算量を減らす等、状況に応じて設定を行うことができる。ここで、本実施の形態では各仮想音源を等間隔となるように想定しているが、各仮想音源の間隔を一律に想定する必要はない。さらに、本実施の形態では、すべての交点について仮想音源を想定するが、例えば任意の位置(交点)について想定するようにしてもよい。
図4は仮想音源の伝搬経路を示す図である。到来時間差Table作成器3では、まず、音速プロファイルのデータを用いて、図4に示すように、音線理論に基づいて、仮想音源Sijから発せられ、音響センサが受信するまでのマルチパス(音波の伝搬経路)を判断する。これを、想定したすべての仮想音源Sijについて判断する。ここで、音線とは音のエネルギ流れの経路(音の空間において伝搬する様子)を示す線である。
図5は到来時間差Tableを視覚的に表す図である。次に、判断した各仮想音源Sijのマルチパスにおいて、音響センサに到来する時刻が早い順に、第1パルス、第2パルス、…、第Nパルスとし、それらのパルスに基づいて、仮想音源Sijについて、伝搬時間理論値t01(i,j),t02(i,j),…,t0N(i,j)を算出する。さらに、第1パルスと第2パルスの到来時間差τ120 (i,j)、第1パルスと第3パルスの到来時間差τ130 (i,j)、…、第1パルスと第Nパルスの到来時間差τ1N0 (i,j)を、次式(1)に基づいて算出する。これを想定したすべての仮想音源Sijについて行い、仮想音源Sijについて到来時間差Tableを作成する。作成した到来時間差Tableは、到来時間差Table作成器3が有する記憶手段(図示せず)にテーブル形式のデータとして記憶される。
次に、到来(観測)したパルスに基づいて音源位置推定装置が音源の位置を推定する動作について説明する。音源から到来したパルスを音響センサがディジタルデータ化し、そのデータを含む信号が入力端子1を介してマルチパス波到来時間差測定器2に入力される。マルチパス波到来時間差測定器2では、到来したパルスが既知の場合、音響センサからの信号が入力端子1を介して入力されると、既知のパルスのデータS1(k)(kは時刻インデクスを表し、k=1,2,…,Kとする)と音響センサからのデータS2(a)(aは時刻インデクスを表し、a=1,2,…,Aとする。S2(a)は複数のパルスを含むデータである)とを用い、次式(2)に基づいて相関結果C(dt)を算出する。さらに、次式(3)に示すように、相関結果C(dt)の自乗値をCP(dt)として算出する。
そして、CP(dt)の値(ピーク点)を結んでいき、包絡線を求める。包絡線の極大となる箇所を判断し、判断した箇所の時刻t1 ,t2 ,…,tn が、各パルスがそれぞれの経路を伝搬して到来した時刻とする。極大となる箇所を判断する方法については、例えば、包絡線を求めた際に算出した各ピーク点の間の傾きに基づいて、その傾きが正から負に変わるピーク点を極大となる箇所とする。
一方、音源からのパルスが未知で、上記のような相関処理を行おうとする場合、入力端子1を介して音響センサから送信された信号に含まれるデータに基づいて、周波数、パルス長等を判断し、音源からのパルスと同様のパルスをレプリカ信号S1(k)(k=1,2,…,K)として作成する。そして、上述した(2)式及び(3)式に基づいてCP(dt)を算出し、包絡線を求めて、各パルスの到来時刻を決定する。
また、音源からのパルスが未知で、上記のような相関処理を行わずに各パルスの到来時刻を決定する場合、音響センサからのデータS2(a)の自乗値の包絡線を算出する。そして、該包絡線の立ち上がりにおいて、あらかじめ定めた閾値を越える箇所の時刻t1 ,t2 ,…,tn が、各パルスがそれぞれの経路を伝搬して到来した時刻とする。また、包絡線における立ち上がり時刻又は立ち下がり時刻を判断し、その中心の時刻をパルスが到来した時刻としてもよい。
図6は到来時間差観測値を表す図である。本実施の形態では、決定した各時刻に基づいて、マルチパス波到来時間差算出器2が、第1パルスと第2パルスの到来時間差τ12、第1パルスと第3パルスの到来時間差τ13、…、第1パルスと第nパルスの到来時間差τ1nを次式(4)に基づいて、到来時間差観測値のデータとして算出する(b=2,3,…,n)。そして、到来時間差観測値のデータ(n−1データ分。これらをベクトルτob=[τ12,τ13,…,τ1n]とした一群のデータとして扱ってもよい)を含む信号をコスト関数算出器4に送信する。
τ1b=tb −t1 …(4)
τ1b=tb −t1 …(4)
図7は差を算出する到来時間差観測値と到来時間差Table値との組み合わせを表す図である。コスト関数算出器4のTable値−観測値差算出器11には、到来時間差観測値のデータτ12,τ13,…,τ1nを含む信号が接続端子9を介して入力される。また、到来時間差Table作成器3に記憶された、仮想音源Sijにおけるパルスの到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)のデータを含む信号が接続端子10を介して入力される。
そして、次式(5)に基づいて、仮想音源Sijにおける到来時間差Table値と到来時間差観測値との差Δτ12(i,j),Δτ13(i,j),…,Δτ1n(i,j)を算出する(ここで、h=2,…,nである)。これを各仮想音源Sij(i=1,2,…,p、j=1,2,…,q)について行う。本実施の形態では、図7に示すように、τ12とτ120 (i,j)、τ13とτ130 (i,j)、…τ1nとτ1n0 (i,j)の組み合わせで演算を行うものとする(本実施の形態ではn+1以上の到来時間差Table値は用いないものとする)。
Δτ1h(i,j)=τ1h0 (i,j)−τ1h …(5)
Δτ1h(i,j)=τ1h0 (i,j)−τ1h …(5)
自乗和算出器12は、Table値−観測値差算出器11が算出した仮想音源Sijにおける到来時間差Table値と到来時間差観測値との差のデータから、次式(6)に基づいてコスト関数Jijを算出する。そして、各仮想音源Sijについて算出したコスト関数Jijのデータを含む信号を、接続端子13を介して音源位置推定器5に出力する。
音源位置推定器5は、各仮想音源Sijのコスト関数Jijに基づいて、その中で最小値のコスト関数Jijを決定する。そして、最小のコスト関数Jijが算出された仮想音源Sijについて設定した位置を、音源の位置として推定する。ここで、各コスト関数Jijは、到来時間差観測値と各仮想音源の到来時間差Table値との一致度を表す。一致度合いが高いほど、到来時間差観測値と到来時間差Table値との差が小さくなり、コスト関数Jijも小さくなる。したがってコスト関数Jijが最小の仮想音源Sijが、観測したパルスを発する音源に最も近いものと推定できる。
以上のように第1の実施の形態によれば、海中の音速プロファイルを考慮し、音波の屈折現象を反映した音線理論に基づいて仮想的に設定した仮想音源Sijの到来時間差Table値を到来時間差Table作成器3が算出して記憶しておき、その到来時間差Table値とマルチパス波到来時間差測定器2が算出した実際の到来時間差観測値との差を算出し、コスト関数算出器4が第1パルスと第2パルス以降の各パルスとの時間差の一致度となるコスト関数Jijを算出し、さらに、音源位置推定器5が最も一致度が高い仮想音源に設定された位置を、音源位置として判断するようにしたので、上述した仮定を満足しない条件下においても、その仮定にとらわれることなく、精度良く音源位置を推定することができる。また、音源深度が既知である等の条件を設定しなくてもよいので、音源深度が未知(条件を充たさない)である場合でも音源位置を推定することができる。
実施の形態2.
図8は差を算出する到来時間差Table値と到来時間差観測値との関係を表す図である。図8(a)は、到来時間差観測値τ12,τ13,…,τ1nと到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)を示す。また、図8(b)は、到来時間差観測値τ12,τ13,…,τ1n及び到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)の組み合わせを時間軸を用いて表したものである。第2の実施の形態における装置は、Table値−観測値差算出器11において、到来時間差Table値と到来時間差観測値との組み合わせを決定するための処理が異なるものの、その装置構成は第1の実施の形態で説明したものと同じである。そこで、本実施の形態では、装置構成については図1及び図2を用い、処理については図8及び9を参照しながら、本実施の形態における音源位置推定装置について説明することにする。
図8は差を算出する到来時間差Table値と到来時間差観測値との関係を表す図である。図8(a)は、到来時間差観測値τ12,τ13,…,τ1nと到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)を示す。また、図8(b)は、到来時間差観測値τ12,τ13,…,τ1n及び到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)の組み合わせを時間軸を用いて表したものである。第2の実施の形態における装置は、Table値−観測値差算出器11において、到来時間差Table値と到来時間差観測値との組み合わせを決定するための処理が異なるものの、その装置構成は第1の実施の形態で説明したものと同じである。そこで、本実施の形態では、装置構成については図1及び図2を用い、処理については図8及び9を参照しながら、本実施の形態における音源位置推定装置について説明することにする。
例えば受信信号のS/Nが低い場合、(ア)、(イ)のような誤検出が発生することがあり、到来時間差観測値と到来時間差Table値を単純に時間差に沿って選択すると、差を算出するための組み合わせが食い違ってしまうことがある。
(ア)雑音を誤ってパルスとして検出する
(イ)存在するはずのパルスを検出できない
差をとる組み合わせに食い違いが生じた結果、例えば真の音源位置に最も近い仮想音源Sijの到来時間差Table値と到来時間差観測値の間でも、伝搬経路(海面及び海底の反射回数)が異なるものどうしの差を算出することになるため、コスト関数が最小とならない可能性がある。本実施の形態は、Table値−観測値差算出器11において、到来時間差観測値と到来時間差Table値との間で、図8(b)の矢印によって示すように、互いに最も近い値どうしを1対1で対応させ、差を算出するための組み合わせを決定するようにしたものである。そして、組み合わせた到来時間差観測値と到来時間差Table値の差Δτ(x)(x=1,2,…,X)を算出する。Xは到来時間差観測値と到来時間差Table値の組み合わせの数を表す。
(ア)雑音を誤ってパルスとして検出する
(イ)存在するはずのパルスを検出できない
差をとる組み合わせに食い違いが生じた結果、例えば真の音源位置に最も近い仮想音源Sijの到来時間差Table値と到来時間差観測値の間でも、伝搬経路(海面及び海底の反射回数)が異なるものどうしの差を算出することになるため、コスト関数が最小とならない可能性がある。本実施の形態は、Table値−観測値差算出器11において、到来時間差観測値と到来時間差Table値との間で、図8(b)の矢印によって示すように、互いに最も近い値どうしを1対1で対応させ、差を算出するための組み合わせを決定するようにしたものである。そして、組み合わせた到来時間差観測値と到来時間差Table値の差Δτ(x)(x=1,2,…,X)を算出する。Xは到来時間差観測値と到来時間差Table値の組み合わせの数を表す。
図9は、Table値−観測値差算出器11による、組み合わせ決定を含めた処理の一例を表すフローチャートである。ここでは、ある仮想音源Sijの到来時間差Table値と到来時間差観測値との組み合わせを決定する手順についてさらに詳しく説明する。ここで図9においては仮想音源の位置に関する添字(i,j)を省略している。
まず、コスト関数算出器4のTable値−観測値差算出器11には、到来時間差観測値のデータτ12,τ13,…,τ1nを含む信号が接続端子9を介して入力される(ここでは、これらのデータの一群をベクトルτob=[τ12,τ13,…,τ1n]として表す)。また、到来時間差Table作成器3に記憶された、仮想音源Sijにおけるパルスの到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)のデータを含む信号が接続端子10を介して入力される(ここでは、これらのデータの一群をベクトルτTable =[τ120 ,τ130 ,…,τ1N0 ]として表す)。
Table値−観測値差算出器11では、まず、τob(1)(=τ12)について、最も近い到来時間差Table値をτTable の要素の中から決定する。つまり、τob(1)との差の絶対値が最も小さくなるτTable の要素を決定する。そして、決定した要素をτnearとし、決定したτTable の要素のインデクス(例えば要素を表す名前、番号等)をインデクスI(1)にデータとして代入する。さらに、τob(1)と最も近い到来時間差Table値との差の絶対値を算出し、Δτ(1)に代入する。以上を式で表すと次式(7)のようになる(ブロックa)。
|τnear−τob(1)|=Δτ(1) …(7)
|τnear−τob(1)|=Δτ(1) …(7)
次に、インデクスを指定するためのポインタxに2を代入する(ブロックb)。そして、τob(x)(ここではx=2)について、前述したように差の絶対値に基づいて、最も近い到来時間差Table値をτTable の要素の中から決定する。そして、決定した要素をτnearとし、決定したτTable の要素のインデクスをI(x)にデータとして代入する。さらに、到来時間差Table値と到来時間差観測値の差の絶対値を算出し、Δτ(x)に代入する(ブロックc)。以上を式で表すと次式(8)のようになる。
|τnear−τob(x)|=Δτ(x) …(8)
|τnear−τob(x)|=Δτ(x) …(8)
次に、インデクスI(x−1)のデータとI(x)のデータとの関係について、I(x−1)=I(x)であるかどうか、つまり、τob(x−1)とτob(x)とについて、最も近いτTable の要素が共通しているかどうかを判断する(分岐d)。
I(x−1)=I(x)であると判断すると、さらにΔτ(x−1)≧Δτ(x)であるかどうか、つまりτob(x−1)とτob(x)とのどちらがよりτnearに近いかどうかを判断する(分岐e)。Δτ(x−1)≧Δτ(x)であれば(τob(x−1)とτob(x)とが同じ又はτob(x)の方がτnearに近ければ)、τob(x)とτnearとを差をとる組み合わせとして採用することとし、Δτ(x−1)に0を代入する(ブロックf)。一方、Δτ(x−1)≧Δτ(x)でなければ(τob(x−1)の方がτnearに近ければ)、τob(x−1)とτnearとを差をとる組み合わせとして採用することとし、Δτ(x)に0を代入する(ブロックg)。そして、xにx+1を代入してxを1増加させる(ブロックh)。
x>n−1であるかどうか(つまり、ベクトルτobの最後(n−1番目)の要素(到来時間差観測値τ1n)について処理を行ったかどうか、すべての到来時間差観測値について処理を行ったかどうか)を判断する(分岐i)。すべての要素について処理が終了していないと判断すると、次の要素について処理を続ける。すべての要素について処理が終了していると判断すると、ある仮想音源Sijにおける処理を終了する。
以上のように決定した組み合わせによりTable値−観測値差算出器11が算出したΔτ(x)(x=1,2,…,n−1)に基づいて、第1の実施の形態と同様に自乗和算出器12がコスト関数Jijを算出する。以上の処理を、コスト関数算出器4(Table値−観測値差算出器11、自乗和算出器12)は、各仮想音源Sijについて行い、算出したコスト関数Jijのデータを含む信号を、接続端子13を介して音源位置推定器5に出力する。音源位置推定器5での音源位置の推定処理は、第1の実施の形態で説明した処理と同様であるので説明を省略する。
以上のように第2の実施の形態によれば、Table値−観測値差算出器11が、互いに最も近い到来時間差Table値と到来時間差観測値との組み合わせを1対1で決定し、その組み合わせに基づいてコスト関数Jij算出のための差を算出するようにしたので、真の音源の位置の近傍において、海面及び海底における反射回数による経路が同じもの同士の到来時間差Table値と到来時間差観測値との組み合わせの差がコストとして計上されることとなり、誤検出に基づく余分な到来時間差観測値や到来時間差Table値はコスト関数Jijの算出に影響しないため、誤検出に対して、ロバストな音源位置を推定することができる。なお、本実施の形態において説明した図9のフローチャートの処理手順は一例であり、本実施の形態のように、互いに最も近い到来時間差観測値と到来時間差Table値との組み合わせを決定してΔτ(x)が算出できるのであれば処理手順は限定しない。
実施の形態3.
図10は差を算出するTable値と到来時間差観測値との関係を表す図である。図10(a)は、到来時間差観測値τ12,τ13,…,τ1nと到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)を示す。また、図10(b)は、到来時間差観測値τ12,τ13,…,τ1n及び到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)の組み合わせを時間軸を用いて表したものである。第3の実施の形態における装置についても、Table値−観測値差算出器11において、到来時間差Table値と到来時間差観測値との組み合わせを決定する処理が異なるものの、その構成は第1の実施の形態で説明したものと同じであるので、装置については図1及び図2を用いて説明する。
図10は差を算出するTable値と到来時間差観測値との関係を表す図である。図10(a)は、到来時間差観測値τ12,τ13,…,τ1nと到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)を示す。また、図10(b)は、到来時間差観測値τ12,τ13,…,τ1n及び到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)の組み合わせを時間軸を用いて表したものである。第3の実施の形態における装置についても、Table値−観測値差算出器11において、到来時間差Table値と到来時間差観測値との組み合わせを決定する処理が異なるものの、その構成は第1の実施の形態で説明したものと同じであるので、装置については図1及び図2を用いて説明する。
例えば、マルチパスの多い環境において、S/Nが低く、マルチパス波の到来時刻推定精度が高くない場合、真の音源位置に近い仮想音源とそうでない仮想音源とのコスト関数の値において差が強調されず、グローバルミニマムとなる可能性がある。そこで、本実施の形態は、図10(b)に示す矢印Aのように、到来時間差観測値に対し、到来時間差Table値の近い値を組み合わせていき、組み合わせができなかった到来時間差Table値についても、矢印Bのように、さらに最も近い到来時間差観測値と組み合わせることで、到来時間差Table値を有効に利用しようとするものである(図10(b)ではτ13とτ130 (i,j),τ140 (i,j))。したがって、本実施の形態は必ずしも1対1の対応になるとは限らない。ここで、第1及び第2の実施の形態と同様に、到来時間差観測値のうち、最も時間差が大きいτ1nと最も近い到来時間差Table値(図10(b)ではτ1(n-1)0 (i,j))よりも大きな時間差を有する到来時間差Table値は採用しないものとする。
図11はTable値−観測値差算出器11による、組み合わせ決定を含めた処理の一例を表すフローチャートである。ここでは、ある仮想音源Sijの到来時間差Table値と到来時間差観測値との組み合わせを決定する手順について説明する(図11では仮想音源の位置に関する添字(i,j)を省略している)。
コスト関数算出器4のTable値−観測値差算出器11には、到来時間差観測値のデータτ12,τ13,…,τ1nを含む信号が接続端子9を介して入力される(ここでは、これらのデータの一群をベクトルτob=[τ12,τ13,…,τ1n]として表す場合もある)。また、到来時間差Table作成器3に記憶された、仮想音源Sijにおけるパルスの到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)のデータを含む信号が接続端子10を介して入力される。
Table値−観測値差算出器11は、まず、到来時間差観測値のインデクスを指定するポインタxに1を代入する(ブロックj)。そして、τob(x)(ここではx=1)について、最も近い到来時間差Table値をτTable の要素の中から決定する。そして、決定した要素をτnearとし、決定したτTable のインデクスをI(x)に代入する。さらに、到来時間差Table値と到来時間差観測値の差を表すベクトルΔτ(x)を前述した(8)式に基づいて算出する(ブロックk)。
次にxにx+1を代入してxを1増加させる(ブロックl)。そして、ベクトルτobを構成するn−1の要素(到来時間差観測値τ12,τ13,…,τ1n)について処理を行ったかどうかを判断する(分岐m)。すべての要素について処理が終了していないと判断すると、次の要素について処理を続ける。
分岐mにおいて、すべての要素について処理が終了していると判断すると、τTable のインデクスの中での最大値I(n−1)の値をImax とする(ブロックn)。図10(b)では、仮想音源Sijの到来時間差Table値におけるτ1(n-1)0 (i,j)のインデクスとなる。そして、τTable のインデクスを指定するポインタyに1を代入する(ブロックo)。さらにそのポインタyの値がベクトルIの要素となっているかどうかを判断する(分岐p)。
ポインタyの値がベクトルIの要素となっていない、つまり、τTable (y)がベクトルτobとの間で、差を算出する組み合わせに採用されていないと判断すると、最も近い到来時間差観測値をベクトルτobの要素の中から決定し、その要素をτnearとする。そして、到来時間差Table値と到来時間差観測値の差を表すベクトルΔτ(x)を算出する(ブロックq)。以上を式で表すと、次式(9)のようになる。その後、xにx+1を代入してxを1増加させる(ブロックr)。
|τnear−τTable(y)|=Δτ(x) …(9)
|τnear−τTable(y)|=Δτ(x) …(9)
そして、yにy+1を代入してyを1増加させる(ブロックs)。そして、分岐p〜sまでの処理をyがImax よりも大きくなったものと判断するまで繰り返す(分岐t)。yがImax よりも大きくなったものと判断すると、ある仮想音源Sijにおいて、到来時間差観測値と到来時間差Table値との組み合わせを決定する処理を終了する。
以上のように決定した組み合わせによりTable値−観測値差算出器11が算出したΔτ(x)(x=1,2,…,X、Xは組み合わせの数)に基づいて、第1の実施の形態と同様に自乗和算出器12がコスト関数Jijを算出する。以上の処理をコスト関数算出器4(Table値−観測値差算出器11、自乗和算出器12)は仮想音源Sij毎に行い、各仮想音源Sijについて算出したコスト関数Jijのデータを含む信号を、接続端子13を介して音源位置推定器5に出力する。音源位置推定器5での音源位置の推定処理は、第1の実施の形態で説明した処理と同様であるので説明を省略する。
以上のように第3の実施の形態によれば、実際に到来したパルスに基づく到来時間差観測値及び到来時間差観測値の最大値に最も近い到来時間差Table値以下の値を有する到来時間差Table値を余すことなく使用し、コスト関数Jijの値を算出し、音源位置の推定を行うようにしたので、一致度を表すコスト関数Jijの値について、実際の音源位置と近い仮想音源とそうでない仮想音源との差を強調し、真の音源位置近傍の仮想音源Sijにおいて算出したコスト関数Jijの勾配を急峻にすることで、精度のよい推定を行うことができる。さらに、真の音源位置から離れた仮想音源においては、互いにずれが大きく対応しない到来時間差観測値、到来時間差Table値の差に基づく値がコスト関数Jijの値に計上することになり、ローカルミニマムの値を大きくすることができる。これにより、マルチパスが多く、到来時刻推定精度が高くない条件の下でもロバストな音源位置の推定を行うことができる。なお、本実施の形態において説明した図11のフローチャートの処理手順は一例であり、本実施の形態のように、到来時間差観測値及び到来時間差Table値を有効に利用してΔτ(x)が算出できるのであれば処理手順は限定しない。
実施の形態4.
図12は本発明の第4の実施の形態に係る音源位置推定装置のブロック図である。図12において図1と同じ符号を付しているものは、上述した実施の形態と同様の動作を行うので説明を省略する。レプリカ相関器16は、観測した音波に基づく値を算出するデータ算出手段であり、レプリカ信号に含まれる既知のパルスのデータと、入力端子1を介して音響センサから送信された信号に含まれるデータとの間で、相関結果であるC(dt)を算出する。図12では、装置外部(例えば記憶手段等)から送信されたレプリカ信号に基づいて相関結果を算出しているが、例えば、上述したように、パルスが未知の場合、音響センサシステムからの受波したデータに基づいてレプリカ信号とするためのデータを作成することもある。また、一度送信されたレプリカ信号に含まれるデータを記憶手段(図示せず)に記憶させておくこともできる。包絡線算出器17は、観測した音波に基づく値を処理する処理手段であり、相関結果の自乗値CP(dt)を算出し、そのピーク点を結ぶことにより、CP(dt)の包絡線CPEを求める。第1パルス到来時刻推定器18は、包絡線の中で最初に閾値を越え、かつ極大となる時刻を求める。これを第1パルス到来時刻推定値t1とする。また、コスト関数算出器19は、コスト関数Jijの値をデータとして算出するが、その算出手順が、後述するようにコスト関数算出器4とは異なる。
図12は本発明の第4の実施の形態に係る音源位置推定装置のブロック図である。図12において図1と同じ符号を付しているものは、上述した実施の形態と同様の動作を行うので説明を省略する。レプリカ相関器16は、観測した音波に基づく値を算出するデータ算出手段であり、レプリカ信号に含まれる既知のパルスのデータと、入力端子1を介して音響センサから送信された信号に含まれるデータとの間で、相関結果であるC(dt)を算出する。図12では、装置外部(例えば記憶手段等)から送信されたレプリカ信号に基づいて相関結果を算出しているが、例えば、上述したように、パルスが未知の場合、音響センサシステムからの受波したデータに基づいてレプリカ信号とするためのデータを作成することもある。また、一度送信されたレプリカ信号に含まれるデータを記憶手段(図示せず)に記憶させておくこともできる。包絡線算出器17は、観測した音波に基づく値を処理する処理手段であり、相関結果の自乗値CP(dt)を算出し、そのピーク点を結ぶことにより、CP(dt)の包絡線CPEを求める。第1パルス到来時刻推定器18は、包絡線の中で最初に閾値を越え、かつ極大となる時刻を求める。これを第1パルス到来時刻推定値t1とする。また、コスト関数算出器19は、コスト関数Jijの値をデータとして算出するが、その算出手順が、後述するようにコスト関数算出器4とは異なる。
図13は包絡線CPEの一例を表す図である。求めた包絡線においてはそれぞれ1の到来時間差観測値が算出される。ここで複数のパルスが重なっているため、通常分離するはずの包絡線が分離できていない場合、上述の第1〜3の実施の形態では、1の包絡線として扱われ、1つの到来時間差観測値しか算出されないことになる。したがって、真の音源位置又はその近傍の仮想音源の到来時間差Table値と到来時間差観測値との間でも対応関係が異なってしまい、算出したコスト関数Jijの値が大きくなり、最小値にならない場合がある。このような場合、本来とは異なる仮想音源の位置が音源位置として推定されてしまうため、正確な音源位置が推定できていないことになってしまう。そこで、本実施の形態は、各到来時間差Table値に対応する相関結果の自乗値の和をコスト関数Jijの値として算出することにより、分離できていない包絡線においても複数の到来時間差Table値との対応を考慮したコスト関数を算出することができる。
次に本実施の形態に係る音源位置推定装置の動作について、図12の構成に基づいて、図14、15を参照しながら説明する。レプリカ相関器16は、入力端子1を介して音響センサから送信された信号に含まれるデータ及び既知のパルスのデータから、第1の実施の形態で説明した(2)式に基づいてレプリカ相関処理を行い、相関結果C(dt)のデータを含む信号を包絡線算出器17に出力する。包絡線算出器17は、(3)式に基づいて相関結果C(dt)の自乗値をCP(dt)として算出する。そして、CP(dt)におけるピーク点を結ぶことにより、包絡線CPEを求め、そのデータを含む信号を第1パルス到来時刻推定器18及びコスト関数算出器19に出力する。ここで、相関出力の自乗値CP(dt)の全てのピーク点を、さらに2次関数フィッティング等を行って結び、包絡線CPEを求めるようにしてもよい。この場合、包絡線CPEを精度よく求められるため、音源位置推定精度を向上させることができる。
図14は包絡線CPEと第1パルス到来時刻推定値t1との関係を表す図である。第1パルス到来時刻推定器18は、包絡線CPEのデータに基づいて、包絡線CPEの中であらかじめ設定した閾値を越えた部分があるかどうかを判断し、越えた部分に基づいて最初の極大点を検出する。そして、検出した極大点における時刻を第1パルスの到来時刻t1として、そのデータを含む信号をコスト関数算出器19に出力する。
図15は包絡線CPEと到来時間差Table値との関係を表す図である。コスト関数算出器19には、包絡線CPEのデータ及び第1パルスの到来時刻推定値t1を含む信号に加え、到来時間差Table作成器3に記憶された、仮想音源Sijにおけるパルスの到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)のデータを含む信号が入力される。そして、包絡線CPEのデータ、第1パルスの到来時刻推定値t1及びパルスの到来時間差Table値のデータに基づいて、図15に示すように、時刻t1+τ120 (i,j),t1+τ130 (i,j),…,t1+τ1N0 (i,j)における包絡線CPEの値を求める。そして、その総和値の逆数を、仮想音源Sijにおけるコスト関数Jijの値とし、そのデータを含む信号を音源位置推定器5に出力する。以上を式で表すと、次式(10)となる。(10)式における包絡線CPEの値CPE(t1+τ1z0 (i,j))(z=1,2,…,N)は、線形補間等の補間処理によって求める。また、t1+τ1z0 (i,j)に最も近い時刻に算出された包絡線による値を用いるようにしてもよい。ここで、τ110 (i,j)=0とする。音源位置推定器5での音源位置の推定処理は、第1の実施の形態で説明した処理と同様であるので説明を省略する。
以上のように第4の実施の形態によれば、コスト関数算出器19が、パルスの到来時間差Table値に基づく仮想音源Sijにおけるコスト関数の値Jijを相関処理がなされた結果である包絡線CPEから直接的に導き出すようにしたので、パルスが重なって受信されたために、本来分離されているべき包絡線が分離されていない場合でも、それを1つの包絡線として扱ってしまうことがなくなるため、真の音源位置及びその近傍に対応する仮想音源Sijのコスト関数の値Jijが大きくならず、最小となる。そのため、重なって受信したパルスを分離できないことによる誤検出が生じた場合に対してロバストな音源位置を推定することができる。
実施の形態5.
図16は本発明の第5の実施の形態に係る音源位置推定装置のブロック図である。図16において図1及び図12と同じ符号を付しているものは、上述した実施の形態と同様の動作を行うので説明を省略する。第1パルス到来時刻幅算出器20は、第1パルス到来時刻に幅をもたせるための開始時刻tstart と終了時刻tend を決定する。
図16は本発明の第5の実施の形態に係る音源位置推定装置のブロック図である。図16において図1及び図12と同じ符号を付しているものは、上述した実施の形態と同様の動作を行うので説明を省略する。第1パルス到来時刻幅算出器20は、第1パルス到来時刻に幅をもたせるための開始時刻tstart と終了時刻tend を決定する。
図17はコスト関数算出器21の詳細構成ブロック図である。図17において、入力端子22には包絡線算出器17から包絡線のデータを含む信号が入力される。また、入力端子23には第1パルス到来時刻幅算出器20から開始時刻tstart と終了時刻tend のデータを含む信号が入力される。さらに、入力端子25には到来時間差Table値のデータを含む信号が入力される。コストベクトル算出器24は、第1パルス到来時刻の幅内でのそれぞれのコスト値を算出し、コストベクトルとする。コスト算出器26は、コストベクトルに基づいて仮想音源Sijにおけるコスト関数の値Jijを算出する。接続端子27はコスト関数の値Jijのデータを含む信号を音源位置推定器5に送信するための端子である。
例えば、S/N比が低く、第1パルスの推定精度が高くない場合、基準となる第1パルスの到来時刻推定値を固定してコスト関数の値Jijを算出すると、第1パルスの到来時刻にずれが生じた場合には、求める包絡線CPEの値CPE(t1+τ1z0 (i,j))もずれるため、コスト関数の値Jijも異なり、音源位置推定精度も劣化する可能性がある。そこで、本実施の形態では、第1パルスの到来時刻に幅を持たせて、その幅で第1パルス到来時刻推定値を変化させながら、コスト関数の値Jijを算出し、音源位置推定の精度を高めようとするものである。
図18は包絡線CPEと第1パルス到来時刻幅との関係を表す図である。次に、第5の実施の形態における音源位置推定装置の動作について、図16、17のブロック図に沿って、図18、19を参照しながら説明する。第1パルス到来時刻幅算出器20は、包絡線CPEのデータに基づいて、包絡線CPEの中であらかじめ設定した閾値を越えた部分をもつ包絡線を第1パルスを含む包絡線であると判断し、越えた部分に基づいて、その閾値を越える開始時刻tstart 及び終了時刻tend を求める。求めた開始時刻tstart 及び終了時刻tend のデータを含む信号をコスト関数算出器21に出力する。
図19は包絡線CPEと到来時間差Table値との関係を表す図である。コスト関数算出器21のコストベクトル算出器24には、包絡線CPEのデータ並びに開始時刻tstart 及び終了時刻tend のデータを含む信号が、それぞれ入力端子22、23を介して入力される。また、到来時間差Table作成器3に記憶された、仮想音源Sijにおけるパルスの到来時間差Table値τ120 (i,j),τ130 (i,j),…,τ1N0 (i,j)のデータを含む信号が接続端子25を介して入力される。
コストベクトル算出器24は、まず、開始時刻tstart 及び終了時刻tend のデータに基づいて、次式(11)に示すΔt(s)を算出する。ここで、fsはサンプリング周波数であり、s=1,2,…,W、W=(tend −tstart )fs+1である。
Δt(s)=(s−1)/fs …(11)
Δt(s)=(s−1)/fs …(11)
そして、包絡線CPEのデータ、第1パルスの到来時刻推定値t1及びパルスの到来時間差Table値のデータに基づいて、図19に示すように、時刻tstart +Δt(s)+τ120 (i,j),tstart +Δt(s)+τ130 (i,j),…,tstart +Δt(s)+τ1N0 (i,j)における包絡線CPEの値を求める(z=2,3,…,N)。あるΔt(s)における包絡線CPEの値の総和に対し、その逆数の値をコストベクトルJVECijの第s要素とする。以上を式で表すと、次式(12)となる。これにより、仮想音源Sijについてのコストベクトルは、JVECij=[J1ij,J2ij,…,Jsij,…JWij]となる。
コスト算出器26は、次式(13)に示すように、コストベクトルJVECijに基づいて、各要素の値の中で最小のものを仮想音源Sijにおけるコスト関数の値Jijとして決定する。そして、各仮想音源Sijにおけるコスト関数の値のデータを含む信号を、接続端子27を介して音源位置推定器5に出力する。音源位置推定器5での音源位置の推定処理は、第1の実施の形態で説明した処理と同様であるので説明を省略する。
Jij=MIN[JVECij]
ここで、i=1,2,…,p、j=1,2,…,qであり、MIN[] は、[]内の要素の最小値をとることを表す。
Jij=MIN[JVECij]
ここで、i=1,2,…,p、j=1,2,…,qであり、MIN[] は、[]内の要素の最小値をとることを表す。
以上のように第5の実施の形態によれば、第1パルス到来時刻幅算出器20が開始時刻tstart及び終了時刻tendを求め、第1パルスの到来時刻推定値にtstart〜tendまでの幅を持たせ、コストベクトル算出器24がその幅の間で第1パルス到来時刻推定値を変化させて、コスト算出器26が仮想音源Sijのコスト関数の値Jijをデータとして算出し、各仮想音源Sijのコスト関数の値Jij中の最小値に基づいて音源位置推定器6が音源位置を推定するようにしたので、精度の高い第1パルスの到来時刻推定値に基づいて音源位置を推定することができる。そのため、第1パルスの到来時刻推定誤差に対してロバストな音源位置を推定することができる。
実施の形態6.
上述の第5の実施の形態では、コストベクトルJVECijの要素の中の最小値を仮想音源Sijにおけるコスト関数の値Jijと決定したが、コストベクトルJVECijの要素の総和の値、平均値等をコスト関数の値Jijとして決定するようにしてもよい。
上述の第5の実施の形態では、コストベクトルJVECijの要素の中の最小値を仮想音源Sijにおけるコスト関数の値Jijと決定したが、コストベクトルJVECijの要素の総和の値、平均値等をコスト関数の値Jijとして決定するようにしてもよい。
また、上述の第4及び第5の実施の形態では、レプリカ信号と音響センサによる受信信号との相関処理結果の自乗値の包絡線CPE及び到来時間差Table値に基づいてコスト関数の値Jijを算出するようにしたが、例えば、受波した音波自体の自乗値の包絡線を用いてコスト関数の値Jijを算出することもできる。本発明は、1つの音響センサが受波した音波による信号を用いて行ってもよいし、複数の音響センサで受波した音波を整相処理した信号を用いて行ってもよい。
さらに、第1〜第5の実施の形態で示した方法を、それぞれ並列に実行し、算出した音源位置の推定結果の中から尤もらしい結果に基づいて、最終的な音源位置を推定するようにすることもできる。
上述の実施の形態では、音源位置推定器5において、推定対象と仮想音源との一致度をコストとして算出し、最も一致度が高い仮想音源の位置を推定対象の位置とした。ただ、本発明はこれに限定されるものではない。例えば、最も一致度が高い仮想音源のコスト関数の値と、その近傍(例えば仮想音源がマトリクス状に配置されている場合、隣接する8つの仮想音源)のコスト関数との値に基づいて、仮想音源間の位置に推定対象の位置を求めることにより、さらに精度の高い位置推定を行うようにしてもよい。
実施の形態7.
上述の実施の形態では、音波が伝搬する媒体を海水としたが、本発明は、液体をはじめとする他の媒体を伝搬する音波を発する音源の位置推定についても適用することができる。
上述の実施の形態では、音波が伝搬する媒体を海水としたが、本発明は、液体をはじめとする他の媒体を伝搬する音波を発する音源の位置推定についても適用することができる。
1 入力端子
2 マルチパス波到来時間差測定器
3 到来時間差Table作成器
4 コスト関数算出器
5 音源位置推定器
6 出力端子
9、10 入力端子
11 Table値−観測値差算出器
12 自乗和算出器
13 接続端子
16 レプリカ相関器
17 包絡線算出器
18 第1パルス到来時刻推定器
19 コスト関数算出器
20 第1パルス到来時刻幅算出器
21 コスト関数算出器
22、23、25 入力端子
24 コストベクトル算出器
26 コスト算出器
27 接続端子
2 マルチパス波到来時間差測定器
3 到来時間差Table作成器
4 コスト関数算出器
5 音源位置推定器
6 出力端子
9、10 入力端子
11 Table値−観測値差算出器
12 自乗和算出器
13 接続端子
16 レプリカ相関器
17 包絡線算出器
18 第1パルス到来時刻推定器
19 コスト関数算出器
20 第1パルス到来時刻幅算出器
21 コスト関数算出器
22、23、25 入力端子
24 コストベクトル算出器
26 コスト算出器
27 接続端子
Claims (12)
- 音源の位置を推定する領域において複数の仮想音源を仮想的に配置したときの、前記領域の環境に応じた、仮想音源から受波地点までの複数の伝搬経路による音波の各到来時刻の理論値に基づく値を各仮想音源に対して算出する工程と、
前記受波地点において観測した音波に基づいて、推定対象の音源からの音波に関する値を算出する工程と、
前記推定対象の音源からの音波に関する値と前記理論値に基づく値から、前記推定対象の音源からの音波による到来時刻と一致又は最も近い到来時刻の仮想音源を決定し、該仮想音源の配置位置を前記推定対象の音源の位置と推定する工程と
を有することを特徴とする音源位置推定方法。 - 音源の位置を推定する領域において複数の仮想音源を仮想的に配置したときの、前記領域の環境に応じた、仮想音源から受波地点までの複数の伝搬経路による音波の到来時刻の理論値に基づいて、最初の音波が到来する時刻と他の伝搬経路による音波が到来する時刻との差の理論値を各仮想音源に対して算出する工程と、
前記受波地点において観測した音波に基づいて、推定対象の音源からの最初の音波が到来した時刻と他の経路により到来した時刻との差の観測値を算出する工程と、
前記推定対象の音源による前記差の観測値と前記各仮想音源による前記差の理論値とにより算出した一致度に基づいて、最も高い一致度が算出された前記仮想音源の配置位置を前記推定対象の音源の位置と推定する工程と
を有することを特徴とする音源位置推定方法。 - 前記一致度を算出するために、前記推定対象の音源による前記差の観測値と前記各仮想音源による前記差の理論値の組み合わせを決定する際に、
前記差の値が互いに最も近いものを1対1で対応させて組み合わせを決定することを特徴とする請求項1記載の音源位置推定方法。
前記差の観測値と前記差の理論値とが互いに最も近いものを1対1で対応させて組み合わせを決定することを特徴とする請求項2記載の音源位置推定方法。 - 前記一致度を算出するために、前記推定対象の音源による前記差の観測値と前記各仮想音源による前記差の理論値の組み合わせを決定する際に、
前記実際の音源による前記差の観測値に対して、前記仮想音源による前記差の理論値で最も近い値をそれぞれ対応させ、
前記差の理論値において組み合わせが決定しなかった値について、前記差の観測値の中で最も近い値を対応させて、組み合わせを決定することを特徴とする請求項2記載の音源位置推定方法。 - 音源の位置を推定する領域において複数の仮想音源を仮想的に配置したときの、前記領域の環境に応じた、各仮想音源から受波地点までの音波の伝搬経路及びそれぞれの前記伝搬経路による音波の到来時刻の理論値を求め、最初の音波が到来する時刻と他の経路による音波が到来する時刻との差の理論値を各仮想音源に対して算出する工程と、
観測した音波に基づいて、推定対象の音源から発せられる音波に基づく値を時系列に算出する工程と、
前記推定対象の音源から発せられる音波が最初に到来する時刻を判断し、該時刻と前記仮想音源における最初の音波が到来する時刻とを一致させたときの、前記差の理論値における前記推定対象の音源から発せられる音波に基づく値から、前記仮想音源と前記推定対象の音源との一致度を算出し、最も高い一致度が算出された前記仮想音源の位置を前記推定対象の音源の位置として推定する工程と
を有することを特徴とする音源位置推定方法。 - 前記推定対象の音源から発せられる音波が最初に到来する時刻としてとり得る時間幅を設定し、前記時間幅内の複数の時刻と前記仮想音源に対して最初の音波が到来する時刻とを一致させたときの前記一致度を前記各仮想音源について算出することを特徴とする請求項5記載の音源位置推定方法。
- 前記推定対象の音源から発せられる音波に基づく値は、推定対象の音源から発せられる音波とあらかじめ設定されたレプリカ用音波との相関関係を表す値とし、前記各仮想音源の前記差の理論値に対応する時刻における前記相関関係を表す値の総和又は平均に基づく値を、前記各仮想音源の前記一致度とすることを特徴とする請求項5又は6記載の音源位置推定方法。
- 複数の音響センサにより観測された、前記推定対象の音源から発せられる音波に基づいて前記推定対象の音源の位置を推定することを特徴とする請求項1〜7のいずれかに記載の音源位置推定方法。
- 音源の位置を推定する領域において複数の仮想音源を仮想的に配置したときの、前記領域の環境に応じた、仮想音源から受波地点までの複数の伝搬経路による音波の到来時刻の理論値に基づいて、最初の音波が到来する時刻と他の伝搬経路による音波が到来する時刻との差の理論値を各仮想音源に対して算出し、データとして記憶する到来時間差理論値算出手段と、
前記受波地点において観測した音波のデータに基づいて、推定対象の音源からの最初の音波が到来した時刻と他の経路により到来した時刻との差の観測値を算出する到来時間差観測値算出手段と、
前記推定対象の音源による前記差の観測値と前記各仮想音源による前記差の理論値とにより算出した一致度をコストとして算出するコスト関数算出手段と、
該コスト関数算出手段が算出した前記コストに基づいて最も高い一致度が算出された前記仮想音源の配置位置を前記推定対象の音源の位置と推定する音源位置推定手段と
を備えたことを特徴とする音源位置推定装置。 - 音源の位置を推定する領域において複数の仮想音源を仮想的に配置したときの、仮想音源から受波地点までの複数の伝搬経路による音波の到来時刻の理論値に基づいて、最初の音波が到来する時刻と他の伝搬経路による音波が到来する時刻との差の理論値を各仮想音源に対して算出し、データとして記憶する到来時間差理論値算出手段と、
観測した音波に基づいて、推定対象の音源から発せられる音波に基づく値のデータを時系列に算出するデータ算出手段と、
前記推定対象の音源から発せられる音波が最初に到来する時刻を推定する到来時刻推定手段と、
前記推定対象の音源から発せられる音波が最初に到来する時刻と各前記仮想音源に対して最初の音波が到来する時刻とに基づいて、前記差の理論値における前記推定対象の音源から発せられる音波に基づく値を判断する処理手段と、
該処理手段が判断した前記推定対象の音源から発せられる音波に基づく値に基づいて、前記推定対象の音源と各前記仮想音源との一致度をコストとして算出するコスト関数算出手段と、
該コスト関数算出手段が算出した前記コストに基づいて最も高い一致度が算出された前記仮想音源の位置を前記推定対象の音源の位置と推定する音源位置推定手段と
を備えたことを特徴とする音源位置推定装置。 - 音源位置推定手段は、最も高い一致度が算出された前記仮想音源の位置を前記推定対象の音源の位置と推定する代わりに、前記最も高い一致度が算出された仮想音源及び該仮想音源の近傍に位置する仮想音源について算出したコストに基づいて、前記推定対象の音源の位置を推定することを特徴とする請求項9又は10記載の音源位置推定装置。
- 観測した音波に基づく数値のデータを含む信号に変換する音響センサと、
請求項9、10又は11に記載の音源位置推定装置とを
備えたことを特徴とするソーナー。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005004243A JP2006194627A (ja) | 2005-01-11 | 2005-01-11 | 音源位置推定方法及び装置並びにソーナー |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005004243A JP2006194627A (ja) | 2005-01-11 | 2005-01-11 | 音源位置推定方法及び装置並びにソーナー |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2006194627A true JP2006194627A (ja) | 2006-07-27 |
Family
ID=36800840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005004243A Pending JP2006194627A (ja) | 2005-01-11 | 2005-01-11 | 音源位置推定方法及び装置並びにソーナー |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2006194627A (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010511172A (ja) * | 2006-11-30 | 2010-04-08 | カルディナーレ、チコッティ、ジュセッペ | 音響波および電磁波を使用して遠隔デバイスを標定するための方法 |
JP2011149864A (ja) * | 2010-01-22 | 2011-08-04 | Oki Electric Industry Co Ltd | 音源位置推定装置 |
JP2013531785A (ja) * | 2010-05-20 | 2013-08-08 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 音声信号を使用した距離推定 |
KR102161246B1 (ko) * | 2020-08-11 | 2020-09-29 | 국방과학연구소 | 해저 고정형 수신기를 갖는 양상태 소나 시스템 및 이를 이용한 해양 환경 기반 음원 배치 방법 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59225374A (ja) * | 1983-06-06 | 1984-12-18 | Hitachi Ltd | 音源位置標定方法とその装置 |
JPH04238284A (ja) * | 1991-01-22 | 1992-08-26 | Oki Electric Ind Co Ltd | 音源位置推定装置 |
JPH06324140A (ja) * | 1993-05-17 | 1994-11-25 | Oki Electric Ind Co Ltd | 音源位置推定装置 |
JPH1062534A (ja) * | 1996-08-26 | 1998-03-06 | Tech Res & Dev Inst Of Japan Def Agency | 目標検出装置 |
JP2004333252A (ja) * | 2003-05-06 | 2004-11-25 | Nippon Telegr & Teleph Corp <Ntt> | 位置推定装置および位置推定方法 |
-
2005
- 2005-01-11 JP JP2005004243A patent/JP2006194627A/ja active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59225374A (ja) * | 1983-06-06 | 1984-12-18 | Hitachi Ltd | 音源位置標定方法とその装置 |
JPH04238284A (ja) * | 1991-01-22 | 1992-08-26 | Oki Electric Ind Co Ltd | 音源位置推定装置 |
JPH06324140A (ja) * | 1993-05-17 | 1994-11-25 | Oki Electric Ind Co Ltd | 音源位置推定装置 |
JPH1062534A (ja) * | 1996-08-26 | 1998-03-06 | Tech Res & Dev Inst Of Japan Def Agency | 目標検出装置 |
JP2004333252A (ja) * | 2003-05-06 | 2004-11-25 | Nippon Telegr & Teleph Corp <Ntt> | 位置推定装置および位置推定方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010511172A (ja) * | 2006-11-30 | 2010-04-08 | カルディナーレ、チコッティ、ジュセッペ | 音響波および電磁波を使用して遠隔デバイスを標定するための方法 |
JP2011149864A (ja) * | 2010-01-22 | 2011-08-04 | Oki Electric Industry Co Ltd | 音源位置推定装置 |
JP2013531785A (ja) * | 2010-05-20 | 2013-08-08 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 音声信号を使用した距離推定 |
KR102161246B1 (ko) * | 2020-08-11 | 2020-09-29 | 국방과학연구소 | 해저 고정형 수신기를 갖는 양상태 소나 시스템 및 이를 이용한 해양 환경 기반 음원 배치 방법 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101106047B1 (ko) | 능동형 소나시스템의 표적거리 오차추정 방법 및 그 능동형 소나시스템 | |
KR20070019987A (ko) | 물체의 위치추정 | |
JP5383374B2 (ja) | 水中探知装置 | |
WO2007145761A2 (en) | Methods and systems for passive range and depth localization | |
CN101482614A (zh) | 声音传播速度建模方法、装置和*** | |
CN108089155A (zh) | 一种深海环境下单水听器声源被动定位方法 | |
JP6075388B2 (ja) | 環境計測システム及び環境計測方法 | |
KR101885991B1 (ko) | 실시간 해양환경 모니터링을 위한 음향 토모그래피 시스템 및 그 방법 | |
RU2634787C1 (ru) | Способ обнаружения локального объекта на фоне распределенной помехи | |
JP2006194627A (ja) | 音源位置推定方法及び装置並びにソーナー | |
JP2007040734A (ja) | 送信モードリコメンドシステム、ソーナーシステム、送信モードリコメンド方法、そのプログラム及び記録媒体 | |
KR101480834B1 (ko) | 다중 경로 음파 전달 모델 및 표적 식별을 이용한 표적 기동분석 방법 | |
JP5757303B2 (ja) | 水中音響測位システム | |
JP3573090B2 (ja) | 水中目標物位置検出装置及び方法 | |
JP2006292435A (ja) | マルチスタティック計測方法及び方式 | |
CN111337881B (zh) | 一种利用螺旋桨噪声的水下目标探测方法 | |
JP4356530B2 (ja) | パルス音の到来時間差推定方法及びその装置 | |
JP2013170817A (ja) | 音源位置推定装置 | |
CN106680824A (zh) | 一种基于垂直机动的水下航行器自主目标被动定位方法 | |
JP5708018B2 (ja) | アクティブソーナー装置 | |
JP2018146354A (ja) | ソーナー画像処理装置、ソーナー画像処理方法およびソーナー画像処理プログラム | |
CN109799480A (zh) | 一种水下通信节点基于多途信道约束的长基线定位方法 | |
Yayu et al. | Research on location of underwater sound source target in deep sea sound field based on bellhop model | |
RU2798390C1 (ru) | Способ пассивного определения координат шумящего объекта | |
WO2023106237A1 (ja) | 空間位置算出装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070131 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090527 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090707 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20091110 |