Partial Fraction Decomposition Over The Reals

In mathematics, partial fractions are used in real-variable integral calculus to find real-valued antiderivatives of rational functions. The partial fraction decomposition of real rational functions is also used for Laplace transforms. This article describes the general method for obtaining the partial fraction decomposition of any real rational function. The proof of the existence of the partial fraction decomposition over an arbitrary field is not given here. For a sketch of the general proof, see partial fraction. For applications of partial fraction decomposition over the reals, see First we state the general result, and then we offer several examples to illustrate the result and show how the method is used in practice.

General result

Let f(x) be any rational function over the real numbers. In other words, suppose there exist real polynomials p(x) and q(x)≠ 0, such that
f(x) = \frac{p(x)}{q(x)}
By removing the leading coefficient of q(x), we may assume without loss of generality that q(x) is monic. By the fundamental theorem of algebra, we can write
q(x) = (x-a_1)^{j_1}\cdots(x-a_m)^{j_m}(x^2+b_1x+c_1)^{k_1}\cdots(x^2+b_nx+c_n)^{k_n}
where a1,..., am, b1,..., bn, c1,..., cn are real numbers with bi2 - 4ci < 0, and j1,..., jm, k1,..., kn are positive integers. The terms (x - ai) are the linear factors of q(x) which correspond to real roots of q(x), and the terms (xi2 + bix + ci) are the irreducible quadratic factors of q(x) which correspond to pairs of complex conjugate roots of q(x). Then the partial fraction decomposition of f(x) is the following.
f(x) = \frac{p(x)}{q(x)} = P(x) + \sum_{i=1}^m\sum_{r=1}^{j_i} \frac{A_{ir}}{(x-a_i)^r} + \sum_{i=1}^n\sum_{r=1}^{k_i} \frac{B_{ir}x+C_{ir}}{(x^2+b_ix+c_i)^r}
Here, P(x) is a (possibly zero) polynomial, and the Air, Bir, and Cir are real constants. There are a number of ways the constants can be found. The most straightforward method is to multiply through by the common denominator q(x). We then obtain an equation of polynomials whose left-hand side is simply p(x) and whose right-hand side has coefficients which are linear expressions of the constants Air, Bir, and Cir. Since two polynomials are equal if and only if their corresponding coefficients are equal, we can equate the coefficients of like terms. In this way, a system of linear equations is obtained which always has a unique solution. This solution can be found using any of the standard methods of linear algebra. However, this is often not the best way to go when computing by pencil and paper, and there are other ways ("tricks") to obtain the constants.

The algorithm

