Пример 1
x(t)=a0+a1*√t
t |
x(t) |
|
t |
x(t) |
(t) |
1 |
1 |
0,68 |
9 |
1 |
70,66 |
2 |
3 |
3,29 |
10 |
3 |
82,21 |
3 |
5 |
5,33 |
11 |
5 |
93,77 |
4 |
7 |
7,05 |
12 |
7 |
105,32 |
5 |
9 |
8,55 |
13 |
9 |
116,88 |
6 |
10 |
9,91 |
14 |
10 |
128,44 |
7 |
11 |
11,16 |
15 |
11 |
139,99 |
11 |
46 |
12,329 |
19 |
46 |
186,22 |
Решение:
Задача решается в отношении неизвестных a0 и a1 .
а) min ∑7t=1 [a0 + a1√t – x(t)]2
сумма квадратичных отклонений (мера близости временного ряда)
Обозначим ∑7t=1 [a0 + a1√t – x(t)]2=f(a0, a1)
Тогда уравнения в нормальной форме будут иметь вид:
∂
F/∂a0
= 2 ∑7t=1[a0
+ a1√t
– x(t)] 1 = 0
∂F/∂a1 = 2 ∑7t=1 [a0 + a1√t – x(t)] √t = 0
Уравнения в каноническом вида:
( ∑7t=11) a0 + (∑7t=1 √t) a1 = ∑7t=1 x(t)
(∑7t=1√t) a0 + (∑7t=1 t) a1 = ∑7t=1 x(t) √t
7 a0 + 13,478 a1 = 46
13,478 a0 + 28 a1 = 101,626
Решая систему методом Крамера, получим
7 13,478
∆ = = 14,344
13,478 28
46 13,478
∆0 = = -81,715
101,626 28
7 46
∆1 = = 91,394
13,478 101,626
a0 = ∆0 / ∆ = -5,697
a1 = ∆1 / ∆ = 6,372
Таким
образом,
=
-5,697 + 6,372√t , для
t=1,2,3,4,5,6,7,11
б) Проводим аналогичные рассуждения для t=9,10,11,12,13,14,15,19
Уравнения в каноническом виде:
7 a0 + 24,161 a1 = 46
24,161 a0 + 84 a1 = 165,788
Решая систему методом Крамера, получим:
7 24,161
∆ = = 4,246
24,161 84
46 24,161
∆0 = = -141,604
165,788 84
7 46
∆1 = = 49,11
24,161 165,788
a0 = ∆0 / ∆ = -33,349
a1 = ∆1 / ∆ = 11,556
Таким образом, = -33,349 + 11,556√t , t = 9,10,11,12,13,14,15,19
t |
x(t) |
(t) |
t |
x(t) |
(t) |
1 |
1 |
0,68 |
9 |
1 |
70,66 |
2 |
3 |
3,29 |
10 |
3 |
82,21 |
3 |
5 |
5,33 |
11 |
5 |
93,77 |
4 |
7 |
7,05 |
12 |
7 |
105,32 |
5 |
9 |
8,55 |
13 |
9 |
116,88 |
6 |
10 |
9,91 |
14 |
10 |
128,44 |
7 |
11 |
11,16 |
15 |
11 |
139,99 |
11 |
46 |
12,329 |
19 |
46 |
186,22 |

(t)