методом наименьших квадратов

Исследование полиномиальной аппроксимации

Составитель: Панкратов Ю.M.

Санкт-Петербург

Лабораторная работа №7

Целью работы является умение интерполировать дискретные данные полиномом и исследование погрешностей интерполяции в зависимости от степени полинома.

1. Теоретические сведения

Метод наименьших квадратов применяют, когда необходимо вывести формулу аппроксимирующей кривой, описывающей некоторую зависимость, полученную в результате, например, эксперимента. Поскольку экспериментальные данные получают, как правило, с некоторой погрешностью, то нет смысла использовать, например, интерполяционный многочлен Лагранжа, всегда проходящий через узлы интерполяции. В этом случае обычно проводят аппроксимирующую кривую, которая не проходит через экспериментальные точки, но в то же время учитывает исследуемую закономерность, сглаживая случайные выбросы результатов эксперимента. В качестве аппроксимирующей функции может использоваться линейная комбинация F(x) любых функций

F(x)= a₀f₀(x)+a₁f₁(x)+ … + a_nf_n(x) , (1)

где f₀(x), f₁(x),…, f_n(x) - набор любых функций, называемых базисными; a₀ ,a₁,…,a_n - набор коэффициентов.

Довольно часто в качестве базисных функций используют комбинацию из возрастающей степенной последовательности аргумента f₀(x)=1, f₁(x)=x, f₂(x)=x²,…, f_n(x)=xⁿ, при этом аппроксимирующая функция будет являться алгебраическим полиномом. Именно этот случай и рассматривается в этой работе при исследовании метода наименьших квадратов.

Пусть некоторая функция задана таблицей значений y_i в узлах x_i (i=0,1,2,…,m), где (m+1) – количество узлов (см. рисунок). Требуется определить коэффициенты a₀ ,a₁,…,a_n полинома

P (x)=a₀+a₁x+ … +a_nxⁿ (2)

таким образом, чтобы сумма S квадратов отклонений полинома P (x) от значений y_i аппроксимируемой функции в заданных точках была бы минимальной, т.е.

(3)

Если степень полинома n = m, то полином (2) будет интерполяционным и пройдет через все узлы и отклонения будут нулевыми. Стремясь понизить степень полинома, принимают n < m, при этом в общем случае аппроксимирующая кривая не будет проходить через заданные точки (x_i, y_i).

Сумма квадратов отклонений зависит от коэффициентов полинома a₀ ,a₁,…,a_n. Поскольку функция S принимает только положительные значения и имеет минимум, вычислим частные производные по всем переменным a₀ ,a₁,…,a_n и приравняем их нулю:

Выполнив суммирование, собрав коэффициенты при каждом a₀ ,a₁,…,a_n и перенеся не содержащие их суммы в правую часть, получим СЛАУ относительно неизвестных a₀ ,a₁,…,a_n:

В матричной форме эту систему можно записать как a×S = b, где

Решение этой системы позволяет найти все коэффициенты a₀ ,a₁,…,a_n и значение полинома в любой точке x вычисляют по формуле (2)

Для полиномиальной аппроксимации (регрессии) в MathCAD’е можно воспользоваться следующими встроенными функциями:

1. regress(x,y,k) – возвращает вектор коэффициентов a полинома, при этом всегда первые три компоненты вектора есть вектор вторых производных и являются параметрами для описываемой ниже функции interp, а остальные компоненты и есть вектор коэффициентов a, где

x - вектор данных аргумента, элементы которого должны быть расположены в порядке возрастания;

y - вектор значений того же размера;

k - степень полинома (целое положительное число).

Функция regress позволяет, вычислив коэффициенты a₀ ,a₁,…,a_n, построить полином в форме (2) и уже его использовать, например, для построения графика.

Например, при значении степени полинома k=2 функция regress возвратит представленное справа значение.

2. interp (s,x,y,t) – возвращает результат полиномиальной регрессии, где

s - вектор вторых производных, созданный, например, предыдущей функцией;

x и y - то же, что и в предыдущей функции;

t - текущее значение аргумента полинома.

Функция interp позволяет непосредственно вычислить значение полинома при любом значении аргумента без построения полинома как такового в форме (2).

3. linfit(x,y,F) – более универсальная функция, возвращающая вектор коэффициентов a линейной комбинации функций (1), где

x и y - то же, что и в предыдущей функции;

F - вектор базисных функций f₀(x), f₁(x),…, f_n(x) в выражении (1). Например, для алгебраического полинома степени n =2 необходимо записать F в следующем виде:

Функция linfit также позволяет непосредственно вычислить значение полинома при любом значении аргумента без построения полинома как такового в форме (2). Например, при значении аргумента, равного 2.3, получим:

2. Задание: апроксимировать табличные значения методом наименьших квадратов полиномами 1, 2, 3 и 4 степеней, при этом для полиномов:

§ 1, 2 и 3 степеней вывести формулы полиномов в форме (2);

§ 1 и 2 степени - использовать функцию regress;

§ 3-й степени - использовать функцию linfit;

§ 4-й степени - использовать функцию interp.

Построить совмещенный график узлов и всех полиномов.

Построить график суммы абсолютных отклонений полиномов во всех узлах в зависимости от степени полинома. Сделать выводы по результатам исследования.

Варианты задания

№	xi	yi	№	xi	yi
1.											16.
2.											17.
3.											18.
4.											19.
5.											20.
6.											21.
7.											22.
8.											23.
9.											24.
10.											25.
11.											26.
12.											27.
13.											28.
14.											29.
15.											30.

3. Пример выполнения задания

4. Литература

1. Бахвалов Н.С. Численные методы. – М., Наука,1975. –632 c.

2. Крылов В.И., Бобков В.В., Монастырный П.И. Вычислительные методы. – М., Наука,1976. –304 с.

3. Волков Е.А. Численные методы: Учебное пособие для вузов. – 2-е изд. – М., Наука, 1987.

4. Копченова Н.В., Марон И.А. Вычислительная математика в примерах и задачах. – М., Наука,1972.

5. Турчак Л.И., Плотников П.В. Основы численных методов: Учеб. пособие.– 2-ое изд., перераб. и доп. -М.: ФИЗМАТЛИТ, 2002. –304 c.

6. Исаков В.Н. Элементы численных методов. – М., Издательский центр «Академия»,2003. –192 c.

---

Дата добавления: 2015-10-24