Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Иохвидов, И. С. Ганкелевы и теплицевы матрицы и формы. Алгебраическая теория

.pdf
Скачиваний:
63
Добавлен:
20.10.2023
Размер:
8.39 Mб
Скачать

170

ПРЕОБРАЗОВАНИЯ

МАТРИЦ

[ГЛ. IV

а все следующие за ним равны нулю.

Поэтому г1 = 4 ( =

г -)- /),

а к1= 1 ( = к) (проверьте это!) — в соответствии с предложением 1°.

2.Если в матрице Т6примера 1 переставить в обратном порядк

строки, то у полученной ганкелевой матрицы

 

а

1

1

1

1

1

 

1

1

1

1

1

1

Я5 = ЛТъ

1

1

1

1

1

1

 

1

1

1

1

1

3

 

1

1

1

1

3

т

 

1

1

1

р

т

б

как легко видеть, гм = 2 ( = г + к). Проверить, что к14 = 3 ( = I)

в соответствии с предложением 2°.

3.Рассмотрим ганкелеву матрицу

0 1 0 1 0

1 0 1 •р 1

0 1 0 1 0

1 0 1 0 0

0 1 0 0 0

с (г, ^-характеристикой (2, 2) (см. упражнение 2 к § 10). В соответ­ ствии с определениями (17.17) и (17.18) для этой матрицы

 

 

Ео = 0,

Е\ =

1

0

= 1,

Ей — Ез = Е\ == 0,

т. е.

=

0 1

2.

матрицы

 

 

 

 

 

 

 

У теплицевой

 

 

0 1

0 1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

1

0 1

 

 

 

Т { = HiJs =

0 1

0 1 0

 

 

 

 

 

 

 

0 0

1

0 1

 

 

 

 

 

 

 

0 0

0

1 0

г1=

2 ( =

r j, А1=

2 ( =

A),

Z1 =

0 ( =

г —’>i) — в соответствии

с предложением 3°. В согласии с этим жещредложеиием у теплице­ вой матрицы

0 1

0

0

0

1

О 1

о о

Г Г = J e ff * = 0 1

0 1

о

1

0

1

0 1

0 1

0 10

' = 2 ( : П), к~ = 0 ( — г - гг), Г = 2 ( = к).

§ 18] ОБРАЩЕНИЕ ТЕПЛИЦЕВЫХ И ГАЫКЕЛЕВЫХ МАТРИЦ 171

4. Для ганкелевой матрицы (ср.

упражнение 3 к § 10)

 

0

4

0

1

Я . =

4

0

1

0

0 10

Vi

 

 

1 0 V* —6

найти все четыре характеристики, о которых идет речь в теореме

17.1.

Ответ,

(г, к) =

(2, 1), (V,

*А) == (3, 0);

(г1, к1, I1) =

(2, 1, 0),

(гн , fcH, Г )

= (2, 0, 1).

5. Найти все четыре характеристики (см. теорему 17.2) теплицевой матрицы

 

 

 

Т2 =

i

2

4i

 

 

 

-V a

t

2 .

 

 

 

 

 

—Ik —Va i

 

 

 

Ответ,

 

(г, к, l) =

 

(1, 0, 1), (r{, к1, ll) = (1,1, 0);

 

 

 

 

 

(r\kl) = ( 2,0), ( г ~ ,Г )= (1, 1).

 

§ 18. Обращение теплицевых и ганкелевых матриц

т.

18.1.

Общеизвестно, что

задача обращения матриц,

е.

отыскания для

задаииой

 

неособенной матрицы А

(|.4

1=54=0)

обратной

матрицы

 

А~г, является одной из

центральных и трудных задач

теории матриц. Важность

ее

решения хотя бы для отдельных классов матриц как

в теоретическом плане,

так и в прикладных задачах (ре­

шение систем линейных уравнений) не вызывает сомне­ ний. К сожалению, несмотря на обширную литературу, посвященную этому вопросу *), проблема во многих ее аспектах требует дальнейшего углубленного исследования.

*) Объем настоящей монографии не позволяет отразить резуль­ таты, относящиеся к численным методам обращения теплицевых и более общих (так называемых блочно-теплицевих) матриц. Эти ре­ зультаты носят по большей части прикладной характер (построение алгоритмов для обращения матриц на ЭВМ) и в большинстве случаев связаны с теорией фильтрации и экстраполяции скалярных и вектор­ ных стационарных случайных процессов, корреляционные матрицы которых как раз и являются тешшцевыми (в векторном случае — блочно-теплицевыми). Для ориентировки читателя укажем здесь на работы Н. Левинсона [53], П. Виттла [55, 56]. Дж. Дурбина [42а], В. Тренча [54], С. Зохара [57], Л. М. Кутикова [34а] и совсем недав­ но появившуюся статью X . Акаике [39а]. Заметим, что во всех этих

172

ПРЕОБРАЗОВАНИЯ МАТРИЦ

[ГЛ. rv

 

В этом параграфе

мы изложим

сперва недавно полу­

ченные

результаты

И. Ц. Гохберга и

А. А. Семенцу-

ла

[17],

а также И. Ц. Гохберга и Н. Я.

Крупника [16]

по

обращению т е п

л и ц е в ы х

матриц Тп (нам удоб­

нее здесь, следуя упомянутым авторам, работать с матри­ цами порядка п -Г 1). Затем используем эти результаты для решения с помощью методов § 17 задачи обраще­ ния г а н к е л е в ы х матриц Нп.

18.2.Центральной в излагаемой ниже теории является

Т е о р е м а 18.1

( Г о х б е р г а и С е м е н ц у -

л а). Если

гпеплицева матрица

Тп = |Ср-qЦ»,q=o такова,

что каждая из систем уравнений

 

 

11

 

 

 

 

)

2 сл -Л =

6ро

(р =

о, 1,.. .,п),

(18.1)

п

q=G

 

 

 

 

 

 

 

 

 

2

ср-чУ<1-п=

 

(р =

0, 1 , . . п)

(18.2)

 

<1=0

 

 

 

 

 

 

 

 

 

разрешима и выполнено условие х0 =j= 0,

то матрица Тп

неособенная и

обратная к ней

матрица Тй1

строится

по формуле *)

 

 

 

 

 

 

 

 

 

