Решатель квадратного корня для уравнений Ляпунова дискретного времени
R = dlyapchol(A,B)
X = dlyapchol(A,B,E)
R = dlyapchol(A,B) вычисляет факторизацию Холеского X = R'*R решения X к матричному уравнению Ляпунова:
A*X*A'- X + B*B' = 0
Все собственные значения A матрица должна находиться в открытом диске блока для R чтобы существовать.
X = dlyapchol(A,B,E) вычисляет факторизацию Холеского X = R'*R из X решение уравнения Сильвестра
A*X*A' - E*X*E' + B*B' = 0
Все обобщенные собственные значения (A,E) должен находиться на открытом диске блока для R чтобы существовать.
dlyapchol использует подпрограммы SLICOT SB03OD и SG03BD.
[1] Бартельс, Р. Х. и Г. В. Стюарт, «Решение матричного уравнения AX + XB = C», Запятая ACM, том 15, № 9, 1972.
[2] Хаммарлинг, С.Дж., «Численное решение стабильного, неотрицательного определённого уравнения Ляпунова», IMA J. Num. Anal., Vol. 2, pp. 303-325, 1982.
[3] Пенцл, Т., «Численное решение обобщенных уравнений Ляпунова», «Достижения в комп. математике», т. 8, стр. 33-48, 1998.