Преобразование Лапласа: ОДУ и системы

Преобразование Лапласа переводит линейное дифференциальное уравнение с постоянными коэффициентами в алгебраическое уравнение относительно образа. Это и есть суть операционного исчисления: вместо интегрирования и подбора частного решения мы делим, складываем дроби и читаем ответ из таблицы. Метод особенно силён, когда заданы начальные условия и в правой части стоят разрывные сигналы, импульсы или экспоненты - то есть в задачах электротехники, механики колебаний и теории управления.

Определение и область сходимости

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

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

Здесь $s = \sigma + i\omega$ - комплексный параметр. Интеграл сходится в полуплоскости $\operatorname{Re}\,s > \sigma_0$, где $\sigma_0$ - абсцисса сходимости (зависит от роста $f(t)$ на бесконечности). Для всех функций, с которыми работает прикладная математика - экспонент, синусов, многочленов, ступенек, дельт - образ существует, и можно не следить за областью сходимости вручную: таблица преобразований и теоремы об операторах уже её учитывают.

Основные свойства

Линейность очевидна из определения: $\mathcal{L}\{\alpha f + \beta g\} = \alpha F(s) + \beta G(s)$. Главное же, ради чего метод применяют к ОДУ, - теорема о производной:

$\mathcal{L}\{f'(t)\} = sF(s) - f(0).$

Производная по $t$ превращается в умножение образа на $s$ с поправкой на начальное значение. Для второй производной формула расширяется:

$\mathcal{L}\{f''(t)\} = s^2 F(s) - s\,f(0) - f'(0).$

В общем виде $\mathcal{L}\{f^{(n)}\} = s^n F(s) - s^{n-1} f(0) - \ldots - f^{(n-1)}(0)$. Именно это свойство и позволяет уравнение с производными свести к алгебраическому: начальные условия зашиваются прямо в уравнение для $X(s)$, отдельной процедуры подбора констант не нужно.

Полезны ещё три оператора: сдвиг по аргументу образа $\mathcal{L}\{e^{at} f(t)\} = F(s - a)$, сдвиг по $t$ с функцией Хевисайда $\mathcal{L}\{u(t-a) f(t-a)\} = e^{-as} F(s)$ и умножение на $t$: $\mathcal{L}\{t\,f(t)\} = -F'(s)$.

Таблица типовых преобразований

Этой короткой таблицы хватает на 90% учебных задач:

Дельта-функция Дирака $\delta(t)$ и функция Хевисайда $u(t)$ - это удобные модели мгновенного импульса и включения сигнала. С их помощью кусочно-заданные правые части записываются одной формулой.

Алгоритм решения линейного ОДУ

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

1. Применить $\mathcal{L}$ к обеим частям. Производные превратятся в $sX(s) - x(0)$, $s^2 X(s) - sx(0) - x'(0)$ и т. д.
2. Алгебраически выразить $X(s)$. Получится дробь вида $X(s) = N(s)/D(s)$, где знаменатель $D(s)$ - характеристический многочлен исходного уравнения.
3. Разложить $X(s)$ на простые дроби методом неопределённых коэффициентов.
4. По таблице найти обратное преобразование каждой дроби - это и есть $x(t)$.

Шаги легко автоматизируются. Чтобы сразу увидеть, как метод раскладывается на твоей задаче, задай уравнение и начальные условия:

Разложение на простые дроби

Когда $X(s)$ получена как рациональная функция, обратное преобразование сводится к декомпозиции знаменателя на сомножители и подбору числителей. Для простого вещественного корня $s = a$ слагаемое имеет вид $A/(s - a)$, и по таблице $A/(s - a) \to A e^{at}$. Для кратного корня кратности $k$ возникает сумма $A_1/(s - a) + A_2/(s - a)^2 + \ldots + A_k/(s - a)^k$, прообраз каждого - $A_j t^{j-1} e^{at}/(j-1)!$. Пара комплексно-сопряжённых корней $a \pm i\omega$ даёт слагаемое $(As + B)/((s - a)^2 + \omega^2)$ - оно прообразуется в комбинацию $e^{at}\sin \omega t$ и $e^{at}\cos \omega t$. Коэффициенты $A_j, B$ находятся либо методом неопределённых коэффициентов (приравнять числители после приведения к общему знаменателю), либо методом подстановки и предельных переходов.

Теорема свёртки

Если изображения $F(s)$ и $G(s)$ перемножить, прообразом будет свёртка функций:

$\mathcal{L}^{-1}\{F(s) G(s)\} = \int_0^t f(\tau) g(t - \tau)\,d\tau.$

Эта теорема нужна, когда в правой части стоит произвольная функция $f(t)$ без табличного образа. Тогда решение однородного уравнения с импульсной правой частью даёт функцию Грина $G(t)$, а ответ на любую правую часть выражается свёрткой $G \ast f$. В электротехнике это классический способ найти отклик линейной цепи на произвольный сигнал через её импульсную характеристику.

Системы ОДУ

Метод напрямую переносится на системы. К вектору $\mathbf{x}(t) = (x_1, \ldots, x_n)^T$ применяется покоординатное преобразование, и система $\mathbf{x}' = A\mathbf{x} + \mathbf{f}(t)$ переходит в

$s\mathbf{X}(s) - \mathbf{x}(0) = A\mathbf{X}(s) + \mathbf{F}(s),$

откуда $\mathbf{X}(s) = (sI - A)^{-1}(\mathbf{x}(0) + \mathbf{F}(s))$. Матрица $(sI - A)^{-1}$ - это резольвента $A$; её элементы - рациональные функции $s$, к которым опять применяются разложение на простые дроби и таблица. Прообраз резольвенты - это матричная экспонента $e^{At}$, и метод Лапласа даёт прямой способ её вычислить, минуя жорданову форму.

Типовые задачи

Гармонический осциллятор $y'' + \omega^2 y = 0$ с $y(0) = y_0$, $y'(0) = v_0$ образуется в $s^2 Y - s y_0 - v_0 + \omega^2 Y = 0$; отсюда $Y(s) = (s y_0 + v_0)/(s^2 + \omega^2)$, и по таблице $y(t) = y_0 \cos \omega t + (v_0/\omega)\sin \omega t$. Если добавить вынуждающую силу $\sin \Omega t$ и сопротивление, получится задача о резонансе и затухании - те же шаги, только больше арифметики при разложении.

RC-цепь $RC\,u' + u = E$ с разряженным конденсатором $u(0) = 0$ переходит в $(RCs + 1)U(s) = E/s$, откуда $U(s) = E/(s(RCs + 1))$. Разложение даёт $u(t) = E(1 - e^{-t/RC})$ - каноническая зарядная кривая, ради которой и считают.

RLC-контур $L y'' + R y' + y/C = f(t)$ при импульсной правой части $f(t) = \delta(t)$ даёт импульсную характеристику, при ступеньке - переходную; для синусоидального сигнала ответ - устойчивая гармоника со сдвигом фазы. Все три случая считаются одним и тем же алгоритмом, отличается только $F(s)$.

Сравнение с методом Фурье

Преобразование Фурье $\hat{f}(\omega) = \int_{-\infty}^{\infty} f(t) e^{-i\omega t} dt$ работает на всей оси и не учитывает начальные условия - это инструмент анализа стационарных сигналов и спектров (близкий объект - характеристическая функция случайной величины, которая по сути есть преобразование Фурье плотности). Лаплас определён на $[0, \infty)$, естественно «съедает» начальные данные и применяется к переходным процессам. Формально, если $f(t) = 0$ при $t < 0$, $\hat{f}(\omega) = F(i\omega)$ - Фурье есть сужение Лапласа на мнимую ось. На практике выбор такой: задача Коши на полуоси - Лаплас, задача о спектре или фильтре - Фурье.

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

• Забыть начальные условия при преобразовании производных. Формула $\mathcal{L}\{f'\} = sF(s)$ верна только при $f(0) = 0$; в общем случае пропадает важное слагаемое.
• Перепутать знак в сдвиге: $\mathcal{L}\{e^{at} f(t)\} = F(s - a)$, а не $F(s + a)$.
• Не довести разложение на простые дроби до конца - оставить квадратичный сомножитель в знаменателе. Прообраз ищется только из табличных строк, а $1/(s^2 + bs + c)$ в таблице нет - нужно выделить полный квадрат.
• Подставить начальные условия в момент $t > 0$ вместо $t = 0$. Все формулы привязаны к $t = 0$; для сдвига начала задачи делайте замену $\tau = t - t_0$.
• Проигнорировать функцию Хевисайда при разрывной правой части. Без $u(t - a)$ нельзя записать включение сигнала и потерять часть решения.

FAQ

Когда метод Лапласа не работает? Если уравнение нелинейно (есть $y^2$, $\sin y$ и т. п.) или коэффициенты переменны и не сводятся к удобной форме - образа не получить в замкнутом виде. Для линейных уравнений с переменными коэффициентами иногда помогает преобразование $\mathcal{L}\{t f\} = -F'(s)$ - оно переводит задачу в новое ОДУ для $F(s)$, но это работает не всегда.

Как считать обратное преобразование, если в таблице нет нужной строки? Сначала пытаются разложить $F(s)$ на простые дроби - почти всегда хватает. Если задача аналитическая, можно применить формулу Бромвича-Меллина $f(t) = (1/2\pi i)\int_{c-i\infty}^{c+i\infty} F(s) e^{st} ds$ с интегрированием по контуру и вычислением вычетов в полюсах $F(s)$.

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

Коротко

Преобразование Лапласа - рабочий инструмент решения линейных ОДУ и систем с постоянными коэффициентами и заданными начальными условиями. Алгоритм всегда один: применить $\mathcal{L}$ к уравнению, выразить $X(s)$, разложить на простые дроби, прочитать прообраз из таблицы. Сила метода - в умении напрямую обработать импульсы, ступеньки, экспоненты и резонансы, а слабость - в нелинейных задачах и переменных коэффициентах, где придётся возвращаться к классическим методам.