lyap

Parâmetros

A, C

matrizes quadradas de reais, C deve ser simétrica

Descrição

X= lyap(A,C,flag) resolve as equações matriciais de tempo contínuo ou de tempo discreto de Lyapunov:

A'*X + X*A = C       ( flag='c' )
A'*X*A - X = C       ( flag='d' )

Perceba que existe uma única solução se e só se um autovalor de A não é um autovalor de -A (flag='c') ou 1 sobre um autovalor de A (flag='d').

Exemplos

A=rand(4,4);C=rand(A);C=C+C';
X=lyap(A,C,'c');
A'*X + X*A -C
X=lyap(A,C,'d');
A'*X*A - X -C

Ver Também

• sylv — equação de Sylvester
• ctr_gram — controllability gramian
• obs_gram — observability gramian