'

х о

0

. . 0

Уо

У- 1 ■•• У-п

 

 

гр-1

_

x i

 

.. . 0

0 У*

■■ ■ У -п +1

 

1 П

х~0х

 

 

 

 

 

 

 

 

 

 

х п

r« - i

■• *0

0

0

■■ Уо

 

 

V

 

0

0

0

 

0

0

0 * „

*,>-1

• • •

 

 

жг

 

У-п

0

0

 

0

0

о о

 

. . .

х „

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У -п и

 

0

0

0 0

0

 

 

 

 

 

.

 

 

 

■ ■ ■ х п

 

У-г

У -2 У - 3

 

У -п 0

0 0

0

. . . 0

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

цевых матриц вида [|Tp_q

где T j (/' = 0, ±

1,

...,

+ (п — 1))

— произвольные

квадратные

матрицы,— обратимости

всех

«усе­

ченных» матриц

|2’„ _ ( ||pi7=a

( k = 0, 1, ..., н

1).

 

 

 

*) Для н е о с о б е н н о й матрицы Тп числа

(х 0,

ад,

..., х п )

и (у», У У-2» ..о

У~п)>как видно пз формул (18.1)

и (18.2), суть эле- ■

§ J8] ОБРАЩЕНИЕ ТЕПЛИЦЕВЫХ И ГАНКЕЛЕВЫХ МАТРИЦ 173

Д о к а з а т е л ь с т в о . При п = 0 утверждения теоремы тривиальны (заметим лишь, что в формуле (18.3) вычитаемое в фигурных скобках в этом случае отсут­ ствует). Итак, пусть п^> 0. Неособенность (обратимость) матрицы Тп докажем от противного.

Пусть

|Тп |= det Тп = 0. Тогда

