JP5035815B2 - 周波数測定装置 - Google Patents

Description

本発明は、ＦＦＴ（Fast Fourier Transform）を用いて周波数を測定する周波数測定装置に関する。

従来、ＦＦＴを用いた周波数測定では、見かけの演算データ数を増加し分解能を高める“零詰め方式”が知られている。この方式によれば、最も正確に周波数を求めることができるものの、処理能力と演算精度で制限され、極端に分解能を増やすことは困難である。

そこで、近似曲線や補間法を用いて周波数を測定することが古くから知られている。例えば、補間法を用いたものとして特許文献１に示すものがある。
特開平１０−２１３６１３号公報

上記した近似曲線や補間法を用いたものは、元のアナログ信号に含まれないような側波帯成分（サイドローブ）を抑圧するためにＡ／Ｄ変換によって得られたデータにハニング窓と称する窓関数をかけてＦＦＴ演算を行い、また複雑な補間式と組み合わせて補間を行うというもので、周波数を求めるのに複雑な過程を必要とする。

本発明は上記問題に鑑みたもので、ＦＦＴを用いた周波数測定を簡易に行うことを目的とする。また、その場合に、精度よく周波数を求めることを他の目的とする。

上記目的を達成するため、請求項１に記載の周波数測定装置では、アナログ信号をディジタルのデータに変換するＡ／Ｄ変換器と、
前記Ａ／Ｄ変換器によって得られたデータに対して、高速フーリエ演算を行う高速フーリエ演算手段と、
前記高速フーリエ演算手段の演算結果をスカラー変換して前記アナログ信号のスペクトラムのデータを求めるベクトル・スカラー変換手段と、を備え、
前記Ａ／Ｄ変換器によって得られたデータにサイドローブを抑圧するための窓関数をかけずに前記高速フーリエ演算が行われ、これにより前記ベクトル・スカラー変換手段によって求められたスペクトラムのデータは、サイドローブが抑圧されていない波形のデータとなるように構成されており、
さらに、前記ベクトル・スカラー変換手段によって求められたスペクトラムのデータの中から、所定次数のピークにおける最大のデータと、周波数軸上でその前後にあるデータとを検出し、さらに前記前後にあるデータのうち大きい方のデータを検出するデータ検出手段と、
前記検出された最大のデータと、前記検出された前後にあるデータのうち大きい方のデータとによる両データの周波数の間において、前記両データの周波数の差を、前記両データのうち前記最大のデータに近い方に按分した値となるように前記前後にあるデータの比により按分したところの周波数を前記アナログ信号の周波数として求める周波数演算手段と、を備えたことを特徴としている。

この発明によれば、周波数測定を簡易に行うことができるとともに、精度よく周波数を求めることができる。

この場合、請求項２に記載の発明のように、上記した両データを、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）とし、前記両データの周波数の差をΔｆとしたとき、前記周波数演算手段は、前記アナログ信号の周波数ｆｘを、ｆｘ＝ｆ（ｎ）＋Δｆ・ｙ（ｎ＋１）／｛ｙ（ｎ）＋ｙ（ｎ＋１）｝あるいはｆｘ＝ｆ（ｎ＋１）−Δｆ・ｙ（ｎ）／｛ｙ（ｎ）＋ｙ（ｎ＋１）｝により求めるようにすることができる。また、請求項３に記載の発明のように、前記周波数演算手段は、前記前後にあるデータのうち大きい方のデータが所定値以下となるとき、あるいは前記前後にあるデータの差がある値以下のときには、前記周波数を求める演算を行わずに、前記最大のデータの周波数を前記アナログ信号の周波数とするものとすることができる。さらに、請求項４に記載の発明のように、前記周波数演算手段は、前記最大のデータに対する前記前後にあるデータの比率がそれぞれ所定値を超えるときには、前記前後にあるデータの周波数の間において、前記前後にあるデータの周波数の差を、前記前後にあるデータのうち大きい方のデータに近い方に按分した値となるように前記前後にあるデータの比により按分したところの周波数を前記アナログ信号の周波数として求めるようにすることもできる。

