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

книги / Линейная алгебра

..pdf
Скачиваний:
34
Добавлен:
12.11.2023
Размер:
13.06 Mб
Скачать

Найти матрицу сопряженного оператора р* в том же базисе.

21. Пусть даны матрица Грама Г базиса е и матрица А линейного опера­ тора в этом базисе:

Г =

(

2

-1

-1

2

 

0

-1

Найти матрицу сопряженного оператора (р* в том же базисе.

22. Построить канонические разложения с ортогональными трансформи­ рующими матрицами для следующих симметрических матриц:

 

/ 0

2

 

2 \

 

/

2 —2 —1 \

/ 2 1 - 2 4

1)

2

3 -1 1,

2)

 

- 2

 

5

2

, 3 )

1 2

 

2

,

 

\ 2 —1

 

3 /

 

\ - 1

 

2

2 /

\ —2 2 —1 /

 

/ 2

 

2 -2

\

/ 1

 

—2

2 4

/ 0

 

2

5 4

4)

 

2

5 - 4

, 5)

 

- 2

 

- 2

4

, 6)

2 - 3

 

- 2

,

 

\ — 2 —4

5 /

 

\

 

2 4 - 2 /

V, 5 —2

 

о /

 

/

3 -2

-4

4

/

 

5 -2

-2 4

/ 11

2

- 8

4

7)

- 2

6 - 2

) , 8)

 

- 2

 

6

0

, 9) (

2 2

10

 

\ -4 -2

3 /

\ - 2

0

4 /

V - 8

Ю -5

/

23 Убедиться в положительной определенности матрицы и найти квадрат-

НЫЙ корень из нее.

/

 

2 - 2

-1 4

/

13

14 -4 4

/ 3

 

2

0 4

1) -2

 

 

5

2

, 2)

14

24 -18

, 3) ( 2

 

4 -2 1 ,

\ -1 2

2 /

\

-4 -18

29 /

\ 0

 

-2

5 /

/

18 6

 

- 6 4

/

44 - 2 2

26 4

/

2

2 -2 4

4)

6

 

21 0

, 5)

-22

29 -4 ) , 6 ) (

2

5

- 4 ,

\ - 6 0

 

16 /

\

26

- 4

53 /

\ -2

 

- 4

5 /

/

 

6 2

- 2 4

/

9 -16

- 8 4

/

5

2 -4 4

7)

2

 

7 0

, 8 )

-16

33 16 , 9)

2 1

- 2

V - 2

 

0

 

5 /

V - 8

16

9 /

\ -4 - 2

5 /

24. Применяя процесс ортогонализации, привести к треугольному виду матрицы из упражнений 22, 23 и следующие матрицы:

25. Ортогональными (унитарными) преобразованиями привести матрицы из упражнений 22-24

а) к треугольному виду; б) к двухдиагональной форме.

26. Ортогональными (унитарными) преобразованиями привести к подоб­ ной матрице почти треугольного вида каждую из следующих матриц:

/

 

9 -1 5

- 3

-3 \

 

 

23

1

 

1

( - 1

4

- 1

- 1

\

- 6

7

5

 

5' У

/

1

1

1

1

\

3)

5

1

- 1

1

f

-3 2

1

1

 

 

V

2

-3 5

1

1

 

to

25

25

0

0

\

25

25

--15

--20

 

0

-15

49

7

 

0

- 2 0

7

1

/

25

25*

0

0

-25*

25

-15

- 2 0

0

-15

49

7

0

- 2 0

7

1

/

625

 

0

375t

500

'

5) 625

О

 

-600

225 - 140*

300 +105*

 

375*

225

- 140»

384 + 480»

288 +140*

 

\

500»

300

+ 105*

288 + 140*

-216 - 480*

)

27.Найти канонический вид В ортогональной матрицы А и ортогональ­

ную матрицу Т

такую, что В = Т” 1А Т

для следующих матриц А :

/

1

2

3

\

l

1

\

1

3

3

2

2

i)

2

1

 

» 2)

1

1

&

3

3

1

2

2

2

 

 

2

1

L&

 

 

 

 

-

2

0 /

 

 

3

3

'

X 2

28. Для следующих матриц построить канонические разложение с уни­ тарной трансформирующей матрицей:

 

