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

哪位高手解答下:
初学MATLAB,要用四阶龙格库塔函数解一个二阶微分方程
(注:不调用ode45函数)
程序运行总是出错
请高手指点修改下:
方程是:y''+4y'+6y=1
y(0)=0,y'(0)=1,0<=x<=3,
程序代码如下:
syms x y z h d n m ff f1 f2
f1=sym('dy=z')
f2=sym('dz=-4*z-6*y+1')
x=0:0.1:2;h=0.1;
y(0)=0;z(0)=1;
for n=1:length(x)-1
k1=-4*z(n)-6*y(n)+1;
k2=-4*[z(n)+h*k1/2]-6*[y(n)+h*z(n)/2]+1;
k3=-4*[z(n)+h*k2/2]-6*[y(n)+h*z(n)/2+h^2*k1/4]+1;
k4=-4*[z(n)+h*k3]-6*[y(n)+h*z(n)+h^2*k2/2]+1;
z(n+1)=z(n)+h*(k1+2*k2+2*k3+k4)/6;
y(n+1)=y(n)+h*z(n)+h^2*(k1+k2+k3)/6;
end
错误提示:
??? Subscript indices must either be real positive integers or logicals.

Error in ==> Untitled34 at 6
y(0)=0;z(0)=1;
不知道这里如何修改,请朋友们帮忙,谢谢拉!!!!

搜索更多相关主题的帖子: 微分方程 解高阶 
2006-11-29 20:49
123tianya
Rank: 1
等 级:新手上路
帖 子:18
专家分:0
注 册:2006-11-29
收藏
得分:0 
谢谢斑竹
结果出来了
不过步长太大
我又修改了一下
还可以
太谢谢拉
2006-11-30 11:50
快速回复:求助:有关解高阶微分方程
数据加载中...
 
   



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

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