FB_CMA_HistArray

Calculates the frequency distribution for single and multi-channel real-valued time series.

The function block FB_CMA_HistArray calculates the frequency distribution (called histogram in the graphic display) of single-channel and also multi-channel real-valued input data. Each channel is considered independently. The frequency distribution of the cyclically incoming data buffer is calculated for every single channel, both individual values and arrays being allowed as input buffers.

The lower and upper limits as well as the number of classes (also called bins) are transferred for the parameterization. The individual class limits are then distributed in identical intervals across the so defined total range, see Histogram. Values that lie below the lower limit and values that lie above the upper limit are each counted in two additional bins.

The return value is a two-dimensional array with unsigned 64-bit integer values. The first index is the number of the channel and the second index is the number of the histogram bin concerned. The counts of the elements with a value below the lower limit value or above the upper limit value are contained in the first or last bin respectively.

If a histogram counter exceeds a value of 2 to the power of 64, approx. 1.8E19, in the current implementation the counter overruns without generating an error message. With a counting step of 100 microseconds, this would happen after 59 million years at the earliest.

Memory properties

The sample quantity N, which is used to calculate the current result, automatically increases with each new incoming data set, i.e. the function block uses all input values since its instantiation. Resetting of the sample quantity to zero (deleting the internal memory of the FB) is provided by a ResetData() method or, if the CallEx() method is used, by the variable bResetData.

Configuration

The initialization parameters specify the limits for counting samples in the regular histogram bins. They can be individually adjusted for each channel with the Configure() method.

NaN occurrence

If a set of input values contains the special constant NaN, no value is added to the statistics for this channel for this time step, i.e. it is treated as indicator for missing values. No NaN values are expected at the output.

Behavior when processing multi-channel input data

When processing several channels (nChannels > 1), there is a possibility of each channel having different return values. In this case, return values can be queried separately on the function block. If the results from one or more channels are impermissible, but not all channels, the value on the function block corresponds to eCM_InfRTime_AmbiguousChannelResults. If the results of all channels are impermissible, then the value on the function block corresponds to eCM_ErrRTime_ErrornousChannelResults.

A list of return values of all channels can be queried using the method GetChannelErrors().

When processing several subchannels (nSubChannels > 0), particular attention must be paid to the formatting of the input and output data. If the input data consist of a multi-channel result of an upstream function block, the value of nChannels must be adopted; further configuration takes place in this case via the parameter nSubChannels.

Sample: In the statistical consideration (e.g. by FB_CMA_Quantiles) of the frequency channels of a multi-channel spectrum (e.g. FB_CMA_MagnitudeSpectrum), the value of nChannels must be identical to the number of input signals; the number of subchannels nSubChannels corresponds to the length of the spectrum.

Sample implementation

A sample implementation is available under the following link: Histogram.

Inputs and outputs

The input and output buffers correspond to the following definition (Shape). The variable parameters are part of the function block input stInitPars.

Versions

Input buffer (MultiArray input stream)
Element type, number of dimensions, dimension sizes

Output buffer (MultiArray output stream)
Element type, number of dimensions, dimension sizes

Standard version
(nSubChannels = 0)

LREAL, 1,
nChannels

ULINT, 2,
nChannels x nBins+2

Standard version for several data sets
(nSubChannels = 0)

LREAL, 2,
nChannels not specified*

ULINT, 2,
nChannels x nBins+2

Multi-channel version
(nSubChannels > 0)

LREAL, 2,
nChannels x nSubChannels

ULINT, 3,
nChannels x nSubChannels x nBins+2

Multi-channel version for several data sets
(nSubChannels > 0)

LREAL, 3,
nChannels x nSubChannels x not specified*

ULINT, 3,
nChannels x nSubChannels x nBins+2

*: The length of this dimension can be selected as desired and can thus adapt itself to the application or to the output buffer of the preceding algorithm.

