Лекция 14. Обращение преобразования Лапласа (формула Меллина). Восстановление оригиналов по известным изображениям. Решение дифференциальных уравнений операционным методом

Теорема 13.5 (Меллина).} Пусть F ( p) ÎC ¥ (Re p > a) и

Пример 13.2.

78

интегрирования замыкается вправо и обходится по часовой стрелке --- в отрицательном направлении)

Указание: можно считать контур интегрирования замкнутым налево и суммировать вычеты в ±iw.

Теоремы свертки и запаздывания.

Теорема. (теорема запаздывания) Если f(t) = 0 при t < 0, то справедлива формула

L[ f (t - t0 )] = e - pt0 L[ f (t)]

где t0 – некоторая точка.

t

Определение. Выражение ò f1 (t) f2 (t - t)dt называется сверткой функций f1(t) и f2(t) и

0

обозначается f1* f2.

Теорема. (теорема свертки) Преобразование Лапласа от свертки равно произведению преобразований Лапласа от функций f1(t) и f2(t) .

· t

F1 ( p)F2 ( p) =· ò f1 (t) f 2 (t - t)dt

0

Теорема. (Интеграл Дюамеля (Дюамель (1797 – 1872) – французский математик)). Если

··

F ( p) = f (t); G( p) = g(t) , то верно равенство

··

Пример. Найти изображение функции sin t . t

79

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

Пусть дано линейное дифференциальное уравнение с постоянными коэффициентами.

обращает исходное уравнение в тождество, значит функция, стоящая в левой части уравнения и функция f(t) имеет (по теореме единственности) одно и то же изображение Лапласа.

80

Этот многочлен зависит от начальных условий. Если эти условия нулевые, то многочлен равен нулю, и формула принимает вид:

x( p) = F ( p) Rn ( p)

Рассмотрим применение этого метода на примерах.

81

·

y = y = e2 x ;

·

Пример. Решить уравнение:

элементарные дроби. Воспользуемся делением многочленов(знаменатель делится без остатка на p – 1):

6p - 6

6p - 6

0

В свою очередь p 2 - 5 p + 6 = ( p - 2)( p - 3)

Получаем: p3 - 6 p 2 +11p - 6 = ( p -1)( p - 2)( p - 3).

Определим коэффициенты А, В и С.

A( p - 2)( p - 3) + B( p -1)( p - 3) + C( p -1)( p - 2) = -6 + p

Ap 2 - 5Ap + 6 A + Bp2 - 4Bp + 3B + Cp 2 - 3Cp + 2C = -6 + p

82

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

Пример. Решить систему уравнений:

Решим полученную систему алгебраических уравнений.

83

Если применить к полученным результатам формулы

Составим систему вспомогательных уравнений:

84

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

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

85

Лекция 15. Системы дифференциальных уравнений. Задача Коши. Теорема существования и единственности решения задачи Коши.

Структура общих решений однородных и неоднородных

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

y'1 = f1(x, y1, y2 , ..., yn), y'2 = f2 (x, y1, y2 , ..., yn),

..............................

y'n = fn (x, y1, y2 , ..., yn),

где x — независимая переменная, а

y1(x), y2(x), ..., yn(x) — неизвестные функции, n — порядок системы.

Обозначив

запишем систему в векторной форме

