← К программе курса 3 семестр

Численные методы I

Первый блок численных методов: корни нелинейных уравнений, интерполяция, квадратурные формулы и вычислительная устойчивость.

Фокус
Строгая теория + вычислительная реализация
Формат
Лекции, практикум, контрольные задачи
Цель
Уровень: самостоятельное решение сложных профильных задач

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}=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. Методика решения типовых задач

  1. Нормировать задачу и оценить масштаб величин до дискретизации.
  2. Выбирать сетку/шаг по требуемой точности, а не по удобству кода.
  3. Сопровождать каждый численный результат оценкой погрешности.
  4. Разделять ошибки модели, метода и машинной арифметики.
  5. Сравнивать прямые и итерационные методы с учетом ресурсов.
  6. Проводить тесты устойчивости к шуму входных данных.

4.1. Формат эталонного решения

  1. Нормировать входные данные и определить масштаб переменных.
  2. Выбрать дискретизацию, обосновать шаг и ожидаемый порядок точности.
  3. Контролировать невязку и апостериорную ошибку на каждом этапе.
  4. Провести исследование сходимости при сгущении сетки.
  5. Сравнить альтернативные алгоритмы по точности и вычислительным затратам.
Оформление полного решения: постановка \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. Экзаменационный минимум и литература

Минимум к экзамену

  • Все базовые определения курса в строгой формулировке.
  • Ключевые теоремы/критерии и условия их применимости.
  • Алгоритм решения типовой задачи каждого раздела.
  • Умение объяснить источник ошибки и устойчивость метода.
  • Интерпретация результата в прикладном контексте.

Рекомендуемая литература

  • Самарский А.А., Гулин А.В. Численные методы; Калиткин Н.Н. Численные методы
  • Материалы семинаров и практикумов кафедры.
  • Набор задач повышенной сложности (подготовка к экзамену).

Тренажер билетов

  1. Дайте строгое определение устойчивости метода и приведите пример неустойчивой схемы.
  2. Выведите оценку погрешности выбранного алгоритма и укажите, что ограничивает порядок.
  3. Покажите, как выбирать шаг/сетку при фиксированном бюджете вычислений.
  4. Разберите ситуацию плохой обусловленности: как диагностировать и как стабилизировать решение.

План повторения перед экзаменом

Эффективная подготовка строится циклом "теория \rightarrow задачи \rightarrow разбор ошибок". Для каждого раздела фиксируйте: формулировку ключевых определений, один эталонный алгоритм решения и типовую ловушку, которая чаще всего приводит к неверному ответу.

  • Сделать короткий one-page summary по каждому разделу с формулами и условиями применимости.
  • Решить минимум 2 задачи базового и 1 задачу повышенного уровня по каждому крупному блоку.
  • Провести устный прогон билета: формулировка теоремы, схема доказательства, прикладная интерпретация.