VAR_INPUT
    stInitPars       : ST_CM_HistArray_InitPars;         // init parameter
    nOwnID           : UDINT;                            // ID for this FB instance
    aDestIDs         : ARRAY[1..cCMA_MaxDest] OF UDINT;  // IDs of destinations for output
    nResultBuffers   : UDINT := 4;                       // number of MultiArrays which should be initialized for results (0 for no initialization)
    tTransferTimeout : LTIME := LTIME#500US;             // timeout checking off during access to inter-task FIFOs
END_VAR

Input parameters

The input parameters of this function block represent initialization parameters and must already be assigned in the declaration of the FB instance! (Alternatively: Init() method). They may only be assigned once. A change at runtime is not possible.

Output parameters

VAR_OUTPUT
    bError         : BOOL;                           // TRUE if an error occurs. Reset by next method call.
    hrErrorCode    : HRESULT;                        // '< 0' = error; '> 0' = info; '0' = no error/info
    ipErrorMessage : I_TcMessage := fbErrorMessage;  // Shows detailed information about occurred errors, warnings and more.
    nCntResults    : ULINT;                          // Counts outgoing results (MultiArrays were calculated and sent to transfer tray).
END_VAR
  • bError: The output is TRUE if an error occurs.
  • hrErrorCode: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes.
  • ipErrorMessage: Contains more detailed information on the current return value. Refer here to the section Error description and information. This special interface pointer is internally secured so that it is always valid/assigned.

Methods

Call():

The method is called each cycle in order to apply the algorithm to the current input data. The function block waits for input data if the method indicates neither new results nor an error. This is a regular behavior in the process of the analysis chain.

  • Return value: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes.
METHOD Call : HRESULT
VAR_OUTPUT
    bNewResult   : BOOL;       // TRUE every time when outgoing MultiArray was calculated and sent to transfer tray.
    bError       : BOOL;       // TRUE if an error occurs.
    hrErrorCode  : HRESULT;    // '< 0' = error; '> 0' = info; '0' = no error/info
END_VAR
  • bError: The output is TRUE if an error occurs.
  • hrErrorCode: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes. This output is identical to the return value of the method.
FB_CMA_HistArray 1:

If a timeout occurs or no MultiArray buffer is available for the result, then neither the input data nor the result data are lost. They are forwarded on the next call.

CallEx() :

The method is called in each cycle in order to calculate the histogram from the input signal. An alternative method is Call().
The histogram evaluation is generally significantly more computationally demanding than the registration of new input values. Therefore a use of the method Callex() can considerably shorten the runtime, depending on the configured parameters, by only calculating statistic results when they are required.

The function block waits for input data if the method indicates neither new results nor an error. This is a regular behavior in the process of the analysis chain.

  • Return value: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes.
METHOD CallEx : HRESULT
VAR_INPUT
    nAppendData  : UDINT;      // count of data buffers which are appended until calculation (1= calculate always)
    bResetData   : BOOL;       // automatic reset of dataset buffer after each calculation
END_VAR
VAR_OUTPUT
    bNewResult   : BOOL;        // TRUE every time when outgoing MultiArray was calculated and sent to transfer tray.
    bError       : BOOL;        // TRUE if an error occurs.
    hrErrorCode  : HRESULT;     // '< 0' = error; '> 0' = info; '0' = no error/info
END_VAR
  • nAppendData: Defines how many input data buffers are to be collected before a calculation is carried out, because several data blocks are preferably added in order to achieve a precise result.
  • bResetData: If set, the internal data buffer is completely deleted after calculation.
  • bError: The output is TRUE if an error occurs.
  • hrErrorCode: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes. This output is identical to the return value of the method.
FB_CMA_HistArray 2:

If a timeout occurs or no MultiArray buffer is available for the result, then neither the input data nor the result data are lost. They are forwarded on the next call.

Configure() :

By calling this method, the histogram arguments can be reconfigured. This allows fine adjustment of the fMinBinned and fMaxBinned parameters during runtime. The corresponding PLC array must be defined as follows.

