1.2 Аппроксимация полиномом первого порядка

Модель первого порядка описывается уравнением вида:

y=a•x+b

Для нахождения коэффициентов а и b составим систему линейных алгебраических уравнений, причем число уравнений в системе равно числу состояний объекта в эксперименте.

Для решения данной системы алгебраических уравнений воспользуемся матричным методом наименьших квадратов. Составим матрицы входных и выходных сигналов:

Получим систему с двумя неизвестными: X ^. A = Y

Транспонируем матрицу Х:

Умножив слева обе части исходной системы на транспонированную матрицу коэффициентов, получим систему, число уравнений в которой равно числу неизвестных, а решение этой системе будет доставлять минимум критерий оптимизации.

X^{T .} X ^. A = X^{T .} Y

Получим систему двух линейных алгебраических уравнений первого порядка:

Найдем главный определитель матрицы:

Найдем вспомогательные определители системы:

Найдем коэффициенты а и b:

Таким образом, получим полином:

у =0.428 . х - 0.198

Для оценки полученного полинома вычислим значения функции и сравним их с экспериментальными данными.

Результаты вычисления сведем в таблицу. таблица 2

Сумма квадратов отклонений:

у_i 2 = 0.174

Ниже приведен проверочный расчет модели объекта первого порядка на ЭВМ в системе MathCad.