Order-RPM map for order analysis

`map = rpmordermap(x,fs,rpm)`

`map = rpmordermap(x,fs,rpm,res)`

`map = rpmordermap(___,Name,Value)`

```
[map,order,rpm,time,res]
= rpmordermap(___)
```

`rpmordermap(___)`

returns
the order-RPM map matrix, `map`

= rpmordermap(`x`

,`fs`

,`rpm`

)`map`

, that results from
performing order analysis on the input vector, `x`

. `x`

is
measured at a set `rpm`

of rotational speeds expressed
in revolutions per minute. `fs`

is the measurement
sample rate in Hz. Each column of `map`

contains
root-mean-square (RMS) amplitude estimates of the orders present at
each `rpm`

value. `rpmordermap`

resamples `x`

to
a constant samples-per-cycle rate and uses the short-time Fourier
transform to analyze the spectral content of the resampled signal.

specifies
options using `map`

= rpmordermap(___,`Name,Value`

)`Name,Value`

pairs in addition to
the input arguments in previous syntaxes.

`rpmordermap(___)`

with no output
arguments plots the order map as a function of rotational speed and
time on an interactive figure.

Order analysis is the study of vibrations in rotating systems
that result from the rotation itself. The frequencies of these vibrations
are often proportional to the rotational speed. The constants of proportionality
are the *orders*.

The rotational speed is usually measured independently and changes with time under most experimental conditions. Proper analysis of rotation-induced vibrations requires resampling and interpolating the measured signal to achieve a constant number of samples per cycle. Through this process, the signal components whose frequencies are constant multiples of the rotational speed transform into constant tones. The transformation reduces the smearing of spectral components that occurs when frequency changes rapidly with time.

The `rpmordermap`

function performs these
steps:

Uses

`cumtrapz`

to estimate the phase angle as the time integral of the rotational speed:$$\varphi (t)={\displaystyle {\int}_{0}^{t}\frac{\mathrm{RPM}(\tau )}{60}\text{\hspace{0.17em}}}d\tau .$$

Uses

`resample`

to upsample and lowpass-filter the signal. This step enables the function to interpolate the signal at nonsampled time points without aliasing of the high-frequency components.`rpmordermap`

upsamples the signal by a factor of 15.Uses

`interp1`

to interpolate the upsampled signal linearly onto a uniform grid in the phase domain. The highest accessible order in a measurement is fixed by the sample rate and the highest rotational speed reached by the system:$${O}_{\mathrm{max}}=\frac{{f}_{s}/2}{\mathrm{max}\left(\frac{\mathrm{RPM}}{60}\right)}.$$

To capture this highest order accurately, it is necessary to sample the signal at twice

*O*_{max}at least. For better results,`rpmordermap`

oversamples by an extra factor of 4. The resulting phase-domain sample rate,*f*, is_{p}$${f}_{p}=4\times 2{O}_{\mathrm{max}}=4\times 2\frac{{f}_{s}/2}{\mathrm{max}\left(\frac{\mathrm{RPM}}{60}\right)}.$$

The default order resolution,

*r*, is$$r=\frac{{f}_{p}}{256}=\frac{4\times 60}{256}\frac{2\times {f}_{s}/2}{\mathrm{max}(\text{RPM})}=\frac{15}{16}\frac{{f}_{s}}{\mathrm{max}(\text{RPM})}.$$

Uses

`spectrogram`

to compute the short-time Fourier transform (STFT) of the interpolated signal. By default, the function divides the signal into*L*-sample segments and windows each of them with a flat top window. There are$${N}_{\text{overlap}}=\mathrm{min}\left(\lceil \frac{{p}_{\text{overlap}}}{100}\times L\rceil ,L-1\right)$$

samples of overlap between adjoining segments, where

*p*_{overlap}is specified using the`'OverlapPercent'`

name-value pair and defaults to 50%. The DFT length is set to*L*. The resolution is related to the sample rate and segment length through$$r=\frac{k{f}_{p}}{L},$$

where

*k*is the equivalent noise bandwidth of the window, as implemented in`enbw`

. Adjust the resolution to differentiate closely spaced orders. Smaller*r*values require larger segment lengths. If you need to attain a given resolution, make sure that your signal has enough samples.

The red dots in the RPM-vs.-time plot at the bottom of the interactive `rpmordermap`

window
correspond to the right edge of each windowed segment. The blue line
in the plot is the RPM signal drawn as a function of time:

[1] Brandt, Anders. *Noise and Vibration Analysis:
Signal Analysis and Experimental Procedures*. Chichester,
UK: John Wiley & Sons, 2011.

`orderspectrum`

| `ordertrack`

| `orderwaveform`

| `rpmfreqmap`

| `rpmtrack`

| `spectrogram`

| `tachorpm`