examp3_14.m

来自「matlab的大量实例代码」· M 代码 · 共 16 行

A=[1 2 3;4 5 6; 7 8 0];  B=A'; C=[1, 5, 4; 5, 6, 7; 4, 7, 9];
A0=kron(A,eye(3))+kron(eye(3),B')

c=[C(1,:) C(2,:) C(3,:)];  x0 = inv(A0)*c'

D=[x0(1:3)'; x0(4:6)'; x0(7:9)']

C0=A*D+D*A'

norm(C0-C)

c=reshape(C,9,1); x0=inv(A0)*c; D=reshape(x0,3,3);

X=lyap(A,A',-C)

norm(A*X+X*A'-C)