Приближение табличных функций по методу наименьших квадратов постановка задачи аппроксимации

• Название: Приближение табличных функций по методу наименьших квадратов постановка задачи аппроксимации
• Дата: 11.10.2022
• Размер: 1.09 Mb.
• Имя файла: priblizhenie_tablichnykh_funkcij_po_mnk.pdf
• Тип: Документы; #727613

ПРИБЛИЖЕНИЕ ТАБЛИЧНЫХ ФУНКЦИЙ
ПО МЕТОДУ НАИМЕНЬШИХ КВАДРАТОВ
Постановка задачи аппроксимации
По результатам экспериментов получена таблица с произвольным расположением аргументов:


,
i
i
x y ,
1,
i
n

. Аналитическое выражение табличной функции может быть неизвестным. На основе этой таблицы требуется найти формулу
 
F
F x

, приближённо описывающую зависимость между экспериментальными данными таблицы. При этом отклонение значений в точках
i
x
,
1,
i
n

, вычисленные по формуле
 
F
F x

, от экспериментальных данных
i
y
должны быть минимальными.
Дано: точки наблюдения


,
i
i
x y
1,
i
n

(
их количество n
).
x
x
1
x
2
…
x
n
y
y
1
y
2
…
y
n
Найти функцию
 
F x и значение
 
*
F x
,


*
1
;
n
x
x x

:
1.
 
i
i
i
F
F x
y


,
1,
i
n
 
2.


,
min


y F
,


1 2
,
,
,
n
y y
y

y
,


1 2
,
,
,
n
F F
F

F
Определение. Полученное функциональное соотношение
 
F x
для приближённого описания зависимости между экспериментальными данными таблицы


,
i
i
x y ,
1,
i
n

называется эмпирической формулой, а функция
 
F x – эмпирической функцией.
Возможным вариантом решения этой задачи является интерполирование. Однако этот способ, требующий обязательного совпадения табличных значений табличной и приближающей функций во всех табличных аргументах мало пригоден. При большом количестве узлов он является неудобным и сложным, ибо потребует отыскания многочлена большой степени.
Кроме того, экспериментальные данные в силу ряда причин могут иметь трудно учитываемые случайные или систематические ошибки.
В этих условиях интерполирование становится сомнительным.
Часто с помощью какой-либо простой функции с проходящим около табличных точек графиком, удаётся добиться эффекта сглаживания ошибок и получить более точное приближение.
Поиск эмпирической формулы
 
F x
начинается с определения класса функций, которые лучше всего отражают связь между табличными данными. Эффективным методом для этого являются графические изображения. На координатной плоскости отмечаются определяемые данной функцией точки, а затем по характеру их расположения подбирается вид приближения из числа известных элементарных функций.

В перечень наиболее часто используемых классов функций входят, например, линейные
1 0
y
a x
a


, полиномиальные
2 2
1 0
y
a x
a x a



логарифмические ln
y
a
x
b


, экспоненциальные
bx
y
ae

, степенные
b
y
ax

, дробно-рациональные
1
y
ax
b


, гиперболические
b
y
a
x
 
и другие функциональные зависимости.
Метод наименьших квадратов для поиска приближённой зависимости
Поиск числовых параметров выбранной эмпирической зависимости
 
F x
сводится к решению задачи поиска минимального значения метрики


 


2 1
,
n
i
i
i
y
F x





y F
y = ax + b
0 10 20 30 40 50 0
5 10 15
y = ax
2
+ bx + c
0 5
10 15 20 25 30 35 0
5 10 15
y = a lnx + b
0 2
4 6
8 10 12 14 0
5 10 15
y = ae
bx
0 10000 20000 30000 40000 0
5 10 15
y = ax
b
0 20000 40000 60000 80000 100000 0
5 10 15 0
1 2
3 4
5 0
5 10 15
y = 1/(ax + b)

Определение.
Квадрат метрики


 


2 2
1
,
n
i
i
i
y
F x





y F
называется среднеквадратическим отклонением.
Задача приближения по методу наименьших квадратов. Требуется определить такие значения параметров
0 1
,
,
n
a a
a
выбранной эмпирической зависимости


