How can i find equilibrium points of 2 ODEs?
조회 수: 97 (최근 30일)
이전 댓글 표시
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!
댓글 수: 0
채택된 답변
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
추가 답변 (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!