между точками.
5.9.5.Сплайновая аппроксимация при большом числе узлов
При большом числе узлов (десятки-сотни и выше) данные представленные точками выглядят нередко не представительно. Например, на рис. 5.21 показан документ, иллюстрирующий сплайновую аппроксимацию функции синуса, представленной 31 отсчетом, но без вывода графика сплайновой функции. Несмотря на равномерное расположение узлов по графику точек невозможно определить, что это функция синуса.
Рис. 5.21 Пример представления функции синуса 31 узловыми точками при равномерном расположении узлов
Рисунок 5.22 отличается от рис. 5.21 только построением сплайновой функции, представленной графическим объектом g1 (на рис. 5.19 он исключен из параметров функции display). После построения графика сплайновой аппроксимирующей функции становится вполне ясным, что точки представляют функцию синуса, которая прекрасно представляется отрезками полиномов сплайн-функции.
Рис. 5.22. Пример сплайновой аппроксимации синусоидальной функции
Здесь полезно обратить внимание на то, что за пределами области узловых точек значения, возвращаемые сплайновой функцией в пакете CurveFitting равны нулю. Так что экстраполяция по ней невозможна (в тоже время функция spline такой возможностью обладает).
5.9.6. Функция реализации метода наименьших квадратов LeastSquares
До сих пор мы рассматривали методы числовой аппроксимации функций или данных, при которых порядок полиномов определялся числом отсчетов функции.
Функция LeastSquares служит для реализации аппроксимации по методу наименьших квадратов. При этом методе происходит статистическая обработка данных (самих по себе или представляющих функцию) исходя из минимума
LeastSquares(xydata, v, opts)
LeastSquares(xdata, ydata, v, opts)
Все входящие в нее параметры были определены выше (см. параметры функции BSplineCurve). Параметр opts задается в форме выражений weight=wlist, curve=f или params=pset.
Следующие примеры иллюстрируют применение функции LeastSquares:
> with(CurveFitting):
LeastSquares([[0,.5], [1,2], [2,4], [3,8]], v);
> LeastSquares([0,1,2,3], [1,2,4,6], v, weight-[1,1,1,10]);
> LeastSquares([0,1,3,5,6], [1,-1,-3,0,5], v, curve=a*v^2+k*v+c);
Наглядную иллюстрацию приближения группы точек кривой (в данном случае представленной полиномом четвертой степени) дает рис. 5.23. Кривая в облаке точек располагается таким образом, что площади квадратов над кривой и под ней в сумме равны нулю.
Рис. 5.23. Графическое представление метода наименьших квадратов
В конце этой главы мы вернемся к реализации метода наименьших квадратов при выполнении регрессионного анализа, построенного па этом методе.
5.9.7. Функция полиномиальной аппроксимации
Функция Polynomial Interpolation реализует полиномиальную интерполяцию и может использоваться в виде:
PolynomialInterpolation(xydata, v)
PolynomialInterpolation(xdata, ydata, v)
Параметры функции были определены выше. Параметр v может быть как именем, так и численным значением. Примеры применения функции представлены ниже
> with(CurveFitting):
PolynomialInterpolation([[0,0], [1,2], [2,4], [3, 3]], z);
> PolynomialInterpolation([0, 2, 5, 8], [2, a, 1, 3], 3);
5.9.8. Функция рациональной аппроксимации
Функция рациональной интерполяции задается в виде:
RationalInterpolation(xydata, z, opts)
RationalInterpolation(xdata, ydata, z, opts)
где необязательный параметр opts задается выражениями method=methodtype или degrees=[d1,d2]. Функция возвращает результат в виде отношения двух полиномов.
Параметр methodtype может иметь значения lookaround или subresultant, задающие учет или пропуск сингулярных точек.
Пример применения функции RationalInterpolation (загрузка пакета опущена, но предполагается):
> xpoints := [0,1,2,3,4,-1]: ypoints := [0, 3, 1, 3, а, 1/11]:
f := RationalInterpolation(xpoints, ypoints, x);
> for i from 1 to 6 do normal(eval(f,x=xpoints[i])-ypoints[i]) end do;
5.9.9. Функция вычисления обычных сплайнов Spline
Функция
Spline(xydata, v, opts)
Spline(xdata, ydata, v, opts)
