Статистические методы и модели. Костин В.Н - 54 стр.

UptoLike

Составители: 

Рубрика: 

54
В данном случае частные производные в выражении (3.13) будут равны
1
0
0
==
j
j
x
b
y
;
1
1
0
j
j
x
b
y
=
; … ;
jk
k
j
x
b
y
0
=
. (3.15)
Тогда система уравнений (3.13) с учетом (3.15) преобразуется к виду:
∑∑
∑∑
===
===
=
=
N
j
k
i
N
j
N
j
k
i
N
j
ji
i
jkjk
jэ
ji
i
jj
jэ
xbxxy
xbxxy
101
101
0
0
000
0
...
00
. (3.16)
Решение системы нормальных уравнений (3.16) целесообразно вести в
матричной форме. С этой целью представим ее в следующем виде:
B
XX
Y
X
TT
0
00
=
, (3.17)
где
0
X матрица входных переменных;
0
T
X транспонированная матрица к матрице
0
X
;
Y матрицастолбец выходного параметра;
B
матрица- столбец коэффициентов регрессии.
Для определения коэффициентов регрессии умножим обе части выра-
жения (3.17) на
1
0
0
XX
T
слева, тогда получим
BXXXXYXXX
TTTT
0
0
0
00
0
0
11
= ,
откуда
YXXXBE
TT
0
0
0
1
= ,
где
Еединичная матрица
или
YXXXB
TT
0
0
0
1
= , (3.18)
         В данном случае частные производные в выражении (3.13) будут равны

                        ∂y j              ∂y j 0          ∂y   0
                             = x j 0 = 1;     = x j1 ; … ; j = x jk .                                (3.15)
                        ∂b0               ∂b1             ∂bk

         Тогда система уравнений (3.13) с учетом (3.15) преобразуется к виду:

                              N                     N           k             0             
                             ∑ y jэ x j 0 − ∑ x j 0 ∑ bi x ji = 0 
                              j =1                  j =1       i =0      
                                                                         
                             ...                                         .                          (3.16)
                              N        0       N 0        k     0        
                             ∑
                             j =1
                                  y jэ x jk − ∑
                                              j =1
                                                   x jk ∑
                                                        i =0
                                                             bi x ji = 0 
                                                                         


      Решение системы нормальных уравнений (3.16) целесообразно вести в
матричной форме. С этой целью представим ее в следующем виде:
                                            0              0          0
                                           X Y = XT X B,
                                                T
                                                                                                     (3.17)
                    0
         где    X       – матрица входных переменных;
                0                                                                                0
               X T – транспонированная матрица к матрице X ;
                 Y – матрица – столбец выходного параметра;
                 B – матрица- столбец коэффициентов регрессии.

         Для определения коэффициентов регрессии умножим обе части выра-
                              −1
                  0 0 
жения (3.17) на  X T X           слева, тогда получим
                         

                                      −1                                          −1
                           0T 0           0
                                                   0 0                                0   0
                          X X
                                         X Y =  X T X 
                                                T
                                                                                       XT X B,
                                                        

откуда
                                                                −1
                                             0T 0                           0
                                      BE =  X X                       XT Y ,
                                                  

         где   Е – единичная матрица

или
                                                               −1
                                            0 0                         0
                                      B =  X T X                  XT Y ,                         (3.18)
                                                   

                                                                                                        54