Вычисление функций от матриц методом суммирования матричных рядов Текст научной статьи по специальности «Математика»

ПРИБОРЫ И СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

УДК 004.056.53

DOI: 10.17586/0021-3454-2017-60-12-1119-1123

ВЫЧИСЛЕНИЕ ФУНКЦИЙ ОТ МАТРИЦ МЕТОДОМ СУММИРОВАНИЯ МАТРИЧНЫХ РЯДОВ

А. И. Коршунов

Военно-морской политехнический институт ВУНЦВМФ „Военно-морская академия им. Н. Г. Кузнецова",

198514, Санкт-Петербург, Россия E-mail: a.i.korshunov@mail.ru

Рассмотрены особенности вычисления в общем виде аналитической функции матричного аргумента, заданной сходящимся бесконечным рядом. Предложен метод вычисления переходной матрицы линейной стационарной системы и других функций от матриц, использующий суммирование матричных рядов. Метод основан на использовании равенства аналитической функции матричного аргумента, заданной бесконечным рядом, сходящимся на спектре матрицы, и полинома от матрицы, совпадающего на спектре матрицы с аналитической функцией. Рассмотрен пример вычисления переходной матрицы линейной стационарной системы с использованием преобразования Лапласа, а также предлагаемого метода. Продемонстрирована большая простота вычисления предлагаемым методом по сравнению с широко используемым в инженерной практике методом, основанным на преобразовании Лапласа.

Ключевые слова: аналитическая функция матричного аргумента, переходная матрица, вычисление, суммирование матричных рядов

Современные ЭВМ позволяют вычислять значение функции от матрицы, определенной бесконечным рядом, практически с любой нужной точностью. Примером служит часто используемая в исследованиях систем автоматического управления [1] переходная матрица системы линейных дифференциальных уравнений с постоянными коэффициентами:

<х fi

H (t) = exp( At) = ^ A1-, i=0 i!

где А — квадратная матрица системы дифференциальных уравнений, t — время.

Однако в аналитических исследованиях во многих случаях необходимо использовать аналитическое выражение Н(/). Обычно переходную матрицу в виде конечного аналитического выражения получают методом приведения матрицы А к диагональной (квазидиагональной) форме или методом преобразования Лапласа [2—5]. В обоих случаях необходимо знать собственные значения матрицы А. Трудоемкость в обоих случаях многократно возрастает при увеличении порядка матрицы. Поэтому разработка менее трудоемких методов вычисления функций от матриц в конечной аналитической форме остается актуальной задачей. Таковым является метод, основанный на суммировании бесконечного матричного ряда, определяющего функцию от матрицы. Его основой служит метод, применяемый при доказательстве сходимости ряда [6]:

/ (А) = / (г)\г=А = а0К + ъ А + а2 А2 +... + атАт +..., (1)

(/(2) = а0 + а12 + а222 +... + ат2т +... — аналитическая, т.е. бесконечное число раз дифференцируемая, функция комплексного переменного 2, заданная рядом с радиусом сходимости р ), и собственных значениях квадратной числовой матрицы а, лежащих внутри круга сходимости ряда:

< р, , =1,2,...,п,

А, — собственные значения матрицы а. Заметим, что функции /(2) и /(А) различаются

тем, что первая принимает комплексные значения, а вторая — матричные, в общем случае комплексные; ее значение определяется подстановкой в тот же полином матрицы А вместо комплексного числа 2. Для квадратной матрицы, как известно [6—8], определены все операции, необходимые при вычислении значения полинома: возведение в степень, умножение на константу и сложение.

В монографии [6, теорема 29] доказано, что матричный ряд (1) сходится и совпадает с многочленом от А: ф(А) степени р < к -1, где к — степень минимального аннулирующего матрицу А многочлена А(2) (многочлена минимальной степени, обращающегося в нулевую

матрицу при подстановке вместо 2 матрицы А), ф(2) = ак^2к-1 +а,к-22к 2 +... + + а0 — полином, заданный на спектре матрицы А. Если известны попарно различные корни полинома А(2) : ^1,^2,...,^ кратности кьк2,...,кг, к1 + к2 +... + кг = к, то коэффициенты полинома ф(2) находятся из системы к уравнений:

а1, а1,

Лф(2)1= Тк1 (2)2' = 0,1,...,к>, ' = 1,2,...,г. (2)

