Численные методы I
Первый блок численных методов: корни нелинейных уравнений, интерполяция, квадратурные формулы и вычислительная устойчивость.
1. Карта курса
- Бисекция, Ньютона, секущих и условия их сходимости.
- Полиномиальная интерполяция и оценка остатка.
- Численное интегрирование и порядок квадратур.
- Погрешности округления и влияние обусловленности.
1.1. Лекционный маршрут и выходные компетенции
- Недели 1-3: постановка вычислительной задачи и типы ошибок.
- Недели 4-6: базовые алгоритмы и теоретическая оценка порядка сходимости.
- Недели 7-10: устойчивость, обусловленность и итерационные критерии остановки.
- Недели 11-14: адаптивные схемы, валидация на тестовых задачах, инженерный отчёт.
2. Теоретический каркас
Численно-вычислительные дисциплины строятся вокруг трех свойств метода: аппроксимация, устойчивость, сходимость. Полноценное решение всегда включает анализ ошибки и вычислительную сложность.
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}=x_k-\frac{f(x_k)}{f'(x_k)}\) |
| Полином Лагранжа | \(P_n(x)=\sum_{i=0}^n y_i l_i(x)\) |
| Квадратура трапеций | \(\int_a^b f\approx \frac{h}{2}(f_0+2\sum_{i=1}^{n-1} f_i+f_n)\) |
4. Методика решения типовых задач
- Нормировать задачу и оценить масштаб величин до дискретизации.
- Выбирать сетку/шаг по требуемой точности, а не по удобству кода.
- Сопровождать каждый численный результат оценкой погрешности.
- Разделять ошибки модели, метода и машинной арифметики.
- Сравнивать прямые и итерационные методы с учетом ресурсов.
- Проводить тесты устойчивости к шуму входных данных.
4.1. Формат эталонного решения
- Нормировать входные данные и определить масштаб переменных.
- Выбрать дискретизацию, обосновать шаг и ожидаемый порядок точности.
- Контролировать невязку и апостериорную ошибку на каждом этапе.
- Провести исследование сходимости при сгущении сетки.
- Сравнить альтернативные алгоритмы по точности и вычислительным затратам.
5. Разбор прикладного кейса
Кейс: выбор метода поиска корня
Для функции с плохой производной около корня сравниваются бисекция и Ньютон. На практике строится гибрид: сначала бисекция для попадания в надежный интервал, затем Ньютон для квадратичного ускорения.
5.1. Углубление кейса
Качественный численный кейс должен завершаться верификацией: совпадает ли наблюдаемый порядок с теоретическим и где начинается деградация точности из-за округления. Без этого численный ответ нельзя считать завершенным результатом.
- Построить таблицу ошибок при нескольких шагах \\(h\\).
- Оценить обусловленность и объяснить чувствительность к шуму входных данных.
- Сравнить минимум два метода по точности на единицу вычислительных затрат.
6. Типичные ошибки
- Слишком крупный шаг дискретизации без оценки влияния на ответ.
- Подмена устойчивости сходимостью.
- Остановка итераций по количеству шагов вместо нормы невязки.
- Игнорирование обусловленности матрицы.
- Отсутствие валидации на контрольных задачах.
6.1. Диагностика ошибок
- Разделять нестабильность метода и плохую обусловленность исходной задачи.
- Не завершать итерации только по лимиту шагов без нормы невязки.
- Проверять влияние округления при очень малых шагах.
- Фиксировать, как шум входа искажает итоговую оценку.
7. Практикум (3 уровня)
Уровень A: базовая техника
- Решить 12-15 стандартных задач с полной записью решения.
- Для каждой задачи указать примененный метод и почему он корректен.
- Проверить 3 задачи альтернативным методом.
Уровень B: продвинутая отработка
- Решить 8 задач с параметрами и анализом вырожденных случаев.
- Оценить погрешность/устойчивость результата на вариациях входа.
- Подготовить короткий отчёт с выводами (1-2 страницы).
Уровень C: мини-проект
- Реализовать вычислительный прототип по теме курса.
- Сравнить минимум 2 метода и обосновать выбор лучшего.
- Подготовить репродуцируемый notebook с графиками и выводами.
8. Экзаменационный минимум и литература
Минимум к экзамену
- Все базовые определения курса в строгой формулировке.
- Ключевые теоремы/критерии и условия их применимости.
- Алгоритм решения типовой задачи каждого раздела.
- Умение объяснить источник ошибки и устойчивость метода.
- Интерпретация результата в прикладном контексте.
Рекомендуемая литература
- Самарский А.А., Гулин А.В. Численные методы; Калиткин Н.Н. Численные методы
- Материалы семинаров и практикумов кафедры.
- Набор задач повышенной сложности (подготовка к экзамену).
Тренажер билетов
- Дайте строгое определение устойчивости метода и приведите пример неустойчивой схемы.
- Выведите оценку погрешности выбранного алгоритма и укажите, что ограничивает порядок.
- Покажите, как выбирать шаг/сетку при фиксированном бюджете вычислений.
- Разберите ситуацию плохой обусловленности: как диагностировать и как стабилизировать решение.
План повторения перед экзаменом
Эффективная подготовка строится циклом "теория \rightarrow задачи \rightarrow разбор ошибок". Для каждого раздела фиксируйте: формулировку ключевых определений, один эталонный алгоритм решения и типовую ловушку, которая чаще всего приводит к неверному ответу.
- Сделать короткий one-page summary по каждому разделу с формулами и условиями применимости.
- Решить минимум 2 задачи базового и 1 задачу повышенного уровня по каждому крупному блоку.
- Провести устный прогон билета: формулировка теоремы, схема доказательства, прикладная интерпретация.