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

Вычислительные технологии

Том 2, № 1, 1997

МЕТОДЫ ВЫЧИСЛЕНИЯ КОЭФФИЦИЕНТОВ ИХП ИНТЕРВАЛЬНЫХ МАТРИЦ *

Е. М. Смагина Томский государственный университет, Россия

А. Н. Моисеев, С. П. Моисеева Анжеро-Судженский филиал Томского государственного педагогического университета, Россия

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

Введение

При решении задач теории автоматического управления (например, при синтезе модального регулятора) часто возникает необходимость в вычислении коэффициентов характеристических полиномов матриц (например, матрицы динамики системы). В случае решения указанных задач для систем с неопределенными (интервальными) параметрами речь идет о вычислении коэффициентов так называемых интервальных характеристических полиномов (ИХП) интервальных матриц. Однако, как будет показано ниже, точное решение этой задачи требует очень больших вычислительных затрат, особенно для интервальных матриц высокого порядка. Поэтому часто говорят о вычислении интервалов, содержащих коэффициенты ИХП интервальной матрицы.

В данной работе предлагается несколько методов расчета подобных интервалов. Представленные методы являются естественными интервальными расширениями таких известных методов вычисления коэффициентов характеристических полиномов числовых матриц, как метод главных миноров, методы Леверье и Фаддеева [1, 2]. Описанные алгоритмы анализируются с позиций трудоемкости и "узости"("точности") получаемых интервальных значений. Настоящая работа является продолжением исследований, начатых в [3].

1. Основные определения. Постановка задачи

Рассмотрим интервальную (п х п)-матрицу [А]. Элементы такой матрицы являются не числами, а интервалами [а^] = [а^] (г,] = 1,п) с известными нижними а^ и верхними

*© Е. М. Смагина, А.Н.Моисеев, С.П.Моисеева, 1997

ац границами :

([«и] . .. [а 1га]' ......... | . (1.1)

[ага1] • • • [агага],

Определение [4]. Под характеристическим полиномом интервальной матрицы [А] понимают семейство характеристических полиномов ^а(Л) вещественных матриц А £ [А]:

|^(Л) = ае1(Л/га - А) = Лп + Р1А"-1 + ■ ■ ■ + рп, А £ [А]}. (1.2)

Коэффициенты pk = Pk(A), k = 1,n являются полилинейными функциями элементов точечных матриц A Е [A], поэтому для каждого из них можно построить объединенное интервальное расширение [4]:

1 I Pk(A) = [min pk(A), maxpk(A)] = p], k = 1,n, (1.3)

Y Ae[A] Ae[A]

Ae[A]

Тогда вместо семейства характеристических полиномов (1.2) можно рассматривать содержащий его интервальный характеристический полином (ИХП)

P[A](A) = An + [pi]An-1 + ••• + [pn]. (1.4)

Пусть для k = 1, n [pk] обозначает естественное интервальное расширение (или просто интервальное расширение) [5] какой-нибудь функции, вычисляющей значение pk (A), на множество всех матриц из [A]. Тогда имеют место включения

[pk] с [pk], k = т;п. (1.5)

Известно [6, 7], что точное значение [pk] можно определить путем вычисления характеристических полиномов 2n точечных матриц, составленных из различных комбинаций верхних и нижних границ элементов интервальной матрицы, и построения на их множестве минимальных интервалов, содержащих значения соответствующих коэффициентов. Такой метод является очень трудоемким и может быть пригодным только для интервальных матриц низкого порядка (так, например, для интервальной матрицы шестого порядка требуется вычислить коэффициенты нескольких миллиардов характеристических полиномов).

Для вычисления интервалов [pk] есть менее трудоемкие алгоритмы, например естественные интервальные расширения обычных методов вычисления коэффициентов характеристического полинома. Но такие алгоритмы дают разные значения p] в зависимости от выбранного метода вычисления коэффициентов и от способа взятия его интервального расширения. Возникает

