无水1324 发表于 2007-10-16 22:01

回复 #30 octopussheng 的帖子

多自由度没有区别，每个一阶方程用一次这个ＲＫ法就可以了

octopussheng 发表于 2007-10-17 08:03

回复 #31 无水1324 的帖子

我是这样想的呀，呵呵，可是小咕说她们老板要她推多自由度的，我当时就晕了！

无水1324 发表于 2007-10-17 12:02

回复 #32 octopussheng 的帖子

那就没有必要推，直接写出来就可以了

咕噜写书不要害人哈？

xdhuang 发表于 2007-12-4 14:44

本人编的4阶ROUGE-KUTTA法程序,非常简单，最能说明问题.

clc;close all;clear all;
%%%%dy/dx=y-x.^2;初值x0=0; y0=1;
%%%%%真实值   y=2+2*x+x.^2-exp(x)
x0=0;
y0=1;
x=x0;
y=y0;
h=0.05;
times=40;

for k=1:times
   
   K1=y-x.^2;
   K2=(y+h/2*K1)-(x+h/2).^2;
   K3=(y+h/2*K2)-(x+h/2).^2;
   K4=(y+h*K3)-(x+h).^2;
   y=y+h/6*(K1+2*K2+2*K3+K4);
   x=x+h;
   y_rk(k,:)=y;
   ytrue(k,:)=2+2*x+x.^2-exp(x);   
end

x=x0+(0:times-1)*h;
plot(x,ytrue); hold on; plot(x,y_rk,'ro');
legend('真实值','4阶R-K估计值')

无水1324 发表于 2007-12-4 17:51

不错！

zhawenhua 发表于 2008-3-28 17:25

谢谢楼主！

zhawenhua 发表于 2008-3-28 17:28

如何上传附件

菜鸟问题，如何上传附件，有几本好书想上传共享。

无水1324 发表于 2008-3-28 17:49

回复 37楼 的帖子

呵呵，建议不要上传书，这样违反论坛规定的。
不过我们可以私下谈哈

cx2000 发表于 2008-4-14 12:49

哈哈
谢谢各位

meggy 发表于 2008-4-16 20:37

回复 2楼 的帖子

正是我要找的
非常感谢啊！:loveliness:

huqingwan 发表于 2008-4-16 21:01

ode45 也可以定步长 tspan=:@)
页: 1 2 [3]
查看完整版本: 有没有ode45的matlab代码