11.4.2. Искусственный пример с двумя входами
Подгонка моделей более чем с одним входным рядом не вызывает принципиальных трудностей, если не считать увеличения числа анализируемых параметров.
Например, в случае двух входов можно записать модель как

с
,
где , , и — стационарные процессы. Для вычисления мы вначале вычислим при конкретных значениях параметров 
, (11.4.3)
и при конкретных значениях параметров 
. (11.4.4)
Тогда шум можно вычислить по формуле
, (11.4.5)
и, наконец, по формуле
. (11.4.6)
Искусственный пример. Ясно, что даже простая ситуация может привести к оцениванию большого числа параметров. В нижеследующем примере с двумя входными переменными и моделями первого порядка с запаздыванием имеется восемь неизвестных параметров. Для того чтобы определить, применима ли для оценок параметров в таких моделях нелинейная процедура наименьших квадратов, описанная в разд. 11.3.2, был проведен эксперимент с искусственными данными.

Рис. 11.8. Данные искусственного примера с двумя входами (ряд ).
Подробности эксперимента описаны в [85]. Данные генерировались моделью, -представление которой имело вид
(11.4.5)
с и . Входные переменные , изменялись в соответствии с рандомизированным -факториальным планом с тремя репликами. Предполагалось, что каждое значение на входе поддерживалось постоянным в течение 5 мин, а наблюдения на выходе делались каждую минуту. Данные представлены на рис. 11.8 и содержатся в виде ряда в сводке временных рядов в конце книги.
Для получения оценок наименьших квадратов использовалась итеративная нелинейная процедура с ограничениями, описанная в гл. 7. Необходимо было только вычислять . Так, для заданных значений параметров значения и можно было получить из
,
.
Они использовались для расчета :
.
Наконец, при заданных значения можно было вычислить по формуле
.
Предполагалось, что входы процесса поддерживаются в типичных условиях некоторое время до начала эксперимента, так что и и, следовательно, можно вычислять вперед с , а с .
Процедура опробовалась дважды, с разными наборами начальных значений. В первом расчете был взят набор параметров, который мог бы предложить человек, знакомый с изучаемым процессом. Во втором в качестве начального значения было взято среднее значение всех наблюдений, а остальные начальные значения были приравнены 0,1. Поэтому вторая попытка описывает значительно худшую ситуацию, чем обычно бывает на практике. В табл. 11.8 показано, что для первого набора начальных значений сходимость наступает после 5 итераций, а табл. 11.9 показывает, что для второго набора сходимость достигнута после 9 итераций. Эти результаты позволяют надеяться, что в реальных обстоятельствах оценивание системы с двумя входами не встретит серьезных трудностей.
Таблица 11.8. Сходимость нелинейного метода наименьших квадратов для системы с двумя входами при использовании предположительных начальных значений
Интерация
|

|

|

|

|

|

|

|

|
Сумма квадратов
|
0
1
2
3
4
5
6
7
8
|
59,19
59,20
59,24
59,35
59,41
59,39
59,39
59,39
59,39
|
10,00
9,07
8,38
9,24
11,90
12,03
12,08
12,07
12,07
|
-7,00
-6,37
-5,35
-3,98
-3,40
-3,52
-3,53
-3,53
-3,53
|
-0,50
-0,58
-0,70
-0,75
-0,75
-0,80
-0,79
-0,79
-0,79
|
-0,50
-0,56
-0,59
-0,55
-0,56
-0,57
-0,56
-0,56
-0,56
|
1,00
1,33
2,03
3,45
5,21
4,99
5,03
5,03
5,03
|
1,00
1,31
1,75
1,95
1,66
1,76
1,77
1,77
1,77
|
0,10
0,24
0,39
0,36
0,22
0,21
0,21
0,21
0,21
|
2046,8
1085,4
621,5
503,5
463,7
461,8
461,8
461,8
461,8
|
Таблица 11.9. Сходимость нелинейного метода наименьших квадратов для системы с двумя входами при использовании экстремальных начальных значений
Интерация
|

|

|

|

|

|

|

|

|
Сумма квадратов
|
0
1
2
3
4
5
6
7
8
9
10
11
|
59,19
59,19
59,22
59,21
59,21
59,31
59,61
59,47
59,41
59,39
59,39
59,39
|
0,10
0,24
1,62
1,80
3,01
6,17
15,83
10,31
11,89
12,07
12,07
12,07
|
0,10
-0,07
-0,29
-0,77
-1,31
-2,82
-3,25
-3,48
-3,41
-3,52
-3,53
-3,53
|
0,10
-1,51
-2,09
-1,75
-1,15
-0,93
-0,70
-0,74
-0,74
-0,79
-0,79
-0,79
|
0,10
1,77
-0,07
0,20
0,91
3,03
8,88
3,52
5,01
5,04
5,03
5,03
|
0,10
1,77
-0,07
0,20
0,91
3,03
8,88
3,52
5,01
5,04
5,03
5,03
|
0,10
-0,28
0,26
-0,10
0,22
1,20
1,64
1,63
1,65
1,76
1,77
1,77
|
0,10
0,15
0,29
0,56
0,72
0,67
0,26
0,23
0,20
0,21
0,21
0,21
|
2496,4
2190,5
1473,6
1016,8
743,1
611,4
534,2
501,9
462,8
461,8
461,8
461,8
|
|