Задача. Выбрать метод вычисления коэффициентов характеристического полинома и способ его интервального расширения, дающие при возможно меньшей трудоемкости наименее широкие интервалы [pk], содержащие коэффициенты интервального характеристического полинома <^[A] (A) интервальной матрицы [A].

Далее будут построены интервальные расширения и произведен их анализ для трех методов вычисления коэффициентов характеристических полиномов матриц: метода главных миноров, метода Леверье и метода Фаддеева [1, 2].

Поскольку метод главных миноров и метод Леверье используют понятие следа матрицы, необходимо ввести это понятие и для интервальных матриц.

Определение. Под следом интервальной матрицы 1г[А] будем понимать интервал, равный сумме всех диагональных элементов этой матрицы:

МЛ] = 5>*]. (1.6)

г=1

2. Метод главных миноров

Метод главных миноров, по сути, вытекает непосредственно из определения характеристического полинома матрицы. Действительно, непосредственное раскрытие детерминанта матрицы Л/п — А, представляющего собой характеристический полином матрицы А, и приведение подобных членов по степеням Л дают следующие аналитические формулы для вычисления коэффициентов характеристического полинома рк, к = 1, п:

Р1 = —1гА, (2.1)

Рк = (—1)к^, к = 2, п — 1, (2.2)

Рп = (—1)п А, (2.3)

где Бк представляет собой сумму всех главных миноров [2] порядка к матрицы А. Если в формулах (2.1)-(2.3) заменить все числовые параметры (то есть элементы матрицы А) и операции над ними на интервальные, то получим интервалы р], к = 1,п, представляющие собой интервальное расширение коэффициентов характеристического полинома, определяемых по методу главных миноров

[Р1] = — ИА], (2.4)

[Рк] = (—1)к [Бк ], к = 27п-Г, (2.5)

[Рп] = (—1)пОв1[А], (2.6)

где каждый интервал [Бк] равен сумме всех главных миноров порядка к матрицы [А].

Очевидно, что в описанном методе, в зависимости от способа вычисления миноров и определителя, границы интервалов [р,], к = 2,п, будут различными. Наименее широкие интервалы получаются, если считать миноры и определитель "точным"методом, используя вычисление 2к2 миноров. Но, несмотря на то, что в формулах (2.4)-(2.6) используются только суммы миноров, даже при точном вычислении этих миноров и определителя интервальной матрицы, мы не всегда сможем получить точные значения коэффициентов ИХП матрицы [А].

3. Метод Леверье

Метод Леверье дает следующую рекуррентную формулу для определения коэффициентов характеристического полинома рк числовой матрицы А:

Рк = —1 МАк)+ р^г(Ак-1) + ••• + рк-^гА}, к =1Гп. (3.1)

к

Запишем интервальное расширение для коэффициентов характеристического полинома, вычисленных по формулам (3.1):

[рк] = - ^(1г[А]к + [й]1г[А]к-1 + ■ ■ ■ + Р-1]ИА]), к = ТЯ (3.2)

к

Как видно из (2.4), (3.2), коэффициенты [р1] и [рТ1] совпадают:

п

[Р1] = [р1] = -ИА] = [ай]. (3.3)

4. Метод Фаддеева

Метод Фаддеева [1] предлагает следующую рекуррентную схему для вычисления коэффициентов характеристического полинома:

А1 = А, р1 = -1гАь В1 = А1 - р5!/п,

A2 = ABi, p2 = - 2trA2, B2 = A2 - p)2In

An — ABn-b Pn — trAn

n

(4.1)

Построим интервальное расширение этого метода. В общем виде формулы для вычисления коэффициентов ИХП методом Фаддеева запишутся следующим образом:

[рк] = - Т МА* ], (4.2)

здесь [Ak] = [A][Bk-i], [Bk-i] = [Ak-i] - [pk-i]1n, [Во] = In, k = 1,n. Как видно из представленных формул, равенство (3.3) можно продолжить:

[pi] = [pi] = [pi] = — tr[A]. (4.3)

К сожалению, нельзя сказать ничего определенного о взаимном расположении интервальных коэффициентов [pk], [pk] и [pk] для k > 1. Их узость зависит от многих факторов: порядка системы, номера коэффициента, порядка умножения интервалов и т.д.

Очевидно, что если при вычислении использовать точное значение произведения матриц (вычисленное как объединенное интервальное расширение указанной бинарной операции на множество всех пар матриц определенных размерностей), то коэффициенты, полученные с помощью методов Леверье и Фаддеева, будут более узкими вследствие многократного вынесения общих множителей и субдистрибутивности операции интервального умножения. Однако нельзя утверждать, что все найденные коэффициенты будут совпадать с истинными значениями коэффициентов ИХП. Кроме того, такой способ расчета потребует больших вычислительных затрат, даже более высоких, чем "точный"способ вычисления коэффициентов.

Далее делается попытка оценить узость интервалов [pk], [pk] и [pk] для k > 1 путем подсчета количества умножений, необходимых для вычисления каждого из них.

5. Анализ методов с точки зрения увеличения ширины интервалов

Основная проблема использования интервальных расширений методов, работающих с точечными параметрами, заключается в "размывании"интервалов, что приводит, в конечном счете, к тому, что результирующие интервальные параметры получаются гораздо шире их истинных значений. Основная причина этого явления заключается в том, что операция умножения субдистрибутивна на множестве правильных ("классических") интервалов (то есть интервалов вида [х, X], где х < х):

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

Итак, оценим количество умножений интервалов, необходимых для вычисления каждого из коэффициентов [рк] и [рк] (к = 2,п). Обозначим через и операцию, результатом которой является интервал, а через Ф(и) — количество умножений интервалов, необходимое для вычисления результата операции и. Например, для вычисления одного элемента произведения [С] = [А][В] интервальных матриц [А] и [В] порядка п необходимо выполнить

умножений. А для вычисления определителя интервальной матрицы [А] порядка п методом разложения по строке —

[а]([Ь] + [с]) С [а][Ь] + [а][с].

(5.1)

Ф([с^- ]) = п

(5.2)

^ к!

к=1

умножений.

Таким образом,

Для получения каждого р], к = 2,п — 1, необходимо вычислить сумму СП главных миноров матрицы, следовательно,

к-1

к-1

(5.5)

Для метода Леверье получаем

к-1

Ф([Рк ]) = Ф(ИА]к) + ^Ф([рг]МА]к-*)

г= 1

£-1

= (к - 1)п х п + ^{Ф(р]) + (к - г - 1)п х п +1} =

г=1

£-1 £-1 = п2£ г +(к - 1) + £ф« =

г=1 г=1 £-1 £-2 = п2£ г + (к - 1) + п^ г +-----+ п2 + 1 =

г=1 г=1

2 - 1) 2 (к - 1)(к - 2) 2 М-

п2—-—;—- + п --2-+-----+ п2 + г

22

(к - г + 1)(к - г) - 1)

г=1

£-1

