Преобразование Лапласа: решаем дифференциальное уравнение

Преобразование Лапласа решает дифференциальное уравнение не интегрированием, а алгеброй: оно переводит производные в умножение на $s$, уравнение становится обычным алгебраическим относительно образа, а начальные условия подставляются автоматически. Вы делите, складываете дроби и читаете ответ из таблицы обратных преобразований. Ниже разберём пошаговый алгоритм и сразу дадим калькулятор, который строит решение по вашим коэффициентам.

Зачем нужен метод преобразования Лапласа

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

Преобразование Лапласа сворачивает всё в одну процедуру. Дифференцирование по времени превращается в умножение на $s$, начальные условия входят в образ сразу, а вместо подбора констант остаётся алгебра дробей. Поэтому метод доминирует в электротехнике, теории управления и механике колебаний - там, где на входе ступеньки, прямоугольные импульсы и экспоненты.

Рассчитать для своей задачи

Что такое преобразование Лапласа

Преобразованием Лапласа функции $f(t)$, заданной при $t \ge 0$, называют функцию комплексной переменной

$$F(s) = \int_0^{\infty} f(t)\, e^{-st}\, dt.$$

Образ $F(s)$ существует там, где интеграл сходится, - обычно в полуплоскости $\operatorname{Re} s > \sigma_0$. На практике интеграл считать не приходится: для всех типовых функций образы давно сведены в таблицу. Достаточно знать несколько строк: $1 \mapsto 1/s$, $e^{at} \mapsto 1/(s-a)$, $\sin \omega t \mapsto \omega/(s^2+\omega^2)$, $\cos \omega t \mapsto s/(s^2+\omega^2)$.

Главное свойство, ради которого всё затевается, - образ производной. Если $f(t) \mapsto F(s)$, то

$$f'(t) \mapsto s F(s) - f(0), \qquad f''(t) \mapsto s^2 F(s) - s f(0) - f'(0).$$

Именно здесь начальные условия $f(0)$ и $f'(0)$ встраиваются в решение без отдельного шага. Подробный вывод этих формул и тонкости с разрывами разобраны в заметке про преобразование Лапласа производной.

Пошаговый алгоритм решения

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

1. Применить преобразование к обеим частям. Каждую производную заменяем по формуле образа производной, подставляя $y(0)$ и $y'(0)$. Уравнение в производных превращается в алгебраическое относительно $Y(s)$.
2. Выразить образ $Y(s)$. Собираем все слагаемые с $Y(s)$ слева, всё остальное - справа, делим. В знаменателе появляется характеристический многочлен исходного уравнения.
3. Разложить на простые дроби. Если $Y(s)$ - дробно-рациональная функция, раскладываем её на сумму элементарных дробей с известными обратными образами.
4. Выполнить обратное преобразование. Каждую простую дробь переводим обратно по таблице, складываем - получаем $y(t)$.
5. Проверить начальные условия. Подставляем $t=0$ в найденную $y(t)$ и её производную: значения обязаны совпасть с заданными.

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

Пример: затухающий осциллятор

Решим уравнение свободных колебаний с затуханием

$$y'' + 2\delta y' + \omega_0^2 y = 0, \qquad y(0) = y_0,\ y'(0) = v_0.$$

Шаг 1 - преобразуем обе части. С учётом образов производных:

$$\bigl(s^2 Y - s y_0 - v_0\bigr) + 2\delta\bigl(s Y - y_0\bigr) + \omega_0^2 Y = 0.$$

Шаг 2 - собираем $Y(s)$:

$$Y(s) = \frac{(s + 2\delta)\, y_0 + v_0}{s^2 + 2\delta s + \omega_0^2}.$$

Знаменатель - характеристический многочлен, его корни $s_{1,2} = -\delta \pm \sqrt{\delta^2 - \omega_0^2}$. Дальше всё зависит от знака дискриминанта $\delta^2 - \omega_0^2$ - он и определяет режим решения.

Три режима решения и корни на плоскости s

Знак подкоренного выражения разбивает задачу на три случая, и каждому отвечает своё расположение корней на комплексной плоскости.

• Колебательное затухание ($\delta < \omega_0$): корни комплексно-сопряжённые $-\delta \pm i\omega_d$, где $\omega_d = \sqrt{\omega_0^2 - \delta^2}$. Решение - синусоида под экспоненциальной огибающей $e^{-\delta t}$.
• Критическое затухание ($\delta = \omega_0$): кратный действительный корень $-\delta$. Решение вида $(y_0 + (v_0 + \delta y_0)t)\,e^{-\delta t}$ - без колебаний, максимально быстрый возврат к нулю.
• Апериодический режим ($\delta > \omega_0$): два разных действительных отрицательных корня, решение - сумма двух затухающих экспонент.

