Решение задачи Коши для линейного неоднородного дифференциального уравнения с частными производными второго порядка. Принцип Дюамеля

Размещено на http://www.allbest.ru/

Тема: Решение задачи Коши для линейного неоднородного дифференциального уравнения с частными производными второго порядка. Принцип Дюамеля



Сначала мы хотим показать, что принцип Дюамеля для дифференциальных уравнений с частными производными является обобщением метода вариации произвольных постоянных для линейных обыкновенных дифференциальных уравнений.

Рассмотрим неоднородное уравнение с постоянными коэффициентами

(1)

где q(t) -- функция, не равная нулю тождественно, и однородное уравнение Ly=0, (2) соответствующее данному неоднородному.

Обозначим 1) через y(t) общее решение уравнения Ly=q(t); 2) через некоторое частное решение этого уравнения; 3) через z(t) - общее решение уравнения (2).

Чтобы найти z(t), достаточно знать корни характеристического уравнения, составленного для дифференциального уравнения Ly=0. А y(t) можно найти, используя метод вариации произвольных постоянных (метод Лагранжа).

Ограничимся случаем, когда уравнение (1) имеет вид

(3)

где k- некоторая положительная постоянная.

Уравнение (3) рассмотрим при t>0, а при t=0 зададим однородные начальные условия. Иными словами, поставим задачу Коши



(I)

Наряду с задачей (I) рассмотрим задачу Коши для однородного уравнения с неоднородными начальными условиями (которые в дальнейшем мы будем называть специальными)

(II)

Сначала решим задачу (II). Характеристическое уравнение

имеет корни Поэтому общее решение уравнения из (II)

(4)

где С₁ и С₂ -- произвольные постоянные.

Используем начальные условия из (II) и получим:



Поэтому решение v(t) задачи Коши (II) имеет вид

(5)

Метод вариации произвольных постоянных позволяет найти общее решение y(t) уравнения из (I) в виде (4), предположив что С₁(t) и С₂(t) -- пока не известные нам функции:

(6)

Для того, чтобы эти функции найти, надо решить систему из двух линейных алгебраических уравнений, неизвестными в которых являются производные функций С_i:

Решением такой системы оказываются функции



Проинтегрируем каждое из полученных равенств и получим

(7)

Под каждым из интегралов в формулах (7) мы понимаем сумму соответствующей первообразной и произвольной постоянной.

Подставим выражения для C₁(t) и C₂(t) из (7) в (6):

(8)

Функция y(t) в формуле (8) представляет собой общее решение уравнения из (I).

Найдем решение задачи Коши (I). Для этого в качестве первообразных для функций C₁(t), C₂(t) рассмотрим интегралы с переменным верхним пределом

соответственно. Тогда частное решение уравнения из (I)

(9)



и его производная

(10)

удовлетворяют однородным начальным условиям, поскольку все интегралы в формулах (9) и (10) обращаются в нуль при t=0.

Следовательно, формула (9) дает решение задачи Коши (I).

Запишем по-другому:

здесь в подынтегральном выражении стоит множитель

где v(t) -- функция из формулы (5), то есть решение задачи Коши (I) получается с помощью решения задачи Коши (II) по формуле

(11)

Формула (11) называется формулой Дюамеля.

Таким образом, можно, получив решение задачи Коши (II) со специальными начальными условиями для однородного уравнения, соответствующего данному неоднородному уравнению, вместо метода вариации произвольных постоянных применить формулу Дюамеля и найти решение задачи Коши (I) для неоднородного уравнения с однородными начальными условиями.

Разъясним, какой физический смысл можно этому придать

Метод импульсов. Интеграл Дюамеля

Рассмотрим бесконечное множество задач с неизвестными функциями

(12)

где функция q(t) отлична от нуля в интервале (0, ), >0, и равна нулю вне этого интервала. Предположим, что

(вскоре мы придадим этому условию определенный физический смысл). По формуле (11)

Пусть переменная t играет роль времени и пусть q(t) -- сила, действующая на некоторую систему, а y(t) -- смещение системы под действием силы q(t).

