Cspline (vx, vy)

Pspline (vx, vy)

Lspline (vx, vy)

Они возвращают вектор коэффициентов вторых производных, который мы будем называть vs. Этот вектор vs обычно используется в функции interp, описанной ниже. Аргументы vx и vy должны быть вещественными векторами одинаковой длины. Значения вектора vx должны быть расположены в порядке возрастания.

Эти три функции отличаются только граничными условиями:

· функция lspline генерирует кривую сплайна, которая приближается к прямой линии в граничных точках;

Рис. 8.3. Кубическая сплайн-интерполяция

· функция pspline генерирует кривую сплайна, которая приближается к параболе в граничных точках.

· функция cspline генерирует кривую сплайна, которая может быть кубическим полиномом в граничных точках.

Чтобы провести кубический сплайн через набор точек (рис. 8.3):

· создайте векторы vx и vy, содержащие координаты x и y, через которые нужно провести кубический сплайн. Элементы vx должны быть расположены в порядке возрастания. (Хотя мы используем имена vx, vy и vz, нет никаких ограничений на имена этих переменных; можно использовать любые имена);

· Вычислите вектор vs: =cspline (vx, vy).

Рис. 8.4. Экстраполяция

Вектор vs содержит вторые производные интерполяционной кривой в рассматриваемых точках;

· чтобы найти интерполируемое значение в произвольной точке x₀, вычислите interp(vs, vx, vy, x₀), где vs, vx, vy - векторы, описанные ранее.

Интерполируемое значение для конкретного x: есть ордината y соответствующей точки сплайна. Для значений x, расположенных перед первой точкой вектора vx, MathCad продолжает сплайн первой из составляющих его кубических парабол. Для значений x, расположенных за последней точкой vx, MathCad продолжает сплайн первой из составляющих его кубических парабол.

Для получения наилучших результатов x должно находится между самыми большими и самыми маленькими значениями vx - маловероятно, что будут полезны значения, вычисленные для х вне этого диапазона. Сплайны предназначены для интерполяции, а не для экстраполяции.

Линейное предсказание (экстраполяция)

Иногда необходимо оценить значения формул в точках, находящихся вне области расположения сетки, на которой заданы значения функции. В MathCad есть функция predict, которая позволяет это сделать.

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

Функция predict использует последние m исходных значений данных, чтобы вычислять коэффициенты предсказания. Как только это сделано она использует последние m точек, чтобы предсказать координаты (m+1)-ой точки, фактически создавая скользящее окно шириной в m точек (рис. 8.4).