2.4. Вычисление коэффициентов регрессионной модели и метод наименьших квадратов.
Метод вычисления
коэффициентов
в (2.10) называют методом наименьших
квадратов (сокращенно - МНК), если при
этом используется критерий наименьших
квадратов
.
Для упрощения
записи будем использовать скалярную
переменную
,
что соответствует
,
но все формулы верны, если
- это массив.
Запишем условия
минимума функции многих переменных
:
,
и определим эти
частные производные, дифференцируя
.
Получаем:

|
(2.12) |
Подставим сюда

,
,
изменим порядок
суммирования и обозначим через
значение
-го
регрессора в
-ой
точке. Получаем систему
уравнений с
неизвестными величинами
.
,

|
(2.13) |
Система является
линейной и решается стандартными
методами. В результате получаем
оптимальные коэффициенты
.
Вводя матрицу
коэффициентов системы (2.13),
,
и векторы
,
,
,
запишем её в компактном виде:

|
(2.14) |
Часто вводят
регрессионную матрицу
c элементами
.
;
транспонированная
матрица:
.
Она полезна при
вычислениях
,
,
где
,
.
Подчеркнем, что
полученная система линейна относительно
коэффициентов
,
а зависимость от входных параметров
может быть
произвольной.
Пример регрессионной
модели.
Пусть имеем
исходных точек,
входных параметра
,
регрессора:
,
,
,
что соответствует регрессионной модели
.
Решив систему трех линейных уравнений,
получим аппроксимирующую функцию с
конкретными числовыми коэффициентами.
Регрессорами могут
быть в частности полиномы или просто
степени
,
например:
,
,
,
,
т.е.
,
это аппроксимирующий полином.
Достоинства
регрессионных моделей: простота,
наглядность, представление данных в
многомерном пространстве, наличие
стандартных статистических методов
для проверки их адекватности. Недостатки:
простота, произвольный выбор регрессоров,
малая область адекватности.
Для аппроксимации
в MathCAD
есть ряд стандартных функций: slope,
intersept,
loess,
regress,
linfit,
genfit.
При выполнении аппроксимации можно
провести предварительное сглаживание
данных и для этого использовать
стандартные функции MathCAD,
например supsmooth.