строки Гр =

= (СР. cp-v

ср- п) матрицы Тп (р = 0, 1, 2,

..., 7г) линейно

зависимы. Поскольку система (18.1) разрешима, а в пер­ вом из ее уравнений (при р = 0) правая часть равна еди­ нице, первая строка Г0 матрицы Тп не может состоять только из нулей: Г 0 Ф 0. Отсюда и из линейной зависи­ мости строк вытекает, что хотя бы одна из остальных строк Г1( Г2, . . ., Г„ выражается линейно через п р е д ы- д у щ и е. Обозначим через m ( > 1) н а и б о л ь ш и й из номеров, для которых

где а э, ах, . . ., a m_x — некоторые

комплексные числа.

Покажем, что m = п. Пусть тп <

«• Поскольку

менты соответственно первого п последнего столбцов обратной мат­ рицы Г”1. Стало быть, в теореме 18.1 релыщето восста новлепил всей

матрицы Г" 1 по этим двум ее столбцам. Тот факт, что матрица Г"1

определяется этими ее столбцами, по-видимому, впервые был уста­ новлен (трансцендентными методами) Г. Бэкстером п И. Хиршманом

[40], которые получили п явную формулу (отличную от (18.3)) для восстановления матрицы Г"1 по этим столбцам. Одпако в их работе

па числа (.тл, ад, ..., хп) н (уо, у- ь

...,

у_п) налагалось дополнитель­

ное ограничение: многочлены

 

 

71

 

 

*(?)= 2 х&’ и

у

= 2 у- £ }

пе должны обращаться в пуль в единичном круге | £

|<1 .

При тех

же ограничениях формулу (18.3) впервые обнаружил

А. А.

Сомен-

цул [38], результаты которого легли в основу соответствующего раз­ дела книги И. Ц. Гохбсрга и И. А. Фельдмана [5а] (гл. III, § 6). В приводимой в тексте редакции теорема 18.1 и ее доказательство заимствованы ив [17]. '

174

ПРЕОБРАЗОВАНИЯ

МАТРИЦ

[ГЛ. IV

то,

произведя замену индексов

р = р' — 1, q =

q' — 1

и возвратясь затем к прежним их обозначениям, получим

т

 

 

c m + l - q 2 a P - l c P ~ q

(? — 1> 2, . . ., п).

(18.4)

Р=1

Всилу разрешимости системы уравнений (18.1) для вся­

кого

ее решения (х0, хх, . . ., хп)

имеем, учитывая, что

2

т + 1

п:

 

 

о — бщ+i, о — 2

еm+1 -g Х п

—“

 

 

<2=0

 

 

 

— 2 СП1+1-<2 —9

2 а Р -1 + 0

 

9 = 0

 

Р= 1

 

тп

 

2

<'П1+1-<? —g

2 ®p-i 2

cp-q—9 —

 

9 = 0

 

 

Р = 1

9 = 0

 

 

 

 

 

III

 

7L

 

 

 

IU

=

С-т+1

2

ap-icp) - 0 + 2

( с,п+1-ч — 2 ар- 1 Ср- 9 I —9 •

 

 

Р = 1

'

9 = 1 '

 

-

Р = 1

В силу

равенств (18.4) и неравенства

х0

ф 0 отсюда сле­

дует, что

 

 

ТП

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cm+l — 2 а Р -1 с Р ’

 

(18.5)

 

 

 

 

р=1

 

 

 

 

Объединяя теперь (18.4)

и (18.5),

получаем

 

 

 

 

т

 

 

 

 

I\ n + l — 2 ® Р -1 Г р ,

Р= 1

аэто противоречит выбору числа т.

Итак, т = п и, следовательно, последняя строка Гп матрицы Тп линейно выражается через предыдущие. Но это в свою очередь противоречит разрешимости системы уравнений (18.2), в которой все правые части, кроме по­ следней, равны нулю, а последняя равна единице.

Таким образом, мы доказали, что матрица Тп обра­ тима. Отсюда вытекает, что решения {ж0, xv . . ., £п} и

§ 18] ОБРАЩЕНИЕ ТЕНЛИЦЕВЫХ И ГАНКЕЛЕВЫХ МАТРИЦ 175

