Решатель дифференциальных уравнений

Дифференциальное уравнение (ДУ) — это уравнение, связывающее функцию с её производными. Обыкновенное дифференциальное уравнение (ОДУ) содержит функцию одной переменной:

$F\left(x, y, y', y'', \ldots, y^{(n)}\right) = 0$

Порядок ДУ — это наивысшая производная, входящая в уравнение. Степень — это степень производной наивысшего порядка (когда уравнение полиномиально относительно производных).

ОДУ первого порядка: $y' = f(x, y)$

ОДУ второго порядка: $y'' + p(x)y' + q(x)y = g(x)$

Решение — это функция $y(x)$, удовлетворяющая уравнению на некотором интервале. Общее решение содержит произвольные постоянные (по одной на каждый порядок). Задача Коши (задача с начальными условиями) задаёт условия вроде $y(x_0) = y_0$ для определения единственного частного решения.

Дифференциальные уравнения моделируют явления реального мира: рост популяций, радиоактивный распад, системы «пружина — масса», электрические цепи, теплопроводность и течение жидкости.

Метод 1: Разделение переменных

Для уравнений вида $\frac{dy}{dx} = f(x)g(y)$:

1. Разделите: $\frac{dy}{g(y)} = f(x)\,dx$
2. Проинтегрируйте обе части: $\int \frac{dy}{g(y)} = \int f(x)\,dx$

Пример: $\frac{dy}{dx} = 2xy$ → $\frac{dy}{y} = 2x\,dx$ → $\ln|y| = x^2 + C$ → $y = Ae^{x^2}$

Метод 2: Интегрирующий множитель (линейное первого порядка)

Для $y' + P(x)y = Q(x)$ умножьте на интегрирующий множитель $\mu(x) = e^{\int P(x)\,dx}$:

$\frac{d}{dx}[\mu(x) \cdot y] = \mu(x) \cdot Q(x)$

Затем проинтегрируйте обе части, чтобы найти $y$.

Пример: $y' + 2y = e^{-x}$. Здесь $P(x) = 2$, поэтому $\mu = e^{2x}$. Умножьте: $(e^{2x}y)' = e^{x}$. Проинтегрируйте: $e^{2x}y = e^x + C$, поэтому $y = e^{-x} + Ce^{-2x}$.

Метод 3: Характеристическое уравнение (постоянные коэффициенты)

Для $ay'' + by' + cy = 0$ решите характеристическое уравнение $ar^2 + br + c = 0$:

Метод 4: Метод неопределённых коэффициентов

Для $ay'' + by' + cy = g(x)$, где $g(x)$ — многочлен, экспонента, синус, косинус или их комбинация:

1. Найдите общее решение однородного уравнения
2. Предположите форму частного решения исходя из $g(x)$
3. Подставьте и найдите коэффициенты
4. Общее решение = однородное + частное

Метод 5: Метод вариации параметров

Общий метод для $y'' + p(x)y' + q(x)y = g(x)$, когда известны однородные решения $y_1, y_2$:

$y_p = -y_1 \int \frac{y_2 g}{W}\,dx + y_2 \int \frac{y_1 g}{W}\,dx$

где $W = y_1 y_2' - y_2 y_1'$ — вронскиан.

Сравнение методов

• Забывают постоянную интегрирования: при разделении переменных постоянную нужно включить до выражения $y$, так как она влияет на окончательную форму решения.
• Неверный интегрирующий множитель: интегрирующий множитель для $y' + P(x)y = Q(x)$ — это $e^{\int P(x)\,dx}$. Убедитесь, что уравнение в стандартном виде (коэффициент при $y'$ должен быть равен 1) перед определением $P(x)$.
• Пропуск случая кратного корня: когда характеристическое уравнение имеет кратный корень $r$, второе решение — $xe^{rx}$, а не снова $e^{rx}$.
• Неверное предположение частного решения: если ваше предположение для $y_p$ уже является решением однородного уравнения, умножьте на $x$ (или $x^2$ при необходимости), чтобы получить допустимую форму.
• Игнорируют начальные условия: общее решение содержит произвольные постоянные. Применяйте начальные условия только после нахождения полного общего решения.