tt.m

来自「一个matlab的将军模型」· M 代码 · 共 10 行

function xini = tt(x0,ts)

A     = [-1     1;...
   		1     -3];
eAT = expm(A*ts);
displacement = step_response(A,inv(A),[0;10],ts);
Xk = transform(x0,eAT,displacement);
xini = Xk;

return