{у-п, у-п+1 , • .

уо} систем

(18.1) и (18.2) соответственно

определяются единственным образом. В частности,

 

х 0 = у 0 =

1^ ~ 1[1

(18.6)

где |Гп_]l |— определитель

усеченной матрицы

Гп_х =

~ 1!ср-д 1р.'т=о-

через в = « ъл R fc=0 матрицу из

правой

Обозначим

части равенства (18.3). Учитывая правила умножения мат­

риц, легко проверить,

что элементы bjk вычисляются по

формуле

 

m in (з, к)

 

bjk = Xq1 2

- y^n. 1+j_sxn+1+s_k) (/, к - 0, 1 ,..., в),

8=0

(18.7)

 

причем для единообразия записей здесь принято хп+1 =

=У -п -1 = 0 .

Преобразуем теперь это выражение для bjk следующим образом. Выделим из суммы в правой части слагаемые, отвечающие значению индекса суммирования s — 0, а в оставшейся сумме сделаем замену индекса s = s' — 1. Тогда получим

b j k = Х 0 ( Х ] У ^ к JZ-n-l+j^n+Hc) 6/-Х, ft-i

0', к = 1, 2, . . ., гг). (18.8)

Кроме того, полагая в формуле (18.7) для bjkиндекс j = 0

=

0), имеем, с учетом равенства (18.6),

 

 

Ьок — У-hi Ьк0 =

хк

= 0, 1,

2, . . ., гг).

(18.9)

Полагая там же

/ ==

гг, находим

 

 

 

к

 

 

 

 

 

Ьпк

2 (Хп-ВУa-к

У-l-Ai+1+s-lc)

 

 

 

8=0

 

 

 

 

 

=

Х 0 ( х п У - к

У—х^-тг+Х—ft “Ь Х п - г У 1 - h

У - 2 х п + 2 - к

“Ь

•■ •"Н ^'n-ft+2j/-2

y~h+lxn-l Н-

У-hXn “Ь

“Ь х п- куо У-i-hXn+i) = хп~к ( к — 0, 1, 2, . . ., гг),

176

ПРЕОБРАЗОВАНИЯ МАТРИЦ

 

1ГЛ.

J V

где исполъзоганы

равенства

(18.6) и я:,1+1 =

0. Этот

ре­

зультат п

аналогичный ему (получающийся

при

к — п)

объединим в виде формул

 

 

 

 

Кк =

хп-к*

bhn = yh. n

= 0, 1, . .

п).

(18.10)

Нам нужно проверить, что Т„В = Е. Покажем сперва, что матрица

