examp4_39.m

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

A=[8,1,6; 3,5,7; 4,9,2]; Q=[16,4,1; 9,3,1; 4,2,1];
X=dlyap(A,Q)

norm(A*X*A'-X+Q)   % 精度验证