Главная страница Случайная страница
Разделы сайта
АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника
Аппроксимация функций по методу наименьших квадратов.
|
|
Пусть на основании экспериментов получена зависимость функции y = φ (x), представленная в виде таблицы.
| x | x1 | x2 | … | x n |
| y | y1 | y2 | … | yn |
Задача состоит в аналитическом представлении функции
y = φ (x) (1)
Возможным вариантом решения является интерполирование с помощью многочлена Ньютона или Лагранжа. Однако, этот способ, требующий обязательного совладения табличных и теоретических значений функции не всегда удобен. При большом количестве узлов интерполяции потребует либо отыскания многочлена большой степени, либо другой «извилистой» функции, график которой проходил бы через все указанные точки. Поэтому рассматривается метод аппроксимации», т.е. приближенного представления искомой функции. Вид функции y = φ (x) устанавливается либо из теоретических соображений, либо на основании расположения на координатной плоскости точек, соответствующих экспериментальным данным. Например, если экспериментальные точки расположены так, как на рис.1, то функцию следует искать в виде y =a0 + a1 x, а на рис.2 в виде y =a0 + a1x + a2.
y •
•
•
•
• •
х
x1 x2 x3 x4 …. xn
Рис. 1
у 
• • •
• • • •
• • •
•
х
Рис. 2.
В общем случае, при выбранном виде функции
y = φ (x, a0, a1, a2, …a k )
следует подобрать параметры a0, a1, a2, …k так, чтобы эта функция в каком-то смысле наилучшим образом описывала рассматриваемый процесс. Наиболее распространенным методом решения данной задачи является метод наименьших квадратов. Этот метод заключается в следующем. Рассмотрим сумму квадратов разностей между экспериментальными значениями и теоретическими значениями функции y = φ (x, a, b, c, …) в соответствующих точках.
(2)
(n – количество точек, в которых известны значения функции, k – количество параметров)
Подберем параметры a0, a1, a2, … так, чтобы эта сумма имела наименьшее значение. Таким образом, задача сведется к нахождению значений параметров a0, a1, a2, …, при которых сумма S(a0, a1, a2, …) имеет минимум.
Из теории функций многих переменных известно, что необходимым условием существования экстремума является равенство нулю частных производных.
(3)
Или в развернутом виде
(4)
Получили систему уравнений относительно a0, a1, a2, ….. Здесь столько уравнений, сколько неизвестных. В каждом отдельном случае исследуется вопрос о существовании решения этой системы и о существовании минимума функции y = φ (x, a0, a1, a2, …a k) при найденных значениях a0, a1, a2, ….ak. В математическом анализе доказывается, что если в качестве аппроксимирующей функции y = φ (x, a0, a1, a2, …a k) берется многочлен
y = φ (x, a0, a1, a2, …, a k ) = a0 + a1x + a2 x2 + …+ a k-1 xk -1 + a k xk, (5)
то решение a0, a1, a2, … системы (4) дает минимум функции y = φ (x, a0, a1, a2, …a k ).
Будем искать аппроксимирующую функцию в виде (5).
Тогда система (4) примет вид
(6)
Раскрывая скобки, получим
(7)
………………………………………………………………….
Получаем систему линейных уравнений для определения неизвестных коэффициентов, a0, a1, a2,...., ak.. Из характера задачи следует, что система имеет единственное решение и что при полученных значениях параметров a0, a1, a2, …a k.
Запишем систему (7) в матричном виде AX = B. Здесь А – матрица системы (матрица, составленная из коэффициентов при неизвестных параметрах a0, a1, a2, …a k), Х – столбец неизвестных, В – столбец свободных членов.
B = 
Решение системы (7) в матричном виде запишется X0 = A -1 B.
|
|