lyap.m

来自「包含大量遗传算法程序」· M 代码 · 共 7 行

function X=lyap(A,B,C)  % 注意应该置于 @sym 目录下
if nargin==2, C=B; B=A'; end
[nr,nc]=size(C); A0=kron(A,eye(nc))+kron(eye(nr),B');
try
   C1=C'; x0=-inv(A0)*C1(:); X=reshape(x0,nc,nr)';
catch, error('singular matrix found.'), end