sylv

Parâmetros

A, B, C

três matrizes de reais de dimensões apropriadas

flag 'c' ou 'd'

Escolha o tipo de equação de Sylvester a resolver.

Descrição

X = sylv(A, B, C, 'c') computa X, solução da equação de "tempo contínuo" de Sylvester A*X + X*B = C.

X = sylv(A, B, C, 'd') computa X, solução da equação de "tempo discreto" de Sylvester A*X*B + X = C.

X = -sylv(-A, B, C, 'd') computes X as solution of the real "discrete time" Sylvester equation A*X*B - X = C.

Exemplos

A = rand(4, 4); C = rand(4, 3); B = rand(3, 3);

// Continuous time Sylvester equation:
X = sylv(A, B, C, 'c');
norm(A*X + X*B - C)

// Modified Discrete time Sylvester equation:
X = sylv(A, B, C, 'd');
norm(A*X*B + X - C)

// Real Discrete time Sylvester equation:
X = -sylv(-A, B, C, 'd');
norm(A*X*B - X - C)

Ver Também

• lyap — equação de Lyapunov