How to find the equilibrium points of dynamics system?
48 views (last 30 days)
Show older comments
Hello! I have this scalar dynamic system.
How to find the equilibrium points in order to check the stability with two methods Lyapunov?
Thanks in advance!
What i do? I solve in matlab this f(x) = 0 and find x=1 x=2 and x=2
0 Comments
Answers (1)
Sam Chak
on 31 Aug 2023
The system has two equilibrium points ( and ). It is possible to evaluate the system's stability through a graphical method. From the vector field, we can observe how the system's state will evolve near these equilibrium points.
[T, X] = meshgrid(-0:5/15:5, 0:3/21:3);
S = - (X - 1).*(X - 2).^2;
L = sqrt(1 + S.^2);
U = 1./L;
V = S./L;
quiver(T, X, U, V, 0.5)
axis tight
xlabel('t'), ylabel('x(t)')
For , it takes an eternity to converge to the equilibrium. For , the trajectory converges in finite time.
fcn = @(t, x) - (x - 1).*(x - 2).^2;
[t, x] = ode45(@(t, x) fcn(t, x), [0 2000], 2.01);
figure(2)
plot(t, x), grid on, xlabel('t'), ylabel('x(t)')
[t, x] = ode45(@(t, x) fcn(t, x), [0 200], 1.99);
figure(3)
plot(t, x), grid on, xlabel('t'), ylabel('x(t)')
0 Comments
See Also
Categories
Find more on Matrix Computations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!