Перейти к содержанию

МатАнПрод:Интегрирование рациональных функций

Материал из Мадока ВТ Вики

Основные понятия

Определение: Многочленом (полиномом) Pn(x) степени n (n0, n) называется функция вида: Pn(x)=a0+a1x+a2x2++anxn, где ai, an0.

Определение: Рациональной функцией (рациональной дробью) называется функция вида Pn(x)Qm(x), где Pn(x) и Qm(x) - многочлены.

Определение: Рациональная функция Pn(x)Qm(x) называется правильной, если степень числителя меньше степени знаменателя: deg(Pn(x))<deg(Qm(x)). В противном случае (если nm) дробь называется неправильной.

Теорема (о делении многочленов с остатком): Если рациональная дробь Pn(x)Qm(x) является неправильной (nm), то существует единственное представление в виде: Pn(x)Qm(x)=Mnm(x)+Nk(x)Qm(x) где Mnm(x) - многочлен (целая часть), а Nk(x)Qm(x) - правильная рациональная дробь (k=deg(Nk(x))<m).

Определение: Число x0 называется корнем многочлена Qm(x), если Qm(x0)=0.

Теорема Безу: Число x0 является корнем многочлена Qm(x) тогда и только тогда, когда Qm(x) делится на (xx0) без остатка, т.е. Qm(x0)=0Qm(x)=(xx0)Qm1(x), где Qm1(x) - многочлен степени m1.

Теорема (о комплексных корнях многочлена с действительными коэффициентами): Если многочлен Qm(x) имеет действительные коэффициенты и число x0=α+iβ (β0) является его корнем, то сопряженное число x¯0=αiβ также является корнем Qm(x). Доказательство: Пусть Qm(x)=c0+c1x++cmxm, где ci. Если Qm(x0)=0, то c0+c1x0++cmx0m=0. Возьмем комплексное сопряжение от обеих частей: c0+c1x0++cmx0m=0 c0+c1x0++cmx0m=0 Так как ci, то ci=ci. Используя свойства сопряжения (a+b=a¯+b¯, ab=a¯b¯), получаем: c0+c1x¯0++cmx¯0m=0. Это означает, что Qm(x¯0)=0, т.е. x¯0 - корень Qm(x).

Основная теорема алгебры: Всякий многочлен степени m1 с действительными (или комплексными) коэффициентами имеет по крайней мере один корень в поле комплексных чисел .

Следствие: Любой многочлен Qm(x) степени m1 с действительными коэффициентами имеет ровно m корней в (с учетом их кратности).

Разложение многочлена на множители

Рассуждение:

  1. Пусть дан многочлен Qm(x) с действительными коэффициентами. По основной теореме алгебры, существует корень x1 такой, что Qm(x1)=0.
  2. По теореме Безу, Qm(x)=(xx1)Qm1(1)(x). 3. Применяя теорему Безу последовательно к Qm1(1)(x), Qm2(2)(x), …, получаем разложение на линейные множители над : Qm(x)=cm(xx1)(xx2)...(xxm) где x1,x2,...,xm - все корни многочлена Qm(x) (с учетом кратности), а cm - старший коэффициент многочлена Qm(x).
  3. Если Qm(x) имеет действительные коэффициенты, то его комплексные корни (β0) входят сопряженными парами. Пусть x0=α+iβ - корень, тогда x¯0=αiβ - тоже корень. В разложении над пара линейных множителей (xx0)(xx¯0) объединяется в один квадратичный множитель с действительными коэффициентами: (xx0)(xx¯0)=(x(α+iβ))(x(αiβ)) =((xα)iβ)((xα)+iβ) =(xα)2(iβ)2=(xα)2+β2 =x22αx+α2+β2 Обозначим p=2α и q=α2+β2. Тогда множитель имеет вид x2+px+q. Дискриминант этого квадратного трехчлена: D=p24q=(2α)24(α2+β2)=4α24α24β2=4β2. Так как β0, то D=4β2<0. Это означает, что квадратный трехчлен x2+px+q не имеет действительных корней и является неприводимым над полем .

  4. Таким образом, любой многочлен Qm(x) с действительными коэффициентами может быть разложен над в произведение своего старшего коэффициента cm, линейных множителей вида (xxk), соответствующих действительным корням xk, и квадратичных множителей вида (x2+pjx+qj) с отрицательным дискриминантом, соответствующих парам комплексно-сопряженных корней. Qm(x)=cm(xxr1)k1(xxrs)ks(x2+p1x+q1)l1(x2+ptx+qt)lt где ki+2lj=m.

