Дифференциальные уравнения
Степенные ряды: когда формулы кончаются
1873. Фердинанд Фробениус. Берлин. Уравнение, которое нельзя решить формулой - только бесконечной суммой. Через 150 лет эти суммы работают внутри GPS-чипов, описывают орбитали атомов водорода и рендерятся в каждой 3D-игре как spherical harmonics. Hermite polynomials, Legendre polynomials, Bessel functions - всё это степенные ряды из одного метода.
- **GPS-антенны и волноводы:** нули функций Бесселя Jₙ(x) определяют допустимые моды распространения сигнала в цилиндрических волноводах. Без них нет расчёта антенн.
- **Квантовая механика:** полиномы Эрмита (орбитали водорода s, p, d), полиномы Лагерра (радиальные волновые функции) - все через метод Фробениуса. scipy.special содержит их как `hermite`, `laguerre`.
- **3D-рендеринг:** Legendre polynomials порождают spherical harmonics - базис для ambient occlusion и environment lighting в Unreal, Unity, Blender.
Предварительные знания
Решение через степенной ряд
У уравнения нет замкнутого решения. Тупик? Нет - перестройка мышления. Если ответ нельзя записать через sin, exp, ln - можно записать через их ряды. Предположение **y = a₀ + a₁x + a₂x² + ... = Σ aₙxⁿ** - это не приближение. Для аналитических уравнений это точный ответ, просто записанный иначе.
Алгоритм: подставить ряд в уравнение, сгруппировать члены при одинаковых степенях x, приравнять каждый коэффициент нулю. Получится рекуррентное соотношение aₙ через aₙ₋₁, aₙ₋₂ - и всё решение автоматически. Точка x₀ называется **обычной**, если коэффициенты P(x) и Q(x) аналитичны в ней. В обычной точке ряд гарантированно сходится.
**Алгоритм метода степенных рядов:** 1. принять y = Σ aₙxⁿ 2. вычислить y', y'' 3. подставить в ОДУ 4. сгруппировать по степеням xⁿ 5. записать рекуррентность для aₙ 6. определить радиус сходимости.
| ОДУ | Вид ряда | Решение |
|---|---|---|
| y'' + y = 0 | Σ aₙxⁿ | sin(x), cos(x) |
| y'' + xy = 0 | Σ aₙxⁿ | Функции Эйри Ai(x), Bi(x) |
| x²y'' + xy' + (x²-ν²)y = 0 | xʳΣ aₙxⁿ (Фробениус) | Функции Бесселя Jν(x) |
| (1-x²)y'' - 2xy' + n(n+1)y = 0 | Σ aₙxⁿ | Многочлены Лежандра Pₙ(x) |
При поиске решения y = Σ aₙxⁿ уравнения y'' - y = 0, какое рекуррентное соотношение получится?
Метод Фробениуса
1873 год. Фердинанд Фробениус - Берлин. Уравнение с особой точкой в нуле: обычный степенной ряд Σ aₙxⁿ расходится прямо в начале. Решение: позволить ряду начинаться не с целой степени. **y = xʳ Σ aₙxⁿ**, где r - дробное или даже комплексное число, определяемое из индициального уравнения.
Точка x = 0 - **регулярная особая**, если P(x) имеет полюс не выше 1-го порядка, Q(x) - не выше 2-го. Подстановка y = xʳΣ aₙxⁿ при n = 0 даёт **индициальное уравнение** - квадратное в r. Его корни r₁ >= r₂ определяют два линейно независимых решения. Если r₁ - r₂ не целое - оба выражаются через xʳΣ aₙxⁿ. Если кратные или разность целая - второе решение содержит логарифм.
Если разность корней r₁ - r₂ = N (целое неотрицательное), второе решение может содержать логарифмический член: y₂ = C·y₁·ln(x) + xʳ²Σ bₙxⁿ. Проверять отдельно подстановкой - иногда C = 0 и логарифм не нужен.
**Hermite polynomials** - орбитали водорода - получаются именно методом Фробениуса для уравнения y'' - 2xy' + 2ny = 0. Квантовая механика: каждый атомный орбиталь s, p, d, f - это рекуррентное соотношение из индициального уравнения.
| Корни r₁, r₂ | Второе решение y₂ | Пример |
|---|---|---|
| r₁ - r₂ ∉ ℤ | xʳ²Σ bₙxⁿ | Бессель ν = 1/3 |
| r₁ = r₂ | y₁·ln(x) + xʳ¹Σ bₙxⁿ | Бессель ν = 0 (Y₀) |
| r₁ - r₂ = N ∈ ℤ⁺ | C·y₁·ln(x) + xʳ²Σ bₙxⁿ (C м.б. 0) | Бессель ν = 1 |
Для уравнения x²y'' + 2xy' + (x² - 2)y = 0 индициальное уравнение r(r-1) + 2r - 2 = 0. Найдите корни r₁, r₂.
Функции Бесселя - ряды в железе
Legendre polynomials описывают освещение в 3D-играх (spherical harmonics, ambient occlusion). Hermite polynomials - орбитали атомов водорода. Bessel functions - GPS-сигналы и рябь на воде. Все три - степенные ряды, полученные методом Фробениуса из одного уравнения каждый. Не «специальные функции» из учебника - реальный физический мир, записанный суммой.
**Уравнение Бесселя** порядка ν: x²y'' + xy' + (x² - ν²)y = 0. Возникает при разделении переменных в цилиндрических координатах - теплопроводность, акустика, волноводы, GPS-антенны. Решение Фробениуса при r = ν:
**Jν(x) = Σ (-1)ᵐ / (m! Γ(m+ν+1)) · (x/2)^(2m+ν)**, m от 0 до ∞. Для ν = 0: J₀(x) ≈ 1 - x²/4 + x⁴/64 - ... похожа на косинус с убывающей амплитудой. scipy.special.jn(0, x) - это оно под капотом.
Функция Бесселя второго рода Yν(x) - второе линейно независимое решение. Yν(x) → -∞ при x → 0+, поэтому в задачах на полном диске (r от 0 до R) Yν отбрасывается: физическое поле обязано быть конечным на оси. В кольцевых областях (r₁ <= r <= r₂) нужны оба решения.
**scipy.special:** `jn(n, x)` - Jₙ(x), `yn(n, x)` - Yₙ(x), `jn_zeros(n, k)` - первые k нулей Jₙ. Для нецелых ν используй `jv(nu, x)` и `yv(nu, x)`. Весь этот API - прямая реализация рядов Фробениуса.
В задаче о температуре внутри цилиндра (0 <= r <= R) с условием T(R,t) = 0, почему Yν(r) нужно отбросить из решения?
Радиус сходимости: где ряд живёт
Ряд Σ aₙ(x - x₀)ⁿ сходится в круге |x - x₀| < R, где **R = 1 / limsup|aₙ|^(1/n)**. Практическое правило: радиус сходимости ряда-решения ОДУ y'' + P(x)y' + Q(x)y = 0 в точке x₀ равен расстоянию от x₀ до **ближайшей особой точки P или Q в комплексной плоскости**. Именно в комплексной - не только на вещественной оси.
Парадокс: функция 1/(1 + x²) вещественная, гладкая на всей вещественной оси - никаких разрывов. Но её ряд Тейлора в нуле сходится только при |x| < 1. Причина: особые точки x = ±i живут в комплексной плоскости. Расстояние от нуля до i равно 1 - и это ровно радиус сходимости.
**Теорема:** если ОДУ с полиномиальными коэффициентами имеет особую точку в x = a, то радиус сходимости ряда в точке x₀ = 0 равен |a|. Для уравнения Бесселя особые точки: x = 0 (регулярная) и x = ∞ - поэтому ряд Jν(x) сходится для всех конечных x (R = ∞).
| Уравнение | Особые точки | Радиус в x₀=0 |
|---|---|---|
| Бесселя x²y''+xy'+(x²-ν²)y=0 | x = 0 (рег.), ∞ | R = ∞ |
| Лежандра (1-x²)y''-2xy'+n(n+1)y=0 | x = ±1 | R = 1 |
| Гаусса x(1-x)y''+(c-(a+b+1)x)y'-aby=0 | x = 0, 1, ∞ | R = 1 |
| Эрмита y''-2xy'+2ny=0 | Нет (целая функция) | R = ∞ |
Ряд для y(x) разложен в точке x₀ = 0. Коэффициенты P(x) и Q(x) имеют особые точки при x = 3 и x = 2+i. Каков радиус сходимости?
Ключевые идеи
- **Степенной ряд y = Σaₙxⁿ** решает ОДУ в обычных точках: подставить, приравнять коэффициенты, получить рекуррентность.
- **Метод Фробениуса y = xʳΣaₙxⁿ** работает в регулярных особых точках. Показатель r - из индициального уравнения.
- **Функции Бесселя Jν(x)** - решение уравнения Бесселя в цилиндрической геометрии. Их нули задают собственные частоты и моды антенн.
- **Радиус сходимости** = расстояние до ближайшей особой точки в комплексной плоскости - не только на вещественной оси.
- **Hermite, Legendre, Bessel** - три имени, три домена применения, один метод: ряды Фробениуса.
Связанные темы
Степенные ряды открывают путь к специальным функциям и прикладным задачам:
- Системы ОДУ — Матричная экспонента - тоже степенной ряд: e^(At) = I + At + (At)²/2! + ...
- Уравнение теплопроводности — Решение в цилиндрических координатах требует Jν(x) из этого урока
- Метод конечных элементов — Базисные функции FEM - полиномы Лежандра/Лагранжа, связанные с рядами
Вопросы для размышления
- Функции Бесселя Jν(x) при больших x ведут себя как затухающие косинусы: Jν(x) ≈ √(2/πx) cos(x - νπ/2 - π/4). Как это согласуется с тем, что уравнение Бесселя - это почти уравнение y'' + y = 0?
- Метод Фробениуса даёт второе решение с логарифмом при кратных корнях. В каких физических задачах логарифмическая особенность в нуле допустима, а в каких нет?
- Почему расстояние до комплексной особой точки ограничивает сходимость вещественного ряда? Объясни на примере 1/(1+x²) и её особых точек x = ±i.
Связанные уроки
- calc-02-series-intro — Сходимость рядов - фундамент метода
- de-03 — Системы ОДУ - предшественник Фробениуса
- de-07 — Теплопроводность в цилиндре требует Jν(x)
- nm-13 — Конечные разности как альтернатива рядам
- de-13 — МКЭ использует полиномы Лежандра из тех же рядов
- calc-16-taylor