24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

444239852

银虫 (正式写手)

应助: 20 (小学生)
金币: 275.9
散金: 40
红花: 1
帖子: 555
在线: 268小时
虫号: 1146976
注册: 2010-11-15
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 用四阶龙格库塔法求解常微分方程组已有5人参与

本人是matlab初学者，自己写了程序但一直不能运行，求助高手帮我写一下程序，重谢！！！！

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 微分方程组.docx

2015-02-01 20:50:51, 33.06 K

» 猜你喜欢

1楼 2015-02-01 20:51:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

7888

木虫 (正式写手)

废都绅士

应助: 9 (幼儿园)
金币: 4251.7
红花: 1
帖子: 434
在线: 328.3小时
虫号: 438241
注册: 2007-09-18
性别: GG
专业: 工程热物理相关交叉领域

【答案】应助回帖

引用回帖:

7楼: Originally posted by 月只蓝 at 2015-02-03 14:36:57
求写出来，供我们学习一下，花费您10分钟，我愿意赞助100金币。...

不好意思今天出差刚回来。特意打开以前老电脑，把程序拷出来，仅供参考。

% begin of matlab function RK4
function [Y,X]=RK4(A,B,C,D,u,Ts,X0,nn)
% Disrcrete an continous system using Runge-Kutta-4th Method
% usually used as:
%       [Y(k),X(k)]=RK4(A,B,C,D,u(k),Ts,X(k-1))
%       [Y(k),X]=RK4(A,B,C,D,u(k),Ts,X)
%       where A B C D IS the state space description of an object
% if nargin<=5;
%    ts=.1,X0=zeros(length(A),1);
% end;
if nargin==7
nn=10;
end;
tss=Ts/nn;
for i=1:nn;
k1=(A*X0+B*u);
k2=(A*(X0+tss*k1/2)+B*u);
k3=(A*(X0+tss*k2/2)+B*u);
k4=(A*(X0+tss*k3)+B*u);
X=X0+(k1+k2*2+k3*2+k4)*tss/6;
X0=X;
end;
Y=C*X+D*u;
% end of matlab function RK4

将以上内容存储成RK4.m文件，然后就可以运行了。