题目

用Matlab能解 常微分方程组 吗?我是要计算其数值解
我是要解这样一个运动方程：
[M]{y2}+[B]{y1}+[C]{y}={f}
其中,y1为y的一阶导数,y2为y的二阶导数.[M],[B],[C]都是6*6矩阵,{y},{f},{y1},{y2}都是6*1的向量.
请问,是否可以用ode45计算,怎么做呢?谢谢

提问时间：2021-02-01

答案

大概这样写：首先先建立2个editor第一个用来描述方程,建立M文件；第二个用来算方程第一个editor：%%%%设y(1)=y,y(2)=dyfunction dy=fun(t,y)global M B C fdy=zeros(2,1);dy(1)=y(2);dy(2)=(f-B.*y(2)-C.*y(1)).*(in...