А — ТпВ — |Лр,, |р, (/=о

теплицева. В самом деле, в силу равенства (18.8) при

р, q = 1, 2, . . ., п имеем

П

 

 

Ярд / 1

^

 

ь—0

п

 

 

 

^ ср^07

2l Ср-$ l^s-1, 7-1 ^0 (ЗД_7 У

я

 

s = l

 

Преобразуя это выражение с учетом формул (18.6), (18.9) и

(18.10),

получаем

 

 

 

 

 

 

 

 

71

 

 

 

 

 

 

— Я0

 

-f- ^ C/j-s-1^>,7-1

"

 

 

 

 

 

 

3 = 0

 

n

 

 

 

 

 

 

 

 

 

 

 

Д’о

q-1 Сp-l-пУ0 4" х 0 bQl 2

Cp-S'X's

 

 

 

 

 

 

s = l

 

71—1

 

 

 

 

 

 

 

 

 

 

 

 

-Т'о

Ьп, 7 -1

2

^P-l-sl/s-n •

 

 

 

 

 

 

 

s = 0

 

Объединяя здесь первое слагаемое с четвертым,

а третье

с пятым,

имеем

 

 

 

 

 

 

п

 

п

 

 

п

 

 

& Р 7 ” 2

^ p - l - S ^ S , q - l 4 “ ^ 0

^ 0 7 2

t ' P - S ^ ' S

 

^ 0 2

^ ( p - l ) - s J /s -7 1 J

3=0

 

6 = 0

 

 

 

 

 

т. e. учитывая

(18.1) и (18.2),

 

 

 

 

 

 

G p q =

^ P - l j Q - l

(.P* 4.

= l i

2 ,

 

ft ) .

 

А это и означает, что матрица А тецлицева.

5 IS) ОБРАЩЕНИЕ ТЕПЛИЦЕ ВЫX И ГАНКЕЛЕВЫХ МАТРИЦ 177

Остается заметить, что в силу равенств (18.9) и (18.10) в сочетании с (18.1) и (18.2) соответственно

 

п

 

п

 

 

Я р о “

2

C p - s ^ s d =

2

=

б р П,

 

8 = 0

 

s = 0

 

 

 

п

 

п

 

 

Я р п =

2 Я р _ вЬ 5П =

2 C p - s l l s - n

— 5 р п

 

s = 0

 

s =0

 

 

т. е. первый и последний столбцы матрицы А совпадают

ссоответствующими столбцами единичной матрицы Е =

||6pgj|p, <2=0. Отсюда в силу теплицевой структуры матрп~ цы А следует полное совпадение этих двух матриц: А = Е-

Теорема доказана.

Условие х0 0 в теореме 18.1 с у щ е с т в е н н о — см. по этому поводу упражнения 3—6 в конце настоящего

-параграфа.

18.3.Легко видеть (см. (18.6)), что при выполнении условий теоремы 18.1 и «усеченная» матрица

 

 

 

 

Т п - г =

II сд - Jp ,q = 0

 

 

 

обратима. Более того, оказывается,

что

 

обратная

к ней

матрица

Тй\ также

может быть построена по решениям

систем уравнений (18.1)

и (18.2).

 

С е м е н ц у л а ) .

Т е о р е м а

18.2

( Г о х б е р г а и

Если

существуют

решения

{х0,

xv

х2, . . .,

х п) и

{у-п,

У-п+1, . . .,

уо)

систем (18.1) и (18.2) и выполняется

условие

х0 ^ 0,

то матрица

Гп_1

=

|ср_9 |рГч==о обра­

тима и обратная к ней строится по формуле.

 

 

 

х 0

0

. .

0

У о

У - 1

У - п +1

 

 

 

X I

Х а

. .

0

0 У о

У —1 1 + 2

 

х п -

1 Х п -2

 

.

х 0

 

 

 

 

У - П

0

. . .

 

0

У - п + 1

У - п

■■ •

.

 

0

 

 

 

0

0

Xп *71-1

' . .

0

хп

. .

У О

Х\

Х2

(18.11)

У - 1

У - 2 ' • У -

0 0

• а;п ,

178

 

 

ПРЕОБРАЗОВАНИЯ МАТРИЦ

 

 

[ГЛ. IV

Д о к а з а т е л ь с т в о .

Обозначим

через

D =

==||^й1л-=о матрицу

из

правой части равенства

(18.11).

Непосредственное

вычисление дает (ср.

(18.7))

 

 

 

 

 

 

min (j, к )

 

 

 

 

 

 

 

 

d j k

= Х 0

 

2

( Z j - s l / s - k

У - n + j - a X n + s - k )

(18.12)

 

 

 

 

6=0

 

 

 

 

 

 

 

для всех /, к = 0,

1,

. .

п — 1.

Отсюда, учитывая ра­

венства

хй =

у0 (см.

(18.6)),

находим

 

 

 

 

&к0

Хк

Хд ХпУк-П1 dofr

i/_(f

Хд У-пХп-к

 

 

 

 

 

 

 

 

 

(к =

0, 1,

. . .,

п -

1),

(18.13)

djk ~ dj-ij h_i “Ь Хд (Х}У^х

У-п+]хп-к)

 

 

 

 

 

 

 

 

 

 

и, к = 1,

2,. .

п -

1)

