Numerical Methods Bracketing Method
Numerical Methods Bracketing Method
Numerical Methods Bracketing Method
Where:
v = downward vertical velocity (m/s),
t = time (s),
g = the acceleration due to gravity (=9.81 m/s2),
cd = a lumped drag coefficient (kg/m),
m = the jumpers mass (kg)
Incremental Search
Value of m that
makes the
function equal to
zero?
Incremental Search
Illustration of a number of
general ways that a root may
occur in an interval prescribed
by a lower bound xl and
an upper bound xu . Parts (a)
and (c) indicate that if both
f (xl ) and f (xu) have the
same sign, either there will
be no roots or there will be an
even number of roots
within the interval. Parts (b)
and (d) indicate that if the
function has different signs at
the end points, there will
be an odd number of roots in
the interval.
Incremental Search
In general, if f (x) is real and continuous in the
interval from xl to xu and f (xl ) and f (xu) have
opposite signs, that is,
Bisection Method
Incremental search method in which the interval
is always divided in half. That is,
xr = (xl + xu) / 2
If a function changes sign over an interval, the
function value at the midpoint is evaluated
The location of the root is then determined as
lying within the subinterval where the sign
change occurs
The process is repeated until the root is known
to the required precision.
FEU Institute of Technology
Bisection Method
Problem
Use bisection method to solve for mass, m, in the
bungee jumper mathematical model. Assume
v=36, t=4, g=9.81 and cd=0.25. Use initial
guesses of xl=50 and xu=200. Tabulate the
values for xl, xu, xr, and the approximate percent
relative error for each iteration. Note also the
signs of f(xl), f(xu) and f(xr). Use a stopping
criterion of s = 0.5%.
Solution
10
FALSE POSITION
(also called the linear interpolation
method)
locates the root by joining f(xl ) and f(xu)
with a line
The intersection of this line with the x-axis
represents an improved estimate of the
root
11
FALSE POSITION
12
Problem
Use false position method to solve the
same (previous) problem approached with
bisection method.
13
Solution
xl
50.0000
xu
200.0000
xr
176.2773
|Ea|%
14
Solution
xl
50.0000
50.0000
xu
200.0000
176.2773
xr
176.2773
162.3828
|Ea|%
8.5566
15
Solution
xl
50.0000
50.0000
50.0000
50.0000
50.0000
50.0000
50.0000
50.0000
xu
200.0000
176.2773
162.3828
154.2446
149.4777
146.6856
145.0501
144.0922
xr
176.2773
162.3828
154.2446
149.4777
146.6856
145.0501
144.0922
143.5311
|Ea|%
8.5566
5.2762
3.1890
1.9035
1.1275
0.6648
0.3909
16
References
Applied Numerical Methods w/MATLAB
for Engineers & Scientist (3rd Edition).
Chapra. 2012.
Numerical Methods for Engineers (6th
Edition). Chapra and Canale. 2010.
17