How can i find equilibrium points of 2 ODEs?

95 次查看（过去 30 天）
Volkan Yangin
Volkan Yangin 2020-11-17
回答： John D'Errico 2020-11-17
Hi,
I have 2 ODEs like this:
syms x1 x2
ode1=- 10*x2 - (9668*sin((pi*((317467494955027*(18*x1 - (2961*x2)/125))/(1125899906842624000*pi) + (11691*atan((317467494955027*(18*x1 - (2961*x2)/125))/(1125899906842624*pi) + 2222272464685189/28147497671065600))/(50*pi) + 2222272464685189/28147497671065600000))/180))/5125 - (9668*sin((pi*((11691*atan((317467494955027*(18*x1 - 180*u + (3339*x2)/125))/(1125899906842624*pi) + 2222272464685189/28147497671065600))/(50*pi) + (317467494955027*(18*x1 - 180*u + (3339*x2)/125))/(1125899906842624000*pi) + 2222272464685189/28147497671065600000))/180))/5125 - 236/1025
ode2= - (4771158*sin((pi*((317467494955027*(18*x1 - (2961*x2)/125))/(1125899906842624000*pi) + (11691*atan((317467494955027*(18*x1 - (2961*x2)/125))/(1125899906842624*pi) + 2222272464685189/28147497671065600))/(50*pi) + 2222272464685189/28147497671065600000))/180))/1665625 - (5380242*sin((pi*((11691*atan((317467494955027*(18*x1 - 180*u + (3339*x2)/125))/(1125899906842624*pi) + 2222272464685189/28147497671065600))/(50*pi) + (317467494955027*(18*x1 - 180*u + (3339*x2)/125))/(1125899906842624000*pi) + 2222272464685189/28147497671065600000))/180))/1665625 - 4956/13325
Is there any way to find the equilibrium points of them?
Thx!

请先登录，再回答此问题。

采纳的回答

John D'Errico
John D'Errico 2020-11-17
You have a system of ODEs. An equilibrium point is where the derivative is zero.
So just use solve, applied to the derivatives, thus solving for the point(s) where dx1/dt==0, dx2/dt==0.
If solve does not find a solution, (as I suspect it may fail here, finding no analytical solution) then use a tool such as vpasolve or fsolve.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by