Во всех трёх случаях корни лежат в левой полуплоскости ($\operatorname{Re} s < 0$) - это и есть условие устойчивости и затухания. Чем дальше корни влево от мнимой оси, тем быстрее гаснет решение.

Рассчитать для своей задачи

Калькулятор выше как раз показывает эту связь: двигаете $\delta$ и $\omega_0$ - и видите, как корни перемещаются по плоскости $s$, а форма решения $y(t)$ меняется от незатухающей синусоиды к плавному спаду.

Уравнения с правой частью и разрывами

Если в правой части стоит не ноль, а внешняя сила $f(t)$, алгоритм не меняется - просто справа добавляется образ $F(s)$:

$$Y(s) = \frac{(s + 2\delta) y_0 + v_0 + F(s)}{s^2 + 2\delta s + \omega_0^2}.$$

Главное преимущество метода раскрывается именно здесь. Ступенчатый вход $1(t)$ даёт $F(s) = 1/s$, единичный импульс $\delta(t)$ - просто $F(s)=1$, экспонента $e^{at}$ - $1/(s-a)$. Никаких сшивок решения на границах интервалов, как в классическом подходе: разрывная правая часть превращается в обычную дробь.

Когда $F(s)$ известно лишь как образ некоторой функции, обратное преобразование произведения $G(s)F(s)$ берут через теорему о свёртке Лапласа: произведению образов отвечает свёртка оригиналов. Это позволяет выписать решение даже для произвольного входа $f(t)$, не зная его явной формулы.

Обратное преобразование и простые дроби

Самый техничный шаг - вернуться от $Y(s)$ к $y(t)$. Для дробно-рациональных образов рецепт стандартный: разложить на простые дроби и перевести каждую по таблице.

Например, образ $Y(s) = \dfrac{s+3}{(s+1)(s+2)}$ раскладывается как

$$Y(s) = \frac{2}{s+1} - \frac{1}{s+2},$$

откуда сразу $y(t) = 2e^{-t} - e^{-2t}$. Для комплексно-сопряжённых корней дробь приводят к виду $\dfrac{s+\delta}{(s+\delta)^2 + \omega_d^2}$ и $\dfrac{\omega_d}{(s+\delta)^2 + \omega_d^2}$ - это образы $e^{-\delta t}\cos\omega_d t$ и $e^{-\delta t}\sin\omega_d t$. Полный список нужных переходов и приёмы их применения собраны в обзоре свойств преобразования Лапласа.

Частые ошибки

• Забыть начальные условия в образе производной. Писать $\mathcal{L}\{y'\} = sY$ вместо $sY - y(0)$ - тогда решение не удовлетворяет условиям задачи.
• Перепутать порядок в образе второй производной. Правильно $s^2 Y - s y(0) - y'(0)$, а не $s^2 Y - y(0) - s y'(0)$: множитель $s$ стоит при $y(0)$.
• Не проверить область сходимости для растущих функций - образ существует лишь правее всех полюсов.
• Ошибиться в разложении на дроби. Коэффициенты удобно находить методом подстановки корней знаменателя; всегда проверяйте обратной сборкой дробей.
• Применять метод к краевой задаче. Преобразование Лапласа рассчитано на начальные условия в нуле, а не на условия на концах отрезка.

FAQ

Чем преобразование Лапласа лучше классического метода? Оно встраивает начальные условия прямо в образ и превращает разрывную правую часть (ступеньки, импульсы) в обычную дробь. Не нужно отдельно искать частное решение и сшивать решение на границах интервалов - всё сводится к алгебре и таблице.

Можно ли так решать уравнения первого порядка? Да, и это даже проще. Для $y' + ay = f(t)$ образ сразу даёт $Y(s) = \dfrac{y(0) + F(s)}{s + a}$, после чего одно обратное преобразование возвращает $y(t)$.

Что делать, если корни знаменателя комплексные? Не раскладывать на дроби с комплексными числами, а выделить полный квадрат в знаменателе: $s^2 + 2\delta s + \omega_0^2 = (s+\delta)^2 + \omega_d^2$. Тогда дробь сразу читается как образ затухающего синуса или косинуса.

Коротко

Преобразование Лапласа решает дифференциальное уравнение, переводя его в алгебраическое относительно образа $Y(s)$: производные становятся множителями $s$, начальные условия входят автоматически, а знаменатель образа оказывается характеристическим многочленом. Алгоритм из пяти шагов - преобразовать, выразить $Y(s)$, разложить на простые дроби, обратно преобразовать, проверить условия - работает одинаково для однородных уравнений, уравнений с разрывной правой частью и систем. Корни знаменателя в левой полуплоскости задают режим затухания, а их мнимая часть - частоту колебаний.