п2 у^ (к - г + 1)(к - г) + к(к - 1) (5

^ 2 2 ' ( ' ; г=1

Коэффициент Метод главных миноров Метод Леверье Метод Фаддеева

2 п=2 2 5 4

2 п=3 6 10 9

3 9 39 45

2 п=4 12 17 16

3 36 67 120

4 40 166 616

2 п=5 20 26 25

3 90 103 175

4 200 256 1075

5 205 510 6475

2 п=6 30 37 36

3 180 147 288

4 600 366 2052

5 1230 730 14400

6 1236 1275 100836

И, наконец, количество умножений для вычисления некоторого коэффициента ИХП методом Фаддеева можно определить следующими рекуррентными формулами:

Ф([Р1]) = 0,

Ф([Рк]) = Ф(1г([А] ■ ([Ак-1] — [Рк-1]/„)}) = п ■ (п + Ф([Ак-1]) + Ф([Рк-1])). (5.7)

Здесь под Ф([Ак-1]) подразумевается количество умножений, необходимых для вычисления одного элемента матрицы [Ак-1 ]. Учитывая что Ф([р5к-1]) = п ■ Ф([Ак-1]), получаем

п + 1

Ф([рк]) = п ■ (п + Ф([рк-1 ]--), к = 2, п. (5.8)

п

В таблице (см. выше) приводятся данные по количеству умножений интервалов для указанных методов применительно к матрицам до шестого порядка. Очевидно, что наименьшее количество умножений интервалов требует метод главных миноров, а наибольшее (среди представленных методов) — метод Фаддеева.

6. Анализ числовых примеров

Представленные в работе методы были протестированы для большого количества примеров интервальных матриц различных порядков. Предлагаем результаты работы описанных методов для интервальных матриц 3-го и 5-го порядков [3].

При вычислении определителей использовалось естественное интервальное расширение разложения определителя по элементам первой строки [7]

п

Бв1[А] = 1)к+1[а1к] ■ [ВД, (6.1)

к=1

где ^^к] — интервальный дополнительный минор, полученный вычеркиванием из матрицы [А] первой строки и к-го столбца.

Итак, рассмотрим интервальную матрицу третьего порядка

( [1, 2] [0,1] [—9, —4]\ [А] = ([—4, —2] 1 2 I . (6.2)

\[—1,2] [3,4] [0,1] )

Ниже приводятся значения коэффициентов ИХП, полученных с использованием различных методов, описанных в работе: методом Леверье

[р1] = [—4, —2], [р2] = [—18, 23], [рз] = [—262.3, 99.7]; (6.3)

методом Леверье с использованием "точного"вычисления произведения матриц

[рх] = [—4, —2], [р2] = [—16, 23], [рз] = [—235.7, 91.7]; (6.4)

методом Фаддеева

[Р1] = [—4, —2], [р2] = [—16.5, 22.5], [рз] = [—242.3, 58.3]; (6.5)

методом Фаддеева с использованием "точного"вычисления произведения матриц

[Р1] = [—4, —2], [р2] = [—14, 22.5], [рз] = [—182.3, 52.3];

(6.6)

методом главных миноров

Ы = [-4,-2], [p2] = [-16, 21], [рз] = [-165, -2];

(6.7)

методом главных миноров с использованием "точного"вычисления всех миноров и определителя матрицы

Можно заметить, что наилучшим, в смысле наименьшей ширины результирующих интервалов, является метод главных миноров, который дает не только самые узкие значения для всех коэффициентов, но и интервалы, достаточно близкие к истинным коэффициентам ИХП. При точном вычислении миноров и определителя интервальной матрицы метод главных миноров в данном примере дает точные значения коэффициентов. Для получения интервалов (6.3)-(6.9) при использованиии компьютера с процессором Intel 80386DX-40 потребовалось следующее время:

(6.3) (метод Леверье) — 0.06 с;

(6.4) (метод Леверье с "точным"умножением матриц) — 25 мин. 40 с;

(6.5) (метод Фаддеева) — 0.05 с;

(6.6) (метод Фаддеева с "точным"умножением матриц) — 8 мин. 32 с;

(6.7) (метод главных миноров) — меньше 0.01 с;

(6.8) (метод главных миноров с "точным"вычислением определителя и миноров) — 0.5

с;

(6.9) (точные значения — перебором матриц) — 1.36 с.

В качестве другого примера рассмотрим матрицу пятого порядка

( 1 2 3 4 7 \

(6.8)

(6.9)

-1 -2 0 2 0 [A] = 0 1 -1 -2 1 0 2 -3 0 2

(6.10)

\-1 1 -5 0 [1,2]/

Значения получаемых интервалов: методом Леверье

[pi] = [0,1], [р2] = [-1,1], [рз] = [37.3, 51.7], [р4] = [-162.5, -40.6], [р5] = [-331.8, 285.3] — за 0.11 с;

методом Фаддеева

[pi] = [0,1], [Р2] = [-2.5, 2], [рз] = [25.2, 64.8], [р4] = [-244.3, 34.2], [р5] = [-995.9, 975.6] — за 0.11 с;

с помощью метода главных миноров

[р1] = [0,1], [р2] = [-2, 2], [рз] = [36, 53],

И = [—131, —84], [р5] = [—32, 34] — за 0.33 с;

методом главных миноров с использованием "точного"вычисления определителя и миноров

ы = [0,1], Ы = [— 2, 2], [рз] = [38,51], И = [—123, —92], [р5] = [—32, 34] — за 11 мин 34 с;

истинные значения коэффициентов ИХП матрицы

Ы = [0,1], Ы = [— 1,1], [рз] = [40,49], [р4] = [—113, —102], [р5] = [—8,10] — за 4 часа 51 мин.

Видно, что метод главных миноров позволяет получить достаточно узкие значения интервальных коэффициентов, близкие к истинным коэффициентам ИХП. Однако в последнем примере значения интервалов для [р2] и [рз], полученные методом Леверье, наиболее узкие. Это связано с тем, что значения Ф([р2]) и Ф([рГ2]), а также Ф([рз]) и Ф([рз]) оказались достаточно близкими; в этом случае увеличение ширины интервалов в значительной степени зависит от порядка умножения и типа интервалов, составляющих матрицу [А].

Анализируя результаты, полученные для большого количества примеров, можно дать следующие рекомендации. Если требуется при минимальных вычислительных затратах получить достаточно "хорошие"(узкие) значения интервалов, содержащих коэффициенты ИХП интервальной матрицы, лучше всего воспользоваться методом главных миноров (раздел 2). Если же вычислительные затраты не так важны, то можно воспользоваться следующим комбинированным методом:

1) вычислить значения интервалов с помощью методов Леверье, Фаддеева и метода главных миноров;

2) выбрать для границ коэффициентов [р2], ... , [рп] те значения, которые обеспечивали бы минимальную ширину каждого интервала.

Выводы

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

Результаты вычислений для конкретных примеров (раздел 6), их анализ, а также анализ алгоритмов с точки зрения увеличения ширины интервалов и трудоемкости (раздел 5) показывают, что достаточно быстрым и качественным методом вычисления интервалов, содержащих коэффициенты ИХП, является метод главных миноров, описанный в разделе 2. Если вычислительные затраты не являются существенными, то можно воспользоваться комбинированным методом, предложенным в разделе 6.

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

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

[1] Гантмахер Ф. Р. Теория матриц. Наука, М', 1988'

[2] Воеводин В. В., Кузнецов Ю. А. Матрицы и вычисления. Наука, М', 1984'

[3] Смагина Е. М., Моисеев А. Н., Моисеева С. П. Анализ двух методов вычисления коэффициентов ИХП интервальных матриц' В "Информатика и процессы управления"' Межвуз' сб' научных работ, Изд-во КГТУ, Красноярск, 1996, вып'З (в печати)'

[4] Дугарова И. В. Применение интервального анализа при проектировании систем управления с неопределенными параметрами: Дис. ... канд. техн. наук. Томск, 1989'

[5] Moore R. E. Interval analysis. Englewood Cliffs, Prentice Hall, N'Y', 1966'

[6] ХлЕБАлин Н. А. Аналитический синтез регуляторов в условиях неопределенности параметров объекта управления: Дис. ... канд. техн. наук. Саратов, 1984'

[7] Дугарова И. В. Определение и анализ интервальных расширений детерминанта интервальной матрицы. Дет в ВИНИТИ Редкол- журн "Изв' РАН, Техн кибернетика", №1З45-В96, 1996'

[8] Алефельд Г., Херцбергер Ю. Введение в интервальные вычисления. Мир, М', 1987'