examp4_46.m

来自「《高等应用数学问题MATLAB求解》的PDF版及习题解答」· M 代码 · 共 9 行

syms t; A=[-3,-1,-1; 0,-3,-1; 1,2,0]; simple(expm(A*t))

[V,J]=jordan(A)   % Jordan 矩阵变换

J1=[exp(-2*t), t*exp(-2*t), 1/2*t^2*exp(-2*t);
            0,   exp(-2*t),       t*exp(-2*t);
            0,           0,         exp(-2*t)];

A1=simple(V*J1*inv(V))