Рациональные дроби как суммы простых слагаемых
В алгебре разложение в частичные дроби или разложение в частичные дроби рациональной дроби (то есть дроби , у которой числитель и знаменатель являются полиномами ) — это операция, заключающаяся в выражении дроби в виде суммы многочлена (возможно, нулевого). ) и одну или несколько дробей с более простым знаменателем. [1]
Важность разложения на частные дроби заключается в том, что оно предоставляет алгоритмы для различных вычислений с рациональными функциями , включая явное вычисление первообразных , [2] разложение в ряд Тейлора , обратные Z-преобразования и обратные преобразования Лапласа . Эта концепция была открыта независимо в 1702 году Иоганном Бернулли и Готфридом Лейбницем . [3]
В символах разложение рациональной дроби в частные дроби вида, где f и g - многочлены, представляет собой ее выражение как![{\textstyle {\frac {f(x)}{g(x)}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {f(x)}{g(x)}}=p(x)+\sum _{j}{\frac {f_{j}(x)}{g_{j}(x )}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
где p ( x ) является многочленом, и для каждого j знаменатель g j ( x ) является степенью неприводимого многочлена ( который не разлагается на многочлены положительных степеней), а числитель f j ( x ) равен многочлен меньшей степени, чем степень этого неприводимого многочлена.
Когда используются явные вычисления, часто предпочтительнее более грубое разложение, которое заключается в замене «неприводимого многочлена» на « полином без квадратов » в описании результата. Это позволяет заменить полиномиальную факторизацию гораздо более простой в вычислении факторизацией без квадратов . Этого достаточно для большинства приложений и позволяет избежать введения иррациональных коэффициентов , когда коэффициенты входных полиномов являются целыми или рациональными числами .
Основные принципы
Позволять
![{\displaystyle R(x)={\frac {F}{G}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
рациональная дробьFGодномерные многочленынеопределенного xПолиномиальная часть
Существуют два полинома E и F 1 такие, что
![{\displaystyle {\frac {F}{G}}=E+{\frac {F_{1}}{G}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg F_{1}<\deg G,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
степеньP![{\displaystyle \градус P}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Это следует непосредственно из евклидова деления F на G , которое утверждает существование E и F 1 таких, что и![{\displaystyle F=EG+F_{1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg F_{1}<\deg G.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Это позволяет предположить на следующих шагах, что![{\displaystyle \deg F <\deg G.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Факторы знаменателя
Если и ![{\displaystyle \deg F <\deg G,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle G=G_{1}G_{2},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
G1G2взаимно простые многочлены , точто![{\displaystyle F_{1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle F_{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {F}{G}}={\frac {F_{1}}{G_{1}}}+{\frac {F_{2}}{G_{2}}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg F_{1}<\deg G_{1}\quad {\text{and}}\quad \deg F_{2}<\deg G_{2}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Это можно доказать следующим образом. Тождество Безу утверждает существование полиномов C и D таких, что
![{\displaystyle CG_{1}+DG_{2}=1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
( по1делитель G1G2.Пусть с — евклидово деление DF . По Установке получаем![{\displaystyle DF=G_{1}Q+F_{1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg F_{1}<\deg G_{1}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle G_{1}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle F_{2}=CF+QG_{2},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{aligned}{\frac {F}{G}}&={\frac {F(CG_{1}+DG_{2})}{G_{1}G_{2}}}= {\frac {DF}{G_{1}}}+{\frac {CF}{G_{2}}}\\&={\frac {F_{1}+G_{1}Q}{G_{1 }}}+{\frac {F_{2}-G_{2}Q}{G_{2}}}\\&={\frac {F_{1}}{G_{1}}}+{\frac {F_{2}}{G_{2}}}.\end{aligned}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg F_{2}<\deg G_{2}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle F=F_{2}G_{1}+F_{1}G_{2},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{aligned}\deg F_{2} &=\deg(FF_{1}G_{2})-\deg G_{1}\leq \max(\deg F,\deg( F_{1}G_{2}))-\deg G_{1}\\&<\max(\deg G,\deg(G_{1}G_{2}))-\deg G_{1}=\ град G_{2}\end{aligned}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Полномочия в знаменателе
Индуктивно используя предыдущее разложение, можно получить дроби вида где G — неприводимый многочлен . Если k > 1 , можно выполнить дальнейшее разложение, используя тот факт, что неприводимый многочлен является многочленом без квадратов , то есть является наибольшим общим делителем многочлена и его производной . Если является производной G , тождество Безу дает полиномы C и D такие, что и, таким образом , евклидово деление на дает полиномы и такие, что и Установка дает![{\displaystyle {\frac {F}{G^{k}}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg F<\deg G^{k}=k\deg G,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle 1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle G'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle CG+DG'=1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle F=FCG+FDG'.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle ФДГ'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle G}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle H_{k}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle Q}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle FDG'=QG+H_{k}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg H_ {k}<\deg G.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle F_{k-1}=FC+Q,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {F}{G^{k}}}={\frac {H_{k}}{G^{k}}}+{\frac {F_{k-1}}{G^ {k-1}}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg H_ {k}<\deg G.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Повторение этого процесса с вместо приводит в конечном итоге к следующей теореме.![{\displaystyle {\frac {F_{k-1}}{G^{k-1}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {F}{G^{k}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Заявление
Теорема . Пусть f и g — ненулевые полиномы над полем K. Запишите g как произведение степеней различных неприводимых многочленов:
![{\displaystyle g=\prod _{i=1}^{k}p_{i}^{n_{i}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Существуют (уникальные) полиномы b и a ij с deg a ij < deg p i такие, что
![{\displaystyle {\frac {f}{g}}=b+\sum _{i=1}^{k}\sum _{j=1}^{n_{i}}{\frac {a_{ij} }{p_{i}^{j}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Если deg f < deg g , то b = 0 .
Единственность можно доказать следующим образом. Пусть d = max(1 + deg f , deg g ) . В совокупности b и aij имеют d коэффициентов . Форма разложения определяет линейную карту векторов коэффициентов в полиномы f степени меньше d . Доказательство существования означает, что это отображение сюръективно . Поскольку два векторных пространства имеют одинаковую размерность, отображение также инъективно , что означает уникальность разложения. Кстати, это доказательство порождает алгоритм вычисления разложения с помощью линейной алгебры .
Если K — поле комплексных чисел , из фундаментальной теоремы алгебры следует, что все p i имеют степень один, а все числители — константы. Когда K является полем действительных чисел , некоторые из pi могут быть квадратичными, поэтому при разложении частных дробей также могут встречаться факторы линейных многочленов по степеням квадратичных многочленов.![{\displaystyle a_{ij}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
В предыдущей теореме можно заменить «различные неприводимые многочлены» на « попарно взаимно простые многочлены, взаимно простые со своей производной». Например, pi могут быть факторами бесквадратной факторизации g . Когда K является полем рациональных чисел , как это обычно бывает в компьютерной алгебре , это позволяет заменить факторизацию вычислением наибольшего общего делителя для вычисления разложения на частичные дроби.
Приложение к символическому интегрированию
В целях символического интегрирования предыдущий результат можно уточнить до
Теорема . Пусть f и g — ненулевые полиномы над полем K. Запишите g как произведение степеней попарно взаимно простых многочленов, не имеющих кратного корня в алгебраически замкнутом поле:
![{\displaystyle g=\prod _{i=1}^{k}p_{i}^{n_{i}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Существуют (уникальные) полиномы b и cij такие, что deg cij < deg p i , что
![{\displaystyle {\frac {f}{g}}=b+\sum _{i=1}^{k}\sum _{j=2}^{n_{i}}\left({\frac {c_ {ij}}{p_{i}^{j-1}}}\right)'+\sum _{i=1}^{k}{\frac {c_{i1}}{p_{i}}} .}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
где обозначает производную
![{\displaystyle X'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle X.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Это сводит вычисление первообразной рациональной функции к интегрированию последней суммы, которая называется логарифмической частью , поскольку ее первообразная представляет собой линейную комбинацию логарифмов.
Существуют различные методы вычисления разложения в теореме. Один простой способ называется методом Эрмита . Во-первых, b немедленно вычисляется путем евклидова деления f на g , что сводится к случаю, когда deg( f ) < deg( g ). Далее, известно deg( cij ) < deg( pi ) , поэтому каждый cij можно записать как многочлен с неизвестными коэффициентами. Приведя сумму дробей в теореме к общему знаменателю и приравняв коэффициенты при каждой степени x в двух числителях, получим систему линейных уравнений , решая которую можно получить искомые (единственные) значения неизвестных коэффициентов. .
Процедура
Учитывая два полинома и , где α n — различные константы и deg P < n , явные выражения для простейших дробей можно получить, предположив, что![{\displaystyle P (х)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle Q(x)=(x-\alpha _{1})(x-\alpha _{2})\cdots (x-\alpha _{n})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {P(x)}{Q(x)}}={\frac {c_{1}}{x-\alpha _{1}}}+{\frac {c_{2}} {x-\alpha _{2}}}+\cdots +{\frac {c_{n}}{x-\alpha _{n}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
c iприравнивания коэффициентовxнеопределенных коэффициентовx.Более прямое вычисление, которое сильно связано с интерполяцией Лагранжа , состоит в записи
![{\displaystyle {\frac {P(x)}{Q(x)}}=\sum _{i=1}^{n}{\frac {P(\alpha _{i})}{Q'( \alpha _{i})}}{\frac {1}{(x-\alpha _{i})}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
остатками/g![{\displaystyle Q'}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle Q}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\tfrac {1}{x-\alpha _{j}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Этот подход не учитывает ряд других случаев, но может быть соответствующим образом модифицирован:
- Если тогда необходимо выполнить евклидово деление P на Q , используя полиномиальное деление в длину , давая P ( x ) = E ( x ) Q ( x ) + R ( x ) с deg R < n . Деление на Q ( x ) дает
![{\ displaystyle \ градус P \ geq \ градус Q,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {P(x)}{Q(x)}}=E(x)+{\frac {R(x)}{Q(x)}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
а затем найдите простейшие дроби для оставшейся дроби (которая по определению удовлетворяет условию deg R < deg Q ). - Если Q ( x ) содержит множители, неприводимые по данному полю, то числитель N ( x ) каждой дроби с таким множителем F ( x ) в знаменателе нужно искать в виде многочлена с deg N < deg F , а не как константа. Например, возьмем следующее разложение по R :
![{\displaystyle {\frac {x^{2}+1}{(x+2)(x-1)\color {Blue}(x^{2}+x+1)}}={\frac {a }{x+2}}+{\frac {b}{x-1}}+{\frac {\color {OliveGreen}cx+d}{\color {Blue}x^{2}+x+1} }.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
- Предположим, Q ( x ) = ( x − α ) r S ( x ) и S ( α ) ≠ 0 , то есть α является корнем Q ( x ) кратности r . При разложении простейших дробей r первых степеней ( x − α ) будут выступать в качестве знаменателей простейших дробей (возможно, с нулевым числителем). Например, если S ( x ) = 1, разложение на простейшие дроби имеет вид
![{\displaystyle {\frac {P(x)}{Q(x)}}={\frac {P(x)}{(x-\alpha )^{r}}}={\frac {c_{1 }}{x-\alpha }}+{\frac {c_{2}}{(x-\alpha )^{2}}}+\cdots +{\frac {c_{r}}{(x-\ альфа )^{r}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Иллюстрация
В примере применения этой процедуры (3 x + 5)/(1 − 2 x ) 2 можно разложить в виде
![{\displaystyle {\frac {3x+5}{(1-2x)^{2}}}={\frac {A}{(1-2x)^{2}}}+{\frac {B}{ (1-2 раза)}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Очистка знаменателей показывает, что 3 x + 5 = A + B (1 − 2 x ) . Разложение и приравнивание коэффициентов при степенях x дает
5 = A + B и 3 x = −2 Bx
Решение этой системы линейных уравнений для A и B дает A = 13/2 и B = −3/2 . Следовательно,
![{\displaystyle {\frac {3x+5}{(1-2x)^{2}}}={\frac {13/2}{(1-2x)^{2}}}+{\frac {- 3/2}{(1-2x)}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Остаточный метод
Что касается комплексных чисел, предположим, что f ( x ) — рациональная правильная дробь, которую можно разложить на
![{\displaystyle f(x)=\sum _{i}\left({\frac {a_{i1}}{x-x_{i}}}+{\frac {a_{i2}}{(x-x_ {i})^{2}}}+\cdots +{\frac {a_{ik_{i}}}{(x-x_{i})^{k_{i}}}}\right).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Позволять
![{\ displaystyle g_ {ij} (x) = (x-x_ {i}) ^ {j-1} f (x),}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
единственности ряда Лоранаa ij( x − x i ) −1g ijxx iвычета![{\displaystyle a_{ij} =\operatorname {Res} (g_{ij},x_{i}).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Это определяется непосредственно формулой
![{\displaystyle a_{ij}={\frac {1}{(k_{i}-j)!}}\lim _{x\to x_{i}}{\frac {d^{k_{i}- j}}{dx^{k_{i}-j}}}\left((x-x_{i})^{k_{i}}f(x)\right),}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
x i![{\displaystyle a_{i1}={\frac {P(x_{i})}{Q'(x_{i})}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f(x)={\frac {P(x)}{Q(x)}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Над реальными событиями
Частные дроби используются в интегральном исчислении с действительными переменными для нахождения действительных первообразных рациональных функций . Разложение действительных рациональных функций в частичные дроби также используется для нахождения их обратных преобразований Лапласа . О применении разложения частичных дробей по действительным числам см.
Общий результат
Пусть – любая рациональная функция над действительными числами . Другими словами, предположим, что существуют действительные полиномиальные функции и такие, что![{\ displaystyle f (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle p (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle q (x) \ neq 0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f(x)={\frac {p(x)}{q(x)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Разделив числитель и знаменатель на старший коэффициент , мы можем без ограничения общности предположить , что . По основной теореме алгебры мы можем написать![{\ displaystyle q (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle q (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle q(x)=(xa_{1})^{j_{1}}\cdots (xa_{m})^{j_{m}}(x^{2}+b_{1 }x+c_{1})^{k_{1}}\cdots (x^{2}+b_{n}x+c_{n})^{k_{n}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
где , , являются действительными числами с , и , являются положительными целыми числами. Члены - это линейные факторы , которые соответствуют действительным корням , и термины - неприводимые квадратичные факторы , которые соответствуют парам комплексно -сопряженных корней .![{\displaystyle a_{1},\dots,a_{m}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle b_{1},\dots,b_{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle c_{1},\dots,c_{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle b_{i}^{2}-4c_{i}<0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle j_{1},\dots,j_{m}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle k_{1},\dots,k_{n}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle (xa_{i})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle q (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle q (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle (x_{i}^{2}+b_{i}x+c_{i})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle q (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ displaystyle q (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Тогда разложение на частичные дроби будет следующим:![{\ displaystyle f (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f(x)={\frac {p(x)}{q(x)}}=P(x)+\sum _{i=1}^{m}\sum _{r=1} ^{j_{i}}{\frac {A_{ir}}{(x-a_{i})^{r}}}+\sum _{i=1}^{n}\sum _{r= 1}^{k_{i}}{\frac {B_{ir}x+C_{ir}}{(x^{2}+b_{i}x+c_{i})^{r}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Здесь P ( x ) — полином (возможно, нулевой), а A ir , B ir и C ir — действительные константы. Константы можно найти несколькими способами.
Самый простой метод — умножить на общий знаменатель q ( x ). Затем мы получаем уравнение полиномов, левая часть которого равна просто p ( x ), а правая часть имеет коэффициенты, которые являются линейными выражениями констант A ir , B ir и C ir . Поскольку два многочлена равны тогда и только тогда, когда равны их соответствующие коэффициенты, мы можем приравнять коэффициенты одинаковых членов. Таким образом получается система линейных уравнений, всегда имеющая единственное решение. Это решение можно найти любым из стандартных методов линейной алгебры . Его также можно найти с пределами (см. пример 5).
Примеры
Пример 1
![{\displaystyle f(x)={\frac {1}{x^{2}+2x-3}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Здесь знаменатель разбивается на два отдельных линейных фактора:
![{\displaystyle q(x)=x^{2}+2x-3=(x+3)(x-1)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Итак, мы имеем разложение на частичные дроби
![{\displaystyle f(x)={\frac {1}{x^{2}+2x-3}}={\frac {A}{x+3}}+{\frac {B}{x-1 }}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Умножение на знаменатель в левой части дает нам полиномиальное тождество
![{\displaystyle 1=A(x-1)+B(x+3)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Подстановка x = −3 в это уравнение дает A = −1/4, а замена x = 1 дает B = 1/4, так что
![{\displaystyle f(x)={\frac {1}{x^{2}+2x-3}}={\frac {1}{4}}\left({\frac {-1}{x+ 3}}+{\frac {1}{x-1}}\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Пример 2
![{\displaystyle f(x)={\frac {x^{3}+16}{x^{3}-4x^{2}+8x}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
После длительного деления имеем
![{\displaystyle f(x)=1+{\frac {4x^{2}-8x+16}{x^{3}-4x^{2}+8x}}=1+{\frac {4x^{ 2}-8x+16}{x(x^{2}-4x+8)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Множитель x 2 − 4 x + 8 неприводим по действительным числам, поскольку его дискриминант (−4) 2 − 4×8 = −16 отрицательен. Таким образом, разложение частных дробей по действительным числам имеет вид
![{\displaystyle {\frac {4x^{2}-8x+16}{x(x^{2}-4x+8)}}={\frac {A}{x}}+{\frac {Bx+ C}{x^{2}-4x+8}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Умножив на x 3 − 4 x 2 + 8 x , мы получим полиномиальное тождество
![{\displaystyle 4x^{2}-8x+16=A\left(x^{2}-4x+8\right)+\left(Bx+C\right)x}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Взяв x = 0, мы видим, что 16 = 8 A , поэтому A = 2. Сравнивая коэффициенты x 2 , мы видим, что 4 = A + B = 2 + B , поэтому B = 2. Сравнивая линейные коэффициенты, мы видим, что — 8 = −4 A + C = −8 + C , поэтому C = 0. В целом,
![{\displaystyle f(x)=1+2\left({\frac {1}{x}}+{\frac {x}{x^{2}-4x+8}\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Дробь можно полностью разложить с помощью комплексных чисел . Согласно основной теореме алгебры каждый комплексный многочлен степени n имеет n (комплексных) корней (некоторые из которых могут повторяться). Вторую фракцию можно разложить на:
![{\displaystyle {\frac {x}{x^{2}-4x+8}}={\frac {D}{x-(2+2i)}}+{\frac {E}{x-(2 -2и)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Умножение на знаменатель дает:
![{\ displaystyle x = D (x- (2-2i)) + E (x- (2 + 2i))}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Приравнивая коэффициенты при x и постоянные (по отношению к x ) коэффициенты обеих частей этого уравнения, получаем систему двух линейных уравнений относительно D и E , решением которой является
![{\displaystyle D={\frac {1+i}{2i}}={\frac {1-i}{2}},\qquad E={\frac {1-i}{-2i}}={ \frac {1+i}{2}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Таким образом, мы имеем полное разложение:
![{\displaystyle f(x)={\frac {x^{3}+16}{x^{3}-4x^{2}+8x}}=1+{\frac {2}{x}}+ {\frac {1-i}{x-(2+2i)}}+{\frac {1+i}{x-(2-2i)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Можно также напрямую вычислить A , D и E с помощью метода вычетов (см. также пример 4 ниже).
Пример 3
Этот пример иллюстрирует почти все «трюки», которые нам, возможно, придется использовать, если не считать обращения к системе компьютерной алгебры .
![{\displaystyle f(x)={\frac {x^{9}-2x^{6}+2x^{5}-7x^{4}+13x^{3}-11x^{2}+12x- 4}{x^{7}-3x^{6}+5x^{5}-7x^{4}+7x^{3}-5x^{2}+3x-1}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
После долгого деления и факторизации знаменателя получим
![{\displaystyle f(x)=x^{2}+3x+4+{\frac {2x^{6}-4x^{5}+5x^{4}-3x^{3}+x^{2 }+3x}{(x-1)^{3}(x^{2}+1)^{2}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Разложение на частичные дроби принимает вид
![{\displaystyle {\frac {2x^{6}-4x^{5}+5x^{4}-3x^{3}+x^{2}+3x}{(x-1)^{3}( x^{2}+1)^{2}}}={\frac {A}{x-1}}+{\frac {B}{(x-1)^{2}}}+{\frac {C}{(x-1)^{3}}}+{\frac {Dx+E}{x^{2}+1}}+{\frac {Fx+G}{(x^{2} +1)^{2}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Умножив на знаменатель в левой части, получим полиномиальное тождество
![{\displaystyle {\begin{aligned}&2x^{6}-4x^{5}+5x^{4}-3x^{3}+x^{2}+3x\\[4pt]={}&A\ влево(x-1\вправо)^{2}\влево(x^{2}+1\вправо)^{2}+B\влево(x-1\вправо)\влево(x^{2}+1 \right)^{2}+C\left(x^{2}+1\right)^{2}+\left(Dx+E\right)\left(x-1\right)^{3}\ влево(x^{2}+1\вправо)+\влево(Fx+G\вправо)\влево(x-1\вправо)^{3}\end{aligned}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Теперь мы используем разные значения x для вычисления коэффициентов:
![{\displaystyle {\begin{cases}4=4C&x=1\\2+2i=(Fi+G)(2+2i)&x=i\\0=AB+CEG&x=0\end{cases}} }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Решая это, мы имеем:
![{\displaystyle {\begin{cases}C=1\\F=0,G=1\\E=AB\end{cases}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Используя эти значения, мы можем написать:
![{\displaystyle {\begin{aligned}&2x^{6}-4x^{5}+5x^{4}-3x^{3}+x^{2}+3x\\[4pt]={}&A\ влево(x-1\вправо)^{2}\влево(x^{2}+1\вправо)^{2}+B\влево(x-1\вправо)\влево(x^{2}+1 \right)^{2}+\left(x^{2}+1\right)^{2}+\left(Dx+\left(AB\right)\right)\left(x-1\right)^ {3}\left(x^{2}+1\right)+\left(x-1\right)^{3}\\[4pt]={}&\left(A+D\right)x^ {6}+\left(-A-3D\right)x^{5}+\left(2B+4D+1\right)x^{4}+\left(-2B-4D+1\right)x ^{3}+\left(-A+2B+3D-1\right)x^{2}+\left(A-2B-D+3\right)x\end{aligned}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Сравниваем коэффициенты при х 6 и х 5 с обеих сторон и имеем:
![{\displaystyle {\begin{cases}A+D=2\\-A-3D=-4\end{cases}}\quad \Rightarrow \quad A=D=1.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Поэтому:
![{\displaystyle 2x^{6}-4x^{5}+5x^{4}-3x^{3}+x^{2}+3x=2x^{6}-4x^{5}+(2B+ 5)x^{4}+(-2B-3)x^{3}+(2B+1)x^{2}+(-2B+3)x}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
что дает нам B = 0. Таким образом, разложение на частичные дроби определяется следующим образом:
![{\displaystyle f(x)=x^{2}+3x+4+{\frac {1}{(x-1)}}+{\frac {1}{(x-1)^{3}} }+{\frac {x+1}{x^{2}+1}}+{\frac {1}{(x^{2}+1)^{2}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Альтернативно, вместо разложения можно получить другие линейные зависимости от коэффициентов, вычисляющих некоторые производные при приведенном выше полиномиальном тождестве. (С этой целью напомним, что производная в точке x = a от ( x − a ) m p ( x ) обращается в нуль, если m > 1, и равна просто p ( a ) для m = 1.) Например, первая производная в точке x = 1 дает![{\ displaystyle x = 1, \ imath }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle 2\cdot 6-4\cdot 5+5\cdot 4-3\cdot 3+2+3=A\cdot (0+0)+B\cdot (4+0)+8+D\cdot 0}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
то есть 8 = 4 B + 8, поэтому B = 0.
Пример 4 (метод остатка)
![{\displaystyle f(z)={\frac {z^{2}-5}{(z^{2}-1)(z^{2}+1)}}={\frac {z^{2 }-5}{(z+1)(z-1)(z+i)(zi)}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Таким образом, f ( z ) можно разложить на рациональные функции, знаменателями которых являются z +1, z −1, z +i, z −i. Поскольку каждый член имеет степень единица, −1, 1, − i и i являются простыми полюсами.
Следовательно, остатки, связанные с каждым полюсом, заданные формулой
![{\displaystyle {\frac {P(z_{i})}{Q'(z_{i})}}={\frac {z_{i}^{2}-5}{4z_{i}^{3 }}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle 1,-1,{\tfrac {3i}{2}},-{\tfrac {3i}{2}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle f(z)={\frac {1}{z+1}}-{\frac {1}{z-1}}+{\frac {3i}{2}}{\frac {1} {z+i}}-{\frac {3i}{2}}{\frac {1}{zi}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Пример 5 (метод ограничения)
Пределы можно использовать для нахождения разложения на частичные дроби. [4] Рассмотрим следующий пример:
![{\displaystyle {\frac {1}{x^{3}-1}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Сначала разложите на множители знаменатель, определяющий разложение:
![{\displaystyle {\frac {1}{x^{3}-1}}={\frac {1}{(x-1)(x^{2}+x+1)}}={\frac { A}{x-1}}+{\frac {Bx+C}{x^{2}+x+1}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Умножив все на и приняв предел при , получим![{\displaystyle x-1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle x\to 1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \lim _{x\to 1}\left((x-1)\left({\frac {A}{x-1}}+{\frac {Bx+C}{x^{2} +x+1}}\right)\right)=\lim _{x\to 1}A+\lim _{x\to 1}{\frac {(x-1)(Bx+C)}{x^ {2}+x+1}}=А.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
С другой стороны,
![{\displaystyle \lim _{x\to 1}{\frac {(x-1)}{(x-1)(x^{2}+x+1)}}=\lim _{x\to 1 }{\frac {1}{x^{2}+x+1}}={\frac {1}{3}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
и поэтому:
![{\displaystyle A={\frac {1}{3}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Умножив на x и приняв предел при , получим![{\displaystyle x\to \infty}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \lim _{x\to \infty }x\left({\frac {A}{x-1}}+{\frac {Bx+C}{x^{2}+x+1}} \right)=\lim _{x\to \infty }{\frac {Ax}{x-1}}+\lim _{x\to \infty }{\frac {Bx^{2}+Cx}{ х^{2}+х+1}}=А+В,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
и
![{\displaystyle \lim _ {x\to \infty }{\frac {x}{(x-1)(x^{2}+x+1)}}=0.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Отсюда следует, что A + B = 0 и т. д .![{\displaystyle B=-{\frac {1}{3}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
При x = 0 получаем и таким образом .![{\displaystyle -1=-A+C,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle C=-{\tfrac {2}{3}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Сложив все вместе, получим разложение
![{\displaystyle {\frac {1}{x^{3}-1}}={\frac {1}{3}}\left({\frac {1}{x-1}}+{\frac { -x-2}{x^{2}+x+1}}\вправо).}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Пример 6 (интегральный)
Предположим, у нас есть неопределенный интеграл :
![{\displaystyle \int {\frac {x^{4}+x^{3}+x^{2}+1}{x^{2}+x-2}}\,dx}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Очевидно, что перед выполнением разложения мы должны выполнить полиномиальное деление в длину и факторизовать знаменатель. Это приведет к:
![{\displaystyle \int \left(x^{2}+3+{\frac {-3x+7}{(x+2)(x-1)}}\right)dx}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
После этого мы теперь можем выполнить разложение на частичные дроби.
![{\displaystyle \int \left(x^{2}+3+{\frac {-3x+7}{(x+2)(x-1)}}\right)dx=\int \left(x^ {2}+3+{\frac {A}{(x+2)}}+{\frac {B}{(x-1)}}\right)dx}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle A(x-1)+B(x+2)=-3x+7}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle A={\frac {-13}{3}}\ ,B={\frac {4}{3}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Включение всего этого обратно в наш интеграл позволяет нам найти ответ:
![{\displaystyle \int \left(x^{2}+3+{\frac {-13/3}{(x+2)}}+{\frac {4/3}{(x-1)}} \right)\,dx={\frac {x^{3}}{3}}\ +3x-{\frac {13}{3}}\ln(|x+2|)+{\frac {4 }{3}}\ln(|x-1|)+C}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Роль полинома Тейлора
Разложение рациональной функции в частные дроби можно связать с теоремой Тейлора следующим образом. Позволять
![{\ displaystyle P (x), Q (x), A_ {1} (x), \ ldots, A_ {r} (x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
быть действительными или комплексными полиномами, предполагают, что
![{\displaystyle Q=\prod _{j=1}^{r}(x-\lambda _{j})^{\nu _{j}},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
удовлетворяет
![{\displaystyle \deg A_{1}<\nu _{1},\ldots,\deg A_{r}<\nu _{r},\quad {\text{and}}\quad \deg(P) <\deg(Q)=\sum _{j=1}^{r}\nu _{j}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Также определите
![{\displaystyle Q_{i}=\prod _{j\neq i}(x-\lambda _{j})^{\nu _{j}}={\frac {Q}{(x-\lambda _ {i})^{\nu _{i}}}},\qquad 1\leqslant i\leqslant r.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Тогда у нас есть
![{\displaystyle {\frac {P}{Q}}=\sum _{j=1}^{r}{\frac {A_{j}}{(x-\lambda _{j})^{\nu _{j}}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
тогда и только тогда, когда каждый многочлен является многочленом Тейлора порядка в точке :![{\displaystyle A_{i}(x)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\tfrac {P}{Q_{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \nu _{i}-1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \lambda _{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle A_{i}(x):=\sum _{k=0}^{\nu _{i}-1}{\frac {1}{k!}}\left({\frac {P }{Q_{i}}}\right)^{(k)}(\lambda _{i})\ (x-\lambda _{i})^{k}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Теорема Тейлора (в вещественном или комплексном случае) затем обеспечивает доказательство существования и единственности разложения на частные дроби, а также характеристику коэффициентов.
Эскиз доказательства
Приведенное выше разложение на частичные дроби подразумевает для каждого 1 ≤ i ≤ r полиномиальное разложение
![{\displaystyle {\frac {P}{Q_{i}}}=A_{i}+O((x-\lambda _{i})^{\nu _{i}}),\qquad {\text {для }}x\to \lambda _{i},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
то же самое происходит и с полиномом Тейлора из-за единственности полиномиального разложения порядка и по предположению .![{\displaystyle A_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\tfrac {P}{Q_{i}}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \nu _{i}-1}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg A_{i}<\nu _{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
И наоборот, если это полиномы Тейлора, приведенные выше разложения в каждом случае выполняются, поэтому мы также имеем![{\displaystyle A_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \lambda _{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle P-Q_{i}A_{i}=O((x-\lambda _{i})^{\nu _{i}}),\qquad {\text{for }}x\to \ лямбда _{i},}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
откуда следует, что многочлен делится на![{\displaystyle P-Q_{i}A_{i}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle (x-\lambda _{i})^{\nu _{i}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
For также делится на , поэтому![{\displaystyle j\neq я,Q_{j}A_{j}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle (x-\lambda _{i})^{\nu _{i}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle P-\sum _{j=1}^{r}Q_{j}A_{j}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
делится на . С![{\displaystyle Q}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \deg \left(P-\sum _{j=1}^{r}Q_{j}A_{j}\right)<\deg(Q)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
тогда у нас есть
![{\displaystyle P-\sum _{j=1}^{r}Q_{j}A_{j}=0,}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
и находим разложение на простейшие дроби при делении на .![{\displaystyle Q}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Дроби целых чисел
Идею простейших дробей можно обобщить и на другие области целости , например, на кольцо целых чисел , где простые числа играют роль неприводимых знаменателей. Например:
![{\displaystyle {\frac {1}{18}}={\frac {1}{2}}-{\frac {1}{3}}-{\frac {1}{3^{2}}} .}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Примечания
- ^ Ларсон, Рон (2016). Алгебра и тригонометрия. Cengage Обучение. ISBN 9781337271172.
- ^ Горовиц, Эллис. «Алгоритмы разложения частичных дробей и интегрирования рациональных функций». Материалы второго симпозиума ACM по символическим и алгебраическим манипуляциям. АКМ, 1971.
- ^ Гросхольц, Эмили (2000). Рост математических знаний . Академические издательства Kluwer. п. 179. ИСБН 978-90-481-5391-6.
- ^ Блюман, Джордж В. (1984). Задачник для первого курса по математическому анализу . Нью-Йорк: Springer-Verlag. стр. 250–251.
Рекомендации
- Рао, КР; Ахмед, Н. (1968). «Рекурсивные методы получения разложения рациональной функции в частные дроби». IEEE Транс. Образование . 11 (2): 152–154. Бибкод : 1968ITEdu..11..152R. дои : 10.1109/TE.1968.4320370.
- Хенрици, Питер (1971). «Алгоритм неполного разложения рациональной функции на простейшие дроби». З. Энджью. Математика. Физ . 22 (4): 751–755. Бибкод :1971ЗаМП...22..751Х. дои : 10.1007/BF01587772. S2CID 120554693.
- Чанг, Фэн-Ченг (1973). «Рекурсивные формулы для разложения рациональной функции в частные дроби с кратными полюсами». Учеб. ИИЭЭ . 61 (8): 1139–1140. дои : 10.1109/PROC.1973.9216.
- Кунг, ХТ; Тонг, DM (1977). «Быстрые алгоритмы разложения частичных дробей». SIAM Journal по вычислительной технике . 6 (3): 582. дои : 10.1137/0206042. S2CID 5857432.
- Юстис, Дэн; Кламкин, М.С. (1979). «О коэффициентах разложения простейших дробей». Американский математический ежемесячник . Том. 86, нет. 6. С. 478–480. JSTOR 2320421.
- Махони, Джей-Джей; Сивазлян, Б.Д. (1983). «Разложение неполных дробей: обзор вычислительной методологии и эффективности». Дж. Компьютер. Прил. Математика . 9 (3): 247–269. дои : 10.1016/0377-0427(83)90018-3 .
- Миллер, Чарльз Д.; Лиал, Маргарет Л.; Шнайдер, Дэвид И. (1990). Основы студенческой алгебры (3-е изд.). Addison-Wesley Educational Publishers, Inc., стр. 364–370. ISBN 0-673-38638-4.
- Вестрайх, Дэвид (1991). «Разложение частичных дробей без оценки производной». IEEE Транс. Цирк. Сист . 38 (6): 658–660. дои : 10.1109/31.81863.
- Кудрявцев, Л.Д. (2001) [1994], «Неопределенные коэффициенты, метод», Энциклопедия математики , EMS Press
- Веллеман, Дэниел Дж. (2002). «Частичные дроби, биномиальные коэффициенты и интеграл нечетной степени сек тета». амер. Математика. Ежемесячно . 109 (8): 746–749. дои : 10.2307/3072399. JSTOR 3072399.
- Слота, Дамиан; Витула, Роман (2005). «Трёхкиричный метод разложения частичных дробей некоторого типа рационального выражения». Вычислительная наука – ICCS 2005 . Лект. Нет. Компьютерные науки. Том. 33516. стр. 659–662. дои : 10.1007/11428862_89. ISBN 978-3-540-26044-8.
- Кунг, Сидни Х. (2006). «Разложение частичных дробей делением». Колл. Математика. Дж . 37 (2): 132–134. дои : 10.2307/27646303. JSTOR 27646303.
- Витула, Роман; Слота, Дамиан (2008). «Разложение некоторых рациональных функций в частные дроби». Прил. Математика. Вычислить . 197 : 328–336. дои : 10.1016/j.amc.2007.07.048. МР 2396331.
Внешние ссылки