Advanced Engineering Mathematics
Advanced Engineering Mathematics
10th Edition
ISBN: 9780470458365
Author: Erwin Kreyszig
Publisher: Wiley, John & Sons, Incorporated
Bartleby Related Questions Icon

Related questions

Question
[8] [Numerical differentiation and round-off error]
In this problem we'll show that numerical differentiation can be unstable
in the presence of round off error.
Consider the 2nd order centered difference approximation to the first
derivative; in exact arithmetic we know that
|ƒ′ (xo) — ƒ (xo + h) — ƒ(xo — h)
_ = — | ≤ 1 42²
M
<
-h²
2h
6
where M = maxx=[a,b] |ƒ"" (x).
Assume now that h can be represented exactly (no floating point error),
but, in contrast, that f(x+h) and f(x−h) have some error in their floating
point representation. Denote these floating point approximations
f(xo+h): Fl(f(xo + h)), f(xoh): Fl(f(xo - h))
and suppose that €1₁, €2 are the floating
f(xo + h) = f(xo + h) + €₁,
(1)
ƒ'(xo)
point errors in each value; that is:
f(xo - h) = f(xo - h) + €2.
(a) Use the triangle inequality and (1) to show that the error in the
centered difference approximation in the presence of round error satisfies:
М
f(xo+h)-f(xo – h)
-
2
| S M N² + 1/
< -h²
2h
6
h'
where € =
|€2- €₁|/2.
(b) As h gets smaller, notice that part of the error from (a) will vanish,
while the other part will increase. Use differential calculus to find h such
that the error derived in part (a) will be minimized.
expand button
Transcribed Image Text:[8] [Numerical differentiation and round-off error] In this problem we'll show that numerical differentiation can be unstable in the presence of round off error. Consider the 2nd order centered difference approximation to the first derivative; in exact arithmetic we know that |ƒ′ (xo) — ƒ (xo + h) — ƒ(xo — h) _ = — | ≤ 1 42² M < -h² 2h 6 where M = maxx=[a,b] |ƒ"" (x). Assume now that h can be represented exactly (no floating point error), but, in contrast, that f(x+h) and f(x−h) have some error in their floating point representation. Denote these floating point approximations f(xo+h): Fl(f(xo + h)), f(xoh): Fl(f(xo - h)) and suppose that €1₁, €2 are the floating f(xo + h) = f(xo + h) + €₁, (1) ƒ'(xo) point errors in each value; that is: f(xo - h) = f(xo - h) + €2. (a) Use the triangle inequality and (1) to show that the error in the centered difference approximation in the presence of round error satisfies: М f(xo+h)-f(xo – h) - 2 | S M N² + 1/ < -h² 2h 6 h' where € = |€2- €₁|/2. (b) As h gets smaller, notice that part of the error from (a) will vanish, while the other part will increase. Use differential calculus to find h such that the error derived in part (a) will be minimized.
Expert Solution
Check Mark
Knowledge Booster
Background pattern image
Recommended textbooks for you
Text book image
Advanced Engineering Mathematics
Advanced Math
ISBN:9780470458365
Author:Erwin Kreyszig
Publisher:Wiley, John & Sons, Incorporated
Text book image
Numerical Methods for Engineers
Advanced Math
ISBN:9780073397924
Author:Steven C. Chapra Dr., Raymond P. Canale
Publisher:McGraw-Hill Education
Text book image
Introductory Mathematics for Engineering Applicat...
Advanced Math
ISBN:9781118141809
Author:Nathan Klingbeil
Publisher:WILEY
Text book image
Mathematics For Machine Technology
Advanced Math
ISBN:9781337798310
Author:Peterson, John.
Publisher:Cengage Learning,
Text book image
Basic Technical Mathematics
Advanced Math
ISBN:9780134437705
Author:Washington
Publisher:PEARSON
Text book image
Topology
Advanced Math
ISBN:9780134689517
Author:Munkres, James R.
Publisher:Pearson,