Характеристическое уравнение: комплексные корни ОДУ

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

Характеристическое уравнение и его дискриминант

Рассмотрим однородное линейное ОДУ второго порядка с постоянными коэффициентами:

$y'' + py' + qy = 0.$

Решение ищут в виде $y = e^{\lambda x}$. Подставив в уравнение и сократив на $e^{\lambda x} \neq 0$, получают характеристическое уравнение:

$\lambda^2 + p\lambda + q = 0.$

Это обычное квадратное уравнение. Его дискриминант $D = p^2 - 4q$ определяет тип корней:

Нас интересует случай $D < 0$. Он возникает, когда $q > p^2/4$, то есть «возвращающая сила» $q$ достаточно велика по сравнению с «затуханием» $p$. Именно этот случай соответствует колебательным процессам - маятнику с малым трением, LC-контуру, упругой балке.

Важно понимать, почему вообще пробуют $y = e^{\lambda x}$. Функция $e^{\lambda x}$ особенная: её производные пропорциональны ей самой, $y' = \lambda e^{\lambda x}$ и $y'' = \lambda^2 e^{\lambda x}$. Это единственный класс функций, для которых все члены уравнения $y'' + py' + qy$ имеют одинаковый «тип», и после подстановки уравнение распадается на множители: $e^{\lambda x}(\lambda^2 + p\lambda + q) = 0$.

Формулы для комплексных корней

При $D < 0$ квадратный корень из дискриминанта мнимый. Решая квадратное уравнение по формуле:

$\lambda_{1,2} = \frac{-p \pm \sqrt{p^2 - 4q}}{2} = \frac{-p \pm \sqrt{-|D|}}{2} = \frac{-p}{2} \pm \frac{\sqrt{|D|}}{2}\,i.$

Обозначают вещественную и мнимую части отдельно:

$$\alpha = -\frac{p}{2}, \qquad \beta = \frac{\sqrt{4q - p^2}}{2} = \frac{\sqrt{|D|}}{2} > 0.$$

Тогда корни записываются кратко: $\lambda_{1,2} = \alpha \pm \beta i$.

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

Переход от комплексных корней к вещественному решению

С формальным записями $C_1 e^{\lambda_1 x} + C_2 e^{\lambda_2 x}$ работать неудобно, потому что $C_1, C_2$ пришлось бы брать комплексными. Используют формулу Эйлера:

$e^{(\alpha \pm \beta i)x} = e^{\alpha x}(\cos \beta x \pm i \sin \beta x).$

После линейной комбинации и переобозначения констант получают вещественную форму общего решения:

$$\boxed{y = e^{\alpha x}\bigl(C_1 \cos \beta x + C_2 \sin \beta x\bigr),}$$

где $C_1, C_2$ - произвольные вещественные константы, определяемые из начальных условий. Это и есть рабочая формула для всех задач с комплексными корнями.

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

На графике видно три характерных случая в зависимости от знака $\alpha$:

• $\alpha < 0$ (т.е. $p > 0$): затухающие колебания - амплитуда убывает как $e^{\alpha x}$;
• $\alpha = 0$ (т.е. $p = 0$): незатухающие колебания - амплитуда постоянна;
• $\alpha > 0$ (т.е. $p < 0$): нарастающие колебания - амплитуда неограниченно растёт.

Физический смысл: что означают alpha и beta

Параметр $\alpha$ отвечает за характер изменения амплитуды: это скорость нарастания или затухания. В механических системах $p$ - коэффициент вязкого трения, поэтому $p > 0$ всегда соответствует физически реальному затуханию. Параметр $\beta$ - это циклическая частота колебаний (в единицах аргумента $x$). Если $x = t$ (время), то $\beta = \omega$ - угловая частота, а период колебаний:

$T = \frac{2\pi}{\beta}.$

