Page 38 - Міністерство освіти і науки України
P. 38

де  матриця  P (N )  буде  пропорційною  коваріаційній  матриці
                            помилок при оцінюванні вектора параметрів        .
                            Позначимо  K    (N )  P (N )X  T (N ).  Тоді  з  (2.54)  дістанемо
                            рівняння
                             P  1 (  N  )  T (  N  )  ( N  )  T (  N  ) 1  ( N  ) 1  X  T (  N  ) X ( N  )

                               P  1 ( N  ) 1  X  T ( N  ) X ( N  ),
                            помноживши, яке спочатку зліва на  P     (N ), а потім справа на
                             P (N   ) 1 , знайдемо
                                        P (N    ) 1  P (N ) P (N )X  T (N )X (N )P (N  ) 1 . (2.56)
                            після  цього  помножимо  обидві  частини  рівняння  (2.56)  на
                             X  T  (N ) :
                             P (N   ) 1 X  T  (N )  P (N )X T (N ) P (N )X  T (N )X (N )P (N  ) 1 X  T (N )

                                            P (N )X  T (N  1 )  X (N )P (N  ) 1 X  T (N ).
                                   З  цього  рівняння  можна  визначити  вектор  корекції
                            (вагомих множників)
                                    K ( N  )  P ( N  ) X  T ( N  )
                                                                                     .   (2.57)
                                      P (  N  ) 1  X  T ( N  1 )  X ( N  P )  (  N  ) 1  X  T (  N  ).  1
                            З урахуванням (2.56), (2.57) маємо
                                   P ( N  )  P (  N  ) 1  P ( N  ) X  T ( N  ) X (  N  P )  (  N  ) 1

                                              P ( N   ) 1  P ( N  ) 1  X  T (  N  )     (2.58)

                                      1   X (  N  P )  (  N  ) 1  X  T (  N  )  1 X ( N  P )  ( N  1 ).

                                   На підставі виразів (2.55), (2.57), (2.58) можна записати
                            загальну  процедуру  (2.49)  для  визначення  середньої
                            квадратичної  оцінки  вектора  параметрів            за  РМНК.
                            Початковими  даними  тут  є  P )0(      I ,    0 (  )  поч   де     -
                            достатньо велике додатне число;  I  - одинична матриця;      поч   -
                            початкове  значення  вектора  параметрів,  яке  в  загальному
                            випадку можна вважати нульовим вектором.
   33   34   35   36   37   38   39   40   41   42   43