examp4_41.m

来自「matlab语言解高等数学中的各种问题」· M 代码 · 共 4 行

A=[8,1,6; 3,5,7; 4,9,2]; Q=[16,4,1; 9,3,1; 4,2,1];
x=lyap(sym(A),-inv(A'),Q*inv(A'))

norm(double(A*x*A'-x+Q))  % 可以证明这样的解没有误差