14.08.2013 Views

Rational functions - Cs.ioc.ee

Rational functions - Cs.ioc.ee

Rational functions - Cs.ioc.ee

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Rational</strong> Expansion Theorem Jaan Penjam<br />

7 Complex numbers vs real numbers<br />

Everything above is valid for the polynomials over complex fields, i.e. all variables and coefficients are<br />

complex numbers. The methods described work also for real numbers, but only if the polynomial Q(z)<br />

has ℓ real roots, where degQ(z) = ℓ. The assumption that a polynomial has as many roots as its degr<strong>ee</strong> is<br />

guarant<strong>ee</strong>d for complex numbers by the Fundamental Theorem of Algebra. This is the reason why you<br />

are strongly suggested to apply complex analysis in this context.<br />

However, it is possible to k<strong>ee</strong>p yourself to real numbers as well, but you n<strong>ee</strong>d to consider possibly<br />

much more complex calculations. For example, to decompose a rational function into partial fractions<br />

within real numbers, you have to consider, that a real polynomial Q(x) can have two types of irreducible<br />

factors:<br />

• x − ρ, where ρ is a real root of Q(x);<br />

• x 2 + px + q, where p and q are real numbers such that p 2 < 4q.<br />

Factors of both types can appear multiple times in the polynomial. Recall the techniques of such decomposition<br />

in notes of your undergraduate course of mathematical analysis, for example pages 111–118 in<br />

http://staff.ttu.<strong>ee</strong>/~janno/ma1chem.pdf.<br />

To use method of Decomposition into Partial Fractions parts of the rational function R(x) = P(x)/Q(x)<br />

are of the form<br />

•<br />

A1 A2<br />

x−ρ + (x−ρ) 2 + ··· + Ak<br />

(x−ρ) k , for a root ρ ∈ R with multiplicity k;<br />

• B1+C1x<br />

x2 B2+C2x<br />

+ +px+q (x2 +px+q) 2 + ··· + Bk+Ckx<br />

(x2 +px+q) k for a factor x2 + px + q with multiplicity k.<br />

All constants Ai,Bi,Ci should be found by solving the system of linear equations similarly to case of<br />

complex numbers (s<strong>ee</strong> Example 6.1).<br />

After decomposition, parts of the form A1<br />

x−ρ can be transformed into the poewer series similarly to<br />

the formula 2:<br />

<br />

1<br />

n + k − 1 xk =<br />

(x − ρ) k k − 1 ρk 1<br />

(−ρ) k (1 − 1<br />

ρ x)k = (−1)k · 1<br />

ρ k ∑ n0<br />

Unfortunately, in general, the parts of the form B+Cx<br />

x 2 +px+q<br />

does not have a simple transformation into<br />

power series. For particular numbers p and q the assistance of an one-line Taylor Series Calculator<br />

may be helpful. S<strong>ee</strong> for example, the Wolfram|Alpha Widgets on http://www.wolframalpha.com/<br />

widgets/view.jsp?id=f9476968629e1163bd4a3ba839d60925. Some first terms of decomposition<br />

of the fraction 1/(x 2 + px + q) looks like:<br />

1<br />

x2 1 px<br />

= −<br />

+ px + q q q2 + (p2 − q)x2 q3 − (p3 − 2pq)x3 q4 mmmmmmmmmmmmmmmmmmmjj − (p5 − 4p 3 q + 3pq 2 )x 5<br />

7<br />

q 6<br />

+ (p4 − 3p2q + q2 )x4 q5 −<br />

+ (p6 − 5p 4 q + 6p 2 q 2 − q 3 )x 6<br />

q 7<br />

− ···

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!