Методы решения. Данная практическая работа выполняется по вариантам

Данная практическая работа выполняется по вариантам. Первый вариант – это метод Ньютона, второй – Лагранжа. Эти полиномы являются степенными.

Известно, что через две точки можно провести одну и только одну прямую, через три – одну и только одну параболу и т.д. Поэтому, через n+1 точку {x_i} можно провести одну и только одну кривую порядка n. Отсюда можно сделать два вывода. Во-первых, чем больше количество точек в заданной сетке, тем выше, в общем случае, будет степень полинома P(x). Именно этим и объясняется осциллирующее поведение полиномов Ньютона и Лагранжа при большом количестве точек – просто их вид становится слишком сложным. Отметим, что для других интерполирующих полиномов это может быть и не так. Например, МНК, независимо от количества точек, дает полином, для которого выполняется условие (2.5.3). Т.е., если в качестве линейно-не­за­ви­си­мых функций взять φ _i(x) = xⁱ, i = 1, 2, …, m, то можно построить, например, кубический полином для любого количества точек (при m = 3). Порядок у него ниже, поэтому он более гладкий. При m = n МНК становится обычным интерполяционным полиномом. Во-вторых, полиномы Ньютона и Лагранжа совпадают, т.е. это просто две формы записи одного и того же полинома, и их можно преобразовать к следующему виду:

где k_i – некоторые константы. Индекс n у полинома указывает на его порядок.

При этом, каждым вариантом необходимо реализовать 6 задач:

1. Вычисление полинома на равномерной сетке;

2. Вычисление полинома на неравномерной сетке;

3. Вычисление первой производной полинома на равномерной сетке;

4. Вычисление первой производной полинома на неравномерной сетке;

5. Вычисление второй производной полинома на равномерной сетке;

6. Вычисление второй производной полинома на неравномерной сетке;

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

(2.5.9)

и подставляется в полином и его производные. Таким образом, получается, что они зависят только от q, а x и {x_i} явным образом в них не входят. Т.е. имеем P(q). Получить его можно самостоятельно, сделав замену (2.5.9) в полиноме P(x).

Выигрыш состоит в том, что не нужно хранить в памяти узлы сетки {x_i}, поэтому ее используется примерно в два раза меньше.