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

book1989

.pdf
Скачиваний:
17
Добавлен:
10.04.2019
Размер:
19.14 Mб
Скачать

при этом для достижения той же точности е потребуется выпол­ нить большее число итераций.

Действительно,

согласно теоремам 1 и 2 число итераций при

малых £ пропорционально

где

£ =

£i =

б/Д для метода (2), (3),

i =

i2 =

(l 0/(1+ 0

Для метода (21) — (23).

Из формулы (23)

находим

 

Пусть, для определенности, 6±< б2, AI<A 2- Тогда получим

|, = 6,м „ b - j v ( l + i ) ( l + A ) > b .

Отношение числа итераций По1'(е) в методе

(2),

(3)

к числу итера­

ций л[,2) (е) в методе (21) — (23) окажется равным

 

 

П? (е)

 

 

 

/(■+£

1 + ? ] > ' ■

л?* («0

V

¥

 

 

 

 

 

Если, например,

 

/2= 0,5/,,

h2=0,5hi, то

получим

62=4б|, Дг=

= 4Д„ и, следовательно,

л!,1’ (s)/n(02) (е) л; 2,5.

 

 

 

З а м е ч а н и е . Существенного

ускорения сходимости

можно добиться пу­

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

 

 

 

Уin-у, Ук

+ А Ук+'А+ А Ук= /■

Tч(*+i)

Ук+1~~ Ук+'/,

т№+О

+ АУкту, + АУк+i —F

с переменными параметрами т^+1’, т^+1', й = 0 , 1 , . . , п— 1. Способ выбора ите­

рационных параметров и оценки погрешности подробно изложены в [30, с. 463]. Отметим лишь, что в случае модельной задачи число итераций, необходимых для

достижения заданной точности е, является величиной О 11п — .

§5. Метод матричной прогонки

1.В в ед ен и е . Матричная прогонка относится к прямым методам

решения разностных уравнений. Она применяется к уравнениям, которые можно записать в виде системы векторных уравнений

^оУо В аУ1

F 0,

Aiiji-x — Ctyt + B i y Ul

= —F it /=1.2,..., N — 1, (1)

Av*/,v-i — СыУм — F N ,

411

где У;— искомые векторы размерности М, Д —заданные

векторы,

 

Аи Bi, Ci— заданные квадратные матрицы порядка М.

 

 

Матричная прогонка представляет собой обобщение обычной

 

прогонки на случай системы векторных уравнений (1). По сравне­

 

нию с другими прямыми методами решения разностных задач ма­

 

тричная прогонка более универсальна, так как позволяет решать

 

уравнения с переменными коэффициентами и не накладывает силь­

 

ных ограничений на вид граничных условий.

Однако применение

 

матричной прогонки к решению двумерных разностных задач стал­

 

кивается

с двумя трудностями: неэкономичностью по числу дей­

 

ствий (т. е. большое время

счета)

и, главным образом, необходи­

 

мостью в больших ресурсах машинной памяти. Если же матрицы

 

Ai, Ви Cf имеют относительно невысокий порядок

(как это бывает

 

при аппроксимации систем одномерных дифференциальных урав­

 

нений), то матричная прогонка ничем не хуже обычной прогонки.

 

Прежде чем излагать алгоритм, покажем на простом прим:ере,

 

каким образом двумерную разностную задачу можно привести к

 

виду (1).

Запись разностного уравнения Пуассона в виде системы век­

2.

торных уравнений. Пусть в прямоугольнике

 

 

 

 

 

 

G= (0< х а< 4 , а =

1, 2}

 

 

 

 

с границей Г требуется найти решение уравнения Пуассона

 

 

 

 

д2и

д'2и

 

 

 

 

 

 

(2)

 

 

 

 

дх\

дх\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

удовлетворяющее условию Дирихле

 

 

 

 

 

 

 

и(хи х2) = ц(хи х2), если

(х„ хг)е Г .

(3)

 

Введем прямоугольную сетку

 

 

 

 

 

 

 

 

 

 

 

Gh=

{хц =

{x'i \

4 я)},

 

 

 

 

где x[l>=

ihlt x{f

= jh2,

i = 0,

1, ....

Nlt j — 0,

1,

..., Nt, htЛГ,= /„

 

h2N2 — l2 и заменим задачу (2),

(3)

разностной схемой

 

 

 

^ £ - 1 , /

^ У ц +

P i'- ц , /

|

V i,f - 1

^ £ / " Е ^ £ ,/ + 1

— —fit i

(4)

 

 

 

‘1

 

+

 

 

 

К

 

 

 

 

 

 

 

 

 

“2

 

 

 

 

 

i

= l , 2......N i - 1,

/ =

1 ,........2 N2-

1,

 

 

 

i/о/

Ро/,

UXxi

P^V,/»

i

 

1, 2, ...

,

1,

(5)

 

 

Уin == pro,

yiNi ==

 

Р£Л/г,

i ;=

1

, 2 , . . . ,

1.

 

 

N,

Р азностн ая

схем а

(4 ),

(5)

 

представляет

собой

систем у

линей­

ных алгебраических уравнений, в которой неизвестными являются

 

значения

yih i= l, 2, ....

—1,

/= 1 ,

2,

...,

Nt-~-1. Число неиз­

 

вестных

равно

числу

уравнений,

т.

е.

(Л\—1) (N2—1). Запишем

 

систему (4), (5) в векторном виде (1).

412

При решении системы (4), (5) матричную прогонку можно про­ водить как по индексу 1, так и по индексу /. Покажем, например, как подготовить систему (4), (5) к виду (1), удобному для приме­ нения прогонки по индексу i. Перепишем систему (4) в виде

Ус-и

( 2Ус/

(/и - 1 - 2У ц + У [.,-+1

 

Уи и !

f

А?

{ а;

А*

) +

А*

Л/’

 

i = l ,

2, . . . , N - 1 ,

/ =

1, 2, ...,

N2—1

 

и учтем граничные условия

 

 

 

 

 

 

У to == У to,

У i'.v, — p.(V,,

t =

1 , 2 , . . .

,

1 •

Тогда получим систему уравнений

 

 

 

 

 

 

Ус-1.1

/ 2Уп

- 2i/a + Ун \

У т л

г

ftp

A;

I, Л2,

Л2

j

 

h\

 

Ы

h\

«/£_!,/

(2у„

Ус,i- i ~ %Уц +

Уи + 1

\

, Уи 1, /

,

“ ч —

 

------------------- ч ----------------

 

 

)

+

ч --------------

 

 

 

/ = 2,3, . ..,

 

2,

 

 

 

Ус-i.N^i

 

Уi .Nr-г — tyj.N r-i

+

Ус+l.Nr

Pi.v,

а;

 

 

 

 

 

А?

 

 

где i = l , 2, ..., JV,—1. Далее, обозначим через £ а единичную матри­ цу порядка — 1 и через Л2— следующую трехдиагональную ма­ трицу того же порядка

2

1

0

0 •••

1

—2

1

0 •••

0

1

—2

0

(6)

 

 

 

 

Ясно, что Л2 представляет собой матрицу оператора второй раз­ ностной производной по направлению хг-

Введем для i= 1,2, ... , N,—1 векторы

Ус(t/iii Усг> • ■•» yc,Nt-i)T■>

(7)

Fc =

|7ii +

,

i»/

, • ■■, /;. 2- , А.л

+

— ■

(8)

 

 

f 2 13

V 2 / , - 1

 

Тогда предыдущую систему уравнений можно записать в вектор­ ном виде

Е аУс-1 (

— Л 2 ^ yi -1— —Eztji+i =

Fi,

h\

U

)

(9)

i = l , 2, ..., Nt— \.

Эту систему уравнений следует дополнить граничными условиями

Уа — Ро, У if , — Рлт,

413

где

P'0— (poi, Рог. •••) I^O'Nj-i)7,

p.v,—

(P N,I P JV,2 • • • Pjf|,rr2- i ) T-

 

Таким образом, разностная схема (4), (5)

записывается в век­

торном виде

(1), где В0 и

— нулевые

матрицы, A{= B i= h^2E2y

Ci — 2И?Ей—Л2, i = l , 2, ...,

Ni— 1.

 

 

 

по

Может оказаться, что N2^>MU т. е. что число точек сетки

направлению х2 гораздо больше числа

точек

по направлению я,

(например

в случае, когда

прямоугольник G сильно растянут в

направлении х2). Тогда выгоднее

пользоваться прогонкой по

ин­

дексу /, так как при этом соответствующие матричные коэффициен­ ты будут иметь порядок А,—1 гораздо меньший, чем N2—1. Соот­ ветствующая система векторных уравнений имеет вид

~ ~ Е1У1-1

( ~~ Ei

АЛ у , +

4 -

Е м +1= — FI,

hl

V hl

 

j

hi

 

j Г= 1I 2,. . N2

1,

Уд== PQ,

УjV== PJVH

где Et— единичная

матрица

порядка

—1, Ad— матрица, анало­

гичная (6) и имеющая порядок А,—1,

 

 

УI = (У\Ь Уги

•••» yNi-i,i)T,

Ро =

(Рю> Р20»• • • >Рл\—1,о)Г|

p,v2 = (PlAf,P2jVj . . . Рл',-1 ,Л'г)Г.

3. Алгоритм матричной прогонки. Пусть задана система урав­ нений (1). Формулы матричной прогонки можно получить так же, как и формулы обычной прогонки (см. п. 7 § 4 ч. I), однако при их выводе надо учитывать, что коэффициенты уравнения (1) непе­ рестановочны. Будем искать решение системы (1) в виде

£/i=ai+1Pi+i+ Pi+i, i= 0, 1, .... N— 1,

(10)

где ai+,— квадратные матрицы того же порядка М, что и порядок матриц Aj, Ви Ci, а р,+1— вектор размерности М. Подставляя у;= = ai+1(/i+I-t-pi+i и yi-i= aiyi + $i=aia,i+lyi+i+ (аф,+1 + Р0 в уравнение

Aiyi-i—Ciyi+Biyi+i= —Fi,

получаем, что это уравнение будет выполнено, если потребовать

(Ада,—С;)ai+1 + Bi+i= 0,

 

(/4;CC;--Ci) Pi+1 = --- (AiPi+fj) .

 

Отсюда приходим к следующим рекуррентным

соотношениям

для определения матриц cci+1 и векторов jii+1:

 

ai+1 =

(Ci — A a£)-1 В,,

(11)

Pi+i =

(Ci Aicci) 1 (.4;Pi -j- F;).

(12)

Здесь (= 1, 2, ..., N—1. Начальные значения a,i и p, задаются в соответствии с уравнением

СоУоА~ ВцУi = —Fц.

414

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

 

 

 

 

 

Уо = Са Вау1-|- С0

F0.

(13)

Сопоставляя (13) с уравнением (10)

при t= 0, получаем

 

 

 

 

 

 

 

 

 

(14)

После того как все коэффициенты

а„ j}f найдены, векторы yit

i = N—1, (V—2,

1,

0,

определяются

последовательно из

урав­

нения (10), начиная с

 

 

Для начала счета надо знать вектор yN,

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

 

АкУп- i

CNyN=

FN,

Ун-i = {Х,хУн+ $N-

 

Отсюда получаем

 

 

 

 

 

 

 

 

 

Ук=={Ся

AN<XN)_1 (-Лифвг+ Fн) •

(15)

Объединяя формулы (10) —(12), (14), (15), приходим к следую­

щему алгоритму матричной прогонки для системы (1):

 

ai+1 — (Ct —Д а£)_1В£)

t =

1, 2, . . . , N — 1, а ^С ^В ^,

(16)

IW = (Ct -

Atat)-1(Л£р£+

Ft),

i =

l,2, . . N, P1=C71F0,

(17)

y i= a,+it/i+i+ Pi+i,

i= N —1, N—2, ..., 1,0, Уп=$w+i-

(18)

При реализации метода матричной прогонки приходится запо­ минать все матрицы оя, i = l , 2, ..., N—1, что ведет в случае матриц больших размеров к необходимости использования внеш­ ней памяти ЭВМ и тем самым к увеличению времени счета.

Кроме того, реализация формул (16) сама по себе требует боль­ шого числа действий. В каждой точке i приходится один раз обра­ тить матрицу и сделать два умножения матриц порядка М, что требует О(АН) арифметических действий. Следовательно, для вы­ числения всех коэффициентов со, i = l , 2, ..., N—1, требуется 0(M 3N) действий. Для модельной задачи, когда M = N = h ~ \ число действий становится величиной О (Л-4). По указанным причинам (большой объем памяти и значительное число арифметических дей­ ствий) матричную прогонку сравнительно редко применяют для решения задач математической физики. Однако в тех случаях, когда матрицы Аи Ви Ct невысокого порядка (небольшое число то­ чек по направлению х2), необходимый объем памяти и число дей­ ствий резко сокращаются и метод можно рекомендовать для прак­ тического использования.

4. Устойчивость матричной прогонки. Так же, как и в случае обычной прогонки, возникает вопрос о численной устойчивости ме­ тода матричной прогонки. Получим здесь достаточные условия устойчивости в виде требований, предъявляемых к матрицам А{, Bit Си i= 0 , 1, ..., N.

Пусть в системе (1) yt и /\— векторы размерности М, Ait Ви Ct— квадратные матрицы порядка М (векторы и матрицы могут быть как вещественными, так и комплексными). Будем рассматривать

415

матрицы Аи Ви Ci как линейные операторы, действующие в М-мер- ном линейном пространстве Н (вещественном или комплексном). Предположим, что в Н определены нормы вектора [Ml и подчинен­ ная ей норма матрицы. При доказательстве устойчивости прогонки нам потребуется следующее известное утверждение.

Л е м м а 1. Если для данной матрицы А существует константа Ч>0 такая, что для любого х ^ Н выполнено неравенство

\\Ах\\^\\х\\,

т>0,

(19)

то матрица А имеет обратную, причем

все собственные

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

сначала, что

числа матрицы А отличны от нуля и, следовательно, существует А~‘. Пусть Я — любое собственное число матрицы А и z — отвечаю­ щий ему собственный вектор, т. е. Az="kz. Согласно условию (19) имеем

WAz\\=\K\\\z\\^\\z\\,

т. е. |Я| ^ у > 0 , и тем самым Я=Д0.

Таким образом, матрица А имеет обратную. Пусть у ^ Н — лю­ бой вектор. Обозначая х=А~'у, получим из условия (19), что

Следовательно,

||Л-‘||

что и требовалось.

Метод прогонки (16) —(18)

будем

называть устойчивым, если

матрицы С,—Л,а, имеют обратные и ||а,||^1 , 1=1, 2, ..., N.

Из устойчивости прогонки следует

однозначная разрешимость

системы (1). Действительно, в этом случае, исходя из рекуррентных формул (18), можно представить решение задачи (1) в явной фор­ ме в виде конечной суммы с коэффициентами, зависящими от он,

Условия Haill^l обеспечивают численную устойчивость счета по формуле (18). Нарушение этих условий не всегда приводит к силь­ ному накоплению погрешности. Однако подробный анализ вычис­ лительной погрешности метода прогонки выходит за рамки данной книги.

Сформулируем теперь теорему об устойчивости матричной про­

гонки.

 

1. Пусть

А{, В,— ненулевые

матрицы,

1 = 1, 2, ...

Т е о р е м а

..., 1V—1 , и пусть существуют матрицы CJ1, 1=0, 1 , ..., N. Если

выполнены неравенства

 

 

 

 

 

 

 

 

 

|CrMil + ]C7lB i |< l f

г = 1, 2,..., iV — 1,

(20)

 

 

||Со1В0|< 1 ,

IIOvVLv [|<

1,

 

(21)

то матричная прогонка устойчива.

 

по

индукции,

что

|!а*1!^1 и

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

Докажем

матрицы

Ci

имеют

обратные,

1=

1, 2,

...,

N. Неравенство

llctill^l

выполнено в силу первого из условий

(21). Предположим,

что Hctill^l для некоторого 1^1. Докажем, что тогда

(С,—Д-аО-1

существует и

||ai+1||^ l .

Поскольку

С;—Aiai= C i{E—С г’Дсо), до­

статочно доказать существование матрицы

— СГМ^а;)-1.

416

 

 

 

 

 

 

 

 

 

 

Пусть i e / / — любой вектор. Тогда

|1(£ — C l l Ai<Xi) х || >

||х ! — j| C 7 1A i(x ix jl

>

 

 

 

 

 

>

1 х | -1 CTMi 11 a, S1 x 1> (1 -

I) C7lA D И .

Отсюда и из условий 1 —\\C'i1Ai\\'^\\C7lBi\\

(см. (20))

получим

1(£ — C l 1 AiCt?) х I >

7,-1x ||,

i = l , 2 .......N — 1,

(22)

где ■Yj= ||C rlfii||>0. Неравенство

^,>0

следует из того,

что С71 —

невырожденная

матрица и

В{ф 0, и

поэтому C ^B t— ненулевая

матрица. Из неравенств (22) и леммы

1

следует существование

матриц, обратных к С,—Л,сс„ г =

1, 2, ..., N—1, и оценки

 

 

I K E - c r M

^ n ^ i c ^ r 1.

 

(23)

Таким образом, матрицы ос1+1, заданные рекуррентным соотно­ шением (16), существуют. Перепишем выражение для ai+l в виде

ai+i = (£ — CrMiGti)-1 (СГ^О-

Отсюда и из оценки (23) получим,что

I a,4i || ||- CrMia,-)-1 IIII t f B i К 1.

Итак, ио индукции доказано, что jlaill^l, г= 1, 2, ..., N. Для завершения доказательства теоремы 1 осталось убедиться в том, что существует матрица, обратная к С.ч-—Л^а». Поскольку ||aNi l ^ l , получим, как и ранее, что

I! СмАман) *Ц (1 —1 СмЛд- |) ;| х ||

для любого х еЯ . Следовательно, неравенство (22) выполняется и

при i— N с константой 7^=1 — ЦСлМ^Ц. Неравенство ^ > 0 выпол­ нено в силу второго из условий (21). Теорема 1 доказана.

З а м е ч а н и е .

Матричная

прогонка

будет

устойчивой

и в том

случае, если

вместо (21) потребовать

 

 

 

 

 

 

 

 

 

 

II C ^ fio lK 1,

 

 

1.

 

 

(24)

 

Доказательство проводится так же, как ив

теореме 1. Надо

заметить только,

что

в случае (24)

выполняются строгие

неравенства Ц0Si] | С

1,

i = l ,

2..........IV, и

 

II

( £ — C ^ A

f j

a . f j ) х ;] О

( 1 — [ ] С у 1 A

v a jV ||) j| х

|j ,

 

 

причем

 

 

 

 

 

 

 

 

 

 

1 - « C ^ M , v a v Ц >

1 _ | l a v [ЩС^Лд, || >

1 - | C ^ . 4 jV Ц > 0 ,

т. е.

1 — ЦС^Л^ад, 1> 0.

 

 

 

 

 

 

 

 

Применим теорему 1 к исследованию устойчивости метода про­ гонки для разностного уравнения Пуассона (см. п. 2). В этом

417

случае система (1) принимает вид (9), причем

 

 

 

Ac = Bi = th2E2,

С1— 2К?Ег— къ,

 

 

i— 1, 2, ..., Лл—1,

Во—Л^=О,

 

 

 

где матрица Л2 определена согласно (6).

 

 

 

 

Условия устойчивости прогонки (20) принимают вид

 

 

||СГ1К 0 ,5 ^ ,

 

1 = 1 , 2 , . . . , ^ - 1

,

 

 

и будут выполнены, если

 

 

 

 

 

 

 

 

\ С с У \ > М у \

 

 

 

(25)

 

 

 

h\

 

 

 

 

 

для любого вектора у размерности Nz—1.

 

 

 

 

Выберем в качестве нормы вектора

 

 

 

 

 

У={Уи Ун,.. .,Унг-i)T

 

 

 

величину ||г/|| — Цу, у), где

 

Л 'а-1

 

 

 

 

 

 

(и, и)=

 

 

 

 

 

 

2 hjUflj.

 

 

 

 

Тогда

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

I Ciyf = ( - \ y ~ k ^ ,

4" У — Aai/ ] =

 

 

 

 

 

l hl

hl

 

J

 

 

 

 

 

 

= 4

fi у ii2—4

y ) +

ii Л.Уf >

hi

i у i2,

 

h[

 

 

h\

 

 

 

поскольку

 

 

 

 

 

 

 

 

 

(А2г/, у) =

— 2

(y7J

2h< 0.

 

 

 

 

 

 

/'=i

 

 

 

 

 

Тем самым условие (25) выполнено и матричная прогонка для системы (4) —(5) устойчива.

§6. Метод редукции

1.Вывод основных формул. Метод редукции является прямым

методом решения системы разностных уравнений, имеющих вид

У -i—Q/i+1/i+i = — j'=l, 2, . . . , N— 1,

(1)

£/о== М-i,

1/JV= (J-21

 

(2)

где у, — искомые векторы размерности М,

pt, р2— заданные век­

торы и С — заданная квадратная

матрица

порядка М. Принципи­

альным отличием данной системы от системы (1) из § 5 является независимость матрицы С от индекса i и равенство коэффициентов

при

и yi+l.

418

 

В основе метода редукции лежит специальный способ исключе­ ния неизвестных из системы (1). Запишем уравнение (1) в точках £—1 и £+1, т. е.

У1 Q/i-l+i/i = Fi-i 1 J/i Cy.i+1-р£/i+2= Fi+!,

и сложим эти уравнения. Тогда получим

£ / t - 2 - ) ~ 2 i / i — С( г / i —! —|—r /i + i ) +

t / i + 2 = —

(FM + f i + i ) .

 

откуда, учитывая, что

 

 

 

1Л—1Ч- Уг+1

1 Т,,

 

 

придем к уравнению

 

 

 

г/;_2—(C2—2E)yi+yi+2= — (Fi^i+ CFi+Fi+i),

(3)

связывающему значения искомого вектора в узлах одинаковой чет­ ности. В частности, если £ — четные, то проведено исключение не­ четных узлов. Далее этот процесс исключения можно продолжить

к 0 0~Т~Т~^~^~Т^В~Т^В~0~Ю ~11~^1Т1ТУв~Тб

к = 1 О

2

в

а

10 12 1*

1Б

к = 3 к-

Рис. 16. Порядок исключения неизвестных в методе редукции

аналогичным образом. При этом необходимо предположить, что чис­ ло узлов N является степенью двойки, N = 2m. Прежде чем перехо­ дить к случаю произвольного т , рассмотрим для наглядности слу­ чай тп=4, т. е. Л£= 16.

Обозначим черезk номер этапа исключения неизвестных. При

k=0 система уравнений совпадает с исходной и содержит значения неизвестных во всех внутренних узлах. На рис. 16 это соответст­ вует верхней горизонтальной черте, где кружочками отмечены но­ мера неизвестных yh входящих в систему. На следующем этапе

419

(6=1) происходит исключение неизвестных с нечетными номерами, в результате чего получаем систему вида (3), содержащую значе­

ния

неизвестных только в четных узлах. Этап k = \ изображен на

рис.

16 второй сверху горизонтальной чертой. Стрелки указывают,

какие неизвестные были исключены. На втором этапе (6 = 2) оста­ ется каждый четвертый узел и на заключительном этапе (й = 3) ос­ тается только одно уравнение, связывающее у8, у0, у1в. Поскольку Уо и у16 заданы (см. (2)), из последнего уравнения можно найти уа. Тем самым начинает осуществляться обратный ход в методе ис­ ключения. Зная уа, можно найти у4 и у12, далее — все неизвестные с четными номерами и, наконец, все остальные неизвестные.

Вернемся к общему случаю, когда N=2m. Согласно (3), в ре­

зультате первого этапа

исключения (k = \ )

получаем систему урав­

нений

 

 

 

 

 

 

 

где

+

 

£=

2, 4, 8,... , 2т - 2 ,

(4)

С(1) = (С<0') 2—2Е,

С<0) = С,

 

(5)

 

 

 

Fl>) = F l-1 + CFt + Fc+1.

 

(6)

По индукции

легко

доказать,

что

на

k-м

этапе исключения,

k = 1, 2, . . ., m, получаем систему

 

 

 

 

 

 

1 -

С * ~ \ + у

. ^

= -

F

t 1',

(7)

 

i = 2*-‘, 3 ■2ft- ‘, . . . , 2m2h~',

 

 

 

 

Уо—у 1>

У»= Р-2>

 

 

 

где матрицы С{к~1> и векторы F

находятся из рекуррентных со­

отношений

 

 

 

 

 

 

 

С(*’ =

(С<4_1))2 — 2Е,

k =

1, 2,. . „ я - 1 ,

 

с (о> =

Cj

 

 

 

 

 

(8)

F\k) =

Ff ^U + c lh-1)F t 1'+

F

t t i-

 

(9)

F f = F c ,

i = 2\ 2-2\ 3-2\ . . . , 2m—2\

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

С*4 и векторов Ff> по формулам (8), (9). Обратный ход состоит в нахождении векторов у; из системы (7), начиная с k — m.

Метод редукции в том виде, как он здесь изложен, не применя­ ется в реальных вычислениях по двум причинам. Во-первых, он не­ экономичен из-за того, что на каждом этапе приходится обращать матрицу Cik) общей структуры. Во-вторых, вычисление правых час­ тей по формулам (9) неустойчиво. В следующих пунктах будет по­ казано, как можно устранить указанные недостатки метода ре­ дукции.

420

Соседние файлы в предмете Численные методы