Linearization

From Wikipedia, the free encyclopedia

Linearization in mathematics and its applications in general refers to finding the linear approximation to a function at a given point. In the study of dynamical systems, linearization is a method for assessing the local stability of an equilibrium point of a system of nonlinear differential equations. This method is used in fields such as engineering, physics, economics, and ecology.

[edit] With a Function y = f(x)

Linearizations of a function are a lines; ones that are usually used for purposes of calculation. Linearization an efficacious method used to approximate the output of function at any x = a based on the value and slope of a function y = f(x) at x = b, given that f(x) is continuous on [a,b] (or [b,a])and that a is close to b. In, short, linearization approximates the output of a function near x = a.

For example, you might know that \sqrt{4} = 2. But without a calculator, what would be a good approximation of \sqrt{4.001} = \sqrt{4 + .001}?

For any given function y = f(x), f(x) can be approximated if it is near a known continuous point. The most basic requisite is that, where La(x) is the linearization of f(x) at x = a, L(a) = f(a). The point-slope form of an equation forms an equation of a line, given a point (H,K) and slope M. The general form of this equation is: yK = M(xH).

Using the point (x,f(x)), La(x) becomes y = f(a) + M(xa). Because continuous functions are locally linear, the best slope to substitute in would be the slope of the line tangent to f(x) at x = a.

While the concept of local linearity applies the most to points arbitrarily close to x = a, those relatively close work relatively well for linear approximations. After all, a linearization is only an approximation. The slope M should be, most accurately, the slope of the tangent line at x = a.

An approximation of f(x) at (x, f(x))
Enlarge
An approximation of f(x) at (x, f(x))

Visually, the accompanying diagram shows the tangent line of f(x) at x. At f(x + h), where h is any small positive or negative value, f(x+h) is very nearly the value of the tangent line at the point (x + h,L(x + h)).

The final equation for the linearization of a function at x = a is:

y = f(a) + f'(a)(x - a)\,

For x = a, f(a) is f(x) at a. The derivative of f(x) is f'(x), and the slope of f(x) at a is f'(a).

Instead of a linear function, perhaps, would x not be better approximated with polynomial? This can be accomplished via Taylor series. This is not the purpose of linearization though.

[edit] Example

To find \sqrt{4.001}, we can use the fact that \sqrt{4} = 2. The linearization of f(x) = \sqrt{x} at x = a is y = \sqrt{a} + \frac{1}{2 \sqrt{a}}(x - a), because the function f'(x) = \frac{1}{2 \sqrt{x}} defines the slope of the function f(x) = \sqrt{x} at x. Plugging in a = 4, the linearization at 4 is y = 2 + \frac{x-4}{4}. In this case x = 4.001, so \sqrt{4.001} is approximately 2 + \frac{4.001-4}{4} = 2.00025. The true value is closer to 2.00024998, which is amazingly accurate.

[edit] With Analysis of Linear Functions

The analysis of linear functions is well defined, but most representations of actual systems are nonlinear. Linearization allows us to use tools for studying linear systems to analyze the behavior of a nonlinear function near a given point. The linearization of a function is the first order term of its Taylor expansion around the point of interest. For a system defined by the equation

\frac{d\bold{x}}{dt} = \bold{F}(\bold{x},t),

the linearized system can be written as

\frac{d\bold{x}}{dt} = D\bold{F}(\bold{x_0},t)  \cdot (\bold{x} - \bold{x_0})

where \bold{x_0} is the point of interest and D\bold{F}(\bold{x_0}) is the Jacobian of \bold{F}(\bold{x}) evaluated at \bold{x_0}.

In stability analysis, one can use the eigenvalues of the Jacobian matrix evaluated at an equilibrium point to determine the nature of that equilibrium. If all the eigenvalues are positive, the equilibrium is unstable; if they are all negative the equilibrium is stable; and if the values are of mixed signs, the equilibrium is a saddle point. Any complex eigenvalues will appear in complex conjugate pairs and indicate spiral (or circular if the real components are zero) around the equilibrium.