Root mean square (RMS)

Collaboration diagram for Root mean square (RMS):

Functions
void	tpt_rms_f32 (f32_t *aResult, const f32_t *aInData, uint32_t aCount)
	Root Mean Square of the elements of a floating-point vector. More...
void	tpt_rms_q15 (q15_t *aResult, const q15_t *aInData, uint32_t aCount)
	Root Mean Square of the elements of a Q15 vector. More...
void	tpt_rms_q31 (q31_t *aResult, const q31_t *aInData, uint32_t aCount)
	Root Mean Square of the elements of a Q31 vector. More...

Detailed Description

Calculates the Root Mean Square of the elements in the input vector. The underlying algorithm is used:

      aResult = sqrt(((aInData[0] * aInData[0] + aInData[1] * aInData[1] +
                ... + aInData[aCount-1] * aInData[aCount-1]) / aCount));
 

There are separate functions for floating point, Q31, and Q15 data types.

Function Documentation

tpt_rms_f32()

void tpt_rms_f32	(	f32_t *	aResult,
		const f32_t *	aInData,
		uint32_t	aCount
	)

Root Mean Square of the elements of a floating-point vector.

Parameters

[out]	aResult	root mean square value returned here
[in]	aInData	points to the input vector
[in]	aCount	number of samples in input vector

Returns

none

tpt_rms_q15()

void tpt_rms_q15	(	q15_t *	aResult,
		const q15_t *	aInData,
		uint32_t	aCount
	)

Root Mean Square of the elements of a Q15 vector.

Parameters

[out]	aResult	root mean square value returned here
[in]	aInData	points to the input vector
[in]	aCount	number of samples in input vector

Returns

none

Scaling and Overflow Behavior

The function is implemented using a 64-bit internal accumulator. The input is represented in 1.15 format. Intermediate multiplication yields a 2.30 format, and this result is added without saturation to a 64-bit accumulator in 34.30 format. With 33 guard bits in the accumulator, there is no risk of overflow, and the full precision of the intermediate multiplication is preserved. Finally, the 34.30 result is truncated to 34.15 format by discarding the lower 15 bits, and then saturated to yield a result in 1.15 format.

tpt_rms_q31()

void tpt_rms_q31	(	q31_t *	aResult,
		const q31_t *	aInData,
		uint32_t	aCount
	)

Root Mean Square of the elements of a Q31 vector.

Parameters

[out]	aResult	root mean square value returned here
[in]	aInData	points to the input vector
[in]	aCount	number of samples in input vector

Returns

none

Scaling and Overflow Behavior

The function is implemented using an internal 64-bit accumulator. The input is represented in 1.31 format, and intermediate multiplication yields a 2.62 format. The accumulator maintains full precision of the intermediate multiplication results, but provides only a single guard bit. There is no saturation on intermediate additions. If the accumulator overflows, it wraps around and distorts the result. In order to avoid overflows completely, the input signal must be scaled down by log2(aCount) bits, as a total of aCount additions are performed internally. Finally, the 2.62 accumulator is right shifted by 31 bits to yield a 1.31 format value.