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

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

Материал из Мадока ВТ Вики
Нет описания правки
Нет описания правки
 
Строка 85: Строка 85:
'''Вывод рекуррентной формулы для <math display="inline">I_k = \int \frac{dt}{(t^2+a^2)^k}</math>:''' <math display="inline">I_k = \int \frac{dt}{(t^2+a^2)^k} = \frac{1}{a^2} \int \frac{a^2+t^2-t^2}{(t^2+a^2)^k} dt</math> <math display="inline">I_k = \frac{1}{a^2} \left( \int \frac{t^2+a^2}{(t^2+a^2)^k} dt - \int \frac{t^2}{(t^2+a^2)^k} dt \right)</math> <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} - \int t \cdot \frac{t}{(t^2+a^2)^k} dt \right)</math> Интегрируем по частям <math display="inline">\int t \cdot \frac{t dt}{(t^2+a^2)^k}</math>: <math display="inline">u=t \implies du=dt</math> <math display="inline">dv = \frac{t dt}{(t^2+a^2)^k} \implies v = \int \frac{t dt}{(t^2+a^2)^k} = \frac{1}{2} \int \frac{d(t^2+a^2)}{(t^2+a^2)^k} = \frac{1}{2(1-k)(t^2+a^2)^{k-1}}</math> <math display="inline">\int t dv = uv - \int v du = \frac{t}{2(1-k)(t^2+a^2)^{k-1}} - \int \frac{dt}{2(1-k)(t^2+a^2)^{k-1}}</math> <math display="inline">\int t \frac{t dt}{(t^2+a^2)^k} = \frac{t}{2(1-k)(t^2+a^2)^{k-1}} - \frac{1}{2(1-k)} I_{k-1}</math> Подставляем обратно в выражение для <math display="inline">I_k</math>: <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} - \left[ \frac{t}{2(1-k)(t^2+a^2)^{k-1}} - \frac{1}{2(1-k)} I_{k-1} \right] \right)</math> <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} \left(1 + \frac{1}{2(k-1)}\right) + \frac{t}{2(k-1)(t^2+a^2)^{k-1}} \right)</math> <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} \frac{2k-2+1}{2(k-1)} + \frac{t}{2(k-1)(t^2+a^2)^{k-1}} \right)</math> <math display="block"> I_k = \frac{2k-3}{2a^2(k-1)} I_{k-1} + \frac{1}{2a^2(k-1)} \frac{t}{(t^2+a^2)^{k-1}} </math> Эта формула позволяет свести вычисление <math display="inline">I_k</math> к <math display="inline">I_{k-1}</math>, затем к <math display="inline">I_{k-2}</math>, и так далее, пока не дойдем до <math display="inline">I_1 = \int \frac{dt}{t^2+a^2} = \frac{1}{a} \arctan(\frac{t}{a}) + C</math>.
'''Вывод рекуррентной формулы для <math display="inline">I_k = \int \frac{dt}{(t^2+a^2)^k}</math>:''' <math display="inline">I_k = \int \frac{dt}{(t^2+a^2)^k} = \frac{1}{a^2} \int \frac{a^2+t^2-t^2}{(t^2+a^2)^k} dt</math> <math display="inline">I_k = \frac{1}{a^2} \left( \int \frac{t^2+a^2}{(t^2+a^2)^k} dt - \int \frac{t^2}{(t^2+a^2)^k} dt \right)</math> <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} - \int t \cdot \frac{t}{(t^2+a^2)^k} dt \right)</math> Интегрируем по частям <math display="inline">\int t \cdot \frac{t dt}{(t^2+a^2)^k}</math>: <math display="inline">u=t \implies du=dt</math> <math display="inline">dv = \frac{t dt}{(t^2+a^2)^k} \implies v = \int \frac{t dt}{(t^2+a^2)^k} = \frac{1}{2} \int \frac{d(t^2+a^2)}{(t^2+a^2)^k} = \frac{1}{2(1-k)(t^2+a^2)^{k-1}}</math> <math display="inline">\int t dv = uv - \int v du = \frac{t}{2(1-k)(t^2+a^2)^{k-1}} - \int \frac{dt}{2(1-k)(t^2+a^2)^{k-1}}</math> <math display="inline">\int t \frac{t dt}{(t^2+a^2)^k} = \frac{t}{2(1-k)(t^2+a^2)^{k-1}} - \frac{1}{2(1-k)} I_{k-1}</math> Подставляем обратно в выражение для <math display="inline">I_k</math>: <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} - \left[ \frac{t}{2(1-k)(t^2+a^2)^{k-1}} - \frac{1}{2(1-k)} I_{k-1} \right] \right)</math> <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} \left(1 + \frac{1}{2(k-1)}\right) + \frac{t}{2(k-1)(t^2+a^2)^{k-1}} \right)</math> <math display="inline">I_k = \frac{1}{a^2} \left( I_{k-1} \frac{2k-2+1}{2(k-1)} + \frac{t}{2(k-1)(t^2+a^2)^{k-1}} \right)</math> <math display="block"> I_k = \frac{2k-3}{2a^2(k-1)} I_{k-1} + \frac{1}{2a^2(k-1)} \frac{t}{(t^2+a^2)^{k-1}} </math> Эта формула позволяет свести вычисление <math display="inline">I_k</math> к <math display="inline">I_{k-1}</math>, затем к <math display="inline">I_{k-2}</math>, и так далее, пока не дойдем до <math display="inline">I_1 = \int \frac{dt}{t^2+a^2} = \frac{1}{a} \arctan(\frac{t}{a}) + C</math>.