The partial fraction decomposition of a real rational function can be obtained by the following procedure:
  1. Remove the leading coefficient of the denominator q(x) if necessary.
  2. Reduce the rational function as much as possible; i.e. cancel any obvious common factors of the numerator and denominator.
  3. Make sure the degree of the numerator p(x) is strictly smaller than the degree of the denominator q(x). If it is not, then use the division algorithm for polynomials and long-divide the denominator into the numerator. You will get a (nonzero) polynomial quotient and a remainder whose degree is now strictly smaller than the degree of the denominator. Take this remainder as your new numerator and continue with the following steps. (Don't forget the polynomial quotient when you're done!)
  4. Factor the denominator q(x) into linear and irreducible quadratic factors. This is not a straightforward task; however, you have some tools to use. If you find any real root a of q(x), then you know that (x - a) is a factor of q(x). Also, if you differentiate q(x) and use the Euclidean algorithm to find the greatest common divisor of q ' (x) and q(x), then any real root a of the greatest common divisor gives you a factor (x - a). You can look for ways to complete the square in different powers of x. Most importantly, you can use the rational root test to determine all possible rational roots of q(x).
  5. Express the partial fraction decomposition of f(x) as the sum of double sums above.
  6. Find the constants Air, Bir, and Cir by whatever methods possible. (These methods are described below.) Once you have found all the constants, you now have the partial fraction decomposition. (Don't forget the polynomial quotient!)

Some examples

Example 1
f(x)=\frac{1}{x^2+2x-3}
Here, the denominator splits into two distinct linear factors:
q(x)=x^2+2x-3=(x+3)(x-1)
so we have the partial fraction decomposition
f(x)=\frac{1}{x^2+2x-3}=\frac{A}{x+3}+\frac{B}{x-1}
Multiplying through by x2 + 2x - 3, we have the polynomial identity
1=A(x-1)+B(x+3)
Substituting x = -3 and x = 1 into this equation gives A = -1/4 and B = 1/4, so that
f(x)=\frac{1}{x^2+2x-3}=\frac{1}{4}\left(\frac{-1}{x+3}+\frac{1}{x-1}\right)
Example 2
f(x)=\frac{x^3+16}{x^3-4x^2+8x}
After long-division, we have
f(x)=1+\frac{4x^2-8x+16}{x^3-4x^2+8x}=1+\frac{4x^2-8x+16}{x(x^2-4x+8)}
Since (-4)2-4(8) = -16 < 0, x2 - 4x + 8 is irreducible, and so
\frac{4x^2-8x+16}{x(x^2-4x+8)}=\frac{A}{x}+\frac{Bx+C}{x^2-4x+8}
Multiplying through by x3 - 4x2 + 8x, we have the polynomial identity
4x^2-8x+16 = A(x^2-4x+8)+(Bx+C)x
Taking x = 0, we see that 16 = 8A, so A = 2. Comparing the x2 coefficients, we see that 4 = A + B = 2 + B, so B = 2. Comparing linear coefficients, we see that -8 = -4A + C = -8 + C, so C = 0. Altogether,
f(x)=1+2\left(\frac{1}{x}+\frac{x}{x^2-4x+8}\right)
The following example illustrates almost all the "tricks" one would need to use short of consulting a computer algebra system. Example 3
f(x)=\frac{x^9-2x^6+2x^5-7x^4+13x^3-11x^2+12x-4}{x^7-3x^6+5x^5-7x^4+7x^3-5x^2+3x-1}
After long-division and factoring, we have
f(x)=x^2+3x+4+\frac{x^5-2x^4+5x^3-5x^2+6x-1}{(x-1)^3(x^2+1)^2}
The partial fraction decomposition takes the form
\frac{x^5-2x^4+5x^3-5x^2+6x-1}{(x-1)^3(x^2+1)^2}=\frac{A}{x-1}+\frac{B}{(x-1)^2}+\frac{C}{(x-1)^3}+\frac{Dx+E}{x^2+1}+\frac{Fx+G}{(x^2+1)^2}
Multiplying through by (x - 1)3(x2 + 1)2 we have the polynomial identity
x^5-2x^4+5x^3-5x^2+6x-1
=A(x-1)^2(x^2+1)^2+B(x-1)(x^2+1)^2+C(x^2+1)^2+(Dx+E)(x-1)^3(x^2+1)+(Fx+G)(x-1)^3
Taking x = 1 gives 4 = 4C, so C = 1. Similarly, taking x = i gives 2 + 2i = (Fi + G)(2 + 2i), so Fi + G = 1, so F = 0 and G = 1 by equating real and imaginary parts. We now have the identity
x^5-2x^4+5x^3-5x^2+6x-1
=A(x-1)^2(x^2+1)^2+B(x-1)(x^2+1)^2+(x^2+1)^2+(Dx+E)(x-1)^3(x^2+1)+(x-1)^3
Taking constant terms gives E = A - B + 1, taking leading coefficients gives A = -D, and taking x-coefficients gives B = 3 - D - 3E. Putting all of this together, E = A - B + 1 = -D - (3 - D - 3E) + 1 = 3E - 2, so E = 1 and A = B = -D. Now,
x^5-2x^4+5x^3-5x^2+6x-1
=A(x-1)^2(x^2+1)^2+A(x-1)(x^2+1)^2+(x^2+1)^2+(-Ax+1)(x-1)^3(x^2+1)+(x-1)^3
Taking x = -1 gives -20 = -8A - 20, so A = B = D = 0. The partial fraction decomposition of f(x) is thus
f(x)=x^2+3x+4+\frac{1}{(x-1)^3}+\frac{1}{x^2+1}+\frac{1}{(x^2+1)^2}

 

<< PreviousWord BrowserNext >>
r. l. stine
election day
charles lenox remond
the wild wild west
rodolphe kreutzer
tri county journal
list of mughal emperors
charles cotesworth pinckney
charles pinckney (father)
charles pinckney
eliza lucas
meramec valley current
dunkard brethren
madden nfl
microhistory
elections canada
yahya ayyash
dorothy loudon
washington missourian
uplink (computer game)
department of conservation
list of alleged al qaeda members
fox trot
ambrosia software
list of world curling men's champions
department of internal affairs
faith evans
triumph (band)
new zealand department of conservation
david michelinie
apostolic christian church
rock and roll machine
big boi
epicyclic gearing
just a game
james l. miller, sr.
new zealand forest service
riverview (ottawa)
list of water sports
ladislaus ii
ladislaus
serin
riduan isamuddin
cor pulmonale