2. Time-Domain Analysis

State-space Review

Formulation of state-space models, state variable selection, matrix representation, and solution methods for LTI systems.

State-Space Review

Hey students! 🎯 Ready to dive into one of the most powerful tools in control engineering? Today we're exploring state-space representation - a mathematical framework that transforms complex control problems into elegant matrix equations. By the end of this lesson, you'll understand how to formulate state-space models, select appropriate state variables, and solve Linear Time-Invariant (LTI) systems. Think of this as learning the "universal language" that modern control engineers use to describe everything from aircraft autopilots to robotic arms! ✈️🤖

Understanding State-Space Fundamentals

State-space representation is like creating a detailed "snapshot" of a dynamic system at any given moment. Instead of dealing with complicated higher-order differential equations, we break down the system into a set of first-order differential equations using something called state variables.

Imagine you're tracking a car's motion 🚗. Traditional methods might use a single second-order equation involving position. But in state-space, we use two first-order equations: one for position and one for velocity. These two variables (position and velocity) completely describe the car's "state" at any time - hence the name!

The beauty of state-space lies in its matrix representation. Every LTI system can be written as:

$$\dot{x}(t) = Ax(t) + Bu(t)$$

$$y(t) = Cx(t) + Du(t)$$

Here, $x(t)$ is our state vector (containing all state variables), $u(t)$ is the input vector, $y(t)$ is the output vector, and A, B, C, D are constant matrices for LTI systems. This compact form can represent systems with multiple inputs and outputs - something that's much messier with traditional transfer functions!

Real-world applications are everywhere. NASA uses state-space models for spacecraft attitude control, where state variables might include angular positions and velocities in three dimensions. Tesla's autopilot system uses state-space representations to track vehicle position, velocity, and acceleration simultaneously.

State Variable Selection and System Formulation

Choosing the right state variables is like picking the best viewpoint to understand a complex machine 🔧. The key principle is that state variables must contain enough information to predict the system's future behavior, given the current state and future inputs.

For mechanical systems, common state variable choices include positions and velocities. Consider a mass-spring-damper system with equation $m\ddot{x} + c\dot{x} + kx = f(t)$. We can choose state variables as $x_1 = x$ (position) and $x_2 = \dot{x}$ (velocity). This gives us:

$$\dot{x_1} = x_2$$

$$\dot{x_2} = -\frac{k}{m}x_1 - \frac{c}{m}x_2 + \frac{1}{m}f(t)$$

For electrical circuits, state variables are typically capacitor voltages and inductor currents - the energy storage elements. In an RLC circuit, you might choose the capacitor voltage and inductor current as states, since these variables cannot change instantaneously and thus "remember" the system's energy state.

The dimension of the state vector equals the system's order, which corresponds to the number of independent energy storage elements. A system with 3 capacitors and 2 inductors (assuming no dependent loops) would have a 5-dimensional state space.

Here's a practical tip: state variables should be linearly independent and minimal. You want just enough variables to capture the system's dynamics without redundancy. Think of it like describing a person's location - you need latitude and longitude, but adding a third coordinate that's just their sum would be redundant and unhelpful.

Matrix Representation and Properties

The state-space matrices tell us everything about system behavior! Let's break down what each matrix represents in our standard form:

The A matrix (system matrix) describes the internal dynamics - how states evolve naturally without external inputs. Its eigenvalues determine system stability: if all eigenvalues have negative real parts, the system is stable. For a simple spring-mass system, the A matrix captures how position affects velocity and how spring/damper forces affect acceleration.

The B matrix (input matrix) shows how external inputs affect each state variable. In our car example, if the input is engine force, B tells us how this force directly impacts velocity (and indirectly affects position through the velocity state).

The C matrix (output matrix) determines what we can measure or observe from the states. If we can only measure position but not velocity, C would be [1 0] for our two-state car system.

The D matrix (feedthrough matrix) represents direct input-to-output coupling. In many physical systems, D is zero because inputs don't instantaneously affect outputs - there's always some dynamic delay.

Controllability and observability are crucial properties determined by these matrices. A system is controllable if we can drive it from any initial state to any final state using appropriate inputs. Mathematically, this requires the controllability matrix $[B \quad AB \quad A^2B \quad ... \quad A^{n-1}B]$ to have full rank.

Observability means we can determine the internal state from output measurements. The observability matrix $[C^T \quad A^TC^T \quad (A^T)^2C^T \quad ... \quad (A^T)^{n-1}C^T]^T$ must have full rank for complete observability.

Modern control systems like the Mars rovers rely heavily on these concepts. Engineers must ensure the rover's navigation system is both controllable (they can steer it where needed) and observable (they can determine its position and orientation from sensor data).

Solution Methods for LTI Systems

Solving state-space equations is like predicting the future behavior of your system 🔮. For LTI systems, we have elegant analytical solutions using matrix exponentials.

The homogeneous solution (when input $u(t) = 0$) is:

$$x(t) = e^{At}x(0)$$

The matrix exponential $e^{At}$ is called the state transition matrix and represents how the system evolves from its initial condition. Computing this involves:

$$e^{At} = I + At + \frac{(At)^2}{2!} + \frac{(At)^3}{3!} + ...$$

For simple 2×2 systems, you can often find closed-form solutions. For larger systems, numerical methods or software tools like MATLAB become essential.

The complete solution including inputs is:

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

This convolution integral shows how the system responds to both initial conditions and input history. The first term represents the natural response, while the integral captures the forced response.

Numerical solution methods are crucial for practical applications. The most common approaches include:

  • Euler's method: Simple but can be inaccurate for stiff systems
  • Runge-Kutta methods: More accurate, widely used in simulation software
  • Matrix exponential computation: Using eigenvalue decomposition or Padé approximation

Modern flight control systems process state-space calculations in real-time using these numerical methods. An F-35 fighter jet's flight computer solves state-space equations hundreds of times per second to maintain stable flight, especially since the aircraft is naturally unstable and requires constant computer correction.

Conclusion

State-space representation transforms complex control problems into manageable matrix equations, providing a unified framework for analyzing multi-input, multi-output systems. We've covered how to select appropriate state variables, construct the characteristic A, B, C, and D matrices, and solve the resulting differential equations. This powerful approach enables modern control applications from autonomous vehicles to spacecraft navigation, making it an essential tool in every control engineer's toolkit.

Study Notes

• State variables: Minimum set of variables needed to completely describe system behavior at any time

• State vector dimension: Equals system order (number of energy storage elements)

• Standard state-space form: $\dot{x} = Ax + Bu$, $y = Cx + Du$

• A matrix: System/plant matrix describing internal dynamics

• B matrix: Input/control matrix showing input effects on states

• C matrix: Output matrix determining measurable quantities

• D matrix: Feedthrough matrix for direct input-output coupling

• State transition matrix: $\Phi(t) = e^{At}$ describes natural system evolution

• Complete solution: $x(t) = e^{At}x(0) + \int_0^t e^{A(t-\tau)}Bu(\tau)d\tau$

• Controllability condition: Controllability matrix $[B \quad AB \quad A^2B \quad ... \quad A^{n-1}B]$ has full rank

• Observability condition: Observability matrix must have full rank

• Stability criterion: All eigenvalues of A matrix must have negative real parts

• Matrix exponential: $e^{At} = I + At + \frac{(At)^2}{2!} + \frac{(At)^3}{3!} + ...$

Practice Quiz

5 questions to test your understanding