Moving RMS

Moving RMS

• 

Libraries:
DSP System Toolbox / Statistics

Description

The Moving RMS block computes the moving root mean square (RMS) of the input signal along each channel independently over time. The block uses either the sliding window method or the exponential weighting method to compute the moving RMS. In the sliding window method, a window of specified length moves over the data sample by sample, and the block computes the RMS over the data in the window. In the exponential weighting method, the block squares the data samples, multiplies them with a set of weighting factors, and sums the weighed data. The block then computes the RMS by taking the square root of the sum. For more details on these methods, see Algorithms.

Examples

Compute Moving RMS of Noisy Step Signal

Compute moving RMS using both the sliding window method and the exponential weighting method.

Open Model

Ports

Input

x — Data input

vector | matrix

The block computes the moving RMS of the data specified at this input port. Specify real- or complex-valued multichannel inputs of the size m-by-n, where m ≥ 1 and n ≥ 1.

When the Allow arbitrary frame length for fixed-size input signals parameter appears and is not selected, and you input a fixed-size signal, the frame length must be a multiple of the hop size (window length − overlap length). In all other cases, the input frame length can be arbitrary.

The block accepts variable-size inputs (frame length changes during simulation). When you input a variable-size signal, the frame length of the signal can be arbitrary.

This port is unnamed until you set Method to Exponential weighting and select the Specify forgetting factor from input port parameter.

Data Types: single | double
Complex Number Support: Yes

lambda — Forgetting factor

positive real scalar in the range (0,1]

The forgetting factor determines how much weight past data is given. A forgetting factor of 0.9 gives more weight to the older data than does a forgetting factor of 0.1. A forgetting factor of 1.0 indicates infinite memory – all previous samples are given an equal weight.

Dependencies

This port appears when you set Method to Exponential weighting and select the Specify forgetting factor from input port parameter.

Data Types: single | double

Output

Port_1 — Moving RMS output

vector | matrix

Moving RMS output, returned as a vector or a matrix. The block computes the moving RMS based on the Method parameter settings using either the sliding window method or the exponential weighting method. For more details, see Algorithms.

This table provides more details on the dimensions of the output signal.

Input Signal	Input Dimensions	Output Dimensions When Allow arbitrary frame length for fixed-size input signals Appears	Output Dimensions When Allow arbitrary frame length for fixed-size input signals Does Not Appear
Fixed-size signal	m-by-n, where m is a multiple of the hop size (window length − overlap length)	(m/hop size)-by-n	m-by-n
Fixed-size signal	m-by-n, where m is not a multiple of the hop size (window length − overlap length)	ceil(m/hop size)-by-n when you select Allow arbitrary frame length for fixed-size input signals If you do not select Allow arbitrary frame length for fixed-size input signals, the block errors.	m-by-n
Variable-size signal	m-by-n	ceil(m/hop size)-by-n	m-by-n

When the output has an upper bound size of ceil(m/hop size)-by-n, during simulation, the size of the first dimension varies within this bound and the size of the second dimension remains constant. For an example that shows this behavior, see Compute Moving RMS of Noisy Step Signal.

Data Types: single | double
Complex Number Support: Yes

Parameters

If a parameter is listed as tunable, then you can change its value during simulation.

Method — Moving RMS method

Sliding window (default) | Exponential weighting

• Sliding window — A window of length Window length moves over the input data along each channel. For every sample the window moves over, the block computes the RMS over the data in the window.

• Exponential weighting — The block multiplies the squares of the samples by a set of weighting factors. The magnitude of the weighting factors decreases exponentially as the age of the data increases, but the magnitude never reaches zero. To compute the RMS, the algorithm sums the weighted data and takes a square root of the sum.

Specify window length — Flag to specify window length

on (default) | off

When you select this check box, the length of the sliding window is equal to the value you specify in Window length. When you clear this check box, the length of the sliding window is infinite. In this mode, the block computes the RMS of the current sample and all the previous samples in the channel.

Dependencies

To enable this parameter, set Method to Sliding window.

Window length — Length of sliding window

4 (default) | positive scalar integer

Specifies the length of the sliding window in samples.

Dependencies

To enable this parameter, set Method to Sliding window and select the Specify window length check box.

Overlap length — Overlap length between windows

3 (default) | nonnegative integer

Specify the overlap length between sliding windows as a nonnegative integer. The value of overlap length varies in the range [0, Window length − 1].

Dependencies

To enable this parameter, set Method to Sliding window and select the Specify window length check box.

Allow arbitrary frame length for fixed-size input signals — Allow arbitrary frame length for fixed-size input signals

off (default) | on

Specify whether fixed-size input signals (whose size does not change during simulation) can have an arbitrary frame length, where the frame length does not have to be a multiple of the hop size. Hop size is defined as Window length − Overlap length. The block uses this parameter setting only for fixed-size input signals and ignores this parameter if the input has a variable-size.

When the input signal is a variable-size signal, the signal can have arbitrary frame length, that is, the frame length does not have to be a multiple of the hop size.

For fixed-size input signals, if you:

