Как решить дифференциальное уравнение с помощью ряда: короткий ответ
Чтобы решить дифференциальное уравнение с помощью ряда, ищи решение в виде степенного ряда: y = a₀ + a₁x + a₂x² + a₃x³ + … Подставь этот ряд и его производные в уравнение, собери слагаемые при одинаковых степенях x и приравняй коэффициенты. Получится рекуррентная формула, по которой числа aₙ находятся одно за другим. Первые из них задают начальные условия.
Зачем решать дифуры рядами, если есть готовые методы?
Знакомая ситуация: открываешь задание, а там уравнение, которое не разделяется, не однородное и не линейное с постоянными коэффициентами. Перебрал все типы из шпаргалки, ни один не подошёл, а в условии написано: «Решите с помощью степенного ряда». Паника понятна: на этот метод в курсе отводят одну пару, и выглядит он совсем не похожим на остальные.
На самом деле идея простая. Многие дифференциальные уравнения, например y″ − xy = 0 или y′ = x² + y², вообще не решаются в элементарных функциях: ни формулой с синусами, ни формулой с экспонентами ответ не записать. Но решение существует, и его можно представить бесконечной суммой степеней x. Ряд даёт ответ с любой нужной точностью, а именно это и требуется в прикладных задачах.
Метод к тому же универсален. Одни и те же шаги работают для ОДУ первого и второго порядка, линейных и нелинейных. Неудобство одно: считать приходится много, и легко запутаться в индексах. Поэтому ниже сначала алгоритм, а потом три полностью разобранных примера, где видно, в каких местах обычно ошибаются. Если разберёшь их по одному на бумаге, остальные задачи из контрольной пойдут по той же схеме, меняются только числа в рекуррентной формуле.
Когда решение ДУ можно записать степенным рядом?
Степенной ряд — это сумма вида a₀ + a₁(x − x₀) + a₂(x − x₀)² + …, где x₀ — точка, вокруг которой ищем решение. Чаще всего берут x₀ = 0, и тогда ряд выглядит как a₀ + a₁x + a₂x² + … Коэффициенты aₙ — обычные числа, их мы и ищем. Начальные условия сразу дают первые два: a₀ = y(0) и a₁ = y′(0).
Решение разлагается в ряд не всегда, но в большинстве учебных задач — да. Если правая часть уравнения y′ = f(x, y) раскладывается в степенной ряд по степеням (x − x₀) и (y − y₀) около начальной точки, то решение задачи Коши с условием y(x₀) = y₀ раскладывается в степенной ряд в некоторой окрестности x₀: это теорема Коши об аналитическом решении. Для линейного уравнения y″ + p(x)y′ + q(x)y = 0 достаточно, чтобы p(x) и q(x) раскладывались в ряды, сходящиеся при |x − x₀| < R: тогда любое решение раскладывается в ряд, сходящийся по крайней мере там же.
Что с областью сходимости? Для линейных уравнений, где при старшей производной стоит единица, а остальные коэффициенты — многочлены, как в первых двух примерах ниже, ряды сходятся при всех x. Для нелинейных уравнений радиус сходимости заранее неизвестен и часто конечен: решение третьего примера уходит в бесконечность примерно при x ≈ 2, и ряд дальше не работает. Это нормально: метод и задуман для приближённого описания решения рядом с заданной точкой.
А если начальные условия заданы не в нуле, например y(2) = 1? Тогда сделай замену t = x − 2. После неё условие переезжает в точку t = 0, а производные по t совпадают с производными по x. Только не забудь заменить x на t + 2 в коэффициентах самого уравнения. Решай уравнение относительно t обычным образом, а в конце вернись к x, подставив t = x − 2. Так что случай x₀ ≠ 0 сводится к уже знакомому и новых приёмов не требует.
Как найти коэффициенты ряда: алгоритм по шагам
Шаг первый: запиши y = a₀ + a₁x + a₂x² + a₃x³ + … и продифференцируй почленно: y′ = a₁ + 2a₂x + 3a₃x² + 4a₄x³ + …, а y″ = 2a₂ + 6a₃x + 12a₄x² + 20a₅x³ + … Шаг второй: подставь всё в уравнение. Шаг третий: раскрой скобки и собери слагаемые при одинаковых степенях x — отдельно при x⁰, при x¹, при x² и так далее.
Шаг четвёртый: приравняй коэффициенты при каждой степени x в левой и правой частях. Так получится цепочка равенств, а из неё рекуррентная формула, выражающая следующий коэффициент через предыдущие. Шаг пятый: подставь a₀ и a₁ из начальных условий и вычисли несколько коэффициентов. Шаг шестой: запиши ответ в виде первых членов ряда и, если получится, общего члена.
Особое внимание сдвигу индексов: после дифференцирования степень x понижается, поэтому коэффициент при xⁿ в y′ равен (n + 1)aₙ₊₁, а в y″ — (n + 2)(n + 1)aₙ₊₂. Именно здесь теряют единицу. Результат можно сверить в калькуляторе дифференциальных уравнений Prepodoff: он умеет решать рядами и показывает ход решения, так что ты увидишь, на каком коэффициенте твоя цепочка разошлась с верной.
Есть и второй путь, через производные. Коэффициенты ряда Тейлора равны aₙ = y⁽ⁿ⁾(0)/n!, то есть a₀ = y(0), a₁ = y′(0), a₂ = y″(0)/2, a₃ = y‴(0)/6. Значения производных находишь, дифференцируя само уравнение и подставляя x = 0. Этот способ удобен для нелинейных уравнений, где ряды неудобно перемножать, а нужны лишь два-три первых члена. Для линейных с многочленами обычно быстрее метод неопределённых коэффициентов.
Пример 1: как решить y′ = y + x при y(0) = 1?
Ищем решение в виде y = a₀ + a₁x + a₂x² + a₃x³ + … Из условия y(0) = 1 получаем a₀ = 1. Производная равна y′ = a₁ + 2a₂x + 3a₃x² + 4a₄x³ + 5a₅x⁴ + … Подставляем в уравнение y′ = y + x и сравниваем коэффициенты при одинаковых степенях x. Слагаемое x в правой части добавляется только к коэффициенту при первой степени.
При x⁰: a₁ = a₀ = 1. При x¹: 2a₂ = a₁ + 1 = 2, значит a₂ = 1. При x²: 3a₃ = a₂ = 1, значит a₃ = 1/3. При x³: 4a₄ = a₃, значит a₄ = 1/12. При x⁴: 5a₅ = a₄, значит a₅ = 1/60. Начиная со второй степени, правило одно: (n + 1)aₙ₊₁ = aₙ.
Ответ: y = 1 + x + x² + x³/3 + x⁴/12 + x⁵/60 + … Проверим. Это линейное уравнение, и его точное решение — y = 2eˣ − x − 1. Разложим: 2eˣ = 2 + 2x + x² + x³/3 + x⁴/12 + x⁵/60 + …, вычтем x + 1 и получим тот же ряд. Подстановкой тоже сходится: y′ = 1 + 2x + x² + x³/3 + x⁴/12 + …, и столько же даёт y + x во всех членах до x⁴ включительно.
Если точного решения нет и сверить не с чем, проверяй иначе: подставь частичную сумму в уравнение и посмотри на остаток. Здесь y′ − y − x для суммы до x⁵ даёт только −x⁵/60: остаток начинается с высокой степени. Значит, все коэффициенты до этой степени найдены верно. Так можно контролировать себя в любом примере, даже когда общий член ряда не угадывается.
Пример 2: как решить y″ − xy = 0 при y(0) = 1, y′(0) = 0?
Это уравнение не решается в элементарных функциях, поэтому ряд здесь единственный разумный путь. Берём y = a₀ + a₁x + a₂x² + …, из условий a₀ = 1 и a₁ = 0. Тогда y″ = 2a₂ + 6a₃x + 12a₄x² + 20a₅x³ + 30a₆x⁴ + …, а xy = a₀x + a₁x² + a₂x³ + a₃x⁴ + …
Приравниваем коэффициенты у y″ и xy. При x⁰: 2a₂ = 0, значит a₂ = 0. При x¹: 6a₃ = a₀ = 1, значит a₃ = 1/6. При x²: 12a₄ = a₁ = 0, поэтому a₄ = 0. При x³: 20a₅ = a₂ = 0. При x⁴: 30a₆ = a₃ = 1/6, значит a₆ = 1/180. Общая формула: (n + 2)(n + 1)aₙ₊₂ = aₙ₋₁ при n ≥ 1.
Ответ: y = 1 + x³/6 + x⁶/180 + … Проверка: y″ = x + x⁴/6 + …, и xy = x + x⁴/6 + …, коэффициенты совпадают до x⁴. Заметь закономерность: ненулевые коэффициенты стоят только при степенях, кратных трём, потому что формула связывает aₙ₊₂ с aₙ₋₁, а шаг по индексу равен трём. Это позволяет сразу видеть нулевые коэффициенты и не считать их.
Обрати внимание на условие n ≥ 1. Слагаемое xy начинается со степени x¹, поэтому при x⁰ в правой части стоит ноль, и уравнение 2a₂ = 0 стоит особняком от общей формулы. Если подставить в формулу n = 0, справа окажется a₋₁, а такого коэффициента нет. Эту первую строку запоминай отдельно: пропуск её даёт неверное a₂ и портит весь ряд.
Пример 3: как решить нелинейное y′ = x² + y² при y(0) = 0?
Здесь общую рекуррентную формулу записать неудобно: y² — это произведение рядов, и коэффициенты складываются из попарных произведений. Поэтому идём по степеням последовательно, находя коэффициенты по очереди и каждый раз проверяя, какие из них уже известны и что получается справа. Из условия a₀ = 0. Берём y = a₁x + a₂x² + a₃x³ + … и подставляем в y′ = x² + y².
При x⁰: a₁ = 0, потому что в правой части свободного члена нет. При x¹: 2a₂ = 0, значит a₂ = 0. При x²: 3a₃ = 1, значит a₃ = 1/3. Теперь y начинается с x³, а y² — с x⁶. Поэтому при x³, x⁴ и x⁵ справа нули, и a₄ = a₅ = a₆ = 0. При x⁶: 7a₇ = (1/3)² = 1/9, значит a₇ = 1/63.
Ответ: y = x³/3 + x⁷/63 + … Проверка: y′ = x² + x⁶/9 + …, и x² + y² = x² + x⁶/9 + … Следующий член получается из условия 11a₁₁ = 2 · (1/3) · (1/63) = 2/189, то есть a₁₁ = 2/2079. Решение этого уравнения нельзя записать элементарной формулой, поэтому ряд здесь и есть рабочий ответ.
Сверим результат вторым способом, через производные. Дифференцируем уравнение: y″ = 2x + 2yy′, затем y‴ = 2 + 2(y′)² + 2yy″. При x = 0 имеем y = 0 и y′ = 0, откуда y″(0) = 0 и y‴(0) = 2. Тогда a₂ = 0 и a₃ = 2/6 = 1/3, что совпадает с коэффициентами, найденными подстановкой. Два независимых способа с одним ответом — хороший признак, что ошибок нет.
Что делать, если запутался в индексах или нужно общее решение?
Чаще всего ошибаются в трёх местах: теряют сдвиг индекса, забывают слагаемое правой части (как x в первом примере) и неверно учитывают начальные условия. Помогает простая привычка: прежде чем искать общий член, выпиши вручную первые четыре-пять коэффициентов. Нужен ли общий член в ответе, зависит от требований твоего преподавателя, часто достаточно нескольких первых членов.
Можно ли получить общее решение рядом? Да. Для уравнения второго порядка оставь a₀ и a₁ произвольными, как C₁ и C₂. В примере 2 получится: y = a₀(1 + x³/6 + x⁶/180 + …) + a₁(x + x⁴/12 + x⁷/504 + …). Для второй скобки: 12a₄ = a₁ и 42a₇ = a₄, отсюда a₄ = a₁/12 и a₇ = a₁/504. Две скобки — это два независимых решения.
Решай задачу сам, а потом сверяйся. Бесплатный калькулятор дифференциальных уравнений Prepodoff (/tools/differential-equations) работает онлайн без регистрации, решает диффуры рядами пошагово, проверяет ответ подстановкой и строит график. Используй его как способ проверить свои коэффициенты и посмотреть ход решения, а не как готовый ответ для списывания. Нужен план подготовки к экзамену по ДУ — в «Моём плане» в Telegram его составит ИИ-тренер.
