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

dy1/dt=998*y1+1998*y2
dy2/dt=-999*y1-1999*y2
初始的y1=1; y2=0;
我编写的程序如下:
function y1y2
y1=1; y2=0;
t1=0.1; t2=0.02;
h1=0.1;h2=0.02;
while(t1<=1)
n=0;
while(n<5)

k11=-999*y1-1999*y2;
k22=-999*y1-1999*(y2+h2*k11/2);
k33=-999*y1-1999*(y2+h2*k22/2);
k44=-999*y1-1999*(y2+h2*k33);
y2=y2+h2*(k11+2*k22+2*k33+k44)/6;
n=n+1;
t2=t2+h2;
end

k1=998*y1+1998*y2;
k2=998*(y1+h1*k1/2)+1998*y2;
k3=998*(y1+h1*k2/2)+1998*y2;
k4=998*(y1+h1*k3)+1998*y2;
y1=y1+h1*(k1+2*k2+2*k3+k4)/6
t1=t1+h1;
w=y1;
y1=y2;
y2=w;
end


可是结果一直都不对,不知道是怎么回事,请高手帮忙看看!不胜感激~

搜索更多相关主题的帖子: 常微分方程 
2007-04-23 18:19
tkovampire
Rank: 1
等 级:新手上路
帖 子:27
专家分:0
注 册:2007-1-16
收藏
得分:0 

那位高手帮帮忙啊,小弟急死了!!!

2007-05-08 21:46
langerxing
Rank: 1
等 级:新手上路
帖 子:5
专家分:0
注 册:2007-5-8
收藏
得分:0 

直接使用ode函数求解就可以

2007-05-10 11:27
快速回复:[求助]解常微分方程组问题
数据加载中...
 
   



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

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