Пусть при t<0 система находилась в состоянии покоя (значит, y(0)=0) и y'(0)=0.

Пусть смещение системы происходит под действием силы q(t), которую мы ввели. В физике интеграл

называется импульсом силы q(t) на промежутке времени (0, t). Будем говорить, что решение задачи Коши (12) есть функция влияния этого импульса. По нашему предположению J(t)=1.

Из интегрального исчисления известно, что если функция (x) непрерывна в интервале [a, b], а функция (x) интегрируема в этом интервале, то интеграл

может быть представлен в виде

где a<c₀<b (обобщенная теорема о среднем).

Рассмотрим как параметр.

Повторим формулу для y(t) и применим к интегралу из этой формулы обобщенную теорему о среднем:



где с₀ -- число, 0< с₀<.

Тогда

Заметим, что когда 0, сила q(t) и импульс этой силы действуют на участке времени, стремящемся к нулю (в то же время по предположению для любого J(t)=1). Поэтому функцию

можно назвать функцией влияния мгновенного импульса интенсивности 1, действующего в момент t=0.

Замечание. Пусть y(t) -- координата точки, движущейся по оси y. Если в момент t=t₀ координата точки не меняется, а скорость получает конечное приращение, то такое воздействие на точку называется толчком

(или мгновенным ударом). Примером служит резкий удар кием по биллиардному шару.

Итак, решение v(t) задачи Коши (II) есть функция влияния мгновенного импульса силы интенсивности 1.

Теперь разъясним физический смысл формулы Дюамеля. Для этого разделим интервал времени (0, t) на n равных частей точками



Длину каждого из участков разбиения обозначим через

Рассмотрим n задач

Каждая из функций z_i представляет собой функцию влияния мгновенного импульса интенсивности Записав z_i как

используем начальные условия и найдем тем же способом, каким мы это уже делали в случае решения задачи (II), что

Вернемся к формуле Дюамеля (11). В нашем случае



Воспользуемся определением интеграла по отрезку и запишем

Напомним, что

q(t) -- непрерывно действующая на систему сила, а

импульс этой силы. Значит, y(t) -- функция влияния импульса непрерывно действующей силы.

Если n достаточно велико, а значит, промежутки времени

малы, можно с большой точностью считать, что

(13)

Эта последняя формула показывает, что функция влияния импульса непрерывно действующей силы есть сумма (суперпозиция, наложение) функций влияния мгновенных импульсов (принцип Дюамеля); как уже было сказано, каждое из слагаемых z_i представляет собой функцию влияния мгновенного импульса интенсивности в момент t=_i.

После перехода к пределу получаем ( формула (11)) точное значение y(t) в виде интеграла (континуальной суммы). Этот интеграл называется интегралом Дюамеля.

В формуле (13) суть метода импульсов.

Если рассматриваются задачи Коши I и II для уравнения Lyy”+p₁y'+p₀y=q(t), то формула Дюамеля имеет вид

где v(t) -- решение задачи

Итак, мы

1) решили задачу Коши I методом вариации произвольных постоянных;

2) полученную для y(t) формулу переписали в виде интеграла Дюамеля;

3) придавая функциям q(t), y(t), v(t) физический смысл, показали, что метод вариации можно рассматривать как метод импульсов.

Формула Дюамеля при этом связывает между собой решение задачи Коши I для неоднородного уравнения с однородными начальными условиями с задачей Коши II для однородного уравнения, соответствующего данному неоднородному, в которой сама функция (v(t)) в начальный момент времени равна нулю, а производная v'(0)=1. Здесь v'(0) -- скорость, которая сообщается системе при t=0.



Принцип Дюамеля для дифференциальных уравнений с частными производными

Вынужденные колебания бесконечной струны. Принцип Дюамеля

Рассмотрим дифференциальное уравнение

Поставим задачу Коши

(14)

Как и в случае обыкновенного дифференциального уравнения, рассматриваем функцию f(x, t) как внешнюю силу, действующую на некоторую динамическую систему. Дифференциальное уравнение в (14), в частности, описывает колебания струны [4], поэтому под f(x, t) можно понимать некоторую внешнюю силу, которая непрерывно (при t>0) действует на струну.

