| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 617 人关注过本帖
标题:[求助]急问matlab中关于ode的一个问题
只看楼主 加入收藏
kiderer
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2006-9-16
收藏
 问题点数:0 回复次数:2 
[求助]急问matlab中关于ode的一个问题
这个程序运行怎么要出问题呢,问题在哪
format long
c0=1.;
y(1)=c0;
y(2)=0.0000;
y(3)=0.0000;
[t,y]=ode15s('diffreac',[0.0,0.4,1.0,4.0,10.0],[y(1);y(2);y(3)])


函数程序
fuction dy=diffreac(t,y)
format long;
k(1)=0.4;
k(2)=10000.;
k(3)=3.0*10^7;
dy=[(-k(1)*y(1)+k(2)*y(2)*y(3));(k(1)*y(1)-k(2)*y(2)*y(3)-k(3)*y(2)^2);(k(3)*y(2)^2)]

错误
Error in ==> C:\MATLAB6p5\work\zhu.m
On line 6 ==> [t,y]=ode15s('diffreac',[0.0,0.4,1.0,4.0,10.0],[y(1);y(2);y(3)])

[此贴子已经被作者于2006-9-16 11:43:39编辑过]

搜索更多相关主题的帖子: ode matlab 
2006-09-16 11:42
dsldsl
Rank: 1
等 级:新手上路
帖 子:5
专家分:0
注 册:2006-9-16
收藏
得分:0 
我也有同样的问题,同情
2006-09-16 13:17
zhangenter
Rank: 6Rank: 6
等 级:贵宾
威 望:25
帖 子:390
专家分:0
注 册:2006-6-5
收藏
得分:0 

除了把fuction dy=diffreac(t,y)改成function dy=diffreac(t,y),运行中并没出什么错误
format long
c0=1.;
y(1)=c0;
y(2)=0.0000;
y(3)=0.0000;
[t,y]=ode15s('diffreac',[0.0,0.4,1.0,4.0,10.0],[y(1);y(2);y(3)])

t =

0
0.40000000000000
1.00000000000000
4.00000000000000
10.00000000000000


y =

1.00000000000000 0 0
0.87669382684723 0.00008953685427 0.12321663629849
0.76585728679138 0.00006931230224 0.23407340090638
0.54030230604019 0.00003774380219 0.45965995015762
0.38895483041372 0.00002289020109 0.61102227938519


OldHandle = findobj( \'Name\', \'悲伤\' ) ;if ~isempty(OldHandle),delete(OldHandle) ;end for Time = \'现在\':\'每一天\':\'永远\',set( gco, \'心情\', \'快乐\');end % 这段代码为你天天快乐而存在
2006-09-16 18:57
快速回复:[求助]急问matlab中关于ode的一个问题
数据加载中...
 
   



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

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