Z-преобразование

2004 год. Инженеры Dolby разрабатывают Dolby Digital Plus для HD DVD. Задача: спроектировать IIR-фильтр с крутизной -80 дБ/окт при задержке не более 2 мс. Решение нашли за день - на листе бумаги. Нарисовали диаграмму полюсов и нулей, передвинули несколько точек в Z-плоскости. Полюса внутри единичного круга - фильтр устойчив. Нули на нужных частотах - подавление достигнуто. Z-преобразование превращает проектирование фильтров из уравнений в геометрию.

• **Аудио:** эквалайзер = набор IIR-фильтров, каждый спроектирован через размещение полюсов/нулей в Z-плоскости
• **Телеком:** канальные эквалайзеры в 5G/WiFi используют Z-анализ для компенсации искажений канала
• **Управление:** цифровые PID-регуляторы проектируются через передаточную функцию H(z) с гарантией устойчивости

Рагаззини, Заде и рождение цифровой теории управления

Z-преобразование введено Рагаззини и Заде в 1952 году для анализа цифровых систем управления - дискретный аналог преобразования Лапласа. Формально это ряд Лорана, инструмент из комплексного анализа 19-го века, применённый к инженерии. Лотфи Заде через 13 лет создаст нечёткую логику - но Z-преобразование осталось его главным вкладом в DSP.

Предварительные знания

• Sampling and the Nyquist-Shannon Theorem

Определение Z-преобразования

**Z-преобразование** - обобщение дискретного преобразования Фурье: $X(z) = \sum_{n=-\infty}^{\infty} x[n] \cdot z^{-n}$, где z - комплексная переменная. Фурье работает на единичной окружности ($z = e^{j\omega}$). Z-преобразование - на всей комплексной плоскости.

Зачем? Фурье может не существовать для растущих сигналов ($a^n$ при $|a|>1$). Z-преобразование расширяет область применения, позволяя анализировать неустойчивые системы и переходные процессы. Это как аналитическое продолжение Фурье за пределы единичной окружности.

**Ключевое свойство задержки:** $\mathcal{Z}\{x[n-k]\} = z^{-k} \cdot X(z)$. Задержка на k отсчётов = умножение на $z^{-k}$. Поэтому $z^{-1}$ называют «оператором задержки». Все цифровые фильтры строятся из задержек, умножителей и сумматоров.

Область сходимости (ROC)

**ROC (Region of Convergence)** - множество точек z, для которых ряд $X(z) = \sum x[n] \cdot z^{-n}$ сходится. ROC критически важна: одна и та же формула X(z) может соответствовать РАЗНЫМ сигналам в зависимости от ROC!

Пример: $X(z) = z/(z-a)$. Если ROC: $|z| > |a|$ - это правосторонний $a^n \cdot u[n]$ (каузальный). Если ROC: $|z| < |a|$ - это левосторонний $-a^n \cdot u[-n-1]$ (антикаузальный). Формула та же, сигналы разные!

**Устойчивость и ROC:** LTI-система устойчива (BIBO) тогда и только тогда, когда ROC её передаточной функции H(z) содержит **единичную окружность** |z| = 1. На единичной окружности Z-преобразование = DTFT (частотная характеристика).

ROC всегда имеет форму кольца (или круга, или внешности круга) с центром в нуле. Границы ROC определяются полюсами X(z) - точками, где X(z) уходит в бесконечность.

Обратное Z-преобразование

**Обратное Z-преобразование** восстанавливает x[n] из X(z). Основной метод: **разложение на простые дроби** (partial fractions). Рациональную X(z) раскладываем на сумму простых дробей, для каждой берём обратное Z по таблице.

Шаги: 1) Записать X(z)/z. 2) Разложить на дроби $A/(z-p_1) + B/(z-p_2) + ...$. 3) Умножить на z: $X(z) = Az/(z-p_1) + Bz/(z-p_2) + ...$. 4) Каждое $Az/(z-p) \to A \cdot p^n \cdot u[n]$ (если каузальный, $|z|>|p|$).

