如何解决非线性ode方程?

2视图(30天)

显示旧的评论

哈姆扎卡里姆 2021年9月14日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1453049-how-to-solve-a-non-linear-ode-equation

评论道: Bjorn Gustavsson 2021年9月16日

你好,

我有以下的颂歌,我想解决

k_r和kaq常量值。

我真的困在我应该采取什么样的措施来实现在matalb方程,我尝试使用odetovectorfield函数线性化方程来解决它采用数值解算器但它不会工作以来方程不是准线性微分方程的函数描述。

我将感激任何指导或帮助

提前谢谢!

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

登录置评。

在回答这个问题。

答案(1)

Bjorn Gustavsson 2021年9月14日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/1453049-how-to-solve-a-non-linear-ode-equation answer_787234

dCAdt你有一个二次方程,解决一个标准的二次方程,你会得到两个解为dCAdt可以解决。金宝搏官方网站如果你比较幸运,你可能最终解决方案解决analyticalls足够简单。金宝搏官方网站根据您的初始条件您可能不得不丢弃其中一个解决方案和选择一个来进行。金宝搏官方网站

HTH

4评论
显示3年长的评论隐藏3年长的评论

哈姆扎卡里姆 2021年9月15日

model.pdf

谢谢你的快速响应Bjorn !

你是对的,这是我试过的第一件事。然而,结果远非我所期望的答案。这就是为什么我想知道还有另一个appraoch来解决这个问题。

我附上一个pdf文件,演示如何获得两个表达式为dCA / dt和下面的代码来解决他们,但我不明白为什么我的CA变量不随时间改变。

                                   清晰的所有
                                  
                                   %
                                  
                                   基米-雷克南= 8.6876平台以及;%反应传递系数(m4.mol - 1. - 1)
                                  
                                   C0 = 116.9;%初始铀浓度(mol.m-3)
                                  
                                   kaq = [1.24 1.73 1.69 1.63 1.57 1.43 e-04 e-04 e-04 e-04 e-04 e-04);%水相的传质系数不同停留时间(米/秒)
                                  
                                   ts_aq = (0.033135652 0.019419467 0.013737619 0.011995912 0.010648619 - 0.009601691);%水相停留时间(s)
                                  
                                   S = 1.25 e04;% A / V (m - 1)
                                  
                                   a1 = (S / 2 * Kr)。* kaq。^ 2;
                                  
                                   a2 = (2 * Kr)。/ kaq;
                                  
                                   a3 = 2 * a2;
                                  
                                   Caq_out1 = 0(1、长度(kaq));
                                  
                                   Caq_out2 = 0(1、长度(kaq));
                                  
                                   为i = 1:长度(kaq)
                                  
                                   tspan = [0 ts_aq (i)];
                                  
                                   dCAdt = @ (t, CA) a1(我)* (a2(我* CA-1 +√(1-a3(我)* CA));
                                  
                                   (t, CA) =数值(dCAdt tspan, C0);
                                  
                                   Caq_out1 (i) = CA(长度(CA));
                                  
                                   结束
                                  
                                   为i = 1:长度(kaq)
                                  
                                   tspan = [0 ts_aq (i)];
                                  
                                   dCAdt = @ (t, CA) a1(我)* (a2 (i) * CA-1-sqrt (1-a3(我)* CA));
                                  
                                   (t, CA) =数值(dCAdt tspan, C0);
                                  
                                   Caq_out2 (i) = CA(长度(CA));
                                  
                                   结束
                                  
                                   Efficiency1 = Caq_out1 / C0;
                                  
                                   Efficiency2 = Caq_out2 / C0;

Bjorn Gustavsson 2021年9月16日

有关单位:你可能已经忘记了正确的反应速率从cm³/ s m³/ s或一些这样的容易忘记的一步。这将给你一个6数量级改变反应速率,但很难获取与量纲分析…

登录置评。

在回答这个问题。

类别

MATLAB 数学数值积分和微分方程常微分方程

找到更多的在常微分方程在帮助中心和文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!