Nonlinear Equations and Continuous Optimization
|
|
- Morris Pearson
- 5 years ago
- Views:
Transcription
1 Nonlinear Equations and Continuous Optimization Sanzheng Qiao Department of Computing and Software McMaster University March, 2014
2 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems of Nonlinear Equations 5 Continuous Optimization 6 Software Packages
3 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems of Nonlinear Equations 5 Continuous Optimization 6 Software Packages
4 Problem setting Find roots f(x) = 0 Often, methods are iterative (roots cannot be found in finite number of steps). Example Compute square roots x 2 A = 0
5 Problem setting Find roots f(x) = 0 Often, methods are iterative (roots cannot be found in finite number of steps). Example Compute square roots x 2 A = 0 Find the side of the square whose area is A
6 Compute square roots Start with a rectangle whose one side is x c, then the other side is A/x c so that its area is A. Make the rectangle more square by setting the new side: x + = 1 ) (x c + Axc 2 Then x c = x + and iterate.
7 Compute square roots Start with a rectangle whose one side is x c, then the other side is A/x c so that its area is A. Make the rectangle more square by setting the new side: x + = 1 ) (x c + Axc 2 Then x c = x + and iterate. A better form (x c Axc ) x + = x c 1 2
8 Compute square roots Three issues to be addressed Initialization (x 0 ) Convergence (x k x?) and rate (how fast?) Termination
9 Initialization Write A in base 4: A = m 4 e, 0.25 m < 1 then A = m 2 e. Now we can assume 4 1 A < 1. Linear interpolation of f(a) = A at A = 0.25, 1.0: p(a) = (1+2A)/
10 Initialization (cont.) Initial error bound: Differentiating 1+2A A 3 with respect to A and then setting the derivative to zero to find the maximum, it can be shown that A (1+2A)/3 0.05
11 Initialization (cont.) Initial error bound: Differentiating 1+2A A 3 with respect to A and then setting the derivative to zero to find the maximum, it can be shown that A (1+2A)/ Initial value: x 0 = (1+2A)/3 Initial error: e
12 Convergence A relation between x k+1 and x k : x k+1 = 1 ) (x 2 k + Axk Denote the error e k = x k A, then the relation between e k+1 and e k : e k+1 = x k+1 A = 1 2 ( x k ) 2 A xk = 1 2 x k e2 k
13 Convergence A relation between x k+1 and x k : x k+1 = 1 ) (x 2 k + Axk Denote the error e k = x k A, then the relation between e k+1 and e k : e k+1 = x k+1 A = 1 2 ( x k ) 2 A xk = 1 2 x k e2 k It can be shown that 0.5 x k 1.0.
14 Convergence (cont.) Since the initial error e , e k e 2 k 1 e2k 0 (0.05)2k We have shown the convergence (e k 0 as k ).
15 Convergence (cont.) Since the initial error e , e k e 2 k 1 e2k 0 (0.05)2k We have shown the convergence (e k 0 as k ). How fast? Rate: quadratic e k+1 cek 2, each iteration doubles the accuracy.
16 Termination Recall: e k (0.05) 2k < 10 2k.
17 Termination Recall: e k (0.05) 2k < 10 2k. When k = 3, e k < When k = 4, e k <
18 Termination Recall: e k (0.05) 2k < 10 2k. When k = 3, e k < When k = 4, e k < Three iterations are enough for IEEE single precision (2 24 ). Four iterations are enough for IEEE double precision (2 53 ).
19 Example Compute 3
20 Example Compute 3 Scale: 3 =
21 Example Compute 3 Scale: 3 = Initial: x 0 = ( )/3 = 2.5/3
22 Example Compute 3 Scale: 3 = Initial: x 0 = ( )/3 = 2.5/3 Iterate: x n+1 = x n (x n 0.75/x n )/2 x 5 = x 4. n x n error < 10 16
23 Example Compute 3 Scale: 3 = Initial: x 0 = ( )/3 = 2.5/3 Iterate: x n+1 = x n (x n 0.75/x n )/2 x 5 = x 4. Scale back: x n x n error < 10 16
24 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems of Nonlinear Equations 5 Continuous Optimization 6 Software Packages
25 Generic algorithm If f(a) f(b) 0 and f(x) is continuous on [a, b], then f(x) has a root on [a, b]. while (b-a)>tol m = (a+b)/2; if f(a)*f(m)<=0 b = m; else a = m; end; end; r = (a + b)/2;
26 Generic algorithm Two problems in the generic algorithm: The while-loop may not terminate.
27 Generic algorithm Two problems in the generic algorithm: The while-loop may not terminate. When a and b are two neighboring floating-point numbers and (b-a)>tol,(a+b)/2 is rounded to either a or b.
28 Generic algorithm Two problems in the generic algorithm: The while-loop may not terminate. When a and b are two neighboring floating-point numbers and (b-a)>tol,(a+b)/2 is rounded to either a or b. Redundant function evaluations.
29 An improved algorithm fa = f(a); while (b-a)>tol + eps*max( a, b ) m = (a + b)/2; fm = f(m); if fa*fm<=0 b = m; else a = m; fa = fm; end; end; r = (a + b)/2;
30 An improved algorithm fa = f(a); while (b-a)>tol + eps*max( a, b ) m = (a + b)/2; fm = f(m); if fa*fm<=0 b = m; else a = m; fa = fm; end; end; r = (a + b)/2; Note: eps*max( a, b ) is about the distance between two consecutive floating-point numbers near max( a, b ). (ulp)
31 Convergence Since b k a k (b a)/2 k, x [a k, b k ], and x k = (a k + b k )/2, we have e k = x k x b k a k 2 = b a 2 k+1 0 In this case, e k+1 0.5e k. Improve accuracy by 1 bit per iteration or 1 decimal digit for every three or so iterations.
32 Convergence In general, linear convergence rate: for some constant c < 1. e k+1 ce k
33 Convergence In general, linear convergence rate: for some constant c < 1. e k+1 ce k Difficulty: Locate the interval [a, b].
34 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems of Nonlinear Equations 5 Continuous Optimization 6 Software Packages
35 Idea The tangent line of f(x) at x c : Set y = 0 y = f(x c )+(x x c )f (x c )
36 Newton s method Newton s method x + = x c f(x c) f (x c )
37 Newton s method Newton s method x + = x c f(x c) f (x c ) Example. Square root problem revisited, find a zero of f(x) = x 2 A x + = x c x2 c A = x c 1 ) (x c Axc 2x c 2
38 Complex case Example f(x) = x 2 + x + 1 (zeros ( 1±i 3)/2) i x i error 0 i i i i i i converge x 6 = x 5
39 Convergence No guarantee of convergence (unlike bisection). For example, f(x) = atan(x), x + = x c (1+x 2 c)atan(x c ) x 0 = 1.5 (> )
40 Convergence No guarantee of convergence (unlike bisection). For example, f(x) = atan(x), x + = x c (1+x 2 c )atan(x c) x 1 =
41 Convergence No guarantee of convergence (unlike bisection). For example, f(x) = atan(x), x + = x c (1+x 2 c )atan(x c) x 2 =
42 Convergence No guarantee of convergence (unlike bisection). For example, f(x) = atan(x), x + = x c (1+x 2 c)atan(x c ) x 3 =
43 Convergence f(x) = atan(x), x + = x c (1+x 2 c)atan(x c ) x 0 = 1.3 ( 1.3 < )
44 Convergence f(x) = atan(x), x + = x c (1+x 2 c)atan(x c ) x 1 =
45 Convergence f(x) = atan(x), x + = x c (1+x 2 c)atan(x c ) x 1 =
46 Convergence f(x) = atan(x), x + = x c (1+x 2 c)atan(x c ) x 1 =
47 Convergence Conditions for convergence (qualitative): x 0 close enough to x f (x) does not change sign near x f(x) is not too nonlinear near x Newton s method is a local method.
48 Convergence Conditions for convergence (qualitative): x 0 close enough to x f (x) does not change sign near x f(x) is not too nonlinear near x Newton s method is a local method. Difficulty: Finding x 0.
49 Hybrid methods Combining bisection and Newton s methods Bracketing interval [a, b], and x c = a or b if x + = x c f(x c )/f (x c ) [a, b] bracketing interval [a, x + ] or [x +, b] else m = (a+b)/2; bracketing interval [a, m] or [m, b] Termination criteria: Any one of (b k a k ) < δ f(x c ) < δ Too many function evaluations
50 Avoiding derivatives Approximation f (x c ) f(x c +δ c ) f(x c ) δ c
51 Avoiding derivatives Approximation f (x c ) f(x c +δ c ) f(x c ) δ c Choice of δ c Example. Secant method (δ c = x x c ) f (x c ) f(x c) f(x ) x c x
52 Secant method x + = x c x c x f(x c ) f(x ) f(x c)
53 Secant method x + = x c x c x f(x c ) f(x ) f(x c) Usually, the convergence rate (if it converges) is (1+ 5)/2 1.6 e k+1 ce 1.6, superlinear, between quadratic and linear. k
54 Zeros of a polynomial Finding the zeros of a polynomial p = x n + c n 1 x n c 1 x 1 + c 0 Many methods were proposed.
55 Zeros of a polynomial Finding the zeros of a polynomial p = x n + c n 1 x n c 1 x 1 + c 0 Many methods were proposed. The eigenvalues of its companion matrix c c 1 C(p) = c 2, det(xi C(p)) = p c n 1
56 Example The zeros of the polynomial x 3 1 are the eigenvalues of
57 Example The zeros of the polynomial x 3 1 are the eigenvalues of One real and two complex conjugate eigenvalues.
58 Note How to compute the eigenvalues of a matrix? Finding the zeros of a polynomial used to be the way of finding the eigenvalues of a matrix A.
59 Note How to compute the eigenvalues of a matrix? Finding the zeros of a polynomial used to be the way of finding the eigenvalues of a matrix A. Text book method: The eigenvalues of a matrix A are the zeros of its characteristic polynomial det(λi A).
60 Note Now, we have efficient and reliable methods for computing eigenvalues of a matrix. QR method, John G.F. Francis and Vera N. Kublanovskaya, late 1950s. We find the zeros of a polynomial by computing the eigenvalues of its companion matrix.
61 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems of Nonlinear Equations 5 Continuous Optimization 6 Software Packages
62 Problem setting f 1 (x 1,..., x n ) = 0 f 2 (x 1,..., x n ) = 0. f n (x 1,..., x n ) = 0 Denote f(x) = 0 f: vector-valued function x: vector
63 Newton s method where s c is the solution of x + = x c + s c f(x c )+J(x c )s c = 0, i.e., s c = J 1 (x c )f(x c ), where J(x c ) is the Jacobian of f at x c : [ ] fi J(x) = = x j f 1 x 1. f n x 1 f 1 x n.. f n x n
64 Example A system of nonlinear equations { x 2 1 x 2 2 = 0 2x 1 x 2 = 1, with starting point x 0 = [ 0 1 ] Solution: x 1 = x 2 = 1/ 2
65 Example f(x) = [ f1 f 2 ] = [ x 2 1 x 2 2 2x 1 x 2 1 ] The Jacobian is and [ ] 2x1 2x J(x) = 2 2x 2 2x 1 J(x 0 ) = [ ]
66 Example Step 1: x 1 = x 0 J 1 (x 0 ) f(x 0 )
67 Example Step 1: x 1 = x 0 J 1 (x 0 ) f(x 0 ) Solving for d 0 in J(x 0 )d 0 = f(x 0 ), we have [ ] 0.5 d 0 = 0.5 Thus x 1 = [ 0 1 ] [ ] = [ ]
68 Example Step 2: J(x 1 ) = x 2 = x 1 J 1 (x 1 ) f(x 1 ) [ ] [, f(x 1 ) = ]
69 Example Solving for d 1 in J(x 1 )d 1 = f(x 1 ), we have [ ] 0.25 d 1 = 0.25 Thus x 2 = [ ] [ ] = [ ]
70 Avoiding derivatives The jth column of J(x) f 1 x j. f n x j = f x j can be approximated by the difference f(x 1,..., x j +δ, x j+1,..., x n ) f(x 1,..., x n ) δ
71 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems of Nonlinear Equations 5 Continuous Optimization 6 Software Packages
72 Problem setting min f(x) or max f(x) x S x S x: vector f(x): objective function and real-valued S: support
73 Problem setting min f(x) or max f(x) x S x S x: vector f(x): objective function and real-valued S: support Find a zero of the gradient f(x) = f(x) x 1. f(x) x n
74 Newton s method View the gradient f(x) as a vector-valued function and apply the Newton s method for solving nonlinear systems. At current x c, find the correction s c : x + = x c + s c where s c is the solution of f(x c )+H(x c )S c = 0. The matrix H(x c ) (the Jacobian of the gradient at x c ) is called the Hessian of f at x c ( 2 f(x c )): H i,j = 2 f x i x j.
75 Example Minimizing f : R 2 R: f(x) = x3 1 3 x 1x x 2. Perform one iteration of Newton s method for minimizing f using the starting point [ ] 0 x 0 =. 1
76 Example Apply the Newton s method for finding a zero of the gradient [ x f(x) = 1 2 ] x2 2 2x 1 x The Hessian [ H(x) = 2 f(x) = 2x 1 2x 2 2x 2 2x 1 ]
77 Example Step 1: x 1 = x 0 2 f(x 0 ) 1 f(x 0 ) [ ] [ 0 0 1/2 = 1 1/2 0 [ ] 1/2 = 1/2 ] [ 1 1 ]
78 Outline 1 Introduction 2 Bisection Method 3 Newton s Method 4 Systems of Nonlinear Equations 5 Continuous Optimization 6 Software Packages
79 Software Packages IMSL zporc, zplrc, zpocc MATLAB roots, fzero NAG c02agf, c02aff NAPACK czero Octave fsolve
80 Summary Issues in an iterative method: Initialization, convergence and rate of convergence, termination. The example of computing square root Bisection method: Numerical termination problem Newton s method: Initial value, convergence problems Newton s method for systems of nonlinear equations, Jacobian matrix Newton s method for minimization, gradient and Hessian.
81 References [1 ] George E. Forsyth and Michael A. Malcolm and Cleve B. Moler. Computer Methods for Mathematical Computations. Prentice-Hall, Inc., Ch 7.
Numerical Integration
Numerical Integration Sanzheng Qiao Department of Computing and Software McMaster University February, 2014 Outline 1 Introduction 2 Rectangle Rule 3 Trapezoid Rule 4 Error Estimates 5 Simpson s Rule 6
More informationContinuous Optimization
Continuous Optimization Sanzheng Qiao Department of Computing and Software McMaster University March, 2009 Outline 1 Introduction 2 Golden Section Search 3 Multivariate Functions Steepest Descent Method
More informationElements of Floating-point Arithmetic
Elements of Floating-point Arithmetic Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Floating-point Numbers Representations IEEE Floating-point Standards Underflow
More informationSolving Ordinary Differential Equations
Solving Ordinary Differential Equations Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Initial Value Problem Euler s Method Runge-Kutta Methods Multistep Methods
More informationElements of Floating-point Arithmetic
Elements of Floating-point Arithmetic Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Floating-point Numbers Representations IEEE Floating-point Standards Underflow
More informationTHE SECANT METHOD. q(x) = a 0 + a 1 x. with
THE SECANT METHOD Newton s method was based on using the line tangent to the curve of y = f (x), with the point of tangency (x 0, f (x 0 )). When x 0 α, the graph of the tangent line is approximately the
More informationNonlinearity Root-finding Bisection Fixed Point Iteration Newton s Method Secant Method Conclusion. Nonlinear Systems
Nonlinear Systems CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Nonlinear Systems 1 / 27 Part III: Nonlinear Problems Not
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 5. Nonlinear Equations
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T Heath Chapter 5 Nonlinear Equations Copyright c 2001 Reproduction permitted only for noncommercial, educational
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More information3.1 Introduction. Solve non-linear real equation f(x) = 0 for real root or zero x. E.g. x x 1.5 =0, tan x x =0.
3.1 Introduction Solve non-linear real equation f(x) = 0 for real root or zero x. E.g. x 3 +1.5x 1.5 =0, tan x x =0. Practical existence test for roots: by intermediate value theorem, f C[a, b] & f(a)f(b)
More informationMotivation: We have already seen an example of a system of nonlinear equations when we studied Gaussian integration (p.8 of integration notes)
AMSC/CMSC 460 Computational Methods, Fall 2007 UNIT 5: Nonlinear Equations Dianne P. O Leary c 2001, 2002, 2007 Solving Nonlinear Equations and Optimization Problems Read Chapter 8. Skip Section 8.1.1.
More informationIntro Polynomial Piecewise Cubic Spline Software Summary. Interpolation. Sanzheng Qiao. Department of Computing and Software McMaster University
Interpolation Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Introduction 2 Polynomial Interpolation 3 Piecewise Polynomial Interpolation 4 Natural Cubic Spline
More informationNonlinearity Root-finding Bisection Fixed Point Iteration Newton s Method Secant Method Conclusion. Nonlinear Systems
Nonlinear Systems CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Nonlinear Systems 1 / 24 Part III: Nonlinear Problems Not all numerical problems
More informationIntro Polynomial Piecewise Cubic Spline Software Summary. Interpolation. Sanzheng Qiao. Department of Computing and Software McMaster University
Interpolation Sanzheng Qiao Department of Computing and Software McMaster University January, 2014 Outline 1 Introduction 2 Polynomial Interpolation 3 Piecewise Polynomial Interpolation 4 Natural Cubic
More information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationCS 450 Numerical Analysis. Chapter 5: Nonlinear Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationSOLVING EQUATIONS OF ONE VARIABLE
1 SOLVING EQUATIONS OF ONE VARIABLE ELM1222 Numerical Analysis Some of the contents are adopted from Laurene V. Fausett, Applied Numerical Analysis using MATLAB. Prentice Hall Inc., 1999 2 Today s lecture
More information15 Nonlinear Equations and Zero-Finders
15 Nonlinear Equations and Zero-Finders This lecture describes several methods for the solution of nonlinear equations. In particular, we will discuss the computation of zeros of nonlinear functions f(x).
More informationNumerical solutions of nonlinear systems of equations
Numerical solutions of nonlinear systems of equations Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan E-mail: min@math.ntnu.edu.tw August 28, 2011 Outline 1 Fixed points
More informationINTRODUCTION, FOUNDATIONS
1 INTRODUCTION, FOUNDATIONS ELM1222 Numerical Analysis Some of the contents are adopted from Laurene V. Fausett, Applied Numerical Analysis using MATLAB. Prentice Hall Inc., 1999 2 Today s lecture Information
More informationMATH 3795 Lecture 12. Numerical Solution of Nonlinear Equations.
MATH 3795 Lecture 12. Numerical Solution of Nonlinear Equations. Dmitriy Leykekhman Fall 2008 Goals Learn about different methods for the solution of f(x) = 0, their advantages and disadvantages. Convergence
More information446 CHAP. 8 NUMERICAL OPTIMIZATION. Newton's Search for a Minimum of f(x,y) Newton s Method
446 CHAP. 8 NUMERICAL OPTIMIZATION Newton's Search for a Minimum of f(xy) Newton s Method The quadratic approximation method of Section 8.1 generated a sequence of seconddegree Lagrange polynomials. It
More informationMATH 350: Introduction to Computational Mathematics
MATH 350: Introduction to Computational Mathematics Chapter IV: Locating Roots of Equations Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2011 fasshauer@iit.edu
More informationp 1 p 0 (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- cant method.
80 CHAP. 2 SOLUTION OF NONLINEAR EQUATIONS f (x) = 0 y y = f(x) (p, 0) p 2 p 1 p 0 x (p 1, f(p 1 )) (p 0, f(p 0 )) The geometric construction of p 2 for the se- Figure 2.16 cant method. Secant Method The
More informationSolution of Nonlinear Equations
Solution of Nonlinear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 14, 017 One of the most frequently occurring problems in scientific work is to find the roots of equations of the form f(x) = 0. (1)
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More informationFinding the Roots of f(x) = 0. Gerald W. Recktenwald Department of Mechanical Engineering Portland State University
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationFinding the Roots of f(x) = 0
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationMATH 350: Introduction to Computational Mathematics
MATH 350: Introduction to Computational Mathematics Chapter IV: Locating Roots of Equations Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2011 fasshauer@iit.edu
More informationComputational Methods CMSC/AMSC/MAPL 460. Solving nonlinear equations and zero finding. Finding zeroes of functions
Computational Methods CMSC/AMSC/MAPL 460 Solving nonlinear equations and zero finding Ramani Duraiswami, Dept. of Computer Science Where does it arise? Finding zeroes of functions Solving functional equations
More informationNonlinear Equations. Chapter The Bisection Method
Chapter 6 Nonlinear Equations Given a nonlinear function f(), a value r such that f(r) = 0, is called a root or a zero of f() For eample, for f() = e 016064, Fig?? gives the set of points satisfying y
More informationIntroductory Numerical Analysis
Introductory Numerical Analysis Lecture Notes December 16, 017 Contents 1 Introduction to 1 11 Floating Point Numbers 1 1 Computational Errors 13 Algorithm 3 14 Calculus Review 3 Root Finding 5 1 Bisection
More informationOptimization and Calculus
Optimization and Calculus To begin, there is a close relationship between finding the roots to a function and optimizing a function. In the former case, we solve for x. In the latter, we solve: g(x) =
More informationMA 8019: Numerical Analysis I Solution of Nonlinear Equations
MA 8019: Numerical Analysis I Solution of Nonlinear Equations Suh-Yuh Yang ( 楊肅煜 ) Department of Mathematics, National Central University Jhongli District, Taoyuan City 32001, Taiwan syyang@math.ncu.edu.tw
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 3 Lecture 3 3.1 General remarks March 4, 2018 This
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 2018 Petros Koumoutsakos, Jens Honore Walther (Last update: April 16, 2018) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationECE580 Exam 1 October 4, Please do not write on the back of the exam pages. Extra paper is available from the instructor.
ECE580 Exam 1 October 4, 2012 1 Name: Solution Score: /100 You must show ALL of your work for full credit. This exam is closed-book. Calculators may NOT be used. Please leave fractions as fractions, etc.
More informationCS 221 Lecture 9. Tuesday, 1 November 2011
CS 221 Lecture 9 Tuesday, 1 November 2011 Some slides in this lecture are from the publisher s slides for Engineering Computation: An Introduction Using MATLAB and Excel 2009 McGraw-Hill Today s Agenda
More informationOptimization. Next: Curve Fitting Up: Numerical Analysis for Chemical Previous: Linear Algebraic and Equations. Subsections
Next: Curve Fitting Up: Numerical Analysis for Chemical Previous: Linear Algebraic and Equations Subsections One-dimensional Unconstrained Optimization Golden-Section Search Quadratic Interpolation Newton's
More informationCLASS NOTES Computational Methods for Engineering Applications I Spring 2015
CLASS NOTES Computational Methods for Engineering Applications I Spring 2015 Petros Koumoutsakos Gerardo Tauriello (Last update: July 2, 2015) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationNumerical Methods. Root Finding
Numerical Methods Solving Non Linear 1-Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More informationNumerical Solution of f(x) = 0
Numerical Solution of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@pdx.edu ME 350: Finding roots of f(x) = 0 Overview Topics covered in these slides
More informationPART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435
PART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435 Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu
More informationAppendix 8 Numerical Methods for Solving Nonlinear Equations 1
Appendix 8 Numerical Methods for Solving Nonlinear Equations 1 An equation is said to be nonlinear when it involves terms of degree higher than 1 in the unknown quantity. These terms may be polynomial
More informationLecture 8. Root finding II
1 Introduction Lecture 8 Root finding II In the previous lecture we considered the bisection root-bracketing algorithm. It requires only that the function be continuous and that we have a root bracketed
More informationNumerical Analysis: Solving Nonlinear Equations
Numerical Analysis: Solving Nonlinear Equations Mirko Navara http://cmp.felk.cvut.cz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office 104a
More informationCS760, S. Qiao Part 3 Page 1. Kernighan and Pike [5] present the following general guidelines for testing software.
CS760, S. Qiao Part 3 Page 1 Testing 1 General Guidelines Kernighan and Pike [5] present the following general guidelines for testing software. Know exactly what you are testing and what results you expect.
More informationChapter 2 Solutions of Equations of One Variable
Chapter 2 Solutions of Equations of One Variable 2.1 Bisection Method In this chapter we consider one of the most basic problems of numerical approximation, the root-finding problem. This process involves
More information6.252 NONLINEAR PROGRAMMING LECTURE 10 ALTERNATIVES TO GRADIENT PROJECTION LECTURE OUTLINE. Three Alternatives/Remedies for Gradient Projection
6.252 NONLINEAR PROGRAMMING LECTURE 10 ALTERNATIVES TO GRADIENT PROJECTION LECTURE OUTLINE Three Alternatives/Remedies for Gradient Projection Two-Metric Projection Methods Manifold Suboptimization Methods
More informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #2 Instructions: This is an open notes exam, i.e., you are allowed to consult any textbook, your class notes, homeworks, or any of the handouts from us.
More informationApplied Mathematics Letters. Combined bracketing methods for solving nonlinear equations
Applied Mathematics Letters 5 (01) 1755 1760 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Combined bracketing methods for
More informationLine Search Methods. Shefali Kulkarni-Thaker
1 BISECTION METHOD Line Search Methods Shefali Kulkarni-Thaker Consider the following unconstrained optimization problem min f(x) x R Any optimization algorithm starts by an initial point x 0 and performs
More informationNUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR)
NUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR) Autumn Session UNIT 1 Numerical analysis is the study of algorithms that uses, creates and implements algorithms for obtaining numerical solutions to problems
More informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #2 Instructions: This is an open notes exam, i.e., you are allowed to consult any textbook, your class notes, homeworks, or any of the handouts from us.
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationChapter 1. Root Finding Methods. 1.1 Bisection method
Chapter 1 Root Finding Methods We begin by considering numerical solutions to the problem f(x) = 0 (1.1) Although the problem above is simple to state it is not always easy to solve analytically. This
More information17 Solution of Nonlinear Systems
17 Solution of Nonlinear Systems We now discuss the solution of systems of nonlinear equations. An important ingredient will be the multivariate Taylor theorem. Theorem 17.1 Let D = {x 1, x 2,..., x m
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationNumerical Analysis Fall. Roots: Open Methods
Numerical Analysis 2015 Fall Roots: Open Methods Open Methods Open methods differ from bracketing methods, in that they require only a single starting value or two starting values that do not necessarily
More informationRoot Finding (and Optimisation)
Root Finding (and Optimisation) M.Sc. in Mathematical Modelling & Scientific Computing, Practical Numerical Analysis Michaelmas Term 2018, Lecture 4 Root Finding The idea of root finding is simple we want
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationMath 273a: Optimization Netwon s methods
Math 273a: Optimization Netwon s methods Instructor: Wotao Yin Department of Mathematics, UCLA Fall 2015 some material taken from Chong-Zak, 4th Ed. Main features of Newton s method Uses both first derivatives
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationIntroduction to Numerical Analysis
Introduction to Numerical Analysis S. Baskar and S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai 400 076. Introduction to Numerical Analysis Lecture Notes
More informationSOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS BISECTION METHOD
BISECTION METHOD If a function f(x) is continuous between a and b, and f(a) and f(b) are of opposite signs, then there exists at least one root between a and b. It is shown graphically as, Let f a be negative
More informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
More informationSolution of nonlinear equations
Chapter 1 Solution of nonlinear equations This chapter is devoted the problem of locating roots of equations (or zeros functions). The problem occurs frequently in scientific work. In this chapter we are
More informationNUMERICAL METHODS USING MATLAB
NUMERICAL METHODS USING MATLAB Dr John Penny George Lindfield Department of Mechanical Engineering, Aston University ELLIS HORWOOD NEW YORK LONDON TORONTO SYDNEY TOKYO SINGAPORE Preface 1 An introduction
More informationSTOP, a i+ 1 is the desired root. )f(a i) > 0. Else If f(a i+ 1. Set a i+1 = a i+ 1 and b i+1 = b Else Set a i+1 = a i and b i+1 = a i+ 1
53 17. Lecture 17 Nonlinear Equations Essentially, the only way that one can solve nonlinear equations is by iteration. The quadratic formula enables one to compute the roots of p(x) = 0 when p P. Formulas
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationOutline. Math Numerical Analysis. Intermediate Value Theorem. Lecture Notes Zeros and Roots. Joseph M. Mahaffy,
Outline Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research
More information8 Numerical methods for unconstrained problems
8 Numerical methods for unconstrained problems Optimization is one of the important fields in numerical computation, beside solving differential equations and linear systems. We can see that these fields
More informationOutline. Scientific Computing: An Introductory Survey. Optimization. Optimization Problems. Examples: Optimization Problems
Outline Scientific Computing: An Introductory Survey Chapter 6 Optimization 1 Prof. Michael. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationNON-LINEAR ALGEBRAIC EQUATIONS Lec. 5.1: Nonlinear Equation in Single Variable
NON-LINEAR ALGEBRAIC EQUATIONS Lec. 5.1: Nonlinear Equation in Single Variable Dr. Niket Kaisare Department of Chemical Engineering IIT Madras NPTEL Course: MATLAB Programming for Numerical Computations
More informationNumerical Methods Lecture 3
Numerical Methods Lecture 3 Nonlinear Equations by Pavel Ludvík Introduction Definition (Root or zero of a function) A root (or a zero) of a function f is a solution of an equation f (x) = 0. We learn
More informationNonlinear Optimization
Nonlinear Optimization (Com S 477/577 Notes) Yan-Bin Jia Nov 7, 2017 1 Introduction Given a single function f that depends on one or more independent variable, we want to find the values of those variables
More informationComputational Methods. Solving Equations
Computational Methods Solving Equations Manfred Huber 2010 1 Solving Equations Solving scalar equations is an elemental task that arises in a wide range of applications Corresponds to finding parameters
More informationINTRODUCTION TO NUMERICAL ANALYSIS
INTRODUCTION TO NUMERICAL ANALYSIS Cho, Hyoung Kyu Department of Nuclear Engineering Seoul National University 3. SOLVING NONLINEAR EQUATIONS 3.1 Background 3.2 Estimation of errors in numerical solutions
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More information1.1: The bisection method. September 2017
(1/11) 1.1: The bisection method Solving nonlinear equations MA385/530 Numerical Analysis September 2017 3 2 f(x)= x 2 2 x axis 1 0 1 x [0] =a x [2] =1 x [3] =1.5 x [1] =b 2 0.5 0 0.5 1 1.5 2 2.5 1 Solving
More informationNumerical Methods in Informatics
Numerical Methods in Informatics Lecture 2, 30.09.2016: Nonlinear Equations in One Variable http://www.math.uzh.ch/binf4232 Tulin Kaman Institute of Mathematics, University of Zurich E-mail: tulin.kaman@math.uzh.ch
More informationcha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO
cha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO ROOTS OF EQUATIONS PT2.1 MOTIVATION Years ago, you learned to use the quadratic formula x = b ± b 2 4ac 2a to solve f(x) = ax 2 + bx + c = 0 (PT2.1) (PT2.2)
More informationSolving Non-Linear Equations (Root Finding)
Solving Non-Linear Equations (Root Finding) Root finding Methods What are root finding methods? Methods for determining a solution of an equation. Essentially finding a root of a function, that is, a zero
More informationLecture 7. Root finding I. 1 Introduction. 2 Graphical solution
1 Introduction Lecture 7 Root finding I For our present purposes, root finding is the process of finding a real value of x which solves the equation f (x)=0. Since the equation g x =h x can be rewritten
More informationE5295/5B5749 Convex optimization with engineering applications. Lecture 8. Smooth convex unconstrained and equality-constrained minimization
E5295/5B5749 Convex optimization with engineering applications Lecture 8 Smooth convex unconstrained and equality-constrained minimization A. Forsgren, KTH 1 Lecture 8 Convex optimization 2006/2007 Unconstrained
More informationECE 595, Section 10 Numerical Simulations Lecture 7: Optimization and Eigenvalues. Prof. Peter Bermel January 23, 2013
ECE 595, Section 10 Numerical Simulations Lecture 7: Optimization and Eigenvalues Prof. Peter Bermel January 23, 2013 Outline Recap from Friday Optimization Methods Brent s Method Golden Section Search
More information, b = 0. (2) 1 2 The eigenvectors of A corresponding to the eigenvalues λ 1 = 1, λ 2 = 3 are
Quadratic forms We consider the quadratic function f : R 2 R defined by f(x) = 2 xt Ax b T x with x = (x, x 2 ) T, () where A R 2 2 is symmetric and b R 2. We will see that, depending on the eigenvalues
More informationLecture 7 Unconstrained nonlinear programming
Lecture 7 Unconstrained nonlinear programming Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University,
More informationTrust Regions. Charles J. Geyer. March 27, 2013
Trust Regions Charles J. Geyer March 27, 2013 1 Trust Region Theory We follow Nocedal and Wright (1999, Chapter 4), using their notation. Fletcher (1987, Section 5.1) discusses the same algorithm, but
More informationChapter 6: Derivative-Based. optimization 1
Chapter 6: Derivative-Based Optimization Introduction (6. Descent Methods (6. he Method of Steepest Descent (6.3 Newton s Methods (NM (6.4 Step Size Determination (6.5 Nonlinear Least-Squares Problems
More informationOptimization. Totally not complete this is...don't use it yet...
Optimization Totally not complete this is...don't use it yet... Bisection? Doing a root method is akin to doing a optimization method, but bi-section would not be an effective method - can detect sign
More informationMATH 3795 Lecture 13. Numerical Solution of Nonlinear Equations in R N.
MATH 3795 Lecture 13. Numerical Solution of Nonlinear Equations in R N. Dmitriy Leykekhman Fall 2008 Goals Learn about different methods for the solution of F (x) = 0, their advantages and disadvantages.
More informationMathematical modelling Chapter 2 Nonlinear models and geometric models
Mathematical modelling Chapter 2 Nonlinear models and geometric models Faculty of Computer and Information Science University of Ljubljana 2018/2019 3. Nonlinear models Given is a sample of points {(x
More informationMath /Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined
Math 400-001/Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined using limits. As a particular case, the derivative of f(x)
More informationOutline. Additional Nonlinear Systems. Abstract. Finding Equilibrium Points Numerically. Newton s Method
Outline Finding Equilibrium Points Numerically Additional Nonlinear Systems James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University June 13, 2017
More informationECE133A Applied Numerical Computing Additional Lecture Notes
Winter Quarter 2018 ECE133A Applied Numerical Computing Additional Lecture Notes L. Vandenberghe ii Contents 1 LU factorization 1 1.1 Definition................................. 1 1.2 Nonsingular sets
More information