Численные методы II
Численные методы II: итерационные решатели, схемы интегрирования ОДУ, устойчивость и контроль погрешности на практике.
1. Карта курса
- Итерационные методы для больших СЛАУ.
- Семейство Рунге-Кутты и многошаговые схемы.
- Жесткие задачи и неявные методы.
- Апостериорный контроль точности.
Рекомендуемый режим освоения: теория \(\rightarrow\) 2-3 задачи вручную
\(\rightarrow\) вычислительная проверка \(\rightarrow\) короткий конспект ошибок.
1.1. Лекционный маршрут и выходные компетенции
- Недели 1-3: постановка вычислительной задачи и типы ошибок.
- Недели 4-6: базовые алгоритмы и теоретическая оценка порядка сходимости.
- Недели 7-10: устойчивость, обусловленность и итерационные критерии остановки.
- Недели 11-14: адаптивные схемы, валидация на тестовых задачах, инженерный отчёт.
2. Теоретический каркас
Численно-вычислительные дисциплины строятся вокруг трех свойств метода: аппроксимация, устойчивость, сходимость. Полноценное решение всегда включает анализ ошибки и вычислительную сложность.
Инженерный критерий: метод хорош только если точность,
устойчивость и время вычисления согласованы с ограничениями задачи.
Практический стандарт: перед запуском крупного эксперимента
проводится grid refinement test и сравнение с эталонным (или аналитическим) решением.
2.1. Строгий минимум раздела
- Для каждого численного метода фиксировать три свойства: аппроксимация, устойчивость, сходимость.
- Разделять погрешность дискретизации, итерационную погрешность и ошибку машинной арифметики.
- Указывать критерий остановки в терминах нормы невязки, а не числа итераций.
- Проводить grid-refinement test и подтверждать порядок сходимости экспериментально.
Требование уровня "отлично": уметь не только выполнить вычисление,
но и письменно обосновать корректность каждого шага с явным указанием условий применимости.
3. Ключевые формулы и зависимости
| Концепт | Формула / интерпретация |
|---|---|
| Глобальная ошибка | \(\|x-x_h\|=O(h^p)\) |
| Итерационный критерий | \(\|x^{k+1}-x^k\|<\varepsilon\) |
| Число обусловленности | \(\kappa(A)=\|A\|\|A^{-1}\|\) |
| Апостериорная оценка | \(\eta_h\approx\|R_h\|\) |
| Итерация Якоби/Зейделя | \(x^{k+1}=Bx^k+c\) |
| Критерий сходимости | \(\rho(B)<1\) |
| Явный Эйлер | \(y_{n+1}=y_n+h f(t_n,y_n)\) |
4. Методика решения типовых задач
- Нормировать задачу и оценить масштаб величин до дискретизации.
- Выбирать сетку/шаг по требуемой точности, а не по удобству кода.
- Сопровождать каждый численный результат оценкой погрешности.
- Разделять ошибки модели, метода и машинной арифметики.
- Сравнивать прямые и итерационные методы с учетом ресурсов.
- Проводить тесты устойчивости к шуму входных данных.
4.1. Формат эталонного решения
- Нормировать входные данные и определить масштаб переменных.
- Выбрать дискретизацию, обосновать шаг и ожидаемый порядок точности.
- Контролировать невязку и апостериорную ошибку на каждом этапе.
- Провести исследование сходимости при сгущении сетки.
- Сравнить альтернативные алгоритмы по точности и вычислительным затратам.
Оформление полного решения: постановка \rightarrow выбор метода \rightarrow
вычисления \rightarrow контроль ошибки \rightarrow интерпретация результата.
5. Разбор прикладного кейса
Кейс: жесткая система химической кинетики
Сравниваются явная и неявная схемы при одинаковом шаге. Явная схема теряет устойчивость, тогда как неявная даёт устойчивое решение при существенно большем шаге и меньшей суммарной стоимости расчета.
Проверка результата: после вычислений обязательно фиксировать
1) численную стабильность, 2) чувствительность к параметрам,
3) интерпретацию в терминах исходной предметной задачи.
5.1. Углубление кейса
Качественный численный кейс должен завершаться верификацией: совпадает ли наблюдаемый порядок с теоретическим и где начинается деградация точности из-за округления. Без этого численный ответ нельзя считать завершенным результатом.
- Построить таблицу ошибок при нескольких шагах \\(h\\).
- Оценить обусловленность и объяснить чувствительность к шуму входных данных.
- Сравнить минимум два метода по точности на единицу вычислительных затрат.
6. Типичные ошибки
- Слишком крупный шаг дискретизации без оценки влияния на ответ.
- Подмена устойчивости сходимостью.
- Остановка итераций по количеству шагов вместо нормы невязки.
- Игнорирование обусловленности матрицы.
- Отсутствие валидации на контрольных задачах.
6.1. Диагностика ошибок
- Разделять нестабильность метода и плохую обусловленность исходной задачи.
- Не завершать итерации только по лимиту шагов без нормы невязки.
- Проверять влияние округления при очень малых шагах.
- Фиксировать, как шум входа искажает итоговую оценку.
7. Практикум (3 уровня)
Уровень A: базовая техника
- Решить 12-15 стандартных задач с полной записью решения.
- Для каждой задачи указать примененный метод и почему он корректен.
- Проверить 3 задачи альтернативным методом.
Уровень B: продвинутая отработка
- Решить 8 задач с параметрами и анализом вырожденных случаев.
- Оценить погрешность/устойчивость результата на вариациях входа.
- Подготовить короткий отчёт с выводами (1-2 страницы).
Уровень C: мини-проект
- Реализовать вычислительный прототип по теме курса.
- Сравнить минимум 2 метода и обосновать выбор лучшего.
- Подготовить репродуцируемый notebook с графиками и выводами.
8. Экзаменационный минимум и литература
Минимум к экзамену
- Все базовые определения курса в строгой формулировке.
- Ключевые теоремы/критерии и условия их применимости.
- Алгоритм решения типовой задачи каждого раздела.
- Умение объяснить источник ошибки и устойчивость метода.
- Интерпретация результата в прикладном контексте.
Рекомендуемая литература
- Hairer, Wanner. Solving ODEs; Самарский А.А. Теория разностных схем
- Материалы семинаров и практикумов кафедры.
- Набор задач повышенной сложности (подготовка к экзамену).
Тренажер билетов
- Дайте строгое определение устойчивости метода и приведите пример неустойчивой схемы.
- Выведите оценку погрешности выбранного алгоритма и укажите, что ограничивает порядок.
- Покажите, как выбирать шаг/сетку при фиксированном бюджете вычислений.
- Разберите ситуацию плохой обусловленности: как диагностировать и как стабилизировать решение.
План повторения перед экзаменом
Эффективная подготовка строится циклом "теория \rightarrow задачи \rightarrow разбор ошибок". Для каждого раздела фиксируйте: формулировку ключевых определений, один эталонный алгоритм решения и типовую ловушку, которая чаще всего приводит к неверному ответу.
- Сделать короткий one-page summary по каждому разделу с формулами и условиями применимости.
- Решить минимум 2 задачи базового и 1 задачу повышенного уровня по каждому крупному блоку.
- Провести устный прогон билета: формулировка теоремы, схема доказательства, прикладная интерпретация.