JP4865582B2 - 回転計と回転数の計測方法 - Google Patents
回転計と回転数の計測方法 Download PDFInfo
- Publication number
- JP4865582B2 JP4865582B2 JP2007030090A JP2007030090A JP4865582B2 JP 4865582 B2 JP4865582 B2 JP 4865582B2 JP 2007030090 A JP2007030090 A JP 2007030090A JP 2007030090 A JP2007030090 A JP 2007030090A JP 4865582 B2 JP4865582 B2 JP 4865582B2
- Authority
- JP
- Japan
- Prior art keywords
- spectrum
- frequency
- search range
- specific
- search
- 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.)
- Expired - Fee Related
Links
- 238000000691 measurement method Methods 0.000 title claims description 3
- 238000001228 spectrum Methods 0.000 claims description 261
- 238000004364 calculation method Methods 0.000 claims description 82
- 238000005259 measurement Methods 0.000 claims description 44
- 230000003595 spectral effect Effects 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 7
- 238000000034 method Methods 0.000 claims description 5
- 238000004880 explosion Methods 0.000 description 13
- 238000012545 processing Methods 0.000 description 13
- 238000002360 preparation method Methods 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 239000000446 fuel Substances 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 239000000567 combustion gas Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000001845 vibrational spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/0097—Electrical control of supply of combustible mixture or its constituents using means for generating speed signals
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/24—Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means
- F02D41/26—Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means using computer, e.g. microprocessor
- F02D41/28—Interface circuits
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/42—Devices characterised by the use of electric or magnetic means
- G01P3/44—Devices characterised by the use of electric or magnetic means for measuring angular speed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/42—Devices characterised by the use of electric or magnetic means
- G01P3/44—Devices characterised by the use of electric or magnetic means for measuring angular speed
- G01P3/48—Devices characterised by the use of electric or magnetic means for measuring angular speed by measuring frequency of generated current or voltage
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/02—Circuit arrangements for generating control signals
- F02D41/14—Introducing closed-loop corrections
- F02D41/1401—Introducing closed-loop corrections characterised by the control or regulation method
- F02D2041/1413—Controller structures or design
- F02D2041/1432—Controller structures or design the system including a filter, e.g. a low pass or high pass filter
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F02—COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
- F02D—CONTROLLING COMBUSTION ENGINES
- F02D41/00—Electrical control of supply of combustible mixture or its constituents
- F02D41/24—Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means
- F02D41/26—Electrical control of supply of combustible mixture or its constituents characterised by the use of digital means using computer, e.g. microprocessor
- F02D41/28—Interface circuits
- F02D2041/286—Interface circuits comprising means for signal processing
- F02D2041/288—Interface circuits comprising means for signal processing for performing a transformation into the frequency domain, e.g. Fourier transformation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mechanical Engineering (AREA)
- Computational Mathematics (AREA)
- Combustion & Propulsion (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Optimization (AREA)
- Discrete Mathematics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Description
本発明は、外部から入力した信号を処理する信号処理装置及び信号探索範囲の追従方法に関するものである。
信号処理装置は、その用途に応じてさまざまな機能が設けられている（例えば、特許文献１参照）。特許文献１の発明は、回転モータの回転子の回転数を測定する信号処理装置（回転計）であり、回転子の回転に起因する物理現象として漏洩磁束を検出し、高速フーリエ変換（ＦＦＴ）により漏洩磁束のスペクトルデータを求めて、そのスペクトルデータの最大値を示す周波数に基づいて回転モータの回転数を求めている。
しかし、このような回転計でＦＦＴ演算されたスペクトルデータには、回転モータの回転数となる主成分（１次成分）のスペクトルだけでなく、０．５次、２次、３次など多数のスペクトルが含まれる。一般的には、スペクトルデータ内の最大値を示すスペクトルが、回転の主成分となって演算されるが、回転モータに対する外乱などの影響や、測定に用いた回転に起因する物理現象として測定された信号の種類によっては、上述の０．５次、２次、３次などのスペクトルが、主成分のスペクトルよりも大きくなってしまい、異なるスペクトルを主成分のスペクトルとして誤って読み取ってしまう問題があった。
そのため、主成分のスペクトルを探索する探索範囲を、上述の０．５次、２次、３次などの周波数成分が含まれないように設定するなどの対策が実施されていた。例えば、探索範囲の下限を主成分の０．６倍とし、上限を主成分の１．８倍として、周波数軸上で主成分（１次成分）のスペクトルと隣り合う０．５次及び２次のスペクトルを除外するようにしている。しかし、このような探索範囲の設定では、回転モータの回転数が変動したときに、その探索範囲内に主成分以外の０．５次成分や２次成分などのスペクトルが含まれてしまい、そのスペクトルを回転モータの回転数に起因する主成分のスペクトルとして誤って読み取ってしまうことがあった。
具体的には、回転モータの主成分が１００Ｈｚであり、主成分の探索範囲が６０〜１８０Ｈｚに設定されているとき、０．５次成分である５０Ｈｚや、２次成分である２００Ｈｚは探索範囲から除外されている。しかし、回転モータの主成分が１００Ｈｚから１５０Ｈｚに変化したとき、探索範囲（６０〜１８０Ｈｚ）内に、変動した主成分（１５０Ｈｚ）の０．５次成分である７５Ｈｚが読み取られ、主成分（１５０Ｈｚ）のスペクトルが０．５次成分（７５Ｈｚ）のそれよりも小さい場合、０．５次成分のスペクトルを主成分として読み取ってしまうことがあった。また、主成分が探索範囲（６０〜１８０Ｈｚ）から外れた場合（例えば、２００Ｈｚ）においても、０．５次などの主成分以外のスペクトルを読み取ってしまう問題があった。
特開平２−２１２６６号公報
本発明の課題は、周波数の変動した特定のスペクトルを正確に読み取ることができる信号処理装置及び信号探索範囲の追従方法を提供することである。
前記課題を解決するために、請求項１の発明は、測定対象から計測された回転運動に起因する時系列の振動の検出結果を入力した信号を、高速フーリエ変換してスペクトルデータを演算するＦＦＴ演算部と、前記測定対象の主成分を特定のスペクトルとして表示して選択可能なスペクトル選択部と、前記スペクトル選択部で選択した前記特定のスペクトルを、前記ＦＦＴ演算部により演算された前記スペクトルデータのうちの所定のスペクトルの探索範囲から逐次探索するスペクトル探索部と、前記スペクトル探索部で探索される前記所定のスペクトルの探索範囲を、前記特定のスペクトルの周波数の２倍よりも小さい周波数を上限とし、前記特定のスペクトルの周波数の０．５倍よりも大きい周波数を下限として、演算する探索範囲演算部と、前記スペクトル探索部によって探索された主成分スペクトルを示す周波数に基づいて、前記測定対象の回転数を演算する回転数演算部と、を含む演算制御部を備え、前記演算制御部は、試運転時には、前記ＦＦＴ演算部、前記スペクトル選択部及び前記探索範囲演算部を動作させ、測定開始時には、前記ＦＦＴ演算部、前記スペクトル探索部、前記探索範囲演算部及び前記回転数演算部を動作させると共に、前記探索範囲演算部を、前記特定のスペクトルの周波数が変動したときに、前記探索範囲を、前記特定のスペクトルの変動した周波数に基づいて演算し、変動した前記特定のスペクトルに追従するように動作させること、を特徴とする回転計である。
請求項２の発明は、請求項１に記載の回転計において、前記スペクトル探索部は、前記探索範囲内における前記スペクトルデータの最大値を示す周波数を前記特定のスペクトルとして探索すること、を特徴とする回転計である。
請求項３の発明は、請求項１に記載の回転計において、前記特定のスペクトルが、前記測定対象の回転周波数のｘ倍の周波数であるときに、前記探索範囲演算部は、前記特定のスペクトルの周波数をｘで除算した周波数に基づいて探索範囲の幅を演算し、演算した前記探索範囲の幅を前記特定のスペクトルの周波数の探索範囲として設定すること、を特徴とする回転計である。
請求項２の発明は、請求項１に記載の回転計において、前記スペクトル探索部は、前記探索範囲内における前記スペクトルデータの最大値を示す周波数を前記特定のスペクトルとして探索すること、を特徴とする回転計である。
請求項３の発明は、請求項１に記載の回転計において、前記特定のスペクトルが、前記測定対象の回転周波数のｘ倍の周波数であるときに、前記探索範囲演算部は、前記特定のスペクトルの周波数をｘで除算した周波数に基づいて探索範囲の幅を演算し、演算した前記探索範囲の幅を前記特定のスペクトルの周波数の探索範囲として設定すること、を特徴とする回転計である。
請求項４の発明は、測定対象の回転数を計測する回転数の計測方法であって、試運転時には、測定対象から計測された回転運動に起因する時系列の振動の検出結果を入力した信号を、高速フーリエ変換してスペクトルデータを演算するＦＦＴ演算ステップと、前記測定対象の主成分を特定のスペクトルとして表示して選択可能なスペクトル選択ステップと、前記ＦＦＴ演算ステップにより演算された前記スペクトルデータのうちの所定のスペクトルの探索範囲を、前記特定のスペクトルの周波数の２倍よりも小さい周波数を上限とし、前記特定のスペクトルの周波数の０．５倍よりも大きい周波数を下限として、演算する探索範囲演算ステップと、を実行し、測定開始時には、測定対象から計測された回転運動に起因する時系列の振動の検出結果を入力した信号を、高速フーリエ変換してスペクトルデータを演算するＦＦＴ演算ステップと、前記測定対象の主成分を特定のスペクトルとして表示して選択可能なスペクトル選択ステップと、前記スペクトル選択ステップで選択した前記特定のスペクトルを、前記ＦＦＴ演算ステップにより演算された前記スペクトルデータのうちの所定のスペクトルの探索範囲から逐次探索するスペクトル探索ステップと、前記スペクトル探索ステップで探索される前記所定スペクトルの探索範囲を、前記特定のスペクトルの周波数の２倍よりも小さい周波数を上限とし、前記特定のスペクトルの周波数の０．５倍よりも大きい周波数を下限として、演算する探索範囲演算ステップと、前記スペクトル探索部によって探索された主成分スペクトルを示す周波数に基づいて、前記測定対象の回転数を演算する回転数演算ステップとを実行すると共に、前記探索範囲演算ステップで、前記特定のスペクトルの周波数が変動したときに、前記探索範囲を、前記特定のスペクトルの変動した周波数に基づいて演算し、変動した前記特定のスペクトルに追従するように実行させること、を特徴とする回転数の計測方法である。
本発明によれば、以下のような効果がある。
（１）選択された特定のスペクトルの周波数に基づいて、特定のスペクトルの探索範囲を演算し、特定のスペクトルの周波数が変動したときに、探索範囲を、変動した特定のスペクトルの周波数に基づいて演算し、変動した特定のスペクトルに追従させるので、特定のスペクトルの周波数が変動したときにおいても、不要なスペクトルを探索範囲から除外することができ、特定のスペクトルを探索ミスすることなく、正確に探索することができる。
（２）スペクトル探索部は、探索範囲内におけるスペクトルデータの最大値を示すスペクトルを特定のスペクトルとして探索するので、探索の判断が容易であり、誤って異なるスペクトルを探索することを防止できる。
（１）選択された特定のスペクトルの周波数に基づいて、特定のスペクトルの探索範囲を演算し、特定のスペクトルの周波数が変動したときに、探索範囲を、変動した特定のスペクトルの周波数に基づいて演算し、変動した特定のスペクトルに追従させるので、特定のスペクトルの周波数が変動したときにおいても、不要なスペクトルを探索範囲から除外することができ、特定のスペクトルを探索ミスすることなく、正確に探索することができる。
（２）スペクトル探索部は、探索範囲内におけるスペクトルデータの最大値を示すスペクトルを特定のスペクトルとして探索するので、探索の判断が容易であり、誤って異なるスペクトルを探索することを防止できる。
（３）探索範囲演算部は、探索範囲の上限を、特定のスペクトルの周波数の２倍よりも小さい周波数として演算し、探索範囲の下限を、特定のスペクトルの周波数の０．５倍よりも大きい周波数として演算するので、特定のスペクトルの０．５次成分及び２次成分を探索範囲から除外することができ、誤って異なるスペクトルを探索するのを防止することができる。
（４）信号処理装置は、測定対象から計測された回転運動に起因する時系列の物理現象の検出結果を信号としているので、容易に、測定対象の回転に起因する物理現象から測定対象の回転数を求めることができる。
（４）信号処理装置は、測定対象から計測された回転運動に起因する時系列の物理現象の検出結果を信号としているので、容易に、測定対象の回転に起因する物理現象から測定対象の回転数を求めることができる。
（５）探索範囲演算部は、特定のスペクトルが、測定対象の回転周波数のｘ倍の周波数であるときに、特定のスペクトルの周波数をｘで除算した周波数に基づいて探索範囲の幅を演算し、演算した探索範囲の幅を特定のスペクトルの周波数の探索範囲として設定するので、特定のスペクトルの周波数の０．５次成分や２次成分などのスペクトルを除外するだけでなく、回転周波数及びその０．５次成分や２次成分などのスペクトルも演算した探索範囲から除外することができ、スペクトル探索部によって正確に特定のスペクトルを探索することができる。
（６）信号処理装置は、測定対象の回転周波数と、特定のスペクトルの周波数との比率を設定し、探索範囲演算部へと入力する設定入力部を備えているので、特定のスペクトルの探索範囲を、特定のスペクトルの周波数及び測定対象の回転周波数に基づいて自動的に演算させることができる。
（６）信号処理装置は、測定対象の回転周波数と、特定のスペクトルの周波数との比率を設定し、探索範囲演算部へと入力する設定入力部を備えているので、特定のスペクトルの探索範囲を、特定のスペクトルの周波数及び測定対象の回転周波数に基づいて自動的に演算させることができる。
本発明は、周波数の変動した特定のスペクトルを正確に読み取ることができる信号処理装置及び信号探索範囲の追従方法を提供するという目的を、変動した特定のスペクトルの周波数に、特定のスペクトルを探索する探索範囲を追従させることにより実現する。
以下、図面等を参照して、本発明の実施例をあげて、さらに詳しく説明する。
図１は、本発明による信号処理装置の実施例１とした回転計を示す図である。図２は、演算制御部で処理されたスペクトルデータを示す図である。図３は、探索範囲演算部により設定された探索範囲を説明する図である。
回転計（信号処理装置）１０は、入力部１１、演算制御部１２、表示部１３、操作部１４、入力アンプ部１５及びフィルタ部１６を有しており、コンプレッサ３０に配置された振動計２０で検出した振動データ（信号）に基づいて、コンプレッサ３０の内部に設けられた回転機構の回転数を測定する測定器である。本実施例では、振動計２０から検出したコンプレッサ３０の振動データに基づいた主成分（１次成分）のスペクトルの周波数ｆ１と、コンプレッサ３０の回転機構の回転周波数は１：１の関係であるとして説明する。
図１は、本発明による信号処理装置の実施例１とした回転計を示す図である。図２は、演算制御部で処理されたスペクトルデータを示す図である。図３は、探索範囲演算部により設定された探索範囲を説明する図である。
回転計（信号処理装置）１０は、入力部１１、演算制御部１２、表示部１３、操作部１４、入力アンプ部１５及びフィルタ部１６を有しており、コンプレッサ３０に配置された振動計２０で検出した振動データ（信号）に基づいて、コンプレッサ３０の内部に設けられた回転機構の回転数を測定する測定器である。本実施例では、振動計２０から検出したコンプレッサ３０の振動データに基づいた主成分（１次成分）のスペクトルの周波数ｆ１と、コンプレッサ３０の回転機構の回転周波数は１：１の関係であるとして説明する。
入力部１１は、振動計２０で検出した時系列の振動データを離散化して、連続的に、回転計１０の演算制御部１２に入力するＡ／Ｄ変換器である。
演算制御部１２は、ＦＦＴ演算部１２ａ、スペクトル選択部１２ｂ、スペクトル探索部１２ｃ、探索範囲演算部１２ｄ及び回転数演算部１２ｅを有し、回転計１０の各部を統括制御する制御回路である。また、後述の操作部１４によって測定時間を設定することができる。
ＦＦＴ演算部１２ａは、入力部１１から連続的に入力した離散化された振動データを、高速フーリエ変換（ＦＦＴ）し、振動のスペクトルデータを演算する演算回路である。
演算制御部１２は、ＦＦＴ演算部１２ａ、スペクトル選択部１２ｂ、スペクトル探索部１２ｃ、探索範囲演算部１２ｄ及び回転数演算部１２ｅを有し、回転計１０の各部を統括制御する制御回路である。また、後述の操作部１４によって測定時間を設定することができる。
ＦＦＴ演算部１２ａは、入力部１１から連続的に入力した離散化された振動データを、高速フーリエ変換（ＦＦＴ）し、振動のスペクトルデータを演算する演算回路である。
スペクトル選択部１２ｂは、測定対象であるコンプレッサ３０の初期動作時の主成分のスペクトル（特定のスペクトル）を選択する回路である。スペクトル選択部１２ｂは、ＦＦＴ演算部１２ａでＦＦＴ演算したスペクトルデータに基づいて、コンプレッサ３０の初期動作時の回転数に起因する主成分のスペクトルを作業者に選択させる。具体的には、図２に示すように、ＦＦＴ演算部１２ａで演算したスペクトルデータを表示部１３に表示し、表示されたスペクトルデータから主成分のスペクトルを作業者に選択させる。また、コンプレッサ３０の初期動作時の回転数が予め判明しているときは、後述の操作部１４から主成分のスペクトルの周波数ｆ１を入力することもできる。
スペクトル探索部１２ｃは、コンプレッサ３０の回転数の測定中に、ＦＦＴ演算部１２ａで演算されるスペクトルデータから、スペクトル選択部１２ｂで選択された主成分のスペクトルを逐次探索する回路である。ここで、主成分のスペクトルは、スペクトルデータの全体から探索されるのではなく、後述の探索範囲演算部１２ｄで演算した探索範囲から探索される。また、主成分のスペクトルの探索は、その探索範囲内で最大値を示すスペクトルを、主成分のスペクトルとして判断している。
探索範囲演算部１２ｄは、スペクトル探索部１２ｃで探索する主成分のスペクトルの探索範囲を演算する回路である。具体的には、演算したスペクトルデータは、図２に示すように、主成分（１次成分）のスペクトルの他に、０．５次成分、２次成分及び３次成分などのスペクトルが検出される。
本実施例では、主成分と隣り合う０．５次成分及び２次成分のスペクトルを、主成分のスペクトルと誤って読み取らないようにするために、探索範囲を０．５次成分及び２次成分のスペクトルを除外するように自動的に設定する。具体的には、図３（ａ）に示すように、コンプレッサ３０の初期動作時の主成分（１次成分）のスペクトルの周波数ｆ１が２ｋＨｚであるときに、０．５次成分及び２次成分のスペクトルはそれぞれ１ｋＨｚ及び４ｋＨｚとして現れる。よって、探索範囲は、この０．５次及び２次成分を含まない範囲を設定する必要があり、例えば、下限の周波数ｆLを周波数ｆ１の４割減（ｆ１の０．６倍
）に設定し（下記の式１参照）、上限の周波数ｆHを周波数ｆ１の８割増（ｆ１の１．８
倍）に設定することによって（下記の式２参照）、０．５次成分及び２次成分のスペクトルを探索範囲から除外している。
本実施例では、主成分と隣り合う０．５次成分及び２次成分のスペクトルを、主成分のスペクトルと誤って読み取らないようにするために、探索範囲を０．５次成分及び２次成分のスペクトルを除外するように自動的に設定する。具体的には、図３（ａ）に示すように、コンプレッサ３０の初期動作時の主成分（１次成分）のスペクトルの周波数ｆ１が２ｋＨｚであるときに、０．５次成分及び２次成分のスペクトルはそれぞれ１ｋＨｚ及び４ｋＨｚとして現れる。よって、探索範囲は、この０．５次及び２次成分を含まない範囲を設定する必要があり、例えば、下限の周波数ｆLを周波数ｆ１の４割減（ｆ１の０．６倍
）に設定し（下記の式１参照）、上限の周波数ｆHを周波数ｆ１の８割増（ｆ１の１．８
倍）に設定することによって（下記の式２参照）、０．５次成分及び２次成分のスペクトルを探索範囲から除外している。
ｆL＝ｆ１−０．４ｆ１＝０．６ｆ１ （式１）
ｆH＝ｆ１＋０．８ｆ１＝１．８ｆ１ （式２）
ｆH＝ｆ１＋０．８ｆ１＝１．８ｆ１ （式２）
よって、探索範囲演算部１２ｄは、スペクトル選択部１２ｂで選択された主成分のスペクトルの周波数ｆ１に基づいて、上記の式１及び式２を用いて探索範囲を演算する。例えば、図３（ａ）に示すように、主成分のスペクトルの周波数ｆ１が２ｋＨｚであるとき、探索範囲の下限の周波数ｆLは、１．２ｋＨｚとなり、上限の周波数ｆHは、３．６ｋＨｚとなる。
また、上記の演算により、探索範囲演算部１２ｄは、コンプレッサ３０の初期動作時の主成分のスペクトルの周波数ｆ１（回転数）が変動したときに、上記の式１及び式２によって、主成分のスペクトルの周波数ｆ１の変動に探索範囲を追従させることができる。例えば、図３（ｂ）に示すように、主成分のスペクトルの周波数ｆ１が２ｋＨｚから２．２
ｋＨｚ（ｆ１‘）に変動したとき、０．５次成分は１．１ｋＨｚとなり、２次成分は４．４ｋＨｚとなり、上述の探索範囲（１．２〜３．６ｋＨｚ）から０．５次成分及び２次成分を除外している。そして、探索範囲演算部１２ｄは、変動した周波数（ｆ１‘＝２．２ｋＨｚ）に基づいて、変動後の探索範囲の下限の周波数ｆL‘及び上限の周波数ｆH‘をｆ１‘に追従させるように、式１及び式２により演算する。ここでは、変動後の下限の周波数ｆL‘が１．３２ｋＨｚとなり、上限の周波数ｆH‘が３．９６ｋＨｚとなる。
変動した主成分のスペクトルの周波数ｆ１‘（２．２ｋＨｚ）がさらに変動したとき（例えば、２．４ｋＨｚ）、その変動した周波数の０．５次成分（１．２ｋＨｚ）及び２次成分（４．８ｋＨｚ）は、ｆ１‘に追従させた探索範囲（１．３２（ｆL‘）〜３．９６
ｋＨｚ（ｆH‘））によって除外される。ここで、図２及び図３のスペクトルデータは、
横軸を周波数ｆ（Ｈｚ）、縦軸をパワースペクトル（ＰＳ）としている。
また、上記の演算により、探索範囲演算部１２ｄは、コンプレッサ３０の初期動作時の主成分のスペクトルの周波数ｆ１（回転数）が変動したときに、上記の式１及び式２によって、主成分のスペクトルの周波数ｆ１の変動に探索範囲を追従させることができる。例えば、図３（ｂ）に示すように、主成分のスペクトルの周波数ｆ１が２ｋＨｚから２．２
ｋＨｚ（ｆ１‘）に変動したとき、０．５次成分は１．１ｋＨｚとなり、２次成分は４．４ｋＨｚとなり、上述の探索範囲（１．２〜３．６ｋＨｚ）から０．５次成分及び２次成分を除外している。そして、探索範囲演算部１２ｄは、変動した周波数（ｆ１‘＝２．２ｋＨｚ）に基づいて、変動後の探索範囲の下限の周波数ｆL‘及び上限の周波数ｆH‘をｆ１‘に追従させるように、式１及び式２により演算する。ここでは、変動後の下限の周波数ｆL‘が１．３２ｋＨｚとなり、上限の周波数ｆH‘が３．９６ｋＨｚとなる。
変動した主成分のスペクトルの周波数ｆ１‘（２．２ｋＨｚ）がさらに変動したとき（例えば、２．４ｋＨｚ）、その変動した周波数の０．５次成分（１．２ｋＨｚ）及び２次成分（４．８ｋＨｚ）は、ｆ１‘に追従させた探索範囲（１．３２（ｆL‘）〜３．９６
ｋＨｚ（ｆH‘））によって除外される。ここで、図２及び図３のスペクトルデータは、
横軸を周波数ｆ（Ｈｚ）、縦軸をパワースペクトル（ＰＳ）としている。
回転数演算部１２ｅは、スペクトル探索部１２ｃによって探索された主成分のスペクトルを示す周波数ｆ１に基づいて、コンプレッサ３０の回転数を演算し、演算結果を表示部１３に表示させる回路である。本実施例では、探索された主成分の周波数ｆ１（Ｈｚ：１秒当たりの振動数）を、分単位の回転数（ｒ／ｍｉｎ）に演算している。
操作部１４は、作業者が回転計１０に対して操作内容を指示するテンキー及び複数のボタンを備えており、主成分のスペクトルの周波数ｆ１が予め判明しているときにその周波数ｆ１値を入力したり、測定時間を設定したりすることができる。
操作部１４は、作業者が回転計１０に対して操作内容を指示するテンキー及び複数のボタンを備えており、主成分のスペクトルの周波数ｆ１が予め判明しているときにその周波数ｆ１値を入力したり、測定時間を設定したりすることができる。
入力アンプ部１５は、振動計２０から検出された振動データを増幅する増幅器である。
フィルタ部１６は、入力アンプ部１５で増幅した振動データに含まれる折り返し雑音を除去するＡＡＦ（アンチ・エイリアシング・フィルタ）である。
振動計２０は、コンプレッサ３０の筐体に固定され、コンプレッサ３０の回転機構の振動を加速度の振動データとして検出する加速度計である。
コンプレッサ３０は、内部に回転機構を備え、その回転機構により気体又は液体を圧縮させて圧力を高める装置である。
フィルタ部１６は、入力アンプ部１５で増幅した振動データに含まれる折り返し雑音を除去するＡＡＦ（アンチ・エイリアシング・フィルタ）である。
振動計２０は、コンプレッサ３０の筐体に固定され、コンプレッサ３０の回転機構の振動を加速度の振動データとして検出する加速度計である。
コンプレッサ３０は、内部に回転機構を備え、その回転機構により気体又は液体を圧縮させて圧力を高める装置である。
次に、回転計１０のコンプレッサ３０の回転数測定の準備について説明する。図４は、回転計１０の測定準備の流れを示すフローである。
作業者により振動計２０がコンプレッサ３０に対して適正に配置され、コンプレッサ３０の回転数の測定準備を開始したら（Ｓ１０１）、作業者は、回転計１０の電源を投入する（Ｓ１０２）。電源が投入されたら、コンプレッサ３０の初期動作時の主成分のスペクトルの周波数ｆ１の値を入力するか否かを確認する（Ｓ１０３）。
初期動作時の主成分のスペクトルの周波数ｆ１の値が予め判明しておらず、入力できない場合（Ｓ１０３：Ｎｏ）、コンプレッサ３０を試運転する旨の指示を作業者に対して表示部１３を介して表示する。コンプレッサ３０が試運転されたら、演算制御部１２は、振動計２０によって検出したコンプレッサ３０の試運転（初期動作）時の振動データを、ＦＦＴ演算部１２ａに入力してスペクトルデータを演算する（Ｓ１０４）。続いて、演算したスペクトルデータを、図２に示すように、表示部１３に表示して、作業者に、コンプレッサ３０の主成分のスペクトル（２ｋＨｚ）を選択させる（Ｓ１０５）。
作業者により振動計２０がコンプレッサ３０に対して適正に配置され、コンプレッサ３０の回転数の測定準備を開始したら（Ｓ１０１）、作業者は、回転計１０の電源を投入する（Ｓ１０２）。電源が投入されたら、コンプレッサ３０の初期動作時の主成分のスペクトルの周波数ｆ１の値を入力するか否かを確認する（Ｓ１０３）。
初期動作時の主成分のスペクトルの周波数ｆ１の値が予め判明しておらず、入力できない場合（Ｓ１０３：Ｎｏ）、コンプレッサ３０を試運転する旨の指示を作業者に対して表示部１３を介して表示する。コンプレッサ３０が試運転されたら、演算制御部１２は、振動計２０によって検出したコンプレッサ３０の試運転（初期動作）時の振動データを、ＦＦＴ演算部１２ａに入力してスペクトルデータを演算する（Ｓ１０４）。続いて、演算したスペクトルデータを、図２に示すように、表示部１３に表示して、作業者に、コンプレッサ３０の主成分のスペクトル（２ｋＨｚ）を選択させる（Ｓ１０５）。
コンプレッサ３０の初期動作時の主成分のスペクトルの周波数ｆ１の値が測定前に予め判明しており、その値を入力する場合（例えば、ｆ１＝２ｋＨｚ、Ｓ１０３：Ｙｅｓ）、演算制御部１２は、作業者に、操作部１４から主成分のスペクトルの周波数ｆ１の値を入力させる（Ｓ１０６）。
コンプレッサ３０の主成分のスペクトルが選択又は入力されたら、演算制御部１２の探索範囲演算部１２ｄは、その主成分のスペクトルの周波数ｆ１に基づいて、測定開始時における主成分のスペクトルを探索する探索範囲（本実施例では１．２ｋＨｚ〜３．６ｋＨｚ）を演算して求め（Ｓ１０７）、測定準備を完了させる（Ｓ１０８）。
コンプレッサ３０の主成分のスペクトルが選択又は入力されたら、演算制御部１２の探索範囲演算部１２ｄは、その主成分のスペクトルの周波数ｆ１に基づいて、測定開始時における主成分のスペクトルを探索する探索範囲（本実施例では１．２ｋＨｚ〜３．６ｋＨｚ）を演算して求め（Ｓ１０７）、測定準備を完了させる（Ｓ１０８）。
次に、回転計１０がコンプレッサ３０の回転数を測定するときの動作について説明する。図５は、コンプレッサ３０の測定開始から終了までの動作の流れを示すフローである。
上述のように測定準備が完了したら、回転計１０は、コンプレッサ３０の回転数の測定を開始する（Ｓ２０１）。まず、コンプレッサ３０を動作させて振動計２０によって連続的に振動データを検出させたら、回転計１０の演算制御部１２は、検出した振動データを入力してＦＦＴ演算部１２ａによりスペクトルデータを演算する（Ｓ２０２）。次に、図３（ａ）に示すように、測定準備時に求めた探索範囲内（１．２ｋＨｚ〜３．６ｋＨｚ）で主成分のスペクトルを、スペクトル探索部１２ｃによって探索する（Ｓ２０３）。主成分のスペクトルを探索したら、演算制御部１２は、探索した主成分のスペクトルの周波数ｆ１［Ｈｚ］を、回転数演算部１２ｅによって、分単位の回転数（ｒ／ｍｉｎ）へと変換して表示部１３に表示する（Ｓ２０４）。
上述のように測定準備が完了したら、回転計１０は、コンプレッサ３０の回転数の測定を開始する（Ｓ２０１）。まず、コンプレッサ３０を動作させて振動計２０によって連続的に振動データを検出させたら、回転計１０の演算制御部１２は、検出した振動データを入力してＦＦＴ演算部１２ａによりスペクトルデータを演算する（Ｓ２０２）。次に、図３（ａ）に示すように、測定準備時に求めた探索範囲内（１．２ｋＨｚ〜３．６ｋＨｚ）で主成分のスペクトルを、スペクトル探索部１２ｃによって探索する（Ｓ２０３）。主成分のスペクトルを探索したら、演算制御部１２は、探索した主成分のスペクトルの周波数ｆ１［Ｈｚ］を、回転数演算部１２ｅによって、分単位の回転数（ｒ／ｍｉｎ）へと変換して表示部１３に表示する（Ｓ２０４）。
次に、演算制御部１２は、継続的に電源が投入された状態であるときに（Ｓ２０５：Ｎｏ）、探索した主成分のスペクトルの周波数ｆ１が、変動しているか否かを確認する（Ｓ２０６）。周波数が変動している場合は（Ｓ２０６：有）、変動した周波数ｆ１‘（例えば、２．２ｋＨｚ）に基づいて探索範囲（例えば、１．３２ｋＨｚ〜３．９６ｋＨｚ）を演算し（Ｓ２０７）、図３（ｂ）に示すように、演算した探索範囲を用いて、ＦＦＴ演算部１２ａで演算したスペクトルデータから主成分のスペクトルを探索し（Ｓ２０２、Ｓ２０３）、回転数を表示する（Ｓ２０４）。
仮に、周波数ｆ１が変動していない場合は（Ｓ２０６：無）、引き続き初期動作時と同じ探索範囲を用いて、ＦＦＴ演算部１２ａで演算したスペクトルデータから主成分のスペクトルを探索し（Ｓ２０２、Ｓ２０３）、回転数を表示する（Ｓ２０４）。
電源が遮断されたら（Ｓ２０５：Ｙｅｓ）、コンプレッサ３０の回転数の測定は終了となる（Ｓ２０８）。
電源が遮断されたら（Ｓ２０５：Ｙｅｓ）、コンプレッサ３０の回転数の測定は終了となる（Ｓ２０８）。
以上より、実施例１の回転計には以下のような効果がある。
（１）回転計１０は、選択された主成分のスペクトルの周波数ｆ１に基づいて、主成分のスペクトルの探索範囲を演算し、主成分のスペクトルの周波数ｆ１が変動したときに、探索範囲を変動した主成分のスペクトルの周波数ｆ１‘に基づいて演算し、変動した主成分のスペクトルに追従させるので、主成分のスペクトルの周波数ｆ１が変動したときにおいても、０．５次成分や２次成分などの不要なスペクトルを探索範囲から除外することができ、主成分のスペクトルの探索ミスを発生させることなく、正確に探索することができる。
（２）スペクトル探索部１２ｃは、探索範囲内におけるスペクトルデータの最大値を示すスペクトルを主成分のスペクトルとして探索するので、探索の判断が容易であり、誤って０．５次成分や２次成分などの異なったスペクトルを探索することを防止できる。
（１）回転計１０は、選択された主成分のスペクトルの周波数ｆ１に基づいて、主成分のスペクトルの探索範囲を演算し、主成分のスペクトルの周波数ｆ１が変動したときに、探索範囲を変動した主成分のスペクトルの周波数ｆ１‘に基づいて演算し、変動した主成分のスペクトルに追従させるので、主成分のスペクトルの周波数ｆ１が変動したときにおいても、０．５次成分や２次成分などの不要なスペクトルを探索範囲から除外することができ、主成分のスペクトルの探索ミスを発生させることなく、正確に探索することができる。
（２）スペクトル探索部１２ｃは、探索範囲内におけるスペクトルデータの最大値を示すスペクトルを主成分のスペクトルとして探索するので、探索の判断が容易であり、誤って０．５次成分や２次成分などの異なったスペクトルを探索することを防止できる。
（３）探索範囲演算部１２ｄは、探索範囲の上限を、主成分のスペクトルの周波数ｆ１の１．８倍の周波数ｆHとして演算し、探索範囲の下限を、主成分のスペクトルの周波数ｆ
１の０．６倍の周波数ｆLとして演算するので、主成分のスペクトルの０．５次成分及び
２次成分を探索範囲から除外することができ、誤って異なるスペクトルを探索するのを防止することができる。
（４）信号処理装置１０は、コンプレッサ３０で計測された回転運動に起因する時系列の振動データを検出しているので、コンプレッサ３０上に配置した振動計２０の振動データから容易に、コンプレッサ３０の回転数を測定することができる。
１の０．６倍の周波数ｆLとして演算するので、主成分のスペクトルの０．５次成分及び
２次成分を探索範囲から除外することができ、誤って異なるスペクトルを探索するのを防止することができる。
（４）信号処理装置１０は、コンプレッサ３０で計測された回転運動に起因する時系列の振動データを検出しているので、コンプレッサ３０上に配置した振動計２０の振動データから容易に、コンプレッサ３０の回転数を測定することができる。
図６は、本発明の実施例２における探索範囲演算部による探索範囲の説明をする図であ
る。実施例２の説明では、前述した実施例１と同様な機能を果たす部分には、同一の符号又は末尾に統一した符号を付して、重複する説明や図面を適宜省略する。
る。実施例２の説明では、前述した実施例１と同様な機能を果たす部分には、同一の符号又は末尾に統一した符号を付して、重複する説明や図面を適宜省略する。
実施例２の信号処理装置も、実施例１と同様に、回転計を例にしており、実施例２の回転計１０−２と実施例１の回転計１０との相違点は、実施例２の測定対象が自動車用の４気筒の４サイクル式のエンジン３０−２となる点と、実施例２の演算制御部１２の探索範囲演算部１２ｄに、探索範囲の演算方法に追加機能が備えられ、探索範囲演算部１２ｄ−２となる点である。
回転計１０−２は、入力部１１、演算制御部１２、表示部１３、操作部１４、入力アンプ部１５及びフィルタ部１６を有しており、エンジン３０−２上に配置された振動計２０で検出した振動データ（信号）に基づいて、エンジン３０−２の回転数を測定する測定器である。
回転計１０−２は、入力部１１、演算制御部１２、表示部１３、操作部１４、入力アンプ部１５及びフィルタ部１６を有しており、エンジン３０−２上に配置された振動計２０で検出した振動データ（信号）に基づいて、エンジン３０−２の回転数を測定する測定器である。
エンジン３０−２は、上述したように、自動車用の４気筒の４サイクル式の内燃機関である。エンジン３０−２は、各気筒内に燃料と空気を混合させた混合気を吸入し、気筒内に設けられたピストンによって吸入した混合気を圧縮し、点火、燃焼、膨張させることによってピストンを往復運動させ、ピストンに連動したクランクシャフトを回転させている。一般に、４サイクルエンジンは、クランクシャフトが２回転する間に１回の爆発（燃焼、膨張）が気筒内で発生する。４気筒の４サイクルエンジンの場合、各気筒の爆発が互い違いに発生するので、エンジン３０−２全体から見ると、クランクシャフト１回転に対して２回の爆発が発生する。ここで、エンジン３０−２の回転数（回転周波数）とは、クランクシャフトの回転数（回転周波数）のことである。
エンジン３０−２上に配置された振動計２０は、前述のエンジン３０−２で発生する爆発の振動を主に検出するが、クランクシャフトの回転に伴う振動も、同じエンジン３０−２で発生する振動であるので微小ではあるが検出する。
エンジン３０−２上に配置された振動計２０は、前述のエンジン３０−２で発生する爆発の振動を主に検出するが、クランクシャフトの回転に伴う振動も、同じエンジン３０−２で発生する振動であるので微小ではあるが検出する。
探索範囲演算部１２ｄ−２に設けられた追加機能は、エンジン３０−２の爆発により発生する振動の主成分のスペクトルの周波数ｆ１と、エンジン３０−２の回転周波数ｆEと
の関係が、１：１以外の比率（本実施例では上述のように２：１）であるときに、その比率を主成分のスペクトルの探索範囲の演算に反映させる機能である。
の関係が、１：１以外の比率（本実施例では上述のように２：１）であるときに、その比率を主成分のスペクトルの探索範囲の演算に反映させる機能である。
具体的には、エンジン３０−２のクランクシャフトが１ｋＨｚで回転しているときに、振動計２０によって検出されるエンジン３０−２の爆発の振動に基づいたスペクトルデータは、２ｋＨｚのスペクトルが主成分として現れ、０．５次成分として１ｋＨｚ、２次成分として４ｋＨｚのスペクトルも現れる。また、クランクシャフトの回転振動に基づいたスペクトルデータは、１ｋＨｚのスペクトルが主成分として現れ、２次、３次及び４次成分として２ｋＨｚ、３ｋＨｚ及び４ｋＨｚのスペクトルも現れる。この爆発及び回転による振動のスペクトルが合成されて、図６に示すようなスペクトルデータとして現れる。
上述したように、振動計２０は、クランクシャフトの回転に伴う振動よりも、爆発による振動を感度よく検出するので、スペクトルデータは、図６に示すように、爆発の振動に基づいたスペクトルは大きく現れ、クランクシャフトの回転に基づいた振動は、爆発のスペクトルよりも小さいスペクトルとして現れる。すなわち、２ｋＨｚのスペクトルは、エンジン３０−２の爆発の主成分及び回転の２次成分を含んでいるため、爆発の０．５次成分及び回転の主成分を含んだ１ｋＨｚのスペクトルよりも振幅が大きくなる傾向となる。したがって、２ｋＨｚのスペクトルを主成分のスペクトルとしてスペクトル探索部１２ｃに探索させるのが、最も探索効率がよく、また、正確に主成分のスペクトルを探索することができる。
前述のように、２ｋＨｚ（ｆ１）のスペクトルを主成分のスペクトルとして探索させる
ときに、実施例１に記載されている式１及び式２の演算で探索範囲を求めると、下限の周波数が１．２ｋＨｚ、上限の周波数が３．６ｋＨｚとなって求められる。しかし、この探索範囲（１．２ｋＨｚ〜３．６ｋＨｚ）には、クランクシャフトの回転に起因する回転周波数ｆE（１ｋＨｚ）のスペクトルの２次成分（２ｋＨｚ）及び３次成分（３ｋＨｚ）が
含まれてしまい、スペクトル探索部１２ｃ−２による主成分のスペクトルの探索ミスを発生させる可能性がある。
ときに、実施例１に記載されている式１及び式２の演算で探索範囲を求めると、下限の周波数が１．２ｋＨｚ、上限の周波数が３．６ｋＨｚとなって求められる。しかし、この探索範囲（１．２ｋＨｚ〜３．６ｋＨｚ）には、クランクシャフトの回転に起因する回転周波数ｆE（１ｋＨｚ）のスペクトルの２次成分（２ｋＨｚ）及び３次成分（３ｋＨｚ）が
含まれてしまい、スペクトル探索部１２ｃ−２による主成分のスペクトルの探索ミスを発生させる可能性がある。
そこで、探索範囲演算部１２ｄ−２は、エンジン３０−２の回転周波数（１ｋＨｚ）に起因した探索範囲の幅Ｆを演算し、その幅Ｆを主成分のスペクトル（２ｋＨｚ）に適用するように演算する。具体的には、探索範囲の下限の周波数ｆL及び上限の周波数ｆHは、以下に示す式３及び式４により演算される。
ｆL＝ｆ１−０．４ｆE （式３）
ｆH＝ｆ１＋０．８ｆE （式４）
ｆH＝ｆ１＋０．８ｆE （式４）
式３及び式４によって、主成分のスペクトルの周波数ｆ１に対して、エンジン３０−２の回転周波数ｆEに基づいて演算された探索範囲の幅Ｆを用いることができ、効率よく正
確に主成分のスペクトルを探索することができる。本実施例では、図６に示すように、主成分のスペクトルの周波数ｆ１が２ｋＨｚであり、エンジン３０−２の回転周波数ｆEは
１ｋＨｚであるので、式３及び式４によって主成分の探索範囲の下限の周波数ｆLが１．
６ｋＨｚと演算され、上限の周波数ｆHが２．８ｋＨｚと演算される。
また、探索範囲演算部１２ｄ−２は、エンジン３０−２の爆発により発生する振動の主成分のスペクトルの周波数ｆ１と、エンジン３０−２の回転周波数ｆEとの比率ｘ（本実
施例ではｘ＝２）の値を作業者によって操作部１４から入力できるようにしており、例えば、以下に示す式５と、上述の式３及び式４によって探索範囲を自動的に演算させることができる。
確に主成分のスペクトルを探索することができる。本実施例では、図６に示すように、主成分のスペクトルの周波数ｆ１が２ｋＨｚであり、エンジン３０−２の回転周波数ｆEは
１ｋＨｚであるので、式３及び式４によって主成分の探索範囲の下限の周波数ｆLが１．
６ｋＨｚと演算され、上限の周波数ｆHが２．８ｋＨｚと演算される。
また、探索範囲演算部１２ｄ−２は、エンジン３０−２の爆発により発生する振動の主成分のスペクトルの周波数ｆ１と、エンジン３０−２の回転周波数ｆEとの比率ｘ（本実
施例ではｘ＝２）の値を作業者によって操作部１４から入力できるようにしており、例えば、以下に示す式５と、上述の式３及び式４によって探索範囲を自動的に演算させることができる。
ｆE＝ｆ１／ｘ （式５）
この場合、実施例１のように測定対象の回転周波数ｆEと、振動計で検出した振動デー
タの主成分のスペクトルの周波数ｆ１とが１：１であるときは、ｘ＝１を入力することで、実施例１と同様な演算結果を求めることができる。
以上のように、本実施例の回転計１０−２は、測定対象の回転運動に起因する物理現象の信号（振動データ）の周波数と、測定対象の回転周波数との関係が特定の比率ｘになっている場合（本実施例ではｘ＝２）に特に用いることができ、周波数ｆ１の探索ミスを発生させることなく、正確に回転数を測定することができる。
タの主成分のスペクトルの周波数ｆ１とが１：１であるときは、ｘ＝１を入力することで、実施例１と同様な演算結果を求めることができる。
以上のように、本実施例の回転計１０−２は、測定対象の回転運動に起因する物理現象の信号（振動データ）の周波数と、測定対象の回転周波数との関係が特定の比率ｘになっている場合（本実施例ではｘ＝２）に特に用いることができ、周波数ｆ１の探索ミスを発生させることなく、正確に回転数を測定することができる。
以上より、実施例２の回転計には以下のような効果がある。
（１）探索範囲演算部１２ｄ−２は、振動計２０で検出した振動データの主成分のスペクトルが、エンジン３０−２の回転周波数ｆEの２倍の周波数であるときに、主成分のスペ
クトルの周波数ｆ１を２で除算した周波数に基づいて探索範囲の幅Ｆを演算し、演算した探索範囲の幅Ｆを主成分のスペクトルの周波数ｆ１の探索範囲として設定するので、主成分のスペクトルの周波数ｆ１の０．５次成分及び２次成分のスペクトルを除外するだけでなく、エンジン３０−２の回転周波数ｆEの２次成分や３次成分のスペクトルも演算した
探索範囲から除外することができ、スペクトル探索部１２ｃによって正確に主成分のスペクトルを探索することができる。
（２）回転計１０−２は、エンジン３０−２の回転周波数ｆEと、主成分のスペクトルの
周波数ｆ１との比率ｘを設定し、探索範囲演算部１２ｄ−２へと入力する操作部１４を備えているので、主成分のスペクトルの探索を、エンジン３０−２の回転周波数ｆEに基づ
いて自動的に演算させることができる。
（１）探索範囲演算部１２ｄ−２は、振動計２０で検出した振動データの主成分のスペクトルが、エンジン３０−２の回転周波数ｆEの２倍の周波数であるときに、主成分のスペ
クトルの周波数ｆ１を２で除算した周波数に基づいて探索範囲の幅Ｆを演算し、演算した探索範囲の幅Ｆを主成分のスペクトルの周波数ｆ１の探索範囲として設定するので、主成分のスペクトルの周波数ｆ１の０．５次成分及び２次成分のスペクトルを除外するだけでなく、エンジン３０−２の回転周波数ｆEの２次成分や３次成分のスペクトルも演算した
探索範囲から除外することができ、スペクトル探索部１２ｃによって正確に主成分のスペクトルを探索することができる。
（２）回転計１０−２は、エンジン３０−２の回転周波数ｆEと、主成分のスペクトルの
周波数ｆ１との比率ｘを設定し、探索範囲演算部１２ｄ−２へと入力する操作部１４を備えているので、主成分のスペクトルの探索を、エンジン３０−２の回転周波数ｆEに基づ
いて自動的に演算させることができる。
（変形例）
以上、説明した実施例に限定されることなく、種々の変形や変更が可能であり、それらも本発明の範囲内である。
（１）各実施例では、回転数を検出する手段として振動計２０で測定した振動データを用いたが、それ以外のものを使用してもよい。例えば、実施例２のエンジン３０−２の回転数は、エンジン３０−２の振動ではなく、音圧計を用いて、エンジン３０−２の燃焼ガスの排気音を検出した音圧データを用いることも可能である。
（２）各実施例では、信号処理装置として回転計を例にしたが、それ以外の測定器などに使用してもよい。
（３）実施例２では、回転計１０−２は、エンジン３０−２の回転数を測定したが、それ以外の測定対象に用いてもよい。例えば、４極の回転モータを測定対象にし、回転モータの回転子の磁気漏洩を磁気センサで検出して、その検出データに基づいて回転数を演算させることも可能である。この場合、磁気センサが検出する磁気漏洩の信号は、回転モータの回転子が４極であるので、１回転当たりに４回検出される。すなわち、磁気センサの信号の周波数は、回転モータの回転周波数の４倍（ｘ＝４）となり、上述の式３〜５によって、探索範囲を演算することができる。
以上、説明した実施例に限定されることなく、種々の変形や変更が可能であり、それらも本発明の範囲内である。
（１）各実施例では、回転数を検出する手段として振動計２０で測定した振動データを用いたが、それ以外のものを使用してもよい。例えば、実施例２のエンジン３０−２の回転数は、エンジン３０−２の振動ではなく、音圧計を用いて、エンジン３０−２の燃焼ガスの排気音を検出した音圧データを用いることも可能である。
（２）各実施例では、信号処理装置として回転計を例にしたが、それ以外の測定器などに使用してもよい。
（３）実施例２では、回転計１０−２は、エンジン３０−２の回転数を測定したが、それ以外の測定対象に用いてもよい。例えば、４極の回転モータを測定対象にし、回転モータの回転子の磁気漏洩を磁気センサで検出して、その検出データに基づいて回転数を演算させることも可能である。この場合、磁気センサが検出する磁気漏洩の信号は、回転モータの回転子が４極であるので、１回転当たりに４回検出される。すなわち、磁気センサの信号の周波数は、回転モータの回転周波数の４倍（ｘ＝４）となり、上述の式３〜５によって、探索範囲を演算することができる。
１０、１０−２ 回転計
１１ 入力部
１２ 演算制御部
１２ａ ＦＦＴ演算部
１２ｂ スペクトル選択部
１２ｃ スペクトル探索部
１２ｄ、１２ｄ−２ 探索範囲演算部
１２ｅ 回転数演算部
１３ 表示部
１４ 操作部
２０ 振動計
３０ コンプレッサ
３０−２ エンジン
１１ 入力部
１２ 演算制御部
１２ａ ＦＦＴ演算部
１２ｂ スペクトル選択部
１２ｃ スペクトル探索部
１２ｄ、１２ｄ−２ 探索範囲演算部
１２ｅ 回転数演算部
１３ 表示部
１４ 操作部
２０ 振動計
３０ コンプレッサ
３０−２ エンジン
Claims (4)
- 測定対象から計測された回転運動に起因する時系列の振動の検出結果を入力した信号を、高速フーリエ変換してスペクトルデータを演算するＦＦＴ演算部と、
前記測定対象の主成分を特定のスペクトルとして表示して選択可能なスペクトル選択部と、
前記スペクトル選択部で選択した前記特定のスペクトルを、前記ＦＦＴ演算部により演算された前記スペクトルデータのうちの所定のスペクトルの探索範囲から逐次探索するスペクトル探索部と、
前記スペクトル探索部で探索される前記所定のスペクトルの探索範囲を、前記特定のスペクトルの周波数の２倍よりも小さい周波数を上限とし、前記特定のスペクトルの周波数の０．５倍よりも大きい周波数を下限として、演算する探索範囲演算部と、
前記スペクトル探索部によって探索された主成分スペクトルを示す周波数に基づいて、前記測定対象の回転数を演算する回転数演算部と、
を含む演算制御部を備え、
前記演算制御部は、
試運転時には、前記ＦＦＴ演算部、前記スペクトル選択部及び前記探索範囲演算部を動作させ、
測定開始時には、前記ＦＦＴ演算部、前記スペクトル探索部、前記探索範囲演算部及び前記回転数演算部を動作させると共に、前記探索範囲演算部を、前記特定のスペクトルの周波数が変動したときに、前記探索範囲を、前記特定のスペクトルの変動した周波数に基づいて演算し、変動した前記特定のスペクトルに追従するように動作させること、
を特徴とする回転計。 - 請求項１に記載の回転計において、
前記スペクトル探索部は、前記探索範囲内における前記スペクトルデータの最大値を示す周波数を前記特定のスペクトルとして探索すること、
を特徴とする回転計。 - 請求項１に記載の回転計において、
前記特定のスペクトルが、前記測定対象の回転周波数のｘ倍の周波数であるときに、
前記探索範囲演算部は、前記特定のスペクトルの周波数をｘで除算した周波数に基づいて探索範囲の幅を演算し、演算した前記探索範囲の幅を前記特定のスペクトルの周波数の探索範囲として設定すること、
を特徴とする回転計。 - 測定対象の回転数を計測する回転数の計測方法であって、
試運転時には、
測定対象から計測された回転運動に起因する時系列の振動の検出結果を入力した信号を、高速フーリエ変換してスペクトルデータを演算するＦＦＴ演算ステップと、
前記測定対象の主成分を特定のスペクトルとして表示して選択可能なスペクトル選択ステップと、
前記ＦＦＴ演算ステップにより演算された前記スペクトルデータのうちの所定のスペクトルの探索範囲を、前記特定のスペクトルの周波数の２倍よりも小さい周波数を上限とし、前記特定のスペクトルの周波数の０．５倍よりも大きい周波数を下限として、演算する探索範囲演算ステップと、
を実行し、
測定開始時には、
測定対象から計測された回転運動に起因する時系列の振動の検出結果を入力した信号を、高速フーリエ変換してスペクトルデータを演算するＦＦＴ演算ステップと、
前記測定対象の主成分を特定のスペクトルとして表示して選択可能なスペクトル選択ステップと、
前記スペクトル選択ステップで選択した前記特定のスペクトルを、前記ＦＦＴ演算ステップにより演算された前記スペクトルデータのうちの所定のスペクトルの探索範囲から逐次探索するスペクトル探索ステップと、
前記スペクトル探索ステップで探索される前記所定スペクトルの探索範囲を、前記特定のスペクトルの周波数の２倍よりも小さい周波数を上限とし、前記特定のスペクトルの周波数の０．５倍よりも大きい周波数を下限として、演算する探索範囲演算ステップと、
前記スペクトル探索部によって探索された主成分スペクトルを示す周波数に基づいて、前記測定対象の回転数を演算する回転数演算ステップと
を実行すると共に、前記探索範囲演算ステップで、前記特定のスペクトルの周波数が変動したときに、前記探索範囲を、前記特定のスペクトルの変動した周波数に基づいて演算し、変動した前記特定のスペクトルに追従するように実行させること、
を特徴とする回転数の計測方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007030090A JP4865582B2 (ja) | 2007-02-09 | 2007-02-09 | 回転計と回転数の計測方法 |
| EP08002190.0A EP1956497A3 (en) | 2007-02-09 | 2008-02-06 | Signal processor and tracking method of signal search range |
| US12/012,823 US7720646B2 (en) | 2007-02-09 | 2008-02-06 | Signal processor and tracking method of signal search range |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007030090A JP4865582B2 (ja) | 2007-02-09 | 2007-02-09 | 回転計と回転数の計測方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2008196876A JP2008196876A (ja) | 2008-08-28 |
| JP4865582B2 true JP4865582B2 (ja) | 2012-02-01 |
Family
ID=39472488
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007030090A Expired - Fee Related JP4865582B2 (ja) | 2007-02-09 | 2007-02-09 | 回転計と回転数の計測方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7720646B2 (ja) |
| EP (1) | EP1956497A3 (ja) |
| JP (1) | JP4865582B2 (ja) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8229754B1 (en) * | 2006-10-23 | 2012-07-24 | Adobe Systems Incorporated | Selecting features of displayed audio data across time |
| JP5180519B2 (ja) * | 2007-05-29 | 2013-04-10 | 富士重工業株式会社 | 車両用音質評価装置、及び、車両用音質評価方法 |
| JP5349980B2 (ja) * | 2009-01-09 | 2013-11-20 | 株式会社小野測器 | 信号処理装置、信号処理装置の制御方法及び信号処理装置用制御プログラム |
| JP5293300B2 (ja) * | 2009-03-16 | 2013-09-18 | 富士電機株式会社 | 回転機の振動監視装置および振動監視方法 |
| JP5375701B2 (ja) * | 2010-03-23 | 2013-12-25 | 新日鐵住金株式会社 | 回転機械の回転数推定方法、装置及びプログラム |
| JP5692109B2 (ja) * | 2012-02-03 | 2015-04-01 | 新日鐵住金株式会社 | 高炉出銑温度測定方法、高炉出銑温度測定システム、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体 |
| JP2015149060A (ja) * | 2014-01-10 | 2015-08-20 | 株式会社神戸製鋼所 | 周期的な運動を行う機構を備える機器のモデルの作成方法、及びそのモデルを用いた機器の故障診断方法 |
| JP6478462B2 (ja) * | 2014-02-06 | 2019-03-06 | 三菱電機株式会社 | 異常検出装置及び方法 |
| JP6539594B2 (ja) * | 2016-01-18 | 2019-07-03 | 株式会社神戸製鋼所 | 回転機異常検出装置および該方法ならびに回転機 |
| US10208697B2 (en) * | 2016-10-26 | 2019-02-19 | GM Global Technology Operations LLC | Detection of irregularities in engine cylinder firing |
| JP7226050B2 (ja) * | 2018-12-27 | 2023-02-21 | 株式会社明電舎 | 設備診断装置 |
| US11148244B2 (en) * | 2019-08-05 | 2021-10-19 | Computational Systems, Inc. | Vibration spectra window enhancement |
| JP7470072B2 (ja) | 2021-02-25 | 2024-04-17 | 株式会社トランストロン | 異常振動検知装置、異常振動検知方法及びプログラム |
| CN113392813B (zh) * | 2021-07-30 | 2022-04-26 | 湖北工业大学 | 精确识别振动信号主频的方法及系统 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS58131518A (ja) | 1982-02-01 | 1983-08-05 | Toyota Motor Corp | トラッキング分析方法 |
| EP0287852A3 (de) | 1987-04-09 | 1989-02-01 | Siemens Aktiengesellschaft | Verfahren und Anordnung zur Drehzahlerfassung einer elektrischen Maschine |
| JPH0763224B2 (ja) | 1988-03-07 | 1995-07-05 | 株式会社小野測器 | 直流モータの回転測定方法 |
| DE4133269A1 (de) | 1991-10-08 | 1993-04-15 | Bosch Gmbh Robert | Verfahren zur messung der drehzahl eines rotierenden teiles |
| GB2283631B (en) * | 1993-11-06 | 1998-04-29 | Roke Manor Research | Radar apparatus |
| GB9501380D0 (en) * | 1995-01-24 | 1995-03-15 | Sun Electric Uk Ltd | Engine and rotary machine analysis |
| DE19547832C2 (de) * | 1995-12-21 | 2002-04-18 | Grundig Ag | Verfahren und Schaltungsanordnung zur Drehzahlermittlung eines mit einem Generator gekoppelten Verbrennungsmotors |
| DE19702234A1 (de) * | 1997-01-23 | 1998-08-06 | Acida Gmbh Aachener Ct Fuer In | Verfahren zur Überwachung und Qualitätsbeurteilung von sich bewegenden und/oder rotierenden Maschinenteilen insbesondere von Maschinenlagern |
| JP2001159586A (ja) | 1999-12-02 | 2001-06-12 | Suzuki Motor Corp | 音評価方法及び装置並びに音評価用プログラムを記憶した記憶媒体 |
| US7031873B2 (en) * | 2002-06-07 | 2006-04-18 | Exxonmobil Research And Engineering Company | Virtual RPM sensor |
| JP2007030090A (ja) | 2005-07-26 | 2007-02-08 | Fujifilm Corp | 微小構造体の形成方法 |
-
2007
- 2007-02-09 JP JP2007030090A patent/JP4865582B2/ja not_active Expired - Fee Related
-
2008
- 2008-02-06 US US12/012,823 patent/US7720646B2/en not_active Expired - Fee Related
- 2008-02-06 EP EP08002190.0A patent/EP1956497A3/en not_active Withdrawn
Also Published As
| Publication number | Publication date |
|---|---|
| EP1956497A2 (en) | 2008-08-13 |
| US20080195350A1 (en) | 2008-08-14 |
| US7720646B2 (en) | 2010-05-18 |
| JP2008196876A (ja) | 2008-08-28 |
| EP1956497A3 (en) | 2014-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4865582B2 (ja) | 回転計と回転数の計測方法 | |
| RU2449252C2 (ru) | Способ обнаружения повреждения опорного подшипника качения двигателя | |
| Randall | State of the art in monitoring rotating machinery-part 1 | |
| US6681634B2 (en) | Bearing defect detection using time synchronous averaging (TSA) of an enveloped accelerometer signal | |
| Feng et al. | An investigation of the orthogonal outputs from an on-rotor MEMS accelerometer for reciprocating compressor condition monitoring | |
| US9016132B2 (en) | Rotating blade analysis | |
| JP2008544131A (ja) | 外部のマイクロホンによるエンジン状態の検知 | |
| US20180371915A1 (en) | Devices and methods for balancing a high-pressure spool of a gas turbine engine | |
| KR20140084159A (ko) | 기어박스의 상태를 모니터링하는 방법 및 시스템 | |
| CN111397877B (zh) | 一种旋转机械拍振故障检测与诊断方法 | |
| JP2006307669A (ja) | 内燃機関のノッキング判定装置 | |
| RU2011147173A (ru) | Способ вибрационной диагностики и прогнозирования внезапного отказа двигателя и носитель | |
| JP4732304B2 (ja) | 回転計及び回転計用プログラム | |
| US20230003756A1 (en) | Determination of rpm based on scanning vibration spectral plot | |
| CN112345827A (zh) | 频谱频率族的图形区分 | |
| JP5349980B2 (ja) | 信号処理装置、信号処理装置の制御方法及び信号処理装置用制御プログラム | |
| US20120330606A1 (en) | Systems and methods for use in monitoring operation of a rotating device | |
| JP2006300036A (ja) | 内燃機関のノッキング判定装置 | |
| CN117686232A (zh) | 一种燃气轮机振动基频实时提取方法、装置及存储介质 | |
| JP2004053383A (ja) | インターナルポンプの軸ねじり振動測定装置 | |
| JP2008032454A (ja) | 振動位相検出装置及び方法 | |
| JP4071161B2 (ja) | 回転機械の劣化診断方法 | |
| EP0709664A1 (en) | Spectral misfire detection system and method therefor | |
| RU2246711C1 (ru) | Способ определения параметров течения в компрессоре и устройство для его осуществления | |
| RU2558170C2 (ru) | Способ определения частоты вынужденных колебаний рабочего колеса в составе ступени турбомашины |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090127 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090324 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090421 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090528 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20090915 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091116 |
|
| A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20091222 |
|
| A912 | Re-examination (zenchi) completed and case transferred to appeal board |
Free format text: JAPANESE INTERMEDIATE CODE: A912 Effective date: 20100122 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110929 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20111110 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141118 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |