我在用数值法求解微分方程时老是出错
我要求解微分方程: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=
=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
>> 原帖由 夜幽书生 于 2007-5-23 21:52 发表 http://www.chinavib.com/forum/images/common/back.gif
我要求解微分方程:
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 多讨论 >> t0=-10;tf=20;
tspan=;
x0=;
=ode23('lab',tspan,x0);
plot(x(:,1),x(:,2));
这样就可以啦 谢谢xun的解答
页:
[1]