АНАЛИЗ МЕТОДОВ РЕШЕНИЯ ДИСКРЕТНОГО УРАВНЕНИЯ СИЛЬВЕСТРА Гнидкина
advertisement
АНАЛИЗ МЕТОДОВ РЕШЕНИЯ ДИСКРЕТНОГО УРАВНЕНИЯ СИЛЬВЕСТРА Гнидкина А.Н. Липецкий государственный технический университет Липецк, Россия METHODS OF SOLVING DISCRETE SYLVESTER EQUATIONS ANALYSIS Gnidkina A.N. Lipetsk State Technical University Lipetsk, Russia В статье рассматривается применение различных методов к решению дискретного уравнения Сильвестра и приводится анализ полученных результатов. Уравнение Сильвестра возникает во многих приложениях, поэтому его решение является важным в задачах линейной алгебры. Рассмотрим дискретное уравнение Сильвестра: АХG – X=С. А – действительная матрица размера m × m, G - действительная матрица размера n × n и С - действительная матрица m × n заданы, а X - искомая матрица m × n. Теорема [1,2]:Дискретное уравнение Сильвестра однозначно разрешимо тогда и только тогда, когда λλ ≠1 ∀λ ∈ σ(A), ∀λ ∈ σ( ). Найдем решение дискретного уравнения Сильвестра при помощи формы Шура. Напомним, что формой Шура матрицы называют матрицу, полученную в результате унитарного подобия, которая имеет треугольный вид. Трансформируем матрицу А к форме Шура[1,2] при помощи ортогональной матрицы U, а G при помощи ортогональной матрицы V: А 0 … 0 ⎡ ⎤ А А … 0 ⎥, A=U AU=⎢ … … … … … … … . ⎢ ⎥ А А ⎦ ⎣А G G … G ⎡ ⎤ 0 G … G ⎥. G = V GV =⎢ … … … … … … … . ⎢ ⎥ 0 G ⎦ ⎣0 Получим нижнетреугольную матрицу A и верхнетреугольную матрицу G. Диагональные блоки G , ..., G матрицы G имеют порядок 1 либо 2. Умножим обе части уравнения Сильвестра слева на U , а справа на V: (U AU)(U XV)(V GV) − (U XV) = U CV. Обозначим через C = U CV и Y = U XV, запишем уравнение в виде: AYG − Y = C. Это уравнение распадается на простейшие матричные уравнения, которые решаются последовательно как линейные системы. Возвращаясь к предыдущему уравнению, восстанавливаем искомую матрицу X: X = UYV . Рассмотрим решение этого уравнения методом произведения Кронекера[2]. Обозначим за матрицу = ⊗A−I : . Заметим, что вектор левой части уравнения Сильвестра равен произведению матрицы на вектор матрицы Х. Тогда уравнение Сильвестра можно записать в виде: vec (AXG − X) = vec (X). vec (X) = vec (С). Откуда можно выразить вектор Х: vec(X) = vec(С). Используя метод произведения Кронекера, могут возникнуть большие погрешности вычисления при обращении матрицы , так как она плохо обусловлена. Также этот метод требует больше памяти для хранения данных при расчетах. Использование формы Шура, является эффективным средством численного решения уравнения Сильвестра, так как основано на вычислении ортогональных матриц. Этот способ решения уравнения Сильвестра более быстрый и удобный, кроме того он является численно устойчивым алгоритмом. Список литературы: 1. Икрамов Х.Д. Несимметричная проблема собственных значений. Численные методы. – М.: Наука. Гл. ред. физ.-мат. лит., 1991. – [13] с.154156. 2. Antoulas A.C.: Approximation of Large-Scale Dynamical Systems. - Department of Electrical and Computer Engineering Rice University Houston, Texas 772511892, USA. – 2004. [6]– с.145-169.
