Учебное пособие: Методы решения краевых задач, в том числе "жестких" краевых задач

B21 ∙ u + B22 ∙ m = t.

Следовательно, искомый вектор m вычисляется по формуле:

m = B12 ∙ (s – B11∙ u).

А искомый вектор n вычисляется через вектор t:

t = B21 ∙ u + B22 ∙ m,

n = t + N ∙ Y*(1←0).

В случае «жестких» дифференциальных уравнений предлагается выполнять поочередное построчное ортонормирование.

Запишем приведенную выше формулу

∙ K(1←0) ∙ =

в виде:


∙ K(1←x2) ∙ K(x2←x1) ∙ K(x1←0) ∙ = .

Эту формулу можно записать в виде разделения левой части на произведение матрицы на вектор:

[ ∙ K(1←x2) ] ∙ { K(x2←x1) ∙ K(x1←0) ∙ } =

[ матрица ] ∙ { вектор } = вектор

Эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

[ ∙ K(1←x2) ] { K(x2←x1) ∙ K(x1←0) ∙ } =

Здесь следует сказать, что подвектор t подвергать преобразованию не нужно, так как невозможно, так как его первоначальное значение не известно. Но подвектор t нам оказывается и не нужен для решения задачи.

Далее запишем:

[[ ∙ K(1←x2) ] ∙ K(x2←x1)] { K(x1←0) ∙ } =

[ матрица ] { вектор } = вектор

Аналогично и эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

[[ ∙ K(1←x2) ] K(x2←x1)] { K(x1←0) ∙ } = .

И так далее.

В результате поочередного ортонормирования получим:

В = ,

= .

Следовательно, искомый вектор m вычисляется по формуле:

m = B12 ∙ (s – B11∙ u).

7 Формула для начала счета методом прогонки С.К.Годунова

К-во Просмотров: 276
Бесплатно скачать Учебное пособие: Методы решения краевых задач, в том числе "жестких" краевых задач