[ÇóÖú] ode45

3Â¥: Originally posted by DIScreet123 at 2015-11-15 18:17:49
ллÄúµÄ´ð¸´ ÎÒÏëÎÊÏÂ
function dx=suijishuru(t,x)
       dx(1)=x(2);
       dx(2)=normrnd(25,1,1,1)-25.*x(1)-x(1).^3-0.5.*x(2);
       dx=dx(;
end
x0=;
=ode45('suijishuru',,x0);
plot(t, ...

2Â¥: Originally posted by feixiaolin at 2015-11-15 16:24:33
% ¶¨ÒåX=£¬
% д³ö΢·Ö·½³Ì
dX=inline(¡­¡­);

% ³õʼÌõ¼þ£¬¿ÉÒÔÊDZäÁ¿
X0=;

% ΢·Ö·½³ÌÇó½â
= ode45(dX, , X0);