Теорема о разложении многочлена над

Теорема: Любой многочлен Qm(x) степени m1 с действительными коэффициентами может быть представлен в виде произведения своего старшего коэффициента cm и множителей вида (xxi) и (x2+pjx+qj), где xi - действительные корни, а x2+pjx+qj - неприводимые над квадратные трехчлены (pj24qj<0), соответствующие парам комплексно-сопряженных корней. Будем считать cm=1 (если нет, можно вынести за скобки). Qm(x)=(xx1)k1(xx2)k2(xxp)kp(x2+p1x+q1)l1(x2+p2x+q2)l2(x2+psx+qs)ls где: * xi - различные действительные корни * ki - кратности действительных корней * pj,qj, pj24qj<0 (квадратные трехчлены неприводимы) * lj - кратности пар комплексно-сопряженных корней * Сумма кратностей: i=1pki+2j=1slj=m=deg(Qm(x)).

Разложение правильной рациональной дроби на простейшие дроби

Лемма 1: Пусть Pn(x)Qm(x) - правильная рациональная дробь (n<m). Пусть a является корнем знаменателя Qm(x) кратности k1, т.е. Qm(a)=0, и Qm(x)=(xa)kQmk(1)(x), где Qmk(1)(a)0. Тогда существует единственное представление: Pn(x)Qm(x)=Ak(xa)k+P~1(x)(xa)k1Qmk(1)(x) где Ak - константа, а P~1(x)(xa)k1Qmk(1)(x) - также правильная рациональная дробь.

Нахождение коэффициента Ak: Умножим обе части равенства на (xa)k: Pn(x)Qmk(1)(x)=Ak+P~1(x)(xa)Qmk(1)(x) Положим x=a. Так как Qmk(1)(a)0, получим: Pn(a)Qmk(1)(a)=Ak+0 (!) Формула для Ak: Ak=Pn(a)Qmk(1)(a)

Доказательство существования и единственности (кратко): Из Pn(x)Qm(x)Ak(xa)k=Pn(x)AkQmk(1)(x)(xa)kQmk(1)(x) Выберем Ak=Pn(a)Qmk(1)(a). Тогда числитель R(x)=Pn(x)AkQmk(1)(x) имеет корень при x=a, так как R(a)=Pn(a)Pn(a)Qmk(1)(a)Qmk(1)(a)=0. Следовательно, R(x)=(xa)P~1(x) по теореме Безу. Тогда: (xa)P~1(x)(xa)kQmk(1)(x)=P~1(x)(xa)k1Qmk(1)(x) Можно показать, что полученная дробь является правильной. Единственность Ak доказывается от противного: предположим два разложения с Ak,1 и Ak,2, приведем к общему знаменателю и подставим x=a, что даст Ak,1=Ak,2. * Равенство A1Q(1)(x)+P~1(x)(xa)=A2Q(1)(x)+P~2(x)(xa) при x=a дает A1Q(1)(a)=A2Q(1)(a), откуда A1=A2, так как Q(1)(a)0.

Следствие (Применение Леммы 1 k-раз): Повторяя эту процедуру k раз для корня a, мы выделим все слагаемые, соответствующие этому корню: Pn(x)Qm(x)=Ak(xa)k+Ak1(xa)k1++A1xa+P^(x)Qmk(1)(x) где Ai, а P^(x)Qmk(1)(x) - правильная рациональная дробь, знаменатель которой уже не имеет корня a (Qmk(1)(a)0).

Лемма 2: Пусть Pn(x)Qm(x) - правильная рациональная дробь (n<m). Пусть (x2+px+q) - неприводимый над множитель (p24q<0) знаменателя Qm(x) кратности k1. То есть Qm(x)=(x2+px+q)kQ~m2k(x), где Q~m2k(x0)0 для корней x0=α±iβ уравнения x2+px+q=0. Тогда существует единственное представление: Pn(x)Qm(x)=Akx+Bk(x2+px+q)k+P~1(x)(x2+px+q)k1Q~m2k(x) (**) где Ak,Bk - константы, а вторая дробь - правильная рациональная дробь.

