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

Материал из testwiki
Перейти к навигации Перейти к поиску

Шаблон:Underlinked

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

Описание метода

Основная идея метода суперпозиции заключается в преобразовании граничной задачи для линейных обыкновенных дифференциальных уравнений к двум или нескольким задачам Коши, которые можно решить одним из методов решения задач Коши, например методом Рунге-Кутта. Это преобразование осуществляется путём представления искомого решения x(t) в виде линейной суммы x(t)=x1(t)+C1x2(t)+...+CNxN(t) нескольких функций x1(t),x2(t),...xN(t), включающей столько неизвестных констант C1,...,CN, сколько недостаёт начальных условий для приведения к задаче Коши. Затем это представление x(t) подставляется в исходное дифференциальное уравнение и в результате получаем систему из N+1 дифференциальных уравнений. При подстановке представления в условия для границ даёт возможность вычислить начальные условия задачи Коши и неизвестные константы C1,...,CN.

Пример

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

f(t,x,dxdt,d2xdt2)=r(t) (1)

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

x(0)=x0,x(1)=x1 (2).

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

x(t)=x1(t)+C1x2(t) (3)

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

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

[f(t,x1,dx1dt,d2x1dt2)r(t)]+C1[f(t,x2,dx2dt,d2x2dt2)]=0

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

f(t,x1,dx1dt,d2x1dt2)=r(t) (4)

f(t,x2,dx2dt,d2x2dt2)=0 (5)

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

x1(0)+C1x2(0)=x0,

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

x1(0)=x0,x2(0)=0 (6 a, b)

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

dx(0)dt=dx1(0)dt+C1dx2(0)dt (7)

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

dx1(0)dt=0,dx2(0)dt=1 (8 a, b)

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

dx(0)dt=C1 (9)

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

x1(1)+C1x2(1)=x1

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

C1=[x1x1(1)]x2(1). (10)

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

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

Литература

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