Question

老师您好，我想请问一下这个方程该怎么用matlab的ode45方程解？

其中Mc和De您就先当做一个常数吧，麻烦您指导一下，您当然可以拒绝我，我这样也算是很唐突的，不好意思。





Answer 1

用ode45求解题主的微分方程组，可以按下列步骤：

1、自定义微分方程组函数，func(t,x)

2、确定x0的初值，x0=[0,0,0,0,0];

3、使用ode45函数求解【t,x】值。

[t,x]=ode45(@(t,x)func(t,x),[0,1.0],x0);

4、使用plot函数绘制t—x(t)，t—dx/dt，t—y(t)，t—dy/dt，t—v(t)数值解的图形

编程运行可以得到如下结果（Mc=10，De=1，F=5）

追问

太感谢您了，可以麻烦您发一下具体matlab程序嘛！太感谢了！