また、請求項５に記載に記載の周波数測定装置では、前記ベクトル・スカラー変換手段によって求められたスペクトラムのデータの中から、所定次数のピークにおける最大のデータと、周波数軸上でその前後にあるデータとを検出するデータ検出手段と、
前記データ検出手段にて検出されたデータのうち前記最大のデータを外した前記前後にあるデータの周波数の間において、前記前後にあるデータの周波数の差を、前記前後にあるデータのうち大きい方のデータに近い方に按分した値となるように前記前後にあるデータの比により按分したところの周波数を前記アナログ信号の周波数として求める周波数演算手段と、を備えたことを特徴としている。

この発明によれば、周波数測定を簡易に行うことができるとともに、中心となるスペクトルのデータに対し、両サイドのスペクトラムのデータが十分小さな値とならない場合において、精度よく周波数を求めることができる。

この場合、請求項６に記載の発明のように、前記最大のデータを、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）、前記前後にあるデータを、周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）とし、各データの周波数の差をΔｆとしたとき、前記周波数演算手段は、前記アナログ信号の周波数ｆｘを、ｆｘ＝ｆ（ｎ−１）＋２Δｆ・ｙ（ｎ＋１）／｛ｙ（ｎ−１）＋ｙ（ｎ＋１）｝あるいはｆｘ＝ｆ（ｎ＋１）−２Δｆ・ｙ（ｎ−１）／｛ｙ（ｎ−１）＋ｙ（ｎ＋１）｝により求めるようにすることができる。

図１に本発明の一実施形態に係る周波数測定装置の構成を示す。この図１に示す周波数測定装置には、入力端子からアナログ信号Ｓ（ｔ）が入力される。このアナログ信号Ｓ（ｔ）は、持続時間の限られた低周波信号（例えば、弦楽器や打楽器など、数ＫＨｚ以下で持続時間が限られた信号）である。アナログ信号Ｓ（ｔ）が減衰振動をする信号である場合、振幅変化が生じ誤差が発生するので、この実施形態では、ＡＬＣ（Automatic Level Control）回路１０により、アナログ信号Ｓ（ｔ）の振幅を一定に保つようにしている。

このＡＬＣ回路１０によって振幅が一定に保たれたアナログ信号Ｓ（ｔ）は、Ａ／Ｄ変換器１１によって所定のサンプリング周波数ｆｓでサンプリングされ、ディジタルのデータに変換され、サンプリングデータメモリ１２に記憶される。

サンプリングデータメモリ１２に記憶されたデータは、ＦＦＴ演算手段１４に入力される。なお、このＦＦＴ演算手段１４、および後述するベクトル・スカラー変換手段１５、ピーク検出手段１７、データ検出手段１８、周波数演算手段１９は、コンピュータを用いてソフト的に演算処理される部分であり、上記したＡＬＣ回路１０は、ＦＦＴ演算手段１４の前段においてＡＬＣ手段１３としてソフト的に行うようになっていてもよい。

ＦＦＴ演算手段１４の演算結果は、ベクトル・スカラー変換手段１５によってスカラー変換され、アナログ信号のスペクトラムのデータとしてスペクトラムデータメモリ１６に記憶される。

このデータは、図２に示すように、横軸にサンプリング周波数ｆｓの１／２、縦軸にスペクトラムの大きさ（マグニチュード）を取った場合、高調波の１次、２次、３次、４次、５次…のピークが現れるデータとなる。また、それぞれのピークにおいて、図３に示すように、周波数ｆ（ｎ−１）、ｆ（ｎ）、ｆ（ｎ＋１）、ｆ（ｎ＋２）…におけるそれぞれのスペクトラムのデータが得られる。なお、図３は、縦軸を圧縮して模式的に表現したものとなっている。また、図では、周波数ｆ（ｎ）におけるスペクトラムのデータをｙ（ｎ）、周波数ｆ（ｎ＋１）におけるスペクトラムのデータをｙ（ｎ＋１）として示している。また、図３に示されるデータは、Ａ／Ｄ変換器１１でのサンプリング周波数を１０ＫＨｚ、サンプリング数Ｎを１０２４とし、スペクトラム間隔Δｆを１０Ｈｚとしたときのものである。

