| 查看: 610 | 回复: 0 | ||
xiaofenger01金虫 (初入文坛)
|
[求助]
Matlab ode45的运用
|
|
下面这段程序中，T是常数，所以比较简单，如果我想让T随着时间变化，比如写成T=50+0.5*tspan，这样T和tspan一样是1*301的向量了，下面函数MassEquations里的k也是T的函数了，这样的话应该如何改写？ 谢谢 function Cha5demo4 T = 220; R = 8.315; k0 = [5.78052E+10 3.92317E+12 1.64254E+4 6.264E+8]; Ea = [124670 150386 77954 111528]; C0 = [1 0 0 0 0]; tspan = [0:1:300]; [t,C] = ode45(@MassEquations,tspan,C0,[],k0,Ea,R,T) plot(t,C(:,1),'r-',t,C(:,2),'k:',t,C(:,3),'b-.',t,C(:,4),'k--'); function dCdt = MassEquations(t,C,k0,Ea,R,T) k = k0.*exp(-Ea/(R*(T+273.15))); rA = -(k(1)+k(2))*C(1); rB = k(1)*C(1)-k(3)*C(2); rC = k(2)*C(1)-k(4)*C(3); rD = k(3)*C(2)-k(4)*C(4); rE = k(4)*C(3); dCdt = [rA; rB; rC; rD; rE]; |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
微分方程ode45求解，最小二乘法优化微分方程参数，程序运行求助
已经有5人回复
- 求助！！对于调用龙格库塔法如何对结果判断反馈
已经有6人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- 使用matlab最优化方法拟合获得多个动力学参数中的问题
已经有4人回复
- 求助：用matlab或其它软件进行如下计算模拟，多谢
已经有3人回复
- matlab编程 ODE45
已经有4人回复
- 求助用matlab三次样条法算一个反应的反应速率
已经有13人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- 请问什么是分段龙格库塔的方法？需要用matlab软件采用这种方法求二阶微分方程的解
已经有3人回复
- 三元二次的三角方程如何用MATLAB求解？？
已经有10人回复
- matlab的lsqnonlin函数怎么用
已经有7人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- 求助利用matlab的ODE45求解微分方程
已经有16人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】MATLAB导入数据精度问题。。。
已经有4人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
10