Main Content

faultBands

Generate fault frequency bands for spectral feature extraction

Description

example

FB= faultBands(F0,N0)generates fault frequency bandsFB, using the fundamental frequency of interestF0and the array of harmonicsN0. For instance, to construct fault bands for an induction motor, the mains frequency of 60 Hz is the fundamental frequency of interest.

example

FB= faultBands(F0,N0,F1,N1)constructs fault frequency bandsFB, using the distance of the first sidebandF1from the fundamental frequencyF0.N1is the array of the sidebands aroundF0. IfF1is not specified, thenfaultBandssetsF1to 10 percent ofF0by default.N1is equivalent to the 'Sidebands' name-value pair. You can use the 'Type' name-value pair to specify separation between successive sidebands.

example

___= faultBands(___,Name,Value)allows you to specify additional parameters using one or more name-value pair arguments.

example

[FB,info] = faultBands(___)also returns the structureinfocontaining information about the generated fault frequency bandsFB.

example

faultBands(___)with no output arguments plots a bar chart of the generated fault frequency bandsFB.

Examples

collapse all

For this example, generate frequency bands for analyzing the signal components around the first 5 harmonics of the mains supply frequency.

With the fundamental frequency of 60 Hz, the frequency of the alternating current in the mains power supply, usefaultBandsto generate the first 5 harmonics of the mains supply.

F0 = 60; N0 = 1:5; FB = faultBands(F0,N0)
FB =5×258.5000 61.5000 118.5000 121.5000 178.5000 181.5000 238.5000 241.5000 298.5000 301.5000

FBis returned as a 5x2 array with default frequency band width of 5% ofF0which is 3 Hz. The first column inFBcontains the values of F - W 2 , while the second column contains all the values of F + W 2 for each harmonic.

对于这个示例,考虑一个感应电动机broken rotor bars. Under normal operation with load, the rotor speed always lags the speed of the magnetic field allowing the rotor bars to cut magnetic lines of force and produce useful torque. This difference is called slip. Considering a slip value of 0.03 in the system with broken rotors, construct frequency bands for sideband components around the fundamental frequency of 60 Hz.

F0 = 60; N0 = 1:2; slip = 0.03; F1 = 2*slip*F0; N1 = 1:3; [FB,info] = faultBands(F0,N0,F1,N1)
FB =12×247.7000 50.7000 51.3000 54.3000 54.9000 57.9000 62.1000 65.1000 65.7000 68.7000 69.3000 72.3000 107.7000 110.7000 111.3000 114.3000 114.9000 117.9000 122.1000 125.1000 ⋮
info =struct with fields:Centers: [1x12 double] Labels: [1x12 string] HarmonicGroups: [1 1 1 1 1 1 2 2 2 2 2 2]

Construct frequency bands for analyzing the signal components around the first three harmonics of the electrical mains supply frequency.

With the fundamental frequency of 60 Hz, the alternating current in the mains power supply, usefaultBandsto visualize the first 3 harmonics of the mains supply.

F0 = 60; N0 = 1:3; faultBands(F0,N0)

Figure contains an axes. The axes with title Fault Frequency Bands contains 6 objects of type patch, text.

From the plot, observe the following:

  • The fundamental frequency, which is also the first harmonic,1F0at 60 Hz

  • The second harmonic,2F0at 120 Hz

  • The third harmonic,3F0at 180 Hz

To better capture the expected variations of the actual system signals around the nominal fault frequencies, set the widths of each band to 10 Hz.

faultBands(F0,N0,'Width',10)

Figure contains an axes. The axes with title Fault Frequency Bands contains 6 objects of type patch, text.

对于这个示例,考虑一个感应电动机static and dynamic rotor eccentricities. Construct and visualize the frequency bands for the 4 sideband components of an induction motor with 4 pole pairs around the fundamental frequency due to the rotor eccentricities.

F0 = 60; N0 = 1; slip = 0.029; polePairs = 4; F1 = 2*F0*(1-slip)/polePairs
F1 = 29.1300
N1 = 0:4; faultBands(F0,N0,F1,N1)
Warning: Truncated or removed negative fault frequency bands.

Figure contains an axes. The axes with title Fault Frequency Bands contains 14 objects of type patch, text.

To avoid truncating negative fault frequency bands, set 'Folding' totrueto fold them onto the positive frequency axis.

