MATLAB的答案

Ode方程组

3次查看(最近30天)

显示旧的评论

" Vossoughian 2020年4月13日

评论道: darova 2020年4月17日

你好,

我试着用龙格库塔来解一个方程组，我做了一个关于RK2(f,h,x0,y0,xfinal)的函数，试着用指定的IC来解如下所示的方程组。有人可以帮助修复代码，因为我得到错误和代码不工作。

                    清晰的所有
                   
                    clc
                   
                    β= 1/3;
                   
                    γ= 1/7;
                   
                    信谊R S我
                   
                    N = S +我+ R;
                   
                    ode1 = -(β*我* S) / N
                   
                    ode2 = -(β*我* S) / N-gamma *我
                   
                    ode3 =γ*我
                   
                    常微分方程= [ode1、ode2 ode3]
                   
                    为j =常微分方程
                   
                    RK2法(j, 0.2, 0 8 e6, 7)
                   
                    结束
                   
                    函数(x, y) = RK2法(f、h、x0, y0, xfinal)
                   
                    f = matlabFunction (f);
                   
                    流('\n x y '）;
                   
                    o = 1;
                   
                    而x0 < = xn
                   
                    流(' \ n % 4.3 f % 4.3 f ', x0, y0);x和y的%值
                   
                    xs (o) = x0;
                   
                    y (o) = y0;
                   
                    k1 = h * f (x0, y0);
                   
                    x1 = x0 + h;
                   
                    k2 = h * f (x1, y0 + k1);
                   
                    日元= y0 + (k1 + k2) / 2;
                   
                    x0 = x1;
                   
                    y0 = y₁;
                   
                    o = o + 1;
                   
                    结束
                   
                    结束

1评论
显示隐藏 None

darova 2020年4月17日

我建议你试着解决使用问题欧拉方法首先

是这样的:

                          %常数I N
                         
                          ds = @(s) -beta*I*s/N;
                         
                          dt = 0.1;
                         
                          s = 0.1;
                         
                          为i = 2:10 0
                         
                          S (i+1) = S (i) + dt*ds(S (i));
                         
                          结束
                         
                          情节(s)