• Select the Allow arbitrary frame length for fixed-size input signals parameter, the frame length of the signal does not have to be a multiple of the hop size. If the input is not a multiple of the hop size, then the output is generally a variable-size signal. Therefore, to support arbitrary input size, the block must also support variable-size operations, which you can enable by selecting the Allow arbitrary frame length for fixed-size input signals parameter.

• Clear the Allow arbitrary frame length for fixed-size input signals parameter, the input frame length must be a multiple of the hop size.

Dependencies

To enable this parameter, set Method to Sliding window and select the Specify window length check box.

Specify forgetting factor from input port — Flag to specify forgetting factor

off (default) | on

When you select this check box, the forgetting factor is input through the lambda port. When you clear this check box, the forgetting factor is specified on the block dialog through the Forgetting factor parameter.

Dependencies

To enable this parameter, set Method to Exponential weighting.

Forgetting factor — Exponential weighting factor

0.9 (default) | positive real scalar in the range (0,1]

The forgetting factor determines how much weight past data is given. A forgetting factor of 0.9 gives more weight to the older data than does a forgetting factor of 0.1. A forgetting factor of 1.0 indicates infinite memory – all previous samples are given an equal weight.

Tunable: Yes

Dependencies

To enable this parameter, set Method to Exponential weighting and clear the Specify forgetting factor from input port check box.

Simulate using — Type of simulation to run

Code generation (default) | Interpreted execution

Specify the type of simulation to run as one of the following:

• Code generation –– Simulate model using generated C code. The first time you run a simulation, Simulink^® generates C code for the block. The C code is reused for subsequent simulations, as long as the model does not change. This option requires additional startup time but provides faster simulation speed than Interpreted execution.

• Interpreted execution –– Simulate model using the MATLAB^®  interpreter. This option shortens startup time but has slower simulation speed than Code generation.

Block Characteristics

Data Types	double | single
Multidimensional Signals	No
Variable-Size Signals	Yes

Algorithms

Sliding Window Method

In the sliding window method, the output for each input sample is the RMS of the current sample and Len – 1 previous samples. Len is the length of the window in samples. To compute the first output sample, the algorithm waits until it receives the hop size number of input samples. Hop size is defined as window length – overlap length. Remaining samples in the window are considered to be zero. As an example, if the window length is 5 and the overlap length is 2, then the algorithm waits until it receives 3 samples of input to compute the first sample of the output. After generating the first output, it generates the subsequent output samples for every hop size number of input samples.

When you do not specify the window length, the algorithm chooses an infinite window length. In this mode, the output is the moving RMS of the current sample and all the previous samples in the channel.

Consider an example of computing the moving RMS of a streaming input data using the sliding window method. The algorithm uses a window length of 4 and an overlap length of 3. With each input sample that comes in, the window of length 4 moves along the data.

Exponential Weighting Method

In the exponential weighting method, the moving RMS is computed recursively using these formulas:

$$\begin{array}{l}{w}_{N,\lambda }=\lambda w_{N-1,\lambda }+1\\ x\_rm{s}_{N,\lambda }=\sqrt{\left(1-\frac{1}{w_{N,\lambda }}\right)x\_rm{s}_{N-1,\lambda }+\left(\frac{1}{w_{N,\lambda }}\right)x^2{}_N}\end{array}$$

• $$x\_rm{s}_{N,\lambda }$$ — Moving RMS at the current sample

• $$x^2{}_N$$ — Square of the current input data sample

• $$x\_rm{s}_{N-1,\lambda }$$ — Moving RMS at the previous sample

• λ — Forgetting factor

• $$w_{N,\lambda }$$ — Weighting factor applied to the current data sample

• $$\left(1-\frac{1}{w_{N,\lambda }}\right)x\_rm{s}_{N-1,\lambda }$$ — Effect of the previous data on the RMS

For the first sample, where N = 1, the algorithm chooses $$w_{N,\lambda }$$ = 1. For the next sample, the weighting factor is updated and used to compute the RMS, as per the recursive equation. As the age of the data increases, the magnitude of the weighting factor decreases exponentially and never reaches zero. In other words, the recent data has more influence on the current RMS than the older data.

The value of the forgetting factor determines the rate of change of the weighting factors. A forgetting factor of 0.9 gives more weight to the older data than does a forgetting factor of 0.1. A forgetting factor of 1.0 indicates infinite memory. All the previous samples are given an equal weight.

Here is an example of computing the moving RMS using the exponential weighting method. The forgetting factor is 0.9.

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using Simulink® Coder™.

Version History

Introduced in R2016b

R2022b: New Overlap length parameter

Starting in R2022b, you can specify the overlap length between sliding windows using the Overlap length parameter.

R2022b: Support for arbitrary input frame length

Moving RMS block supports input signals with arbitrary frame lengths when the:

• Input signal is a fixed-size signal (frame length does not change during simulation) and you select the Allow arbitrary frame length for fixed-size input signals parameter (if enabled).

• Input signal is a variable-size signal (frame length changes during simulation).

When this block supports an arbitrary frame length input signal, the input frame length does not have to be a multiple of the hop size.