Convert between partial fraction expansion and ratio of two polynomials
The inputs to
residue are vectors of coefficients
of the polynomials
b = [bn ... b1 b0] and
= [am ... a1 a0]. The outputs are the residues
= [rm ... r2 r1], the poles
p = [pm ... p2 p1],
and the polynomial
k. For most textbook problems,
Find the partial fraction expansion of the following ratio of polynomials F(s) using
b = [ 5 3 -2 7]; a = [-4 0 8 3]; [r,p,k] = residue(b,a)
r = -1.4167 -0.6653 1.3320 p = 1.5737 -1.1644 -0.4093 k = -1.2500
This represents the partial fraction expansion
Convert the partial fraction expansion back to polynomial coefficients using
[b,a] = residue(r,p,k)
b = -1.2500 -0.7500 0.5000 -1.7500 a = 1.0000 -0.0000 -2.0000 -0.7500
The result is normalized for the leading coefficient in the denominator and represents
The ratio of polynomials F(s) and its partial fraction expansion is
Find the partial fraction expansion of F(s) using
residue. When the degree of the numerator is greater than the degree of the denominator, output
k is a vector that represents the coefficients of a polynomial in s.
b = [2 0 0 1 0]; a = [1 0 1]; [r,p,k] = residue(b,a)
r = 0.5000 - 1.0000i 0.5000 + 1.0000i p = 0.0000 + 1.0000i 0.0000 - 1.0000i k = 2 0 -2
k represents the polynomial
b— Coefficients of numerator polynomialvector of numbers
Coefficients of the polynomial in the numerator, specified as a vector of numbers representing the coefficients of the polynomial in descending powers of s.
Complex Number Support: Yes
a— Coefficients of denominator polynomialvector of numbers
Coefficients of the polynomial in the denominator, specified as a vector of numbers representing the coefficients of the polynomial in descending powers of s.
Complex Number Support: Yes
r— Residues of partial fraction expansioncolumn vector of numbers
Residues of partial fraction expansion, returned as a column vector of numbers.
p— Poles of partial fraction expansioncolumn vector of numbers
Poles of partial fraction expansion, returned as a column vector of numbers.
k— Direct termrow vector of numbers
Direct term, returned as a row vector of numbers that specify the coefficients of the polynomial in descending powers of s.
Consider the fraction F(s) of two polynomials b and a of degree n and m, respectively
The fraction F(s) can be represented as a sum of simple fractions
This sum is called the partial fraction expansion of F. The values rm,...,r1 are the residues, the values pm,...,p1 are the poles, and k(s) is a polynomial in s. For most textbook problems, k(s) is 0 or a constant.
The number of poles
n = length(a)-1 = length(r) = length(p)
The direct term vector is empty if
length(k) = length(b)-length(a)+1
p(i) = ... = p(i+j-1) is a pole of multiplicity
then the expansion includes terms of the form
residue first obtains the poles using
Next, if the fraction is nonproper, the direct term
deconv, which performs polynomial long
residue determines the residues
by evaluating the polynomial with individual roots removed. For repeated
resi2 computes the residues at the repeated
Numerically, the partial fraction expansion of a ratio of polynomials represents an ill-posed problem. If the denominator polynomial, a(s), is near a polynomial with multiple roots, then small changes in the data, including roundoff errors, can result in arbitrarily large changes in the resulting poles and residues. Problem formulations making use of state-space or zero-pole representations are preferable.
 Oppenheim, A.V. and R.W. Schafer. Digital Signal Processing. Prentice-Hall, 1975, p. 56.