Решатель уравнения Бернулли

О Решатель уравнения Бернулли

Решатель уравнения Бернулли справляется с одним из самых известных нелинейных дифференциальных уравнений первого порядка — уравнением Бернулли y' + P(x)y = Q(x)yⁿ — и превращает классический учебный вывод в интерактивное пошаговое руководство. Он линеаризует уравнение с помощью подстановки v = y¹⁻ⁿ, строит интегрирующий множитель μ(x) и накладывает результирующую кривую в замкнутом виде на численное решение RK4 и поле направлений, чтобы вы могли видеть все детали сразу.

Что такое дифференциальное уравнение Бернулли?

Представленное Якобом Бернулли в 1695 году, уравнение Бернулли — это ОДУ первого порядка вида

Когда n = 0, уравнение уже является линейным; когда n = 1, оно является уравнением с разделяющимися переменными. Для любого другого действительного n уравнение нелинейно, но классическая подстановка v = y¹⁻ⁿ преобразует его в линейное ОДУ относительно v, которое можно решить стандартным методом интегрирующего множителя.

Шестиэтапный метод Бернулли

Начиная с y' + P(x)y = Q(x)yⁿ:

1. Разделить на yⁿ: \( y^{-n}y' + P(x)\,y^{1-n} = Q(x) \).
2. Сделать подстановку v = y¹⁻ⁿ: заметим, что \( v' = (1-n)y^{-n}y' \), следовательно \( y^{-n}y' = v'/(1-n) \).
3. Линеаризовать: \( v' + (1-n)P(x)\,v = (1-n)Q(x) \) — линейное ОДУ первого порядка относительно v.
4. Найти интегрирующий множитель: \( \mu(x) = \exp\!\left(\int (1-n)P(x)\,dx\right) \), тогда \( (\mu v)' = \mu(1-n)Q(x) \).
5. Решить для v(x): \( v(x) = \frac{1}{\mu(x)}\left[\mu(x_0)v_0 + \int_{x_0}^{x}\mu(t)(1-n)Q(t)\,dt\right] \).
6. Сделать обратную подстановку: \( y(x) = v(x)^{1/(1-n)} \).

Когда участвующие интегралы являются элементарными, вы получаете чистое решение в замкнутом виде; когда нет — калькулятор вычисляет их численно по методу Симпсона для построения кривой решения.

Автоматическая обработка частных случаев

Практический пример — n = 2, логистического типа

Рассмотрим y' + y/x = x·y² с начальным условием y(1) = 1. Здесь P(x) = 1/x, Q(x) = x и n = 2, следовательно 1 − n = −1.

1. Подстановка v = y⁻¹ = 1/y. Тогда v' = −y⁻²y' и уравнение принимает вид v' − (1/x)v = −x.
2. Интегрирующий множитель: μ(x) = exp(∫−1/x dx) = 1/x.
3. (μ·v)' = μ·(−x) = −1. Интегрируем: (1/x)·v = −x + C, т.е. v = −x² + Cx.
4. Применяем НУ: при x = 1, v = 1/1 = 1, значит 1 = −1 + C ⇒ C = 2. Следовательно v(x) = −x² + 2x.
5. Обратная подстановка: y(x) = 1/v(x) = 1/(2x − x²) = 1/(x(2 − x)).

Решение в замкнутом виде y = 1/(x(2−x)) имеет вертикальные асимптоты при x = 0 и x = 2 — именно это поле направлений делает очевидным с первого взгляда.

Как пользоваться этим калькулятором

1. Заполните конструктор уравнений. Введите P(x) и Q(x) в синие поля, а показатель степени n — в маленькое поле в верхнем индексе. Макет повторяет стандартную форму y' + P(x)y = Q(x)yⁿ.
2. Установите начальное условие (x₀, y₀) и диапазон построения [x min, x max]. Диапазон должен содержать x₀.
3. Нажмите «Решить». Калькулятор определит, является ли случай частным (n = 0 или n = 1), и покажет соответствующий вывод. В противном случае он выполнит полную шестиэтапную подстановку Бернулли с формулами в MathJax.
4. Изучите график. Оранжевая кривая — это численное решение RK4. Синяя пунктирная кривая — это замкнутая форма, вычисленная через интегрирующий множитель. Поле стрелок показывает y' повсюду, так что вы можете оценить и другие решения.
5. Скопируйте CSV с точками выборки, если хотите импортировать траекторию в другую программу.

Советы, ловушки и особые случаи

• Нецелое n требует y > 0. Решатель отмечает комбинации типа n = 1/2 при y₀ ≤ 0, где yⁿ было бы комплексным.
• y₀ = 0 часто является сингулярным. Любое уравнение Бернулли с Q ≠ 0 и n > 0 имеет тривиальное решение y ≡ 0, которое обычно не является той ветвью, которая вам нужна.
• Избегайте разрывов P(x) вблизи x₀. Выражения типа 1/x требуют x₀ ≠ 0; решатель проверяет это перед запуском.
• Большие показатели степени (|n| > 20) отклоняются для предотвращения переполнения. На практике уравнения Бернулли с такими большими n почти никогда не встречаются в реальных задачах.
• Вертикальные асимптоты. Если RK4 расходится, попробуйте сузить диапазон x до той стороны от x₀, где решение остается конечным.

Где встречаются уравнения Бернулли

• Динамика популяций — логистическое уравнение y' = ry(1 − y/K) является замаскированным уравнением Бернулли (n = 2 после перегруппировки).
• Химическая кинетика — автокаталитические реакции часто подчиняются закону y' ∝ y − y².
• Электрические цепи — некоторые цепи RL с нелинейными резисторами приводят к форме Бернулли.
• Механика жидкости — уравнения пограничного слоя после приведения подобия.
• Эпидемиологические модели — доля восприимчивых в модели SIR может быть сведена к форме Бернулли.
• Экономический рост — модель Солоу — Суона с постоянной нормой сбережений является уравнением Бернулли с n = α.

Часто задаваемые вопросы

Что такое дифференциальное уравнение Бернулли?

Уравнение Бернулли — это ОДУ первого порядка вида y' + P(x)y = Q(x)yⁿ, где P и Q — непрерывные функции, а n — любое действительное число. Это классический пример нелинейного ОДУ, которое можно преобразовать в линейное с помощью подстановки v = y¹⁻ⁿ.

Как работает подстановка v = y¹⁻ⁿ?

Умножьте исходное уравнение на y⁻ⁿ, чтобы каждый член с y стал y¹⁻ⁿ или y⁻ⁿy'. Установка v = y¹⁻ⁿ дает v' = (1−n)y⁻ⁿy'. Подстановка преобразует уравнение Бернулли в v' + (1−n)P(x)v = (1−n)Q(x), которое линейно относительно v и решается с помощью интегрирующего множителя.

Что происходит при n = 0 или n = 1?

При n = 0 уравнение уже является линейным первого порядка, поэтому подстановка не требуется. При n = 1 метод Бернулли привел бы к делению на 1 − n = 0, поэтому мы обрабатываем этот случай отдельно: уравнение сводится к y' = (Q(x) − P(x))·y, которое является уравнением с разделяющимися переменными с решением в замкнутом виде y = y₀·exp(∫(Q−P) dx).

Всегда ли уравнения Бернулли можно решить в замкнутом виде?

В принципе да, но результирующие интегралы, включающие интегрирующий множитель, могут не иметь элементарных первообразных. В таких случаях калькулятор вычисляет их численно по методу Симпсона и строит кривую решения. Сам метод всегда сводит ОДУ Бернулли к квадратурам.

Почему отрицательные y и нецелые n вызывают проблемы?

Если n не является целым числом, yⁿ определяется как exp(n·ln y) и является действительным только при y > 0. Использование отрицательного y привело бы к получению комплексного числа. Решатель отмечает эту ситуацию и запрашивает y₀ > 0 или целый показатель степени, чтобы решение оставалось действительным.

Что показывает поле направлений?

Поле направлений — это сетка крошечных отрезков касательных, угол наклона которых равен y' в данной точке (x, y). Любая кривая решения вынуждена следовать этим касательным, поэтому поле направлений позволяет увидеть качественную форму всех решений сразу, при этом начальное условие выделяет конкретную кривую.

Дополнительная литература

• Уравнение Бернулли — Википедия
• Интегрирующий множитель — Википедия
• Логистическое уравнение — Википедия
• Поле направлений — Википедия (англ.)