m_chap2_ex2_sub.m

来自「simulink建模与仿真 基础入门源代码」· M 代码 · 共 9 行

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