28.3 Произведения многочленов

: conv (a, b) ¶

: conv (a, b, shape) ¶

Выполняет свёртку двух векторов a и b.

Когда a и b являются коэффициентными векторами двух многочленов, свёртка представляет собой коэффициентный вектор произведения многочленов.

Размер результата определяется необязательным аргументом shape, который принимает следующие значения

shape = "full"

Возвращает полную свёртку. (по умолчанию) Результат — вектор длиной, равной length (a) + length (b) - 1.

shape = "same"

Возвращает центральную часть свёртки того же размера, что и a.

shape = "valid"

Возвращает только части, не включающие обрезаемые нулями края. Размер результата — max (size (a) - size (b) + 1, 0).

См. также: deconv, conv2, convn, fftconv.

: C = convn (A, B) ¶

: C = convn (A, B, shape) ¶

Возвращает n-мерную свёртку A и B.

Размер результата определяется необязательным аргументом shape, который принимает следующие значения

shape = "full"

Возвращает полную свёртку. (по умолчанию)

shape = "same"

Возвращает центральную часть свёртки того же размера, что и A. Центральная часть свёртки начинается с индексов floor ([size(B)/2] + 1).

shape = "valid"

Возвращает только части, не включающие обрезаемые нулями края. Размер результата — max (size (A) - size (B) + 1, 0).

См. также: conv2, conv.

: b = deconv (y, a) ¶

: [b, r] = deconv (y, a) ¶

Выполняет де свёртку двух векторов (деление многочленов).

[b, r] = deconv (y, a) решает относительно b и r таким образом, чтобы y = conv (a, b) + r.

Если y и a — коэффициентные векторы многочленов, b будет содержать коэффициенты многочлена частного, а r — многочлен остатка наименьшей степени.

См. также: conv, residue.

: conv2 (A, B) ¶

: conv2 (v1, v2, m) ¶

: conv2 (…, shape) ¶

Возвращает 2-мерную свёртку A и B.

Размер результата определяется необязательным аргументом shape, который принимает следующие значения

shape = "full"

Возвращает полную свёртку. (по умолчанию)

shape = "same"

Возвращает центральную часть свёртки того же размера, что и A. Центральная часть свёртки начинается с индексов floor ([size(B)/2] + 1).

shape = "valid"

Возвращает только части, не включающие обрезаемые нулями края. Размер результата — max (size (A) - size (B) + 1, 0).

Когда третий аргумент является матрицей, возвращает свёртку матрицы m по вектору v1 в направлении столбцов и по вектору v2 в направлении строк.

См. также: conv, convn.

: q = polygcd (b, a) ¶

: q = polygcd (b, a, tol) ¶

Находит наибольший общий делитель двух многочленов.

Это эквивалентно многочлену, полученному умножением всех общих корней. Вместе с де свёрткой вы можете упростить отношение двух многочленов.

Допуск tol по умолчанию равен sqrt (eps).

Внимание: Этот алгоритм числово неустойчив и не должен использоваться с большими многочленами.

Пример кода:

polygcd (poly (1:8), poly (3:12)) - poly (3:8)
⇒ [ 0, 0, 0, 0, 0, 0, 0 ]
deconv (poly (1:8), polygcd (poly (1:8), poly (3:12))) - poly (1:2)
⇒ [ 0, 0, 0 ]

См. также: poly, roots, conv, deconv, residue.

: [r, p, k, e] = residue (b, a) ¶

: [b, a] = residue (r, p, k) ¶

: [b, a] = residue (r, p, k, e) ¶

Первый вызов вычисляет разложение в частичные дроби для частного многочленов, b и a.

Частное определяется как

B(s)    M       r(m)        N
---- = SUM ------------- + SUM k(i)*s^(N-i)
A(s)   m=1 (s-p(m))^e(m)   i=1

где M — количество полюсов (длина векторов r, p и e), вектор k — многочлен порядка N-1, представляющий прямой вклад, а вектор e задаёт кратность полюса m-го остатка.

Например,

b = [1, 1, 1];
a = [1, -5, 8, -4];
[r, p, k, e] = residue (b, a)
   ⇒ r = [-2; 7; 3]
   ⇒ p = [2; 2; 1]
   ⇒ k = [](0x0)
   ⇒ e = [1; 2; 1]

что представляет собой следующее разложение в частичные дроби

s^2 + s + 1       -2        7        3
   ------------------- = ----- + ------- + -----
   s^3 - 5s^2 + 8s - 4   (s-2)   (s-2)^2   (s-1)

Второй вызов выполняет обратную операцию и вычисляет воссозданное частное многочленов, b(s)/a(s), из разложения в частичные дроби; представленного остатками, полюсами и прямым многочленом, заданным r, p и k, и кратностью полюса e.

Если кратность e не задана явно, кратность определяется функцией mpoles.

Например:

r = [-2; 7; 3];
p = [2; 2; 1];
k = [1, 0];
[b, a] = residue (r, p, k)
   ⇒ b = [1, -5, 9, -3, 1]
   ⇒ a = [1, -5, 8, -4]

where mpoles is used to determine e = [1; 2; 1]

В качестве альтернативы кратность может быть определена явно, например,

r = [7; 3; -2];
p = [2; 1; 2];
k = [1, 0];
e = [2; 1; 1];
[b, a] = residue (r, p, k, e)
   ⇒ b = [1, -5, 9, -3, 1]
   ⇒ a = [1, -5, 8, -4]

что представляет собой следующее разложение в частичные дроби

-2        7        3         s^4 - 5s^3 + 9s^2 - 3s + 1
----- + ------- + ----- + s = --------------------------
(s-2)   (s-2)^2   (s-1)          s^3 - 5s^2 + 8s - 4

См. также: mpoles, poly, roots, conv, deconv.