Скачиваний:
12
Добавлен:
13.06.2022
Размер:
1.7 Mб
Скачать

31 Тема 02: Построение эмпирических статистических моделей ХТП

Для определения (подгонки) коэффициентов (параметров) модели необходимо, чтобы критерий МНК стал наименьшим.

Задача определения коэффициентов модели сводится к реализации одного из алгоритмов оптимизации для определения минимума критерия МНК:

 

 

n

 

 

 

 

2

 

 

i

i

 

 

 

 

 

min

 

y

p

y

 

 

 

 

 

a a

доп

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

a

доп

 

- допустимая область изменения параметров рода.

a

- ограничения первого

Именно так решается задача параметрической идентификации для нелинейных моделей.

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

32 Тема 02: Построение эмпирических статистических моделей ХТП

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

Cr

0;

Cr

0; ...

Cr

0

a

a

a

 

 

 

m

 

0

 

1

 

 

 

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

.

a

0

, a

,...a

m

 

1

 

Однако, как показывает практика, решение систем нелинейных уравнений ничуть не проще, чем прямое решение задачи оптимизации с поиском минимума критерия МНК.

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

33 Тема 02: Построение эмпирических статистических моделей ХТП

Определение выборочных (эмпирических) коэффициентов регрессии для линейных по параметрам моделей (произвольное число входных переменных

x

x

s

s 1,...r

r 1

 

 

 

 

 

Таблица проведения экспериментальных исследований в этом случае имеет вид:

p

n

1

2

n

x

 

1

x

 

11

x

21

 

xn1

x

2

 

 

x

 

 

12

x

22

 

xn2

x

r

 

x

 

1r

x

2r

 

xnr

y

э

 

y1э

y2э

ynэ

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

34 Тема 02: Построение эмпирических статистических моделей ХТП

Для линейных или линеаризованных по параметрам моделей необходимо представить критерий МНК в виде:

n

 

m

 

 

2

 

 

 

 

 

 

 

 

Cr a j j x yi

 

 

i 1

j 0

 

 

 

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

(СЛАУ):

Cr

 

n

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

y

 

 

x

0

 

2

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

j

 

j

i

 

i

 

0

 

i

 

 

i 1

j 0

 

 

 

 

 

 

 

 

 

 

 

 

Cr

 

n

m

 

 

 

x

y

 

 

 

x 0

 

2

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

a1

 

j

 

j

i

 

i

1

 

i

 

 

 

i 1

j 0

 

 

 

 

 

 

 

 

 

 

 

 

............................................................

 

Cr

 

n

m

 

 

 

x

y

 

 

 

 

x 0

 

2

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

am

 

j

 

j

i

 

i

 

m

 

i

 

 

i 1

j 0

 

 

 

 

 

 

 

 

 

 

 

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

35 Тема 02: Построение эмпирических статистических моделей ХТП

Перегруппировав члены в последней системе уравнений, можно записать СЛАУ в виде:

m

 

n

 

 

x

 

 

x

n

 

 

x y

 

 

 

a

j

 

 

 

 

 

 

 

( j, u 0,1,...m)

j

u

u

i

 

 

i

 

i

 

i

 

j 0

 

i 1

 

 

 

 

 

 

i 1

 

 

 

 

 

И, если ввести в рассмотрение информационную матрицу

 

uj

j 0,1,...m; u 0,1,...m :

 

 

 

n

 

xi j xi

 

uj u

j 0,1,...m; u 0,1,...m

i 1

то она окажется квадратной, симметричной и значения её элементов зависят

только от входных переменных и конкретного вида функций

 

j

x

 

 

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

36 Тема 02: Построение эмпирических статистических моделей ХТП

В матричном виде информационную матрицу

 

можно представить в виде

 

произведения транспонированной и исходной матрицы входных переменных

 

 

:

T

 

 

 

Матрица, зависящая от входных переменных, имеет вид:

 

 

x

x

...

 

 

x

 

 

 

0

1

 

1

1

 

 

 

 

m

1

 

 

 

 

 

 

 

 

 

 

 

x

 

x

 

...

 

 

 

x

 

 

 

 

 

 

0

2

 

1

2

 

 

 

m

2

 

 

 

 

 

 

 

n m 1

...

 

...

 

... ...

 

 

 

 

x

 

 

x

 

...

 

 

 

x

 

 

 

m

 

 

 

0

n

 

1

n

 

 

 

n

 

Соответственно правую часть рассматриваемой СЛАУ можно записать:

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

37 Тема 02: Построение эмпирических статистических моделей ХТП

 

 

n

 

x y

 

 

b

 

 

 

 

(u 0,1,...m)

i

u

 

u

i

 

 

 

i 1

 

 

 

 

или в матричном виде:

T

 

 

b

y

 

В результате СЛАУ, решаемая для определения коэффициентов эмпирической модели, может быть представлена:

m

 

 

 

 

 

 

 

 

uj

a

j

b

(u 0,1,...m)

 

 

u

 

j 0

 

 

 

 

 

 

или в матричном виде:

a b

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

38 Тема 02: Построение эмпирических статистических моделей ХТП

Если для определения коэффициентов использовать метод обратной матрицы, то получится:

 

1

 

1

 

 

 

a

 

b

 

1

 

 

 

 

 

,

 

 

 

 

1

a

 

 

b ,

 

1

 

 

 

a

b

 

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

39 Тема 02: Построение эмпирических статистических моделей ХТП

Матричная формула для определения коэффициентов линейной регрессии (параметров эмпирической модели):

ˆ

 

 

 

T

 

 

1

 

 

T

y

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

• сформировать вектор наблюдений

y

 

 

и вычислить его компоненты (только для

линеаризованных моделей);

• сформировать и рассчитать компоненты матрицы, зависящей от входных переменных

;

 

 

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

40 Тема 02: Построение эмпирических статистических моделей ХТП

• транспонировать матрицу

 

T

 

 

;

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

• перемножить транспонированную матрицу

на исходную матрицу

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

• выполнить обращение информационной матрицы -

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

• умножить полученную обратную матрицу на матрицу

 

 

 

;

 

 

 

 

 

 

• умножить полученный результат на вектор наблюдений

коэффициенты регрессии

ˆ

a

и

y получить выборочные

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

Соседние файлы в папке Лекции ХТП