この実施形態では、上記したＡ／Ｄ変換器１１、ＦＦＴ演算手段１４、ベクトル・スカラー変換手段１５およびスペクトラムデータメモリ１６として、背景技術の欄で示した特開平１０−２１３６１３号公報に記載のものと同じものを用いることができる。しかし、特開平１０−２１３６１３号公報に記載のもののように、ＦＦＴ演算手段１４の前段に窓関数演算手段を用いていない。このため、ベクトル・スカラー変換手段１５によってスカラー変換されたスペクトラムのデータは、サイドローブが抑圧されていない波形のデータとなる。すなわち、アナログ信号の周波数がスペクトラム間隔Δｆの整数倍の場合には、データの取り込み区間とアナログ信号の周期が一致し、該当周波数に一本のスペクトル（メインローブ）が立つ波形のデータとなるが、アナログ信号の周波数がスペクトル間隔の整数倍でない場合には、メインローブ以外にサイドローブが立った波形のデータとなる。したがって、図３では、サイドローブを示すデータとして、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）と、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）が示されている。

ここで、アナログ信号の周波数は、その波形におけるピーク値ｙｍとなるときの周波数ｆｘである。この実施形態では、上記した周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）と、スペクトラム間隔Δｆすなわち周波数ｆ（ｎ＋１）−ｆ（ｎ）から、下記の数式１を用いてアナログ信号の周波数ｆｘを求める。
（数１）
ｆｘ＝ｆ（ｎ）＋Δｆ・ｙ（ｎ＋１）／｛ｙ（ｎ）＋ｙ（ｎ＋１）｝
この数式１は、以下のようにして導かれる。ｙ（ｎ）とｙ（ｎ＋１）を結ぶ線上において、周波数ｆｘとなる点がｘ１：ｘ２で与えられるとし、ｘ１：ｘ２がｙ（ｎ＋１）：ｙ（ｎ）となる、つまりｙ（ｎ）とｙ（ｎ＋１）を結ぶ線上においてｙ（ｎ）とｙ（ｎ＋１）のうち大きい方のデータｙ（ｎ）に近い方に周波数ｆｘの点を示すようにｙ（ｎ）とｙ（ｎ＋１）で按分（比例配分）する。そして、ｘ１：ｘ２＝ｙ（ｎ＋１）：ｙ（ｎ）およびΔｆ＝ｘ１＋ｘ２から、ｘ１＝Δｆ・ｙ（ｎ＋１）／｛ｙ（ｎ）＋ｙ（ｎ＋１）｝が得られ、それにｆ（ｎ）を加えることにより、周波数ｆｘが求められる。つまり、上記した数式１の意味するところは、周波数ｆ（ｎ）と周波数ｆ（ｎ＋１）の間において、スペクトラム間隔Δｆを、ｙ（ｎ）とｙ（ｎ＋１）のうち大きいデータ（図３の場合ｙ（ｎ））に近い方に按分した値となるようにｙ（ｎ）とｙ（ｎ＋１）の比により按分したところの周波数ｆｘをアナログ信号の周波数とするものである。

そこで、上記のようにしてアナログ信号の周波数ｆｘを求めるため、この実施形態では、図１に示すように、まず、ピーク検出手段１７により、図２に示すピークが検出される。この場合、ピークは、スペクトラムデータメモリ１６に記憶されたスペクトラムデータを周波数の低い方から順に比較していくことによりデータが増加から減少に変化することでピークを求めることができる。最初のピーク検出においては１次のピークが検出される。次に、データ検出手段１８により、まず１次のピークについてスペクトラムデータメモリ１６に記憶されたスペクトラムデータの中から最大のデータが検出される。この最大のデータは、図３の場合、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）である。さらに、このデータ検出手段１８により、周波数軸上で最大のデータの前後にあるデータのうち大きい方のデータが検出される。図３の場合、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）である。

このようにして検出された最大のデータおよびその前後にあるデータのうち大きい方のデータ、さらに最大のデータの周波数とスペクトラム間隔Δｆとにより、周波数演算手段１９において、数式１を用いてアナログ信号の周波数ｆｘが求められる。

また、１次のピークにおける周波数が測定された後、２次以降のピークにおける周波数についても、上記したピーク検出手段１７、データ検出手段１８および周波数演算手段１９の処理を繰り返すことにより求められる。