/ 3 -* 0 \

/ 0 * 0 \

/

3 2 + 2* 0 \

1)

*

3 0 , 2)

-* 0 0

, 3) ( 2 - 2*

1 0

 

\ 0 0 4 /

\ 0 0 4 /

V

0

0 4 /

29. Убедиться в положительной определенности матрицы и найти ква­ дратный корень из нее:

 

/ 3

-*

0 \

/

3

2 —*

0 \

/ 3

*

0 \

1)

»

3 0 , 2 )

\

2 + *

7

0 , 3)

-*

3 0 ) .

 

\ 0

0

4 /

0

0

4 /

\ 0

0

9 /

30. Убедиться, что матрица является нормальной и построить для нее каноническое разложение с унитарной трансформирующей матрицей:

ч 0 O'г)(:i I -!)• з»(

31.Построить QA-разложения для матриц из упражнений 22 и 24.

32.Решить систему А Х = Ь при матрице Л, заданной ее QAразложением, и столбце Ь свободных членов:

33. Построить сингулярное разложение матриц:

7

1

2 ^

 

 

11 10

4 - 1

8

5

1 \

5

5

- 2

 

 

 

 

10

4 - 1

8

5

1

3

3

6

 

 

 

 

2

8

7

4

7

5

 

 

 

 

^ \

 

-7

 

 

1

7

2

 

- 8

-4

-7

-5 /

 

 

 

^

- 2

9

3

3 \,

 

 

/(

6

3

0

15

-3

3 \

7

7

- 1

 

 

 

 

6

3

0

9

9

-9

5

5

7

 

5)

 

2

5

4

3

3

15

з 9

1

\^ - 2

-5

-4

- з

15

3 /

з )

 

 

Указание. Одно из ненулевых характеристических чисел матриц А * А следует искать среди чисел 36, 144, 324.

34. Построить полярные разложения следующих матриц:

 

4

-1

4

- 2

 

1

1

1

1

 

1

4

2

 

1

1

1

1

1)

4

> 3)

2

2

-1

1

1

1

1

 

-2

 

 

- 2

2

1

 

1

1

1

1

 

 

 

35.Построить скелетные разложения для матриц из упражнения 33.

36.Используя скелетные разложения, построить псевдообратные ма­ трицы для следующих матриц:

2

/

1

2

3

1

2

3)

2

2

1

1

1

V

-3

-4

-4

-2

1

 

 

 

 

 

2

 

 

 

 

 

2

6)

-3

-2

-3

-5

2

 

2

 

 

 

 

 

1

■ 9)

О

-1

-2

-5

3

 

37. Используя сингулярные разложения, построить псевдообратные ма­ трицы для матриц из упражнения 33 и следующих матриц:

/

12

12

6 \

/ 18 18

9

\

1)

2

-1

-2

 

4

-2

-4

 

4

-2

-4

 

8

-4

- 8

 

 

4

-2

-4

/

\

8 -4

- 8

/

38.Для матриц из упражнения 33 построить псевдообратные матрицы методом Гревилля.

39.Методом наименьших квадратов решить системы линейных уравне­

ний:

 

XI

2

3

=

18,

1)

XI

 

+ХЗ

=

18,

2xi

2

 

=

18,

 

3xi

2

-Хз

=

36.

 

XI

2

+ 2хз

=

1,

з)

XI

2

+ХЗ

=

- 1,

s-XI

2

=

1,

 

-XI

2

=

- 1.

 

—Х

— хз

 

 

Xi

2

 

=

1,

 

XI

—Х

+ 2хз

5)

2

=

- 1,

XI

2

3

=

2 ,

 

 

2xi

2

3

=

3.

 

XI

 

+Х4

=

4,

 

 

2

+Х4

=

8 ,

 

 

Х

=

 

 

хз

4

4.

Г

*1

+X2

+ 2 хз =1,

1

2X1

-X 2

+x3

==

o,

I

- * i

+2X2

+x3

=

0 ,

н1

*1

+2X2

+X3

=:

2 .

f

 

X2

+X3

=

1,

 

Xl

4) <

 

+X3

=

0 ,

XI

-X 2

—2x3

=

- 1,

,

