<<
>>

6. Метод Галеркина-Петрова для нелинейных уравнений

Приближения и системы Галеркина. Пусть Е — сепарабель- ное нормированное пространство с базисом и F — отображение Е в Е*. Метод Галеркина приближенного решения уравнения F(u) = 0 заключается в следующем.

Сначала в пространстве Е задается базис, т. е. линейно независимая система векторов фьф2,фз, -ч обладающая тем свойством, что всякий вектор из Е представляется единственным образом в виде и = а*(к)ф*. При помощи этого базиса строится последовательность конечномерных подпространств (?„), где Еп — л-мерное подпространство, натянутое на векторы ф і, фг,..., ф„.

Приближением Галеркина решения уравнения F(u) = 0 называется вектор ип ? Е„, т.е. ^

(21)

*=i

удовлетворяющий системе уравнений

(*=1,2,..., л). (22)

*=1

Эта система, решение которой определяет приближение Галеркина к„, называется системой Галеркина. Отметим, что, несмотря на то, что система (22) совпадает с системой (13), вопрос о ее разрешимости решается здесь иначе. В четвертом разделе этой главы использовались свойства функционала /, а здесь при изучении системы (22) будут использованы свойства отображения F.

Связь с проекционными методами. Пусть Рп — оператор проектирования Е на Е„ и Р* — сопряженный оператор, который, как известно, проектирует Е* на л-мерное подпространство Е*. Проекционный метод приближенного решения уравнения F(u) = 0, где F — отображение Е в ?*, заключается в том, что данное уравнение заменяется уравнением в конечномерном пространстве

P;iF(Pnu)= 0, (23)

причем решение последнего уравнения называется приближенным ре-шением исходного уравнения. Покажем, что уравнение (23) эквивалентно системе (22). Действительно, если h — произвольный вектор из ?*, то уравнение (23) эквивалентно уравнению

О = {P*nF{Pnu), h) = (F(Pnu), Pnh) =

= (р (І «*ф*) , і Р.Ф») = і Р. {р (І я*ф*). ф»). 1=1 1=1 *=1

