题主给出的微分方程组,用MATLAB求解的话,可以使用ode函数得到其数值解。
求解过程:
第一步,根据微分方程组,自定义其函数,odefun(t,x)
第二步,确定t的范围,如TSPAN=[0 1];
第三步,确定初始条件,确定Y0的初值,即Y0=[100,20];
第四步,使用ode45函数,得到其数值解,即
[t,x] = ode45(@odefun,TSPAN,Y0);
第五步,使用plot函数,绘出t—x(t)和t—y(t)的曲线图,即
plot(t,x(:,1),t,x(:,2));
运行程序后,可以得到如下结果。