Операционное исчисление: как преобразование Лапласа решает дифференциальные уравнения

Задача знакомая до зубной боли: решить $y^{\prime \prime }+y=t$. Классический путь — общее решение однородного уравнения, угадывание частного (тут подходит $t$), сборка констант из начальных условий. Работает, пока правая часть — многочлен или экспонента. Подайте внутрь прямоугольный импульс или пилу — и подбор превращается в мучение. Операционное исчисление предлагает другой маршрут: перевести всю задачу в комплексную плоскость $p$, где дифференцирование становится умножением на $p$, а само уравнение — линейным алгебраическим, по духу крошечной СЛАУ. Решаем его там и возвращаемся назад по таблице. Пять строк выкладок вместо страницы подбора — за это метод любят электротехники и механики.

Оригинал и изображение: договор о переводе#

Оригинал — функция $f(t)$ при $t\ge 0$, кусочно-непрерывная (что такое непрерывность, вспоминаем отдельно) и растущая не быстрее экспоненты: $|f(t)|\le M{e}^{s_{0}t}$. Изображение — интеграл Лапласа от неё; $p=s+i\omega$ — комплексное число, с которым работаем по правилам комплексных чисел:

В этой полуплоскости $F(p)$ — аналитическая функция, то есть ведёт себя как настоящая «хорошая» функция из теории аналитических функций: дифференцируемая, без особенностей. Пара «оригинал ↔ изображение» взаимно однозначна: потеряв $F(p)$, оригинал восстановим, и наоборот. Линейность бесплатна: изображение суммы равно сумме изображений, константа выносится. Базовые картинки считаются прямо по определению: $1\to 1/p$, $e^{at}\to 1/(p-a)$ — и уже из них, плюс свойства ниже, собирается вся таблица. Полезная картинка для интуиции: интеграл Лапласа раскладывает сигнал по экспонентам $e^{pt}$, спрашивая у каждого $p$ — сколько в $f(t)$ такой экспоненты? Производная при дифференцировании превращается в умножение именно потому, что $\frac{d}{dt}{e}^{pt}=p{e}^{pt}$: экспоненты — собственные функции оператора дифференцирования. Отсюда и название «операционное исчисление»: дифференцирование есть оператор, и на экспонентах он действует как умножение на число.

Три свойства-двигателя: производная, смещение, свёртка#

Главный двигатель — дифференцирование оригинала. Интегрируем по частям ${\int }_0^{\infty }{f}^{\prime }(t)e^{-pt}dt$ и получаем связь с $F(p)$; применяем формулу дважды — добираемся до второй производной. Именно здесь метод покупает свою магию:

Второй двигатель — теорема смещения: умножение оригинала на экспоненту сдвигает аргумент изображения. Она объясняет, почему в таблице всё «с $a$» получается из базовых строк простой заменой $p\to p-a$:

Третий двигатель — свёртка. Произведению изображений отвечает свёртка оригиналов, а не произведение функций. Для рациональных изображений со сложными корнями знаменателя обратный переход часто быстрее через вычеты, но свёртка нужна там, где произведение уже собрано:

Решаем y'' + y = t: пять шагов операционного метода#

Задача Коши: $y^{\prime \prime }+y=t$, $y(0)=0$, $y^{\prime }(0)=0$. Обозначим $Y(p)=L\{y\}$ и переведём каждый член уравнения в плоскость $p$.

Структура ответа видна уже в $Y(p)$: знаменатель левой части $p^2+1$ переезжает в оригинал как свободные колебания $\sin t$, а полюс правой части $1/p^2$ — как прямой тренд $t$. Общее правило: однородная динамика живёт в полюсах левой части, вынужденный режим — в полюсах правой. Поэтому операционный метод не просто решает задачу — он показывает, из каких режимов собран ответ.

Ненулевые начальные условия меняют немного: слагаемые $py(0)+y^{\prime }(0)$ просто переезжают в правую часть алгебраического уравнения. Мини-пример на закрепление, без картинок: $y^{\prime }-y=e^{2t}$, $y(0)=0$. Изображения: $L\{y^{\prime }\}=pY$ (начальное значение нулевое), $L\{e^{2t}\}=1/(p-2)$, и уравнение мгновенно становится $(p-1)Y=1/(p-2)$. Значит $Y=1/((p-1)(p-2))$; разложение на простые дроби даёт $1/(p-2)-1/(p-1)$, а оригинал $y=e^{2t}-e^t$. Проверка: $y^{\prime }=2e^{2t}-e^t$, $y^{\prime }-y=(2e^{2t}-e^t)-(e^{2t}-e^t)=e^{2t}$ ✓, $y(0)=0$ ✓. Вся задача — три строки; сравните с вариацией постоянных. Системы линейных уравнений с постоянными коэффициентами переводятся в системы линейных алгебраических уравнений для изображений и решаются теми же приёмами. Дроби с комплексными корнями знаменателя ($p^2+{\omega }^2$ при желании разбирается на $(p-i\omega )(p+i\omega )$) — вход в мир вычетов.

Полюса изображения кодируют поведение оригинала — этим удобно пользоваться для самоконтроля. Полюс $F(p)$ в точке $p=-a$ отвечает слагаемому $e^{-at}$: затухание со скоростью $a$. Пара комплексных полюсов $-b\pm i\omega$ отвечает колебаниям $e^{-bt}\cos \omega t$: вещественная часть задаёт скорость затухания, мнимая — частоту. Полюс в нуле — константа или степенной рост. Поэтому по виду $Y(p)$ инженер читает будущее решения, ещё не возвращаясь к оригиналу: устойчива ли система (полюса слева от мнимой оси), будет ли дребезг, каков масштаб времени. Это тот же язык, на котором говорят передаточные функции в теории управления.

Ошибки, которые режут баллы#

• Забыть начальные условия: писать $L\{y^{\prime \prime }\}=p^{2}Y$ при ненулевых $y(0)$, $y^{\prime }(0)$ — правильная формула $p^{2}Y-py(0)-y^{\prime }(0)$.
• Перепутать направление сдвига: $e^{at}f(t)$ даёт $F(p-a)$, а не $F(p+a)$.
• Перемножать оригиналы вместо изображений: произведению $F\cdot G$ отвечает свёртка оригиналов, а не произведение $f\cdot g$.
• Искать «частное решение» в операционном методе — его роль играет разложение изображения на простые дроби.
• Не указывать полуплоскость сходимости: $1/(p-a)$ — это $e^{at}$ при $\mathrm{Re}p>a$; без области ответ не полный.
• Складывать изображения, не глядя на области: сумма $F+G$ живёт в пересечении полуплоскостей, и это влияет на оригинал.
• Не проверять ответ подстановкой в исходное уравнение — самая дешёвая страховка в этом разделе.

Итог маршрута: оригинал → изображение → алгебра → простые дроби → оригинал. Держите в голове пять шагов — перевод, алгебра, дроби, возврат, проверка — и метод перестаёт быть магией, превращаясь в конвейер. Полная таблица изображений, список свойств и формулы обратного перехода собраны в шпаргалке по ТФКП; следующий уровень владения методом — изображения импульсных функций и разложения через вычеты. Одна и та же механика потом всплывёт в уравнениях с запаздыванием и в теории управления — там без операционного исчисления никуда.