а это уравнение в силу произвольности (3, эквивалентно системе (22).

Таким образом, метод Галеркина решения рассматриваемого нами уравнения F(u) = 0 совпадает с проекционным методом решения этого уравнения. Отметим еще, что если Р — отображение Ех в Еу, где Ех и Еу — нормированные пространства, то проекционный метод решения уравнения Р(и) = 0 заключается в следующем. Задаются две последовательности подпространств (Ех'^) и (Еу"^) (где п — указатель размерности), объединения которых соответственно плотны в Ех и Еу, а также последовательности проекторов (Р„) и (Q„), где РпЕх — (Ех'^), QnEy = (Еу"^). Затем уравнение Р(и) — 0 заменяется уравнением

QaP(P„u) = О, (24)

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

Если Ех — Еу — Н, где Н — гильбертово пространство, то проекционный метод называется методом Бубнова-Галеркина; в общем случае проекционный метод носит название метода Галеркина-Петрова.

В заключение этого пункта отметим, что в силу (24), если F — отображение Е в Е, то приближения Галеркина для уравнения F(u) = 0 найдутся из уравнения в конечномерном подпространстве

P„F{P„u) = 0. (25)

Разрешимость систем Галеркина. Пусть F: Е -> Е* — монотонный непрерывный оператор, где Е — вещественное сепарабельное нормированное пространство, удовлетворяющий на сфере ||u|| = г > 0 условию

(F(u),u)>0 (||u|| = г). (26)

В этих условиях можно показать разрешимость системы (22). В силу эквивалентности систем (22) и (23) достаточно показать, что уравнение (23) имеет решение. Положим Р„и = w € Е„- Тогда в силу моно-тонности F и конечномерности Е,„ On(iv) = P*F(W) есть непрерывное отображение Е„ в Е* — Р*Е*, т.е. непрерывное отображение л-мерного пространства в п-мерное. Отождествляя Еп и Е*, находим, что Ф„ — непрерывное отображение л-мерного пространства в себя. Далее, так как при IMI = Г

(Ф„Н, W> = (p:F(w), w) = (F(w), W) > 0,

то на сфере ІИІ = г будем иметь |j<С>„(w)11 > 0. Отсюда находим, что уравнение Фп(и') = 0 имеет решение, принадлежащее шару ||w|| < г.

Лемма 7.

Пусть монотонный непрерывный оператор F: Е —> Е*, где Е — вещественное сепарабельное нормированное пространство, удовлетворяет на сфере ||м|| = г > 0 условию (26).

Тогда система Галеркина (22) разрешима при любом п и приближения Галеркина и„ удовлетворяют неравенству ||м„|| < г.

Сходимость метода Галеркина-Петрова. Здесь будем предполагать, что и„ — приближения Галеркина решения уравнения F(u) = 0, где F — отображение нормированного пространства Е в Е'. В данном случае приближения ип находятся из системы (22) или из уравнения (23).

Лемма 8. Пусть приближения Галеркина и„, удовлетворяющие системам (22), существуют при любом п, причем ||м„|| < г. Тогда если F — ограниченный оператор, то последовательность (F(un)) сходится Е -слабо к нулю.

Пусть Е — сепарабельное нормированное пространство с базисом фьфг.фз,---; (Е„) — последовательность подпространств, рассмотренных в п. 6.1; (Р„) — последовательность проекторов (РпЕ = Еп) и (Р*) — последовательность сопряженных проекторов. Как было отмечено в п. 6.2, Р* проектирует Е* на п -мерное подпространство Е*.

Пусть h — произвольный вектор из Е*, hn = P*h и М") — h-h„. Го-ворят, что последовательность {?,*} предельно плотна в Е*, если для всякого h?E* М") -> 0 при /?->«>.

Лемма 9. Пусть приближения Галеркина ип, удовлетворяющие уравнениям (25), существуют при любом п и ||м„|| < г = const. Тогда если F: Е —> Е есть ограниченный оператор и последовательность (Е*) предельно плотна в Е*, то последовательность (F(un)) сходится слабо к нулю.

В дальнейшем используется следующее определение.

Определение 17. Отображение G: Е -> Е" называется равномерно монотонным, если (G(u) — G(v), и — v) > ||u - v||у(||м - v||), где ty(t) — возрастающая вещественная функция, обращающаяся в нуль в нуле.

Справедлива следующая теорема сходимости.

Теорема 23. Пусть Е — рефлексивное вещественное банахово пространство с базисом {ф*}, а непрерывный равномерно монотонный и ограниченный оператор F: Е -> Е* удовлетворяет неравенству (F(u), и) > 0, если ||u|| > г > 0.

Тогда приближения Галеркина и„ существуют при любом п и сходятся к единственному решению ио уравнения F(u) = 0.

В условиях теоремы 23 метод Галеркина-Петрова можно использовать для отыскания приближенного решения уравнения F(u) = 0. Назовем вектор и„ при п = N из (21) N-м приближением к точному решению ио¦ Согласно теореме 23 для любого є > 0 существует номер N такой, что

— "oil < Е- Таким образом, используя метод Галеркина-Петрова, мы можем найти приближенное решение идг с любой наперед заданной точностью.

<< | >>
Источник: Агошков, Валерий Иванович. Методы решения задач математической физики:. 2002

Еще по теме 6. Метод Галеркина-Петрова для нелинейных уравнений:

  1. 3.2. Метод Ньютона для системы нелинейных уравнений
  2. 3.3. Метод итерации для нелинейной системы уравнений
  3. Методы численного интегрирования нелинейных уравнений движения
  4. Методы построения уравнений движения геометрически нелинейных стержневых механических систем
  5. Глава 7Методы решения нелинейных уравнений
  6. 17) Метод Фурье решения начально-краевых задач для однородного волнового уравнения (уравнение теплопроводности) с однородными краевыми условиями.
  7. Прямое численное интегрирование нелинейных уравнений движения
  8. Тема 5 Решение систем нелинейных уравнений (СНУ).
  9. 3.1. Отделение корней нелинейного уравнения.
  10. Модуль прямого численного интегрирования уравнений движения геометрически нелинейных стержневых систем
  11. 1.7. Решение систем нелинейных уравнений и задач оптимизации
  12. Решение краевой задачи для линейного дифференциального уравнения второго порядка методом прогонки
  13. Петров Г.В.. Краткий справочник по философии Составитель Петров Г.В. Псков: ПГПИ,2004. - 52 с., 2004
  14. 1.4 Нелинейная диэлектрическая проницаемость и методы её исследования
  15. Методика исследования сегнетоэлектрических материалов методом нелинейной диэлектрической спектроскопии
  16. Метод нелинейной диэлектрической спектроскопии
  17. Учет нелинейной зависимости сил инерции от перемещений в методах прямого численного интегрирования