我在用数值法求解微分方程时老是出错

夜幽书生

我要求解微分方程:
dx/dt=-0.15*y(t)
dy/dt=-0.1*x(t)
x0=10000,y0=5000
代码:帮忙看一下
function xdot=a(t,x)
xdot=zeros(2,1);
xdot(1)=-0.15*x(2);
xdot(2)=-0.10*x(1);
存为a.m


>> t0=-10;tf=20;
x0=[10000 5000]
[t,x]=ode23('a',t0,tf,x0);
plot(x(:,1),x(:,2))

x0 =

       10000        5000

Warning: Obsolete syntax.  Use ode23(fun,tspan,y0,...) instead.
(Type "warning off MATLAB:odearguments:ObsoleteSyntax" to suppress this warning.)
> In C:\MATLAB6p5\toolbox\matlab\funfun\private\odearguments.m at line 39
  In C:\MATLAB6p5\toolbox\matlab\funfun\ode23.m at line 154
>>

eight

原帖由 夜幽书生 于 2007-5-23 21:52 发表

                               
登录/注册后可看大图

我要求解微分方程:
dx/dt=-0.15*y(t)
dy/dt=-0.1*x(t)
x0=10000,y0=5000
代码:帮忙看一下
function xdot=a(t,x)
xdot=zeros(2,1);
xdot(1)=-0.15*x(2);
xdot(2)=-0.10*x(1);
存为a.m


>> t0= ...

请看看 matlab 帮助中的例子，或者搜索一下本版，有 n 多讨论

xun

>> t0=-10;tf=20;
tspan=[t0,tf];
x0=[10000 5000];
[t,x]=ode23('lab',tspan,x0);
plot(x(:,1),x(:,2));
这样就可以啦

夜幽书生

谢谢xun的解答