Математика / Начала анализа
Неявная средняя точка: шаг и сохранение структуры
Неявная средняя точка вычисляет производную в среднем старого и нового состояний. Новый уровень входит в правую часть, поэтому шаг находят решением уравнения. Схема второго порядка обратима по времени при согласованном решении.
Формула
q₀=1, v₀=0, h=0,2. Новые координаты найдены совместно.
Обозначения
- $yₙ,yₙ₊₁$
- старое и искомое новое состояние, единица y
- $tₙ$
- начальный момент, единица t
- $h$
- полный шаг, единица t
- $f$
- правая часть дифференциального уравнения, единица y / единица t
Условия применения
Требуется гладкая задача Коши и возможность решить неявное уравнение для нового состояния. Выбирают ветвь решения, непрерывно связанную со старым состоянием при h→0. Для гамильтоновых свойств нужны соответствующая система и достаточно точное решение шага.
Ограничения
A-устойчивость на линейном тесте не означает сильного подавления жёсткого затухания: множитель стремится к −1. Нелинейное уравнение шага может иметь несколько корней или плохо сходиться. При неточном решении симметрия и сохранение инвариантов могут нарушаться.
Подробное объяснение
Средняя точка здесь определяется не прогнозным полушагом, а арифметическим средним двух принятых состояний. Поэтому неизвестное yₙ₊₁ находится внутри функции f. Даже при одном обращении к f в записи схемы её реализация может требовать многих итераций решения уравнения.
При замене h на −h и перестановке старого и нового состояний формула возвращается к себе. Такая симметрия объясняет обратимость временного шага при согласованном выборе корня. Явный прогноз середины не обладает тем же свойством, несмотря на похожее название.
Для y′=λy получается множитель (1+z/2)/(1−z/2). Он совпадает с неявной трапецией на этом тесте. Для нелинейного f среднее значений на концах и значение в среднем состоянии различаются, поэтому две неявные схемы не тождественны.
У линейного осциллятора q′=v, v′=−q изменение q²+v² равно произведению сумм координат на их приращения. Подстановка средних наклонов даёт взаимное сокращение, так что квадратичная энергия сохраняется. Для общего нелинейного потенциала вместо точной энергии исследуют структуру и долгосрочное поведение ошибки.
Как пользоваться формулой
- Составьте уравнение нового состояния с правой частью в среднем времени и среднем старого и нового состояний.
- Выберите начальное приближение к новому состоянию и решите уравнение с контролем невязки, а не только числа итераций.
- Проверьте ветвь, связанную с малым шагом, и сравните невязку с требуемой точностью временной дискретизации.
- Для динамики измеряйте инварианты и фазу отдельно; уменьшите шаг и проверьте обратный переход с отрицательным шагом.
Историческая справка
Неявная средняя точка принадлежит симметричным методам Рунге — Кутты и интерпретируется как одноузловая гауссова схема. В современной теории численного интегрирования её выделяют потому, что временная симметрия и симплектичность дают иной набор свойств, чем у явных схем того же второго порядка.
В курсах геометрической интеграции конца XX–XXI века метод сравнивают с Эйлером, Верле и трапецией на задачах механики. Лекции Хайрера показывают сохранение геометрической структуры, а современные исследования различают точное сохранение исходного инварианта и сохранение структуры, объясняющей долгосрочное поведение. Поэтому историческая учебная роль метода связана не только с порядком ошибки, но и с качеством воспроизведения движения.
Историческая линия формулы
Неявная средняя точка — одношаговая гауссова схема второго порядка из семейства Рунге — Кутты. Имя геометрического узла определяет место вычисления правой части. Оно не объединяет метод с явной средней точкой, у которой иная таблица стадий и иные свойства.
Пример
Для безразмерного осциллятора q′=v, v′=−q возьмём q₀=1, v₀=0 и h=0,2. Уравнения шага: q₁=1+0,1v₁ и v₁=−0,1(1+q₁). Подстановка даёт 1,01q₁=0,99, поэтому q₁=99/101≈0,980198020, v₁=−20/101≈−0,198019802. Проверка энергии: q₁²+v₁²=(9801+400)/10201=1, как в начале. Точные cos(0,2)≈0,980066578 и −sin(0,2)≈−0,198669331 отличаются: сохранение энергии не устраняет фазовую ошибку. Все числа относятся к одному решению линейной системы шага, а не к последовательной подстановке уже обновлённой координаты. Шаг безразмерен.
Частая ошибка
Не заменяйте среднее состояний явным полушагом Эйлера. Не вычисляйте координату и скорость последовательно с разными средними. Совпадение с трапецией на линейном тесте не распространяется на нелинейный случай. Сохранение энергии осциллятора не означает точной фазы.
Практика
Задачи с решением
Затухание
Условие. y′=−y, y₀=1, h=0,5. Найдите новый уровень.
Решение. y₁=1−0,25(1+y₁); 1,25y₁=0,75.
Ответ. y₁=0,6.
Нелинейная ветвь
Условие. y′=y², y₀=1, h=0,2. Выберите корень, связанный с h→0.
Решение. y₁=1+0,05(1+y₁)²; y₁²−18y₁+21=0. Корни 9±2√15; выбираем меньший.
Ответ. 9−2√15≈1,254033308.
Дополнительные источники
- https://www.unige.ch/~hairer/preprints/dundee.html
- https://www.unige.ch/~hairer/preprints/iavernaro-mazzia.html
Связанные формулы
Математика
Явный метод средней точки для задачи Коши
Явный метод средней точки делает пробный полушаг по Эйлеру, оценивает производную в середине и использует её на полном шаге. Это двухстадийная схема Рунге — Кутты второго порядка.
Математика
Неявная трапециевидная схема Кранка — Николсон для ОДУ
Неявная трапециевидная схема усредняет производные на концах шага, включая неизвестный новый уровень. Она имеет второй порядок и A-устойчивость, но при крупном шаге затухание может сопровождаться сменой знака.
Физика
Метод Стёрмера — Верле для гармонического осциллятора
Схема Стёрмера — Верле обновляет координату осциллятора по двум соседним моментам и хорошо сохраняет структуру обратимого движения при малом шаге.