| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1320 人关注过本帖
标题:[求助]我在用数值法求解微分方程时老是出错
只看楼主 加入收藏
夜幽书生
Rank: 1
等 级:新手上路
帖 子:19
专家分:0
注 册:2007-5-11
收藏
 问题点数:0 回复次数:3 
[求助]我在用数值法求解微分方程时老是出错

我要求解微分方程:
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
>>


搜索更多相关主题的帖子: 微分方程 数值 求解 
2007-05-23 21:42
abingchem
Rank: 6Rank: 6
等 级:贵宾
威 望:24
帖 子:716
专家分:0
注 册:2004-12-30
收藏
得分:0 

应该没有什么问题:
[t,x]=ode23(@a,[t0,tf],[10000 50000]);

version:7.1


此人已消失
2007-05-24 22:17
夜幽书生
Rank: 1
等 级:新手上路
帖 子:19
专家分:0
注 册:2007-5-11
收藏
得分:0 
谢谢你的解答,但是那明明显示有错啊
2007-05-25 12:24
快速回复:[求助]我在用数值法求解微分方程时老是出错
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.023311 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved