龙格库塔法求解一阶常微分方程组 的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])
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
东莞大凡
2024-11-14 广告
标定板认准大凡光学科技,专业生产研发厂家，专业从事光学影像测量仪,光学投影测量仪.光学三维测量仪,光学二维测量仪,光学二维测量仪,光学三维测量仪,光学二维测量仪.的研发生产销售。东莞市大凡光学科技有限公司创立于 2018 年，公司总部坐落于...
点击进入详情页
本回答由东莞大凡提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