(18.14)

(ср.

аналогичное] вычисление для bjk в доказательстве тео­

ремы 18.1).

 

 

 

 

 

 

 

 

 

для

Далее из равенства (18.12) при j = п — 1 имеем

к =

0,1,

. . .,

п

1

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

dn-1, к =

Х0 2

{Хп-1-sya-к

У-1-rTn+s-k)

 

 

 

 

 

 

6=0

 

 

 

 

 

 

 

 

 

 

=

Хд1 (Хп-! У-к —

y_xXn_k +

Xn-a2/i-k — УгХп+г-к +

• . •

 

 

 

•••+ Хп-кУ~1 У-кХп-1 + Xn-^цУд У1-кхп) —

— Xn-k-i *01ХпУ-к-1-

Выпишем этот результат вместе с аналогичным ему, по­ лучающимся из (18.12) при к = п — 1, отдельно:

dji-i, к = xn-k-i

Хд хпу-к^ ‘, dkfn_i =

2/-п+л+1

Хд У-nxk+i

 

= 0, 1,

..., п -

1). (18.15)

Убедимся теперь, что первый и последний столбцы матрицы Т n_iZ> = 1 ё,к|[£~;iLo совпадают с соответствую­ щими столбцами единичной матрицы Е = |8jk ||£"к=0.

В самом деле, учитывая (18.13), (18.1), (18.2) и (18.6),

§ 181 ОБРАЩЕНИЕ ТЕПЛИЦЕВЫХ И ГАНКЕЛЕВЫХ МАТРИЦ 179

имеем

 

 

 

 

П-1

 

 

 

п—1

ej0 2

fy-ic dko 2

fy-k'X'k

^0

2 Cj-кУк-п

k=0

k=0

 

k=0

n

 

 

 

 

— 2 ^i~kx k cj-nxn x0

xtl. 2 ^j—кУк-п 4 ~ x0 хп^]—пУо—

k=o

 

 

k=a

 

= Й,'n

X n X n & in =

Й11

 

 

JJ0

•b0 x nujn —

u70

 

 

(/ = О, 1, .. . , n — 1).

Аналогично (с привлечением

еще формул

(18.15))

по­

лучаем

n—1

 

 

n —1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n-i —

2

ci-k

n-i — 2

ci-k (У-n+k+i

y^nxk+i)

 

 

 

ft=o

 

 

Jc=o

n

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

== 2

-с}+1-кУк-п

xo У-п 2

ci+i-kxk ~

 

 

 

 

 

k=l

 

 

 

k=l

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

— 2 ^J+l-кУк—п

 

x o У-п 2 Cf+i—k^t4~Д'О У-ПС]+ 0

 

 

k = 0

 

 

 

 

k=0

 

 

 

 

 

fij+l,n —

У-n^j+l, о —

n -i

( / — 0, 1) . .

П

1).

Для завершения доказательства теоремы осталось по­

казать,

что

матрица

1 ejk ll£"k=0 теплицева.

Но из (18.14),

(18.15)

и (18.13) следует,

что

 

 

 

 

 

 

71—1

 

 

 

 

 

 

 

 

 

e j k ~

2 c j ~ s ^ s k = =

 

 

 

 

 

 

 

 

 

s=0

n —1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cj d-ok

4" 2

cJ-e [^s-l, k-l +

Ж01 (х аУ-к У-п+$хп-к)]

=

 

 

 

 

8=1

 

 

 

 

 

 

 

 

 

71—2

 

 

 

 

n—1

 

 

 

 

2 C i - 1 ~ S d S l f t _ l 4 “ C j d g h + X g У - К 2 C j - s X s

 

 

 

 

 

 

 

 

 

 

s= l

 

П—1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xo Xn—k 2

У^-п 1

 

71— 1

 

 

 

 

 

 

 

S=1

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2

 

d a , k-l

C j ^ n d n - 1 , k-l 4 “ c j

4 “

 

 

 

8=0

Соседние файлы в папке книги из ГПНТБ