'''Вывод:''' Интеграл от любой рациональной функции <math display="inline">\int \frac{P_n(x)}{Q_m(x)} dx</math> всегда выражается через элементарные функции (многочлены, логарифмы, арктангенсы). ## Интегрирование некоторых иррациональных функций
'''Вывод:''' Интеграл от любой рациональной функции <math display="inline">\int \frac{P_n(x)}{Q_m(x)} dx</math> всегда выражается через элементарные функции (многочлены, логарифмы, арктангенсы).  
 
Здесь <math display="inline">R(u,v)</math> обозначает рациональную функцию своих аргументов, т.е. отношение двух многочленов от <math display="inline">u</math> и <math display="inline">v</math>.
 
<span id="интегралы-вида-int-rleftx-leftfracaxbcxdrightr_1-dots-leftfracaxbcxdrightr_nright-dx"></span>
=== 1. Интегралы вида <math display="inline">\int R\left(x, \left(\frac{ax+b}{cx+d}\right)^{r_1}, \dots, \left(\frac{ax+b}{cx+d}\right)^{r_n}\right) dx</math> ===
 
где <math display="inline">r_1, \dots, r_n \in \mathbb{Q}</math> - рациональные числа, <math display="inline">a, b, c, d \in \mathbb{R}</math>, <math display="inline">ad-bc \neq 0</math>. Пусть <math display="inline">r_i = p_i/q_i</math>. Находим <math display="inline">h = \text{НОК}(q_1, \dots, q_n)</math> - наименьшее общее кратное знаменателей <math display="inline">q_i</math>. Используем '''подстановку''': <math display="block"> t^h = \frac{ax+b}{cx+d} </math> Из этого равенства выражаем <math display="inline">x</math> как рациональную функцию от <math display="inline">t</math>. Тогда <math display="inline">dx</math> также будет выражаться как произведение рациональной функции от <math display="inline">t</math> на <math display="inline">dt</math>. Все дробные степени <math display="inline">(\frac{ax+b}{cx+d})^{p_i/q_i} = (t^h)^{p_i/q_i} = t^{h p_i / q_i}</math> станут целыми степенями <math display="inline">t</math>, так как <math display="inline">h</math> делится на каждый <math display="inline">q_i</math>. В результате интеграл сводится к интегралу от рациональной функции от <math display="inline">t</math>.
 
