从系统的一分之三返回不正确的图阶常微分方程采用数值

1视图(30天)

显示旧的评论

迈克尔Temesgen 2020年5月2日

评论道: 迈克尔Temesgen2020年5月8日

我使用这个系统的颂歌,它代表的生物膜的增长在一个装运箱(S -浓度限制衬底,u -浮游生物量浓度/免费的细菌,w -区域生物量密度aherent /附着细菌):

产生以下图左边:

我有附加常量变量和所有信息从我的初始值(S) u和w以及我的时间。我采用数值但我似乎无法得到上面的曲线图和得到一个警告(无法满足集成公差没有减少步长低于允许的最小值)。本文使用缩放变量所以我不确定如果我也必须这么做,或者使用实际的初始值(S, u和w条形,ubar, wbar如本文所示)。我附上了链接我用的纸,以防我遗漏了任何重要的信息需要解决我的问题。下面是我的matlab代码在运行文件我用们初始值和数值。

                        函数dxdt =测试仪(t, x)
                       
                        s = x (1);
                       
                        u = x (2);
                       
                        w = x (3);
                       
                        % %常数变量(正确)
                       
                        wmax = 2.78 * 10 ^ 6;
                       
                        γ= 0.5;
                       
                        β= 0.1;
                       
                        s01 = 2.09 * 10 ^ 6;
                       
                        α= 3000;
                       
                        σ= 6/1;
                       
                        一个= 9 * 10 ^ 7;
                       
                        m = 1.66;
                       
                        D = 0.23;
                       
                        k = 0;
                       
                        千瓦= 0;
                       
                        % %方程(正确)
                       
                        W = (W / wmax);
                       
                        fs = ((m * s) / (a + s));
                       
                        GW = ((1 w) / (1.1 - w));
                       
                        % %微分方程(正确)
                       
                        dxdt (1) = (D * (s01-s)) - (((u * fs) +(σ* w * fs)) /γ);
                       
                        dxdt (2) = (u * (fs-D-k)) +(σβ* * w) +(σ* w * fs * (1-GW))(((α* u) * (1 w));
                       
                        dxdt (3) = (w * ((fs * GW) -beta-kw)) +(((α* u) * (1 w)) /σ);
                       
                        dxdt = dxdt ';
                       
                        结束

这是我的运行文件:

                        % %运行文件
                       
                        s01 = 2.09 * 10 ^ 6;
                       
                        一个= 9 * 10 ^ 7;
                       
                        wmax = 2.78 * 10 ^ 6;
                       
                        γ= 0.5;
                       
                        % %初始条件(不确定是否正确)
                       
                        s0 = (s01);
                       
                        情况= 10 ^ 5 *γ*;
                       
                        w0 = 0
                       
                        y0 = (s0) uoh w0);
                       
                        tspan = [0 160];
                       
                        (t, x) =数值(@tester tspan, y0);
                       
                        % %图的密度与时间(incorrect-doesnt匹配纸)
                       
                        情节(t) x (: 1), t, x (:, 2), t, x (:, 3))
                       
                        传奇(“年代”,“u”,' w ')