(57)【要約】【課題】フーリエ変換によるピークパワースペクトル検出において、測定対象とする時間範囲の分解能を高めつつ、周波数分解能を高める。【解決手段】時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する方法は、信号を所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング段階と、必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、基準サンプル数に達するように、デジタルデータ列にゼロデータを付加する付加段階と、ゼロデータが付加されたデジタルデータ列に対して高速フーリエ変換を行うＦＦＴ段階と、を有する。高速フーリエ変換の結果において最大値を示す周波数をピークパワー周波数として検出する。【選択図】図７
2012247304.tif 000002
【特許請求の範囲】
【請求項１】
時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する方法であって、
前記信号を前記所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング段階と、
必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、該基準サンプル数に達するように、前記デジタルデータ列にゼロデータを付加する付加段階と、
前記ゼロデータが付加された前記デジタルデータ列に対して高速フーリエ変換を行うＦＦＴ段階と、
を有し、
前記高速フーリエ変換の結果において最大値を示す周波数を前記ピークパワー周波数として検出する方法。
【請求項２】
前記基準サンプル数は、前記サンプリング段階で得られる前記デジタルデータ列のサンプル数の２倍以上である、請求項１に記載の方法。
【請求項３】
前記信号は、所定の周波数帯域を有し、前記所定の周波数帯域が低周波側に変換されるように、前記サンプリング段階の前に、前記信号に対して周波数変換を行う段階を有する、請求項１または２に記載の方法。
【請求項４】
時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する装置であって、
前記信号を前記所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング手段と、
必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、該基準サンプル数に達するように、前記デジタルデータ列にゼロデータを付加するゼロ付加手段と、
前記ゼロデータが付加された前記デジタルデータ列に対して高速フーリエ変換を行うＦＦＴ手段と、
を有する装置。
【請求項５】
前記基準サンプル数は、前記サンプリング手段で得られる前記デジタルデータ列のサンプル数の２倍以上である、請求項４に記載の装置。
【請求項６】
前記サンプリング手段は、前記信号に対してアナログ／デジタル変換を行うＡ／Ｄ変換器と、前記所望の時間幅に応じ、前記Ａ／Ｄ変換器から出力されるデジタルデータ列を切り出す信号切り出しゲートと、を有する、請求項４または５に記載の装置。
【請求項７】
前記信号は所定の周波数帯域を有しており、
さらに、前記サンプリング手段の前に、前記所定の周波数帯域が低周波側に変換されるように前記信号に対して周波数変換を行う周波数変換手段を有する、請求項４乃至６のいずれか１項に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、時間に応じて変化する信号の中から所望の時間幅内の部分を切り出し、その切り出された部分における信号のピークパワースペクトルを検出する方法と装置とに関し、特に、検出の高速化を図りつつ、検出されるスペクトルの周波数精度を向上させた方法及び装置に関する。
【背景技術】
【０００２】
所定の周波数帯域幅を有して時々刻々と変化する信号の中から所望の時間幅の部分のみを切り出し、その切り出された部分内での信号のピークパワースペクトルを検出することは、後述するように、例えばドップラー計測などで広く行われている。ここでいう所望の時間幅は、信号全体の継続時間に比べた場合には十分に短く、例えば、その開始時刻が与えられているかあるいは開始時刻を自動検出できるものである。また、検出されるピークパワースペクトルは、切り出された部分内での信号の周波数スペクトルの中で一番強いパワースペクトルを有するものである。パワースペクトルの検出には、一般的にフーリエ変換が用いられ、計算量の観点から、中でも、ＦＦＴ（高速フーリエ変換：Fast Fourier Transform）が使用される。
【０００３】
ところで、ある信号を媒質中に発射（送信）し、その媒質中から反射した信号（エコー信号）を受信したとき、送信した時刻をスタート時刻とすると、スタート時刻から反射によるエコー信号を受信するまでの時間は、媒質中でのその信号の伝搬距離すなわち対象物までの距離に依存する。したがって、信号の伝搬速度（音波であれば音速）を考慮して、スタート時刻から所望の時間経過した時点で受信されるエコー信号は、その時間経過に応じた距離の位置にある対象物からのエコー信号となる。スタート時刻からの時間（すなわち距離）に応じてエコー信号の強度を映像として表示するようにした応用機器が、例えば、レーダーや魚群探知機である。
【０００４】
信号を反射する対象物が相対的に移動しているとすると、反射エコーでは、ドップラー効果により、送信した周波数からの周波数変位が観測される。この周波数差、すなわち媒質に対して送信したときの周波数と反射エコーで観測される周波数との差をドップラー周波数という。ドップラー計測では、一定の周波数の信号を媒質中に発射し、媒質中の対象物によって反射された信号（反射エコー）を受信し、その周波数解析を行って、ドップラー周波数を決定するものである。ドップラー計測を行うものとして、例えば、ドップラーレーダー、ドップラーソナー、ＡＤＣＰ（音波ドップラー流速プロファイラー；acoustic Doppler current profiler）、潮流計、音波レーダー、超音波ドップラー速度計などがあり、これらは、ドップラー周波数に基づいて、対象物との相対速度を決定する。
【０００５】
送信回路によって一定の送信周波数ｆtxの信号を生成して、この信号を送受波器から媒質に発射する。媒質中を対象物（反射体）が移動しているとして、ドップラー周波数ｆdopは、対象物の移動速度をＶ（ｍ／ｓ）、媒質中での信号の伝搬速度をＣ（ｍ／ｓ）とすると、媒質中の信号の伝搬速度が対象物の速度よりも十分に速い場合（Ｃ≫Ｖ）には、次の近似式で表される。
【０００６】
【数１】
2012247304.tif 000003
【０００７】
ここでθは、対象物から送受波器を見たときの方向と対象物の移動方向とがなす角度であり、同じ方向であれば（対象物が送受波器にむかってまっすぐ進んでいるときには）、θ＝０となる。したがって、ドップラー周波数ｆdopは、対象物が近付いているときには正（＋）となり、遠ざかるときには負（−）となる。
【０００８】
エコー信号の周波数解析を行ってドップラー周波数ｆdopを求めるが、この周波数解析には、一般的に、フーリエ変換の手法が使われている。フーリエ変換として、パワースペクトル解析を行う場合にはＤＦＴ（離散フーリエ変換；Discrete Fourier Transform）が用いられ、あるいは信号の標本数（サンプル数）が２のべき乗個であるときには、高速演算処理が可能なＦＦＴ（高速フーリエ変換：Fast Fourier Transform）が通常用いられている。ＤＦＴもＦＦＴも複素乗算演算を反復するものであるが、サンプル数をＮとすると、ＤＦＴでの複素乗算数回数はＮ×Ｎ回（すなわちＮ2回）であり、Ｎ＝２nとすると、ＦＦＴでの複素乗算回数は、２n×ｌｏｇ2 ２n回（すなわちＮ×ｌｏｇ2 Ｎ回）であり、Ｎが大きくなるほど、ＦＦＴはＤＦＴに比べて高速演算ができる。
【０００９】
ところで、ドップラー計測では、対象物の位置と速度の両方に関する情報を取得したい場合が多い。例えば、海洋において深さごとの潮流の方向と速度とを求めて流速プロファイルを得るためには、海中に対して斜め方向に音波パルスを放射して反射エコーを取得し、周波数解析を行う。この際、音波パルスの放射から所定の時間の経過後にある時間幅で反射エコーを切り出し、切り出された部分に対して周波数解析を行うことになる。音波パルスの発射からの所定の時間によって、どの深度の流速を計測しようとするかが定まることになり、反射エコーを切り出す時の時間幅によって、どの程度の深度範囲（厚さ）にわたる平均として流速が算出されるかが定まることになる。例えば、切り出す時の時間幅が大きければ、海中における厚みの大きい層での平均的な流速が求められることになる。これに対し、時間幅が小さければ、薄い層での流速を求めることが可能になって、結果として、細かな流速プロファイルが得られることになる。
【００１０】
しかしながら、ドップラー周波数を求めるための周波数解析としてのフーリエ変換の問題点として、周波数分解能と時間（距離）分解能とが相反する関係にある、ということが挙げられる。ここでいう周波数分解能は、ドップラー周波数を求める際の分解能であり、これは速度の分解能に比例する。時間分解能とは、反射エコーを切り出す時にどこまで時間幅を小さくできるかということであり、これは、ドップラー計測における距離あるいは位置の分解能に比例する。時間分解能が粗ければ、上述した流速プロファイルの場合であれば、深度方向の厚みが大きな層でのおおざっぱな平均的な流速しか得られないことになる。
【００１１】
離散フーリエ変換（ＤＦＴ）は、次の式で表される。
【００１２】
【数２】
2012247304.tif 000004
【００１３】
ＤＦＴまたはＦＦＴでのサンプリング周波数をｆft、フーリエ変換の基本波の周期をＴ0とすると、周期の逆数は、フーリエ変換での基本周波数ｆ0と呼ばれ、
ｆ0＝１／Ｔ0                        ：式１−４
と表される。この基本周波数ｆ0は、ＤＦＴやＦＦＴでの最小分解周波数すなわち周波数の分解能となる。サンプルされた信号点の総数すなわちサンプル数をＮとすると、
Ｔ0＝Ｎ／ｆft                        ：式１−５
が当然に成立し、周波数分解能である最小分解周波数ｆ0は、
ｆ0＝１／Ｔ0＝ｆft／Ｎ                ：式１−６
と表されることになる。逆に、所望の周波数分解能ｆ0が与えられたとして、必要なサンプル数Ｎは、
Ｎ＝ｆft／ｆ0                        ：式１−７
と表されることになる。フーリエ変換では、基本波の周期Ｔ0と等しい時間幅の信号から等間隔でＮ個のサンプルを取り出すから、Ｔ0が信号の切り出し時間（切り出しの時間幅）ということになる。
【００１４】
ＦＦＴで得られるパワースペクトルの例が図１に示されている。ＦＦＴのパワースペクトルは、サンプル数Ｎに対し、０から（Ｎ／２）−１までの一定の周波数間隔で周波数ごとの値１４，１５が与えられるスペクトルとして表される。ここで周波数間隔は上記のｆ0と一致し、ｆ0＝１／Ｔ0である。
【００１５】
式１−４に示すように、周波数分解能ｆ0は、信号の切り出し時間幅Ｔ0の逆数であるから、周波数分解能を上げる（ｆ0を低くする）ためには、信号の切り出し時間幅Ｔ0を長くしなければならない。その結果、周波数分解能を２倍細かくすると、時間分解能すなわち位置（距離）の分解能は２倍粗くなって劣化する。
【００１６】
対象物との相対的な速度によって、ドップラー周波数ｆdopは変化する。このドップラー周波数ｆfopが変化する範囲をΔｆpとおくことにすると、近づく方向と遠ざかる方向での移動が同等に起こるものと仮定すれば（この過程は、潮流計などの用途において妥当な仮定である）、ドップラー信号の周波数範囲ΔｆPの中心周波数ｆc、つまりドップラー周波数ｆdopがゼロのときの受信信号の周波数は、送信周波数ｆtxと一致する（ｆc＝ｆtx）。
【００１７】
以下の条件に従い、超音波のパルス状の信号を海中に送信し、移動物からの反射エコーのドップラー周波数を測定する場合を考える。
【００１８】
海水中の音波伝搬速度Ｃ：Ｃ＝１５００（ｍ／ｓ）
送信周波数ｆtx：ｆtx＝１２０ｋＨｚ
検出最大速度（水平方向）Ｖ：Ｖ＝１５（ｍ／ｓ）
また、音波は水平方向から斜め方向（θ＝６０°）に送受信するものとする。式１−１により、
【００１９】
【数３】
2012247304.tif 000005
【００２０】
となり、検出最大速度が１５ｍ／ｓであるとすれば、ドップラー信号として観測される周波数範囲Δｆpは、１２０±１．２（ｋＨｚ）となる。
【００２１】
測定速度範囲すなわち検出最大速度は、通常、予め決められているので、ドップラー信号についても、その周波数範囲が予め定められて測定する。
【００２２】
ドップラー計測による速度精度を０．１ｍ／ｓとすると、それに必要な周波数分解能ｆ0は１２Ｈｚとなる。ドップラー計測において検出対象とする位置の範囲（位置分解能）を７．５ｍとすると、信号の切り出し時間幅は、海水中で音波が１ｍを往復する時間である２／１５００（ｍ／ｓ）を乗ずると、７．５×２／１５００＝０．０１となって、１０ｍｓとなる。
【００２３】
しかしながら、周波数分解能ｆ0である１２Ｈｚを式１−４に代入すると、Ｔ0＝１／ｆ0＝１／１２≒８３．３ｍｓとなり、８３．３ｍｓは、海水中で音波が１ｍを往復する時間２／１５００（ｍ／ｓ）で割ることにより距離に換算すると、８３．３（ｍｓ）×１５００（ｍ／ｓ）／２≒６２．５（ｍ）となる。
【００２４】
すなわち、速度精度を０．１ｍ／ｓとするために周波数分解能ｆ0を１２Ｈｚとすると、信号の切り出し時間幅として８３ｍｓが必要であり、これは距離に換算すれば６３ｍに相当し、位置（距離）に関する分解能を７．５ｍとする要求を満足することできない。逆に、位置に関する分解能が７．５ｍである信号でＦＦＴを行うとすると、式１−４より周波数分解能ｆ0が１００Ｈｚとなり、速度精度から定められた周波数分解能１２Ｈｚの要求を満足できない。
【００２５】
このようにフーリエ変換には、周波数分解能と、測定対象に関する時間分解能とが両立しない、という基本的な問題があり、これは、ドップラー計測の場合には、周波数分解能あるいは速度分解能と、距離あるいは位置の分解能とが両立しないことを意味する。
【００２６】
図２は、一般的なドップラー計測機器の構成の一例を示している。ここでは、信号が音響信号であり、信号が伝搬する媒質が水（あるいは海水）であるとする。音響信号の送受波器１は、送受信切替回路２を介し、送信回路３の出力と受信増幅器４の入力に接続している。送信回路３は、送信周波数ｆtxの信号を生成するものである。受信増幅器４の出力には、受信信号を中間周波数信号に周波数変換する変調器５が設けられており、変調器５には、局部発振回路６から局部発振周波数ｆlocの信号が供給されている。変調器５の出力すなわち中間周波数信号は、アナログフィルタ７を介してアナログ／デジタル変換器（Ａ／Ｄ）１０に供給されてデジタル信号に変換され、このデジタル信号は、所定のサンプル数（すなわち時間幅）で信号を切り出す信号切り出しゲート１１に供給される。信号切り出しゲート１１で切り出された信号は、次にＦＦＴ処理部１３に送られてＦＦＴによる周波数解析の対象となる。ＦＦＴ処理部１３から、計算されたパワースペクトルが出力される。
【００２７】
このドップラー計測機器では、水中に音波パルスを発射する時には、送信周波数ｆtxの信号を送信回路３から送受信切替え回路２を介して送受波器１に供給する。その後、直ちに送受信切替え回路２を受信側とし、媒質（水中）からの反射エコーを送受波器１で受信して電気信号である受信信号に変換し、送受信切替え回路２を介して受信信号を受信増幅器４に供給する。その結果、受信信号は受信増幅器４において信号増幅され、変調器５に対して供給されることになる。
【００２８】
図３は、受信信号や中間周波数信号などの周波数間の関係を示している。受信増幅器４から出力される信号自体は、図３において符号５２で示されるような周波数帯域を有するものとする。受信信号の帯域は、符号５２で示す周波数帯域の一部であり、検出対象となるドップラー信号は、符号５１で示される周波数範囲Δｆpの中にある。周波数範囲Δｆpの中心、すなわち、ドップラー周波数がゼロの時の周波数ｆc（４７）は、送信周波数ｆtxと同じである（ｆc＝ｆtx）。局部発振回路６は、周波数ｆcよりも中間周波数ｆmidだけ低い局部発振周波数ｆloc（４５）の局部発振信号を生成し（すなわちｆc＞ｆloc）、変調器５に供給する。その結果、変調器５では、局部発振信号を受信信号に乗算することにより、受信信号を中間周波数ｆmidの中間周波数信号に変換する。
【００２９】
受信信号の周波数帯域のうち、トップラー計測のために着目する周波数範囲Δｆpの上限を上限周波数ｆmax（４９）とし、下限を下限周波数ｆmin（４８）とすると、フーリエ変換すべき周波数範囲Δｆp（５１）は、Δｆp＝ｆmax−ｆminで表される。
【００３０】
変調器５における中間周波数ｆmidの信号の生成は、下記の式で表されるものであり、２つの変換周波数から、ローパスフィルタかバンドパスフィルタを使用して、低い方の周波数を得る。
【００３１】
【数４】
2012247304.tif 000006
【００３２】
ドップラー信号の周波数範囲Δｆpは、中間周波数ｆmidを中心に、
【００３３】
【数５】
2012247304.tif 000007
【００３４】
の範囲に変換されることとなり、検出対象の周波数範囲Δｆpは、周波数変換しても保たれる。上記の式は、ｆloc＜ｆcの時の処理を示しているが、ｆloc＞ｆcの時も、複号（符号（±））が逆になるだけで、同様に処理がなされる。
【００３５】
図３において、右側部分が受信信号と局部発振信号との関係を表しており、ここでｆmid（４６）は、ｆmid＝ｆc−ｆlocの関係を満たしている。一方、図３の左側部分では、中間周波数信号に変換後のスペクトルを示しており、ここでは中間周波数信号の中心周波数はｆmid（５０）となっている。また、上述の周波数帯域５２も周波数帯域５４に変換され、検出対象の周波数範囲Δｆp（５１）は、中心周波数がｆmidである周波数範囲Δｆp（５３）に変換されている。
【００３６】
図４は、潮流における流速プロファイル測定を行うために、観測船からθ＝６０°の方向（水平方向から６０°下を向いた方向）に対して音波パルスを発射して反射エコーを受信したときの、受波波形の一例を示している。図において横軸は時刻を表している。ここでは反射エコー１６は、中間周波数信号に変換されたのちの信号として示されている。時刻０の直後の振幅の大きな信号１７は、送信した音波パルス自体を送受波器が受信したことによるものである。また、振幅の大きい波形１８は、海底によって反射されたエコー（海底エコー）を示している。海底エコーにも、観測船の船速に応じ、±Δｆp／２の周波数範囲内でのドップラーシフトが含まれる。
【００３７】
ここで中間周波数ｆmid（すなわち周波数変換された信号の中心周波数）は５ｋＨｚであるとし、アナログ／デジタル変換でのサンプリング周波数及びＦＦＴでのサンプリング周波数ｆftは２０ｋＨｚであるとする。このとき、フーリエ変換でのサンプル数Ｎは、式１−７より、Ｎ＝ｆft／ｆ0＝２００００／１２≒１６６７となるが、ＦＦＴの使用を前提とすると、サンプル数は２のべき乗としなければならないので、サンプル数Ｎ＝２11＝２０４８とする。サンプル数Ｎが２０４８であるとき、周波数分解能は、式１−６よりｆ0＝２００００／２０４８＝９．７６６Ｈｚとなり、上述した速度精度の０．１ｍ／ｓに対応する周波数分解能１２Ｈｚよりも細かく、要求事項を満足している。
【００３８】
図４において、符号２０は、音波パルスの送信と同時に取り込みを開始するとして、ＦＦＴに必要な２０４８個の信号サンプルを取り込むのに要する時間Ｔ0を表している。この２０４８個の信号サンプルに対してＦＦＴを行った結果を図５に示す。
【００３９】
図５に示すスペクトルは、ＦＦＴステップで５１２番目に位置するピーク３８と、５７４番目あたりに位置するピーク３９とを顕著なピークとして含んでおり、強度としては、ピーク３８の方が大きい。ＦＦＴステップ数に周波数分解能ｆ0を乗じたものが実際の周波数であり、ここでは、ｆ0＝２００００／２０４８＝９．７６６Ｈｚであるから、ピーク３８の周波数は、９．７６６Ｈｚ×５１２＝５０００Ｈｚとなる。上述のように、中間周波数信号における中心周波数ｆmidは５ｋＨｚであって、これは送信周波数の信号を周波数変換したものと一致するから、ピーク３８は、送信波を直接検出したピークに対応する。また、同様の考察を行うと、ピーク３９は、海底エコーに相当するものであり、船速に応じて送信周波数がドップラーシフトした周波数となっている。ここに示した例では、サンプリング周波数ｆft＝２０ｋＨｚであってサンプル数Ｎ＝２０４８としているから、式１−５からＴ0＝１０２．４ｍｓとなる。これは、距離に換算すれば、送受波器の位置から海底までをまとめてＦＦＴの対象としていることになる。その結果、最もパワーが大きな送信信号自体がピークパワースペクトルとして検出されたことになる。つまり、所望する時間範囲以外に強い信号があるＦＦＴ入力信号では、その所望する時間範囲内でのピークパワーとなるスペクトルあるいは周波数を求めることはできない。
【００４０】
図４において、音波パルスの送信時刻を始点として時間ｔ1（２２）だけシフトさせ、ピーク１８だけを含む時間範囲Ｔ1（２１）の信号だけを切り出してＦＦＴを行うことも考えられる。ちょうど１２８個のサンプルが切り出されてＦＦＴを行ったとすると、周波数分解能ｆ0は、ｆ0＝２００００／１２８≒１５６Ｈｚであり、これは、ＦＦＴステップが１つずれれば１５６Ｈｚ異なるという、粗い精度である。この周波数分解能は、要求されている分解能１２Ｈｚと比べて極めて粗く、要求に応えることができない。
【００４１】
なお、説明のため、図６ではスペクトルが連続する実線で描かれているが、実際にはＦＦＴによる解析結果であるので、図１に示したものと同様に、ＦＦＴステップごとの離散値として表されるものである。
【００４２】
従来のドップラー計測では、所望の時間間隔内でピークパワーとなる周波数を求める場合、まず、距離（時間）分解能を優先させて、受信信号から必要な範囲Ｔ1を切り出し、周波数分解能が粗いままでＦＦＴによる周波数解析を実行する。解析結果において最も高い値となったサンプリング点、次に大きな値となったサンプリング点、また次に大きな値となったサンプリング点などから、スペクトルにおけるピークの形状を仮定したカーブフィッティングなどの手法によって補間を行うことにより、真のピークパワー周波数を推定する方法などを用いていた。補間によるので、真のピークパワー周波数は、ＦＦＴにおける周波数分解能よりも細かい分解能で算出できる。
【００４３】
しかしながら、このような従来のやり方では、補間によってピークパワー周波数を求めるための回路や推定アルゴリズムが複雑になり、その結果、ＦＦＴ以外の処理にも演算パワーを費やす必要が生じ、演算速度の低下やコストアップの要因となっていた。また、ＦＦＴの後で周波数推定を行うため、測定精度が悪化しやすいという問題もある。
【００４４】
特許文献１には、ＦＦＴにおける演算ポイント数を減らしつつ、周波数分解能を高めるようにした信号処理方法が示されている。
【先行技術文献】
【特許文献】
【００４５】
【特許文献１】ＷＯ２００６／０４３５１１
【発明の概要】
【発明が解決しようとする課題】
【００４６】
上述したように、ドップラー計測などにフーリエ変換を適用する場合、周波数分解能を向上させることと測定対象とする範囲の時間（あるいは距離）についての分解能を向上させることとが両立しない、というフーリエ変換の基本的な問題に起因する問題が生じる。この問題を回避するために、所望の短い時間幅で信号を切り出した上でＦＦＴを行い、そのＦＦＴ結果に対して各種の補間方法や推定方法を適用して、ＦＦＴによって離散化した周波数間隔よりも細かくピークパワーの周波数を求めることも行われているが、これは、構成の肥大化やアルゴリズムの複雑化をもたらし、そのために演算速度の低下やコストアップがもたらされる。さらに、これらの補間方法や推定方法を用いた場合には、必ずしも所望の周波数精度が得られない、という問題点もある。
【００４７】
本発明の目的は、フーリエ変換における上述した基本的な問題を解決し、測定対象とする時間範囲の分解能を高めつつ、周波数分解能を高めることができる検出方法及び装置を提供することにある。
【課題を解決するための手段】
【００４８】
本発明の検出方法は、時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する方法であって、信号を所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング段階と、必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、この基準サンプル数に達するように、デジタルデータ列にゼロデータを付加する付加段階と、ゼロデータが付加されたデジタルデータ列に対して高速フーリエ変換を行うＦＦＴ段階と、を有し、高速フーリエ変換の結果において最大値を示す周波数をピークパワー周波数として検出する。
【００４９】
本発明の検出装置は、時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する装置であって、信号を所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング手段と、必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、この基準サンプル数に達するように、デジタルデータ列にゼロデータを付加するゼロ付加手段と、ゼロデータが付加されたデジタルデータ列に対して高速フーリエ変換を行うＦＦＴ手段と、を有する。
【発明の効果】
【００５０】
本発明では、粗い周波数分解能でのＦＦＴの結果から補間や推定によってピークパワー周波数を求めるのではなく、所望の周波数分解能が得られるようにゼロデータを付加してからＦＦＴを行うので、ＦＦＴ後の補間や推定などを行うことなく、所望の周波数分解能でＦＦＴからピークパワー周波数を直接求めることができるようになる。本発明は、複雑な回路やアルゴリズムを必要としないので、シンプルな構成で、低コストで、精度良く、かつ高速に、ピークパワー周波数を求めることができる。
【図面の簡単な説明】
【００５１】
【図１】ＦＦＴによる周波数スペクトル解析結果を示すグラフである。
【図２】ドップラー計測機器の構成の一例を示すブロック図である。
【図３】受信信号の中間周波数への変換などを説明するスペクトル図である。
【図４】斜め方向に海中に音響信号を送波したときのその反射エコーの受信波形の一例を示す波形図である。
【図５】送信信号を含んだ反射エコーの周波数スペクトル解析結果の一例を示す図である。
【図６】入力サンプル数が少ない時のＦＦＴによる周波数スペクトル解析結果を示す図である。
【図７】本発明の実施の一形態のピークパワースペクトル検出装置の構成を示すブロック図である。
【図８】ＦＦＴへの入力信号を説明する波形図である。
【発明を実施するための形態】
【００５２】
図７は、本発明の実施の一形態のピークパワースペクトル検出方法の実施に適したドップラー計測機器の構成を示している。このドップラー計測機器は、本発明に基づくピークパワースペクトル検出装置の一例である。
【００５３】
図７に示したドップラー計測機器は、図２に示したドップラー計測機器において、信号切り出しゲート１１とＦＦＴ処理部１３との間に、ゼロ付加部１２を挿入したものである。信号切り出しゲート１１からは、測定対象の時間範囲に対応して信号をサンプリングしデジタル化して得られるデジタルデータ列が出力されるが、ゼロ付加部１２は、必要とされる周波数分解能を満たすためにＦＦＴにおいて必要なサンプル数を基準サンプル数として、この基準サンプル数に達するように、デジタルデータ列にゼロデータを付加するものである。ゼロ付加部１２での処理の詳細については後述する。
【００５４】
図７に示すドップラー計測機器では、図２に示したものと同様に、変調器５及び局部発振回路６を使用して、送受波器１で受信した信号における上述の周波数範囲Δｆpを低周波数側に周波数変換している。これは、高い周波数成分を検出するためにＦＦＴを行うためには、サンプリング周波数を高くする必要があり、サンプリング周波数が高い分だけサンプル数が膨大となり過ぎるからである。ドップラー計測の場合、例えば、１２０ｋＨｚ±１．２（ｋＨｚ）の周波数範囲内でドップラー周波数を検出すればよく、検出対象とすべき周波数範囲Δｆpが決まっているから、ＦＦＴにおいてサンプリング周波数を低下させ、かつ、必要な分解能の範囲内でサンプル数を減らすために、本実施形態では、受信信号を元の周波数よりも低い周波数の中間周波数信号に周波数変換している。前述の式２−１〜２−３などを用いて示したように、周波数変換後の中間周波数信号でも、周波数範囲Δｆpは保たれている。
【００５５】
変調器５の出力から低周波側の中間周波数成分だけを取り出すことと、後段のアナログ／デジタル変換においてエリアシングが生じないようにすることのために、アナログフィルタ７が設けられている。アナログ／デジタル変換器（Ａ／Ｄ）１０は、ＦＦＴの実行が可能なサンプリング周波数でアナログフィルタ７の出力をサンプリングしてアナログ／デジタル変換する。これにより、Ａ／Ｄ変換器１０から、中間周波数信号における瞬時値を表すデジタル値からなるデジタルデータ列が出力されることになる。信号切り出しゲート１１は、所望する時間範囲のデジタルデータ列を切り出す。この切り出されたデータのサンプル数を切り出しサンプル数とする。ゼロ付加部１２は、要求されている周波数分解能に対応するＦＦＴのサンプル数を基準サンプル数として求め、基準サンプル数に達するように、基準サンプル数と切り出しサンプル数との差に相当する数のゼロデータ（値がゼロであるサンプル）を切り出されたデジタルデータ列に付加する。そして、ゼロデータが付加されてサンプル数が基準サンプル数とされたデジタルデータ列に対し、ＦＦＴ処理部１３がＦＦＴを実行し、ピークパワースペクトルを算出する。算出されたピークパワースペクトルは、信号切り出しゲート１１において切り出された範囲に対するものであって、信号におけるその時間範囲でのピークパワー周波数を示している。また、ゼロデータが付加されたことにより、信号においてＦＦＴの対象となる時間幅が実質的に拡大されており、この時間幅の逆数が周波数分解能に該当することから、周波数分解能が向上している。このように、本実施形態によれば、受信信号の中での所望の時間範囲（あるいは距離範囲）のピークパワースペクトルを、ＦＦＴの演算だけで、所望の周波数精度を満たして検出することができる。
【００５６】
以下、数値例を挙げて、この実施形態における周波数分解能について説明する。
【００５７】
背景技術の欄で述べたものと同じ条件で、超音波のパルス状の信号を海中に送信し、移動物からの反射エコーのドップラー周波数を測定する場合を考える。すなわち、海水中の音波伝搬速度ＣをＣ＝１５００（ｍ／ｓ）とし、送信周波数ｆtxをｆtx＝１２０ｋＨｚとし、検出最大速度（水平方向）ＶをＶ＝１５（ｍ／ｓ）とし、音波を水平方向から斜め方向（θ＝６０°）に送受信するものとする。
【００５８】
このとき、ドップラー信号は１２０±１．２（ｋＨｚ）の範囲であり、速度精度として０．１ｍ／ｓを要求すれば周波数分解能ｆ0は１２Ｈｚとなる。また、ドップラー計測において検出対象とする位置の範囲（位置分解能）を７．５ｍとし、周波数変換後の信号の中心周波数すなわち中間周波数ｆmid（５０）を５ｋＨｚとし、アナログ／デジタル変換でのサンプリング周波数及びＦＦＴのサンプリング周波数ｆftを２０ｋＨｚとする。
【００５９】
周波数分解能ｆ0から必要なＴ0を求めると、Ｔ0＝１／ｆ0＝１／１２≒８３．３ｍｓとなる。
【００６０】
海中に音波パルスを発射して受信信号を観測し、その受信信号を中間周波数信号に変換した後の波形を示す図４において、上述したように、符号１６は反射エコーを表し、符号１７は送信信号自体を表し、符号１８は海底エコーを表している。海底エコー１８の周波数を解析し、式１−１に基づき、ドップラー周波数から船の速度を求める場合を考える。ここでは、送信信号などの不要な信号が含まれないようにするために、海底エコー１８の部分だけを切り出してＦＦＴを行い、海底からの反射波に関してドップラー周波数を求める。図８は、図４に示す信号波形から海底エコーだけを切り出した信号１９を示している。ここでは、海底エコー１９の部分として、音波パルスの送信から時間ｔ2（２６）の経過後に時間幅Ｔ2（２４）を切り出すものとする。ここで切り出された信号は、サンプル数としては１００個である。
【００６１】
時間幅Ｔ2の１００個のサンプルのみを用いてフーリエ変換を行った場合の周波数分解能ｆ0は、ｆ0＝２００００Ｈｚ／１００＝２００Ｈｚとなり、要求されている周波数分解能である１２Ｈｚを全く満たさない。
【００６２】
ここで、要求されている周波数分解能である１２Ｈｚをフーリエ変換が満たすために必要なサンプル数Ｎを求めると、式１−７から、Ｎ＝ｆft／ｆ0＝２００００／１２≒１６６７となる。１６６７は２のべき乗数ではないので、１６６７を超えて最も１６６７に近い２のべき乗数を考えると、２11＝２０４８が見つかる。そこで、ＦＦＴにおけるサンプル数Ｎとして２０４８を採用すると、式１−６より、周波数分解能ｆ0は、ｆ0＝２００００Ｈｚ／２０４８＝９．７６６Ｈｚとなって、要求事項を満足している。
【００６３】
そこでこの実施形態では、ＦＦＴをサンプル数Ｎ＝２０４８で実行するものとし、海底エコーの信号１９の分のサンプル数は１００なので、不足した１９４８個のサンプルについては、ゼロ値からなるサンプルを挿入する。図８においてＴ0（２３）は、２０４８個のサンプルに対応する時間であり、ＦＦＴ解析の対象となる時間幅である。
【００６４】
ここでＴ0＝Ｔ2＋Ｔ3の関係が成り立っており、測定したい範囲の信号１９を切り出してＦＦＴ解析に必要なサンプル数Ｎを満たすように付加されるゼロ値信号の継続時間が、Ｔ3に対応することになる。時間幅Ｔ3の挿入位置は、時間幅Ｔ2すなわち切り出された信号１９の前でも後でもよい。あるいは、信号１９を挟むように前後に分割してゼロ値信号を挿入する時間幅を配置してＦＦＴ解析の時間幅Ｔ0（２７）を設定してもよい。
【００６５】
ここで示す例では、音波パルスの送信から信号の切り出しを開始するまでの経過時間ｔ2（図８での符号２６）と時間幅Ｔ2（図８での符号２４）を予め定めているが、海底エコーのような強度の大きな反射エコーを検出対象とする場合であれば、切り出しまでの経過時間ｔ2も切り出しの時間幅Ｔ2も自動検出して決定することも可能である。深度ごとの潮流の流速を求める場合などには、自動検出を行うことができないのでｔ2及びＴ2について予め設定しておくことになる。
【００６６】
ここでゼロ付加部１２によるゼロ値データの付加について説明する。
【００６７】
一般にＦＦＴではその対象とするサンプル数は２のべき乗である必要がある。その一方で、ＦＦＴによる解析を行いたい対象データにおけるサンプル数が２のべき乗となっていないことはよくあることである。そこで、サンプル数が２のべき乗となるように、ゼロデータのサンプルを付加することがゼロパディングなどと称して行われている。その場合、とにかくＦＦＴを行えるだけのサンプル数となればよいのであるから、対象データのサンプル数をＮsとして、２m-1＜Ｎs≦２mを満たす自然数ｍが存在する時に、２m−Ｎs個のゼロ値データを付加することになる。それに対して本実施形態では、単にＦＦＴの実行が可能になるようにゼロ値データを付加するのではなく、ＦＦＴにおけるサンプリング周波数ｆftを考慮して所望の周波数分解能が得られるように、ゼロ値データを付加する。したがって、付加されるゼロ値データの数も、単にＦＦＴを実行可能になるようにする場合に比べ、はるかに多くなることがある。上述した例では、海底エコーを切り出してサンプル数が１００になったとして、単にＦＦＴを実行するためだけであれば２7＝１２８であるので２８個のゼロ値データを付加すればよいだけであるが、所望の周波数分解能を満たすようにするために、１９４８個もののゼロ値データを付加しているのである。
【００６８】
以上、本発明の好ましい実施形態について、反射エコーにおけるドップラー周波数の検出を例に挙げて説明したが、本発明の適用範囲はこれに限られるものではない。検出対象とすべき周波数範囲を想定できる信号から短い時間範囲の部分を切り出してＦＦＴによる解析を行う場合に、本発明は広く一般的に適用可能である。例えば、バンドパスフィルター等のフィルターを通して所望する周波数帯域の中の信号へと変換し、その後、短い時間信号を切りだして周波数検出を行うことに対しても、本発明は適用可能である。
【符号の説明】
【００６９】
１    送受波器
２    送受信切替回路
３    送信回路
４    受信増幅器
５    変調器
６    局部発振回路
７    アナログフィルタ
１０    アナログ／デジタル変換器（Ａ／Ｄ）
１１    信号切り出しゲート
１２    ゼロ付加部
１３    ＦＦＴ処理部
【図１】
2012247304.tif 000008
【図２】
2012247304.tif 000009
【図３】
2012247304.tif 000010
【図４】
2012247304.tif 000011
【図５】
2012247304.tif 000012
【図６】
2012247304.tif 000013
【図７】
2012247304.tif 000014
【図８】
2012247304.tif 000015