matlab求微分方程的系数,如何利用matlab求解矩阵系数的二阶微分方程

M=[2,0;0 1 ];                      %质量矩阵

K=[6 -2;-2 4];                   %刚度矩阵

a=0;b=0;

C=a*K+b*M;

dt=0.28;

t=0:dt:2.8;

ft0=zeros(length(K),length(t));

for i=1:length(t)

ft0(1,i)=10;  %在节点4的竖直方向加大小为200N的阶跃力

end

dsp=zeros(length(K),length(t));                                         % 位移

vel=zeros(length(K),length(t));                                             % 速度

acc=zeros(length(K),length(t));                                          % 加速度

%--------------------------------------------------------------------------

%  (2) Newmark

alpha=0.3; beta=0.6;                                        % 稳定条件

acc(:,1)=inv(M)*(ft0(:,1)-K*dsp(:,1)-C*vel(:,1));   % 计算初始加速度 (t=0)

ekk=K+M/(alpha*dt^2)+C*beta/(alpha*dt);  % 计算有效刚度矩阵</