Главная страница

Методы интерполяции и аппроксимации


Скачать 93.46 Kb.
НазваниеМетоды интерполяции и аппроксимации
Дата31.01.2022
Размер93.46 Kb.
Формат файлаdocx
Имя файлаInterp_app.docx
ТипДокументы
#347959
страница11 из 13
1   ...   5   6   7   8   9   10   11   12   13

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


Суть метода наименьших квадратов заключается в нахождении таких значений хi, при которых сумма квадратов отклонений (ошибок) ei=yifi(x) будет стремиться к минимуму


2
n


e
i

n

( yi
f(xi
))2 min . (9)

i1

i1 x

Т.к. каждое значение xiв общем случае «сопровождается» соответствующим коэффициентом аi(i = 0, 1, 2, …, n), то задача сводится к нахождению данных коэффициентов. Введем обозначение функции

n

F(a, a, ...,a) ( y f(x))2.

(10)

0 1 n i ii1

Тогда, на основе обращения в точке минимума функции F в нуль ее производных, для определения вышеупомянутых коэффициентов составляется нормальная система:

dF



da0

dF

da

0;
0;

1

...



dF

dan

 0.

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


В ряде экспериментов данные распределяются таким образом, что оказывается возможным описать их изменение линейной зависимостью (линейным уравнением) (рис. 7)

P(x)=ax+b. (11)

Формулы для расчета коэффициентов aи bопределяются по методу наименьших квадратов (9), подставив (11) в (10)

n

F ( yi

i1

  • a xi

b)2 min . (12)



Рис. 7. Линейная аппроксимация

Для решения (12) составляется система из двух уравнений с двумя неизвестными

dF

da

dF

0;

(13)

0.

db

Подставляя в (13) формулу (12), получаем

dF

n

2 ( y

  • a x

b) 1 0,

db



i i

i1
(14)

dF

da

2

n

( yi

i1

и

  • a xi

  • b) xi

0.

an x2 bn x

n

 (x
y)

i

i1

n

i

i1

n

i

i1

i

, (15)



a xi

  • nb yi

i1

i1

Решая полученную систему (15) методом подстановки, получаем формулы для нахождения коэффициентов aи b:

n n n

n (xi yi) xi yi

ai1 i1 i1 , (16)

n n 2

n x2 x

i

i1

i

i1 

n n

y a x

n n

yx2

n n

x (x

y)

i i

i i

i i i

b i1 i1 i1 i1 i1 i1 . (17)

n n n 2

n x2 x



1   ...   5   6   7   8   9   10   11   12   13


написать администратору сайта