Частичное дробное разложение символьных выражений

Найдите частичное разложение дробей одномерных и многомерных выражений.

Во-первых, найдите частичное разложение дробей одномерных выражений. Для выражений с одной переменной можно опустить указание переменной.

syms x
partfrac(x^2/(x^3 - 3*x + 2))

ans =
5/(9*(x - 1)) + 1/(3*(x - 1)^2) + 4/(9*(x + 2))

Найдите частичное разложение дроби многомерного выражения относительно конкретной переменной.

syms a b
partfrac(a^2/(a^2 - b^2),a)

ans =
b/(2*(a - b)) - b/(2*(a + b)) + 1

partfrac(a^2/(a^2 - b^2),b)

ans =
a/(2*(a + b)) + a/(2*(a - b))

Если вы не задаете переменную, то partfrac вычисляет разложение частичной дроби относительно переменной, определяемой как symvar.

symvar(a^2/(a^2 - b^2),1)
partfrac(a^2/(a^2 - b^2))

ans =
b
 
ans =
a/(2*(a + b)) + a/(2*(a - b))

Режимы факторизации

Выберите конкретный режим факторизации при помощи FactorMode вход.

Найдите частичное разложение дроби, не задавая режим факторизации. По умолчанию, partfrac использует факторизацию по рациональным числам. В этом режиме, partfrac содержит числа в их точной символической форме.

syms x
f = 1/(x^3 + 2);
partfrac(f,x)

ans =
1/(x^3 + 2)

Повторите разложение с числовым факторизацией по вещественным числам. В этом режиме, partfrac множит знаменатель в линейные и квадратичные неприводимые полиномы с действительными коэффициентами. Этот режим преобразует все числовые значения в числа с плавающей запятой.

partfrac(f,x,'FactorMode','real')

ans =
0.2099868416491455274612017678797/(x + 1.2599210498948731647672106072782) -...
(0.2099868416491455274612017678797*x - 0.52913368398939982491723521309077)/(x^2 -...
1.2599210498948731647672106072782*x + 1.5874010519681994747517056392723)

Повторите разложение с факторизацией по комплексным числам. В этом режиме, partfrac уменьшает квадратичные полиномы в знаменателе до линейных выражений с комплексными коэффициентами. Этот режим преобразует все числа в плавающие точки.

partfrac(f,x,'FactorMode','complex')

ans =
0.2099868416491455274612017678797/(x + 1.2599210498948731647672106072782) +...
(- 0.10499342082457276373060088393985 - 0.18185393932862023392667876903163i)/...
(x - 0.62996052494743658238360530363911 - 1.0911236359717214035600726141898i) +...
(- 0.10499342082457276373060088393985 + 0.18185393932862023392667876903163i)/...
(x - 0.62996052494743658238360530363911 + 1.0911236359717214035600726141898i)

Найдите частичное разложение дроби этого выражения с помощью полного режима факторизации. В этом режиме, partfrac множит знаменатель в линейные выражения, сводя квадратичные полиномы к линейным выражениям со сложными коэффициентами. Этот режим сохраняет числа в их точной символической форме.

pfFull = partfrac(f,x,'FactorMode','full')

pfFull =
2^(1/3)/(6*(x + 2^(1/3))) +...
(2^(1/3)*((3^(1/2)*1i)/2 - 1/2))/(6*(x + 2^(1/3)*((3^(1/2)*1i)/2 - 1/2))) -...
(2^(1/3)*((3^(1/2)*1i)/2 + 1/2))/(6*(x - 2^(1/3)*((3^(1/2)*1i)/2 + 1/2)))

Аппроксимируйте результат с числами с плавающей запятой при помощи vpa. Потому что выражение не содержит никаких символьных параметров, кроме переменной x, результат такой же, как и в режиме комплексной факторизации.

vpa(pfFull)

ans =
0.2099868416491455274612017678797/(x + 1.2599210498948731647672106072782) +...
(- 0.10499342082457276373060088393985 - 0.18185393932862023392667876903163i)/...
(x - 0.62996052494743658238360530363911 - 1.0911236359717214035600726141898i) +...
(- 0.10499342082457276373060088393985 + 0.18185393932862023392667876903163i)/...
(x - 0.62996052494743658238360530363911 + 1.0911236359717214035600726141898i)

В комплексном режиме, partfrac множители только те выражения в знаменателе, коэффициенты которых могут быть преобразованы в числа с плавающей запятой. Показать это путем замены 2 в f с символьной переменной и найти частичное разложение дроби в комплексном режиме. partfrac возвращает выражение без изменений.

syms a
f = subs(f,2,a);
partfrac(f,x,'FactorMode','complex')

ans =
1/(x^3 + a)

Когда вы используете полный режим факторизации, partfrac выражения факторов в знаменателе символически. Таким образом, partfrac в режиме полного факторизации выражение становится фактором.

partfrac(1/(x^3 + a), x, 'FactorMode', 'full')

ans =
1/(3*(-a)^(2/3)*(x - (-a)^(1/3))) -...
((3^(1/2)*1i)/2 + 1/2)/(3*(-a)^(2/3)*(x + (-a)^(1/3)*((3^(1/2)*1i)/2 + 1/2))) +...
((3^(1/2)*1i)/2 - 1/2)/(3*(-a)^(2/3)*(x - (-a)^(1/3)*((3^(1/2)*1i)/2 - 1/2)))

Полный режим факторизации возвращает root

В режиме полного факторизации, partfrac представляет коэффициенты, использующие root когда математически невозможно найти коэффициенты как точные символьные числа. Показать это поведение.

syms x
s = partfrac(1/(x^3 + x - 3), x, 'FactorMode','full')

s =
symsum(-((6*root(z^3 + z - 3, z, k)^2)/247 +...
         (27*root(z^3 + z - 3, z, k))/247 +...
          4/247)/(root(z^3 + z - 3, z, k) - x), k, 1, 3)

Аппроксимируйте результат с числами с плавающей запятой при помощи vpa.

vpa(s)

ans =
0.1846004942289254798185772017286/(x - 1.2134116627622296341321313773815) +...
(- 0.092300247114462739909288600864302 + 0.11581130283490645120989658654914i)/...
(x + 0.60670583138111481706606568869074 - 1.450612249188441526515442203395i) +...
(- 0.092300247114462739909288600864302 - 0.11581130283490645120989658654914i)/...
(x + 0.60670583138111481706606568869074 + 1.450612249188441526515442203395i)

Числители и знаменатели частичного разложения дробей

Верните вектор числителей и вектор знаменателей разложения частичной дроби.

Сначала найдите частичное дробное разложение выражения.

syms x
P = partfrac(x^2/(x^3 - 3*x + 2), x)

P =
5/(9*(x - 1)) + 1/(3*(x - 1)^2) + 4/(9*(x + 2))

Частичное разложение дроби является суммой дробей. Используйте children функция для возврата вектора, содержащего условия этой суммы. Затем используйте numden извлечь числители и знаменатели терминов.

[N,D] = numden(children(P))

N =
[ 5, 1, 4]
 
D =
[ 9*x - 9, 3*(x - 1)^2, 9*x + 18]

Восстановите частичное разложение дроби из векторов числителей и знаменателей.

P1 = sum(N./D)

P1 =
1/(3*(x - 1)^2) + 5/(9*x - 9) + 4/(9*x + 18)

Проверьте, что восстановленное выражение, P1, эквивалентно исходному частичному разложению дроби, P.

isAlways(P1 == P)

ans =
  logical
     1