| 查看: 533 | 回复: 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]; |
» 猜你喜欢
0703 物理化学调剂
已经有5人回复
化学工程321分求调剂
已经有7人回复
0703化学调剂 290分有科研经历,论文在投
已经有7人回复
308求调剂
已经有3人回复
本人考085602 化学工程 专硕
已经有10人回复
焦虑
已经有9人回复
0856求调剂
已经有3人回复
283求调剂
已经有11人回复
深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿)
已经有8人回复
0703化学调剂
已经有7人回复
找到一些相关的精华帖子,希望有用哦~
微分方程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人回复
科研从小木虫开始,人人为我,我为人人













回复此楼
点击这里搜索更多相关资源