METHOD Configure : HRESULT
VAR_INPUT
    pArg     : POINTER TO LREAL;  // pointer to 2-dimensional array (LREAL) of arguments
    nArgSize : UDINT;             // size of arguments buffer in bytes 
END_VAR

The input buffers correspond to one of the following definitions (input shape). The variable parameters are part of the function block input stInitPars. The two parameters to be configured per channel and subchannel are [fMinBinned, fMaxBinned].

Versions

Input buffer (MultiArray input stream)
Element type, number of dimensions, dimension sizes

Identical configuration of all channels and subchannels

LREAL, 1,
2

Channel-specific configuration
(nSubChannels = 0)

LREAL, 2,
nChannels x 2

Subchannel-specific configuration
(nSubChannels > 0)

LREAL, 2,
nSubChannels x 2

Channel and subchannel-specific configuration
(nSubChannels > 0)

LREAL, 3,
nChannels x nSubChannels x 2

METHOD Init : HRESULT
VAR_INPUT
    stInitPars     : ST_CM_HistArray_InitPars;         // init parameter
    nOwnID         : UDINT;       // ID for this FB instance
    aDestIDs       : ARRAY[1..cCMA_MaxDest] OF UDINT;  // IDs of destinations for output
    nResultBuffers : UDINT := 4;  // number of MultiArrays which should be initialized for results (0 for no initialization)
END_VAR

Init():

This method is not usually necessary in a Condition Monitoring application. It offers an alternative to the function block initialization. The Init() method may only be called during the initialization phase of the PLC. It cannot be used at runtime. You are referred to the use of an FB_init method or the attribute 'call_after_init' (see TwinCAT PLC reference). In addition, this facilitates the function block encapsulation.

The input parameters of the function block instance may not be assigned in the declaration if the initialization is to take place using the Init() method.

  • Return value: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes.

ResetData():

The method deletes all data records that have already been added, see Memory property of the function block. If the Call() method is called again after a ResetData(), the internal memory must be replenished in order to calculate a valid result.

  • Return value: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes.
METHOD ResetData : HRESULT
VAR_INPUT
END_VAR

Alternatively the automatic reset in the method CallEx() can be used.

PassInputs():

As long as an FB_CMA_Source instance is called and signal data are thus transferred to a target block, all further function blocks of the analysis chain have to be called cyclically as explained in the API PLC Reference.
Sometimes it is useful not to execute an algorithm for a certain time. For example, some algorithms should be executed only after prior training or configuration. The function block must be called cyclically, but it is sufficient for the data arriving at the function block to be forwarded in the communication ring. This is done using the PassInputs() method in place of the Call() method. The algorithm itself is not called here, and accordingly no result is calculated and no output buffer generated.

  • Return value: If an error occurs, a corresponding error code of the type HRESULT is output. Possible values are described in the List of error codes.
METHOD PassInputs : HRESULT
VAR_INPUT
END_VAR

GetChannelErrors():

The method enables the querying of a list of the channel-specific return values when processing several channels (nChannels > 1). A call is useful in the case that the return value of the function block corresponds to one of the values eCM_InfRTime_AmbiguousChannelResults or eCM_ErrRTime_ErrornousChannelResults.

  • Return value: Information on the reading process of the list of error codes. The value is set to TRUE if the query was successful, otherwise to FALSE.
    METHOD GetChannelErrors : BOOL
VAR_IN_OUT
    aChannelErrors : ARRAY[*] OF HRESULT;
END_VAR
  • aChannelErrors: Error list of the type HRESULT of the length nChannels.

Similar function blocks

The FB_CMA_Quantiles function block calculates the quantiles of input value distributions.

The FB_CMA_MomentCoefficients function block calculates the statistical moment coefficients: average value, standard deviation, skew and kurtosis.

Requirements

Development environment

Target platform

PLC libraries to include

TwinCAT v3.1.4022.25

PC or CX (x86, x64)

Tc3_CM, Tc3_CM_Base

FB_CMA_HistArray 3:

Limited functional scope already available with CM 3.1. See section Compatibility.