`Y '=`F(x,`Y ).

Решением системы называется вектор-функция `Y , которая определена и нерерывно дифференцируема на интервале (a, b) и удовлетворяет системе, т.е. для всех x0О (a, b) справедливо

`Y '(x) =`F(x,`Y (x)).

Задачей Коши (задачей с начальными условиями) называется следующая задача: найти такое решение`Y (x) системы `Y '=`F(x,`Y ), что `Y (x0) =`Y 0,

где x0 — заданное число, а `Y 0 — заданный вектор.

Интегральной кривой системы называется кривая в (n+1) -мерном пространстве Rn+1x,y, заданная уравнением `Y =`Y (x), где `Y (x) - решение системы.Таким образом, решить задачу Коши — это значит найти интегральную кривую, проходящую через заданную точку пространства Rn+1x,y.

·Система ОДУ первого порядка, разрешённая относительно производных

(нормальная система -го порядка)

·ОДУ -го порядка, разрешённое относительно старшей производной

86

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

Теорема 15.1(теорема существования и единственности решения задачи Коши.)

Если вектор-функция `F(x,`Y (x)) и ее частные производные по переменным yi , i = 1, 2, ..., n, непрерывны в области G пространства Rn+1x,y, то на некотором интервале (x0 -h, x0+h) существует единственное решение системы

`Y '(x) =`F(x,`Y (x)) ,

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

`Y (x0) =`Y 0,

т.е. через каждую точку области G проходит единственная интегральная кривая системы.

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

Опишем алгоритм решения задачи Коши для уравнения второго порядка y'' + a1 y' + a2 y = f(x), y(x0)=y0, (y)'(x0)=y0,1.

Будем искать решение задачи в виде y(x)= c1(x) y1(x) + c2(x) y2(x),

где y1(x), y2(x) — линейно независимые решения однородного уравнения y'' + a1 y' + a2 y = 0.

Вычислим y'(x), y''(x) и подставим полученные выражения в уравнение. Вычислим первую производную

y'(x)= (c1'(x) y1(x) + c2(x)' y2(x)) + (c1(x) y1'(x) + c2(x) y2'(x)), положим c1'(x) y1(x) + c2(x)' y2(x) = 0

и тогда y'(x)= c1(x) y1'(x) + c2(x) y2'(x), y''(x)= (y'(x))'= (c1(x) y1'(x) + c2(x) y2'(x))'=

=c1'(x) y1'(x) + c2(x)' y2'(x) + c1(x) y1''(x) + c2(x) y2''(x).

Подставив y(x) и ее производные в уравнение, получим: y'' + a1 y' + a2 y =

=c1'(x) y1'(x) + c2(x)' y2'(x) + c1(x) y1''(x) + c2(x) y2''(x) +

+a1(c1(x) y1'(x)+c2(x) y2'(x)) + a2(c1(x) y1(x)+c2(x) y2(x)) =

=c1(x)( y1''(x)+a1 y1'(x)+a2 y1(x)) + c2(x)( y2''(x)+a1 y2'(x)+a2 y2(x)) +

+c1'(x) y1'(x) + c2(x)' y2'(x) = 0 + 0 + c1'(x) y1'(x) + c2(x)' y2'(x) = f(x), при условии c1'(x) y1(x) + c2(x)' y2(x) = 0.

Тогда неизвестные функции c1(x) и c2(x) являются решениями системы линейных дифференциальных уравнений c1'(x) y1'(x) + c2(x)' y2'(x) = f(x),

c1'(x) y1(x) + c2(x)' y2(x) = 0

с известными y1(x) и y2(x).

Эта система легко разрешима относительно c1(x) и c2(x):

87

c1'(x) = f(x)y2(x)/(y1'(x)y2(x)-y1(x)y2'(x)), c1'(x)= f(x)y1(x)/(y1(x)y2'(x)-y1'(x)y2(x)).

Вычислив интегралы в правой части системы, получим

Произвольные константы C1 и C2 определяются из начальных условий.

ПРИМЕР 1. Решение методом вариации задачи Коши для линейного неоднородного дифференциального уравнения 2-го порядка с постоянными коэффициентами.

Заметим, что разрешимость системы дифференциальных уравнений для

c1'(x) и c2'(x) и однозначная разрешимость системы начальных условий для произвольных констант C1 и C2 гарантированы линейной независимостью y1(x) и y2(x), (y1'(x)y2(x)-y1(x)y2'(x))№0 для линейно независимых y1(x) и y2(x).

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

Решение задачи Коши ищем в виде

y(x)= c1(x) y1(x) + c2(x) y2(x) + ... + cn(x) yn(x),

где y1(x), y2(x), ..., yn(x) — линейно независимые решения соответствующего однородного уравнения.

Неизвестные функции c1(x) , c2(x), ..., cn(x)

находим как решения линейной системы дифференциальных уравнений c1'(x) y1(x) + c2(x)' y2(x) + ... + cn'(x) yn(x) = 0

c1'(x) y1'(x) + c2'(x) y2'(x) + ... + cn'(x) yn'(x) = 0, c1'(x) y1''(x) + c2'(x) y2''(x) + ... + cn'(x) yn''(x) = 0,

.................

c1'(x) y1(n-1)(x) + c2'(x) y2(n-1)(x) + ... + cn'(x) yn(n-1)(x) = f(x),

которая в силу линейной независимости y1(x), y2(x), ..., yn(x) разрешима относительно ci'(x). Вычислив ci(x) = Fi(x) + Ci находим произвольные постоянные Ci из начальных условий и тогда искомое решение уравнения имеет вид

y(x)= F1(x) y1(x) + F2(x) y2(x) + ...+ Fn(x) yn(x) + C1y1(x) + C2 y2(x) +...+ Cnyn(x).

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

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

записать характеристическое уравнение;

найти все корни характеристического уравнения l1, l2, ... , ln; найти фундаментальную систему решений y1(x), y2(x), ..., yn(x));

представить искомое решение задачи Коши в виде линейной комбинации

88

y(x)= c1(x)y1(x) + c2(x)y2(x) + ... + cn(x)yn(x),

с неизвестными функциями c1(x), c2(x), ..., cn(x); составить и решить систему для c1 (x), c2(x), ..., cn(x);

подставить вычисленные ci(x) = Fi(x) + Ci в выражение для решения и записать для него начальные условия;

найти из начальных условий значения констант Ci и записать искомое решение.

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

Для отыскания общего решения линейного неоднородного дифференциального уравнения с постоянными коэффициентами следует

найти общее решение соответствующего однородного уравнения (записать характеристическое уравнение, найти все корни характеристического уравнения l1, l2, ... , ln, записать фундаментальную систему решений y1(x), y2(x), ..., yn(x));

найти методом вариации произвольных постоянных любое частное решение неоднородного уравнения yч(x);

записать выражение для общего решения y(x)= c1 y1(x) + c2 y2(x) + ... + cn yn(x) + yч(x).

Доказательство формулы Лиувилля-Остроградского для линейной системы обыкновенных дифференциальных уравнений

Пусть вектор-функции — решения линейной системы ОДУ. Введем матрицу следующим образом

Тогда . Воспользуемся тем, что - решения системы ОДУ, то есть .