faultBands(F0,N0,F1,N1,'Folding',真正的)

Figure contains an axes. The axes with title Fault Frequency Bands contains 18 objects of type patch, text.

Observe that the sideband frequencies1F0-3F1and1F0-4F1are now visible on the positive axis.

Input Arguments

collapse all

基频of interest, specified as a positive scalar.faultBandsconstructs the fault frequency bands around the fundamental frequencyF0. For instance, to construct fault bands for a faulty induction motor, the mains frequency of 60 Hz is the fundamental frequency of interest. Similarly, to generate fault bands for a faulty gear train, the input shaft frequency is the fundamental frequency.

You can specifyF0in either hertz or orders.

Harmonics of the fundamental frequency, specified as a vector of positive integers. Specify fault bands around the fundamental frequencyF0and its harmonics byN0.N0is equivalent to the 'Harmonics' name-value pair with a default value of 1.

Distance of the first sideband from the fundamental frequency, specified as a positive scalar. IfF1is not specified, thenfaultBandsassumes a value of 10 percent of the fundamental frequency forF1.

Sidebands of the fundamental frequency and its harmonics, specified as a vector of nonnegative integers.N1is equivalent to the 'Sidebands' name-value pair with a default value of 0.

Name-Value Pair Arguments

Specify optional comma-separated pairs ofName,Valuearguments.Nameis the argument name andValueis the corresponding value.Namemust appear inside quotes. You can specify several name and value pair arguments in any order asName1,Value1,...,NameN,ValueN.

Example:...,'Harmonics',[1,3,5]

Harmonics of the fundamental frequency to be included, specified as the comma-separated pair consisting of 'Harmonics' and a vector of positive integers. The default value is 1. Specify 'Harmonics' when you want to construct the frequency bands with more harmonics of the fundamental frequency.

Sidebands around the fundamental frequency and its harmonics to be included, specified as the comma-separated pair consisting of 'Sidebands' and a vector of nonnegative integers. The default value is 0. Specify 'Sidebands' when you want to construct the frequency bands with sidebands around the fundamental frequency and its harmonics.

Width of the frequency bands centered at the nominal fault frequencies, specified as the comma-separated pair consisting of 'Width' and a positive scalar. The default value is 5 percent of the fundamental frequency. Avoid specifying 'Width' with a large value so that the fault bands do not overlap.

Separation value between successive sidebands, specified as the comma-separated pair consisting of 'Type' and either'additive'or'multiplicative'. Specify 'Type' as:

  • 'additive', to set the separation between successive sidebands toF1.

  • 'multiplicative', to set the separation between successive sidebands proportional to both the harmonic order and the sideband value.

Logical value specifying whether negative nominal fault frequencies have to be folded about the frequency origin, specified as the comma-separated pair consisting of 'Folding' and eithertrueorfalse. If you set 'Folding' totrue, thenfaultBandsfolds the negative nominal fault frequencies about the frequency origin by taking their absolute values such that the folded fault bands always fall in the positive frequency intervals. The folded fault bands are computed as [ max ( 0 , | F | W 2 ) , | F | + W 2 ] , whereWis the 'Width' name-value pair andFis one of the nominal fault frequencies.

Output Arguments

collapse all

Fault frequency bands, returned as anNx2array, whereNis the number of fault frequencies.FBis returned in the same units asF0, in either Hertz or orders. The generated fault bands, [ F W 2 , F + W 2 ] , are centered depending on the sideband specification as follows:

  • If you do not specify the sidebands, then the fault bands are centered at F = n 0 F 0 , where the integern0ranges through the elements of the array of harmonics,N0.

  • If you specify sidebands usingN1or the 'Sidebands' name-value pair, then fault bands are centered at:

    • F = n 0 F 0 ± n 1 F 1 , when 'Type' is specified as'additive'. Here, the integern1ranges through the elements of the array of sidebands,N1.

    • F = n 0 ( F 0 ± n 1 F 1 ) , when 'Type' is specified as'multiplicative'.

关于故障频带的信息inFB, returned as a structure with the following fields:

  • Centers— Center fault frequencies

  • Labels— Labels describing each frequency

  • HarmonicGroups— Harmonic group numbers equal to the harmonic order of each frequency band to be able to identify fault bands associated with the nominal fault frequency F = n 0 F 0 , where the integern0ranges through the elements of the array of harmonics,N0

Introduced in R2019b