СВОЙСТВА Z-ПРЕОБРАЗОВАНИЯ [2].

Важнейшим свойством z-преобразования является свойство его единственности. Любая последовательность s(k) однозначно определяется z-изображением в области его сходимости, и наоборот, однозначно восстанавливается по z-изображению.

Без углубления в теорию, можно констатировать, что все свойства ДПФ действительны и для z-преобразования. Отметим некоторые из них.

Линейность: Если s(k) = a·x(k)+b·y(k), то S(z) = aX(z)+bY(z). Соответственно, z-преобразование допустимо только для анализа линейных систем и сигналов, удовлетворяющих принципу суперпозиции.

Задержка на n тактов: y(k) = x(k-n).

Y(z) = y(k) z^k = x(k-n) z^k =zⁿ x(k-n) z^k-n = zⁿ x(m) z^m = zⁿ X(z).

Соответственно, умножение z-образа сигнала на множитель zⁿ вызывает сдвиг сигнала на n тактов дискретизации.

Преобразование свертки. При выполнении нерекурсивной цифровой фильтрации односторонними операторами фильтров:

s(k) = h(n) y(k-n), k = 0, 1, 2, …

Z-преобразование уравнения свертки:

S(z) = h(n) y(k-n) z^k = h(n) zⁿ y(k-n) z^k-n =

= h(n) zⁿ y(k-n) z^k-n = H(z) Y(z).

Таким образом, свертка дискретных функций отображается произведением z-образов этих функций. Аналогично, для z-преобразования могут быть доказаны все известные теоремы о свойствах z-образов, что вполне естественно, т.к. при z=exp(-jw) эти свойства полностью эквивалентны свойствам спектров функций.

Разложение сигналов на блоки последовательной свертки. Z-преобразование позволяет производить разложение сигналов и функций, например передаточных функций фильтров, на короткие составляющие операции свертки, для чего достаточно приравнять z-полином к нулю, найти его корни a_i, и переписать полином в виде произведения двучленов:

S(z) = a₀(z-a₁)(z-a₂)...,

где а₀- последний отсчет сигнала (коэффициент при старшей степени z).

Но произведению в z-области соответствует свертка в координатной области, и при обратном преобразовании двучлены (z-a_i) превращаются в двухточечные диполи {-a_i,1}, а сигнал длиной N представляется сверткой (N-1) диполей:

s_k= a₀{-a₁,1}_*{-a₂,1}_*{-a₃,1}_* ...

Пример. s_k = {1.4464, -2.32, 3.37, -3, 1}. S(z) = z⁴-3z³+3.37z²-2.32z+1.4464. a₀ = 1.

Корни полинома S(z): a₁ = 0.8+0.8j, a₂ = 0.8-0.8j, a₃ = 0.7+0.8j, a₄ = 0.7-0.8j,

S(z) = (z-0.8-0.8j)(z-0.8+0.8j)(z-0.7-0.8j)(z-0.7+0.8j).

Корни полинома представлены на z-плоскости на рис. 8.1.1. Корни полинома комплексные и четыре двучлена в координатной области также будут комплексными. Но они являются сопряженными, и для получения вещественных функций следует перемножить сопряженные двучлены и получить биквадратные блоки:

S(z) = (z²-1.4z+1.13)(z²-1.6z+1.28).

При переходе в координатную область: s_k = {1.13, -1.4, 1} _* {1.28, -1.6, 1}.

Таким образом, исходный сигнал разложен на свертку двух трехчленных сигналов (функций).

Дифференцирование. Если имеем s(k) « S(z), то z-образ функции s'(k) можно найти, продифференцировав S(z), что бывает полезно для вычисления обратного z-преобразования функций S(z) с полюсами высокого порядка:

s'(k) « z dS(z)/dz.

8.4. ОБРАТНОЕ Z-ПРЕОБРАЗОВАНИЕ [43]

Методы преобразования. Обратное z-преобразование позволяет восстанавливать дискретную функцию по ее z-образу. Оно широко используется, например, при определении импульсных характеристик рекурсивных цифровых фильтров. В символической форме:

x(k) = TZ^-1[X(z)].

На практике X(z) в процессе расчетов обычно выражается через отношение двух многочленов от z:

X(z) = (b₀+b₁z+b₂z² + …+ b_Nz^N ) / (a₀+a₁z+a₂z² + …+ a_Mz^M ) = (8.4.1)

