Построение и анализ сплайн-интерполяции: кубические сплайны vs B-сплайны на неравномерной сетке

Тема относится к вычислительной математике и численному анализу. Главная сложность здесь связана с неравномерностью шага сетки, при которой формулы классической трехточечной прогонки усложняются, а выбор узлового вектора для B-сплайнов начинает напрямую влиять на обусловленность систем линейных уравнений.

Постановка задачи и теоретические основы сплайнов

В теоретической части необходимо строго определить оба математических объекта на отрезке с сеткой x_0 < x_1 < ... < x_n. Для кубического сплайна выпишите требования непрерывности самой функции, а также ее первой и второй производных во внутренних узлах. Обязательно разберите постановку краевых условий, поскольку именно они замыкают систему уравнений. Рассмотрите естественные краевые условия с нулевой второй производной на концах отрезка, зажатые края с заданными первыми производными и периодические условия.

Для аппарата B-сплайнов приведите рекуррентную формулу Кокса — де Бора. Четко разграничьте точки интерполяции данных и элементы расширенного узлового вектора t_j, кратность которых на границах обычно принимают равной порядку сплайна k = 4 для достижения совпадения крайних точек кривой с граничными узлами интерполяции.

Алгоритмическая реализация и вычислительная сложность

Практическая часть требует программной реализации обоих методов без использования готовых высокоуровневых функций подгонки кривых. В алгоритме кубических сплайнов вывод сводится к решению системы линейных уравнений с трехдиагональной матрицей относительно моментов сплайна, то есть вторых производных. При неравномерном шаге сетки h_i = x_{i+1} - x_i диагональное преобладание сохраняется, что гарантирует корректность применения метода прогонки с вычислительной сложностью порядка O(n).

Для B-сплайнов построение интерполянта требует формирования ленточной матрицы интерполяции размера n на n с полушириной ленты, зависящей от степени базиса. Покажите, как алгоритм де Бора вычисляет значения базисных функций за время O(k^2) для заданной точки, и сопоставьте число арифметических операций при заполнении и решении систем для обеих моделей.

Свойства локальности и устойчивость к шуму

Центральный эксперимент работы должен продемонстрировать различие в локальности сплайнов. При смещении одного входного значения классический кубический сплайн меняет форму на всем расчетном интервале, хотя амплитуда возмущения экспоненциально затухает по мере удаления от измененной точки. B-сплайн третьей степени по построению финитен, поэтому локальное возмущение точки строго ограничено четырьмя соседними интервалами сетки.

Для проверки устойчивости к экспериментальным погрешностям сгенерируйте аналитическую тестовую функцию, например функцию Рунге f(x) = 1 / (1 + 25x^2) или синусоиду с затуханием, и наложите аддитивный гауссов шум с нулевым математическим ожиданием. Сравните среднеквадратическую погрешность аппроксимации на плотной сетке проверки при варьировании дисперсии шума от 1% до 10% от амплитуды сигнала.

Что преподаватель оценивает при проверке

• Математическую строгость вывода сеточных уравнений с учетом переменного шага h_i
• Корректность выбора узлового вектора для B-сплайнов по алгоритмам де Бора или метода усреднения узлов
• Собственную программную реализацию алгоритмов решения разреженных СЛАУ, а не вызов библиотечной интерполяции одной строкой
• Наглядность графического материала с визуализацией базисных функций, графиков интерполянтов и карт погрешностей
• Анализ поведения графиков вблизи границ интервала при возникновении краевых осцилляций