Кроме задачи (14) рассмотрим ряд задач Коши (при различных t=) для однородного уравнения:

(15)

где f(x, ) - значение функции f(x, t) из (14) при t=>0.

Чтобы подчеркнуть, что функция v зависит от , будем писать v=v(x, t, ). При каждом 0 решение задачи (15) находится по формуле Даламбера [4]. При этом, поскольку v=0 при t=, выражение для v(x, t, ) содержит только одно слагаемое:

(16)

Решение v(x, t, ) продолжим тождественным нулем при t. v(x, t, ) - функция влияния мгновенного импульса силы интенсивности f(x, ). Оказывается, что решение u(x, t) задачи (16) ищется по формуле Дюамеля:

(17)

(подынтегральная функция в нашем случае находится по формуле (16)). Интеграл в формуле (17) называется интегралом Дюамеля. Если эта формула верна, то, как и в случае с обыкновенным дифференциальным уравнением, она показывает, что u(x, t) - функция влияния непрерывно действующей силы- есть суперпозиция функций влияния мгновенных импульсов (принцип Дюамеля).

Проверим, что u(x, t) из (17) является решением задачи Коши (14). Предположим, что дифференцирование под знаком интеграла в (17) законно. Имеем



здесь v(x, t, t)=0 по первому из начальных условий задачи (15). Далее

здесь

по второму из начальных условий задачи (15).

Поэтому

поскольку v(x, t, ) - решение однородного дифференциального уравнения из формулы (15). Следовательно, u(x, t) -- решение неоднородного дифференциального уравнения из формулы (14). Функция u(x, t) удовлетворяет нулевым начальным условиям, поскольку



Формула (17) дает математическое выражение принципа Дюамеля. Повторим физическую интерпретацию этого принципа: результат воздействия на систему (на частицу) непрерывно действующей силы эквивалентен воздействию континуальной (то есть непрерывной) суммы последовательных толчков.

Обобщенный принцип суперпозиции

Ранее мы уже упоминали о принципе суперпозиции для линейного дифференциального уравнения.

В случае, когда i меняется от 1 до бесконечности, справедлив обобщенный принцип суперпозиции, если только вычисление производных от ряда фигурирующих в операторе Lu, можно совершить при помощи почленного дифференцирования этого ряда.

Обобщенный принцип суперпозиции справедлив и в том случае, когда линейное дифференциальное уравнение Lu=0 имеет бесконечно много решений вида U(x, t, ), где - некоторый параметр. Интеграл

где - произвольная ограниченная функция, является решением уравнения Lu=0, если производные, входящие в линейный дифференциальный оператор Lu, можно вычислить при помощи дифференцирования под знаком интеграла.

С математической точки зрения принцип Дюамеля можно рассматривать как следствие принципа суперпозиции, и поэтому он применим к задаче Коши только для линейного дифференциального уравнения (обыкновенного или с частными производными).

Отметим некоторые случаи, в которых возможно применение принципа Дюамеля. задача коши интеграл дюамель

1. Уравнение

с однородными начальными условиями, если x=(x₁,…,x_n) и

2. Уравнение

с однородными начальными условиями, если x=(x₁,…,x_n) и L - произвольный линейный дифференциальный оператор по x и t с переменными коэффициентами, который может содержать производную по t первого порядка, но не содержит производной по t более высокого порядка.

3. Задача Коши для уравнения теплопроводности в случае, когда xR¹:



4. Задача Коши

если xRⁿ, L(u) - произвольный линейный дифференциальный оператор по x c переменными коэффициентами.

5. Задача Коши для систем дифференциальных уравнений первого порядка:

[ Эту же задачу можно сформулировать в матричной форме. Для этого

Обозначим



Тогда система запишется в матричной форме:

Начальные условия на u_i(x, t) перейдут при этом в начальное условие для вектор-функции u(x, t): u(x, 0)=0].

Размещено на Allbest.ru

...