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

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

Материал из Мадока ВТ Вики
Нет описания правки
Нет описания правки
 
(не показаны 3 промежуточные версии этого же участника)
Строка 1: Строка 1:
== Основные понятия ==
'''Определение:''' Многочленом (полиномом) <math display="inline">P_n(x)</math> степени <math display="inline">n</math> (<math display="inline">n \ge 0</math>, <math display="inline">n \in \mathbb{Z}</math>) называется функция вида: <math display="inline">P_n(x) = a_0 + a_1 x + a_2 x^2 + \dots + a_n x^n</math>, где <math display="inline">a_i \in \mathbb{R}</math>, <math display="inline">a_n \neq 0</math>.
'''Определение:''' Многочленом (полиномом) <math display="inline">P_n(x)</math> степени <math display="inline">n</math> (<math display="inline">n \ge 0</math>, <math display="inline">n \in \mathbb{Z}</math>) называется функция вида: <math display="inline">P_n(x) = a_0 + a_1 x + a_2 x^2 + \dots + a_n x^n</math>, где <math display="inline">a_i \in \mathbb{R}</math>, <math display="inline">a_n \neq 0</math>.


Строка 29: Строка 30:


'''Теорема:''' Любой многочлен <math display="inline">Q_m(x)</math> степени <math display="inline">m \ge 1</math> с действительными коэффициентами может быть представлен в виде произведения своего старшего коэффициента <math display="inline">c_m</math> и множителей вида <math display="inline">(x-x_i)</math> и <math display="inline">(x^2+p_j x+q_j)</math>, где <math display="inline">x_i</math> - действительные корни, а <math display="inline">x^2+p_j x+q_j</math> - неприводимые над <math display="inline">\mathbb{R}</math> квадратные трехчлены (<math display="inline">p_j^2 - 4q_j < 0</math>), соответствующие парам комплексно-сопряженных корней. Будем считать <math display="inline">c_m=1</math> (если нет, можно вынести за скобки). <math display="inline">Q_m(x) = (x-x_1)^{k_1} (x-x_2)^{k_2} \dots (x-x_p)^{k_p} (x^2+p_1x+q_1)^{l_1} (x^2+p_2x+q_2)^{l_2} \dots (x^2+p_s x+q_s)^{l_s}</math> где: * <math display="inline">x_i \in \mathbb{R}</math> - различные действительные корни * <math display="inline">k_i \in \mathbb{N}</math> - кратности действительных корней * <math display="inline">p_j, q_j \in \mathbb{R}</math>, <math display="inline">p_j^2 - 4q_j < 0</math> (квадратные трехчлены неприводимы) * <math display="inline">l_j \in \mathbb{N}</math> - кратности пар комплексно-сопряженных корней * Сумма кратностей: <math display="inline">\sum_{i=1}^{p} k_i + 2 \sum_{j=1}^{s} l_j = m = \deg(Q_m(x))</math>.
'''Теорема:''' Любой многочлен <math display="inline">Q_m(x)</math> степени <math display="inline">m \ge 1</math> с действительными коэффициентами может быть представлен в виде произведения своего старшего коэффициента <math display="inline">c_m</math> и множителей вида <math display="inline">(x-x_i)</math> и <math display="inline">(x^2+p_j x+q_j)</math>, где <math display="inline">x_i</math> - действительные корни, а <math display="inline">x^2+p_j x+q_j</math> - неприводимые над <math display="inline">\mathbb{R}</math> квадратные трехчлены (<math display="inline">p_j^2 - 4q_j < 0</math>), соответствующие парам комплексно-сопряженных корней. Будем считать <math display="inline">c_m=1</math> (если нет, можно вынести за скобки). <math display="inline">Q_m(x) = (x-x_1)^{k_1} (x-x_2)^{k_2} \dots (x-x_p)^{k_p} (x^2+p_1x+q_1)^{l_1} (x^2+p_2x+q_2)^{l_2} \dots (x^2+p_s x+q_s)^{l_s}</math> где: * <math display="inline">x_i \in \mathbb{R}</math> - различные действительные корни * <math display="inline">k_i \in \mathbb{N}</math> - кратности действительных корней * <math display="inline">p_j, q_j \in \mathbb{R}</math>, <math display="inline">p_j^2 - 4q_j < 0</math> (квадратные трехчлены неприводимы) * <math display="inline">l_j \in \mathbb{N}</math> - кратности пар комплексно-сопряженных корней * Сумма кратностей: <math display="inline">\sum_{i=1}^{p} k_i + 2 \sum_{j=1}^{s} l_j = m = \deg(Q_m(x))</math>.
== Разложение правильной рациональной дроби на простейшие дроби ==
'''Лемма 1:''' Пусть <math display="inline">\frac{P_n(x)}{Q_m(x)}</math> - правильная рациональная дробь (<math display="inline">n < m</math>). Пусть <math display="inline">a \in \mathbb{R}</math> является корнем знаменателя <math display="inline">Q_m(x)</math> кратности <math display="inline">k \ge 1</math>, т.е. <math display="inline">Q_m(a)=0</math>, и <math display="inline">Q_m(x) = (x-a)^k Q_{m-k}^{(1)}(x)</math>, где <math display="inline">Q_{m-k}^{(1)}(a) \neq 0</math>. Тогда существует единственное представление: <math display="inline">\frac{P_n(x)}{Q_m(x)} = \frac{A_k}{(x-a)^k} + \frac{\tilde{P}_1(x)}{(x-a)^{k-1} Q_{m-k}^{(1)}(x)}</math> где <math display="inline">A_k \in \mathbb{R}</math> - константа, а <math display="inline">\frac{\tilde{P}_1(x)}{(x-a)^{k-1} Q_{m-k}^{(1)}(x)}</math> - также правильная рациональная дробь.
'''Нахождение коэффициента <math display="inline">A_k</math>:''' Умножим обе части равенства на <math display="inline">(x-a)^k</math>: <math display="inline">\frac{P_n(x)}{Q_{m-k}^{(1)}(x)} = A_k + \frac{\tilde{P}_1(x) (x-a)}{Q_{m-k}^{(1)}(x)}</math> Положим <math display="inline">x=a</math>. Так как <math display="inline">Q_{m-k}^{(1)}(a) \neq 0</math>, получим: <math display="inline">\frac{P_n(a)}{Q_{m-k}^{(1)}(a)} = A_k + 0</math> '''(!) Формула для <math display="inline">A_k</math>:''' <math display="inline">A_k = \frac{P_n(a)}{Q_{m-k}^{(1)}(a)}</math>
'''Доказательство существования и единственности (кратко):''' Из <math display="inline">\frac{P_n(x)}{Q_m(x)} - \frac{A_k}{(x-a)^k} = \frac{P_n(x) - A_k Q_{m-k}^{(1)}(x)}{(x-a)^k Q_{m-k}^{(1)}(x)}</math> Выберем <math display="inline">A_k = \frac{P_n(a)}{Q_{m-k}^{(1)}(a)}</math>. Тогда числитель <math display="inline">R(x) = P_n(x) - A_k Q_{m-k}^{(1)}(x)</math> имеет корень при <math display="inline">x=a</math>, так как <math display="inline">R(a) = P_n(a) - \frac{P_n(a)}{Q_{m-k}^{(1)}(a)} Q_{m-k}^{(1)}(a) = 0</math>. Следовательно, <math display="inline">R(x) = (x-a) \tilde{P}_1(x)</math> по теореме Безу. Тогда: <math display="inline">\frac{(x-a) \tilde{P}_1(x)}{(x-a)^k Q_{m-k}^{(1)}(x)} = \frac{\tilde{P}_1(x)}{(x-a)^{k-1} Q_{m-k}^{(1)}(x)}</math> Можно показать, что полученная дробь является правильной. Единственность <math display="inline">A_k</math> доказывается от противного: предположим два разложения с <math display="inline">A_{k,1}</math> и <math display="inline">A_{k,2}</math>, приведем к общему знаменателю и подставим <math display="inline">x=a</math>, что даст <math display="inline">A_{k,1} = A_{k,2}</math>. * Равенство <math display="inline">A_1 Q^{(1)}(x) + \tilde{P}_1(x)(x-a) = A_2 Q^{(1)}(x) + \tilde{P}_2(x)(x-a)</math> при <math display="inline">x=a</math> дает <math display="inline">A_1 Q^{(1)}(a) = A_2 Q^{(1)}(a)</math>, откуда <math display="inline">A_1=A_2</math>, так как <math display="inline">Q^{(1)}(a) \neq 0</math>.
'''Следствие (Применение Леммы 1 k-раз):''' Повторяя эту процедуру <math display="inline">k</math> раз для корня <math display="inline">a</math>, мы выделим все слагаемые, соответствующие этому корню: <math display="inline">\frac{P_n(x)}{Q_m(x)} = \frac{A_k}{(x-a)^k} + \frac{A_{k-1}}{(x-a)^{k-1}} + \dots + \frac{A_1}{x-a} + \frac{\hat{P}(x)}{Q_{m-k}^{(1)}(x)}</math> где <math display="inline">A_i \in \mathbb{R}</math>, а <math display="inline">\frac{\hat{P}(x)}{Q_{m-k}^{(1)}(x)}</math> - правильная рациональная дробь, знаменатель которой уже не имеет корня <math display="inline">a</math> (<math display="inline">Q_{m-k}^{(1)}(a) \neq 0</math>).
'''Лемма 2:''' Пусть <math display="inline">\frac{P_n(x)}{Q_m(x)}</math> - правильная рациональная дробь (<math display="inline">n < m</math>). Пусть <math display="inline">(x^2+px+q)</math> - неприводимый над <math display="inline">\mathbb{R}</math> множитель (<math display="inline">p^2-4q < 0</math>) знаменателя <math display="inline">Q_m(x)</math> кратности <math display="inline">k \ge 1</math>. То есть <math display="inline">Q_m(x) = (x^2+px+q)^k \tilde{Q}_{m-2k}(x)</math>, где <math display="inline">\tilde{Q}_{m-2k}(x_0) \neq 0</math> для корней <math display="inline">x_0 = \alpha \pm i\beta</math> уравнения <math display="inline">x^2+px+q=0</math>. Тогда существует единственное представление: <math display="inline">\frac{P_n(x)}{Q_m(x)} = \frac{A_k x + B_k}{(x^2+px+q)^k} + \frac{\tilde{P}_1(x)}{(x^2+px+q)^{k-1} \tilde{Q}_{m-2k}(x)}</math> ('''**''') где <math display="inline">A_k, B_k \in \mathbb{R}</math> - константы, а вторая дробь - правильная рациональная дробь.
'''Нахождение коэффициентов <math display="inline">A_k, B_k</math>:''' Умножим обе части равенства на <math display="inline">(x^2+px+q)^k</math>: <math display="inline">\frac{P_n(x)}{\tilde{Q}_{m-2k}(x)} = A_k x + B_k + \frac{\tilde{P}_1(x) (x^2+px+q)}{\tilde{Q}_{m-2k}(x)}</math> Пусть <math display="inline">x_0 = \alpha+i\beta</math> - один из комплексных корней уравнения <math display="inline">x^2+px+q=0</math>. Подставим <math display="inline">x=x_0</math>: <math display="inline">\frac{P_n(x_0)}{\tilde{Q}_{m-2k}(x_0)} = A_k x_0 + B_k + 0</math> Обозначим комплексное число <math display="inline">R = \frac{P_n(x_0)}{\tilde{Q}_{m-2k}(x_0)} = \text{Re}(R) + i \text{Im}(R)</math>. <math display="inline">A_k (\alpha+i\beta) + B_k = \text{Re}(R) + i \text{Im}(R)</math> <math display="inline">(A_k \alpha + B_k) + i (A_k \beta) = \text{Re}(R) + i \text{Im}(R)</math> Приравнивая действительные и мнимые части, получаем систему для <math display="inline">A_k, B_k</math>: <math display="inline">\begin{cases} A_k \alpha + B_k = \text{Re}(R) \\ A_k \beta = \text{Im}(R) \end{cases}</math> Так как <math display="inline">\beta \neq 0</math>, система имеет единственное решение: <math display="inline">A_k = \frac{\text{Im}(R)}{\beta}</math> <math display="inline">B_k = \text{Re}(R) - A_k \alpha = \text{Re}(R) - \alpha \frac{\text{Im}(R)}{\beta}</math> Так как <math display="inline">P_n, \tilde{Q}_{m-2k}</math> имеют действительные коэффициенты, и если <math display="inline">x_0</math> - корень, то <math display="inline">\bar{x}_0</math> - тоже, можно показать, что <math display="inline">A_k, B_k</math> всегда получаются действительными числами.
'''Доказательство существования и единственности (кратко):''' Аналогично случаю действительного корня. Рассмотрим разность <math display="inline">P_n(x) - (A_k x + B_k) \tilde{Q}_{m-2k}(x)</math>. Подстановка <math display="inline">x=x_0</math> дает <math display="inline">P_n(x_0) - (A_k x_0 + B_k) \tilde{Q}_{m-2k}(x_0) = P_n(x_0) - \frac{P_n(x_0)}{\tilde{Q}_{m-2k}(x_0)} \tilde{Q}_{m-2k}(x_0) = 0</math>. Значит <math display="inline">x_0</math> - корень этого полинома. Так как полином имеет действительные коэффициенты, то <math display="inline">\bar{x}_0</math> тоже корень. Следовательно, полином делится на <math display="inline">(x-x_0)(x-\bar{x}_0) = x^2+px+q</math>. <math display="inline">P_n(x) - (A_k x + B_k) \tilde{Q}_{m-2k}(x) = (x^2+px+q) \hat{P}(x)</math>. Это позволяет переписать исходную дробь в виде ('''**'''). Можно показать, что остаточная дробь правильная. Единственность доказывается аналогично случаю действительных корней, приравнивая два разложения и подставляя <math display="inline">x=x_0</math>. * Равенство <math display="inline">(A_1x+B_1)\tilde{Q}(x) + \tilde{P}_1(x)(x^2+px+q) = (A_2x+B_2)\tilde{Q}(x) + \tilde{P}_2(x)(x^2+px+q)</math> при <math display="inline">x=x_0</math> дает <math display="inline">(A_1x_0+B_1)\tilde{Q}(x_0) = (A_2x_0+B_2)\tilde{Q}(x_0)</math>. Так как <math display="inline">\tilde{Q}(x_0) \neq 0</math>, имеем <math display="inline">A_1x_0+B_1=A_2x_0+B_2</math>. Приравнивая действительные и мнимые части, получаем <math display="inline">A_1=A_2</math> и <math display="inline">B_1=B_2</math>.
'''Следствие (Применение Леммы 2 k-раз):''' Повторяя эту процедуру <math display="inline">k</math> раз для множителя <math display="inline">(x^2+px+q)</math>, мы выделим все слагаемые, соответствующие этой паре корней: <math display="inline">\frac{P_n(x)}{Q_m(x)} = \frac{A_k x + B_k}{(x^2+px+q)^k} + \frac{A_{k-1} x + B_{k-1}}{(x^2+px+q)^{k-1}} + \dots + \frac{A_1 x + B_1}{x^2+px+q} + \frac{\hat{P}(x)}{\tilde{Q}_{m-2k}(x)}</math> где <math display="inline">A_i, B_i \in \mathbb{R}</math>, а <math display="inline">\frac{\hat{P}(x)}{\tilde{Q}_{m-2k}(x)}</math> - правильная рациональная дробь, знаменатель которой уже не имеет корней <math display="inline">x_0, \bar{x}_0</math>.
<span id="общая-теорема-о-разложении-на-простейшие-дроби"></span>
=== Общая теорема о разложении на простейшие дроби ===
'''Теорема:''' Любую правильную рациональную дробь <math display="inline">\frac{P_n(x)}{Q_m(x)}</math>, знаменатель которой разложен на множители над <math display="inline">\mathbb{R}</math> вида <math display="inline">Q_m(x) = c_m (x-x_1)^{k_1} \dots (x-x_p)^{k_p} (x^2+p_1x+q_1)^{l_1} \dots (x^2+p_s x+q_s)^{l_s}</math> можно единственным образом представить в виде суммы простейших дробей: <math display="block">
\frac{P_n(x)}{Q_m(x)} =
  \sum_{i=1}^{p} \left( \frac{A_{i,k_i}}{(x-x_i)^{k_i}} + \frac{A_{i,k_i-1}}{(x-x_i)^{k_i-1}} + \dots + \frac{A_{i,1}}{x-x_i} \right) + \sum_{j=1}^{s} \left( \frac{B_{j,l_j} x + C_{j,l_j}}{(x^2+p_j x+q_j)^{l_j}} + \frac{B_{j,l_j-1} x + C_{j,l_j-1}}{(x^2+p_j x+q_j)^{l_j-1}} + \dots + \frac{B_{j,1} x + C_{j,1}}{x^2+p_j x+q_j} \right)
