lyap

연속시간 랴푸노프 방정식 풀기

구문

X = lyap(A,Q)

X = lyap(A,B,C)

X = lyap(A,Q,[],E)

설명

lyap를 사용하여 랴푸노프 방정식의 특수 형식과 일반 형식을 풉니다. 랴푸노프 방정식은 안정성 이론, 시스템의 RMS(제곱평균제곱근) 동작 연구와 같은 제어의 여러 분야에서 대두됩니다.

X = lyap(A,Q)는 랴푸노프 방정식 $$AX+X{A}^T+Q=0$$의 해를 반환합니다. 여기서 A와 Q는 동일한 크기의 정사각 행렬을 나타냅니다. Q가 대칭 행렬인 경우 해 X도 대칭 행렬입니다.

X = lyap(A,B,C)는 실베스터 방정식 $$AX+XB+C=0$$의 해를 반환합니다. 여기서 입력값 A는 m×m 행렬이고, 입력값 B는 n×n 행렬이며, C와 X는 모두 m×n 행렬입니다.

X = lyap(A,Q,[],E)는 일반화된 랴푸노프 방정식 $$AX{E}^T+EX{A}^T+Q=0$$을 풉니다. 여기서 Q는 대칭 행렬입니다. 이 구문에는 빈 대괄호 []을 사용해야 합니다. 괄호 안에 값을 넣으면 함수가 오류를 반환합니다.

예제

랴푸노프 방정식 풀기

이 예제에서는 랴푸노프 방정식을 푸는 방법을 보여줍니다.

$\mathit{AX}+{\mathit{XA}}^{\mathit{T}}+\mathit{Q}=0$,

여기서

$\mathit{A}=\left[\begin{array}{cc}1& 2\\ -3& -4\end{array}\right]$ 및 $\mathit{Q}=\left[\begin{array}{cc}3& 1\\ 1& 1\end{array}\right]$입니다.

$\mathit{A}$ 행렬은 안정적이고 $\mathit{Q}$ 행렬은 양의 정부호입니다.

행렬을 정의합니다.

A = [1 2; -3 -4];  
Q = [3 1; 1 1];

랴푸노프 방정식을 풀기 위해 lyap 함수를 사용합니다.

X = lyap(A,Q)

X = 2×2

    6.1667   -3.8333
   -3.8333    3.0000

이 함수는 대칭 행렬 $\mathit{X}$를 반환합니다. $\mathit{X}$가 양의 정부호 행렬인지 확인하기 위해 고유값을 계산할 수 있습니다.

eig(X)

ans = 2×1

    0.4359
    8.7308

해는 양의 정부호 행렬입니다.

실베스터 방정식 풀기

이 예제에서는 실베스터 방정식을 푸는 방법을 보여줍니다.

$\mathit{AX}+\mathit{XB}+\mathit{C}=0$,

여기서 $\mathit{A}=5$, $\mathit{B}=\left[\begin{array}{cc}4& 3\\ 3& 3\end{array}\right]$, $\mathit{C}=\left[\begin{array}{cc}2& 1\end{array}\right]$입니다.

행렬을 정의합니다.

A = 5;
B = [4 3; 3 3];
C = [2 1];

lyap 함수를 사용하여 A, B, C의 이러한 값에 대한 실베스터 방정식을 풉니다.

X = lyap(A,B,C)

X = 1×2

   -0.2063   -0.0476

결과는 1×2 행렬입니다.

일반화된 랴푸노프 방정식 풀기

이 예제에서는 일반화된 랴푸노프 방정식 ${\mathit{AXE}}^{\mathit{T}}+{\mathit{EXA}}^{\mathit{T}}+\mathit{Q}=0$을 푸는 방법을 보여줍니다.

난수 값으로 구성된 행렬 $\mathit{A}$와 $\mathit{E}$를 생성합니다.

rng(0)
A = rand(2)

A = 2×2

    0.8147    0.1270
    0.9058    0.9134

E = randn(2)

E = 2×2

    0.3188   -0.4336
   -1.3077    0.3426

복소수 값으로 구성된 대칭 $\mathit{Q}$ 행렬을 생성합니다.

Q = complex(randn(2),randn(2));
Q = Q*Q'

Q = 2×2 complex

  15.6642 + 0.0000i   5.6211 + 4.1271i
   5.6211 - 4.1271i  16.9265 + 0.0000i

일반화된 랴푸노프 방정식을 풉니다.

X = lyap(A,Q,[],E)

X = 2×2 complex

  -2.0000 + 0.0000i  19.6841 - 3.6552i
  19.6841 + 3.6552i  20.9934 + 0.0000i

해 $\mathit{X}$는 2×2 대칭 행렬입니다.

입력 인수

A, Q, B, C, E — 입력 행렬

행렬

입력 행렬로, 다음 크기의 행렬로 지정됩니다.

• A, Q, E — m×m 정사각 행렬

• B — n×n 행렬

• C — m×n 행렬

데이터형: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | fi
복소수 지원 여부: 예

출력 인수

X — 해

행렬

해로, 행렬로 반환됩니다. 랴푸노프 방정식의 경우 해 X는 m×m 정사각 행렬입니다. 실베스터 방정식의 경우 해 X는 C와 동일한 크기의 행렬입니다.

제한 사항

연속 랴푸노프 방정식은 A의 고유값 $${\alpha }_1,{\alpha }_2,\mathrm{...},{\alpha }_n$$과 B의 고유값 $${\beta }_1,{\beta }_2,\mathrm{...},{\beta }_n$$이 모든 쌍 (i,j)에 대해 $${\alpha }_i+{\beta }_j\ne 0$$을 충족할 경우 유일한 해를 갖습니다.

이 조건을 위반하는 경우 lyap는 다음 오류 메시지를 생성합니다.

Solution does not exist or is not unique.

알고리즘

lyap는 랴푸노프 방정식의 경우 SLICOT 루틴 SB03MD 및 SG03AD를, 실베스터 방정식의 경우 SB04MD(SLICOT) 및 ZTRSYL(LAPACK)을 사용합니다.

참고 문헌

[1] Bartels, R. H., and G. W. Stewart. “Algorithm 432 [C2]: Solution of the Matrix Equation AX + XB = C [F4].” Communications of the ACM 15, no. 9 (September 1972): 820–26. https://doi.org/10.1145/361573.361582.

[2] Barraud, A. “A Numerical Algorithm to solveA^{T}XA - X = Q.” IEEE Transactions on Automatic Control 22, no. 5 (October 1977): 883–85. https://doi.org/10.1109/TAC.1977.1101604.

[3] Hammarling, S. J. “Numerical Solution of the Stable, Non-Negative Definite Lyapunov Equation Lyapunov Equation.” IMA Journal of Numerical Analysis 2, no. 3 (1982): 303–23. https://doi.org/10.1093/imanum/2.3.303.

[4] Penzl, Thilo. “Numerical Solution of Generalized Lyapunov Equations.” Advances in Computational Mathematics 8, no. 1 (January 1, 1998): 33–48. https://doi.org/10.1023/A:1018979826766.

[5] Golub, G., S. Nash, and C. Van Loan. “A Hessenberg-Schur Method for the Problem AX + XB= C.” IEEE Transactions on Automatic Control 24, no. 6 (December 1979): 909–13. https://doi.org/10.1109/TAC.1979.1102170.

버전 내역

R2006a 이전에 개발됨

참고 항목

covar | dlyap