m_chap2_ex2_1_sub.m

simulink 动态仿真实例m文件

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);