numerical_differentiation

numerical_differentiation - Num erical Differentiation...

This preview shows pages 1–9. Sign up to view the full content.

Numerical Differentiation ° Estimate the derivatives (slope, curvature, etc.) of a function by using the function values at only a set of discrete points ° Ordinary differential equation (ODE) ° Partial differential equation (PDE) ° Represent the function by Taylor polynomials or Lagrange interpolation ° Evaluate the derivatives of the interpolation polynomial at selected (unevenly distributed) nodal points

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
0 ( ) ( ) Differentiation lim ( ) ( ) ( ) Difference ( ) x y f(x) dy f x x f x dx x x x y f x x f x dy x x x x dx ∆ → = +∆ - = +∆ - +∆ - = +∆ - A smaller step (or h) results in a smaller error x
x i-3 x i-2 x i-1 x i x i+1 x i+2 x i+3 h x = Evenly distributed points along the x-axis Distance between two neighboring points is the same, i.e. h. x 1 x 2 x 3 Unevenly distributed points along the x-axis

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Lagrange Interpolation ° 1st-order Lagrange polynomial (Two (Two-Point Point ) ) ( ) ( ) ( ) ( ) ( ) ( 1 0 1 0 0 1 0 1 1 1 0 0 1 x f x x x x x f x x x x x f x L x f L x f - - + + - - = + + = ° Second-order Lagrange polynomial (Three (Three- Point ) ) x ( f ) x x )( x x ( ) x x )( x x ( ) x ( f ) x x )( x x ( ) x x )( x x ( ) x ( f ) x x )( x x ( ) x x )( x x ( ) x ( f 2 1 2 0 2 1 0 1 2 1 0 1 2 0 0 2 0 1 0 2 1 2 - - - - + + - - - - + + + - - - - =
Formula ± Lagrange interpolation polynomial for unequally spaced data ) ( ) )( ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( 1 i 1 i i 1 i 1 i i i 1 i 1 i x f x x x x x f x L x f x L x f x L x f - + + + + + + + - - - - - - - = + + + + = ) ( ) )( ( ) )( ( ) ( ) )( ( ) )( ( ) )( ( 1 i i 1 i 1 i 1 i i 1 i i 1 i i 1 i i 1 i 1 i 1 i 1 i i 1 i x f x x x x x x x x x f x x x x x x x x x x x x + + + + - + + - + + - + + - + + - - - - - - + + - - - - + +

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Lagrange Interpolation L 2 (x)f(x 2 ) L 0 (x)f(x 0 ) L 1 (x)f(x 1 ) x 0 x 1 x 2
Formula ± Lagrange interpolation polynomial for unequally spaced data ) )( ( ) )( ( ) )( ( ) ( ) )( ( ) )( ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( 1 i i 1 i i 1 i 1 i i 1 i 1 i i 1 i 1 i i 1 i 1 i 1 i i i 1 i 1 i x x x x x x x x x x x x x f x x x x x x x x x f x f x L x f x L x f x L x f - - - - - - + + - - - - = + + = + + - + + - + + - - + + - + + + + - - ± First derivative ) )( ( ) ( i 1 i 1 i 1 i i 1 i 1 i x x x x x f - - + + + + - + + - + + ) )( ( ) ( ) )( ( ) ( ) )( ( ) ( ) ( i 1 i 1 i 1 i i 1 i 1 i 1 i i 1 i i 1 i 1 i i 1 i 1 i i 1 i 1 i i 1 i x x x x x x x 2 x f x x x x x x x 2 x f x x x x x x x 2 x f x f - - - - + + - - - - + + - - - - = + + + - + + + - + + + + - + + - + + - - + + -

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Lagrange polynomial (Two (Two-Point ) Point ) ) ( ) ( ) ( ) ( ) ( ) ( 1 0 0 1 1 1 0 0 1 x f x x x f x x x f x L x f L x f - - + - - = + = 0 1 1 0 x x x x
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}

What students are saying

• As a current student on this bumpy collegiate pathway, I stumbled upon Course Hero, where I can find study resources for nearly all my courses, get online help from tutors 24/7, and even share my old projects, papers, and lecture notes with other students.

Kiran Temple University Fox School of Business ‘17, Course Hero Intern

• I cannot even describe how much Course Hero helped me this summer. It’s truly become something I can always rely on and help me. In the end, I was not only able to survive summer classes, but I was able to thrive thanks to Course Hero.

Dana University of Pennsylvania ‘17, Course Hero Intern

• The ability to access any university’s resources through Course Hero proved invaluable in my case. I was behind on Tulane coursework and actually used UCLA’s materials to help me move forward and get everything together on time.

Jill Tulane University ‘16, Course Hero Intern