Discrete-Time Integrator (Using Simulink)

Perform discrete-time integration of a signal.

Library

Discrete

Description

The Discrete-Time Integrator block can be used in place of the Integrator block when constructing a purely discrete system.

The Discrete-Time Integrator block allows you to:

Define initial conditions on the block dialog box or as input to the block.
Output the block state.
Define upper and lower limits on the integral.
Reset the state depending on an additional reset input.

These features are described below.

Integration Methods

The block can integrate using these methods: Forward Euler, Backward Euler, and Trapezoidal. For a given step k, Simulink updates y(k) and x(k+1). T is the sampling period (delta T in the case of triggered sampling time). Values are clipped according to upper or lower limits. In all cases, x(0)=IC (clipped if necessary):

Forward Euler method (the default), also known as Forward Rectangular, or left-hand approximation. For this method, 1/s is approximated by T/(z-1). This gives us y(k) = y(k-1) + T * u(k-1).

Let x(k) = y(k), then we have:

x(k+1) = x(k) + T * u(k)   (clip if necessary)
  y(k) = x(k)

With this method, input port 1 does not have direct feedthrough.

Backward Euler method, also known as Backward Rectangular or right-hand approximation. For this method, 1/s is approximated by
T*z/(z-1). This gives us y(k) = y(k-1) + T * u(k).

Let x(k) = y(k-1), then we have:

x(k+1) = y(k)
  y(k) = x(k) + T * u(k)     (clip if necessary)

With this method, input port 1 has direct feedthrough.

Trapezoidal method. For this method, 1/s is approximated by
T/2*(z+1)/(z-1). This gives us y(k) = y(k-1) + T/2 * (u(k) + u(k-1)).

When T is fixed (equal to the sampling period), let
x(k) = y(k-1) + T/2 * u(k-1), then we have:

x(k+1) = y(k) + T/2 * u(k)       (clip if necessary)
  y(k) = x(k) + T/2 * u(k)       (clip if necessary)

Here, x(k+1) is the best estimate of the next output. It isn't quite the state, in the sense that x(k) != y(k).

When T is variable (that is, obtained from the triggering times), we have:

x(k+1) = y(k)
  y(k) = x(k) + T/2 * (u(k) + u(k-1)) (clip if necessary)

With this method, input port 1 has direct feedthrough.

The block icon reflects the selected integration method, as this figure shows.

Defining Initial Conditions

You can define the initial conditions as a parameter on the block dialog box or input them from an external signal:

To define the initial conditions as a block parameter, specify the Initial condition source parameter as internal and enter the value in the Initial condition parameter field.
To provide the initial conditions from an external source, specify the Initial condition source parameter as external. An additional input port appears under the block input, as shown in this figure.

Using the State Port

In two known situations, you must use the state port instead of the output port:

When the output of the block is fed back into the block through the reset port or the initial condition port, causing an algebraic loop. For an example of this situation, see the bounce model.
When you want to pass the state from one conditionally executed subsystem to another, which may cause timing problems. For an example of this situation, see the clutch model.

You can correct these problems by passing the state through the state port rather than the output port. Although the values are the same, Simulink generates them at slightly different times, which protects your model from these problems. You output the block state by selecting the Show state port check box.

By default, the state port appears on the top of the block, as shown in this figure.

Limiting the Integral

To prevent the output from exceeding specifiable levels, select the Limit output check box and enter the limits in the appropriate parameter fields. Doing so causes the block to function as a limited integrator. When the output is outside the limits, the integral action is turned off to prevent integral wind up. During a simulation, you can change the limits but you cannot change whether the output is limited. The output is determined as follows:

When the integral is less than the Lower saturation limit and the input is negative, the output is held at the Lower saturation limit.
When the integral is between the Lower saturation limit and the Upper saturation limit, the output is the integral.
When the integral is greater than the Upper saturation limit and the input is positive, the output is held at the Upper saturation limit.

To generate a signal that indicates when the state is being limited, select the Show saturation port check box. A saturation port appears below the block output port, as shown in this figure.

The signal has one of three values:

1 indicates that the upper limit is being applied.
0 indicates that the integral is not limited.
-1 indicates that the lower limit is being applied.

When the Limit output option is selected, the block has three zero crossings: one to detect when it enters the upper saturation limit, one to detect when it enters the lower saturation limit, and one to detect when it leaves saturation.

Resetting the State

The block can reset its state to the specified initial condition based on an external signal. To cause the block to reset its state, select one of the External reset choices. A trigger port appears below the block's input port and indicates the trigger type, as shown in this figure.

Select rising to trigger the state reset when the reset signal has a rising edge. Select falling to trigger the state reset when the reset signal has a falling edge. Select either to trigger the reset when either a rising or falling signal occurs.

The reset port has direct feedthrough. If the block output is fed back into this port, either directly or through a series of blocks with direct feedthrough, an algebraic loop results. To resolve this loop, feed the block state into the reset port instead. To access the block's state, select the Show state port check box.

Choosing All Options

When all options are selected, the icon looks like this.

Data Type Support

A Discrete-Time Integrator block accepts and outputs real signals of type double.

Parameters and Dialog Box

Integrator method

The integration method. The default is ForwardEuler.

External reset

Resets the states to their initial conditions when a trigger event (rising, falling, or either) occurs in the reset signal.

Initial condition source

Gets the states' initial conditions from the Initial condition parameter (if set to internal) or from an external block (if set to external).

Initial condition

The states' initial conditions. Set the Initial condition source parameter value to internal.

Limit output

If checked, limits the states to a value between the Lower saturation limit and Upper saturation limit parameters.

Upper saturation limit

The upper limit for the integral. The default is inf.

Lower saturation limit

The lower limit for the integral. The default is -inf.

Show saturation port

If checked, adds a saturation output port to the block.

Show state port

If checked, adds an output port to the block for the block's state.

Sample time

The time interval between samples. The default is 1.

Characteristics


Direct Feedthrough	Yes, of the reset and external initial condition source ports
Sample Time	Discrete
Scalar Expansion	Of parameters
States	Inherited from driving block and parameter
Vectorized	Yes
Zero Crossing	One for detecting reset; one each to detect upper and lower saturation limits, one when leaving saturation

[ Previous | Help Desk | Next ]