Введение в численные методы. Дулов Е.Н. - 51 стр.

UptoLike

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

51
( ) ( ) ( )
( )
ΦΦ+Φ=ΦΦ= ,,2,
,
2
fffffs
(8.3.2)
Наложим условие экстремума как равенство нулю всех частных производных:
( )
( )
0
,,
2
2
=
ΦΦ
+
Φ
=
jjj
fs
ααα
,
nj ,..,0=
(8.3.3)
В свою очередь
Φ
можно представить как:
=
=Φ
n
k
kk
0
ϕα
(8.3.4)
Тогда:
( )
( )
j
j
f
f
ϕ
α
,
,
=
Φ
(8.3.5)
Чтобы найти
( )
j
α
ΦΦ ,
, удобно представить
в виде
таблицы:
( ) ( )
( ) ( )
..
..
,,
..,,
11111010
10100000
ϕϕααϕϕαα
ϕϕααϕϕαα
(8.3.6)
Тогда, дифференцируя, к примеру, по
1
α
, можно видеть, что результатом будет
нулевой вклад всех столбцов и строк, кроме 1-х:
( )
( )
=
=
ΦΦ
n
i
jii
j
0
,2
,
ϕϕα
α
(8.3.7)
Подставляя (8.3.5) и (8.3.7) в (8.3.3), получаем СЛАУ:
( )
( )
j
n
i
jii
f
ϕϕϕα
,,
0
=
=
,
nj ,..,0=
(8.3.8)
В матричном виде:
( ) ( )
( ) ( )
( )
( )
=
.
.
,
,
.
.
..
..
,,
..,,
1
0
1
0
1110
1000
ϕ
ϕ
α
α
ϕϕϕϕ
ϕϕϕϕ
f
f
F=Γ
α
ˆ
F
1
ˆ
Γ=
α
(8.3.9)
Матрица
Γ
ˆ
называется матрицей Грама.
В итоге задача аппроксимации для полиномов сводится к решению СЛАУ или к
нахождению обратной матрицы и ее умножению на вектор.
    s 2 = (Φ − f , Φ − f ) = ( f , f ) − 2(Φ , f ) + (Φ , Φ )                                                   (8.3.2)
    Наложим условие экстремума как равенство нулю всех частных производных:
    ∂s 2      ∂ (Φ , f ) ∂ (Φ , Φ )
         = −2           +           = 0 , j = 0,.., n                                                           (8.3.3)
    ∂α j        ∂α j        ∂α j

    В свою очередь Φ можно представить как:
               n
    Φ = ∑α kϕ k                                                                                                 (8.3.4)
           k =0


    Тогда:
    ∂ (Φ , f )
               = ( f ,ϕ j )                                                                                     (8.3.5)
      ∂α j

                                        ∂ (Φ , Φ )
                                                                            (Φ , Φ ) = ∑∑α α (ϕ ,ϕ )
                                                                                    n     n
    Чтобы найти                                    , удобно представить                                           в виде
                                           ∂α j
                                                                                                i   j   i   j
                                                                                    j =0 i =0


таблицы:
     α 0α 0 (ϕ 0 ,ϕ 0 ) α 0α 1 (ϕ 0 ,ϕ1 ) .              .
                                                            
     α 0α 1 (ϕ 0 ,ϕ1 ) α 1α 1 (ϕ1 ,ϕ1 )                     
                                                                                                              (8.3.6)
              .                            .
                                                            
    
             .                                           . 

    Тогда, дифференцируя, к примеру, по α1 , можно видеть, что результатом будет
нулевой вклад всех столбцов и строк, кроме 1-х:
    ∂ (Φ , Φ )
               = 2∑α i (ϕ i ,ϕ j )
                    n
                                                                                                                (8.3.7)
       ∂α j       i =0


    Подставляя (8.3.5) и (8.3.7) в (8.3.3), получаем СЛАУ:

    ∑α (ϕ ,ϕ ) = ( f ,ϕ ) ,
     n

           i       i   j            j       j = 0,.., n                                                         (8.3.8)
    i =0


    В матричном виде:
     (ϕ 0 ,ϕ 0 )      (ϕ 0 ,ϕ1 )       .   .  α 0   ( f ,ϕ 0 )
                                                              
     (ϕ 0 ,ϕ1 )       (ϕ1 ,ϕ1 )                α 1   ( f ,ϕ1 )
     .                                 .       .  =  . 
                                                              
     .
                                           .  .   . 

    Γˆ α = F

    α = Γˆ −1 F                                                                                                 (8.3.9)
    Матрица Γ̂ называется матрицей Грама.
    В итоге задача аппроксимации для полиномов сводится к решению СЛАУ или к
нахождению обратной матрицы и ее умножению на вектор.
                                                                       51