—21 1 —21 1

Во многих случаях минимальный аннулирующий многочлен совпадает с характеристическим полиномом матрицы, например, при ее простых собственных значениях. Однако степень аннулирующего многочлена А( 2) меньше степени характеристического полинома [7, 8]

для с(2) = ёе1;(2Е - А), Е = {1| — единичная матрица, (2Е - А) — характеристическая матрица матрицы А.

В общем случае [7, 8]:

А( 2) = , (3)

где Оп-1( 2) — наибольший общий делитель всех миноров (и-1)-го порядка характеристической матрицы гЕ - А .

Для сравнения предлагаемого метода с наиболее часто применяемым в инженерной практике методом, использующим преобразование Лапласа, вычислим переходную матрицу:

0 1

ао а1

а0, а1 = со^

— X = АХ, (4)

Н(г) = ехр(Аг), А = системы дифференциальных уравнений;

—г

т

X = [х, х2 ] — вектор-столбец.

При заданных начальных условиях X(0) = Хо решение системы (4), как известно [2—5], имеет вид

X (г) = Н (г) Хо, (5)

где H (t) =

= exp( At).

" ИцЦ ) /12(0"

Ь>1 ( ) МО _

Следуя методу преобразования Лапласа, получим:

X (р) = (рЕ - V

откуда обратным преобразованием Лапласа находим

Н(/) = 1-Х {(рЕ - Л)-1}.

(6)

(7)

Процедуру вычисления обратной характеристической матрицы, трудоемкость которой многократно возрастает с увеличением п — порядка матрицы А, можно исключить. Для этого следует решить систему линейных уравнений.

(pE - A) X (p) = Xo, =

x10, x20,..., xj 0,..., xn0

(8)

¡1

что дает х(р) = ^ Ьц(Р)хцо, % 0) = {Ьт/(Р)}, * =1 п.

М

В рассматриваемом примере эта система

Рх1(Р) - х2(Р) = x10, 1 а0 х1 (Р) + (Р + а1) х2 (Р) = х201

имеет решение:

. (р + Д1)х10 + х20 -а0 хЛр) + РХ20 2

Х1(Р)= —-——, Х2(р) = 01^/ ^20, Б(р) = р2 + а1р + ^0.

Отсюда при различных собственных значениях матрицы А ( ^ ^ ^2) следует

hu{t)=l-1 ¡Pia. [= 1 D(p) J

^^^, ^t) = l-1

X1 - X2

_X.it Xot I

h2i(t) = L-1 J-2V ^ = -«0 £:r-T-, h22(t) = L-1 1 P

D( p)

X1 - X 2

1 j = в X1t - ^ D(p)j А1 -2 ' ' = X^ - X2eX2t

D(p) J X1 -2

V2 2

« /4 - « — корни характеристического уравнения p

+ « p + Ü0 = 0 .

В случае А1— X2— X = -«1 /2 («2 = 4«0) предельный переход при X2 ^ X1 дает

