m_chap2_ex2_sub.m

simulink建模与仿真 基础入门源代码

function xdot=f(t,x)

m=1*eye(3,3);
k=1*[3 -1 0;-1 2 -1;0 -1 3];

a=[zeros(3,3),eye(3,3);-inv(m)*k,zeros(3,3)];   b=[zeros(3,1);inv(m)*[1;0;0]];

xdot=a*x+b*sin(3*t);