В ОДУ вида $y'' + 2\delta y' + \omega_0^2 y = 0$ (уравнение свободных колебаний с затуханием) параметры напрямую соответствуют физике: $\delta$ - коэффициент затухания, $\omega_0$ - собственная частота. Формулы перехода: $p = 2\delta$, $q = \omega_0^2$, поэтому $\alpha = -\delta$ и $\beta = \sqrt{\omega_0^2 - \delta^2}$ - уменьшенная (реальная) частота колебаний.

Отсюда следует практически важное условие: комплексные корни (колебательный режим) существуют тогда и только тогда, когда $\delta < \omega_0$, то есть затухание меньше собственной частоты. Это называют режимом недодемпфирования (underdamping). Если $\delta = \omega_0$, система критически демпфирована (кратный корень); если $\delta > \omega_0$ - передемпфирована (два вещественных корня, никаких колебаний). Калькулятор выше наглядно показывает все три режима: уменьшайте $q$ при фиксированном $p$ - и кривая из колебательной перейдёт в монотонную.

Пример решения задачи

Разберём стандартную постановку: найти общее решение уравнения $y'' + 2y' + 5y = 0$.

Шаг 1. Составляем характеристическое уравнение: $\lambda^2 + 2\lambda + 5 = 0$.

Шаг 2. Вычисляем дискриминант: $D = 4 - 20 = -16 < 0$. Корни комплексные.

Шаг 3. Находим вещественную и мнимую части:

$$\alpha = -\frac{2}{2} = -1, \qquad \beta = \frac{\sqrt{|-16|}}{2} = \frac{4}{2} = 2.$$

Корни: $\lambda_{1,2} = -1 \pm 2i$.

Шаг 4. Записываем общее решение:

$$y = e^{-x}(C_1 \cos 2x + C_2 \sin 2x).$$

Проверка: $\alpha = -1 < 0$, значит колебания затухающие. Период $T = 2\pi/2 = \pi$.

Убедиться в правильности ответа можно, подставив обратно в уравнение. Берём $y = e^{-x}(C_1 \cos 2x + C_2 \sin 2x)$, вычисляем $y'$ и $y''$, складываем $y'' + 2y' + 5y$ и проверяем, что результат тождественно равен нулю. Это занимает время, но при любых сомнениях - стоит.

Задача Коши с начальными условиями

Теперь добавим условия: $y(0) = 1$, $y'(0) = 0$.

Из $y(0) = e^0(C_1 \cdot 1 + C_2 \cdot 0) = C_1 = 1$ сразу получаем $C_1 = 1$.

Для $C_2$ нужна производная:

$$y' = e^{-x}\bigl[(-C_1 + 2C_2)\cos 2x + (-2C_1 - C_2)\sin 2x + \ldots\bigr].$$

Точнее, по правилу произведения:

$$y' = -e^{-x}(C_1 \cos 2x + C_2 \sin 2x) + e^{-x}(-2C_1 \sin 2x + 2C_2 \cos 2x).$$

При $x = 0$: $y'(0) = -C_1 + 2C_2 = 0$, откуда $C_2 = C_1/2 = 1/2$.

Ответ: $y = e^{-x}\!\left(\cos 2x + \frac{1}{2}\sin 2x\right)$.

Это можно записать в амплитудно-фазовой форме: $y = A e^{-x} \cos(2x - \varphi)$, где $A = \sqrt{C_1^2 + C_2^2} = \sqrt{1 + 1/4} = \sqrt{5}/2$ и $\tan\varphi = C_2/C_1 = 1/2$. Такая запись удобна, когда нужно сразу видеть начальную амплитуду колебаний.

Связь с уравнениями высшего порядка и системами

Метод характеристического уравнения распространяется на ОДУ любого порядка $n$: подстановка $y = e^{\lambda x}$ превращает его в алгебраическое уравнение степени $n$. Среди его корней комплексные всегда идут парами (теорема о комплексно-сопряжённых корнях для уравнений с вещественными коэффициентами). Каждой паре $\alpha \pm \beta i$ отвечает двумерное пространство решений - точно такой же «блок» $e^{\alpha x}(C_k \cos \beta x + C_{k+1} \sin \beta x)$.