-XI

-X 2

=

1.

(

Xl

—X2

+X3

=

6 6 ,

 

 

Xl

 

=

-33,

 

Xl

+X2

+ 2 хз

=

-33,

' I 1. Si

+2X2

+3x3

=

6 6 .

 

Г Xl

+X2

+X3

+X4

=

2 ,

8 ) {

Xl

+X2

+Хз

+X4

=

3,

1 Xl

+X2

+X3

+X4

=

4.

40.Для систем из предыдущего упражнения найти нормальные решения, решая соответствующие одну или две системы с невырожденными матри­ цами.

41.Для систем из упражнения 39 найти общие и нормальные псевдореше­ ния, решая эти системы в матричном виде с применением псевдообратных матриц.

42. Найти общее и нормальное псевдорешения системы А Х — b и их проекции Х к , к = 1, г на подпространства правых сингулярных векторов при матрице А , заданной ее сингулярным разложением, и столбце b сво­ бодных членов:

 

 

1

1

1

1

а) А =|

I

1

1

- 1

- 1

2

1

- 1

1 - 1

 

 

1

- 1

- 1

1

6 = (216, 0, 0, —216)т,

<Г1= 18, <Т2

= 12,

<71 := <72 = 6,

II ь

<73 ==0;

2)

\ — 18, <72 = 6 , <7з= 0 ;

О

<тЛ

 

3

0

0

0 \ /

<Т\

0

0

<0 ^ 4

0

1

2

2 11

0

<Т2

0

 

0

2

1

-2

0

0

с73

 

0

2

-2

1/ V 0

0

0

!, 0, о, -—162)т,

II

<71 =

С»

,

^2

18,

 

<*2

II II

3

II

II b

о

0;

2) <71 =

43. Методом регуляризации найти нормальные решения систем из упр.

39 и следующих систем линейных уравнений:

 

XI

+32

=

з,

Г

XI

+2X2

=

3,

 

 

1)

XI

 

=

 

I

2xi

- Х 2

=

1,

 

 

- Х 2

1,

2)

 

 

 

XI

+2X2

=

5.

1 *1

+S2

=

2.

 

 

 

XI

+Х2

=

 

Г

XI

 

+хз

=

0,

з)

2xi

- Х 2

=

з,

I

 

+8X2

+6х3

=

 

4)

- * 1

1,

 

2xi

+Z2

=

5.

1

XI

-7X2

—5хз

=

-1.

44. Пользуясь нормой ||А||оо, найти число обусловленностей коо матриц:

(

5

-3,31

\

2 .

е ф 0.

* \

6

-9,97

у ’

П - 1 3 0 7

45.

Оценить возможное изменение решений систем

 

 

 

 

 

Г xi

2x2

==

1 + 6 ,

Г 5xi

—3,31x2

=

1,69

+ 6 ,

 

2xi

+ 4 , 01x2

=

2 + 6 ,

у

6xi

3 , 97x2

=

2,07

+ 6 ,

при изменении е

и ei

в пределах 0< е

< 0.01, 0 < 6i < 0.03.

 

 

46.

В эмпирической формуле b = xiai + хгаг + хзаз найти коэффициенты

xi, Х2, хз по результатам наблюдений:

 

 

 

 

 

 

 

ai

0>2

аз

ь

 

ai

а2

аз

Ь

 

 

 

 

1

0

0

0

 

1

1

-1

12

 

 

 

1.

0

1

0

0

2.

-1

1

1

12

 

 

 

 

0

0

1

0

 

1

-1

1

12

 

 

 

 

1

1

1

4

 

1

1

-1

12

 

 

 

 

ai

2

аз

6

 

ai

2

аз

Ь

 

 

 

 

 

а

 

 

а

9

 

 

 

3.

1

-1

0

12

4.

1

-1

0

 

 

 

2

-1

1

-12

-1

2

1

0

 

 

 

 

0

1

1

12

 

2

-3

-1

-9

 

 

 

 

1

-1

1

-24

 

0

1

1

0

 

 

 

47.

В эмпирической формуле p (t) = хо + x\ t + хгt2 найти коэффициенты

xi, Х2, хз по результатам наблюдений:

 

 

 

 

 

 

1.

t

0

1

2

3

P (i)

1

4

9

16

 

3.

t

-1

0

1

to

p (t)

5

1

-1

-1

 

5.

t

-1

0

1

2

p it)

4

1

0

1

 

t

0

-1

-2

-3

 

2.

1

0

1

4

 

_ еШ _

 

t

0

 

1

2

3

4.

1.1

 

3.9

9.1

16.1

p it)

 

t

0

 

-1

-2

-3

6.

1.1

-0.1

1.1

3.9

p(f)

 

 

 

 

 

48. Найти многочлен второй степени p (t) = хо + x\t + хгt2, приближаю­ щий с наименьшей квадратичной погрешностью функцию /(t), заданную таблицей:

1.

t

0

1

2

3

t

-1 0

1

2

1

3

3

1

2 .

4 2

2

4

 

49. Найти многочлен третьей степени p (t) = хо + x\t + хг<2 + хз*3, при­ ближающий с наименьшей квадратичной погрешностью функцию /(<), за­ данную таблицей:

1.

t

0

l

2

<

-1

0

1

m

l

l

-1

2 .

4

2

-2

 

 

Глава 7

Квадратичные формы

7.1.Определение квадратичной ф орм ы

Действительной квадратичной ф орм ой о т п действитель­

ных переменны х ж2, хп называют многочлен с действи­ тельными коэффициентами а*у, каждый член aijXiXj которого имеет вторую степень, т.е. многочлен вида

 

п

п

 

/( * 1 , Х2,

хп) =

ai)xixj

(7.1)

<=1 j =1

Переменные жх, ж2, •••>хп можно рассматривать как координаты вектора х = (жх, ж2, •••> £п)Т в некотором фиксированном базисе действительного n-мерного пространства Х п. Поэтому квадратич­ ную форму (7.1) можно рассматривать как числовую функцию век­ торного аргумента х из Х п с действительными значениями.

Пусть /(ж 1 , ж2, ..., хп) - квадратичная форма (7.1). Коэффициент при ж? в ней обозначен через а,», а коэффициент при произведении XiXj - через а,-у. Будем считать, что

aij =

(7.2)

Этого всегда можно достичь, сложив члены а,уж,жу и ау,жуж,- и разде­ лив сумму на два. При такой договоренности квадратичную форму (7.1) можно записать следующим образом:

пп

/ ( Ж1 , ж2, а?п) = ^ Х ^ а‘>Ж|Жу =ацж ? + а22*2 + --- +

»=1 ;= 1

+flnn^n + 2ах2Х1Ж2 + 2ахзЖ1Жз + ... +

+ 2а1ПЖ1Жп + ... + 2ап_ 1>пжп_ 1Ж„.

Ее также можно записать в виде

пп

/(Ж1 , ®2j •••>жп) ==

У ! aijxixj =

 

»=1 J=1

17*

с матрицей А . Применим к ее переменным жх, жг, •••>жп невырожден­ ное линейное преобразование переменных

I *1 = ЯпУ1 + <7122/2+

+?1пУп,

 

(7.9)

[ Хп = Яп1У1 + ?г»2У2+

АЧппУп j

или, что то же самое,

 

х = Qy.

(7.10)

Подставив в форму (7.8) вместо каждой переменной ®i, жг, . .., х„ выражение через переменные j/i, J/2? ••-5Уп по формулам (7.10), получим квадратичную форму

9{Уи У2, ■■■, Уп) = (Qy)r A(Qy) = yr QTAQy.

(7.11)

Это означает, что матрицей квадратичной формы д(уи У2> •••, Уп) является матрица

В = QTAQ.

(7.12)

Таким образом, при линейном преобразовании переменных с матри­ цей Q квадратичная форма с матрицей А преобразуется в квадратич­ ную форму с матрицей В — QTAQ. При этом, если Q — невырожден­ ная матрица, то из соотношения (7.12) следует, что ранг матрицы В равен рангу матрицы А } т.е. ранг квадратичной формы не менлетсл при невырожденном линейном преобразовании ее переменных. В то же время определители матриц В и А связаны соотношением

\B\ = \A\-\Q\\

(7.13)

Напомним, что переход к новым переменным можно рассматри­ вать как переход к новому базису в линейном пространстве. Поэтому квадратичную форму (7.11) в новых переменных j/i, t/2>. . уп можно рассматривать как квадратичную форму (7.8) в новом базисе.

Примечание. Не следует смешивать правило QTAQ изменения матрицы А квадратичной формы при переходе от одного базиса к другому с правилом Q~lAQ (см. п. 3.4) изменения матрицы А ) рас­ сматриваемой как линейный оператор, при переходе от одного базиса к другому. Матрица QTAQ в общем случае не является матрицей того же оператора. Совпадение матриц QTAQ и Q~lAQ будет лишь в случае ортогональной матрицы Q, т.е. при QT = Q” 1, или, что

то же самое, когда переход совершается от отортонормированного базиса к ортонормированному базису.

Если квадратичная форма /( х i, Х2, хп) некоторым невыро­ жденным линейным преобразованием переменных переведена в ква­ дратичную форму <7(2/1, 2/2, •••, Уп), то формы / и д называют экви ­ валентными. Так как матрицы А и В таких форм связаны соот­ ношением В = QTAQ> где Q — невырожденная матрица, то ранги этих форм одинаковые. Эквивалентность квадратичных форм обо­ значают знаком ~ Она обладает свойствами:

1)

рефлексивности, т.е. / ~ /,

 

2)

симметричности, т.е. если /

~ <7, то д ~ /,

3)

транзитивности, т.е. если /

~ у, а д ~ Л, то / ~ Л.

Пусть квадратичная форма /(х 1 , Х2, ..., х„) некоторым невыро­ жденным преобразованием переменных переведена в квадратичную форму

я {у и 2/2, •••, Уп) = Ciyl + С2у\ + . . . + cny l ,

(7.14)

содержащую только квадраты переменных. Квадратичную форму д называют каноническим видом квадратичной ф орм ы / . Ма­ трицей формы (7.14) является диагональная матрица

( Cl

0

(7.15)

0 Сп

Поскольку ранги у форм / и д одинаковые, то число ненулевых диаго­ нальных элементов в матрице (7.15) равно рангу квадратичной фор­ мы / .

Базис пространства, в котором квадратичная форма имеет ка­ нонический вид, называют каноническим базисом эт о й квадра­ тичной ф орм ы . Очевидно, что две квадратичные формы эквива­ лентны между собой тогда и только тогда, когда они могут быть приведены к одинаковому каноническому виду.

7.4.П риведение квадратичной ф орм ы к каноническом у виду

Впредыдущем параграфе мы видели, что если в квадратичной форме подвергнуть переменные линейному преобразованию, то по­

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

Оказывается верным следующее утверждение.

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

Такое приведение квадратичной формы легко достигается мето­ дом Лагранжа. Основная идея этого метода состоит в последователь­ ном выделении в квадратичной форме полных квадратов по каждой переменной. При этом, если в квадратичной форме нет членов с ква­ дратами переменных, т.е. если в ней все коэффициенты ац = 0, то сначала добиваются, чтобы в квадратичной форме появились члены с квадратами переменных. Так, если все ац = 0, но a,j ф 0, то совер­ шив невырожденное преобразование переменных

*.• = Vi + У], Xj = у,- - у,-, хк = Ук при

к Ф i,j ;

получим, что член 2aijXiXj квадратичной формы примет вид

2<4j{yi + Vj)(yi - Vj) = 2aijyf - 2dijyj,

т. е. в квадратичной форме появятся члены с квадратами переменных у, - и уу, причем они не могут сократиться с другими членами формы, так как в каждый другой ее член входит у* при к ф i ,j . Таким образом, можно считать, что в квадратичной форме

n n

есть члены с квадратами переменных.

Пусть,

например, в форме /

есть член с квадратом переменной

x i, т.е. в

ней коэффициент ац

ф 0. Соберем в форме / все члены,

содержащие xi, и дополним их сумму до полного квадрата. Тогда получим

f{x 1 ,

X2,

xn) = (a n ij + 2a12iCiX2 + ... + 2alna:ixn) +

 

 

+ a22®2 + 2023X2X3 + ... + a

 

=

- — (011X1 + 0, 2X2 + ... + ai„x „)2 + g(x2,

x„),