Разложение рациональной функции в ряд — различия между версиями
(→Общий алгоритм) |
(→Общий алгоритм) |
||
Строка 29: | Строка 29: | ||
==Общий алгоритм== | ==Общий алгоритм== | ||
# Привести дробь <tex>\dfrac{P(z)}{Q(z)}</tex> к такому виду, чтобы степень числителя была меньше степени знаменателя. Если <tex>\deg(P) > \deg(Q)</tex>, то можем записать <tex>G(z)=\dfrac{P(z)}{Q(z)} = R(z)+\dfrac{P_0(z)}{Q(z)},</tex> где <tex>\deg(P_0) < \deg(Q)</tex>. | # Привести дробь <tex>\dfrac{P(z)}{Q(z)}</tex> к такому виду, чтобы степень числителя была меньше степени знаменателя. Если <tex>\deg(P) > \deg(Q)</tex>, то можем записать <tex>G(z)=\dfrac{P(z)}{Q(z)} = R(z)+\dfrac{P_0(z)}{Q(z)},</tex> где <tex>\deg(P_0) < \deg(Q)</tex>. | ||
+ | # Отыскать корни уравнения <tex>Q(z)=0</tex> и разбить знаменатель на множители вида <tex>(z_s−z)^{k_s}</tex> (здесь <tex>z_s</tex> — корень кратности <tex>k_s<tex>). | ||
# Разобьем знаменатель <tex>Q(z)</tex> на множители <tex>Q(z) = (z_k-z)^{k_s} *\ldots</tex>, где <tex>z_1, z_2, \ldots, z_s</tex> - корни уравнения <tex>Q(z) = 0</tex>. При этом, <tex>k_1+k_2+\ldots+k_s=\deg (Q)</tex>. После разбиения знаменателя на множители получим: <tex>G(z)=\dfrac{P(z)}{(z1-z)^k1 *...(zs-z)^ks}</tex> (k1, ks - сделать индексами) | # Разобьем знаменатель <tex>Q(z)</tex> на множители <tex>Q(z) = (z_k-z)^{k_s} *\ldots</tex>, где <tex>z_1, z_2, \ldots, z_s</tex> - корни уравнения <tex>Q(z) = 0</tex>. При этом, <tex>k_1+k_2+\ldots+k_s=\deg (Q)</tex>. После разбиения знаменателя на множители получим: <tex>G(z)=\dfrac{P(z)}{(z1-z)^k1 *...(zs-z)^ks}</tex> (k1, ks - сделать индексами) | ||
# Приведем G(z) к сумме дробей, знаменатели которых будут иметь вид (zj−z)^kj, а числители — полиномы Pj(z), причем deg Pj(z)<kj. <tex>G(z)=\dfrac{P(z)}{(z1-z)^k1 *...(zs-z)^ks} = \sum\limits \dfrac{Pj(z)}{(zj-z)^kj}</tex>. Найдем Pj(z) с помощью [[Разложение рациональной функции в ряд#Метод неопределенных коэффициентов|метода неопределенных коэффициентов]]. | # Приведем G(z) к сумме дробей, знаменатели которых будут иметь вид (zj−z)^kj, а числители — полиномы Pj(z), причем deg Pj(z)<kj. <tex>G(z)=\dfrac{P(z)}{(z1-z)^k1 *...(zs-z)^ks} = \sum\limits \dfrac{Pj(z)}{(zj-z)^kj}</tex>. Найдем Pj(z) с помощью [[Разложение рациональной функции в ряд#Метод неопределенных коэффициентов|метода неопределенных коэффициентов]]. |
Версия 14:22, 28 мая 2017
Определения
Определение: |
Рациональная функция — это функция вида:
, |
Рациональные производящие функции получаются при решении линейных рекуррентных соотношений. По этой причине актуальной является задача о разложении рациональной функции в ряд по степеням переменной .
Чтобы разложить дробь в ряд, необходимо разбить её на сумму элементарных дробей.
Определение: |
Элементарными дробями будем называть дроби вида:
, |
Затем, элементарные дроби сможем разложить в ряд, пользуясь формулами преобразования производящих функций и таблицей производящих функций.
Общий алгоритм
- Привести дробь к такому виду, чтобы степень числителя была меньше степени знаменателя. Если , то можем записать где .
- Отыскать корни уравнения и разбить знаменатель на множители вида (здесь — корень кратности на множители , где - корни уравнения . При этом, . После разбиения знаменателя на множители получим: (k1, ks - сделать индексами)
- Приведем G(z) к сумме дробей, знаменатели которых будут иметь вид (zj−z)^kj, а числители — полиномы Pj(z), причем deg Pj(z)<kj. метода неопределенных коэффициентов. . Найдем Pj(z) с помощью
Метод неопределенных коэффициентов
- Записать сумму дробей, знаменатили которых будут иметь вид (zs−z)ks, а числители — полиномы с неопределёнными коэффициентами, имеющие степень ks−1.
- Сложить выписанные дроби и сгруппировать слагаемые в числителе по степеням z.
- Прировнять полученные выражения с неопределёнными коэффициентами к соответсвующим коэффициентам полинома P(z), составив, таким образом, систему линейных уравнений.
- Решить систему и получить значения неопределённых коэффициентов.