How do I numerically compute a state-space equation?

• hkBattousai
In summary, the conversation discusses simulating a linear system in a computer by sending input values and receiving corresponding output values. The basic idea is to iterate the equation for the state vector and use it in the equation for the output. To ensure accuracy and convergence, it is recommended to use a numerical method with adaptive step size control, such as the Runge-Kutta method. This allows for a balance between truncation error, rounding error, and stability. One example is the embedded Runge-Kutta method, which uses two approximations with different orders to calculate the step size. There are various textbooks available on the subject of numerical methods for ordinary differential systems.

hkBattousai

Suppose that, I have a linear system which is analytically defined as below:

$\frac{dx}{dt} = Ax + Bu \, ... \, (I)$
$y = Cx + Du \, ... \, (II)$
A, B, C, D are matrices defining the system,
u is input, y is output.

I want to simulate this system in computer (not by using Matlab or any other libraries/tools) by sending input values and receiving corresponding output values. How do I do this? What is the basic idea?

Do I have to iterate the equation (I) by calculating dx/dt and using it to calculate the value of x, then use it in equation (II) to find the output?

What is x vector for in the first place? Why do we calculate it?

Yes, you should use something like a "Runge-Kutta" method to solved the differential equation for x, then use it in the equation for y.

As to "what is x vector for", that would depend entirely upon the problem the system of equations is from!

HallsofIvy said:
Yes, you should use something like a "Runge-Kutta" method to solved the differential equation for x, then use it in the equation for y.

As to "what is x vector for", that would depend entirely upon the problem the system of equations is from!

Is this scenario correct?:

Suppose that, input sample rate is 1 samples/second. I use Runge-Kutta method with a time step (h) which is very low compared to the sample rate and integer divisor of it; for example h=0.1 seconds. I receive an input sample, then apply it Runge Kutta method 10 times (10 = 1/0.1), then receive the next input and so on...

In this case, can we guaranty that the accumulating error of Runge Kutta (or any other derivation method) converges? In other words, will the error remain small, or it will grow up in time?

Could it be by solving it with the equation below?:

$x(t) = e^{At}x(t-T) + \int_{t-T}^{t} e^{A(t - \tau)}Bu(\tau)d\tau$

"T" is the sampling period of the input signal. To calculate the state vector between any two input samples, I will just solve this equation. Would this method work? Even if it works, the integration term requires input signal to be a continuous function; but it originally is discreet. How do I use it as a continuous function? Would interpolation with Legendre polynomials work for estimating the input signal "u "between u(t-T) and u(t)?

hkBattousai said:
Suppose that, input sample rate is 1 samples/second. I use Runge-Kutta method with a time step (h) which is very low compared to the sample rate and integer divisor of it; for example h=0.1 seconds. I receive an input sample, then apply it Runge Kutta method 10 times (10 = 1/0.1), then receive the next input and so on...

In this case, can we guaranty that the accumulating error of Runge Kutta (or any other derivation method) converges? In other words, will the error remain small, or it will grow up in time?

All numerical methods must produce a numerical solution that converges to the real solution as the step size goes towards zero. However, finding a fixed step size that balances truncation error, rounding error and stability is often difficult unless the system being integrated is very smooth. For simple integration methods, like Euler, the requirement for a stable solution may easily drive the step size so small that the solution globally suffers from severe rounding errors and requires much more computational effort than if done using methods of higher order.

In order to find the "right" step size, you would in practice normally employ so-called adaptive step size control when using Runge-Kutta and similar methods which works by estimating and monitoring the truncation error in the numerical solution. For instance, using RK4, a popular fourth-order method, the truncation error should be proportional with h4 so from the same initial state you can calculate the next state using different h (say, h and h/2) and compare the normed difference in state to a set limit. If error is above limit you reduce h, if well below, you can increase h. Most ODE solver libraries implements something like this principle and allows you to specify the relative or absolute error the solver should strive to maintain.

As a further example, one RK variant, called embedded RK, allows you to calculate a step from the same state using two approximations with different orders but only with a minimum of additional computational effort, making adaptive step size control particular attractive for those methods. One common example of such a method is the Runge-Kutta-Feldberg 7/8 (aka RKF78) method.

There are plenty of good textbooks on the subject if you like to read more. I can personally recommend [1], although it is a bit old. There may be more appropriate reference if you are mostly interested in the state-space solution to a given control system.

[1] Numerical Methods for Ordinary Differential Systems, Lambert, Wiley 1991.

Last edited:

1. What is a state-space equation?

A state-space equation is a mathematical model that describes the behavior of a dynamic system with respect to its internal state variables and external inputs. It is commonly used in the field of control theory and can be represented in either continuous or discrete form.

2. How do I convert a system of differential equations into a state-space equation?

To convert a system of differential equations into a state-space equation, you can follow these steps:
1. Identify the state variables of the system
2. Write the differential equations in terms of these state variables
3. Express the equations in matrix form by grouping the variables and their derivatives
4. The resulting matrix equation is the state-space form of the system.

3. What is the purpose of numerically computing a state-space equation?

Numerically computing a state-space equation allows us to simulate the behavior of a system over time and predict its future states. This is useful for designing controllers and analyzing the stability and performance of a system.

4. What are the common methods for numerically computing a state-space equation?

The most commonly used methods for numerically computing a state-space equation are the Euler method, Runge-Kutta methods, and the state-transition matrix method. These methods use numerical integration techniques to approximate the behavior of the system over time.

5. How do I validate the accuracy of my numerically computed state-space equation?

One way to validate the accuracy of a numerically computed state-space equation is to compare the results with analytical solutions, if available. Additionally, you can vary the time step used in the numerical computation and observe the effect on the results. The smaller the time step, the more accurate the results will be.