Page 277 - MATLAB an introduction with applications
P. 277

262 ———  MATLAB: An Introduction with Applications

                   Differentiating Eq.(5.1) w.r.t. x  gives
                                           i
                                  ∂ F  =− b +  A x
                                         i ∑
                                   x ∂  i   j  ij  j
                   or in vector notation

                                  ∇= −  b +  Ax                                                      ...(5.2)
                                   F
                   where ∇F is the gradient of F.
                   The gradient along u when the motion takes place along the line x = x  + su, where s is the distance moved
                                                                           0
                   is given by
                              ∇
                                 F  0 x + su  = − b +  A x +  su ) = ∇ F  0 x  +  s Au
                                             ( 0
                   If the change in the gradient sAu is perpendicular to a vector V, then
                                   VAu =  0                                                        ...(5.2A)
                                 T
                   The directions of u and V are said to be mutually conjugate.


                    5.3 NEWTON’S METHOD

                   Newton’s method is a gradient method and can be conveniently used to optimize functions with several
                   parameters. Let the function to be optimized be ()Ux , where  x  is the vector of parameters x , x , …, x .
                                                                                                1
                                                                                                  2
                                                                                                        n
                   The function  Ux can be expanded in the Taylor’s series about a point x as
                                ()
                                                                                *
                                              n  ∂ Ux * )
                                                   (
                                Ux =    (x * ) + ∑     (x −  x * i  )
                                  () U
                                                        i
                                              i= 1  ∂ X i
                   Newton’s method uses only two terms in the series. Expressing in concise form, the above series can be
                   written as
                                  ( ) Ux +
                                Ux =    (  * ) g − T  (x * )(x −  x * )                             ...(5.3)
                   where  (gx * ) is the vector of first derivatives given by
                                        ∂ Ux * ) ∂ Ux * )   T
                                                  (
                                           (
                                  ( gx * ) =   ,     ,.... 
                                          x ∂  1  x ∂  2  
                                   *
                   The minimum  ()Ux obtained by setting
                                  ∂ U  =  0
                                   x ∂  i

                   In Eq.(5.3) which yields the set of equations
                                       +
                                   ( gx * ) J  (x * )(x −  x * ) =  0                               ...(5.4)
   272   273   274   275   276   277   278   279   280   281   282