advanced27 min studyLesson 351 in Mathematics

State-space models and input–output maps

Separate internal state, forcing and measurement while deriving a causal response from the matrix exponential.

On this page

What you’ll understand

  • Derive the forced response of a finite-dimensional continuous-time linear system.
  • Distinguish a transfer matrix from the full initial-state response.
  • Change state coordinates without changing the measured input–output behavior.

A dynamical model needs more than a differential equation for what happens when left alone. It must say how an actuator enters, what a sensor measures, and which internal quantities carry memory. A state-space model keeps these roles separate. The same input–output behavior can have many state coordinate descriptions, and a transfer function can even hide internal modes completely. We begin with finite-dimensional real linear time-invariant systems, explicitly separating the initial condition from the external input before making any stability or design claim.

State, actuator and sensor have different jobs

Write the state equation as ẋ=Ax+Bu and output equation as y=Cx+Du. The state has n components, the input m, and the output p; the matrix dimensions must follow those choices. A describes autonomous motion, B maps actuator commands into state derivatives, C selects measured combinations, and D permits an instantaneous direct input contribution. A physical position and velocity may form a useful state, but state variables need not equal sensor readings. A single sensor can reveal information about several states through their different effects over time.

This is a continuous-time model. A sampled discrete-time model instead uses x_{k+1}=A_d x_k+B_d u_k, and its stability region is the open unit disk rather than the left half-plane. Under a held input over a sampling interval h, A_d=e^{Ah} and B_d is the integral of e^{Aτ}B from zero to h. Replacing A_d by A without a discretization argument mixes units and changes the dynamics. Linearization of a nonlinear plant also has a region of validity; a local approximation does not automatically describe saturated actuators or large excursions.

x˙=Ax+Bu,y=Cx+Du,x(t)=eAtx0+∫0teA(t−τ)Bu(τ) dτ\dot x=Ax+Bu,\quad y=Cx+Du,\qquad x(t)=e^{At}x_0+\int_0^t e^{A(t-\tau)}Bu(\tau)\,d\tau
Variation of constants separates the natural response from the causal forcing response. Locally integrable inputs suffice for this finite-dimensional formula.↔ Scroll wide equations horizontally.

A transfer matrix assumes a zero initial state

With zero initial state, taking Laplace transforms gives G(s)=C(sI−A)^{-1}B+D wherever the resolvent exists and the transforms converge. The impulse response away from time zero is Ce^{At}B; direct feedthrough contributes D times the input immediately, represented by a delta term in an impulse description. A nonzero initial state adds C(sI−A)^{-1}x₀ to the transformed output. It is therefore incorrect to use G alone to predict every response unless the initial state is zero or its natural response has been separately included.

Transfer poles need not include every eigenvalue of A. A mode may be unreachable through B or invisible through C, so algebraic cancellation can remove it from the transfer description. Such a cancellation is not evidence that the internal state is harmless. An unmeasured unstable mode can grow while the measured forced response looks stable. Later lessons identify the controllable and observable part, which is the minimal realization of the transfer behavior. For design, one must decide whether internal stability or only a particular measured response is being analyzed.

G(s)=C(sI−A)−1B+D,Y(s)=G(s)U(s)+C(sI−A)−1x0G(s)=C(sI-A)^{-1}B+D,\qquad Y(s)=G(s)U(s)+C(sI-A)^{-1}x_0
The second term is the initial-state contribution. Omitting it changes the problem rather than simplifying the same response.↔ Scroll wide equations horizontally.

Coordinates change descriptions, not behavior

For an invertible real matrix T, define z=Tx. Then the transformed matrices are TAT^{-1}, TB, CT^{-1}, and D. Substitution proves that the measured y is unchanged for the same input and appropriately transformed initial state. Eigenvalues, controllability rank and observability rank are invariant under this similarity transformation. Individual entries and numerical conditioning are not. Scaling a velocity coordinate from meters per second to millimeters per second can produce very different matrix magnitudes while describing exactly the same physical trajectories.

A useful workflow checks dimensions, specifies initial conditions, derives the response, and only then asks whether a proposed coordinate choice makes the important structure easier to inspect. A transfer description is particularly convenient for frequency response; a state description is particularly convenient for feedback, estimation and hidden-mode analysis. Neither is a universal replacement for the other. When reducing a model, preserving its zero-state transfer map does not preserve every possible initial-state response, because discarded internal variables may have carried memory that the chosen transfer experiment never excited.

Include both forcing and initial memory

For ẋ=−2x+3u, y=x+u, x(0)=1 and constant input u=2, find x(t), y(t) and the zero-state transfer function.

  1. The constant equilibrium under u=2 is x=3. The natural deviation starts at 1−3=−2, so x(t)=3−2exp(−2t).

  2. The output adds the direct input: y(t)=5−2exp(−2t). In particular y(0)=3, while its limiting value is 5.

  3. The zero-state transfer is G(s)=3/(s+2)+1=(s+5)/(s+2)G(s)=3/(s+2)+1=(s+5)/(s+2). Multiplying it by 2/s would omit the additional initial-state term 1/(s+2).

The result

The full output contains direct feedthrough, a forced equilibrium and a decaying initial deviation. The transfer function alone represents only the zero-initial-state experiment.

Explore the model

Interactive lab

Interactive lab · oscillator

A state trajectory carries more than one measured value

Pause and Reset, then set damping ratio ζ=0 and natural frequency ω₀=1. Predict x(0)=1 and velocity 0, then Step to t=0.5: x=cos(0.5) and velocity=−sin(0.5). For a fresh comparison Reset and set ζ=1,ω₀=1; the position is (1+t)exp(−t). These are unforced two-state oscillator responses with a fixed initial condition, not an arbitrary A,B,C,D editor.

A state trajectory carries more than one measured valueUnit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Current readings: Displacement x(t): 1; Velocity v(t): 0; Regime: underdamped. All parameters are adjustable with the labeled controls below.02.557.510−1.15−0.5700.571.15Time · sDisplacementunderdamped
Displacement x(t)Current time
Displacement x(t)
1
Velocity v(t)
0
Regime
underdamped
t = 0 s

Unit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA state includes the memory needed to predict future motion; a single measured position is only one component of that state.

Check your understanding

What extra information is needed to predict output from G(s) when the initial state is nonzero?

Try it yourself

For A=[[0,1],[−2,−3]], B=(0,1)ᵀ, C=(1,0), D=0, derive the scalar zero-state transfer function and explain what changes if the initial position is nonzero.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Eliminate the velocity state to obtain a second-order differential equation for the measured position.

Compare with the solution

The equations give ẏ=x₂ and ÿ=−2y−3ẏ+u. With zero initial state, (s2+3s+2)Y=U(s^2+3s+2)Y=U, so G(s)=1/((s+1)(s+2)). A nonzero initial position or velocity contributes the homogeneous solution of ÿ+3ẏ+2y=0, determined by those initial values. That contribution must be added separately; changing the initial position does not change G.

Take it with you

  • A,B,C,D describe distinct internal, actuation, sensing and direct-feedthrough roles.
  • A transfer matrix excludes nonzero-initial-state memory.
  • Similarity changes state coordinates while preserving the corresponding input–output experiment.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.