= x(0) + x(1)z + x(2)z² + … (8.4.1')

Самые распространенные методы обратного преобразования из этой формы X(z):

• Преобразование интегрированием по контуру (метод вычетов).
• Метод разложения на элементарные дроби.
• Метод разложения в степенной ряд.

Метод разложения в степенной ряд наиболее прост и пригоден для выполнения на компьютерах, но он не дает решения в аналитической форме. При задании большого числа точек обратного преобразования требуется также следить за возможным нарастанием числовых ошибок вследствие рекурсии его алгоритма.

Два первых метода позволяют получать результаты в аналитическом виде, но требуют вычисления полюсов функции X(z), что может представлять трудности при высоком порядке функции. При высоких порядках полюсов потребуется также дифференцирование соответствующих порядков.

Преобразование интегрированием по контуру относится к числу математически строгих методов. Оно выполняется интегрированием по произвольному замкнутому контуру C, расположенному в области сходимости и окружающему все особые точки (нули и полюсы) z-образа. Интегрирование удобнее выполнять над полюсами, расположенными внутри контура, включающего центр системы координат, т.е. в символике z^-1. В этой символике мы и будем рассматривать данных параграф. Контурный интеграл обратного преобразования:

s_k = (1/2pj) . (8.4.2)

Согласно теореме Коши о вычетах, интеграл (8.4.2) равен сумме вычетов (Res) подынтегральной функции относительно всех полюсов этой функции, лежащих внутри контура интегрирования. Каждый вычет связан с определенным полюсом p_k:

Res[F(z), p_k] = [(z-p_k) F(z)] при z=p_k. (8.4.3)

где F(z) = z^k-1 S(z), m – порядок полюса в точке p_k. Для простого полюса:

Res[F(z), p_k] = (z-p_k) F(z) = (z-p_k) z^k-1 S(z) при z=p_k. (8.4.3')

Пример.Z-образ функции: X(z) = z² / (z-0.5)(z-1)².

x(k) = Res[F(z), p₁] + Res[F(z), p₂]. F(z) = z^k-1 X(z) = z^k+1 / (z-0.5)(z-1)².

Функция F(z) имеет простой полюс p₁ = 0.5 и полюс второго порядка p₂ = 1.

Res[F(z), 0.5] = (z-0.5) z^k+1 / (z-0.5)(z-1)² = z^k+1 / (z-1)² |_z=0.5 = 0.5 (0.5)^k / (0.5)² = 2(0.5)^k.

Res[F(z), 1] = [(z-1)² z^k+1 / (z-0.5)(z-1)²] = [(z-0.5)(k+1)z^k-z^k+1] / (z-0.5)² |_z=1= 2(k-1).

Результат: x(k) = 2[(k-1) + (0.5)^k].

Преобразование разложением на дроби. В этом методе z-образ (8.4.1) раскладывается на рациональные простые дроби с последующим почленным обратным преобразованием с помощью таблицы. Наиболее просто это выполняется, если функция S(z) может быть разложена по степеням z в символике z^-1, т.е. представлена в следующем виде:

S(z) = s(0) + s(1) z^-1+ s(2) z^-2 + …

Соответственно, в выражении (8.4.1) отношение многочленов также должно быть в символике z^-1. Если полюсы S(z) первого порядка и N = M, то (8.4.1) можно разложить на следующую сумму:

S(z) = B₀ + C₁/(1-p₁z^-1) + C₂/(1-p₂z^-2) + … + C_M/(1-p_Mz^-M) =

B₀ + C₁z/(z-p₁) + C₂z/(z-p₂) + … + C_Mz/(z-p_M) = B₀ + C_kz/(z-p_k). (8.4.4)

B₀ = b_N / a_N.

где С_k – коэффициенты элементарных дробей, которые являются вычетами функции S(z).

Для вычисления коэффициентов C_k умножим левую и правую стороны выражения (8.4.4) на (z-p_k)/z и положим z=p_k, при этом в правой части за счет множителя (z-p_k)=0 при z=p_k обнуляются все члены суммы кроме члена с C_k данного полюса, а в левой остается произведение S(z)(z-p_k)/z, что позволяет вычислить значения C_k:

C_k = S(z)(z-p_k)/z |_z=pk (8.4.5)

Если в (8.4.1) N < M, то значение B₀ равно нулю. Если функция S(z) в точке z=p_k имеет полюс m-ного порядка, то коэффициент C_k заменяется суммой коэффициентов:

D_i /(z-p_k)ⁱ, (8.4.6)

D_i = [ X(z) (z-p_k)^m/z], при z=p_k. (8.4.7)

Пример.Повторим пример преобразования данным способом z-образа функции X(z) = z² / (z-0.5)(z-1)², использованного в предыдущем примере. Функция имеет простой полюс p₁ = 0.5 и полюс второго порядка p₂ = 1.

X(z) = Cz/(z-0.5) + D₁z/(z-1) + D₂z/(z-1)².

С = z/(z-1)² = 0.5/(0.5-1)² = 2.

D₁ = [(z-1)² X(z)/z] = [z / (z-0.5)] |_z=1= -2.

D₂ = (z-1)² X(z)/z = z/(z-0.5) |_z=1= 2.

X(z) = 2z/(z-0.5) + D₁z/(z-1) + D₂z/(z-1)².

Обратное преобразование каждой простой дроби выполним по таблице 8.2.1.

Результат: x(k) = 2(0.5)^k -2 +2k = 2[(k-1) + (0.5)^k]. Результат аналогичен методу вычетов.

Если z-изображение имеет вид дробно-рациональной функции, то разложение на простые дроби с последующим применением таблицы соответствий обычно труда не представляет. Так, например:

S(z) = (b₀+b₁z^-1+b₂z^-2) / (1-az^-1) = b₀/(1-az^-1) + b₁z^-1/(1-az^-1) + b₂z^-2/(1-az^-1).

По таблице соответствия:

X(z) = 1/(1-az^-1) → x(k) = a^k.

Отсюда, с учетом линейности преобразования и свойства задержки:

x(k) = b₀ a^k + b₁ a^k-1 + b₂ a^k-2.

При преобразовании функций со знаменателями более высоких порядков предварительно следует найти полюса функции. Например, для многочлена второго порядка с полюсами p₁ и p₂:

S(z) = 1/(1-a₁z^-1+a₂z^-2) = 1/[(1-p₁z^-1)(1-p₂z^-1).

Представим S(z) в виде суммы дробей с неизвестными коэффициентами b₁ и b₂:

S(z) = b₁/(1-p₁z^-1)+b₂/(1-p₂z^-1) = (b₁- b₁p₂z^-1+b₂-b₂p₁z^-1)/[(1-p₁z^-1)(1-p₂z^-1).

При равенстве знаменателей в этих двух выражениях должны быть равны и числители:

(b₁ + b₂) – (b₁ p₂+b₂ p₁)z^-1 = 1,

а это обеспечивается равенством коэффициентов при одинаковых степенях z. Отсюда получаем систему уравнений:

b₁ + b₂ = 1.

b₁ p₂+b₂ p₁ = 0.

Решая эту систему уравнений, находим значения коэффициентов b₁ и b₂, подставляем коэффициенты в S(z), выраженное в виде суммы дробей, и по таблице соответствия переводим дроби во временные функции.

Метод степенных рядов. Выражение (8.4.1) можно разложить непосредственно в степенной ряд (8.4.1') путем деления в столбик, для чего числитель и знаменатель функции выражаются предварительно через нарастающий или уменьшающийся показатель степени z. Обратное z-преобразование степенного ряда очевидно.

Пример нарастающей степени z.X(z) = (1+2z+z²) / (1-z+0.4z²).

1 + 2z + z² | 1 – z + 0.4z²

1 – z + 0.4z² 1 + 3z + 3.6z² + 2.4z³ + 0.96z⁴ + … Ряд может быть бесконечным.

3z + 0.6z²

3z – 3z² + 1.2z³

3.6z² – 1.2z³

3.6z² – 3.6z³ + 1.44z⁴

2.4z³ – 1.44z⁴

2.4z³ – 2.4z⁴ + 0.96z⁵

0.96z⁴ – 0.96z⁵

0.96z⁴ – 0.96z⁵ + 0.384z⁶, и т.д.

Обратное преобразование выполняется путем идентификации коэффициентов степеней при z^k с k-отсчетами функции: x(k) = {1, 3, 3.6, 2.4, 0.96, …}.

Пример уменьшающихся номеров степени z.X(z) = (1+2z+z²) / (1-z+0.4z²) → (деление на z^N числителя и знаменателя полинома) → (z^-2+2z^-1+1) / (z^-2-z^-1+0.4).

z^-2 + 2z^-1 + 1 | z^-2 – z^-1 + 0.4

z^-2 – z^-1 + 0.4 1 + 3z + 3.6z² + 2.4z³ + 0.96z⁴ + … Результат тот же.

3z^-1 + 0.6

3z^-1 – 3 + 1.2z

3.6 – 1.2z

3.6 – 3.6z + 1.44z²

2.4z – 1.44z²

2.4z – 2.4z² + 0.96z³

0.96z² – 0.96z³

0.96z² – 0.96z³ + 0.384z⁴, и т.д.

Метод деления полинома (8.4.1) можно выполнять рекурсивно:

x(0) = b₀ / a₀,

x(1) = (b₁ – x(0) a₁) / a₀,

x(2) = (b₂ – x(1) a₁ – x(0) a₂) / a₀,

…

x(n) = (b_n – (x(n-i) a_i) /a₀, n = 1, 2, 3, …