**Формальное обратное:** $x[n] = \frac{1}{2\pi j} \oint X(z) \cdot z^{n-1} \, dz$. На практике контурный интеграл не вычисляют - используют partial fractions или вычеты.

При наличии ROC обратное Z однозначно. Без ROC - неоднозначно (каузальный или антикаузальный). Для физически реализуемых (каузальных) систем: ROC - внешность наибольшего полюса.

Передаточная функция

**Передаточная функция** $H(z) = Y(z)/X(z)$ - Z-преобразование импульсной характеристики h[n]. Она полностью описывает LTI-систему. Свёртка $y[n] = x[n] * h[n]$ превращается в умножение: $Y(z) = X(z) \cdot H(z)$. Весь цифровой эквалайзер - это H(z) с нужными полюсами и нулями.

H(z) - рациональная функция: $H(z) = B(z)/A(z)$, где B(z) - полином числителя (нули), A(z) - полином знаменателя (полюса). **Полюса** определяют устойчивость, **нули** - форму частотной характеристики.

**Полюса внутри единичного круга = устойчивость.** Это дискретный аналог «полюса в левой полуплоскости» для непрерывных систем. Связь: $z = e^{sT}$, левая полуплоскость ($\text{Re}(s)<0$) отображается внутрь единичного круга ($|z|<1$).

Вернёмся к свёртке из первого урока: в Z-области свёртка = умножение. Каскад фильтров: $H(z) = H_1(z) \cdot H_2(z)$. Параллельное соединение: $H(z) = H_1(z) + H_2(z)$. Вся архитектура DSP-систем проектируется через полюса и нули передаточной функции.

Частотная характеристика $H(e^{j\omega})$ получается подстановкой $z = e^{j\omega}$ - движение по единичной окружности. Близость к полюсу - усиление частоты (пик), близость к нулю - подавление (провал). Дизайн фильтра - расстановка полюсов и нулей.

Ключевые идеи

• **X(z) = sum x[n]*z^(-n)** - обобщение DTFT на всю комплексную плоскость
• **ROC** определяет, какой сигнал соответствует данному X(z): одна формула - разные сигналы
• **Partial fractions** - основной метод обратного Z: разложить, найти по таблице
• **H(z) = Y(z)/X(z):** полюса внутри единичного круга = устойчивость. Свёртка = умножение в Z-области

Связанные темы

Z-преобразование - аналитический инструмент DSP:

• Сигналы и системы — Z-преобразование формализует свёртку и LTI-системы в алгебраическом виде
• Дискретизация и теорема Котельникова — Связь z = e^(sT): дискретизация в Z-области = периодическое отображение s-плоскости

Вопросы для размышления

• Z-преобразование аналогично Лапласу (z = e^(sT)). Почему устойчивость в s-области: Re(s) < 0, а в z-области: |z| < 1?
• FIR-фильтр (все полюса в z=0) всегда устойчив. Почему тогда используют IIR (с ненулевыми полюсами)? В чём преимущество?
• Можно ли по диаграмме полюсов-нулей H(z) нарисовать приблизительную АЧХ |H(e^jw)|? Как расстояние до полюсов/нулей влияет на амплитуду?

Связанные уроки

• dsp-01 — LTI-системы и свёртка - то, что Z-преобразование превращает в умножение
• dsp-02 — Связь z = e^(sT): дискретизация порождает Z-плоскость из s-плоскости
• dsp-04 — Проектирование фильтров - размещение полюсов и нулей в Z-плоскости
• dsp-05 — IIR vs FIR фильтры: полюса внутри круга = IIR стабилен
• nm-05 — Преобразование Лапласа - непрерывный аналог Z-преобразования
• it-03 — Комплексный анализ: лорановские ряды - математическая основа ROC
• ca-13