Операционный метод решения дифференциального уравнения

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

В чём суть операционного метода

Идея операционного исчисления восходит к Хевисайду и формализована через преобразование Лапласа. Каждой функции-оригиналу $f(t)$, определённой при $t \ge 0$ и не растущей быстрее экспоненты, ставится в соответствие функция-изображение комплексного аргумента:

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

Обозначают это соответствие так: $f(t) \risingdotseq F(p)$, или короче $F(p) = L\{f(t)\}$. Главное свойство, ради которого всё затевается, - операция дифференцирования в пространстве оригиналов превращается в умножение на $p$ в пространстве изображений. Это и есть «алгебраизация»: трудная операция (взятие производной) заменяется лёгкой (умножение на число-параметр).

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

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

Ключевое свойство: изображение производной

Центральная формула операционного метода - изображение производной. Если $f(t) \risingdotseq F(p)$, то для первой производной справедливо:

$$f'(t) \risingdotseq p F(p) - f(0).$$

Для второй производной формула наращивается:

$$f''(t) \risingdotseq p^2 F(p) - p f(0) - f'(0).$$

В общем случае для $n$-й производной:

$$f^{(n)}(t) \risingdotseq p^n F(p) - p^{n-1} f(0) - p^{n-2} f'(0) - \dots - f^{(n-1)}(0).$$

Обратите внимание: начальные условия $f(0)$, $f'(0)$ и так далее входят в формулу сразу. Именно поэтому операционный метод особенно удобен для задачи Коши - не нужно сначала искать общее решение с произвольными постоянными, а потом отдельно определять эти постоянные из начальных данных. Условия учитываются автоматически на этапе перехода к изображению.

Алгоритм решения по шагам

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

1. Переход к изображениям. Применяем преобразование Лапласа к обеим частям уравнения. Каждую производную заменяем по формуле изображения производной с учётом начальных условий, а правую часть - её изображением по таблице.
2. Решение алгебраического уравнения. После замены получаем алгебраическое уравнение относительно изображения $X(p)$. Решаем его как обычное уравнение: выражаем $X(p)$ в явном виде.
3. Подготовка к обратному переходу. Полученное $X(p)$ - обычно дробно-рациональная функция. Раскладываем её на простейшие дроби (метод неопределённых коэффициентов), чтобы каждое слагаемое нашлось в таблице соответствий.
4. Обратное преобразование. По таблице (или формуле обращения) переходим от $X(p)$ к оригиналу $x(t)$. Это и есть искомое решение задачи Коши.

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

Ключевой технический навык здесь - разложение на простейшие дроби. Именно на нём чаще всего застревают: без правильного разложения обратное преобразование «не находится» в таблице.

Таблица основных соответствий

Минимальный набор оригиналов и изображений, который покрывает большинство учебных задач:

$$\begin{aligned} 1 &\risingdotseq \frac{1}{p}, \\ t^n &\risingdotseq \frac{n!}{p^{n+1}}, \\ e^{at} &\risingdotseq \frac{1}{p-a}, \\ \sin \omega t &\risingdotseq \frac{\omega}{p^2 + \omega^2}, \\ \cos \omega t &\risingdotseq \frac{p}{p^2 + \omega^2}. \end{aligned}$$

Дополняют их две теоремы смещения. Теорема затухания: если $f(t) \risingdotseq F(p)$, то $e^{at} f(t) \risingdotseq F(p-a)$ - умножение оригинала на экспоненту сдвигает изображение по аргументу. Теорема запаздывания работает с функцией Хевисайда и нужна, когда правая часть «включается» не сразу.

Именно таблица превращает обратное преобразование из интеграла в узнавание: разложив $X(p)$ на простейшие дроби, мы просто сопоставляем каждое слагаемое строке таблицы.

Разбор типового примера

Решим задачу Коши $x'' + x = 0$ с условиями $x(0) = 0$, $x'(0) = 1$ операционным методом.

Переходим к изображениям. Пусть $x(t) \risingdotseq X(p)$. Тогда:

$$\bigl(p^2 X(p) - p \cdot 0 - 1\bigr) + X(p) = 0.$$

Это уже алгебраическое уравнение. Решаем его относительно $X(p)$:

$$p^2 X(p) - 1 + X(p) = 0 \quad\Rightarrow\quad X(p)\,(p^2 + 1) = 1 \quad\Rightarrow\quad X(p) = \frac{1}{p^2 + 1}.$$

По таблице $\dfrac{1}{p^2 + 1}$ - это изображение синуса с $\omega = 1$. Следовательно:

$$x(t) = \sin t.$$

Проверка элементарна: $\sin 0 = 0$, $(\sin t)' = \cos t$, $\cos 0 = 1$ - оба начальных условия выполнены. Вся «дифференциальная» часть работы свелась к одному делению и одному обращению по таблице.

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

Когда операционный метод особенно выгоден

Метод даёт наибольший выигрыш в нескольких ситуациях. Во-первых, при задаче Коши с заданными начальными условиями - они встраиваются автоматически, не приходится решать систему на произвольные постоянные. Во-вторых, при разрывной или импульсной правой части (включение источника, удар, ступенька): функции Хевисайда и Дирака естественно живут в операционном исчислении, тогда как классические методы требуют сшивки решений на интервалах. В-третьих, при системах линейных уравнений - изображения превращают систему в линейную алгебраическую относительно $X_1(p), X_2(p), \dots$, и её решают теми же приёмами.

Ограничения тоже стоит помнить. Операционный метод в стандартном виде работает для линейных уравнений с постоянными коэффициентами и для функций, удовлетворяющих условиям существования изображения (заданы при $t \ge 0$, экспоненциальный рост). Для нелинейных уравнений или переменных коэффициентов прямой переход к изображению не даёт алгебраического уравнения.

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

• Забыли начальные условия в формуле производной. Самая частая ошибка - написать $f'(t) \risingdotseq p F(p)$ без слагаемого $-f(0)$. Тогда задача Коши решается неверно: условия просто теряются.
• Перепутан порядок начальных данных у второй производной. В формуле $p^2 F(p) - p f(0) - f'(0)$ при $f(0)$ стоит множитель $p$, а $f'(0)$ идёт без множителя. Обратный порядок - типичная описка.
• Неверное разложение на простейшие дроби. Если знаменатель имеет кратные или комплексные корни, наивное разложение даёт неверные коэффициенты, и оригинал «не находится» в таблице.
• Применение к уравнению с переменными коэффициентами. Операционный метод в базовом виде на них не рассчитан: умножение на $t$ даёт дифференцирование изображения, и алгебраического уравнения не получится.
• Путаница оригинала и изображения при ответе. Финальный ответ - это оригинал $x(t)$, а не изображение $X(p)$. Остановиться на $X(p)$ значит не довести решение до конца.

FAQ

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

Обязательно ли знать преобразование Лапласа, чтобы пользоваться методом? Да, операционный метод - это прикладная форма преобразования Лапласа. На практике достаточно владеть таблицей соответствий, формулой изображения производной и теоремами смещения; вычислять интеграл преобразования каждый раз не нужно.

Что делать, если изображения нет в таблице? Разложить $X(p)$ на простейшие дроби - после разложения каждое слагаемое почти всегда попадает в таблицу. Если мешают комплексные корни знаменателя, их объединяют в блоки вида $\dfrac{p}{p^2 + \omega^2}$ и $\dfrac{\omega}{p^2 + \omega^2}$, дающие косинус и синус.

Коротко

Операционный метод решения дифференциального уравнения переводит задачу в пространство изображений по Лапласу, где дифференцирование становится умножением на $p$, а начальные условия встраиваются автоматически. Алгоритм из четырёх шагов - перейти к изображению, решить алгебраическое уравнение для $X(p)$, разложить результат на простейшие дроби и вернуться к оригиналу $x(t)$ - превращает решение задачи Коши в почти механическую процедуру. Метод особенно силён при задаче Коши, разрывной правой части и системах линейных уравнений с постоянными коэффициентами.