</math> где <math display="inline">A_{i,r}, B_{j,t}, C_{j,t}</math> - действительные коэффициенты. ''(Примечание: В конспекте коэффициенты для квадратичных множителей обозначались <math display="inline">A_j x + B_j</math>, здесь использовано <math display="inline">B_j x + C_j</math> для избежания конфликта с <math display="inline">A_{i,r}</math> для линейных множителей, что является стандартной практикой).''
Интегрирование рациональной функции сводится к интегрированию многочлена (если дробь неправильная) и суммы простейших дробей.
ivabus@celerrime-x ➜ ~ $ pandoc -f markdown -t mediawiki /Users/ivabus/Downloads/Матан\ №4.md  -o -
<span id="теорема-о-разложении-правильной-рациональной-дроби"></span>
=== Теорема о разложении правильной рациональной дроби ===
'''Теорема:''' Пусть <math display="inline">\frac{P_n(x)}{Q_m(x)}</math> - правильная рациональная дробь (<math display="inline">n < m</math>), и знаменатель <math display="inline">Q_m(x)</math> (со старшим коэффициентом, равным 1) разложен на множители над <math display="inline">\mathbb{R}</math>: <math display="inline">Q_m(x) = (x-x_1)^{k_1} (x-x_2)^{k_2} \dots (x-x_p)^{k_p} (x^2+p_1x+q_1)^{l_1} (x^2+p_2x+q_2)^{l_2} \dots (x^2+p_s x+q_s)^{l_s}</math> где <math display="inline">x_i \in \mathbb{R}</math>, <math display="inline">k_i \in \mathbb{N}</math>, <math display="inline">p_j, q_j \in \mathbb{R}</math>, <math display="inline">p_j^2 - 4q_j < 0</math>, <math display="inline">l_j \in \mathbb{N}</math>. Тогда существует '''единственное''' разложение дроби <math display="inline">\frac{P_n(x)}{Q_m(x)}</math> в сумму простейших дробей вида: <math display="block">
\frac{P_n(x)}{Q_m(x)} = \sum_{i=1}^{p} \sum_{r=1}^{k_i} \frac{A_{i,r}}{(x-x_i)^r} + \sum_{j=1}^{s} \sum_{t=1}^{l_j} \frac{B_{j,t} x + C_{j,t}}{(x^2+p_j x+q_j)^t}
</math> где <math display="inline">A_{i,r}, B_{j,t}, C_{j,t} \in \mathbb{R}</math> - некоторые константы.
''(Примечание: В оригинальном конспекте использовалась несколько иная индексация сумм: <math display="inline">\sum_{i=1}^{p} \sum_{j=1}^{k_i} \frac{A_{i,j}}{(x-x_i)^{k_i-j+1}} + \sum_{i=1}^{s} \sum_{j=1}^{l_i} \frac{B_{i,j}x+C_{i,j}}{(x^2+p_i x+q_i)^{l_i-j+1}}</math>. Представленная выше форма является стандартной и эквивалентной.)''
'''Замечание:''' Если исходная дробь <math display="inline">\frac{P_n(x)}{Q_m(x)}</math> неправильная (<math display="inline">n \ge m</math>), то сначала необходимо выделить целую часть (многочлен) с помощью деления “уголком”: <math display="inline">\frac{P_n(x)}{Q_m(x)} = M_{n-m}(x) + \frac{N_k(x)}{Q_m(x)}</math>, где <math display="inline">k < m</math>. Затем правильную дробь <math display="inline">\frac{N_k(x)}{Q_m(x)}</math> раскладывают на простейшие.
<span id="интегрирование-простейших-дробей"></span>
=== Интегрирование простейших дробей ===
'''Определение:''' Простейшими дробями I, II, III, IV типов называются соответственно дроби вида: I. <math display="inline">\frac{A}{x-a}</math> II. <math display="inline">\frac{A}{(x-a)^k}</math>, где <math display="inline">k \in \mathbb{N}, k \ge 2</math>. III. <math display="inline">\frac{Ax+B}{x^2+px+q}</math>, где <math display="inline">p^2-4q < 0</math>. IV. <math display="inline">\frac{Ax+B}{(x^2+px+q)^k}</math>, где <math display="inline">p^2-4q < 0</math>, <math display="inline">k \in \mathbb{N}, k \ge 2</math>.
Интегралы от простейших дробей:
I. <math display="inline">\int \frac{A}{x-a} dx = A \int \frac{d(x-a)}{x-a} = A \ln|x-a| + C</math>
<ol start="2" style="list-style-type: upper-roman;">
<li><p><math display="inline">\int \frac{A}{(x-a)^k} dx = A \int (x-a)^{-k} d(x-a) = A \frac{(x-a)^{-k+1}}{-k+1} + C = \frac{A}{(1-k)(x-a)^{k-1}} + C</math></p></li>
<li><p><math display="inline">\int \frac{Ax+B}{x^2+px+q} dx</math> Знаменатель <math display="inline">x^2+px+q</math> имеет комплексные корни, т.к. <math display="inline">D=p^2-4q<0</math>. Выделим в числителе производную знаменателя <math display="inline">(x^2+px+q)' = 2x+p</math>: <math display="inline">Ax+B = \frac{A}{2}(2x+p) - \frac{Ap}{2} + B = \frac{A}{2}(2x+p) + \left(B - \frac{Ap}{2}\right)</math> Тогда интеграл разбивается на два: <math display="inline">\int \frac{\frac{A}{2}(2x+p)}{x^2+px+q} dx + \int \frac{B - \frac{Ap}{2}}{x^2+px+q} dx</math> Первый интеграл: <math display="inline">\frac{A}{2} \int \frac{d(x^2+px+q)}{x^2+px+q} = \frac{A}{2} \ln(x^2+px+q) + C_1</math> (модуль не нужен, т.к. <math display="inline">x^2+px+q > 0</math> при <math display="inline">D<0</math>). Второй интеграл: Выделим полный квадрат в знаменателе: <math display="inline">x^2+px+q = \left(x + \frac{p}{2}\right)^2 + q - \frac{p^2}{4} = \left(x + \frac{p}{2}\right)^2 + a^2</math>, где <math display="inline">a^2 = q - \frac{p^2}{4} > 0</math>. <math display="inline">\left(B - \frac{Ap}{2}\right) \int \frac{dx}{(x+p/2)^2 + a^2} = \left(B - \frac{Ap}{2}\right) \frac{1}{a} \arctan\left(\frac{x+p/2}{a}\right) + C_2</math> Объединяя: <math display="inline">\int \frac{Ax+B}{x^2+px+q} dx = \frac{A}{2} \ln(x^2+px+q) + \frac{B - Ap/2}{\sqrt{q-p^2/4}} \arctan\left(\frac{x+p/2}{\sqrt{q-p^2/4}}\right) + C</math></p></li>
<li><p><math display="inline">\int \frac{Ax+B}{(x^2+px+q)^k} dx</math>, где <math display="inline">k \ge 2</math>. Аналогично типу III, выделяем производную знаменателя в числителе: <math display="inline">\int \frac{\frac{A}{2}(2x+p)}{(x^2+px+q)^k} dx + \int \frac{B - \frac{Ap}{2}}{(x^2+px+q)^k} dx</math> Первый интеграл: <math display="inline">\frac{A}{2} \int (x^2+px+q)^{-k} d(x^2+px+q) = \frac{A}{2} \frac{(x^2+px+q)^{-k+1}}{-k+1} + C_1 = \frac{A}{2(1-k)(x^2+px+q)^{k-1}} + C_1</math>. Второй интеграл: <math display="inline">\left(B - \frac{Ap}{2}\right) \int \frac{dx}{((x+p/2)^2 + a^2)^k} = \left(B - \frac{Ap}{2}\right) \int \frac{dt}{(t^2+a^2)^k}</math>, где <math display="inline">t=x+p/2</math>, <math display="inline">a^2=q-p^2/4</math>. Обозначим <math display="inline">I_k = \int \frac{dt}{(t^2+a^2)^k}</math>. Этот интеграл вычисляется с помощью рекуррентной формулы.</p></li></ol>
'''Вывод рекуррентной формулы для <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> всегда выражается через элементарные функции (многочлены, логарифмы, арктангенсы).


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

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