上記した実施形態によれば、従来の近似曲線や補間法を用いたもののように窓関数をかけてＦＦＴ演算を行いサイドローブを抑圧することを行っていないので、窓関数演算が不要であり、また数式１に示すように、求める周波数ｆｘの前後の２つのスペクトラムデータの比によりスペクトラム間隔Δｆを按分して周波数ｆｘを求めているので、極めて簡単に周波数ｆｘを求めることができる。

なお、上記した周波数測定装置を用い、発振器の信号についてその周波数を測定したところ、入力周波数＝発振器の表示周波数＋サンプリング周波数オフセット（実験データの平均から求めたもの）＝２００２．９５＋０．０４８＝２００２．９９９８に対して、２０００＋１０×３６８．８５／（８６１．９０＋３６８．８５）＝２００２．９９６９５３０７７（単位は全てＨｚ）の結果が得られ、極めて精度よく周波数を測定することができた。

このように高精度で簡易な方法で周波数を測定できる周波数測定装置は、弦楽器や打楽器などの波形解析やプロトン磁力計による磁場測定など、持続時間の限られた信号の周波数測定を行う場合に有効である。

なお、上記した実施形態において、アナログ信号が減衰しない一定振幅の信号であって、それを測定する場合には、ＡＬＣ回路１０、ＡＬＣ手段１３をなくしてもよい。

また、数式１では、２つのスペクトラムのうち周波数の低い方ｆ（ｎ）を基準として周波数ｆｘを求めるものを示したが、周波数の高い方ｆ（ｎ＋１）を基準として周波数ｆｘを求めるようにしてもよい。この場合、周波数ｆｘは数式２を用いて求めることができる。
（数２）
ｆｘ＝ｆ（ｎ＋１）−Δｆ・ｙ（ｎ）／｛ｙ（ｎ）＋ｙ（ｎ＋１）｝
また、アナログ信号の周波数がスペクトラム間隔Δｆの整数倍となる場合には、一本のスペクトルが立つ波形のデータとなり、その前後の周波数のスペクトルのデータは理論的にはゼロであるが実際には演算誤差等によってゼロにはならないので、そのデータを用いると周波数の測定精度が低下することになる。この場合、演算誤差等によって、前後の周波数のスペクトルのデータはそれぞれ極めて小さな値となるか、あるいは両者がほぼ等しい値になる。したがって、周波数演算手段１９では、前後の周波数のスペクトルのデータのうち大きい方のデータが所定値（スレッショルド値）以下となるとき、あるいは前後の周波数のスペクトルのデータの差がほぼ等しいとき（すなわち差がある値以下のとき）に、上記した周波数演算を行わずに、最大のデータの周波数をｆｘとして出力し、前後の周波数のスペクトルのデータのうち大きい方のデータが所定値より大きく、かつ前後の周波数のスペクトルのデータの差がある値より大きい場合に、上記した周波数演算を行って周波数ｆｘを出力するようにするのが好ましい。

また、上記した実施形態では、サイドローブを示すデータとして、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）と、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）を用い、それによりスペクトラム間隔Δｆを按分してアナログ信号の周波数ｆｘを求めるものを示したが、サイドローブを示すデータとして３つ以上を用い、所定の関係式を用いてアナログ信号の周波数ｆｘ、すなわちメインローブの周波数ｆｘを求めるようにしてもよい。

