Математика / Начала анализа
Метод Тейлора второго порядка для задачи Коши
Метод Тейлора второго порядка добавляет к шагу Эйлера член с полной производной правой части. Для скалярного уравнения y′=f(t,y) вторая производная решения равна fₜ+fᵧf, вычисленной в начале шага.
Формула
Для y′=y², y₀=1, h=0,2 показано разложение одного шага.
Обозначения
- $yₙ,tₙ$
- приближённое состояние и время начала шага, единицы y и t
- $h$
- шаг времени, единица t
- $fₙ$
- f(tₙ,yₙ), первая производная решения, единица y / единица t
- $fₜ$
- частная производная f по времени при фиксированном y, единица y / единица t²
- $fᵧ$
- частная производная f по состоянию при фиксированном t, 1 / единица t
Условия применения
Решается скалярная гладкая задача Коши y′=f(t,y). Производные fₜ и fᵧ известны и определены в узле. Для второго глобального порядка требуется достаточная гладкость решения и устойчивое накопление ошибки на фиксированном интервале.
Ограничения
Численное дифференцирование шумной правой части добавляет собственную ошибку. Формула не имеет встроенного контроля шага и остаётся явной, поэтому жёсткость ограничивает устойчивость. Для системы fᵧ становится матрицей Якоби, и нужно матрично-векторное произведение.
Подробное объяснение
Разложение точного решения за малое время даёт y(t+h)=y(t)+hy′(t)+h²y″(t)/2+O(h³). Первый член после состояния совпадает с Эйлером. Второй учитывает изменение наклона внутри шага и поэтому исправляет главный недостающий вклад кривизны.
Для y′=f(t,y(t)) нельзя получать y″ только частным дифференцированием по t. По правилу цепочки y″=fₜ+fᵧy′=fₜ+fᵧf. В автономном уравнении fₜ=0, но второй член обычно остаётся: изменение состояния само меняет производную.
Все части приращения имеют единицу y. В частности, fᵧ имеет единицу обратного времени, а произведение fᵧf — единицу y/t², как и fₜ. Такая проверка быстро выявляет потерянный множитель f или лишнюю степень шага.
На линейном тесте y′=λy множитель шага равен 1+z+z²/2, z=hλ. Он совпадает с явными RK2, однако на нелинейном уравнении конечные шаги могут отличаться. Метод Тейлора использует производные в старом узле; методы Рунге — Кутты добиваются того же порядка дополнительными вычислениями f.
Как пользоваться формулой
- Запишите правую часть как функцию двух независимых аргументов и найдите частные производные по времени и состоянию.
- Вычислите f, fₜ и fᵧ в одном исходном узле; составьте полную вторую производную как fₜ+fᵧf.
- Прибавьте к старому состоянию линейное и квадратичное приращения, сохраняя полный шаг h в обоих слагаемых.
- Повторите на более мелкой сетке и сравните при одинаковом конечном времени; отдельно проверьте устойчивость модели.
Историческая справка
Метод продолжает классическую идею разложения гладкой функции Тейлора, применённую непосредственно к решению дифференциального уравнения. В численном анализе XIX–XX веков такой путь дал естественный ориентир для определения порядка: сколько членов точного приращения воспроизводит дискретный переход.
Современное изложение Fundamentals of Numerical Computation начинает вывод методов Рунге — Кутты именно со второго члена Тейлора и правила цепочки. Оно показывает причину перехода к многостадийным схемам: аналитические старшие производные быстро усложняются, а дополнительные обращения к f могут воспроизвести нужные члены без их явного вычисления. Поэтому метод Тейлора сохраняет учебную и проверочную ценность даже там, где основной расчёт выполняют другим решателем.
Историческая линия формулы
Название указывает на использование разложения Тейлора для решения ОДУ. Конкретная формула второго порядка следует из правила цепочки для f(t,y(t)). Это способ временного интегрирования, который следует отличать от общего утверждения о разложении произвольной функции.
Пример
В безразмерной задаче y′=y², y(0)=1 возьмём h=0,2. Здесь f=y²=1, fₜ=0 и fᵧ=2y=2. Значит, y″=fᵧf=2. Шаг Тейлора: y₁=1+0,2·1+(0,2²/2)·2=1,24. Точное решение y(t)=1/(1−t) даёт y(0,2)=1,25, абсолютная ошибка 0,01. Только Эйлер дал бы 1,2 с ошибкой 0,05. Ненулевое y″ возникает целиком из изменения y: отсутствие явного t в правой части не делает решение линейным. Проверка единиц проводится на безразмерной учебной модели; для физического y′=ay² коэффициент a должен иметь единицу 1/(y·t). Проверку шага проводят до особенности t=1.
Частая ошибка
Не считайте fₜ полной второй производной решения. Не забывайте множитель f в цепочке fᵧf. Производные берут в исходном узле, а не в уже обновлённом состоянии. Высокий порядок разложения не снимает ограничение явной схемы по устойчивости. Сначала проверяйте правило цепочки.
Практика
Задачи с решением
Зависимость от времени
Условие. Для y′=t+y, y(0)=1, h=0,1 найдите шаг.
Решение. f=1, fₜ=1, fᵧ=1; y″=2. Поэтому y₁=1+0,1+0,01=1,11.
Ответ. 1,11.
Линейное затухание
Условие. Для y′=−2y, y₀=1, h=0,25 найдите шаг.
Решение. f=−2, fᵧ=−2; y″=4. y₁=1−0,5+(0,25²/2)·4=0,625.
Ответ. 0,625.
Калькулятор
Расчёт: Новое y
Один скалярный шаг; f, fₜ и fᵧ задаются в исходном узле. Гладкость и устойчивость проверяют отдельно.
- Один скалярный шаг; f, fₜ и fᵧ задаются в исходном узле. Гладкость и устойчивость проверяют отдельно.
Дополнительные источники
- https://fncbook.com/rk/
- https://fncbook.com/euler/
Связанные формулы
Математика
Формула Тейлора с остаточным членом
Формула Тейлора с остаточным членом: формула f(x)=\sum_{k=0}^{n}\frac{f^{(k)}(a)}{k!}(x-a)^k+R_n(x),\quad R_n(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!}(x-a)^{n+1} помогает разобрать функцию через пределы, производные, интегралы или ряды. В тексте есть условия, пример, ошибки и проверка результата.
Математика
Явный метод Эйлера: решение задачи Коши по шагам
Явный метод Эйлера приближённо решает уравнение y′ = f(t,y), продолжая решение по наклону в начале каждого шага. При гладком решении глобальная ошибка имеет первый порядок по шагу.
Математика
Классический метод Рунге — Кутты четвёртого порядка
Классический RK4 вычисляет четыре наклона за шаг и объединяет их с весами 1:2:2:1. Для гладкой задачи Коши он даёт четвёртый порядок общей ошибки без явного вычисления высших производных.