0 1
,
, ,
n
F x a a
a
, при которых квадрат метрики


2
,

y F (среднеквадратическое отклонение) будет наименьшим:






2 2
0 1
1
,
,
,
,
min
n
i
i
n
i
y
F x a a
a






y F
Эта задача является задачей поиска минимума функции многих переменных
0 1
,
,
n
a a
a
:






2 0
1 0
1 1
,
,
,
,
,
,
min
n
n
i
i
n
i
S x a a
a
y
F x a a
a





Необходимое условие экстремума функции многих переменных – равенство нулю всех частных производных функции:
1,
i
n
 
:
0
i
S
a



Поскольку










2 0
1 0
1 0
1 1
1
,
,
,
,
,
,
2
,
,
,
0
n
n
i
i
n
i
n
i
i
n
i
i
i
i
i
y
F x a a
a
F x a a
a
S
y
F x a a
a
a
a
a















, параметры
0 1
,
,
n
a a
a
определяются из системы уравнений:






0 1
0 1
1
,
,
,
,
,
,
0
n
i
n
i
i
n
i
i
F x a a
a
y
F x a a
a
a






,
1,
i
n

(1)
Линейная аппроксимация
Пусть на основании графического изображения в качестве эмпирической функции выбирается линейная функция
 
1 0
F x
a x a


Система (1) примет вид:












1 0
1 0
1 0
1 0
1 0
1 1
0 0
n
i
i
i
i
n
i
i
i
i
F a x
a
y
a x
a
a
F a x
a
y
a x
a
a







 








 




,








1 0
1 1
0 1
1 0
0
n
i
i
i
n
i
i
i
i
y
a x
a
y
a x
a
x





  





 




Отсюда получается система линейных уравнений для определения параметров
0 1
,
a a
линейной аппроксимации:
2 1
0 1
1 1
1 0
1 1
n
n
n
i
i
i
i
i
i
i
n
n
i
i
i
i
a
x
a
x
x y
a
x
na
y



















(2)
Задача 1
Методом наименьших квадратов определить параметры линейной эмпирической функции
 
1 0
F x
a x a


по экспериментальным данным, представленным в таблице. Найти
 
6
F
x
y
1 1,6 2
1 3
3 4
3,3 5
5,4
Решение
Количество точек наблюдения n = 5.
Для определения параметров линейной зависимости
 
1 0
F x
a x a


составляется вспомогательная таблица:
x
y
x
2
xy
1 1,6 1 1,6 2
1 4
2 3
3 9
9 4 3,3 16 13 5 5,4 25 27

15 14
55
53
Полученные значения сумм подставляются в систему линейных уравнений для определения параметров
0 1
,
a a
линейной аппроксимации:
1 0
1 0
55 15 52,8 15 5
14, 3
a
a
a
a







Комментарий к решению системы. Точные методы решения этой системы линейных уравнений известны: решение системы с помощью обратной матрицы, с помощью формул
Крамера, метод исключения неизвестных – метод Гаусса. В данном случае проще применить метод Гаусса. Действительно, если умножить второе уравнение системы на (-3) и сложить его с первым, то легко определяется параметр
1
a
:

1 0
1 0
55 15 52,8 45 15 42, 9
a
a
a
a






 

,
1 0
1 55 15 52,8 10 9, 9
a
a
a






,
1 0
1 55 15 52,8 0,99
a
a
a






При подстановке
1 0, 99
a

в любое уравнение исходной системы получается
0 0,11
a
 
Искомая эмпирическая функция имеет вид:
 
0,99 0,11
F x
x


 
6 0,99 6 0,11 5,83
F

 

Реализация решения примера в MS Excel
Услуга Мастера диаграмм Построение линии тренда реализует метод наименьших квадратов для поиска коэффициентов эмпирической функции и построения её графика.
Параметры линейной зависимости определяются с помощью Мастера функций: с помощью функции ЛИНЕЙН из категории Статистические функции.
Квадратическая
(параболическая,
полиномиальная
2-й
степени)
аппроксимация
Пусть на основании графического изображения в качестве эмпирической функции выбирается линейная функция
 
2 2
1 0
F x
a x
a x a



Система (1) примет вид:


















