Метод суперпозиции

Материал из Википедии — свободной энциклопедии
Перейти к навигации Перейти к поиску

Метод суперпозиции — метод решения краевой задачи для линейных обыкновенных дифференциальных уравнений путём преобразования в задачу Коши.

Описание метода[править | править код]

Основная идея метода суперпозиции заключается в преобразовании граничной задачи для линейных обыкновенных дифференциальных уравнений к двум или нескольким задачам Коши, которые можно решить одним из методов решения задач Коши, например методом Рунге-Кутта. Это преобразование осуществляется путём представления искомого решения в виде линейной суммы нескольких функций , включающей столько неизвестных констант , сколько недостаёт начальных условий для приведения к задаче Коши. Затем это представление подставляется в исходное дифференциальное уравнение и в результате получаем систему из дифференциальных уравнений. При подстановке представления в условия для границ даёт возможность вычислить начальные условия задачи Коши и неизвестные константы .

Пример[править | править код]

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

(1)

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

(2).

Для приведения краевой задачи к задаче Коши недостаёт одного условия, поэтому представим решение в виде

(3)

с одной неизвестной константой .

Подставив это разложение в (1) получаем:

В этом уравнении оба слагаемых должны быть равны нулю.

(4)

(5)

Первое граничное условие в (2) принимает вид:

,

отсюда вытекает:

(6 a, b)

Начальные условия для производной найдем путём дифференцирования (3) в точке 0:

(7)

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

(8 a, b)

Из (6) получаем:

(9)

Граничное условие во второй точке имеет вид:

Из этого уравнения получаем:

. (10)

Итак, мы получили все начальные данные для задачи Коши. Граничная задача (1), (2) решается следующим образом:

  1. Интегрируем уравнение (4) с начальными условиями (6 a), (8 a) от 0 до 1. Получаем .
  2. Интегрируем уравнение (5) с начальными условиями (6 b), (8 b) от 0 до 1. Получаем .
  3. По формуле (10) вычисляем константу , которая в силу (9) является недостающим начальным значением.
  4. По формуле (3) вычисляем решение исходной задачи.

Литература[править | править код]

  • Ц. На Вычислительные методы решения прикладных граничных задач. - М.: Мир, 1982.