chap5_section53.m

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

A=[0,0,1,0;0,0,0,1;-6,4,-0.2,0.1;2,-2,0.05,-0.05];

T=0.02;
F0=eye(4,4);

for j=1:500
    F=F0+(A^j)*(T^j)/prod(1:j);
    F0=F;
end
F