Если пара является кратным корнем кратности $m$, решение усложняется: каждый «блок» умножается на полином степени $m-1$:

$$e^{\alpha x}\bigl[(C_1 + C_2 x + \ldots + C_m x^{m-1})\cos \beta x + (D_1 + D_2 x + \ldots + D_m x^{m-1})\sin \beta x\bigr].$$

На практике в задачах второго порядка кратные комплексные корни не встречаются (это потребовало бы специального выбора $p$ и $q$), но в задачах четвёртого порядка (балки, уравнение Эйлера-Бернулли) - уже возможны.

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

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

• Забыть проверить знак дискриминанта. Формулу $y = e^{\alpha x}(C_1 \cos \beta x + C_2 \sin \beta x)$ применяют только при $D < 0$. При $D > 0$ и $D = 0$ формулы другие.
• Перепутать $\alpha$ и $\beta$. Помните: $\alpha = -p/2$ - вещественная часть корня (влияет на амплитуду), $\beta = \sqrt{|D|}/2$ - мнимая часть (влияет на частоту).
• Подставить $\lambda$ в $y = e^{\lambda x}$ целиком вместо перехода к вещественной форме. Получается формально правильно, но работать с комплексными $C_1, C_2$ на практике труднее.
• Ошибиться в производной при задаче Коши. При дифференцировании $e^{\alpha x}(C_1 \cos \beta x + C_2 \sin \beta x)$ нужно применить правило произведения аккуратно - по обеим функциям от $x$.
• Путать $\beta$ с угловой частотой уравнения. В уравнении $y'' + \omega_0^2 y = 0$ (без члена $py'$) имеем $p = 0$, $q = \omega_0^2$, поэтому $\beta = \omega_0$ - это совпадение только при нулевом затухании.

FAQ

Что такое характеристическое уравнение ОДУ? Характеристическое уравнение получают, подставив $y = e^{\lambda x}$ в однородное линейное ОДУ с постоянными коэффициентами и сократив на $e^{\lambda x}$. Для $y'' + py' + qy = 0$ оно имеет вид $\lambda^2 + p\lambda + q = 0$. Его корни определяют, из каких функций складывается общее решение.

Почему при комплексных корнях общее решение остаётся вещественным? Комплексно-сопряжённые корни $\alpha \pm \beta i$ дают пару функций $e^{(\alpha+\beta i)x}$ и $e^{(\alpha-\beta i)x}$. Их линейные комбинации по формуле Эйлера сводятся к $e^{\alpha x}\cos \beta x$ и $e^{\alpha x}\sin \beta x$ - вещественным функциям. Произвольные вещественные $C_1, C_2$ обеспечивают полный набор решений.

Как связаны комплексные корни и физика колебаний? Мнимая часть $\beta$ задаёт частоту колебаний, а вещественная часть $\alpha$ - их характер. При $\alpha < 0$ колебания затухают (диссипативная система), при $\alpha = 0$ - гармонические (идеальный колебательный контур или маятник без трения), при $\alpha > 0$ - нарастают (неустойчивая система).

Коротко

Комплексные корни характеристического уравнения $\lambda_{1,2} = \alpha \pm \beta i$ возникают при $D = p^2 - 4q < 0$. Вещественная часть $\alpha = -p/2$ отвечает за нарастание или затухание амплитуды, мнимая часть $\beta = \sqrt{4q-p^2}/2$ - за частоту колебаний. Общее решение записывается по формуле $y = e^{\alpha x}(C_1 \cos \beta x + C_2 \sin \beta x)$ с произвольными вещественными $C_1, C_2$, которые определяются из начальных условий.