Нахождение коэффициентов Ak,Bk: Умножим обе части равенства на (x2+px+q)k: Pn(x)Q~m2k(x)=Akx+Bk+P~1(x)(x2+px+q)Q~m2k(x) Пусть x0=α+iβ - один из комплексных корней уравнения x2+px+q=0. Подставим x=x0: Pn(x0)Q~m2k(x0)=Akx0+Bk+0 Обозначим комплексное число R=Pn(x0)Q~m2k(x0)=Re(R)+iIm(R). Ak(α+iβ)+Bk=Re(R)+iIm(R) (Akα+Bk)+i(Akβ)=Re(R)+iIm(R) Приравнивая действительные и мнимые части, получаем систему для Ak,Bk: {Akα+Bk=Re(R)Akβ=Im(R) Так как β0, система имеет единственное решение: Ak=Im(R)β Bk=Re(R)Akα=Re(R)αIm(R)β Так как Pn,Q~m2k имеют действительные коэффициенты, и если x0 - корень, то x¯0 - тоже, можно показать, что Ak,Bk всегда получаются действительными числами.

Доказательство существования и единственности (кратко): Аналогично случаю действительного корня. Рассмотрим разность Pn(x)(Akx+Bk)Q~m2k(x). Подстановка x=x0 дает Pn(x0)(Akx0+Bk)Q~m2k(x0)=Pn(x0)Pn(x0)Q~m2k(x0)Q~m2k(x0)=0. Значит x0 - корень этого полинома. Так как полином имеет действительные коэффициенты, то x¯0 тоже корень. Следовательно, полином делится на (xx0)(xx¯0)=x2+px+q. Pn(x)(Akx+Bk)Q~m2k(x)=(x2+px+q)P^(x). Это позволяет переписать исходную дробь в виде (**). Можно показать, что остаточная дробь правильная. Единственность доказывается аналогично случаю действительных корней, приравнивая два разложения и подставляя x=x0. * Равенство (A1x+B1)Q~(x)+P~1(x)(x2+px+q)=(A2x+B2)Q~(x)+P~2(x)(x2+px+q) при x=x0 дает (A1x0+B1)Q~(x0)=(A2x0+B2)Q~(x0). Так как Q~(x0)0, имеем A1x0+B1=A2x0+B2. Приравнивая действительные и мнимые части, получаем A1=A2 и B1=B2.

Следствие (Применение Леммы 2 k-раз): Повторяя эту процедуру k раз для множителя (x2+px+q), мы выделим все слагаемые, соответствующие этой паре корней: Pn(x)Qm(x)=Akx+Bk(x2+px+q)k+Ak1x+Bk1(x2+px+q)k1++A1x+B1x2+px+q+P^(x)Q~m2k(x) где Ai,Bi, а P^(x)Q~m2k(x) - правильная рациональная дробь, знаменатель которой уже не имеет корней x0,x¯0.

Общая теорема о разложении на простейшие дроби

Теорема: Любую правильную рациональную дробь Pn(x)Qm(x), знаменатель которой разложен на множители над вида Qm(x)=cm(xx1)k1(xxp)kp(x2+p1x+q1)l1(x2+psx+qs)ls можно единственным образом представить в виде суммы простейших дробей: Pn(x)Qm(x)=i=1p(Ai,ki(xxi)ki+Ai,ki1(xxi)ki1++Ai,1xxi)+j=1s(Bj,ljx+Cj,lj(x2+pjx+qj)lj+Bj,lj1x+Cj,lj1(x2+pjx+qj)lj1++Bj,1x+Cj,1x2+pjx+qj) где Ai,r,Bj,t,Cj,t - действительные коэффициенты. (Примечание: В конспекте коэффициенты для квадратичных множителей обозначались Ajx+Bj, здесь использовано Bjx+Cj для избежания конфликта с Ai,r для линейных множителей, что является стандартной практикой).

Интегрирование рациональной функции сводится к интегрированию многочлена (если дробь неправильная) и суммы простейших дробей. ivabus@celerrime-x ➜ ~ $ pandoc -f markdown -t mediawiki /Users/ivabus/Downloads/Матан\ №4.md -o -

Теорема о разложении правильной рациональной дроби

Теорема: Пусть Pn(x)Qm(x) - правильная рациональная дробь (n<m), и знаменатель Qm(x) (со старшим коэффициентом, равным 1) разложен на множители над : Qm(x)=(xx1)k1(xx2)k2(xxp)kp(x2+p1x+q1)l1(x2+p2x+q2)l2(x2+psx+qs)ls где xi, ki, pj,qj, pj24qj<0, lj. Тогда существует единственное разложение дроби Pn(x)Qm(x) в сумму простейших дробей вида: Pn(x)Qm(x)=i=1pr=1kiAi,r(xxi)r+j=1st=1ljBj,tx+Cj,t(x2+pjx+qj)t где Ai,r,Bj,t,Cj,t - некоторые константы.

(Примечание: В оригинальном конспекте использовалась несколько иная индексация сумм: i=1pj=1kiAi,j(xxi)kij+1+i=1sj=1liBi,jx+Ci,j(x2+pix+qi)lij+1. Представленная выше форма является стандартной и эквивалентной.)

Замечание: Если исходная дробь Pn(x)Qm(x) неправильная (nm), то сначала необходимо выделить целую часть (многочлен) с помощью деления “уголком”: Pn(x)Qm(x)=Mnm(x)+Nk(x)Qm(x), где k<m. Затем правильную дробь Nk(x)Qm(x) раскладывают на простейшие.

Интегрирование простейших дробей

Определение: Простейшими дробями I, II, III, IV типов называются соответственно дроби вида: I. Axa II. A(xa)k, где k,k2. III. Ax+Bx2+px+q, где p24q<0. IV. Ax+B(x2+px+q)k, где p24q<0, k,k2.

Интегралы от простейших дробей:

I. Axadx=Ad(xa)xa=Aln|xa|+C

  1. A(xa)kdx=A(xa)kd(xa)=A(xa)k+1k+1+C=A(1k)(xa)k1+C

  2. Ax+Bx2+px+qdx Знаменатель x2+px+q имеет комплексные корни, т.к. D=p24q<0. Выделим в числителе производную знаменателя (x2+px+q)=2x+p: Ax+B=A2(2x+p)Ap2+B=A2(2x+p)+(BAp2) Тогда интеграл разбивается на два: A2(2x+p)x2+px+qdx+BAp2x2+px+qdx Первый интеграл: A2d(x2+px+q)x2+px+q=A2ln(x2+px+q)+C1 (модуль не нужен, т.к. x2+px+q>0 при D<0). Второй интеграл: Выделим полный квадрат в знаменателе: x2+px+q=(x+p2)2+qp24=(x+p2)2+a2, где a2=qp24>0. (BAp2)dx(x+p/2)2+a2=(BAp2)1aarctan(x+p/2a)+C2 Объединяя: Ax+Bx2+px+qdx=A2ln(x2+px+q)+BAp/2qp2/4arctan(x+p/2qp2/4)+C

  3. Ax+B(x2+px+q)kdx, где k2. Аналогично типу III, выделяем производную знаменателя в числителе: A2(2x+p)(x2+px+q)kdx+BAp2(x2+px+q)kdx Первый интеграл: A2(x2+px+q)kd(x2+px+q)=A2(x2+px+q)k+1k+1+C1=A2(1k)(x2+px+q)k1+C1. Второй интеграл: (BAp2)dx((x+p/2)2+a2)k=(BAp2)dt(t2+a2)k, где t=x+p/2, a2=qp2/4. Обозначим Ik=dt(t2+a2)k. Этот интеграл вычисляется с помощью рекуррентной формулы.

Вывод рекуррентной формулы для Ik=dt(t2+a2)k: Ik=dt(t2+a2)k=1a2a2+t2t2(t2+a2)kdt Ik=1a2(t2+a2(t2+a2)kdtt2(t2+a2)kdt) Ik=1a2(Ik1tt(t2+a2)kdt) Интегрируем по частям ttdt(t2+a2)k: u=tdu=dt dv=tdt(t2+a2)kv=tdt(t2+a2)k=12d(t2+a2)(t2+a2)k=12(1k)(t2+a2)k1 tdv=uvvdu=t2(1k)(t2+a2)k1dt2(1k)(t2+a2)k1 ttdt(t2+a2)k=t2(1k)(t2+a2)k112(1k)Ik1 Подставляем обратно в выражение для Ik: Ik=1a2(Ik1[t2(1k)(t2+a2)k112(1k)Ik1]) Ik=1a2(Ik1(1+12(k1))+t2(k1)(t2+a2)k1) Ik=1a2(Ik12k2+12(k1)+t2(k1)(t2+a2)k1) Ik=2k32a2(k1)Ik1+12a2(k1)t(t2+a2)k1 Эта формула позволяет свести вычисление Ik к Ik1, затем к Ik2, и так далее, пока не дойдем до I1=dtt2+a2=1aarctan(ta)+C.

Вывод: Интеграл от любой рациональной функции Pn(x)Qm(x)dx всегда выражается через элементарные функции (многочлены, логарифмы, арктангенсы).