JP6550476B2 - 信号を分析するための方法およびその方法を実行するための装置 - Google Patents
信号を分析するための方法およびその方法を実行するための装置 Download PDFInfo
- Publication number
- JP6550476B2 JP6550476B2 JP2017560963A JP2017560963A JP6550476B2 JP 6550476 B2 JP6550476 B2 JP 6550476B2 JP 2017560963 A JP2017560963 A JP 2017560963A JP 2017560963 A JP2017560963 A JP 2017560963A JP 6550476 B2 JP6550476 B2 JP 6550476B2
- Authority
- JP
- Japan
- Prior art keywords
- amplitude
- spectrum
- highest
- point
- frequency
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 77
- 230000003595 spectral effect Effects 0.000 claims description 237
- 238000001228 spectrum Methods 0.000 claims description 170
- 238000012937 correction Methods 0.000 claims description 110
- 238000011156 evaluation Methods 0.000 claims description 19
- 230000036962 time dependent Effects 0.000 description 8
- 238000000605 extraction Methods 0.000 description 7
- 238000002485 combustion reaction Methods 0.000 description 6
- 230000001419 dependent effect Effects 0.000 description 6
- 238000009434 installation Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000001143 conditioned effect Effects 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
- G01R23/167—Spectrum analysis; Fourier analysis using filters with digital filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L19/00—Details of, or accessories for, apparatus for measuring steady or quasi-steady pressure of a fluent medium insofar as such details or accessories are not special to particular types of pressure gauges
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
-
- 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
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/12—Analogue/digital converters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/12—Analogue/digital converters
- H03M1/124—Sampling or signal conditioning arrangements specially adapted for A/D converters
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Engines (AREA)
Description
− 所定の窓関数による分析される信号の乗算、およびその後の、特に離散または高速フーリエ変換の結果である、分析される信号のスペクトルが提供され、
− 最も高いスペクトル線およびそれに直接隣接する2番目に高いスペクトル線が、少なくとも1つの極大値の領域内のスペクトル内で決定され、
− 最も高いスペクトル線の振幅およびそれに直接隣接する2番目に高いスペクトル線の振幅が決定され、
− 最も高いスペクトル線の振幅および2番目に高いスペクトル線の振幅に依存する少なくとも1つの値、特に、最も高いスペクトル線の振幅の2番目に高いスペクトル線の振幅に対する比、ならびに/または最も高いスペクトル線の振幅の最も高いスペクトル線および2番目に高いスペクトル線の振幅の平均値に対する比が決定され、
− 最も高いスペクトル線の周波数に対する周波数補正値が、所定の窓関数のために作成された周波数補正特性曲線を使用して、決定された値、特に比から決定され、
− 最も高いスペクトル線の周波数が、周波数補正値によって補正される。
− 2つ以上の、特に正弦波の、参照信号が提供され、
− 各々の参照信号が窓関数によって乗算され、
− 窓関数によって乗算された各々の参照信号が、特に離散または高速フーリエ変換にかけられ、
− 各々の場合において、最も高いスペクトル線およびそれに直接隣接する2番目に高いスペクトル線が、極大値の領域内の各々の変換されたスペクトルから決定され、
− 最も高いスペクトル線の振幅および2番目に高いスペクトル線の振幅が、各々の変換されたスペクトルに対して決定され、
− 最も高いスペクトル線の振幅および2番目に高いスペクトル線の振幅に依存する少なくとも1つの値、特に、最も高いスペクトル線の振幅の2番目に高いスペクトル線の振幅に対する比、ならびに/または最も高いスペクトル線の振幅の最も高いスペクトル線および2番目に高いスペクトル線の振幅の平均値に対する比が、各々の変換されたスペクトルに対して決定され、
− 決定された値、特に比が、参照信号の関連周波数と一緒に値のペアとして格納され、
− 特性曲線が値のペアを基に作成される、という手順によって作成される。
− 特に正弦波の、既知の周波数の参照信号が提供され、
− fn=flineL+n/N*Δf、ここで、n=0…Nで分解能Δfの周波数を有するN個のさらなる参照信号を得るために、参照信号が、任意のスペクトル線の周波数flineLとそれに隣接するスペクトル線の周波数との間で変動する、という手順によって提供される複数の参照信号が実現され得る。
本発明に従う方法を使用して条件付けされる周波数値は、例えば、設備、特にガスタービンの状態に関する情報のため、メッセージ収集のため、または設備の動作に対する影響のために使用され得る。
− 各々の場合において、最も高いスペクトル線およびそれに直接隣接する2番目に高いスペクトル線が、複数の極大値の領域内の、特にすべての極大値の領域内のスペクトル内で決定され、
− 最も高いスペクトル線の振幅および2番目に高いスペクトル線の振幅が各々の極大値に対して決定され、
− 最も高いスペクトル線の振幅および2番目に高いスペクトル線の振幅に依存する少なくとも1つの値、特に、最も高いスペクトル線の振幅の2番目に高いスペクトル線の振幅に対する比、ならびに/または最も高いスペクトル線の振幅の、最も高いスペクトル線および2番目に高いスペクトル線の振幅の平均値に対する比が、各々の極大値に対して算出され、
− 最も高いスペクトル線の周波数に対する周波数補正値が、各々の極大値に対する周波数補正特性曲線を使用して、算出された値、特に比から決定され、
− 最も高いスペクトル線の周波数が、各々の極大値に対するそれぞれの周波数補正値によって補正される、ということによって特徴付けられる。
− 分析される信号を記録するための少なくとも1つのセンサ、特に、少なくとも1つの圧力センサおよび/または少なくとも1つの振動センサに接続されるように設計され、
− 少なくとも1つのセンサによって記録された信号を獲得するように、および特に、それをフィルタリングする、および/またはデジタル化するように設計され、
− 信号を所定の窓関数によって乗算するように設計され、
− 信号のスペクトルを得るために、信号を、特に高速または離散フーリエ変換にかけるように設計される、
データ獲得および評価デバイスを備え、
データ獲得および評価デバイスは、少なくとも1つのセンサによって獲得される信号を分析するための本発明に従う方法を実行するように設計および構成される。
A1/Avg12(n/N)=A1/Avg12(−n/N)
A1/Avg12(n/N)=A1/Avg12(n/N+1)
f(A1,A2,A3,n/N)=f(A1,A2,A3,−n/N)
f(A1,A2,A3,n/N)=f(A1,A2,A3,n/N+1)
A1/Avg12≒(Max(A1/Avg12)−Min(A1/Avg12))/0.5*n/N+Max(A1/Avg12)
および範囲n/N=0〜n/N=0.5については、以下が当てはまる。
A1/Avg12≒(Max(A1/Avg12)−Min(A1/Avg12))/0.5*n/N+Max(A1/Avg12)
F1:y=0.8139x+1.4056が得られ、ここで、x=n/Nであり、
n/N=0〜n/N=0.5の範囲にわたって第2の線形フィット関数、
F2:y=−0.7971x+1.4022が得られ、ここで、x=n/Nである。
F3:y=0.925x2−0.0328x+1.0008、ここで、x=n/Nであり、
さらなる二次多項式フィット関数は、n/N=0.5〜n/N=1の範囲について得られ、
F4:y=0.9666x2−1.8925x+1.9265であり、ここで、x=n/Nである。
n/N=0でMin(A1setpoint/A1)=1.0
n/N=0.5でMax(A1setpoint/A1)=1.22
F5:y=0.7081x2−0.7083x+1.3979、ここで、x=n/Nである。
n/N=−0.5*(A1/Avg12−Max(A1/Avg12)/((Max(A1/Avg12)−Min(A1/Avg12))
および、範囲n/N=0.5〜n/N=1で範囲n/N=−0.5〜n/N=0での関数の周期性を考慮して、すなわち“MaxIsRight”=TRUEの場合は、
n/N=0.5*(A1Avg12−Max(A1/Avg12)/((Max(A1/Avg12)−Min(A1/Avg12))
へと転置される。
F6: y=−1.2496x+1.7527、ここで、x=A1/Avg12であり、
以下が範囲n/N=−0.5〜n/N=0で得られ、
F7:y=1.2496x−1.7527、ここで、x=A1/Avg12である。
fncorr=flineL+n/N*Δf
と推定され、信号成分の実際の周波数によく一致する、最も高いスペクトル線S1の補正された周波数、
fncorr=10.26Hz
が得られる。
A1corrMax=Avg12measured*F5
ここで、
F5:y=0.7081x2−0.7073x+1.3979であり、ここで、x=n/Nである。
A1corrMax=2.98、
それは2.97という実際の振幅によく一致する。
A1corrMax>A1
を有する値を生み出す。
ΔA1corrRed/ΔA1corrMax=Actual2−3/Setpoint2−3、および
(A1corrRed−A1)/(A1corrMax−A1)=Actual2−3/Setpoint2−3
A1corrRed=A1+(A1corrMax−A1)*(Actual2−3/Setpoint2−3)である。
Actual2−3=0.7、および
Setpoint2−3=1.14である。
精度=ΔA1corrRed/ΔA1corrMax
精度=0.035/0.06=62%
A2 振幅、2番目に高いスペクトル線S2の振幅
A3 振幅、3番目に高いスペクトル線S3の振幅
Avg12 振幅A1、A2の平均値
f0、f30、f50、f80、f100 周波数
P1、P2、P3 ピーク
S1 スペクトル線、最も高いスペクトル線
S2 スペクトル線、2番目に高いスペクトル線
S3 スペクトル線、3番目に高いスペクトル線
Claims (20)
- 信号を分析するための方法であって、
2つ以上の、正弦波の、異なる周波数の参照信号を使用して、所定の窓関数のために周波数補正特性曲線が作成され、
前記周波数補正特性曲線が、
− 2つ以上の、正弦波の、異なる周波数の参照信号が提供され、
− 各々の参照信号が前記窓関数によって乗算され、
− 前記窓関数によって乗算された各々の参照信号が、離散または高速フーリエ変換にかけられ、
− 複数の前記参照信号の各々の場合において、振幅が最も高いスペクトル上の点(S 1 )およびそれに直接隣接する前記振幅が2番目に高いスペクトル上の点(S 2 )が、振幅の極大値を含む領域内の各々の変換されたスペクトルから決定され、
− 前記振幅が最も高いスペクトル上の点(S 1 )の前記振幅(A 1 )および前記振幅が2番目に高いスペクトル上の点(S 2 )の前記振幅(A 2 )が、各々の変換されたスペクトルに対して決定され、
− 前記振幅が最も高いスペクトル上の点(S 1 )の前記振幅(A 1 )の前記振幅が2番目に高いスペクトル上の点(S 2 )の前記振幅(A 2 )に対する比、ならびに/または前記振幅が最も高いスペクトル上の点(S 1 )の前記振幅(A 1 )の、前記振幅が最も高いスペクトル上の点および前記振幅が2番目に高いスペクトル上の点(S 1 、S 2 )の前記振幅(A 1 、A 2 )の平均値(Avg 12 )に対する比が、各々の変換されたスペクトルに対して決定され、
− 前記決定された比が、前記比を決定するために使用された前記参照信号の周波数と一緒に値のペアとして格納され、
− 特性曲線が前記値のペアを基に作成される、という手順によって、前記所定の窓関数のために作成され、
前記所定の窓関数による、分析される前記信号の乗算、およびその後の、離散または高速フーリエ変換の結果である、分析される前記信号のスペクトルが提供され、
− 振幅が最も高いスペクトル上の点(S1)およびそれに直接隣接する振幅が2番目に高いスペクトル上の点(S2)が、少なくとも1つの振幅の極大値を含む領域内の前記スペクトル内で決定され、
− 前記振幅が最も高いスペクトル上の点の振幅(A1)およびそれに直接隣接する振幅が2番目に高いスペクトル上の点の振幅(A2)が決定され、
− 前記振幅が最も高いスペクトル上の点の前記振幅(A1)の前記振幅が2番目に高いスペクトル上の点の前記振幅(A2)に対する比、ならびに/または前記振幅が最も高いスペクトル上の点の前記振幅(A1)の、前記振幅が最も高いスペクトル上の点および振幅が2番目に高いスペクトル上の点の振幅(A1、A2)の平均値(Avg12)に対する比、が決定され、
− 前記振幅が最も高いスペクトル上の点(S1)の周波数に対する周波数補正値が、前記所定の窓関数のために作成された周波数補正特性曲線を使用して、前記決定された比から決定され、
− 前記振幅が最も高いスペクトル上の点(S1)の前記周波数が、前記周波数補正値によって補正されることを特徴とする、方法。 - 前記分析される信号が複数の周波数及び振幅の信号成分を含む場合において、
− 前記スペクトルに現れる極大値の各々において、前記振幅が最も高いスペクトル上の点(S1)およびそれに直接隣接する前記振幅が2番目に高いスペクトル上の点(S2)が、複数の振幅の極大値を含む領域内のスペクトル内で決定され、
− 前記振幅が最も高いスペクトル上の点(S1)の前記振幅(A1)および前記振幅が2番目に高いスペクトル上の点(S2)の前記振幅(A2)が、前記極大値の各々に対して決定され、
− 前記振幅が最も高いスペクトル上の点(S1)の前記振幅(A1)の前記振幅が2番目に高いスペクトル上の点(S2)の前記振幅(A2)に対する比、ならびに/または前記振幅が最も高いスペクトル上の点(S1)の前記振幅(A1)の、前記振幅が最も高いスペクトル上の点および前記振幅が2番目に高いスペクトル上の点(S1、S2)の前記振幅(A1、A2)の平均値(Avg12)に対する比が、前記極大値の各々に対して算出され、
− 前記振幅が最も高いスペクトル上の点(S1)の周波数に対する周波数補正値が、前記極大値の各々に対する前記周波数補正特性曲線を使用して、前記算出された比から決定され、
− 前記振幅が最も高いスペクトル上の点(S1)の前記周波数が、前記極大値の各々に対するそれぞれの周波数補正値によって補正される、ことを特徴とする、請求項1に記載の方法。 - 複数の前記参照信号が、
− 正弦波の、既知の周波数の参照信号が提供され、
− fn=flineL+n/N*Δf、
ここで、n=0…Nで分解能Δfの周波数を有するN個のさらなる参照信号を得るために、前記参照信号が、任意のスペクトル上の点の周波数(flineL)とそれに隣接するスペクトル上の点の周波数(flineL+1)との間で変動する、
という手順によって提供されることを特徴とする、請求項1に記載の方法。 - 前記参照信号の前記スペクトルから決定される比であって、前記振幅が最も高いスペクトル上の点(S1)の前記振幅(A1)の前記振幅が2番目に高いスペクトル上の点(S2)の前記振幅(A2)に対する比、ならびに/または前記振幅が最も高いスペクトル上の点(S1)の前記振幅(A1)の、前記振幅が最も高いスペクトル上の点および前記振幅が2番目に高いスペクトル上の点(S1、S2)の前記振幅(A1、A2)の平均値(Avg12)に対する比が、変動した周波数(fn)に対する、または比n/Nに対する、または前記周波数のオフセットn/N*Δfに対する関数として求められることを特徴とする、請求項3に記載の方法。
- 少なくとも1つのフィット関数(F1、F2)が、前記参照信号の前記スペクトルから決定される前記比であって、前記振幅が最も高いスペクトル上の点(S 1 )の前記振幅(A 1 )の前記振幅が2番目に高いスペクトル上の点(S 2 )の前記振幅(A 2 )に対する比、ならびに/または前記振幅が最も高いスペクトル上の点(S 1 )の前記振幅(A 1 )の、前記振幅が最も高いスペクトル上の点および前記振幅が2番目に高いスペクトル上の点(S 1 、S 2 )の前記振幅(A 1 、A 2 )の平均値(Avg 12 )に対する前記比の値に適用され、前記少なくとも1つのフィット関数(F1、F2)が、周波数補正特性曲線として使用される、ことを特徴とする、請求項4に記載の方法。
- 振幅補正値および/または補正された振幅が、異なる周波数の2つ以上の参照信号を使用して作成された振幅補正特性曲線を使用して、分析される前記信号からの前記少なくとも1つの振幅の極大値を含む領域内の前記振幅が最も高いスペクトル上の点(S1)に対する前記周波数補正値または前記補正された周波数から決定される、ことを特徴とする、請求項1から5のいずれか一項に記載の方法。
- 振幅補正値および/または補正された振幅が、異なる周波数の2つ以上の参照信号を使用して作成された振幅補正特性曲線を使用して、分析される前記信号からの前記少なくとも1つの振幅の極大値を含む領域内の前記振幅が最も高いスペクトル上の点(S1)に対する前記周波数補正値または前記補正された周波数から決定され、
前記振幅が最も高いスペクトル上の点(S1)の前記振幅の所定の設定値(A1setpoint)の、前記振幅が最も高いスペクトル上の点および振幅が2番目に高いスペクトル上の点(S1、S2)の前記振幅(A1、A2)の平均値(Avg12)に対する比が、各々の参照信号のスペクトルから決定され、前記決定された比が、前記比を決定するために使用された前記参照信号の周波数と一緒に値のペアとして格納され、特性曲線が前記値のペアを基に作成される、という手順によって、前記振幅補正特性曲線が作成されることを特徴とする、請求項3から5のいずれか一項に記載の方法。 - 前記参照信号の前記スペクトルから決定された前記比が、変動した周波数(fn)に対する、前記比n/Nに対する、または前記周波数のオフセットn/N*Δfに対する関数として求められ、フィット関数(F5)が、前記参照信号の前記スペクトルから決定された前記比の値に適用され、前記フィット関数が、振幅補正関数として使用される、ことを特徴とする、請求項7に記載の方法。
- 前記振幅が最も高いスペクトル上の点(S1)に直接隣接する振幅が3番目に高いスペクトル上の点(S3)および前記振幅が3番目に高いスペクトル上の点の振幅(A3)が、各々の場合において前記振幅が最も高いスペクトル上の点および前記振幅が2番目に高いスペクトル上の点(S1、S2)が決定されたその振幅の極大値を含む領域内の前記参照信号の前記スペクトルから決定され、前記振幅が2番目に高いスペクトル上の点(S2)の前記振幅(A2)と前記振幅が3番目に高いスペクトル上の点(S3)の前記振幅(A3)との間の差が決定されることを特徴とする、請求項7または8に記載の方法。
- 前記参照信号の前記スペクトルから決定された前記差が、変動した周波数(fn)に対する、または前記比n/Nに対する、または前記周波数のオフセットn/N*fにに対する関数として求められ、少なくとも1つのフィット関数(F8)が、前記参照信号の前記スペクトルから決定された前記差の値に適用され、前記フィット関数(F8)が、さらなる振幅補正特性曲線として使用される、ことを特徴とする、請求項9に記載の方法。
- 前記振幅が最も高いスペクトル上の点(S1)に直接隣接する前記振幅が3番目に高いスペクトル上の点(S3)および前記振幅が3番目に高いスペクトル上の点の前記振幅(A3)が、前記振幅が最も高いスペクトル上の点および前記振幅が2番目に高いスペクトル上の点(S1、S2)が決定されたその振幅の極大値を含む領域内の分析される前記信号の前記スペクトル内で決定され、前記振幅が2番目に高いスペクトル上の点(S2)の前記振幅(A2)と前記振幅が3番目に高いスペクトル上の点(S3)の前記振幅(A3)との間の差が、実際の差の値(Actual2−3)として決定されることを特徴とする、請求項10に記載の方法。
- 前記振幅が2番目に高いスペクトル上の点(S2)の前記振幅(A2)と前記振幅が3番目に高いスペクトル上の点(S3)の前記振幅(A3)との間の前記差の設定差の値(SetpointA2−A3)が、前記フィット関数(F8)から決定され、前記実際の差の値(ActualA2−A3)と比較され、
前記設定差の値が、前記振幅が最も高いスペクトル上の点及び前記振幅が2番目に高いスペクトル上の点(S1、S2)の振幅(A1、A2)の平均値(Avg12)に前記フィット関数(F5)を乗じて求めた補正された振幅に、前記比n/Nに窓関数ごとに異なる所定の定数を乗じた値の絶対値を乗じて得られることを特徴とする、請求項11に記載の方法。 - 前記振幅補正値および/または前記補正された振幅が、前記実際の差の値(ActualA2−A3)が前記設定差の値(SetpointA2−A3)から外れるところまで減少させられることを特徴とする、請求項12に記載の方法。
- 前記減少させられた振幅補正値の前記振幅補正値に対する比および/または前記減少させられた補正された振幅の前記補正された振幅に対する比が形成され、精度値として決定されることを特徴とする、請求項13に記載の方法。
- 前記振幅が最も高いスペクトル上の点(S1)の前記振幅の所定の設定値(A1setpoint)の前記振幅が最も高いスペクトル上の点(S1)の前記振幅のそれぞれの実際の値(A1actual)に対する比が、前記参照信号の前記スペクトルから決定され、前記決定された比が、変動した周波数(fn)に対する、または前記比n/Nに対する、または前記周波数のオフセットn/N*fにに対する関数として求められることを特徴とする、請求項3から14のいずれか一項に記載の方法。
- 前記振幅が最も高いスペクトル上の点(S1)の前記振幅の前記所定の設定値(A1setpoint)と前記振幅が最も高いスペクトル上の点(S1)の前記振幅の前記それぞれの実際の値(A1actual)の比から、前記比の最小値および/または最大値が決定され、前記最大値が、前記減少させられた振幅補正値の妥当性のチェックに使用されることを特徴とする、請求項15に記載の方法。
- 前記参照信号が、N=10ステップで変動することを特徴とする、請求項1から16のいずれか一項に記載の方法。
- ハミング窓関数またはテューキー窓関数またはハニング窓関数または矩形窓関数が、前記窓関数として使用されることを特徴とする、請求項1から17のいずれか一項に記載の方法。
- − 分析される信号を記録するための少なくとも1つのセンサに接続されるように設計され、
− 前記少なくとも1つのセンサによって記録された信号を獲得するように設計され、
− 前記信号を所定の窓関数によって乗算するように設計され、
− 前記信号のスペクトルを得るために、前記信号を、高速または離散フーリエ変換にかけるように設計される、
データ獲得および評価デバイスを備え、
前記データ獲得および評価デバイスが、前記少なくとも1つのセンサによって獲得された信号を分析するための請求項1から18のいずれか一項に記載の方法を実行するように設計および構成される、請求項1から18のいずれか一項に記載の信号を分析するための方法を実行するための装置。 - 前記データ獲得および評価デバイスに接続される、分析される信号を記録するための少なくとも1つのセンサを備える、請求項19に記載の装置。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP15173342.5A EP3109646A1 (de) | 2015-06-23 | 2015-06-23 | Verfahren zur analyse eines signals sowie vorrichtung zur durchführung des verfahrens |
EP15173342.5 | 2015-06-23 | ||
PCT/EP2016/061402 WO2016206876A1 (de) | 2015-06-23 | 2016-05-20 | Verfahren zur analyse eines signals sowie vorrichtung zur durchführung des verfahrens |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018523101A JP2018523101A (ja) | 2018-08-16 |
JP6550476B2 true JP6550476B2 (ja) | 2019-07-24 |
Family
ID=53491293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017560963A Active JP6550476B2 (ja) | 2015-06-23 | 2016-05-20 | 信号を分析するための方法およびその方法を実行するための装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10234491B2 (ja) |
EP (2) | EP3109646A1 (ja) |
JP (1) | JP6550476B2 (ja) |
CN (1) | CN107787455B (ja) |
WO (1) | WO2016206876A1 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110166813B (zh) * | 2019-05-22 | 2021-06-08 | 南京巨鲨显示科技有限公司 | 一种显示器多窗口图像独立校正方法及系统 |
FR3097960B1 (fr) * | 2019-06-27 | 2021-06-18 | Safran Aircraft Engines | Procédé d’acquisition de données pour la détection d’un endommagement d’un palier |
CN113009432B (zh) * | 2020-02-28 | 2024-05-31 | 加特兰微电子科技(上海)有限公司 | 提高测量精度、提高目标探测精度的方法、装置及设备 |
CN113484544B (zh) * | 2021-08-25 | 2024-01-12 | 浙江上风高科专风实业股份有限公司 | 一种面向无线振动速度传感器的中低频频谱矫正方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0619390B2 (ja) * | 1985-08-08 | 1994-03-16 | 横河・ヒユ−レツト・パツカ−ド株式会社 | デイジタル・フ−リエ変換の後処理方法 |
WO2003102610A2 (de) * | 2002-05-31 | 2003-12-11 | Technische Universität München | Verfahren und vorrichtung zur erfassung im frequenzbereich basierend auf einer zeitbereichsmessung |
JP5035815B2 (ja) * | 2004-07-05 | 2012-09-26 | 学校法人中部大学 | 周波数測定装置 |
DE102006022056A1 (de) * | 2006-02-20 | 2007-08-30 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Spektralanalyse zur zuverlässigeren Bestimmung von Vitalparametern |
EP1843230A1 (de) * | 2006-04-07 | 2007-10-10 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zum Erfassen von Werten eines Signals |
JP4260831B2 (ja) * | 2006-09-20 | 2009-04-30 | 三菱電機株式会社 | 車載用周波数変調レーダ装置 |
US7751277B2 (en) * | 2008-03-17 | 2010-07-06 | Pgs Geophysical As | Method for interpolating seismic data by anti-alias, anti-leakage Fourier transform |
EP2799838B1 (en) * | 2011-12-27 | 2016-07-27 | Canon Kabushiki Kaisha | Information signal generating method |
WO2014203708A1 (ja) * | 2013-06-17 | 2014-12-24 | アルプス電気株式会社 | 信号周波数算出方法 |
-
2015
- 2015-06-23 EP EP15173342.5A patent/EP3109646A1/de not_active Withdrawn
-
2016
- 2016-05-20 CN CN201680037016.7A patent/CN107787455B/zh active Active
- 2016-05-20 JP JP2017560963A patent/JP6550476B2/ja active Active
- 2016-05-20 EP EP16727330.9A patent/EP3281021B1/de active Active
- 2016-05-20 US US15/578,901 patent/US10234491B2/en active Active
- 2016-05-20 WO PCT/EP2016/061402 patent/WO2016206876A1/de active Application Filing
Also Published As
Publication number | Publication date |
---|---|
EP3109646A1 (de) | 2016-12-28 |
EP3281021A1 (de) | 2018-02-14 |
CN107787455B (zh) | 2019-07-19 |
WO2016206876A1 (de) | 2016-12-29 |
EP3281021B1 (de) | 2021-03-17 |
US10234491B2 (en) | 2019-03-19 |
JP2018523101A (ja) | 2018-08-16 |
CN107787455A (zh) | 2018-03-09 |
US20180172743A1 (en) | 2018-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6550476B2 (ja) | 信号を分析するための方法およびその方法を実行するための装置 | |
CN1796756B (zh) | 缸内压的检测 | |
RU2514061C2 (ru) | Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток | |
JP5422401B2 (ja) | レゾルバ信号変換装置及び方法 | |
JP4965613B2 (ja) | レゾルバインタフェース、レゾルバからの入力信号を監視する方法、および差動位置センサインタフェース | |
US9952065B2 (en) | Position sensor device to determine a position of a moving device | |
CN112577746B (zh) | 一种转速波动下滚动轴承包络阶次谱故障特征的提取方法 | |
JP4476321B2 (ja) | Rdコンバータ及び角度検出装置 | |
CN108956117B (zh) | 电磁振动分量的去除方法、旋转机械诊断方法及装置 | |
JP6506681B2 (ja) | 空気流量測定装置 | |
EA022630B1 (ru) | Устройство и способ для анализа вибраций машины, имеющей вращающуюся часть | |
CN112105907B (zh) | 用于监测齿轮系统的方法和设备 | |
CN111835262B (zh) | 半导体装置、电机控制系统及误差检测方法 | |
CA2955531A1 (en) | Vibration monitoring systems | |
JP7391341B2 (ja) | 電気角取得システム、電気角取得方法および電気角取得プログラム | |
US20210140928A1 (en) | Rubbing detection device for rotary machine and method for detecting rubbing of rotary machine | |
US20040153780A1 (en) | Method and device for the correction of the dynamic error of a sensor | |
JP5773932B2 (ja) | 電動モータの回転情報検出方法、電動モータの回転情報検出装置、電動モータ制御装置 | |
JP5325266B2 (ja) | 回転体のアンバランス量測定方法および装置 | |
JP4951060B2 (ja) | 渦流量計からの測定信号を処理するための方法 | |
US8412478B2 (en) | Device for determining an error induced by a high-pass filter and associated error correction method | |
CN110537082B (zh) | 振动检测装置及异常判定系统 | |
JP6760612B2 (ja) | 位相計測装置およびこの位相計測装置を適用した機器 | |
JP2005214932A5 (ja) | ||
JP2005214932A (ja) | 信号処理装置、この信号処理装置を用いた電圧測定装置及び電流測定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20180618 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180702 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20181002 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190104 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190307 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20190603 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190701 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6550476 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
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 |