'''Пример:''' <math display="inline">\int \frac{1+3\sqrt{x}+\sqrt[3]{x}}{5-\sqrt[6]{x}+\sqrt[3]{x^2}} dx</math> Здесь <math display="inline">\frac{ax+b}{cx+d} = x</math>. Степени: <math display="inline">1/2, 1/3, 1/6, 2/3</math>. Знаменатели: <math display="inline">2, 3, 6, 3</math>. <math display="inline">h = \text{НОК}(2, 3, 6) = 6</math>. Подстановка: <math display="inline">x = t^6</math>. Тогда <math display="inline">dx = 6t^5 dt</math>. <math display="inline">\sqrt{x} = t^3</math>, <math display="inline">\sqrt[3]{x} = t^2</math>, <math display="inline">\sqrt[6]{x} = t</math>, <math display="inline">\sqrt[3]{x^2} = x^{2/3} = (t^6)^{2/3} = t^4</math>. <math display="inline">\int \frac{1+3t^3+t^2}{5-t+t^4} (6t^5 dt) = 6 \int \frac{t^5(1+3t^3+t^2)}{t^4-t+5} dt</math> Получили интеграл от рациональной функции.
 
<span id="интегралы-вида-int-rx-sqrtax2bxc-dx"></span>
=== 2. Интегралы вида <math display="inline">\int R(x, \sqrt{ax^2+bx+c}) dx</math> ===
 
где <math display="inline">a \neq 0</math>, <math display="inline">b^2-4ac \neq 0</math>. Такие интегралы берутся с помощью '''подстановок Эйлера'''.
 
# '''Первая подстановка Эйлера:''' Если <math display="inline">a > 0</math>, полагаем: <math display="inline">\sqrt{ax^2+bx+c} = \pm \sqrt{a} x + t</math> (выбираем один из знаков) Например, <math display="inline">\sqrt{ax^2+bx+c} = \sqrt{a} x + t</math>. Возводим в квадрат: <math display="inline">ax^2+bx+c = ax^2 + 2t\sqrt{a} x + t^2</math> <math display="inline">bx+c = 2t\sqrt{a} x + t^2</math> <math display="inline">x(b - 2t\sqrt{a}) = t^2 - c \implies x = \frac{t^2-c}{b-2t\sqrt{a}}</math> Отсюда <math display="inline">x</math> и <math display="inline">dx</math> выражаются рационально через <math display="inline">t</math>. <math display="inline">\sqrt{ax^2+bx+c}</math> также выражается рационально через <math display="inline">t</math> (из подстановки). Интеграл сводится к интегралу от рациональной функции от <math display="inline">t</math>.
# '''Вторая подстановка Эйлера:''' Если <math display="inline">c > 0</math>, полагаем: <math display="inline">\sqrt{ax^2+bx+c} = xt \pm \sqrt{c}</math> (выбираем один из знаков) Например, <math display="inline">\sqrt{ax^2+bx+c} = xt + \sqrt{c}</math>. Возводим в квадрат: <math display="inline">ax^2+bx+c = x^2t^2 + 2xt\sqrt{c} + c</math> <math display="inline">ax^2+bx = x^2t^2 + 2xt\sqrt{c}</math> Делим на <math display="inline">x</math> (при <math display="inline">x\neq 0</math>): <math display="inline">ax+b = xt^2 + 2t\sqrt{c}</math> <math display="inline">x(a-t^2) = 2t\sqrt{c} - b \implies x = \frac{2t\sqrt{c}-b}{a-t^2}</math> Снова <math display="inline">x, dx, \sqrt{ax^2+bx+c}</math> выражаются рационально через <math display="inline">t</math>.
# '''Третья подстановка Эйлера:''' Если квадратный трехчлен <math display="inline">ax^2+bx+c</math> имеет действительные корни <math display="inline">x_1, x_2</math> (т.е. <math display="inline">b^2-4ac > 0</math>), то <math display="inline">ax^2+bx+c = a(x-x_1)(x-x_2)</math>. Полагаем: <math display="inline">\sqrt{ax^2+bx+c} = t(x-x_1)</math> (или <math display="inline">t(x-x_2)</math>) <math display="inline">\sqrt{a(x-x_1)(x-x_2)} = t(x-x_1)</math> Возводим в квадрат: <math display="inline">a(x-x_1)(x-x_2) = t^2(x-x_1)^2</math> Делим на <math display="inline">(x-x_1)</math> (при <math display="inline">x \neq x_1</math>): <math display="inline">a(x-x_2) = t^2(x-x_1)</math> <math display="inline">ax - ax_2 = t^2 x - t^2 x_1</math> <math display="inline">x(a-t^2) = ax_2 - t^2 x_1 \implies x = \frac{ax_2 - t^2 x_1}{a-t^2}</math> Снова <math display="inline">x, dx, \sqrt{ax^2+bx+c}</math> выражаются рационально через <math display="inline">t</math>.
 
