時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する方法であって、前記信号を前記所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング段階と、必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、該基準サンプル数に達するように、前記デジタルデータ列にゼロデータを付加する付加段階と、前記ゼロデータが付加された前記デジタルデータ列に対して高速フーリエ変換を行うFFT段階と、を有し、前記高速フーリエ変換の結果において最大値を示す周波数を前記ピークパワー周波数として検出する方法。前記基準サンプル数は、前記サンプリング段階で得られる前記デジタルデータ列のサンプル数の2倍以上である、請求項1に記載の方法。前記信号は、所定の周波数帯域を有し、前記所定の周波数帯域が低周波側に変換されるように、前記サンプリング段階の前に、前記信号に対して周波数変換を行う段階を有する、請求項1または2に記載の方法。時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する装置であって、前記信号を前記所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング手段と、必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、該基準サンプル数に達するように、前記デジタルデータ列にゼロデータを付加するゼロ付加手段と、前記ゼロデータが付加された前記デジタルデータ列に対して高速フーリエ変換を行うFFT手段と、を有する装置。前記基準サンプル数は、前記サンプリング手段で得られる前記デジタルデータ列のサンプル数の2倍以上である、請求項4に記載の装置。前記サンプリング手段は、前記信号に対してアナログ/デジタル変換を行うA/D変換器と、前記所望の時間幅に応じ、前記A/D変換器から出力されるデジタルデータ列を切り出す信号切り出しゲートと、を有する、請求項4または5に記載の装置。前記信号は所定の周波数帯域を有しており、さらに、前記サンプリング手段の前に、前記所定の周波数帯域が低周波側に変換されるように前記信号に対して周波数変換を行う周波数変換手段を有する、請求項4乃至6のいずれか1項に記載の装置。上述したように、ドップラー計測などにフーリエ変換を適用する場合、周波数分解能を向上させることと測定対象とする範囲の時間あるいは距離についての分解能を向上させることとが両立しない、というフーリエ変換の基本的な問題に起因する問題が生じる。この問題を回避するために、所望の短い時間幅で信号を切り出した上でFFTを行い、そのFFT結果に対して各種の補間方法や推定方法を適用して、FFTによって離散化した周波数間隔よりも細かくピークパワーの周波数を求めることも行われているが、これは、構成の肥大化やアルゴリズムの複雑化をもたらし、そのために演算速度の低下やコストアップがもたらされる。さらに、これらの補間方法や推定方法を用いた場合には、必ずしも所望の周波数精度が得られない、という問題点もある。本発明の目的は、フーリエ変換における上述した基本的な問題を解決し、測定対象とする時間範囲の分解能を高めつつ、周波数分解能を高めることができる検出方法及び装置を提供することにある。本発明の検出方法は、時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する方法であって、信号を所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング段階と、必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、この基準サンプル数に達するように、デジタルデータ列にゼロデータを付加する付加段階と、ゼロデータが付加されたデジタルデータ列に対して高速フーリエ変換を行うFFT段階と、を有し、高速フーリエ変換の結果において最大値を示す周波数をピークパワー周波数として検出する。本発明の検出装置は、時間変化する信号における所望の時間幅の範囲内でのピークパワー周波数を検出する装置であって、信号を所望の時間幅内においてサンプリングしてデジタルデータ列とするサンプリング手段と、必要とされる周波数分解能を満たすために高速フーリエ変換において必要なサンプル数を基準サンプル数として、この基準サンプル数に達するように、デジタルデータ列にゼロデータを付加するゼロ付加手段と、ゼロデータが付加されたデジタルデータ列に対して高速フーリエ変換を行うFFT手段と、を有する。本発明では、粗い周波数分解能でのFFTの結果から補間や推定によってピークパワー周波数を求めるのではなく、所望の周波数分解能が得られるようにゼロデータを付加してからFFTを行うので、FFT後の補間や推定などを行うことなく、所望の周波数分解能でFFTからピークパワー周波数を直接求めることができるようになる。本発明は、複雑な回路やアルゴリズムを必要としないので、シンプルな構成で、低コストで、精度良く、かつ高速に、ピークパワー周波数を求めることができる。図7は、本発明の実施の一形態のピークパワースペクトル検出方法の実施に適したドップラー計測機器の構成を示している。このドップラー計測機器は、本発明に基づくピークパワースペクトル検出装置の一例である。図7に示したドップラー計測機器は、図2に示したドップラー計測機器において、信号切り出しゲート11とFFT処理部13との間に、ゼロ付加部12を挿入したものである。信号切り出しゲート11からは、測定対象の時間範囲に対応して信号をサンプリングしデジタル化して得られるデジタルデータ列が出力されるが、ゼロ付加部12は、必要とされる周波数分解能を満たすためにFFTにおいて必要なサンプル数を基準サンプル数として、この基準サンプル数に達するように、デジタルデータ列にゼロデータを付加するものである。ゼロ付加部12での処理の詳細については後述する。図7に示すドップラー計測機器では、図2に示したものと同様に、変調器5及び局部発振回路6を使用して、送受波器1で受信した信号における上述の周波数範囲Δfpを低周波数側に周波数変換している。これは、高い周波数成分を検出するためにFFTを行うためには、サンプリング周波数を高くする必要があり、サンプリング周波数が高い分だけサンプル数が膨大となり過ぎるからである。ドップラー計測の場合、例えば、120kHz±1.2kHzの周波数範囲内でドップラー周波数を検出すればよく、検出対象とすべき周波数範囲Δfpが決まっているから、FFTにおいてサンプリング周波数を低下させ、かつ、必要な分解能の範囲内でサンプル数を減らすために、本実施形態では、受信信号を元の周波数よりも低い周波数の中間周波数信号に周波数変換している。前述の式2-1~2-3などを用いて示したように、周波数変換後の中間周波数信号でも、周波数範囲Δfpは保たれている。変調器5の出力から低周波側の中間周波数成分だけを取り出すことと、後段のアナログ/デジタル変換においてエリアシングが生じないようにすることのために、アナログフィルタ7が設けられている。アナログ/デジタル変換器A/D10は、FFTの実行が可能なサンプリング周波数でアナログフィルタ7の出力をサンプリングしてアナログ/デジタル変換する。これにより、A/D変換器10から、中間周波数信号における瞬時値を表すデジタル値からなるデジタルデータ列が出力されることになる。信号切り出しゲート11は、所望する時間範囲のデジタルデータ列を切り出す。この切り出されたデータのサンプル数を切り出しサンプル数とする。ゼロ付加部12は、要求されている周波数分解能に対応するFFTのサンプル数を基準サンプル数として求め、基準サンプル数に達するように、基準サンプル数と切り出しサンプル数との差に相当する数のゼロデータ値がゼロであるサンプルを切り出されたデジタルデータ列に付加する。そして、ゼロデータが付加されてサンプル数が基準サンプル数とされたデジタルデータ列に対し、FFT処理部13がFFTを実行し、ピークパワースペクトルを算出する。算出されたピークパワースペクトルは、信号切り出しゲート11において切り出された範囲に対するものであって、信号におけるその時間範囲でのピークパワー周波数を示している。また、ゼロデータが付加されたことにより、信号においてFFTの対象となる時間幅が実質的に拡大されており、この時間幅の逆数が周波数分解能に該当することから、周波数分解能が向上している。このように、本実施形態によれば、受信信号の中での所望の時間範囲あるいは距離範囲のピークパワースペクトルを、FFTの演算だけで、所望の周波数精度を満たして検出することができる。以下、数値例を挙げて、この実施形態における周波数分解能について説明する。背景技術の欄で述べたものと同じ条件で、超音波のパルス状の信号を海中に送信し、移動物からの反射エコーのドップラー周波数を測定する場合を考える。すなわち、海水中の音波伝搬速度CをC=1500m/sとし、送信周波数ftxをftx=120kHzとし、検出最大速度水平方向VをV=15m/sとし、音波を水平方向から斜め方向θ=60°に送受信するものとする。このとき、ドップラー信号は120±1.2kHzの範囲であり、速度精度として0.1m/sを要求すれば周波数分解能f0は12Hzとなる。また、ドップラー計測において検出対象とする位置の範囲位置分解能を7.5mとし、周波数変換後の信号の中心周波数すなわち中間周波数fmid50を5kHzとし、アナログ/デジタル変換でのサンプリング周波数及びFFTのサンプリング周波数fftを20kHzとする。周波数分解能f0から必要なT0を求めると、T0=1/f0=1/12≒83.3msとなる。海中に音波パルスを発射して受信信号を観測し、その受信信号を中間周波数信号に変換した後の波形を示す図4において、上述したように、符号16は反射エコーを表し、符号17は送信信号自体を表し、符号18は海底エコーを表している。海底エコー18の周波数を解析し、式1-1に基づき、ドップラー周波数から船の速度を求める場合を考える。ここでは、送信信号などの不要な信号が含まれないようにするために、海底エコー18の部分だけを切り出してFFTを行い、海底からの反射波に関してドップラー周波数を求める。図8は、図4に示す信号波形から海底エコーだけを切り出した信号19を示している。ここでは、海底エコー19の部分として、音波パルスの送信から時間t226の経過後に時間幅T224を切り出すものとする。ここで切り出された信号は、サンプル数としては100個である。時間幅T2の100個のサンプルのみを用いてフーリエ変換を行った場合の周波数分解能f0は、f0=20000Hz/100=200Hzとなり、要求されている周波数分解能である12Hzを全く満たさない。ここで、要求されている周波数分解能である12Hzをフーリエ変換が満たすために必要なサンプル数Nを求めると、式1-7から、N=fft/f0=20000/12≒1667となる。1667は2のべき乗数ではないので、1667を超えて最も1667に近い2のべき乗数を考えると、211=2048が見つかる。そこで、FFTにおけるサンプル数Nとして2048を採用すると、式1-6より、周波数分解能f0は、f0=20000Hz/2048=9.766Hzとなって、要求事項を満足している。そこでこの実施形態では、FFTをサンプル数N=2048で実行するものとし、海底エコーの信号19の分のサンプル数は100なので、不足した1948個のサンプルについては、ゼロ値からなるサンプルを挿入する。図8においてT023は、2048個のサンプルに対応する時間であり、FFT解析の対象となる時間幅である。ここでT0=T2T3の関係が成り立っており、測定したい範囲の信号19を切り出してFFT解析に必要なサンプル数Nを満たすように付加されるゼロ値信号の継続時間が、T3に対応することになる。時間幅T3の挿入位置は、時間幅T2すなわち切り出された信号19の前でも後でもよい。あるいは、信号19を挟むように前後に分割してゼロ値信号を挿入する時間幅を配置してFFT解析の時間幅T027を設定してもよい。ここで示す例では、音波パルスの送信から信号の切り出しを開始するまでの経過時間t2図8での符号26と時間幅T2図8での符号24を予め定めているが、海底エコーのような強度の大きな反射エコーを検出対象とする場合であれば、切り出しまでの経過時間t2も切り出しの時間幅T2も自動検出して決定することも可能である。深度ごとの潮流の流速を求める場合などには、自動検出を行うことができないのでt2及びT2について予め設定しておくことになる。ここでゼロ付加部12によるゼロ値データの付加について説明する。一般にFFTではその対象とするサンプル数は2のべき乗である必要がある。その一方で、FFTによる解析を行いたい対象データにおけるサンプル数が2のべき乗となっていないことはよくあることである。そこで、サンプル数が2のべき乗となるように、ゼロデータのサンプルを付加することがゼロパディングなどと称して行われている。その場合、とにかくFFTを行えるだけのサンプル数となればよいのであるから、対象データのサンプル数をNsとして、2m1<Ns≦2mを満たす自然数mが存在する時に、2m-Ns個のゼロ値データを付加することになる。それに対して本実施形態では、単にFFTの実行が可能になるようにゼロ値データを付加するのではなく、FFTにおけるサンプリング周波数fftを考慮して所望の周波数分解能が得られるように、ゼロ値データを付加する。したがって、付加されるゼロ値データの数も、単にFFTを実行可能になるようにする場合に比べ、はるかに多くなることがある。上述した例では、海底エコーを切り出してサンプル数が100になったとして、単にFFTを実行するためだけであれば27=128であるので28個のゼロ値データを付加すればよいだけであるが、所望の周波数分解能を満たすようにするために、1948個もののゼロ値データを付加しているのである。以上、本発明の好ましい実施形態について、反射エコーにおけるドップラー周波数の検出を例に挙げて説明したが、本発明の適用範囲はこれに限られるものではない。検出対象とすべき周波数範囲を想定できる信号から短い時間範囲の部分を切り出してFFTによる解析を行う場合に、本発明は広く一般的に適用可能である。例えば、バンドパスフィルター等のフィルターを通して所望する周波数帯域の中の信号へと変換し、その後、短い時間信号を切りだして周波数検出を行うことに対しても、本発明は適用可能である。