龙格-库塔数值

25日视图(30天)

显示旧的评论

阿尔瓦罗·MªZumalacarregui德尔珈朵 2021年3月1日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/759194-runge-kutta-ode45

回答: Meysam Mahooti 2021年5月5日

答:接受 KALYAN ACHARJYA

我想解决两个江户系统龙格库塔方法,我写了这段代码,但我什么都不要,有人知道为什么这是错的吗?如果我使用命令数值越好

函数[]= runge_kutta_sis (f, g, x0, y0, a、b、h)

t = a: h: b;n =长度(t);

x = (x0);y = (y0);

为我= 1:n - 1

k1 = h * f (x(我),(我)、t (i));

l1 = h * g (x(我),(我)、t (i));

k2 = h * f (x(我)+ k1/2 y (i) + l1/2, t(我)+ h / 2);

l2 = h * g (x(我)+ k1/2 y (i) + l1/2, t(我)+ h / 2);

k3 = h * f (x(我)+ k2/2 y (i) + l2/2, t(我)+ h / 2);

l3 = h * g (x(我)+ k2/2 y (i) + l2/2, t(我)+ h / 2);

k4 = h * f (x + k3(我),(我)+ l3, t (i) + h);

l4 = h * g (x + k3(我),(我)+ l3, t (i) + h);

(i + 1) = x (i) + (1/6) * (k1 + 2 * k2 + 2 * k3 + k4);

(i + 1) = y (i) + (1/6) * (l1 + 2 * l2 + 2 * l3 + l4);

结束

情节(t t, x,, y)

我写这篇文章在命令窗口中但是我看不到任何图

> > f = @ (x, y, t) 10 - 0.1 * x * y;

> > g = @ (x, y, t) 20 - 0.2 * x * y;

> > runge_kutta_sis (0 f, g, 500600年,20岁的0.5)

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

KALYAN ACHARJYA 2021年3月1日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/759194-runge-kutta-ode45 answer_636064

小修改: 函数代码(保存在同一工作目录,文件名runge_kutta_sis.m)

                             函数runge_kutta_sis (f, g, x0, y0, a、b、h)
                            
                             t = a: h: b;
                            
                             n =长度(t);
                            
                             x = 0(1、长度(n));
                            
                             y = 0(1、长度(n));
                            
                             x (1) = x0;
                            
                             y (1) = y0;
                            
                             为我= 1:n - 1
                            
                             k1 = h * f (x(我),(我)、t (i));
                            
                             l1 = h * g (x(我),(我)、t (i));
                            
                             k2 = h * f (x(我)+ k1/2 y (i) + l1/2, t(我)+ h / 2);
                            
                             l2 = h * g (x(我)+ k1/2 y (i) + l1/2, t(我)+ h / 2);
                            
                             k3 = h * f (x(我)+ k2/2 y (i) + l2/2, t(我)+ h / 2);
                            
                             l3 = h * g (x(我)+ k2/2 y (i) + l2/2, t(我)+ h / 2);
                            
                             k4 = h * f (x + k3(我),(我)+ l3, t (i) + h);
                            
                             l4 = h * g (x + k3(我),(我)+ l3, t (i) + h);
                            
                             (i + 1) = x (i) + (1/6) * (k1 + 2 * k2 + 2 * k3 + k4);
                            
                             (i + 1) = y (i) + (1/6) * (l1 + 2 * l2 + 2 * l3 + l4);
                            
                             结束
                            
                             情节(t t, x,, y);
                            
                             网格在;

主要的脚本:

                             f = @ (x, y, t) 10 - 0.1 * x * y;
                            
                             g = @ (x, y, t) 20 - 0.2 * x * y;
                            
                             runge_kutta_sis (0 f, g, 500600年,20岁的0.5)

输出: