
book1989
.pdfгде У;— искомые векторы размерности М, Д —заданные |
векторы, |
|
|||||||||||
Аи 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 = —Fь 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 |
Fу 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- ‘, . . . , 2m—2h~', |
|
|
||||
|
|
Уо—у 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