[問題] 非線性ODE求解已回收
各位先進大家好,
小弟在解非線性ODE上遇到問題,
想解的題目如下:
X'=a*X/Y*Y'+b*(X-c)
其中a.b.c都是常數,
X、Y都是t的函數,
已知Y(t)= sin(t)+2*cos(t)
要求X(t)
如果要用 m file 以ode45求解,
該怎麼假設會比較好?
我所知道的如下:
function dx=odedT(t,x)
% X(t)設為 x(1)
% dX(t)設為 x(2)
% Y(t)設為 x(3)
% dY(t)設為 x(4)
dx=[a*x(1)/x(3)*x(4)+b*(x(1)-c);
x(1)'=x(2);
sin(t)+2*cos(t);
x(3)'=x(4)];
再以ode45去解
不過這樣都解不出來,
請問有哪個環節出錯嗎??
感謝!!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.118.2.37
→
03/05 10:09, , 1F
03/05 10:09, 1F
→
03/05 10:09, , 2F
03/05 10:09, 2F
推
03/10 18:17, , 3F
03/10 18:17, 3F
討論串 (同標題文章)