Systems are most easily modeled using systems of linear constant coefficient differential equations. They can be represented either as a set of state-space equations or as a transfer function in the Laplace domain.
Electrical and Mechanical Systems
Eletrical Systems
Mechanical Systems
Electro-Mechanical Equivalence
It turns out that electrical and mechanical systems are analogous to each other. In other words, given an electrical system, we can convert it into a mechanical system and vice versa. Capacitors act like springs as energy storage, resistors act like dampers which dissipate energy, and inductors act like inertial masses which resist movement. These are clear from their force/voltage differential equations (in the Laplace domain) in Table 1. Under these analogies, forces are like voltages, currents are like velocities, and charge is like position.
Linearization
Because non-linear systems often have dynamics which are complicated to analyze, a standard trick to make them simpler is to linearize them.
Definition 5
Using Definition 5, we can see that around our operating point, we have
State-Space Equations
Definition 6
System variables are variables which depend on either the input or the system's internal state.
Definition 7
We can easily go from State-Space Equations to a transfer function via the Unilateral Laplace transform. After taking the Laplace Transform of both sides of Equation 2, Equation 3,
Phase Variable Form
We can also derive state space equations from their transfer functions. First, we assume that the transfer function comes from the LCCDE
meaning our transfer function will be of the form
Using this intermediary variable, we can now let
Converting this back to the time-domain,
When we do control in State-Space Control, this makes it easier to place the system poles where we want them to be.
Time Domain Solution
Notice that
Combining these two equations, we see that
Notice that Equation 7 is broken into two pieces.
Definition 8
The zero-input response is how the system will behave when no input is supplied.
Definition 9
Controllability
Definition 10
By the Cayley-Hamilton Theorem (see Cayley-Hamilton),
Definition 11
The controllability matrix is
Theorem 1
Observability
Definition 12
Definition 13
The observability matrix is
A theorem analogous to Theorem 1 exists for observability.
In electrical systems, there are three basic components: resistors, capacitors, and inductors. See Table 1 for their Laplace domain relationships. At an electrical node, ∑V=0 by Kirchoff’s Voltage Law, and at an electrical junction, ∑Iin=∑Iout by Kirchoff’s Current Law.
In mechanical systems, there are also three basic components: dampers, springs, and masses. There are also rotational counterparts. See Table 1 for their Laplace domain relationships. At a massless node, ∑F=0 by Newton’s 2nd law. Because we consider dampers and springs are massless, the force at two ends of a damper or spring must be equal. In rotational systems, we can also have a gear train. Rotational impedances are reflected through gear trains by multiplying by (Nsource2Ndest2).
Linearization is when a nonlinear system f(x) is approximated by the first two terms of its Taylor series about a particular operating point.
f(x0+δx)≈f(x0)+∇x∣x0+δxδx
f(x)−f(x0)=δf(x)≈∇x∣x0+δxδx(1)
Equation 1 will hold so long as δx is small enough to be within the linear regime (i.e where the Taylor Series expansion is a good approximation). If f is a multi-variable equation, then Equation 1 becomes
δf(x,u,…)≈∇x∣x0+δxδx+∇u∣u0+δuδu+⋯
The state variables of a system are the smallest set of linear independent system variables that can uniquely determine all the other system variables for all t>0.
One can think of the state variables x as capturing the internal dynamics of the system. The dynamics are described by matrices A (the state-evolution matrix) and B (the input matrix)
dtdx=Ax+Bu
where u is the input to the system. Sometimes the states are not directly observable, but instead the sensor in Figure 2 only provides a linear combination of the states determined by the output matrix C and the feedforward matrix D. Together, Equation 2 and Equation 3 are the state-space equations of the system.
It is possible that ∃M<N such that ∀k≥M,bk=0. In other words, the numerator can have fewer terms than the denominator. We now introduce an intermediary variable X so
Together, Equation 5, Equation 6 are known as Phase Variable Form. Notice that the characteristic polynomial of the A matrix when it is in phase variable form is
Δ(s)=sn+∑i=0N−1aNaisi.
For transfer functions, the time domain solution for a particular input is given by L−1{H(s)U(s)}. How do we do the same for state-space equations? Equation 2 is a inhomogenous, first-order vector ordinary differential equation. If it was a scalar homogenous ODE, then we know the solution would be x(t)=x(0)eat, so for our vector case, let us first define
eAt=∑k=0∞k!1Ak
using the Taylor Series expansion. With this definition, we can solve Equation 2 using integrating factors. If we let e−At be our integrating factor, then multiplying it to both sides of Equation 2 gives
The zero-state response is how the system response to an input when its initial state is x(0)=0. It is the convolution of the input with eAtBu(t)u(t) where u(t) is the unit step.
x(t)=∫0teA(t−τ)Bu(τ)dτ
A system is controllable if for any initial state x0, we can reach a new state xf in finite time with no constraints on the input u.
Let us assume that we have a controllable system and we want to reach the state 0 from x0, and we reach it at time tf. Then using Equation 7,
Notice that if C is invertible, then we can find the c which will recover −x0, but if it is not invertible, then we may not be able to do this.
If Cis invertible, then the system is controllable.
A system is observable if for any initial state x0, we can determine x0 from u(t) and y(t)over a finite time interval.
O=CCA⋮CAn−1.
If Ois invertible, then the system is observable.
Sometimes systems have a time-delay in them. This is equivalent to placing a system before the plant with impulse response δ(t−T) since x(t)∗δ(t−T)=x(t−T). In the Laplace domain, this is the same as the transfer function e−sT as shown in Figure 3.