サイドローブを示すデータを３つ以上用いてメインローブの周波数ｆｘを求める場合の一例について以下説明する。例えば、アナログ信号の周波数がスペクトラム間隔Δｆの整数倍と接近した時で、演算誤差やサンプリング誤差により、図４に示すように、中心となるスペクトルのデータｙ（ｎ）に対し、両サイドのスペクトラムのデータｙ（ｎ＋１）、ｙ（ｎ−１）が十分小さな値とならないことがある。また、ＡＬＣ回路が無くてサンプリング期間中の最初と最後の信号の振幅変化が大きい場合も、図４に示すような関係になる場合がある。このような場合に、数式１あるいは数式２を用いて周波数ｆｘを求めると、その誤差が大きくなるので、この場合には、３つのサイドローブを示すデータのうち、中心となるサイドローブのデータを外し、両サイドのサイドローブのデータである、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）と周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）により、周波数の間隔を２Δｆとして、数式３あるいは数式４により、周波数ｆｘを求める。
（数３）
ｆｘ＝ｆ（ｎ−１）＋２Δｆ・ｙ（ｎ＋１）／｛ｙ（ｎ−１）＋ｙ（ｎ＋１）｝
（数４）
ｆｘ＝ｆ（ｎ＋１）−２Δｆ・ｙ（ｎ−１）／｛ｙ（ｎ−１）＋ｙ（ｎ＋１）｝
この例の場合、以下のようにして実現することができる。データ検出手段１８は、ベクトル・スカラー変換手段１５によって求められたスペクトラムのデータの中から、所定次数のピークにおける最大のデータ、すなわち周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）と、周波数軸上でその前後にある両データ、すなわち周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）と周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）を検出する。周波数演算手段１９は、検出されたデータのうち、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）を外し、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）と周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）とスペクトラム間隔Δｆとにより、数式３あるいは数式４を用いてアナログ信号の周波数ｆｘを求める。この場合、周波数演算手段１９において、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）に対し、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）と周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）の比率がそれぞれ所定値（例えば１／１０００）を超えるときに、上記した数式３あるいは数式４を用いて周波数ｆｘを求め、それ以外のときには、上記した数式１あるいは数式２を用いて周波数ｆｘを求めるようにすることができる。

また、３つよりも多い数のサイドローブのデータを用いて周波数ｆｘを求めるようにしてもよい。例えば、５つのサイドローブを示すデータを用い、そのうち中心となるサイドローブのデータを外し、両サイドの４つのサイドローブのデータにより、周波数ｆｘを求めるようにしてもよい。この場合、内側２つのサイドローブに対しては、数式３あるいは数式４を用いて周波数ｆｘ’を求め、外側２つのサイドローブに対しては、周波数の間隔を４Δｆとし数式３あるいは数式４と同様の考えた方で周波数ｆｘ’’を求め、それぞれの周波数のいずれか１つを所定の要件で選択することによって周波数ｆｘを求める、あるいは両者を加重平均などすることによって周波数ｆｘを求めるようにしてもよい。

また、上記した例に対し、数式１、２を用いずに、数式３、４を用いて周波数ｆｘを求めるようにしてもよい。この場合、データ検出手段１８は、ベクトル・スカラー変換手段１５によって求められたスペクトラムのデータの中から、最大のデータを外して、周波数軸上でその前後にある両データ、すなわち周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）と周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）を検出し、周波数演算手段１９は、検出された、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）と周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）とスペクトラム間隔Δｆとにより、数式３あるいは数式４を用いてアナログ信号の周波数ｆｘを求める。

本発明の一実施形態における周波数測定装置の構成を示す図である。 アナログ信号のスペクトラムのデータを示す図である。 図２に示すそれぞれのピークの波形の詳細を示す図である。 本発明の他の実施形態におけるピークの波形の詳細を示す図である。

符号の説明

１０…ＡＬＣ回路、１１…Ａ／Ｄ変換器、１２…サンプリングデータメモリ、１４…ＦＦＴ演算手段、１５…ベクトル・スカラー変換手段、１６…スペクトラムデータメモリ、１７…ピーク検出手段、１８…データ検出手段、１９…周波数演算手段。

Claims (6)

1. アナログ信号をディジタルのデータに変換するＡ／Ｄ変換器と、
   前記Ａ／Ｄ変換器によって得られたデータに対して、高速フーリエ演算を行う高速フーリエ演算手段と、
   前記高速フーリエ演算手段の演算結果をスカラー変換して前記アナログ信号のスペクトラムのデータを求めるベクトル・スカラー変換手段と、を備え、
   前記Ａ／Ｄ変換器によって得られたデータにサイドローブを抑圧するための窓関数をかけずに前記高速フーリエ演算が行われ、これにより前記ベクトル・スカラー変換手段によって求められたスペクトラムのデータは、サイドローブが抑圧されていない波形のデータとなるように構成されており、
   さらに、前記ベクトル・スカラー変換手段によって求められたスペクトラムのデータの中から、所定次数のピークにおける最大のデータと、周波数軸上でその前後にあるデータとを検出し、さらに前記前後にあるデータのうち大きい方のデータを検出するデータ検出手段と、
   前記検出された最大のデータと、前記検出された前後にあるデータのうち大きい方のデータとによる両データの周波数の間において、前記両データの周波数の差を、前記両データのうち前記最大のデータに近い方に按分した値となるように前記両データの比により按分したところの周波数を前記アナログ信号の周波数として求める周波数演算手段と、を備えたことを特徴とする周波数測定装置。
