%特征值法syms tA = [0 1; -2 -3];[P D] = eig(A);Q = inv(P);expD = exp(D * t) .* (eye(2));eAt = simple(P * expD * Q)%预解矩阵Fssyms sFs = inv(s*eye(2)-A);%预解矩阵xfdxwz反变换得到矩阵指数函数eAt = ilaplace(Fs,s,t);%化简矩阵指数函数eAt = simple(eAt)%有限项展开法[P D] = eig(A);D = eig(A);%生成xsdtn特矩阵for i = 1 : 1 : 2 vand(:,i)=D.^(i-1);enda = (inv(vand) * exp(D * t));a0 = a(1,:);a1 = a(2,:);eAt = simple(a0 * eye(2) + a1 * A)
运算结果如下: