Control Systems Assignment Solutions
Control Systems Assignment Solutions
The stability of a nonlinear system can be determined using Lyapunov's method by constructing a suitable Lyapunov function that is positive definite and demonstrating that its time derivative along system trajectories is negative definite (or at least non-positive). For instance, in the system \(\dot{x}_1 = -x_2 -x_1(1-x_1^2-x_2^2)\), \(\dot{x}_2 = x_1 -x_2(1-x_1^2-x_2^2)\), the Lyapunov function \( V = \frac{1}{2}(x_1^2 + x_2^2) \) shows \( \dot{V} < 0 \) for the specific region \( \|x\|_2 \leq r^2 \), indicating local stability in that region .
A system is globally asymptotically stable if every solution trajectory tends toward the equilibrium point as time approaches infinity, regardless of initial conditions. In contrast, a system is locally asymptotically stable if only trajectories originating within a certain neighborhood of the equilibrium move towards it over time. Thus, global stability implies stronger robustness to initial conditions compared to local stability .
A stable node is characterized by all eigenvalues of the system's Jacobian matrix at the equilibrium point having negative real parts, indicating all perturbations decay over time and the system returns to equilibrium. Conversely, a saddle node features eigenvalues with both positive and negative real parts, leading to some diverging trajectories, resulting in instability .
When eigenvalues of a Jacobian matrix have both positive and negative parts, it indicates the presence of a saddle node. A saddle node implies that the system is unstable because some trajectories diverge due to the directions associated with positive real parts of the eigenvalues, despite other directions being stable .
The Jacobian matrix evaluated at equilibrium points gives the linear approximation of the system at those points. Analyzing the eigenvalues of the Jacobian matrix helps determine stability: if all real parts of eigenvalues are negative, the point is a stable node; otherwise it is an unstable or saddle node. For the given system, at \((0,0)\), eigenvalues \(\lambda = -1 \pm j\sqrt{3}/2\) indicate stability. At \((\sqrt{6},0)\) and \((-\sqrt{6},0)\), eigenvalues \(\lambda = -2, 1\) indicate saddle nodes, implying instability at these points .
Lyapunov's stability theory is used to show that the system is stable in the sense of Lyapunov by constructing a Lyapunov function \( V = a(1 - \cos x_1) + 0.5x_2^2 \). Differentiating \( V \) along system trajectories gives \( \dot{V} = -bx_2^2 \), which is non-positive, indicating that \( V \) doesn't increase over time, thus confirming stability .
Linearization simplifies a nonlinear system around equilibrium points, providing an approximation that captures the dynamics near those points. By evaluating the Jacobian matrix and examining its eigenvalues for the linearized version, insights are gained into the stability and dynamic responses, such as oscillations or exponential divergence, thus predicting behavior for small perturbations .
Phase portraits graphically represent trajectories of the system in state space, showing how the system behaves near equilibrium points. They help illustrate the stability and nature of equilibria (e.g., nodes, saddles) and provide a visual tool to understand the evolution of system states over time. Different configurations of a phase portrait provide insights into the system's dynamics and how perturbations might evolve .
The equilibrium points are found where both derivatives \( \dot{x}_1 \) and \( \dot{x}_2 \) equal zero. From \( \dot{x}_1 = x_2 = 0 \) and \( \dot{x}_2 = -x_1 + x_3 - \frac{1}{6} \cdot 0 = 0 \), resulting in \( x_1 = 0 \) and \( x_1 = \pm \sqrt{6} \). Thus, the equilibrium points are \( (0, 0) \), \( (\sqrt{6}, 0) \), and \( (-\sqrt{6}, 0) \).
The function \( f(x) = \sin(x) \; \text{sgn}(x) \) is locally Lipschitz because it is not continuously differentiable at \( x = 0 \). The derivatives from the positive and negative sides do not match: \( \lim_{x \to 0^+} f'(x) = 1 \neq \lim_{x \to 0^-} f'(x) = -1 \). This discontinuity prevents it from being globally Lipschitz .