'''Пример:''' <math display="inline">\int \frac{dx}{x \sqrt{x^2-x+1}}</math> Здесь <math display="inline">a=1>0, c=1>0, b^2-4ac = (-1)^2 - 4(1)(1) = -3 < 0</math>. Можно использовать 1-ю или 2-ю подстановку Эйлера. Используем 2-ю: <math display="inline">\sqrt{x^2-x+1} = xt + \sqrt{c} = xt + 1</math>. <math display="inline">x^2-x+1 = (xt+1)^2 = x^2t^2 + 2xt + 1</math> <math display="inline">x^2-x = x^2t^2 + 2xt</math> Делим на <math display="inline">x</math> (<math display="inline">x \neq 0</math>): <math display="inline">x-1 = xt^2 + 2t</math> <math display="inline">x - xt^2 = 1+2t \implies x(1-t^2) = 1+2t \implies x = \frac{1+2t}{1-t^2}</math> <math display="inline">dx = \left( \frac{2(1-t^2) - (1+2t)(-2t)}{(1-t^2)^2} \right) dt = \frac{2-2t^2+2t+4t^2}{(1-t^2)^2} dt = \frac{2t^2+2t+2}{(1-t^2)^2} dt = \frac{2(t^2+t+1)}{(1-t^2)^2} dt</math> <math display="inline">\sqrt{x^2-x+1} = xt+1 = \frac{1+2t}{1-t^2} t + 1 = \frac{t+2t^2 + 1-t^2}{1-t^2} = \frac{t^2+t+1}{1-t^2}</math> Подставляем в интеграл: <math display="inline">\int \frac{1}{x \sqrt{x^2-x+1}} dx = \int \frac{1}{\frac{1+2t}{1-t^2} \cdot \frac{t^2+t+1}{1-t^2}} \cdot \frac{2(t^2+t+1)}{(1-t^2)^2} dt</math> <math display="inline">= \int \frac{(1-t^2)^2}{(1+2t)(t^2+t+1)} \cdot \frac{2(t^2+t+1)}{(1-t^2)^2} dt = \int \frac{2}{1+2t} dt</math> <math display="inline">= \ln|1+2t| + C</math> Возвращаемся к <math display="inline">x</math>. Из <math display="inline">\sqrt{x^2-x+1} = xt+1</math>, получаем <math display="inline">xt = \sqrt{x^2-x+1}-1</math>, т.е. <math display="inline">t = \frac{\sqrt{x^2-x+1}-1}{x}</math>. Итоговый ответ: <math display="inline">\ln\left|1 + 2 \frac{\sqrt{x^2-x+1}-1}{x}\right| + C = \ln\left|\frac{x+2\sqrt{x^2-x+1}-2}{x}\right| + C</math>.
 
''(Примечание: Пример в конспекте был <math display="inline">\int \frac{\sqrt{x^2-x+1}-1}{x\sqrt{x^2-x+1}} dx</math>. Это несколько другой интеграл, но подстановка та же самая.)''


[[Категория:МатАнПрод]]
[[Категория:МатАнПрод]]

Текущая версия от 16:35, 14 апреля 2025

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

Определение: Многочленом (полиномом) 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 всегда выражается через элементарные функции (многочлены, логарифмы, арктангенсы).