2 2
1 0
2 2
1 0
1 0
2 2
1 0
2 2
1 0
1 1
2 2
1 0
2 2
1 0
1 2
0 0
0
n
i
i
i
i
i
i
n
i
i
i
i
i
i
n
i
i
i
i
i
i
F a x
a x
a
y
a x
a x
a
a
F a x
a x
a
y
a x
a x
a
a
F a x
a x
a
y
a x
a x
a
a
























 


















,












2 2
1 0
1 2
2 1
0 1
2 2
2 1
0 1
1 0
0 0
n
i
i
i
i
n
i
i
i
i
i
n
i
i
i
i
i
y
a x
a x
a
y
a x
a x
a
x
y
a x
a x
a
x







 






  






 




y = 0,99x - 0,11
R² = 0,8368 0
1 2
3 4
5 6
0 1
2 3
4 5
6

Отсюда получается система линейных уравнений для определения параметров
0 1
2
,
,
a a a
параболической аппроксимации:
4 3
2 2
2 1
0 1
1 1
1 3
2 2
1 0
1 1
1 1
2 2
1 0
1 1
1
n
n
n
n
i
i
i
i
i
i
i
i
i
n
n
n
n
i
i
i
i
i
i
i
i
i
n
n
n
i
i
i
i
i
i
a
x
a
x
a
x
x y
a
x
a
x
a
x
x y
a
x
a
x
na
y









































Экспоненциальная аппроксимация
Пусть на основании графического изображения в качестве эмпирической функции выбирается экспоненциальная функция
 
bx
F x
ae

Для применения метода наименьших квадратов экспоненциальная функция линеаризуется: ln ln
bx
F
ae

, ln ln
F
a bx


Метод наименьших квадратов применяется к линейной функции
1 0
ˆ
F
a X
a


, где
0
ln
a
a

,
1
a
b

,
ˆ
ln
F
F

,
X
x

Задача 2
Методом наименьших квадратов определить параметры экспоненциальной эмпирической функции
 
bx
F x
ae

по экспериментальным данным, представленным в таблице. Найти минимальную величину


2
,

y F . Найти
 
6
F
x
y
1 2,1 2
6,2 3
13,3 4
24,6 5
117,5
Решение
Количество точек наблюдения n = 5.
Для применения метода наименьших квадратов экспоненциальная функция линеаризуется: ln ln
bx
F
ae

, ln ln
F
a bx



Метод наименьших квадратов применяется для определения коэффициентов линейной функции
1 0
ˆ
F
a X
a


, аппроксимирующей экспериментальные данные


,
i
i
X Y , где
i
i
X
x

, ln
i
i
Y
y

,
0
ln
a
a

,
1
a
b

,
ˆ
ln
F
F

,
1,
i
n

Для определения параметров линейной зависимости
1 0
ˆ
F
a X
a


составляется вспомогательная таблица:
x
y
X = x
Y=lny
X
2
XY
F
(y – F)
2
1 2,1 1
0,7 1
0,7 2,13 0,00 2
6,2 2
1,8 4
3,6 4,82 1,90 3
13,3 3
2,6 9
7,8 10,90 5,76 4
24,6 4
3,2 16 12,8 24,65 0,00 5
117,5 5
4,8 25 23,8 55,73 3815,82

15
163,7
15
13,1
55
48,8
3823,493
Полученные значения сумм подставляются в систему линейных уравнений для определения параметров
0 1
,
a a
линейной аппроксимации:
1 0
1 0
55 15 48,8 15 5
13,1
a
a
a
a







Решением системы являются значения
1 0, 94
a

и
0 0, 20
a
 
. Линеаризованная эмпирическая функция имеет вид: ˆ
0, 94 0, 20
F
x


Определение параметров экспоненциальной эмпирической зависимости производится из ранее принятых соотношений:
0
ln
a
a

,
1
a
b

Откуда
0 0,20 0,82
a
a
e
e




,
1 0,94
b
a


. Искомая экспоненциальная эмпирическая функция будет иметь вид:
 
0,82 0,94
bx
x
F x
ae
e



Логарифмическая аппроксимация
Пусть на основании графического изображения в качестве эмпирической функции выбирается логарифмическая функция ln
F
a
x b


Метод наименьших квадратов применяется к линейной функции
1 0
ˆ
F
a X
a


, где
0
a
b

,
1
a
a

,
ˆ
F
F

, ln
X
x

Степенная аппроксимация
Пусть на основании графического изображения в качестве эмпирической функции выбирается степенная функция
 
b
F x
ax


Для применения метода наименьших квадратов степенная функция линеаризуется:
 
ln ln
b
F x
ax

,
 
ln ln ln
F x
a b
x


Метод наименьших квадратов применяется для определения коэффициентов линейной функции
1 0
ˆ
F
a X
a


, аппроксимирующей экспериментальные данные


,
i
i
X Y , где ln
i
i
X
x

, ln
i
i
Y
y

,
0
ln
a
a

,
1
a
b

,
ˆ
ln
F
F

,
1,
i
n

Для определения параметров линейной зависимости
1 0
ˆ
F
a X
a


составляется вспомогательная таблица:
x
y
X =lnx
Y=lny
X
2
XY
F
(y – F)
2
1 2,1 0,0
-0,5 0,0 0
2,13 0,00 2
6,2 0,7 0,6 0,5 0,407423 4,82 1,90 3
13,3 1,1 1,7 1,2 1,852699 10,90 5,76 4
24,6 1,4 2,3 1,9 3,164054 24,65 0,00 5
117,5 1,6 2,9 2,6 4,704647 55,73 3815,82

15
163,7
4,8
7,0
6,2
10,1
3823,493
Полученные значения сумм подставляются в систему линейных уравнений для определения параметров
0 1
,
a a
линейной аппроксимации:
1 0
1 0
6, 2 4,8 10,1 4,8 5
7, 0
a
a
a
a







Решением системы являются значения
1 2,14
a

и
0 0, 65
a
 
. Линеаризованная эмпирическая функция имеет вид: ˆ
2,14 0, 65
F
x


Определение параметров степенной эмпирической зависимости производится из ранее принятых соотношений:
0
ln
a
a

,
1
a
b

Откуда
0 0,65 0, 52
a
a
e
e




,
1 2,14
b
a


. Искомая степенная эмпирическая функция будет иметь вид:
 
2,14 0,52
b
F x
ax
x


Вопрос выбора наилучшей эмпирической зависимости
Если для таблицы


,
i
i
x y можно указать несколько классов эмпирических функций, то сначала из каждого класса методом наименьших квадратов ищется наилучшая функция, а затем из них выбирается та, которая даёт наименьшее среднеквадратическое отклонени е


,

y F .

Решение задачи 5 из Типового расчёта №2
Задание 5. Приближение табличных функций по методу наименьших
квадратов
Найти функцию, являющуюся наилучшим приближением к данной табличной функции по методу наименьших квадратов.
x y
1 2,6 2 6,4 3 19 4 38,9 5 85,8 6 322,3 7 796,9
Для исследования в MS Excel использовать
1) линейную функцию;
2) полиномиальную функцию 2-й степени;
3) экспоненциальную функцию;
4) логарифмическую функцию;
5) степенную функцию.
Выбрав наиболее подходящую зависимость, найти прогнозное значение
 
8
F
Решение
Порядок работы:
1.
Ввести двумерный массив данных


,
i
i
x y ,
1, 7
i

в таблицу MS Excel.
2.
Сформировать массивы предполагаемых эмпирических зависимостей, используя абсолютные ссылки на параметры функций (пусть ячейки с этими параметрами пока не заполнены).
3.
Сформировать суммы квадратов разностей значений эмпирической функции и соответствующих эмпирических данных.

4.
С помощью Надстройки MS Excel Поиск решения провести минимизацию целевых функций - сумм квадратов разностей, используя в решении метод ОПГ.
5.
Из полученных минимальных значений сумм квадратов выбрать наименьшее значение.
Соответствующую функцию можно принимать за искомую эмпирическую зависимость, по которой вычисляется прогнозное значение.
6.
Построить графики полученных эмпирических функций.

-400
-200 0
200 400 600 800 1000 1
2 3
4 5
6 7
F_линейн
F_поли y
F_эксп
F_лог
F_степ