Комплекс лабораторных работ по курсу "Моделирование систем". Давыдов Р.В - 8 стр.

UptoLike

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

=
Pnnn
P
P
ххх
ххх
ххх
х
,...,,
...................
,...,,
,...,,
10
21202
11101
(8)
Квадрат разности для 3-го опыта запишется так:
2
1100
2
)...()(
PjPjjj
j
j
xbxbxbyyy =
(9)
Подставляя зависимость (9) в выражение (6), получим:
=
=
n
j
PjPjjj
xbxbxbybS
1
2
1100
]...[)( min (10)
Необходимым условием минимума функции S (b ) является равенство нулю
ее частных производных по искомым параметрам (поскольку функция S (b ) явля-
ется квадратической, то эти условия выделяют единственную точку минимума).
0
0
=
b
S
,
0
1
=
b
S
, … ,
0=
P
b
S
или
=
=
=
=
=
=
n
j
PjPjPjjj
n
j
jPjPjjj
n
j
jPjPjjj
хxbxbxby
хxbxbxby
хxbxbxby
1
1100
1
11100
1
01100
0))(...(2
....................................................................
0))(...(2
0))(...(2
Запишем эту систему в виде, удобном для анализа,
=+++
=+++
=+++
∑∑
∑∑
∑∑
== ==
== ==
== = =
n
j
n
j
n
j
n
j
jPjPjPPjjPjj
n
j
n
j
n
j
n
j
jjPjjPjjj
n
j
n
j
n
j
n
j
jjPjjPjjJ
yxxbxxbxxb
yxxxbxbxxb
yxxxbxxbxb
11 11
2
1100
11 11
11
2
11100
11 1 1
00101
2
00
...
............................................................................
...
...
(11)
Полученная система линейных алгебраических уравнений содержит столько
уравнений, сколько в нее входит неизвестных параметров Б. В теории метода сис-
тему (11) принято называть системой нормальных уравнений.
8
                                                 ⎛ х01 , х11 ,..., х P1 ⎞
                                                 ⎜                      ⎟
                                                 ⎜ х    , х  ,..., х P2 ⎟
                                             х = ⎜ 02 12                                                (8)
                                                   ................... ⎟
                                                 ⎜                      ⎟
                                                 ⎜ х , х ,..., х ⎟
                                                 ⎝ 0 n 1n            Pn ⎠


     Квадрат разности для 3-го опыта запишется так:
                                         ∧
                            ( y j − y j ) 2 = ( y j − b0 x0 j − b1 x1 j − ... − bP x Pj ) 2             (9)
     Подставляя зависимость (9) в выражение (6), получим:
                                  n
                      S (b) = ∑ [ y j − b0 x0 j − b1 x1 j − ... − bP x Pj ]2 → min                     (10)
                                  j =1


     Необходимым условием минимума функции S ( b ) является равенство нулю
ее частных производных по искомым параметрам (поскольку функция S ( b ) явля-
ется квадратической, то эти условия выделяют единственную точку минимума).
                                    ∂S       ∂S                             ∂S
                                        = 0,     = 0,                 ,         =0
                                    ∂b0      ∂b1                            ∂bP
     или
                       ⎧2 n ( y − b x − b x − ... − b x )(− х ) = 0
                       ⎪ ∑  j =1
                                   j      0 0j        1 1j             P Pj          0j

                       ⎪ n
                       ⎪⎪2∑ ( y j − b0 x0 j − b1 x1 j − ... − bP x Pj )(− х1 j ) = 0
                        ⎨ j =1
                        ⎪....................................................................
                        ⎪ n
                        ⎪2∑ ( y − b x − b x − ... − b x )(− х ) = 0
                        ⎪⎩ j =1 j         0 0j        1 1j             P Pj          Pj



     Запишем эту систему в виде, удобном для анализа,
                    ⎧b n x 2 + b n x x + ... + b n x x = n x y
                    ⎪ 0∑  j =1
                                0J      1∑ 0 j 1j
                                          j =1
                                                                   P ∑ 0 j Pj
                                                                      j =1
                                                                                      ∑ j =1
                                                                                              0j   j

                    ⎪ n                         n                     n                n
                    ⎪⎪b0 ∑ x0 j x1 j + b1 ∑ x12j + ... + bP ∑ x1 j x Pj = ∑ x1 j y j
                     ⎨ j =1                    j =1                  j =1             j =1             (11)
                     ⎪............................................................................
                     ⎪ n                         n                          n           n
                     ⎪b ∑ x x + b ∑ x x + ... + b ∑ x 2 = ∑ x y
                     ⎪⎩ 0 j =1 0 j Pj        1
                                               j =1
                                                    1 j Pj              P
                                                                           j =1
                                                                                Pj
                                                                                      j =1
                                                                                             Pj  j



     Полученная система линейных алгебраических уравнений содержит столько
уравнений, сколько в нее входит неизвестных параметров Б. В теории метода сис-
тему (11) принято называть системой нормальных уравнений.


                                                        8