книги / Течения газа около тупых тел. Метод расчета и анализ течений А. Н. Любимов, В. В. Русанов
.pdfДля уравнений газовой динамики, в декартовых координатах имеем
Ц/ = |
(х^ф! + |
ХгОД + |
Хз»Фз) + &1С V (х ^ ) 2 + (х2ф2)2 + (х3ф3)2, |
||
где ‘0,1,2,з = 0, |
#4 = -г |
1, й5 = |
— 1. |
Очевидно, |
что для /= 1, 2, 3 щ всегда может |
обратиться в нуль при ср =4=О, а для |
I = 4 и 5 — только если и2 + V2 + ю2 ^ с2 . |
||||
Геометрически возможность обращения в нуль |
соответствует возможности построе |
ния стационарной характеристической поверхности соответствующего семейства. Из структуры матриц А, В у С в системе (3.15) следует, что одно из собственных
значений (т$ непосредственно выражается через элементы а5б, 655, сб5 точно так жег как в скалярном случае 6г7 выражается через а, Ь, с. Поэтому при 6 = 0 условие- (6.54) не обеспечивает устойчивости и необходимо вводить в схему параметры о*2 и о^. Выбор соответствующего 60 производится либо по (6.52), либо по (6.59). Строго говоря, вместо со0, определенного выше, в этих формулах должна стоять величина
шах {х2| Ь5б| + *з 1*551} = та х 1 + ЛгУ + г+г^и; | + ХдГ"1 1V)|} = < Щ
и соответственно вместо 60 — величина 60.
Для остальных собственных значений нижнюю границу можно определить толь ко по аналогии, пользуясь одной из тех же формул, но уже с со 0, а не с со0. Величины
60и 60 можно брать, вообще говоря, различными, тем самым полагая сг2 и сг3 разными в разных уравнениях. Это, однако, ведет к усложнению схемы и более удобно значе ния а2 и 0*3 задавать одними и теми же во всех уравнениях.
При 60 =^= 0 шаг по времени т определяется по-прежнему из условия (6.61), где-
со0 должно удовлетворять условиям (6.52) или (6.59), |
которые имеют вид: |
||
где |
$0> <?“ о» |
|
(6.62) |
|
|
|
|
О — 22'т |
1 |
д = |
2 / |
^ |
(2/ — I)2,7-1 2а — 1 * |
|
|
для формулы (6.52) и (? = 0.25, д = 5 для формулы (6.59). |
о*2х2 = 0*3X3. Подставляя |
||
Будем задавать о*2 и сг3 так, чтобы ог2//г2 = &3/к3 и 60 = |
в (6.62) 60 = о*2//г2тя, со0 = А™/к2т:п, получим окончательную формулу для определе ния соо и тп на слое I = Vх:
(6.63>
(6.64>
Коэффициент дх < 1 введен для обеспечения «запаса устойчивости».
§ 7. Устойчивость метода прогонки
В этом параграфе приведено полное исследование устойчивости алгоритма решения системы (5.12) методом прогонки. Построение алгоритма и устойчивость прямой про гонки были исследованы в [15] для системы с неособенными матрицами ЪтЛу,. Опытрасчетов показал, что это условие в нестационарном многомерном течении, как пра вило, неоднократно нарушается в процессе решения, что приводит к невозможности продолжать счет или потере точности. Алгоритм, изложенный в § 5 и исследованный ниже, не требует неособенности матричных коэффициентов системы (5 .12) и основан лишь на предположениях, необходимых для корректности системы.
Для исследования устойчивости применяется ряд понятий матричного исчисле ния, в частности свойства пучка матриц. Для удобства читателя необходимые сведе ния по теории матриц помещены в Приложении.
1. Алгебраическая инвариантность переноса граничных условий. Рассмотрим краевую задачу для системы разностных уравнений, аналогичной (5.12), но с линей ными граничными условиями:
= 2о» |
(7.1а) |
Дл+7* Х’ к+1 + Ьк+1/яХ к = %+у,, |
(7.16) |
УмХм = км, |
(7.1 в) |
к = 0, 1 , . . . , Л Г - 1 , |
|
где а, Ь, |л0, чм — матрицы размеров г X г, г X г, гг X г, г2 X г соответственно, при чем г1 + г2 = г.
Будем предполагать, что система (7.1) имеет единственное решение, т. е. что ее определитель
[А0 О О О О
&1/, Й1Д |
0 |
0 |
0 |
|
|
О |
Ьуг |
аа/а . . . |
О |
0 |
^ 2) |
0 |
0 |
0 |
ЪМ-у. Дм-1/* |
|
|
0 |
0 |
0 |
0 |
\м |
|
отличен от нуля. Отсюда следует (и это мы существенно используем в дальнейшем), что ранг матрицы, составленной из любого количества столбцов матрицы (7.2), равен их числу.
Покажем, что при любом т = 0,1,..., М существует соотношение вида
— 2т» |
(7.3) |
где ранг р,т равен гх, являющееся следствием уравнения |
(7.1а) и уравнений (7.16) |
при к = 0 ,1 ,..., т — 1 и такое, что всякое иное соотношение вида (7.3) получается из (7.3) умножением на неособенную матрицу порядка гг.
Наше утверждение очевидно при т = 0. Пусть т > 0. Опишем прежде всего однозначно определенный способ нахождения |Хт и §т . Для этого рассмотрим первые г1 + тг уравнений системы
= 2о, |
|
|
Дк+7*^к+1 ”Ь Ък+угХь — Як+1/,, |
к — 0, 1 , ... , т — 1. |
(7.4) |
Прямоугольная матрица, составленная |
из коэффициентов при |
неизвестных |
Х0,..., Х т _г, имеет тг столбцов, и ее ранг равен тг, так как остальные уравнения системы (7.1) не содержат Х0,..., Выберем среди отличных от нуля миноров тг-то порядка этой матрицы тот, который имеет наибольшее абсолютное значение, а если таких несколько, то выберем из них тот, сумма номеров строк которого мини мальна. Выразим затем из соответствующих уравнений Х 01\Х ц..., Х т_1 через Х т и правые части системы и подставим их значения в оставшиеся уравнения. Мы по лучим вместо (7.4) эквивалентную систему
Х& 4- МкХт = Укч |
к = 0, 1 , . . . , т 1 , |
(7.5а) |
|
= 2т, |
(7.56) |
где щ — квадратные матрицы порядка г, а матрица Дт |
имеет ранг г19 причем Дт и |
1гт определены однозначно. |
|
Любое соотношение вида |
|
Н’то-^т = ёту |
(7-6) |
являющееся следствием (7.4), есть также следствие эквивалентной системы (7.5), Рассмотрим какую-нибудь из строк (7.6). Она имеет вид
ат-^т == Рт»
где а т — матрица-строка с г элементами.
С другой стороны, произвольная линейная комбинация уравнений системы (7.5) имеет вид:
771—1 |
771—1 |
771—1 |
|
2 |
^к^к “Ь ( 2 |
^к^к “Ь ТРт) ^-т = 2 ^кУк “Ь Тёту |
(^•^) |
7с=о |
К*= 0 |
к=о |
|
где ъ'к и у — матрицы-строки с г и ?\ элементами.
Сравнивая (7.6) и (7.7), получаем, что все ук = 0 и ат = ур,т , Рт = у#т , откуда следует, что
И'т = |
Г'тН'т» |
ё т = |
^ т ё т • |
(^ •^ ) |
Очевидно, что аналогичное доказательство может быть проведено для соотноше |
||||
ния вида чтХ т = кт, т = М, |
М — 1,..., |
О, |
являющегося |
следствием (7.16) и |
(7.1в). |
|
|
|
|
Геометрически соотношения (7.6) можно интерпретировать как задание в про странстве Н Т г измерений г — гг = г2-мерной гиперплоскости Пт , в которой дол жен лежать конец вектора Х ш. Очевидно, что эта гиперплоскость не зависит от выбора матрицы Гт и, следовательно, граничное условие (7.1а), заданное при пг = 0, ин дуцирует посредством системы (7.16) инвариантную последовательность гиперпло скостей Пт . Строки матриц (хт можно интерпретировать как векторы сопряженного пространства Н* , и тогда последовательность р,т определяет в Н* последовательность подпространств Ь т размерности г, натянутых на строки р,т . Подпространства Ь т также инвариантны относительно способа вычисления р,т , так как умножение р,т олева на неособенную матрицу Гт меняет лишь базис в Ьт , оставляя само Ьт неиз менным.
2. Метод прогонки. Методы решения системы (7.1) с помощью переноса гранич ного условия (7.1а) называются методами прогонки. Процесс решения методом про гонки состоит из двух этапов. Сначала выполняется прямая прогонка, при которой вычисляются матричные коэффициенты р,т и правые части соотношений (7.6) для т = 1, 2,..., М . Последнее из них (при т = М) совместно с (7.1в) образует систему г уравнений для компонент вектора Х м• Разрешимость этой системы вытекает из предположенной разрешимости всей системы (7.1), что легко показать, используя алгоритм исключения. После определения Хм выполняется второй этап решения — обратная прогонка, состоящая в последовательном вычислении Х м+ь Хм+г*--*» •Х’о» для чего используются соотношения (7.6) и уравнения (7.16).
Существует целый ряд различных методов прогонки как для общих систем, так и для матриц а, Ь, обладающих известной спецификой [165—170]. Из доказанного выше следует, что множество значений Х т , удовлетворяющих (7.6) для данного т , не зависит от способа вычисления прогоночных соотношений, иначе говоря, что все методы йрогонки алгебраически эквивалентны.
Опишем один общий способ вычисления \хт и ^т , частный случай которого при гг = 1 применен в § 5 г.Пусть при некотором т известно соотношение (7.6). Присоеди-1
1 Очевидно, что способ вычисления рт п ут , использованный придоказательстве инвариантно сти в п. 1, имеет только теоретическое значение.
ним к нему уравнения (7.1) при к = т : |
|
\^тХт = |
|
Ът+Ч^Хт“Ь Ят+чХт +1 = ЯтН-**. |
(7.9) |
Ранг матрицы коэффициентов левой части равен числу ее строк г + г1э в противном случае определитель матрицы (7.2) был бы равен нулю. Пусть ранг матрицы
имеющей г столбцов, равен р ^ г. Тогда, исключив Х т из (7.9), мы получим г + + Г1 — Р > Г1 независимых соотношений между Х т+1. По доказанному выше их число равно точно и, следовательно, р = гх.
Если матрица Ът+ч, неособенная, то результат исключения легко выписывается
явном виде: |
|
Рт^771+V*^т+г/>Xт+1 = Рт^т+1/*^т+1/*— ёгп» |
(7.10) |
Очевидно, что (7.10) есть искомое прогоночное соотношение для к = т + |
1 . |
Если матрица &ж+./а особенная, то априори нельзя указать, какие именно строки матрицы %т образуют отличный от нуля минор, с помощью которого следует произ водить исключение Х т . Для построения алгоритма в этом случае необходимо либо использовать специфику матриц а, Ъ, либо применять общий метод исключения с вы бором главного элемента, автоматически обеспечивающий выбор минора, не равного
нулю [166]. |
= 1 оказывается возможным несколько видоизменить |
формулу |
Только при |
||
(7.10) так, что она |
остается справедливой и при особенной Ьто+1/а. Пусть Бе1; Ьт+г/л= А. |
|
Предполагая сначала, что А =^= 0, положим 6^+»/. = А~1ргп+1/„ где |
Рт+*/* |
транспонированная матрица алгебраических дополнений элементов Ъп&>. Подстав
ляя в (7.10), получим после умножения на А |
|
Р'тпРт-!-*,2ат+1/«^т+1 = 1АтРт+1/*Л'т+,/г — А^ж. |
(7.11) |
Покажем, что если г1= 1, то формула (7.11) остается справедливой и при А = 0 . Для этого прежде всего заметим, что если А = 0, то: 1) ранг Ьт+у, равен г — 1 и
2)прогоночное соотношение при к = т + 1 определяется независимо от (7.6). Первое очевидно, так как иначе р < г, что невозможно. Поэтому ранг рт+>/2равен
единице и рж+1/3= Еж+у* Ут+у2, где ?7П,у3— правый, а ут+уа — левый собственный век тор матрицы &тп+1/2, соответствующие ее нулевому собственному значению. Умножая (7.16) при к = т на ут+ч2слева, получим
Тт+1 |
т+У:Х 7П+! = Тт+Уз^т+У*. |
(7.12) |
Вектор ут+Уавт+у, не нулевой в |
силу неособенности^ матрицы |
(7.2), и (7.12) есть |
следствие (7.1а) и (7.16) при к = |
1 , 2,..., т. Поэтому (7.12) есть прогоночное соотно |
шение, определенное независимо от (7.6). |
Ст+у2, Гт+у._, получим |
|
Подставляя в (7.11) А = 0 и рт+1/2= |
|
|
(И’ТП^т+Уг) Ттп+у.^т+у*-^т+У« |
(Ит^т+Уг) Т'т+У*^т+,,1» |
(7.13) |
т. е. действительно это соотношение лишь множителем (р-шСт+у.) отличается от (7.12). Этот множитель не равен нулю, так как иначе система |х7Л^7П+1/1 = 0, &т+у,Ст+»;» = 0 имела бы ненулевое решение, что невозможно, так как ранг р ее матрицы равен г.
Формула (7.11) и была применена для вычисления прогоночных соотношений в алгоритме решения системы (5.12). Использование специфики матриц а и Ьпозволило записать выражения для рт+1 и #тП через рт , ^т , а7,.+,/, и Ът+цл в явном виде, что значительно сократило время решения задачи.
Процесс обратной прогонки, как и прямой, не определяется однозначным обрат зом. В самом деле, пусть найдено Хш+1, и рассмотрим (7.9) как систему относительно Х т. Из способа получения цт и %т следует, что эта система совместна, хотя и содер жит больше уравнений, чем неизвестных (г компонент вектора Х т). Поэтому для вы числения Х тможно использовать любые г комбинаций уравнений (7.9), разрешимые относительно Х т. Пусть Т — прямоугольная матрица размера г X (гх + г), умноже ние которой слева на (7.9) выделяет такие г комбинаций:
) х т + |
г ( |
° ) 1 га+1 = |
г ( |
|
0т-|-У* ) |
I |
) |
I ^т+у* |
|
или |
|
|
|
(7.14) |
(ТгМ'т Ч" Т2^т-{-У2) Х т -(- Т^(1т-^.у2Х тх.1 — Т1^т -(- ТгЯлг-р/з, |
где Т1'в.Т2 — матрицы, состоящие из гг первых и г последних столбцов Т соответствен но. Разрешая (7.14) относительно Х т , получим
где |
Х т = Лт-±.1{лХ т+1 Ч" Ут+У21 |
(7.15) |
||
|
|
|
|
|
Лт+Чг “ |
(^хИ'т Ч~ Тфт-\-У*) |
1 Т |
(7.15а) |
|
Ут+'/ж= (?1 |
+ |
^ 2^т+У2)П (Т^т + |
^ 2Лт+У2). |
(7.156) |
Коэффициенты Лт+чя и |
Ут + |
выражаются только через |
р,т , #т , ат+у2, Ьт+у„ |
я ы+«/. и поэтому они могут быть определены заранее, в процессе прямой прогонки. Хотя С/т+у. и Ут+Чжзависят от выбора матрицы Г, значение Х т формально от
него не зависит, в силу алгебраической совместности системы (7.9).
3. Устойчивость прогонки. Несмотря на алгебраическую эквивалентность спо собов вычисления прогоночных коэффициентов в прямой и функций Х т в обратной прогонке, не все способы в одинаковой мере пригодны для практического применения. Как в любом численном алгоритме, в алгоритме прогонки существенную роль играет его устойчивость. Устойчивость прогонки складывается из устойчивости вычисления матриц \1т и $тв прямой и векторов Х т в обратной прогонке. В силу линейности урав нений устойчивость можно определить как ограниченность норм \*>т,8 т,Х т линейной комбинацией норм правых частей системы (7.1) с постоянными коэффициентами. При этом устойчивость решения Х т , т. е. корректность краевой задачи (7.1), имеет смысл рассматривать независимо от устойчивости способа его получения, в данном случае метода прогонки.
В соответствии с принятой терминологией будем говорить, что задача (7.1) кор ректна, если ее решение Х т единственно и устойчиво относительно правых частей,
иначе говоря, если выполнено неравенство |
|
ш ах||Х т ||< С 0Х, |
(7.16) |
ТП |
|
где |
|
К = т а х {||’#0II, )| ки ||, т а х ||я*+./г ||}. |
|
к |
|
Будем говорить также, что метод прогонки устойчив, если |
|
||р 'т||^С ’1» |
(7.17а) |
|1 * » |< с ,а:. |
(7.176) |
Если М фиксировано и матрицы ц 0, а^Чж Ьл+у„ Ум. постоянны, то очевидно, что для корректности задачи (7.1) необходимо и достаточно, чтобы определитель (7.2) не равнялся нулю. Любой метод прогонки в этом случае будет устойчив, если он факти чески реализуем, т. е. если в алгоритме вычисления |1т , §т и Х т не встречается де ления на нуль. Конкретное построение таких алгоритмов прогонки описано выше.
Если система (7.1) есть разностная.аппроксимация краевой задачи для системы дифференциальных уравнений, то матрицы а, Ъсуть функции некоторой координаты
5 и параметра к (шага сетки): |
ак+«/а == а(Ък+г/%, |
к){ |
Ьк+Ъ = Ь(Ъкт1/г, к); &+*/. = |
— (к + 1/2) к. Число М зависит от к |
и М->-оо при к |
0. |
В определение корректности |
задачи в этом случае вводится дополнительное требование равномерного по М (или по к) выполнения неравенств (7.16). Аналогичное требование равномерного выполне ния (7.17) вводится и в определение устойчивости прогонки. Иначе говоря, величин** Со, С1г С2 не должны зависеть от М.
Отличие от нуля определителя (7.2) при любом к остается необходимым условием корректности задачи (7.1), однако достаточным оно уже не является. Необходимые и достаточные условия корректности для произвольных матриц а, Ъ были найдены В. С. Рябеньким [171]. Из структуры системы (3.15) следует, что матрицы а, Ь, оп ределенные формулами (5.13), приводятся одновременно к диагональному виду с действительными диагональными элементами. Иначе говоря, пучок матриц оа + + со Ъпри любых | и к имеет простую структуру и все его собственные значения дей ствительны (см. Приложение). Для этого случая необходимые и достаточные условия корректности, приведенные в [171], формулируются следующим образом.
1)Матрицы а(^, к) и Ь(^, к) суть непрерывные функции своих аргументов, и соб
ственные значения ог: (ог пучка оа + |
(&Ь при любых 6 и к ^ к0удовлетворяют ус |
|||||
ловиям: |
|
|
|
|
|
|
|
|
I | > 9011 |
| при |
г = 1 , 2, |
|
(7.18) |
|
|
|ш Г|<?о|°Г| |
при |
1 = г г + 1 ..... г, |
|
|
где (/0< 1 и каждое |
собственное значение считается столько раз, какова |
его крат |
||||
ность. |
удовлетворяют при любом к ^ |
к0условиям |
|
|||
|
2) Матрицы |х о |
|
||||
|
ОеЬ{ Ц о^о^О . |
ОеЬ{гдх<?Г}=^0, |
(7.19) |
|||
<?о = |
(?+(0, к) — матрица размера г |
X г1? |
столбцы которой суть правые |
собствен |
||
ные |
векторы пучка |
аа(0, к) + со&(0, /г), |
отвечающие |
собственным |
значениям |
|
|
(/= 1,2,3,..., гх). Аналогично (?Г= 0~(1» Щ — матрица размера г X г2, столб |
цы которой суть правые собственные векторы пучка сга(1 , к) + ю 6(1 , А), отвечающие
<т7 : со/ (/ = |
+ 1 ,..., г). |
Очевидно, что (7.19) аналогично рассмотренному в § 4 требованию допустимости граничных условий для системы гиперболических дифференциальных уравнений.
Критерии устойчивости различных методов прогонки всегда содержат, как не обходимые, условия корректности краевой задачи, а в ряде случаев и некоторые до полнительные условия на коэффициенты системы (7.1), вызванные спецификой ме тода. Естественно считать методы прогонки, не требующие для устойчивости этих дополнительных условий, в некотором смысле оптимальными.
Для корректной эадачи всегда возможно построить метод прогонки, для которого выполнено (7.17а). В самом деле, для этого достаточно, последовательно вычисляя матрицы |хш, любым способом вякий раз нормировать их так, чтобы элементы их были ограничены. Однако устойчивость обратной прогонки обеспечивается не при всяком способе такой нормировки. Существо дела состоит в том, что базис в Ьт , образован ный строками нормированной матрицы р,т , не должен быть близок к вырожденному {равномерно по М). Аналитически это условие можно сформулировать так: отноше ние определителя Грама, построенного на строках \кт1 к произведению сумм квадра тов их компонент должно быть ограничено снизу положительной константой, не за
висящей от М: |
|
(М'щИ'т) |
— /п \ Л |
иц.дГИлтдГ-! л »л р |
^ |
Ю |
|
равномерно по М. Здесь обозначает транспонированную матрицу и
IIIV*!2=2 |
г1- |
1 - 1 |
|
Необходимость выполнения неравенства (7.20) очевидна, так как в противном слу |
|
чае левая часть его может оказаться настолько малой, |
что в пределах погрешности |
вычислений строки матрицы рт будут зависимы. Простейшим способом нормировки, удовлетворяющим условию ограниченности элементов рт , является деление каждой строки матрицы цт на ее норму как вектора. Однако этот способ обеспечивает выпол нение (7.20) только если гг = 1.
Покажем это на примере системы с постоянными матрицами а, Ь, причем пусть матрица Ь — неособенная. Тогда имеет место (7.10), нормируя которое с комощыо некоторой матрицы Гт , получим
М*т+1Хт+1 — #771+1»
где |
Н'ТЛ+Х^ Г'7П|Х77г&1Я, |
#тп+1 = |
Г7п[^ттг^ ^тп+#/* ^тпЦт• |
|
(7.21) |
|
|
||||
При указанном |
способе нормировки |
матрица |
Гт — диагональная, |
и ее |
элементы |
суть обратные значения норм строк |
матрицы |
р,т+1 = ^ Ъ ^ а . Очевидно, |
что тогда |
||
каждая строка |
будет при возрастании тп приближаться к одному и тому же век |
тору — собственному вектору матрицы Ъ~га, отвечающему ее максимальному по мо дулю собственному значению. Следовательно, еслигх ^> 1,то базис, образованный строками \хт1 будет становиться вырожденным при возрастании т.
Обеспечить практически выполнение (7.20) при гх > 1 можно, если потребовать, помимо ограниченности элементов [хт , ортогональности ее строк. Тогда левая часть (7.20) будет всегда равна единице, т. е. своему максимально возможному значению.
Метод ортогонализации базиса в основном пространстве Н г был |
предложен и приме |
|
нен в ряде задач С. К. Годуновым [167]. Упомянутый ранее |
метод |
исключения с |
выбором главного элемента [166] также обеспечивает выполнение |
условия (7.20). |
В отличие от метода ортогонализации, в нем не требуется вычислять квадратные кор ни и исключена возможность случайной потери точности. Однако оба метода доволь но трудоемки, и их применение оправдано лишь в тех случаях, когда либо затрудни
тельно учесть специфику матриц, |
либо |
быстрое получение численного результата |
важнее, чем затраты машинного времени. |
|
|
Выполнение неравенства (7.20) само по себе еще не обеспечивает ограниченности |
||
норм $т. В самом деле, обращаясь снова |
к формуле (7.21), видим, что если ||Гт ||^>1, |
|
то ||#т || может возрастать, хотя |
|| р,т || |
остается ограниченной. Это может про |
изойти, например, в случае, когда все собственные числа Ь"1 по модулю меньше еди ницы, т. е. если не выполнено условие корректности (7.18).
Доказательство ограниченности $т и устойчивости обратной прогонки для кор ректной задачи и при условии выполнения неравенства (7.18), (7.19) требует более де тального изучения структуры прогоночных соотношений. Чтобы выяснить существо дела, не загромождая изложение сложными выкладками, мы предположим, что пра вые и левые собственные векторы пучка сга(^, к) + со &(^, к) не зависят от ^ и к. Это значит (см. Приложение), что существуют постоянные неособенные матрицы Р и (), приводящие а и 6 к диагональному виду:
*< «. + . * ) < ? - { |
ГаЕг — со2 |
0 |
1 |
(7.22) |
„ |
МЙ1_ |
з 0 }. |
Матрицы Е г и Е2 — единичные, а 2 и й — диагональные порядков гх и г2 соответст венно.
Собственные значения пучка оа + юЬ выражаются через диагональные элементы
2 и & следующим образом: |
|
|
|
б,+ :со| = |
о,:1, |
|
2 = 1 , 2 , . . . , ^ , |
сГ:со7 = |
1 :ог, |
|сог|< д 0, |
/ = 1, . . г2. |
Величины ог и со г могут зависеть от |, Н1но при этом неравенства (7.23) не нарушают ся.
Разобьем () на две прямоугольные матрицы |
и , составленные из |
правых соб |
|
ственных векторов пучка, отвечающих собственным значениям с г |
1 и |
1 : со 1 соот |
|
ветственно. Тогда (? = {()+, (?“}, где ()+ имеет размер гх г1? а |
— размер г X г2. |
Точно так же разобьем матрицу Ф = Р -1, а матрицы Р и Т = |
^_1 разобьем аналогич |
||||
но по строкам: |
|
( Р + |
|
|
|
Ф = {ф+, ф-}; |
Р = |
Т = |
|
|
|
| р _ |
|
|
|||
Очевидно, что |
|
|
Р±ф+ = о, |
|
|
Р+ф+ = Еъ |
Р~Ф- = |
|
|
||
Т +<?+ = Ег, |
Ф“(Г = |
Е,, |
= |
о, |
(7.24) |
а = ф+Т+— Ф-ЙЧ'-, |
Ъ= —ф+2 Т - + Ф -Т-. |
Выберем в пространстве Нг векторов X новый базис, образованный правыми соб ственными векторами пучка, т. е. столбцами матрицы (?. Пусть ЗС — вектор, образо ванный компонентами X в новом базисе, тогда
Х = <2&= |
+ <2~ЗГ, |
(7.25) |
ЗС = (Г1Х = т х , |
ЗС± = Т±Х. |
|
Здесь ЗС+ и ЗС~ — векторы, составленные из гг первых и г2последних компонент ЗС. Любое соотношение вида (7.6) записывается в новом базисе так:
+ |
==**,, |
(7-26) |
где
«А* = И™@+» А п = Н'тф .
Применив это к соотношениям (7.1а) п (7.1в) системы (7.1), получим
•А+^о + А "А ' = л , |
А |
А |
+ |
А А = |
(7.27) |
|
По условию корректности (7.19) матрицы А = |
| а0(?+и |
неособенные и, |
||||
следовательно, (7.1а) и (7.1в) можно умножить слева на (М+У1 и (Жм)"1. |
||||||
Подставив Х т = (){!Ств (7.16) и |
умножив |
эти уравнения слева на Р, получим |
||||
окончательно запись (7.1) в новом базисе: |
|
|
(7.28а) |
|||
|
•2? о - И ^ = |
|
||||
|
— 2 к+1/а«й?^ + |
«27^+1= |
5&л+уц |
(7.286) |
||
|
|
|
= |
^к+,л» |
||
|
|
|
(7.28в) |
|||
|
—■ЪиЗОи + |
«2?м = |
|
|||
где |
А = - ( Х Г А |
|
2 лг = - № ) -1^ , |
|
||
|
|
|
||||
$о = ( |
ЗС^1* = -Р^Як+у,, |
№\м = (Л/>м) |
^лг. |
Так как матрицы Р хъ(2 неособенные, тосистемы (7.1) и.(7.28) эквивалентны. Со отношения (7.26), полученные из (7.6) преобразованием базиса,.являются прогоночными соотношениями системы (7.28) и инвариантны относительно способа их получе ния. Имея это в виду, покажем, что любое соотношение (7.26) может быть записано так:
|
|
|
|
= |
|
|
(7.29) |
Очевидно, (7.29) верно при т = |
0. Исключая |
3?т, ЗСт из (7.29) и (7.286), при |
|||||
к = т получим |
|
|
|
|
|
|
|
|
|
<2?т+1 — |
771+1 «27т+1 = |
$ т + Ь |
|
|
|
где |
|
И^п+1 = 2т+|/*^т^т+,;/ |
|
(7.30а) |
|||
|
|
|
|||||
|
Эт-г1 = |
|
|
+ 31т+'1*+ 2т+1/2И^7П5?т+1/2. |
(7.306) |
||
Таким образом, (7.29) определено рекуррентно при всех т. |
|
|
|||||
Применяя (7.30а) последовательно, получим |
|
|
|
||||
|
ТУ*П= |
|
|
|
• . . Йт _здйт _|Д. |
(7.31) |
|
Так как из (7.23) следует, |
что |
||2 т+1/а|| ^ д0 и ||йт+у3|К д „ , |
то ||И^т | К д»™||й^0||. |
||||
Следовательно, ||РРт || равномерно по М ограничены и |
0 при т->- оо. |
||||||
Далее, |
|
|
|
|
|
|
|
I! З о и у, + |
р т+7^ т ^ т |
+у, II < II 31^Н II + до II И^оII №&*ЧшII < |
с ъж, |
||||
где ЭС = т а х {||#0||, |
||3^м||, |
т а х |
||5ГЛ+*д||} и Съ не зависят от М. |
|
|||
Следовательно, |
|
к |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 2 ? т +1 | | < д о | | 2 и + с 5^ |
|
|
|||
и |
|
|
|
• + д Г 1) с5ж< $ ц% || + |
|
д0Г 1 сьж. (7.32) |
|
II 5?т II < дот и% | ■+ (1 + |
до + . . |
(1 + |
Итак, ||$ т || равномерно по М ограничены и прямая прогонка устойчива в новом ба
зисе.
В силу инвариантности прогоночных соотношений (7.6) отличается от (7.29)
только умножением на некоторую неособенную матрицу Гт , т. е. |
|
цт = Тт (У+ - И ^ П , ёт = Гт9т. |
(7.33) |
Из (7.33) следует, что если нормировка цт удовлетворяет условию (7.20), то выполне но и (7.176), т. е. прямая прогонка устойчива в исходном базисе. В самом деле, из
ограниченности ||р,т || следует |
ограниченность ||Гт ||, а следовательно, и ограничен |
|
ность е т. |
|
из ]Ут ->• 0 следует, что подпространство векторов X , |
Заметим, что при т |
оо |
удовлетворяющих условию \1т X = 0, стремится к подпространству, натянутому на собственные векторы ()~. В самом деле, предельное подпространство определяется
системой уравнений |
&+ = 'У+Х — 0, общее решение которой есть X = |
где |
||
— произвольный |
вектор с г2 составляющими. |
|
||
После окончания прямой прогонки получаем систему уравнений для Хм» |
|
|||
|
ЦмХм = 1 м {Т + — Игм хР~}Хм = ёл/ 1 |
|
||
|
ум Х = ^ м ) " 1{— |
|
+ 'Г'} Х м = Ллг. |
|
Матрицу этой системы можно записать так: |
|
|
||
/ |
I _ ] Гз' 0 |
1 ! |
Е1 ~ ^ м \х р |
|
Ь л , / _ 1 0№ ) - » /\ - г м |
Е г / |
* |
Еслд выполнено условие (7.20), то определитель Гм равномерно по М ограничен сни зу. и поскольку элементы Гм ограничены, то ||Гм11 равномерно по АТ ограничена свер-
{ Н'Л/ ) 1
Чм I
ограничена и, следовательно,
|Х м |< С ,Я . |
(7.34) |
Нам осталось показать устойчивость обратной прогонки, т. е. что при надлежа - щем выборе матрицы Т из (7.18), (7.19) следует (7.16).
Применяя последовательно (7.15), получим:
Х т= |
М-1-^М + Ут, М-2^ДГ-V* + • • • + пУт+У* + |
^т+*/а |
|
где |
^ 1п.т+к = ^ т+У»Ит+У* •••IIтп+к+у*« |
|
|
Если |
|
||
ЦУт| т+*1< Л в*+1. |
(7.35) |
||
|
|||
где и у <; 1 не зависят от т и /с, то |
|
||
II Хм К |
Я {Ям~т\\Хм\\ + (1 + д+-..+ Чь м )тах.| |
1}< |
|
|
< ЩС*К + (1 - д)-1т а х | У*+у, ||}. |
|
|
|
к |
|
Следовательно, обратная прогонка будет устойчива, если выполнено (7.35) и ||Уц+'и\\^С7К. Покажем, что это действительно имеет место при подходящем выборе мат рицы Т. Для этого выберем Т так, чтобы г линейных комбинаций состояли из: 1) гх прргоночных соотношений и 2) г2линейных комбинаций строк разностных уравнений,
т. е.
где ТХ1 и Т22 — матрицы размеров |
гх X гх и г2 X г. |
|
|
|
|
|
|
|||||||||
Используя (7.15) и (7.24), получим: |
|
|
|
|
|
|
|
|
|
|||||||
|
Ит+7* -- |
' |
( ТиРтп V-1 |
Г 0 |
| |
|
|
|
|
|
||||||
|
1 |
Т22^771+7* ] |
|
1 7'22а т + У , |
) |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
•\г |
|
|
|
|Г ц |1 т |
I - 1 |
| 2*115т |
1 |
|
|
|
|
|
||
|
|
* я»+7* — 1 Т 22 Ьтп+у, |
] |
1 Г 22ЗСт+'/г |
1 |
|
|
|
|
|
||||||
|
( ГиШ |
1 |
|
\ |
|
ТцТт |
|
- |
ГпГ т И Ц |
( |
|
|
|
|
||
|
н/. ) |
|
1 — Га2Ф+2 т+7, |
|
Т22Ф- |
/ |
1 |
|
|
|
|
|||||
|
|
|
|
|
[ ГиГщО |
М |
|
|
|
|
|
|
|
|
||
|
|
1 | Г и Г . 0 | | |
I— |
|
|
Т |
ф - I * |
|
|
|||||||
| |
о |
о |
|
о |
|
|
| |
|
(7.36) |
|
||||||
|
0 |
|
|
|
|
|||||||||||
1ГИат ../ ./ “ |
1 |
|
|
0 |
Я ./1 7 ..Ф * |
|
|
|
а! |
’ |
|
|
||||
|
|
|
|
А |
|
Е , |
|
|
- И |
Ц |
- |
Ч |
о |
ПИ-7*К |
0 |
|
|
С/т+7. - |
- |
^ |
|
|
\ _ 2 ’ИФ+2т+7, |
|
ГМФ -/ |
1ГМФ+ |
-Г„Ф - |
|
у _ * - * ( |
- " ' • П Р ь |
1. |
(7.37) |
* 1-Г„Ф*2т+71 |
Гпф-| \Т гхЖт^ы\ |
|