2. 前記両データを、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）とし、前記両データの周波数の差をΔｆとしたとき、前記周波数演算手段は、前記アナログ信号の周波数ｆｘを、ｆｘ＝ｆ（ｎ）＋Δｆ・ｙ（ｎ＋１）／｛ｙ（ｎ）＋ｙ（ｎ＋１）｝あるいはｆｘ＝ｆ（ｎ＋１）−Δｆ・ｙ（ｎ）／｛ｙ（ｎ）＋ｙ（ｎ＋１）｝により求めることを特徴とする請求項１に記載の周波数測定装置。
3. 前記周波数演算手段は、前記前後にあるデータのうち大きい方のデータが所定値以下となるとき、あるいは前記前後にあるデータの差がある値以下のときには、前記周波数を求める演算を行わずに、前記最大のデータの周波数を前記アナログ信号の周波数とすることを特徴とする請求項１または２に記載の周波数測定装置。
4. 前記周波数演算手段は、前記最大のデータに対する前記前後にあるデータの比率がそれぞれ所定値を超えるときには、前記前後にあるデータの周波数の間において、前記前後にあるデータの周波数の差を、前記前後にあるデータのうち大きい方のデータに近い方に按分した値となるように前記前後にあるデータの比により按分したところの周波数を前記アナログ信号の周波数として求めることを特徴とする請求項１または２に記載の周波数測定装置。
5. アナログ信号をディジタルのデータに変換するＡ／Ｄ変換器と、
   前記Ａ／Ｄ変換器によって得られたデータに対して、高速フーリエ演算を行う高速フーリエ演算手段と、
   前記高速フーリエ演算手段の演算結果をスカラー変換して前記アナログ信号のスペクトラムのデータを求めるベクトル・スカラー変換手段と、を備え、
   前記Ａ／Ｄ変換器によって得られたデータにサイドローブを抑圧するための窓関数をかけずに前記高速フーリエ演算が行われ、これにより前記ベクトル・スカラー変換手段によって求められたスペクトラムのデータは、サイドローブが抑圧されていない波形のデータとなるように構成されており、
   さらに、前記ベクトル・スカラー変換手段によって求められたスペクトラムのデータの中から、所定次数のピークにおける最大のデータと、周波数軸上でその前後にあるデータとを検出するデータ検出手段と、
   前記データ検出手段にて検出されたデータのうち前記最大のデータを外した前記前後にあるデータの周波数の間において、前記前後にあるデータの周波数の差を、前記前後にあるデータのうち大きい方のデータに近い方に按分した値となるように前記前後にあるデータの比により按分したところの周波数を前記アナログ信号の周波数として求める周波数演算手段と、を備えたことを特徴とする周波数測定装置。
6. 前記最大のデータを、周波数ｆ（ｎ）のときのスペクトラムのデータｙ（ｎ）、前記前後にあるデータを、周波数ｆ（ｎ−１）のときのスペクトラムのデータｙ（ｎ−１）、周波数ｆ（ｎ＋１）のときのスペクトラムのデータｙ（ｎ＋１）とし、各データの周波数の差をΔｆとしたとき、前記周波数演算手段は、前記アナログ信号の周波数ｆｘを、ｆｘ＝ｆ（ｎ−１）＋２Δｆ・ｙ（ｎ＋１）／｛ｙ（ｎ−１）＋ｙ（ｎ＋１）｝あるいはｆｘ＝ｆ（ｎ＋１）−２Δｆ・ｙ（ｎ−１）／｛ｙ（ｎ−１）＋ｙ（ｎ＋１）｝により求めることを特徴とする請求項５に記載の周波数測定装置。

Images