Re: [問題] 請問非線性聯立微分方程求解已回收
※ 引述《Ak1 (阿強)》之銘言:
: 請教版上各位前輩
: 有一非線性系統
: x1'=x2-x1(x1^2+x2^2-1)
: x2'=-x1-x2(x1^2+x2^2-1)
: 初始值為x1(0)=x2(0)=0
: 想用m file以ODE45求解,請問該怎麼做比較好
: 以往都只解單一非線性方程..
以前在研究所修了"非線性控制".那時後在課本裡面好像有看到這一題
這題大概有三種類型的解.
1.零解(就是你給的那個條件就會這樣)
2.初始條件同號大於零
3.初始條件同號小於零
初始條件異號似乎沒有解答
--------------------------------------------
function pttex137
clc
clear
[t x] = ode45(@pttxxxx,[0 10],[0.001;0.001]);
figure(1)
plot(t,x(:,1),t,x(:,2))
legend('x1','x2')
xlabel('time')
ylabel('x')
[t1 x1] = ode15s(@pttxxxx,[0 10],[-0.001;-0.001]);
figure(2)
plot(t1,x1(:,1),t1,x1(:,2))
legend('x1','x2')
xlabel('time')
ylabel('x')
function f = pttxxxx(t,x)
f = zeros(2,1);
f(1) = x(2)-x(1).*(x(1).^2+x(2).^2-1);
f(2) = x(1)-x(2).*(x(1).^2+x(2).^2-1);
----------------------------------------------
本問題沒有 STIFF 的性質.
--
1.MATLAB programming 2.ASPEN process simulation package
3.FORTRN programming 4.Advance Engineering Mathematics
5.Process Control Theory
6.Chemical Engineering Basic Theory(Kinetic.thermodynamics.transport)
7.Numerical Method and Analysis
8.MATLAB Toolbox.Simulink system basic design
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 113.160.148.219
推
06/01 13:41, , 1F
06/01 13:41, 1F
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 2 之 3 篇):