龙格库塔法求解一阶常微分方程组 的Matlab代码 20
方程组：dCa/dt=1.65E-3*(14.93-Ca)-1.70E13*exp(-9064.23/T)*Ca;dT/dt=4.93-0.01655*T+1.117E1...
方程组：
dCa/dt=1.65E-3*(14.93-Ca)-1.70E13*exp(-9064.23/T)*Ca ;
dT/dt=4.93-0.01655*T+1.117E13*exp(-9064.23/T)*Ca;
对于给定初值（Ca=0.8,T=300)，用龙格库塔法求方程组数值解。 展开
dCa/dt=1.65E-3*(14.93-Ca)-1.70E13*exp(-9064.23/T)*Ca ;
dT/dt=4.93-0.01655*T+1.117E13*exp(-9064.23/T)*Ca;
对于给定初值（Ca=0.8,T=300)，用龙格库塔法求方程组数值解。 展开
1个回答
展开全部
clc
f=@(y,t)([1.65e-3*(14.93-y(1))-1.70e13*exp(-9064.23/y(2))*y(1);
4.93-0.01655*y(2)+1.117e13*exp(-9064.23/y(2))*y(1)]);
[Ca T]=ode45(f,[0 1],[.8 300])
f=@(y,t)([1.65e-3*(14.93-y(1))-1.70e13*exp(-9064.23/y(2))*y(1);
4.93-0.01655*y(2)+1.117e13*exp(-9064.23/y(2))*y(1)]);
[Ca T]=ode45(f,[0 1],[.8 300])
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
