JP2022169984A - Measurement device, measurement method and program - Google Patents
Measurement device, measurement method and program Download PDFInfo
- Publication number
- JP2022169984A JP2022169984A JP2021075770A JP2021075770A JP2022169984A JP 2022169984 A JP2022169984 A JP 2022169984A JP 2021075770 A JP2021075770 A JP 2021075770A JP 2021075770 A JP2021075770 A JP 2021075770A JP 2022169984 A JP2022169984 A JP 2022169984A
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- digital signal
- band
- distribution data
- waveform
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 78
- 238000000691 measurement method Methods 0.000 title claims abstract description 4
- 238000005070 sampling Methods 0.000 claims abstract description 47
- 238000000605 extraction Methods 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 description 25
- 239000000284 extract Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000000034 method Methods 0.000 description 6
- 230000005856 abnormality Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
Description
本発明は、計測装置、計測方法及びプログラムに関する。 The present invention relates to a measuring device, measuring method and program.
振動の計測結果に基づいて軸受の異常を診断する装置(振動の計測装置)が、特許文献1に開示されている。図11は、計測装置の構成例を示す図である。図11に例示された計測装置は、計測対象(例えば、軸受)の振動の卓越周波数を計測する。計測装置は、ローパスフィルタと、アナログ・デジタル変換部と、高速フーリエ変換部とを備える。
A device (vibration measuring device) for diagnosing bearing abnormalities based on vibration measurement results is disclosed in
図12は、計測処理の例を示す図である。ローパスフィルタには、計測対象の振動の波形信号が入力される。ローパスフィルタは、入力された波形信号において、カットオフ周波数を超える帯域を制限する。このため、カットオフ周波数を超える卓越周波数を計測することはできない。図12では、カットオフ周波数は、所定のナイキスト周波数である。アナログ・デジタル変換部は、入力された波形信号(アナログ信号)に対してサンプリング処理を実行することよって、波形信号(デジタル信号)を生成する。高速フーリエ変換部は、波形信号(デジタル信号)に対して、高速フーリエ変換を実行する。 FIG. 12 is a diagram illustrating an example of measurement processing. A waveform signal of vibration to be measured is input to the low-pass filter. A low-pass filter limits the band exceeding the cutoff frequency in the input waveform signal. Therefore, it is not possible to measure dominant frequencies exceeding the cutoff frequency. In FIG. 12, the cutoff frequency is the predetermined Nyquist frequency. The analog/digital converter generates a waveform signal (digital signal) by performing sampling processing on the input waveform signal (analog signal). The fast Fourier transform unit performs fast Fourier transform on the waveform signal (digital signal).
このような計測装置では、計測対象の振動の卓越周波数(損傷周波数)が予め特定されている必要がある。なぜなら、標本化定理により、計測対象の振動の卓越周波数がナイキスト周波数以下となるように、アナログ・デジタル変換部のサンプリング周波数が予め定められる必要があるからである。サンプリング周波数は、ナイキスト周波数の2逓倍である。このため、アナログ・デジタル変換部のサンプリング周波数は損傷周波数に対して十分に高い必要がある。しかしながら、サンプリング周波数が十分に高いアナログ・デジタル変換部のコストは非常に高い。アナログ・デジタル変換部のコストの増加を抑制するためには、波形信号においてナイキスト周波数を超える卓越周波数が計測可能であることが望ましい。 In such a measuring device, the dominant frequency (damage frequency) of the vibration to be measured must be specified in advance. This is because, according to the sampling theorem, the sampling frequency of the analog/digital converter must be determined in advance so that the dominant frequency of the vibration to be measured is equal to or lower than the Nyquist frequency. The sampling frequency is twice the Nyquist frequency. Therefore, the sampling frequency of the analog/digital converter must be sufficiently high with respect to the damage frequency. However, the cost of an analog-to-digital converter with a sufficiently high sampling frequency is very high. In order to suppress an increase in the cost of the analog-to-digital converter, it is desirable to be able to measure dominant frequencies exceeding the Nyquist frequency in the waveform signal.
上記事情に鑑み、本発明は、波形信号においてナイキスト周波数を超える卓越周波数を計測することが可能である計測装置、計測方法及びプログラムを提供することを目的としている。 In view of the above circumstances, it is an object of the present invention to provide a measuring apparatus, a measuring method, and a program capable of measuring dominant frequencies exceeding the Nyquist frequency in waveform signals.
本発明の一態様は、入力されたアナログ信号の帯域を制限するフィルタと、帯域が制限された前記アナログ信号に基づく第1デジタル信号と帯域が制限された又は制限されていない前記アナログ信号に基づく第2デジタル信号とを、サンプリング周波数に基づいて生成するサンプリング実行部と、前記第1デジタル信号の周波数分布データと前記第2デジタル信号の周波数分布データとを、前記サンプリング周波数に基づくナイキスト周波数以下の帯域について生成する高速フーリエ変換部と、前記第1デジタル信号の周波数分布データにおけるピーク波形の有無と前記第2デジタル信号の周波数分布データにおけるピーク波形の有無との差分又は共通点を周波数ごとに抽出し、前記ナイキスト周波数を基準として、前記差分又は前記共通点に基づいてピーク波形の周波数を折り返す抽出部とを備える計測装置である。 One aspect of the present invention is a filter for band-limiting an input analog signal; a sampling execution unit that generates a second digital signal based on a sampling frequency; and a frequency distribution data of the first digital signal and a frequency distribution data of the second digital signal that are equal to or lower than the Nyquist frequency based on the sampling frequency. A fast Fourier transform unit that generates a band and extracts, for each frequency, a difference or a common point between the presence or absence of a peak waveform in the frequency distribution data of the first digital signal and the presence or absence of a peak waveform in the frequency distribution data of the second digital signal. and an extraction unit that folds the frequency of the peak waveform based on the difference or the common point with the Nyquist frequency as a reference.
本発明の一態様は、計測装置が実行する計測方法であって、入力されたアナログ信号の帯域を制限するフィルタステップと、帯域が制限された前記アナログ信号に基づく第1デジタル信号と帯域が制限された又は制限されていない前記アナログ信号に基づく第2デジタル信号とを、サンプリング周波数に基づいて生成するサンプリング実行ステップと、前記第1デジタル信号の周波数分布データと前記第2デジタル信号の周波数分布データとを、前記サンプリング周波数に基づくナイキスト周波数以下の帯域について生成する高速フーリエ変換ステップと、前記第1デジタル信号の周波数分布データにおけるピーク波形の有無と前記第2デジタル信号の周波数分布データにおけるピーク波形の有無との差分又は共通点を周波数ごとに抽出し、前記ナイキスト周波数を基準として、前記差分又は前記共通点に基づいてピーク波形の周波数を折り返す抽出ステップとを含む計測方法である。 One aspect of the present invention is a measurement method performed by a measurement device, comprising: a filter step for limiting the band of an input analog signal; a sampling execution step of generating a second digital signal based on said analog signal that is limited or unlimited based on a sampling frequency; frequency distribution data of said first digital signal and frequency distribution data of said second digital signal; a fast Fourier transform step of generating a band below the Nyquist frequency based on the sampling frequency; Presence or absence of a peak waveform in the frequency distribution data of the first digital signal; an extracting step of extracting a difference or a common point between presence and absence for each frequency, and folding back the frequency of the peak waveform based on the difference or the common point with the Nyquist frequency as a reference.
本発明の一態様は、上記の計測装置としてコンピュータを機能させるためのプログラムである。 One aspect of the present invention is a program for causing a computer to function as the above measuring device.
本発明により、波形信号においてナイキスト周波数を超える卓越周波数を計測することが可能である。 The present invention makes it possible to measure dominant frequencies above the Nyquist frequency in a waveform signal.
本発明の実施形態について、図面を参照して詳細に説明する。
(第1実施形態)
図1は、第1実施形態における、計測装置1aの構成例を示す図である。計測装置1aは、周波数を計測する装置である。計測装置1aは、切替部10と、ローパスフィルタ11と、アナログ・デジタル変換部12と、高速フーリエ変換部13と、計測部14aと、記憶装置15と、出力部16とを備える。なお、ローパスフィルタ11は、バンドパスフィルタでもよい。
Embodiments of the present invention will be described in detail with reference to the drawings.
(First embodiment)
FIG. 1 is a diagram showing a configuration example of a
計測装置1aの各機能部のうちの一部又は全部は、CPU(Central Processing Unit)等のプロセッサが、不揮発性の記録媒体(非一時的な記録媒体)を有する記憶部(例えば、記憶装置15)に記憶されたプログラムを実行することにより、ソフトウェアとして実現される。プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ROM(Read Only Memory)、CD-ROM(Compact Disc Read Only Memory)等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置などの非一時的な記録媒体である。
Some or all of the functional units of the
計測装置1aの各機能部のうちの一部又は全部は、例えば、LSI(Large Scale Integrated circuit)、ASIC(Application Specific Integrated Circuit)、PLD(Programmable Logic Device)又はFPGA(Field Programmable Gate Array)等を用いた電子回路(electronic circuit又はcircuitry)を含むハードウェアを用いて実現されてもよい。
Some or all of the functional units of the
センサ(不図示)又は測定装置(不図示)は、計測対象の振動(例えば、音)を測定する。計測対象とは、例えばポンプの回転軸の軸受である。センサは、例えば加速度センサである。測定装置は、例えば騒音計である。センサ又は測定装置は、計測された振動の波形信号を生成する。 A sensor (not shown) or measuring device (not shown) measures the vibration (eg, sound) to be measured. The object to be measured is, for example, the bearing of the rotary shaft of the pump. The sensor is, for example, an acceleration sensor. The measuring device is, for example, a sound level meter. A sensor or measurement device produces a waveform signal of the measured vibration.
切替部10は、センサ又は測定装置から、計測された振動の波形信号(アナログ信号)を取得する。切替部10は、ローパスフィルタ11及びアナログ・デジタル変換部12に、例えば所定周期で交互に波形信号を出力する。すなわち、切替部10は、波形信号の出力先を、例えば所定周期で交互に切り替える。所定周期は、例えば、計測対象の損傷周波数に基づいて定められる。所定周期は、例えば、測定対象の軸受における回転軸の回転周期に基づいて定められる。
The
ローパスフィルタ11は、入力された波形信号の帯域を制限する機能部である。ローパスフィルタ11は、波形信号が切替部10から入力された場合、カットオフ周波数以下の帯域の波形信号を、アナログ・デジタル変換部12に出力する。以下では、カットオフ周波数は、一例として、ナイキスト周波数に等しい。なお、カットオフ周波数は、計測対象の振動特性に応じて着目された周波数の帯域(例えば、損傷周波数を含む帯域)に合わせて定められてもよい。周波数の帯域が絞られることによって、表示処理及び演算処理を軽量化することが可能である。
The low-
ナイキスト周波数は、アナログ・デジタル変換部12が実行するアナログ・デジタル変換処理におけるサンプリング周波数の半分の周波数である。換言すれば、サンプリング周波数は、ナイキスト周波数の2逓倍である。
The Nyquist frequency is half the sampling frequency in the analog-to-digital conversion processing executed by the analog-to-
アナログ・デジタル変換部12には、例えば所定周期で、波形信号が切替部10及びローパスフィルタ11から交互に入力される。すなわち、第1期間において、波形信号が、アナログ・デジタル変換部12にローパスフィルタ11から入力される。第1期間とは異なる第2期間において、波形信号が、アナログ・デジタル変換部12に切替部10から入力される。
Waveform signals are alternately input to the analog-to-
アナログ・デジタル変換部12は、入力された波形信号(アナログ信号)に対して、アナログ・デジタル変換処理を実行する。すなわち、アナログ・デジタル変換部12(サンプリング実行部)は、入力された波形信号に対して、所定のサンプリング周波数でサンプリング処理を実行する。ここで、サンプリング周波数は、一例として100Hzである。したがって、ナイキスト周波数は、一例として50Hz(=100/2Hz)である。アナログ・デジタル変換部12は、生成された波形信号(デジタル信号)を、高速フーリエ変換部13に出力する。
The analog/
アナログ・デジタル変換部12にローパスフィルタ11から波形信号が入力されている場合、高速フーリエ変換部13は、ローパスフィルタ11によって帯域が制限された波形信号に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13は、フィルタによって帯域が制限された波形信号の周波数分布データ(以下「制限周波数データ」という。)を、計測部14aに出力する。
When a waveform signal is input from the low-
アナログ・デジタル変換部12に切替部10から波形信号が入力されている場合、高速フーリエ変換部13は、ローパスフィルタ11によって帯域が制限されていない波形信号(エイリアシングの影響を含む波形信号)に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13は、フィルタによって帯域が制限されていない波形信号の周波数分布データ(以下「非制限周波数データ」という。)を、計測部14aに出力する。
When a waveform signal is input from the switching
アナログ・デジタル変換部12にローパスフィルタ11から波形信号が入力されている場合、計測部14aは、制限周波数データを高速フーリエ変換部13から取得する。計測部14aは、制限周波数データを記憶装置15に記録する。
When a waveform signal is input from the low-
アナログ・デジタル変換部12に切替部10から波形信号が入力されている場合、計測部14aは、非制限周波数データを高速フーリエ変換部13から取得する。計測部14aは、非制限周波数データを記憶装置15に記録する。
When the waveform signal is input from the switching
計測部14a(抽出部)は、制限周波数データにおけるピーク波形と非制限周波数データにおけるピーク波形との有無の差分を、周波数ごとに抽出する。計測部14aは、非制限周波数データにおいて差分として抽出されたピーク波形の周波数(卓越周波数)を計測する。計測部14aは、ナイキスト周波数を基準として、差分として抽出されたピーク波形の周波数を折り返す。計測部14aは、制限周波数データと、折り返された周波数のピーク波形のデータとを、出力部16に出力する。
The measuring unit 14a (extracting unit) extracts, for each frequency, the difference between the peak waveform in the restricted frequency data and the peak waveform in the non-restricted frequency data. The measurement unit 14a measures the frequency (dominant frequency) of the peak waveform extracted as the difference in the unlimited frequency data. The measurement unit 14a repeats the frequency of the peak waveform extracted as the difference with reference to the Nyquist frequency. The measurement unit 14 a outputs the limit frequency data and the peak waveform data of the folded frequency to the
記憶装置15は、制限周波数データと非制限周波数データとを記憶する。記憶装置15は、折り返された周波数のピーク波形のデータを記憶してもよい。 The storage device 15 stores limited frequency data and non-limited frequency data. The storage device 15 may store the peak waveform data of the folded frequency.
出力部16は、制限周波数データと、折り返された周波数のピーク波形のデータとを、表示装置等の外部装置(不図示)に出力する。出力部16(表示部)は、制限周波数データに基づく周波数分布(周波数スペクトル)の画像と、折り返された周波数のピーク波形の画像とを、並べて表示してもよい。
The
次に、計測処理について説明する。
図2は、第1実施形態における、計測処理の例を示す図である。計測部14aは、制限周波数データ100におけるピーク波形と非制限周波数データ200におけるピーク波形との有無の差分を、周波数ごとに抽出する。図2では、制限周波数データ100におけるピーク波形は、一例として、ピーク波形101とピーク波形102とである。図2では、非制限周波数データ200におけるピーク波形は、一例として、ピーク波形101とピーク波形102とピーク波形103とである。
Next, measurement processing will be described.
FIG. 2 is a diagram showing an example of measurement processing in the first embodiment. The measurement unit 14a extracts the difference between the peak waveform in the restricted
計測部14aは、非制限周波数データ200において差分として抽出されたピーク波形103の周波数(卓越周波数)を計測する。図2では、ピーク波形103の周波数は、一例として30Hzである。計測部14aは、ナイキスト周波数(50Hz)を基準として、差分として抽出されたピーク波形103の周波数を折り返す。すなわち、計測部14aは、ナイキスト周波数を鏡として、ピーク波形103を反転する。この結果、折り返された周波数は、70Hz(=50+(50-30)Hz)である。計測部14aは、制限周波数データ100と、折り返された周波数(70Hz)のピーク波形104のデータとを、出力部16に出力する。
The measurement unit 14 a measures the frequency (dominant frequency) of the
次に、計測装置1aの動作例を説明する。
図3は、第1実施形態における、計測装置1aの動作例を示すフローチャートである。切替部10は、波形信号(アナログ信号)をセンサ等(不図示)から取得する(ステップS101)。計測部14aは、アナログ・デジタル変換部12のサンプリング周波数(例えば、100Hz)に基づいて、ナイキスト周波数を導出する(ステップS102)。
Next, an operation example of the measuring
FIG. 3 is a flow chart showing an operation example of the measuring
第1期間において、ローパスフィルタ11は、入力された波形信号に対して、ローパスフィルタ処理を実行する(ステップS103)。アナログ・デジタル変換部12は、サンプリング周波数に基づいて、帯域が制限された波形信号(アナログ信号)に対して、サンプリング処理を実行する。これによって、アナログ・デジタル変換部12は、帯域が制限された波形信号(アナログ信号)を、帯域が制限された波形信号(デジタル信号)に変換する(ステップS104-1)。高速フーリエ変換部13は、ナイキスト周波数以下の帯域について、帯域が制限された波形信号(デジタル信号)に対して高速フーリエ変換を実行する。これによって、高速フーリエ変換部13は、制限周波数データを生成する(ステップS105-1)。
In the first period, the low-
第1期間とは異なる第2期間において、アナログ・デジタル変換部12は、サンプリング周波数に基づいて、帯域が制限されていない波形信号(アナログ信号)に対して、サンプリング処理を実行する。これによって、アナログ・デジタル変換部12は、帯域が制限されていない波形信号(アナログ信号)を、帯域が制限されていない波形信号(デジタル信号)に変換する(ステップS104-2)。高速フーリエ変換部13は、ナイキスト周波数以下の帯域について、帯域が制限されていない波形信号(デジタル信号)に対して高速フーリエ変換を実行する。これによって、高速フーリエ変換部13は、非制限周波数データを生成する(ステップS105-2)。
In a second period different from the first period, the analog-to-
計測部14aは、制限周波数データにおけるピーク波形と非制限周波数データにおけるピーク波形との有無の差分を、周波数ごとに抽出する(ステップS106)。計測部14aは、ナイキスト周波数を基準として、差分のピーク波形の周波数を折り返す(ステップS107)。計測部14aは、サンプリング周波数以下の帯域について、周波数分布を生成する(ステップS108)。出力部16は、サンプリング周波数以下の帯域について、周波数分布を出力する。すなわち、出力部16は、ナイキスト周波数の2逓倍以下の帯域について、周波数分布を出力する(ステップS109)。
The measurement unit 14a extracts, for each frequency, the difference between the peak waveform in the restricted frequency data and the peak waveform in the non-restricted frequency data (step S106). The measurement unit 14a repeats the frequency of the differential peak waveform with reference to the Nyquist frequency (step S107). The measurement unit 14a generates a frequency distribution for bands below the sampling frequency (step S108). The
以上のように、ローパスフィルタ11は、入力されたアナログ信号(波形信号)の帯域を、ナイキスト周波数以下の帯域となるように制限する。アナログ・デジタル変換部12(サンプリング実行部)は、ナイキスト周波数以下の帯域となるように帯域が制限されたアナログ信号に基づく第1デジタル信号を、サンプリング周波数に基づいて生成する。アナログ・デジタル変換部12は、帯域が制限されていないアナログ信号に基づく第2デジタル信号を、サンプリング周波数に基づいて生成する。高速フーリエ変換部13は、第1デジタル信号の周波数分布データ(制限周波数データ100)と、第2デジタル信号の周波数分布データ(非制限周波数データ200)とを、ナイキスト周波数以下の帯域について生成する。
As described above, the low-
計測部14a(抽出部)は、第1デジタル信号の周波数分布データ(制限周波数データ100)におけるピーク波形の有無と、第2デジタル信号の周波数分布データ(非制限周波数データ200)におけるピーク波形の有無との差分を、周波数ごとに抽出する。計測部14aは、ナイキスト周波数を基準として、差分に基づいてピーク波形103の周波数を折り返す。ここで、計測部14aは、制限周波数データ100との差分として非制限周波数データ200から抽出されたピーク波形103の周波数を、ナイキスト周波数を基準として折り返す。
The measurement unit 14a (extraction unit) determines whether or not there is a peak waveform in the frequency distribution data (limit frequency data 100) of the first digital signal and whether or not there is a peak waveform in the frequency distribution data (unlimited frequency data 200) of the second digital signal. , is extracted for each frequency. The measurement unit 14a repeats the frequency of the
これによって、波形信号においてナイキスト周波数を超える卓越周波数を計測することが可能である。すなわち、エイリアシング(折り返し雑音)の現象を利用して、周波数分布をサンプリング周波数まで計測することが可能である。ここで、計測対象(軸受等)の周波数(損傷周波数)が予め特定されている必要がない。すなわち、取り扱われる周波数領域が予め特定されている必要がない。サンプリング周波数を不要に高くする必要がないので、計測装置のコストの増加を抑制することが可能である。サンプリング周波数を不要に高くする必要がないので、計測対象の周波数(損傷周波数)を計測するためのデータ容量の増加を抑制することが可能である。 This makes it possible to measure dominant frequencies exceeding the Nyquist frequency in the waveform signal. That is, the phenomenon of aliasing (aliasing noise) can be used to measure the frequency distribution up to the sampling frequency. Here, the frequency (damage frequency) of the object to be measured (bearing, etc.) need not be specified in advance. That is, there is no need to specify the frequency range to be handled in advance. Since there is no need to increase the sampling frequency unnecessarily, it is possible to suppress an increase in the cost of the measuring device. Since there is no need to increase the sampling frequency unnecessarily, it is possible to suppress an increase in data capacity for measuring the frequency of the measurement object (damage frequency).
(第2実施形態)
第2実施形態では、計測装置がローパスフィルタを備えない点と、計測装置がハイパスフィルタを備える点とが、第1実施形態との差分である。第2実施形態では、第1実施形態との差分を中心に説明する。
(Second embodiment)
The difference between the second embodiment and the first embodiment is that the measuring device does not have a low-pass filter and the measuring device has a high-pass filter. 2nd Embodiment demonstrates centering around the difference with 1st Embodiment.
図4は、第2実施形態における、計測装置1bの構成例を示す図である。計測装置1bは、周波数を計測する装置である。計測装置1bは、切替部10と、アナログ・デジタル変換部12と、高速フーリエ変換部13と、計測部14bと、記憶装置15と、出力部16と、ハイパスフィルタ17とを備える。なお、ハイパスフィルタ17は、バンドパスフィルタでもよい。
FIG. 4 is a diagram showing a configuration example of the measuring
ハイパスフィルタ17は、入力された波形信号の帯域を制限する機能部である。ハイパスフィルタ17は、波形信号が切替部10から入力された場合、ナイキスト周波数以上の帯域の波形信号を、アナログ・デジタル変換部12に出力する。ハイパスフィルタ17のカットオフ周波数は、例えば、ナイキスト周波数に等しい周波数である。
The high-
アナログ・デジタル変換部12には、例えば所定周期で、波形信号が切替部10及びハイパスフィルタ17から交互に入力される。すなわち、第1期間において、波形信号が、アナログ・デジタル変換部12にハイパスフィルタ17から入力される。第1期間とは異なる第2期間において、波形信号が、アナログ・デジタル変換部12に切替部10から入力される。
Waveform signals are alternately input to the analog/
アナログ・デジタル変換部12にハイパスフィルタ17から波形信号が入力されている場合、高速フーリエ変換部13は、ハイパスフィルタ17によって帯域が制限された波形信号に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13は、制限周波数データを計測部14bに出力する。
When a waveform signal is input from the high-
アナログ・デジタル変換部12に切替部10から波形信号が入力されている場合、高速フーリエ変換部13は、ハイパスフィルタ17によって帯域が制限されていない波形信号に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13は、非制限周波数データを計測部14bに出力する。
When a waveform signal is input from the switching
次に、計測処理について説明する。
図5は、第2実施形態における、計測処理の例を示す図である。計測部14bは、制限周波数データ400におけるピーク波形と非制限周波数データ200におけるピーク波形との有無の差分を、周波数ごとに抽出する。図5では、制限周波数データ400におけるピーク波形は、一例として、ピーク波形103である。図5では、非制限周波数データ200におけるピーク波形は、一例として、ピーク波形101とピーク波形102とピーク波形103とである。
Next, measurement processing will be described.
FIG. 5 is a diagram showing an example of measurement processing in the second embodiment. The measuring unit 14b extracts the difference between the peak waveform in the restricted
計測部14bは、非制限周波数データ200において差分(相違点)として抽出されないピーク波形103の周波数(卓越周波数)を計測する。図5では、ピーク波形103の周波数は、一例として30Hzである。計測部14bは、ナイキスト周波数(50Hz)を基準として、差分として抽出されないピーク波形103(共通点として抽出されたピーク波形103)の周波数を折り返す。すなわち、計測部14aは、ナイキスト周波数を鏡として、ピーク波形103を反転する。この結果、折り返された周波数は、70Hz(=50+(50-30)Hz)である。計測部14bは、非制限周波数データ200と、折り返された周波数(70Hz)のピーク波形104のデータとを、出力部16に出力する。
The measurement unit 14b measures the frequency (predominant frequency) of the
次に、計測装置1bの動作例を説明する。
図6は、第2実施形態における、計測装置1bの動作例を示すフローチャートである。ステップS201からステップS202までの各動作は、図3に示されたステップS101からステップS102までの各動作と同様である。
Next, an operation example of the measuring
FIG. 6 is a flow chart showing an operation example of the measuring
第1期間において、ハイパスフィルタ17は、入力された波形信号に対して、ハイパスフィルタ処理を実行する(ステップS203)。ステップS204-1からステップS205-1までの各動作は、図3に示されたステップS104-1からステップS105-1までの各動作と同様である。また、ステップS204-2からステップS205-2までの各動作は、図3に示されたステップS104-2からステップS105-2までの各動作と同様である。
In the first period, the high-
ステップS206の動作は、図3に示されたステップS106の動作と同様である。計測部14bは、ナイキスト周波数を基準として、差分として抽出されないピーク波形103の周波数を折り返す(ステップS207)。ステップS208からステップS209までの各動作は、図3に示されたステップS108からステップS109までの各動作と同様である。
The operation of step S206 is the same as the operation of step S106 shown in FIG. The measurement unit 14b repeats the frequency of the
以上のように、ハイパスフィルタ17は、入力されたアナログ信号(波形信号)の帯域を、ナイキスト周波数以上の帯域となるように制限する。アナログ・デジタル変換部12(サンプリング実行部)は、ナイキスト周波数以上の帯域となるように帯域が制限されたアナログ信号に基づく第1デジタル信号を生成する。アナログ・デジタル変換部12は、帯域が制限されていないアナログ信号に基づく第2デジタル信号を生成する。計測部14b(抽出部)は、第2デジタル信号の周波数分布データ(非制限周波数データ200)との共通点として第1デジタル信号の周波数分布データ(制限周波数データ400)から抽出されたピーク波形103の周波数を、ナイキスト周波数を基準として折り返す。
As described above, the high-
これによって、波形信号においてナイキスト周波数を超える卓越周波数を計測することが可能である。計測対象(軸受等)の周波数(損傷周波数)が予め特定されている必要がない。すなわち、取り扱われる周波数領域が予め特定されている必要がない。サンプリング周波数を不要に高くする必要がないので、計測装置のコストの増加を抑制することが可能である。サンプリング周波数を不要に高くする必要がないので、計測対象の周波数(損傷周波数)を計測するためのデータ容量の増加を抑制することが可能である。 This makes it possible to measure dominant frequencies exceeding the Nyquist frequency in the waveform signal. The frequency (damage frequency) of the object to be measured (bearing, etc.) need not be specified in advance. That is, there is no need to specify the frequency range to be handled in advance. Since there is no need to increase the sampling frequency unnecessarily, it is possible to suppress an increase in the cost of the measuring device. Since there is no need to increase the sampling frequency unnecessarily, it is possible to suppress an increase in data capacity for measuring the frequency of the measurement object (damage frequency).
(第3実施形態)
第3実施形態では、ローパスフィルタとハイパスフィルタとを計測装置が並列に備える点が、第1実施形態との差分である。第3実施形態では、第1実施形態との差分を中心に説明する。
(Third embodiment)
The difference between the third embodiment and the first embodiment is that the measuring device includes a low-pass filter and a high-pass filter in parallel. 3rd Embodiment demonstrates centering around the difference with 1st Embodiment.
図7は、第3実施形態における、計測装置1cの構成例を示す図である。計測装置1cは、周波数を計測する装置である。計測装置1cは、切替部10と、ローパスフィルタ11と、アナログ・デジタル変換部12と、高速フーリエ変換部13と、計測部14aと、記憶装置15と、出力部16と、ハイパスフィルタ17とを備える。なお、ローパスフィルタ11とハイパスフィルタ17とのうちの少なくとも一方は、バンドパスフィルタでもよい。
FIG. 7 is a diagram showing a configuration example of the measuring device 1c in the third embodiment. The measuring device 1c is a device for measuring frequency. The measurement device 1c includes a
切替部10は、センサ又は測定装置から、波形信号(アナログ信号)を取得する。切替部10は、ローパスフィルタ11及びハイパスフィルタ17に、例えば所定周期で交互に波形信号を出力する。すなわち、切替部10は、波形信号の出力先を、例えば所定周期で交互に切り替える。
The switching
ローパスフィルタ11は、波形信号が切替部10から入力された場合、ナイキスト周波数以下の帯域の波形信号を、アナログ・デジタル変換部12に出力する。ハイパスフィルタ17は、波形信号が切替部10から入力された場合、ナイキスト周波数以上の帯域の波形信号を、アナログ・デジタル変換部12に出力する。
When a waveform signal is input from the switching
アナログ・デジタル変換部12には、例えば所定周期で、波形信号がローパスフィルタ11及びハイパスフィルタ17から交互に入力される。すなわち、第1期間において、波形信号が、アナログ・デジタル変換部12にローパスフィルタ11から入力される。第1期間とは異なる第2期間において、波形信号が、アナログ・デジタル変換部12にハイパスフィルタ17から入力される。
Waveform signals are alternately input from the low-
アナログ・デジタル変換部12にローパスフィルタ11から波形信号が入力されている場合、高速フーリエ変換部13は、ローパスフィルタ11によって帯域が制限された波形信号に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13は、第1の制限周波数データを計測部14cに出力する。
When a waveform signal is input from the low-
アナログ・デジタル変換部12にハイパスフィルタ17から波形信号が入力されている場合、高速フーリエ変換部13は、ハイパスフィルタ17によって帯域が制限された波形信号に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13は、第2の制限周波数データを計測部14cに出力する。
When a waveform signal is input from the high-
計測部14cは、アナログ・デジタル変換部12にローパスフィルタ11から波形信号が入力されている場合、第1の制限周波数データを高速フーリエ変換部13から取得する。計測部14cは、第1の制限周波数データを記憶装置15に記録する。計測部14cは、アナログ・デジタル変換部12にハイパスフィルタ17から波形信号が入力されている場合、第2の制限周波数データを高速フーリエ変換部13から取得する。計測部14cは、第2の制限周波数データを記憶装置15に記録する。
The
計測部14cは、第1の制限周波数データにおけるピーク波形と第2の制限周波数データにおけるピーク波形との有無の差分を、周波数ごとに抽出する。計測部14cは、第2の制限周波数データにおいて差分として抽出されたピーク波形の周波数(卓越周波数)を計測する。計測部14cは、ナイキスト周波数を基準として、差分として抽出されたピーク波形の周波数を折り返す。計測部14cは、制限周波数データと、折り返された周波数のピーク波形のデータとを、出力部16に出力する。
The measuring
次に、計測処理について説明する。
図8は、第3実施形態における、計測処理の例を示す図である。計測部14cは、制限周波数データ100におけるピーク波形と制限周波数データ400におけるピーク波形との有無の差分を、周波数ごとに抽出する。図8では、制限周波数データ100におけるピーク波形は、一例として、ピーク波形101とピーク波形102とである。図8では、制限周波数データ400におけるピーク波形は、一例として、ピーク波形103である。
Next, measurement processing will be described.
FIG. 8 is a diagram showing an example of measurement processing in the third embodiment. The measuring
計測部14cは、制限周波数データ400において差分として抽出されたピーク波形103の周波数(卓越周波数)を計測する。図8では、ピーク波形103の周波数は、一例として30Hzである。計測部14cは、ナイキスト周波数(50Hz)を基準として、差分として抽出されたピーク波形103の周波数を折り返す。すなわち、計測部14cは、ナイキスト周波数を鏡として、ピーク波形103を反転する。この結果、折り返された周波数は、70Hz(=50+(50-30)Hz)である。計測部14cは、制限周波数データ100と、折り返された周波数(70Hz)のピーク波形104のデータとを、出力部16に出力する。
The
次に、計測装置1cの動作例を説明する。
図9は、第3実施形態における、計測装置1cの動作例を示すフローチャートである。第1期間において、ステップS301からステップS305-1までの各動作は、図3に示されたステップS101からステップS105-1までの各動作と同様である。第2期間において、ステップS303-2からステップS305-2までの各動作は、図6に示されたステップS203からステップS205-1までの各動作と同様である。
Next, an operation example of the measuring device 1c will be described.
FIG. 9 is a flow chart showing an operation example of the measuring device 1c in the third embodiment. In the first period, each operation from step S301 to step S305-1 is the same as each operation from step S101 to step S105-1 shown in FIG. In the second period, each operation from step S303-2 to step S305-2 is the same as each operation from step S203 to step S205-1 shown in FIG.
計測部14cは、第1の制限周波数データにおけるピーク波形と第2の制限周波数データにおけるピーク波形との有無の差分を、周波数ごとに抽出する(ステップS306)。計測部14cは、ナイキスト周波数を基準として、第2の制限周波数データにおける差分のピーク波形の周波数を折り返す(ステップS307)。ステップS308からステップS309までの各動作は、図3に示されたステップS308からステップS309までの各動作と同様である。
The
以上のように、ローパスフィルタ11は、入力されたアナログ信号(波形信号)の帯域を、ナイキスト周波数以下の帯域となるように制限する。ハイパスフィルタ17は、入力されたアナログ信号(波形信号)の帯域を、ナイキスト周波数以上の帯域となるように制限する。アナログ・デジタル変換部12(サンプリング実行部)は、ナイキスト周波数以下の帯域となるように帯域が制限されたアナログ信号に基づく第1デジタル信号を生成する。アナログ・デジタル変換部12は、ナイキスト周波数以上の帯域となるように帯域が制限されたアナログ信号に基づく第2デジタル信号を生成する。計測部14c(抽出部)は、第1デジタル信号の周波数分布データ(制限周波数データ100)との差分として第2デジタル信号の周波数分布データ(制限周波数データ400)から抽出されたピーク波形103の周波数を、ナイキスト周波数を基準として折り返す。
As described above, the low-
第3実施形態ではフィルタの数が多いので、騒音の周波数が広帯域に分布する環境でも、波形信号においてナイキスト周波数を超える卓越周波数を計測することが可能である。 Since the number of filters in the third embodiment is large, even in an environment where noise frequencies are distributed over a wide band, it is possible to measure dominant frequencies exceeding the Nyquist frequency in the waveform signal.
(第4実施形態)
第4実施形態では、ローパスフィルタの系統とハイパスフィルタの系統との両方に波形信号が並列に分配される点が、第3実施形態との差分である。第4実施形態では、第3実施形態との差分を中心に説明する。
(Fourth embodiment)
The fourth embodiment differs from the third embodiment in that waveform signals are distributed in parallel to both the low-pass filter system and the high-pass filter system. 4th Embodiment demonstrates centering around the difference with 3rd Embodiment.
図10は、第4実施形態における、計測装置1dの構成例を示す図である。計測装置1dは、周波数を計測する装置である。計測装置1dは、ローパスフィルタ11と、アナログ・デジタル変換部12-1と、アナログ・デジタル変換部12-2と、高速フーリエ変換部13-1と、高速フーリエ変換部13-2と、計測部14aと、記憶装置15と、出力部16と、ハイパスフィルタ17と、分配部18とを備える。なお、ローパスフィルタ11とハイパスフィルタ17とのうちの少なくとも一方は、バンドパスフィルタでもよい。
FIG. 10 is a diagram showing a configuration example of a
分配部18は、入力された波形信号を、ローパスフィルタ11の系統とハイパスフィルタ17の系統との両方に(並列に)分配する。ローパスフィルタ11は、ナイキスト周波数以下の帯域の波形信号を、アナログ・デジタル変換部12-1に出力する。アナログ・デジタル変換部12-1には、波形信号がローパスフィルタ11から入力される。高速フーリエ変換部13-1は、ローパスフィルタ11によって帯域が制限された波形信号に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13-1は、制限周波数データ100を計測部14dに出力する。
The
ハイパスフィルタ17は、ナイキスト周波数以上の帯域の波形信号を、アナログ・デジタル変換部12-2に出力する。アナログ・デジタル変換部12-2には、波形信号がハイパスフィルタ17から入力される。高速フーリエ変換部13-2は、ハイパスフィルタ17によって帯域が制限された波形信号に対して、高速フーリエ変換処理を実行する。高速フーリエ変換部13は、制限周波数データ400を計測部14dに出力する。
The high-
計測部14dは、制限周波数データ100におけるピーク波形と制限周波数データ400におけるピーク波形との有無の差分を、周波数ごとに抽出する。計測部14dは、制限周波数データ400において差分として抽出されたピーク波形103の周波数(卓越周波数)を計測する。計測部14dは、ナイキスト周波数を基準として、差分として抽出されたピーク波形103の周波数を折り返す。計測部14dは、制限周波数データ100と、折り返された周波数のピーク波形103のデータとを、出力部16に出力する。
The measuring
以上のように、分配部18は、入力された波形信号を、ローパスフィルタ11の系統とハイパスフィルタ17の系統との両方に(並列に)分配する。これによって、波形信号の入力系統が切り替えられる必要がないので、短時間で変動する波形信号(定常的でない波形信号)においてナイキスト周波数を超える卓越周波数を計測することが可能である。
As described above, the
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。 Although the embodiment of the present invention has been described in detail with reference to the drawings, the specific configuration is not limited to this embodiment, and design and the like are included within the scope of the gist of the present invention.
例えば、計測装置は、卓越周波数に基づいて計測対象の異常を診断してもよい。 For example, the measuring device may diagnose abnormalities in the measurement target based on the dominant frequency.
1a,1b,1c,1d…計測装置、10…切替部、11…ローパスフィルタ、12…アナログ・デジタル変換部、13…高速フーリエ変換部、14a,14b,14c,14d…計測部、15…記憶装置、16…出力部、17…ハイパスフィルタ、18…分配部、100…制限周波数データ、101…ピーク波形、102…ピーク波形、103…ピーク波形、104…ピーク波形、200…非制限周波数データ、300…計測周波数データ、400…制限周波数データ
1a, 1b, 1c, 1d... Measuring
Claims (6)
帯域が制限された前記アナログ信号に基づく第1デジタル信号と帯域が制限された又は制限されていない前記アナログ信号に基づく第2デジタル信号とを、サンプリング周波数に基づいて生成するサンプリング実行部と、
前記第1デジタル信号の周波数分布データと前記第2デジタル信号の周波数分布データとを、前記サンプリング周波数に基づくナイキスト周波数以下の帯域について生成する高速フーリエ変換部と、
前記第1デジタル信号の周波数分布データにおけるピーク波形の有無と前記第2デジタル信号の周波数分布データにおけるピーク波形の有無との差分又は共通点を周波数ごとに抽出し、前記ナイキスト周波数を基準として、前記差分又は前記共通点に基づいてピーク波形の周波数を折り返す抽出部と
を備える計測装置。 a filter that limits the band of the input analog signal;
a sampling execution unit that generates a first digital signal based on the band-limited analog signal and a second digital signal based on the band-limited or unband-limited analog signal based on a sampling frequency;
a fast Fourier transform unit that generates the frequency distribution data of the first digital signal and the frequency distribution data of the second digital signal for a band equal to or lower than the Nyquist frequency based on the sampling frequency;
A difference or a common point between the presence or absence of a peak waveform in the frequency distribution data of the first digital signal and the presence or absence of a peak waveform in the frequency distribution data of the second digital signal is extracted for each frequency, and the Nyquist frequency is used as a reference, A measurement device comprising: an extractor that folds back the frequency of the peak waveform based on the difference or the common point.
前記抽出部は、前記第1デジタル信号の周波数分布データとの差分として前記第2デジタル信号の周波数分布データから抽出されたピーク波形の周波数を、前記ナイキスト周波数を基準として折り返す、
請求項1に記載の計測装置。 The sampling execution unit generates a first digital signal based on the analog signal whose band is limited to a band equal to or lower than the Nyquist frequency, and a second digital signal based on the analog signal whose band is not limited. death,
The extracting unit folds the frequency of the peak waveform extracted from the frequency distribution data of the second digital signal as a difference from the frequency distribution data of the first digital signal with reference to the Nyquist frequency,
The measuring device according to claim 1.
前記抽出部は、前記第2デジタル信号の周波数分布データとの共通点として前記第1デジタル信号の周波数分布データから抽出されたピーク波形の周波数を、前記ナイキスト周波数を基準として折り返す、
請求項1に記載の計測装置。 The sampling execution unit generates a first digital signal based on the analog signal whose band is limited to a band equal to or higher than the Nyquist frequency, and a second digital signal based on the analog signal whose band is not limited. death,
The extraction unit repeats the frequency of the peak waveform extracted from the frequency distribution data of the first digital signal as a common point with the frequency distribution data of the second digital signal, with the Nyquist frequency as a reference.
The measuring device according to claim 1.
前記抽出部は、前記第1デジタル信号の周波数分布データとの差分として前記第2デジタル信号の周波数分布データから抽出されたピーク波形の周波数を、前記ナイキスト周波数を基準として折り返す、
請求項1に記載の計測装置。 The sampling execution unit generates a first digital signal based on the analog signal whose band is limited to a band equal to or lower than the Nyquist frequency, and the analog signal whose band is limited to a band equal to or higher than the Nyquist frequency. and a second digital signal based on
The extracting unit folds the frequency of the peak waveform extracted from the frequency distribution data of the second digital signal as a difference from the frequency distribution data of the first digital signal with reference to the Nyquist frequency,
The measuring device according to claim 1.
入力されたアナログ信号の帯域を制限するフィルタステップと、
帯域が制限された前記アナログ信号に基づく第1デジタル信号と帯域が制限された又は制限されていない前記アナログ信号に基づく第2デジタル信号とを、サンプリング周波数に基づいて生成するサンプリング実行ステップと、
前記第1デジタル信号の周波数分布データと前記第2デジタル信号の周波数分布データとを、前記サンプリング周波数に基づくナイキスト周波数以下の帯域について生成する高速フーリエ変換ステップと、
前記第1デジタル信号の周波数分布データにおけるピーク波形の有無と前記第2デジタル信号の周波数分布データにおけるピーク波形の有無との差分又は共通点を周波数ごとに抽出し、前記ナイキスト周波数を基準として、前記差分又は前記共通点に基づいてピーク波形の周波数を折り返す抽出ステップと
を含む計測方法。 A measurement method performed by a measurement device,
a filter step that limits the band of the input analog signal;
a sampling performing step of generating a first digital signal based on the band-limited analog signal and a second digital signal based on the band-limited or unband-limited analog signal based on a sampling frequency;
a fast Fourier transform step of generating the frequency distribution data of the first digital signal and the frequency distribution data of the second digital signal for a band below the Nyquist frequency based on the sampling frequency;
A difference or a common point between the presence or absence of a peak waveform in the frequency distribution data of the first digital signal and the presence or absence of a peak waveform in the frequency distribution data of the second digital signal is extracted for each frequency, and the Nyquist frequency is used as a reference, and an extracting step of folding the frequency of the peak waveform based on the difference or the common point.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021075770A JP2022169984A (en) | 2021-04-28 | 2021-04-28 | Measurement device, measurement method and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021075770A JP2022169984A (en) | 2021-04-28 | 2021-04-28 | Measurement device, measurement method and program |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022169984A true JP2022169984A (en) | 2022-11-10 |
Family
ID=83944881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021075770A Pending JP2022169984A (en) | 2021-04-28 | 2021-04-28 | Measurement device, measurement method and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2022169984A (en) |
-
2021
- 2021-04-28 JP JP2021075770A patent/JP2022169984A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10634554B2 (en) | Integrated vibration measurement and analysis system | |
RU2383894C2 (en) | Method and device for detection and measurement of frequency disturbances of rotor rotation speed | |
JP2015505058A (en) | A system for acquiring vibration signals of rotary motors | |
JP2011232334A (en) | Axial vibration monitoring for detecting shaft misalignment in turbomachinery train | |
JP2014077791A (en) | Method of indicating correlation between multiple signals, and test measurement device | |
JP2022169984A (en) | Measurement device, measurement method and program | |
Samoylov et al. | Estimation to efficiency of the using of anti-alias filter in the A/D interface of instrumentation and control systems | |
JP2009177835A (en) | Device for signal-processing time discrete values | |
Katayama et al. | Low distortion sine wave generator with simple harmonics cancellation circuit and filter for analog device testing | |
JP2013160534A (en) | Inclination measuring apparatus and inclination measuring system with the same | |
JP2001074547A (en) | Characteristic display method, characteristic display system, and record medium storing characteristic display program | |
JP5867209B2 (en) | Sound removal apparatus, sound inspection apparatus, sound removal method, and sound removal program | |
Szary et al. | Signal conditioning perspectives on pyroshock measurement systems | |
Giannoccaro et al. | A digital analysis of the experimental accelerometers data used for buildings dynamical identification | |
JP7196504B2 (en) | Acoustic property measurement system, acoustic property measurement method, and program | |
JP5191184B2 (en) | Downsampling apparatus, downsampling method, and program thereof | |
JP7351193B2 (en) | Acoustic property measurement system, acoustic property measurement method, and acoustic property measurement program | |
JP2006081168A (en) | Audio signal processing apparatus | |
Yanagida et al. | Low-distortion one-tone and two-tone signal generation using AWG over full Nyquist region | |
JP7056739B2 (en) | Wave source direction estimator, wave source direction estimation method, and program | |
JP6505546B2 (en) | Sound effect giving device | |
WO2024057364A1 (en) | Vibration detection device, abnormality detection assistance system, and vibration detection method | |
JP2020022055A (en) | Sound output device and sound output system | |
JP4950801B2 (en) | Upsampling apparatus, upsampling method, and program thereof | |
JP4229231B2 (en) | Analog-to-digital conversion method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20231206 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240515 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240521 |