微分方程组ODE法求解。
显示 更早的评论
有如下微分方程组，之前有人推荐我用mathematica解，结果可能方程太复杂，一直显示running，运行不出来。所以想用matlab解一下试试。可是本人接触matlab时间不久，还写不了这个代码，所以求个ode解法的代码。:'(
方程组如图所示，附件中也有方程组和我写的mathematica的代码。
采纳的回答
gegif
2022-11-19
syms t
y6=107163*(0.24*(51749*t)^0.2+0.35)^2
y5=0.0526*log(0.0001955*t*y6+1)
y6=matlabFunction(y6);
y5=matlabFunction(y5);
odefunc=@(t,y)[111.732*y5(t)*y(2)-(111.732*y5(t)*y(1))/47500;
(111.732*y5(t)*y(1)+0.0000936*y(3))-(y(2)*(111.732*y5(t)+0.00243))/9175;
(0.00243*y(2)+0.326*111.732*y5(t))*y(4)-(y(3)*(0.0000936+0.326*111.732*y5(t))/(3.9765*10^12 ));
(1+0.326*111.732*y5(t)*y(3))-(y(4)*0.326*111.732*y5(t))/1329500
]
tspan=[2,2.9]; %时间区间
t0=[ 0.0000000155705,0,0.000032,0];
[t,y]=ode45(odefunc,tspan,t0);
plot(t,y(:,1))
方程如果都能够显示表达的话，matlab求解也是挺方便的。mathematica之所以一直busy，应该跟时间区间的选取有关，时间区间大于2.9的时候，matlab也会给出警告
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息
