1 / 36

Lecture 10 - Nonlinear gradient techniques and LU Decomposition

Lecture 10 - Nonlinear gradient techniques and LU Decomposition. CVEN 302 June 24, 2002. Lecture’s Goals. Nonlinear Gradient technique LU Decomposition Crout’s technique Doolittle’s technique Cholesky’s technique. Nonlinear Equations.

Download Presentation

Lecture 10 - Nonlinear gradient techniques and LU Decomposition

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Lecture 10 - Nonlinear gradient techniques and LU Decomposition CVEN 302 June 24, 2002

  2. Lecture’s Goals • Nonlinear Gradient technique • LU Decomposition • Crout’s technique • Doolittle’s technique • Cholesky’s technique

  3. Nonlinear Equations • The nonlinear equations can be solved using a gradient technique. • The minimization technique calculates a positive scalar value and use a gradient to find the zero of multiple functions.

  4. Minimization algorithm • Calculate the square function. • h(x) = [S (f(x))2] • Calculate a scalar value • z0 = h(x) • Calculate the gradient • dx = - dh/dx

  5. Minimization algorithm • Multiple loops to convergence xnew = xold + dx ; z1 = h(xnew ); dif = z1 - z0; if dif > 0 dx = dx/2 xnew = xold + dx else end loop endif

  6. Program FFMIN • The program is adapted from the book to do a minimization of scalar and uses a gradient technique to find the roots.

  7. Example of the 2-D Problem f1(x,y) = x2 + y2 - 1 f2(x,y) = x2 - y

  8. Example of the 2-D Problem The gradient function: h(x,y) =[( x2 + y2 - 1)2 +( x2 - y) 2] The derivative of the function: dh=[-(4(x2 + y2-1)x + 4( x2 - y)x) -(4(x2 +y2 - 1)y - 2( x2 - y))]

  9. Example of the 3-D Problem f1(x,y,z) = x2 + 2y2 + 4z2 - 7 f2(x,y,z) = 2x2 + y3 + 6z2 - 10 f3(x,y,z) = xyz + 1

  10. Example of the 3-D Problem The gradient function: h(x,y,z) = [ (x2 + 2y2 + 4z2 - 7)2 + (2x2 + y3 + 6z2 - 10)2 + (xyz + 1) 2]

  11. End of material on Exam 1 Exam 1 Chapter 1 through 5 Monday July 3, 2002 open book and open notes

  12. Chapter 6 LU Decomposition of Matrices

  13. LU Decomposition A modification of the elimination method, called the LU decomposition. The technique will rewrite the matrix as the product of two matrices. A = LU

  14. LU Decomposition The technique breaks the matrix into a product of two matrices, L and U, L is a lower triangular matrix and U is an upper triangular matrix.

  15. LU Decomposition There are variation of the technique using different methods. • Crout’s reduction (U has ones on the diagonal) • Doolittle’s method( L has ones on the diagonal) • Cholesky’s method ( The diagonal terms are the same value for the L and U matrices)

  16. Decomposition

  17. LU Decomposition Solving Using the LU decomposition [A]{x} = [L][U]{x} = [L]{[U]{x}} = {b} Solve [L]{y} = {b} and then solve [U]{x} = {y}

  18. LU Decomposition The matrices are represented by

  19. Equation Solving What is the advantage of breaking up one linear set into two successive ones? • The advantage is that the solution of triangular set of equations is trivial to solve.

  20. Equation Solving • First step - forward substitution

  21. Equation Solving • Second step - back substitution

  22. LU Decomposition (Crout’s reduction) Matrix decomposition

  23. LU Decomposition (Doolittle’s method) Matrix decomposition

  24. Cholesky’s method Matrix is decomposed into: where, lii = uii

  25. LU Decomposition (Crout’s reduction) Matrix decomposition

  26. Crout’s Reduction The method alternates from solving from the lower triangular to the upper triangular

  27. Crout’s Reduction

  28. General formulation of Crout’s These are the general equations for the component of the two matrices

  29. Example The matrix is broken into a lower and upper triangular matrices.

  30. LU Decomposition (Doolittle’s method) Matrix decomposition

  31. Doolitte’s method The method alternates from solving from the upper triangular to the lower triangular

  32. General formulation of Doolittle’s The problem is reverse of the Crout’s reduction, starting with the upper triangular matrix and going to the lower triangular matrix.

  33. LU Programs • There are two programs • LU_factor - the program does a Doolittle decomposition of a matrix and returns the L and U matrices • LU_solve uses an L and U matrix combination to solve the system of equations.

  34. Example • The matrix is broken into a lower and upper triangular matrices.

  35. Summary • Nonlinear scalar gradient method uses a simple step to find the crossing terms. • Setup of the LU decomposition techniques.

  36. Homework • Check the Homework webpage

More Related