-1.5000 0.5000 -0.50002.0000 0.5000 -4.5000 0.50002.0000 1.0000 -2.0000 -2.0000>> A=[-2 0.5 -0.5 0.5;0 -1.5 0.5 -0.5;2 0.5 -4.5 0.5;2 1 -2 -2;];[V J]=jordan(sym(A))V =[ 0, 1/2, 1/2, -1/4][ 0, 0, 1/2, 1][ 1/4, 1/2, 1/2, -1/4][ 1/4, 1/2, 1, -1/4]J =[ -4, 0, 0, 0][ 0, -2, 1, 0][ 0, 0, -2, 1][ 0, 0, 0, -2]十四数值法:>> A=[3,-6,-4,0,5;1,4,2,-2,4;-6,3,-6,7,3;-13,10,0,-11,0;0,4,0,3,4];B=[3,-2,1;-2,-9,2;-2,-1,9];C=-[-2,1,-1;4,1,2;5,-6,1;6,-4,-4;-6,6,-3];X=lyap(A,B,C),norm(A*X+X*B+C)X =4.0569 14.5128 -1.5653-0.0356 -25.0743 2.7408-9.4886 -25.9323 4.4177-2.6969 -21.6450 2.8851-7.7229 -31.9100 3.7634