МО = ^ + Л М0 =^ = „22(0 = ^ Л

Следуя предлагаемому методу и учитывая, что при простых корнях характеристического уравнения А1, X2 минимальный аннулирующий многочлен А(z) имеет порядок k—n—2,

определим коэффициенты имеющего порядок k—1—1 полинома ф(z) = а^ + а0, заданного на

спектре матрицы At (А^, X2t), из системы уравнений (2):

a1X1t + а0 = в14,

a1X2t + а0 = в Решение этой системы уравнений дает:

X2t

в -в12г

а1 (11 - 12 )г; а0 =

в V - в V " 0 г '

(11 - -12 )г -а^г -а^

Х1гв1 2 - X 2геХ1

(1 - к 2 )г

; Н (г) = ехр( Аг) = ф( Аг) = а1 Аг + а0 Е =

к^в12 - X2гвХ1

(1 -12 )г

1 0 0 1

11в12г -12 в1г к1 - 12 в 1 - в 2

12г

-а

11 - 12

11г в 1 - в

11 - 12 11в11 -^в1 11 - 12

Это совпадает с результатом, полученным преобразованием Лапласа. В случае 11 =12 =1 определяем Бп-1(г) по вычисленным минорам порядка п—1=2—1=1 характеристической матрицы:

г + «1 -ао

гЕ - А =

" г -1 "

, или

а0 ^ + а1

1

Очевидно, что как наибольший общий множитель миноров г+а1, —а0, 1, г — Вп-1=1. Согласно формуле (3), многочлен Д(г) совпадает с характеристическим полиномом и, следовательно, имеет порядок к=2, а полином ф(г) имеет порядок £-1=2-1=1, т.е. ф(2)=а1^+а0. Система уравнений для определения коэффициентов а1 и а0, согласно (2), имеет вид:

а11г + а0 = в1г, |

Хг

откуда

а1 = в

а1 = в1г, а0 = (1 - 1г )в1г.

Переходная матрица в этом случае получена с учетом 1=-а1/2 в виде:

Н (г) = в

—а0г —а1

+ (1 - 1г )в

1 0 0 1

1+^ г I в

1г

—а^в

гв

а

1г

1 г в

что совпадает с результатом, полученным предельным переходом.

Как показал рассмотренный пример, в случае простых корней матрицы предлагаемый метод менее трудоемок, он не требует таких трудноформализуемых операций, как обратное преобразование Лапласа. Несколько усложняется использование метода при наличии кратных собственных значений матрицы вследствие необходимости вычисления алгебраических дополнений элементов характеристической матрицы и определения их наибольшего общего делителя.

г

список литературы

1. Коршунов А. И. Численный метод построения реакции автоматической системы // Изв. вузов СССР. Приборостроение. 1985. № 29. С. 16—22.

2. Бромберг П. В. Матричные методы в теории линейного и импульсного регулирования. М.: Наука, 1967. 324 с.

3. Сю Д., Мейер А. Современная теория автоматического управления и ее применение. М.: Машиностроение, 1972. 544 с.

4. Современная теория систем управления / Под ред. К. Т. Леондеса. М.: Наука, 1970. 511 с.

5. Деруссо П., Рой Р., Клоуз Ч. Пространство состояний в теории управления. М.: Наука, 1970. 620 с.

6. Понтрягин Л. С. Обыкновенные дифференциальные уравнения. М.: Наука, 1965. 332 с.

7. Ланкастер П. Теория матриц. М.: Наука, 1977. 280 с.

8. Гантмахер Ф. Р. Теория матриц. М.: Наука, 1988. 552 с.

Сведения об авторе

Анатолий Иванович Коршунов — д-р техн. наук, профессор; Военно-морской политехнический институт ВУНЦ ВМФ „Военно-морская академия им. Н. Г. Кузнецова", кафедра радиоэлектроники; E-mail: a.i.korshunov@mail.ru

Рекомендована кафедрой радиоэлектроники

Поступила в редакцию 21.03.17 г.

Ссылка для цитирования: Коршунов А. И. Вычисление функций от матриц методом суммирования матричных рядов // Изв. вузов. Приборостроение. 2017. Т. 60, № 12. С. 1119—1123.

CALCULATION OF FUNCTIONS FROM MATRIX BY METHOD OF SUMMING OF MATRIX SERIES

The features of computing of an analytic function of a matrix argument given by a convergent infinite series are considered in the general form. A method is proposed for calculating the transition matrix of a linear stationary system and other functions of matrices, using summation of matrix series. The method is based on the use of the equality of the analytic function of the matrix argument given by an infinite series that converges on the spectrum of the matrix, and the polynomial from the matrix that coincides on the spectrum of the matrix with the analytic function. An example of the calculation of the transition matrix of a linear stationary system using the Laplace transform and the proposed method is considered. The comparison with the method widely used in engineering practice, based on the Laplace transform, demonstrates a much greater simplicity of calculation by the proposed method.

Keywords: analytical functions of matrix argument, transitional matrix, calculation, summation of matrix series

For citation: Коrshunov А. I. Calculation of functions from matrix by method of summing of matrix series. Journal of Instrument Engineering. 2017. Vol. 60, N 12. P. 1119—1123 (in Russian).

DOI: 10.17586/0021-3454-2017-60-12-1119-1123

А. I. Korshunov

N. G. Kuznetsov Naval Academy, Naval Polytechnic Institute, 198514, St. Petersburg, Russia E-mail: a.i.korshunov@mail.ru

Anatoly I. Korshunov

Data on author

Dr. Sci., Professor; N. G. Kuznetsov Naval Academy, Naval Polytechnic Institute, Department of Radio-Electronics; E-mail: a.i.korshunov@mail.ru