Основная идея, известная теперь как Z-преобразование, была известна Лапласу и была повторно представлена в 1947 году У. Гуревичем [1] [2] и другими как способ обработки систем управления выборками данных, используемых с радаром. Это дает удобный способ решения линейных разностных уравнений с постоянным коэффициентом . Позже это было названо «z-преобразованием» Рагаццини и Заде в контрольной группе выборочных данных в Колумбийском университете в 1952 г. [3] [4]
Идея, содержащаяся в Z-преобразовании, также известна в математической литературе как метод производящих функций, который можно проследить еще в 1730 году, когда он был введен де Муавром вместе с теорией вероятностей. [7]
С математической точки зрения Z-преобразование также можно рассматривать как ряд Лорана, в котором рассматриваемая последовательность чисел рассматривается как разложение (Лорана) аналитической функции.
Определение [ править ]
Z-преобразование может быть определено как одностороннее или двустороннее преобразование. [8]
Двустороннее Z-преобразование [ править ]
Двусторонний или двусторонний Z-преобразование дискретного сигнала времени является формальным степенным рядом , определенный как
( Уравнение 1 )
где - целое число и , как правило, комплексное число :
где - величина , - мнимая единица , - комплексный аргумент (также называемый углом или фазой ) в радианах .
Одностороннее Z-преобразование [ править ]
В качестве альтернативы, в случаях, когда определено только для , одностороннее или одностороннее Z-преобразование определяется как
( Уравнение 2 )
При обработке сигналов это определение может использоваться для оценки Z-преобразования единичной импульсной характеристики причинной системы с дискретным временем .
Важным примером одностороннего Z-преобразования является функция , генерирующая вероятность , где компонент - это вероятность того, что дискретная случайная величина примет значение , а функция обычно записывается как в терминах . Свойства Z-преобразований (см. Ниже) имеют полезные интерпретации в контексте теории вероятностей.
Обратное Z-преобразование [ править ]
Обратное Z-преобразование
( Уравнение 3 )
где C - замкнутый путь против часовой стрелки, охватывающий начало координат и полностью находящийся в области конвергенции (ROC). В случае, когда ROC является причинным (см. Пример 2 ), это означает, что путь C должен охватывать все полюса .
Частный случай этого контурного интеграла возникает, когда C - единичная окружность. Этот контур можно использовать, когда ROC включает единичную окружность, что всегда гарантируется, когда она устойчива, то есть когда все полюса находятся внутри единичной окружности. С этим контуром обратное Z-преобразование упрощается до обратного дискретного преобразования Фурье или ряда Фурье периодических значений Z-преобразования по единичной окружности:
Z-преобразование с конечным диапазоном значений n и конечным числом равномерно распределенных значений z может быть эффективно вычислено с помощью алгоритма БПФ Блустейна . Дискретным временем преобразования Фурье (ДВПФ) -не следует путать с дискретным преобразованием Фурье (ДПФ) -это частный случай такой Z-преобразования , получаемого путем ограничения г лежать на единичной окружности.
Область конвергенции [ править ]
Область конвергенции (ROC) - это набор точек на комплексной плоскости, для которых суммирование Z-преобразования сходится.
Пример 1 (без ROC) [ править ]
Пусть x [n] = (0,5) n . Разложив x [n] на интервал (−∞, ∞), получим
Глядя на сумму
Следовательно, не существует значений z , удовлетворяющих этому условию.
Пример 2 (причинный ROC) [ править ]
ROC показан синим цветом, единичный круг в виде серого пунктирного круга и круг | z | = 0,5 отображается пунктирным черным кружком.
Пусть (где u - ступенчатая функция Хевисайда ). Разложив x [n] на интервал (−∞, ∞), получим
Глядя на сумму
Последнее равенство возникает из бесконечного геометрического ряда, и равенство выполняется только в том случае, если | 0.5 z −1 | <1, который можно переписать через z как | z | > 0,5. Таким образом, ОКР | z | > 0,5. В этом случае ROC представляет собой комплексную плоскость с «пробитым» диском радиуса 0,5 в начале координат.
Пример 3 (антикаузальная ОКР) [ править ]
ROC показан синим цветом, единичный круг в виде серого пунктирного круга и круг | z | = 0,5 отображается пунктирным черным кружком.
Пусть (где u - ступенчатая функция Хевисайда ). Разложив x [n] на интервал (−∞, ∞), получим
Глядя на сумму
Опять же, используя бесконечный геометрический ряд , равенство выполняется только в том случае, если | 0.5 −1 z | <1, который можно переписать через z как | z | <0,5. Таким образом, ОКР | z | <0,5. В этом случае ROC представляет собой диск с центром в начале координат и радиусом 0,5.
Что отличает этот пример от предыдущего, так это только ROC. Это сделано намеренно, чтобы продемонстрировать, что одного результата преобразования недостаточно.
Вывод из примеров [ править ]
Примеры 2 и 3 ясно показывают, что Z-преобразование X (z) для x [n] уникально тогда и только тогда, когда задается ROC. Создание графика полюс-ноль для причинно-следственного и антикаузального случая показывает, что ROC для любого случая не включает полюс, находящийся на 0,5. Это распространяется на случаи с несколькими полюсами: ROC никогда не будет содержать полюсов.
В примере 2 причинная система дает ROC, который включает | z | = ∞, в то время как антикаузальная система в примере 3 дает ROC, который включает | z | = 0.
ROC показан синим кольцом 0.5 <| z | <0,75
В системах с несколькими полюсами возможно наличие ROC, не содержащего ни | z | = ∞ ни | z | = 0. ROC создает круговую полосу. Например,
имеет полюса 0,5 и 0,75. ROC будет 0,5 <| z | <0,75, что не включает ни начало координат, ни бесконечность. Такая система называется системой смешанной причинности, поскольку она содержит причинный член (0.5) n u [ n ] и антикаузальный член - (0.75) n u [- n −1].
Стабильность системы также можно определить, зная РПЦ в одиночку. Если ROC содержит единичную окружность (т. Е. | Z | = 1), то система устойчива. В приведенных выше системах причинная система (Пример 2) устойчива, поскольку | z | > 0,5 содержит единичный круг.
Предположим, нам предоставлено Z-преобразование системы без ROC (т.е. неоднозначное x [n] ). Мы можем определить уникальный x [n], если нам нужно следующее:
Стабильность
Причинно-следственная связь
Для устойчивости ROC должен содержать единичный круг. Если нам нужна причинная система, тогда ROC должен содержать бесконечность, а функция системы будет правосторонней последовательностью. Если нам нужна антикаузальная система, тогда ROC должен содержать начало координат, а функция системы будет левосторонней последовательностью. Если нам нужна как стабильность, так и причинность, все полюса системной функции должны находиться внутри единичного круга.
Затем можно найти уникальный x [n] .
Свойства [ править ]
Свойства z-преобразования
Область времени
Z-домен
Доказательство
ROC
Обозначение
Линейность
Содержит ROC 1 ∩ ROC 2
Расширение времени
с
Децимация
ohio-state.edu или ee.ic.ac.uk
Временная задержка
с и
ROC, кроме z = 0, если k > 0, и z = ∞, если k <0
Время вперед
с
Двустороннее Z-преобразование:
Одностороннее Z-преобразование: [9]
Первое отличие назад
с x [ n ] = 0 для n <0
Содержит пересечение ROC X 1 (z) и z ≠ 0
Первая разница вперед
Обратное время
Масштабирование в z-области
Комплексное сопряжение
Реальная часть
Мнимая часть
Дифференциация
ОКР, если рационально;
ROC возможно без границы, если иррационально [10]
Свертка
Содержит ROC 1 ∩ ROC 2
Взаимная корреляция
Содержит пересечение ОКР и
Накопление
Умножение
-
Теорема Парсеваля
Теорема о начальном значении : если x [ n ] причинно, то
Теорема об окончательном значении : если полюса ( z −1) X ( z ) находятся внутри единичной окружности, то
Таблица общих пар Z-преобразований [ править ]
Здесь:
- ступенчатая функция единицы (или Хевисайда) и
представляет собой единичную импульсную функцию с дискретным временем (см. дельта-функцию Дирака, которая является версией с непрерывным временем). Две функции выбираются вместе, так что функция единичного шага представляет собой накопление (промежуточный итог) единичной импульсной функции.
Сигнал,
Z-преобразование,
ROC
1
1
все z
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
, для положительного целого числа [10]
18
, для положительного целого числа [10]
19
20
21 год
22
Связь с рядами Фурье и преобразованием Фурье [ править ]
Основная статья: Дискретное преобразование Фурье
Для значений в области , известной как единичный круг , мы можем выразить преобразование как функцию единственной действительной переменной ω путем определения . А двустороннее преобразование сводится к ряду Фурье :
( Уравнение 4 )
которое также известно как преобразование Фурье с дискретным временем (ДВПФ) последовательности. Этот 2 π -периодической функция является периодическим суммированием из преобразования Фурье , что делает его широко используемым инструментом анализа. Чтобы понять это, позвольте быть преобразованием Фурье любой функции, выборки которой на некотором интервале T равны последовательности x [ n ]. Тогда ДВПФ последовательности x [ n ] можно записать следующим образом.
( Уравнение 5 )
Когда T имеет единицы секунд, имеет единицы герц . Сравнение двух серий показывает, что это нормализованная частота с единицами радиан на выборку . Значение ω = 2 π соответствует Гц. И теперь, с заменой Eq.4 может быть выражено через преобразование Фурье, X (•) :
( Уравнение 6 )
При изменении параметра T отдельные члены уравнения (5) перемещаются дальше друг от друга или сближаются вдоль оси f. Однако в уравнении 6 центры остаются на расстоянии 2 π друг от друга, а их ширина расширяется или сжимается. Когда последовательность x ( nT ) представляет собой импульсную характеристику системы LTI , эти функции также известны как ее частотная характеристика . Когда последовательность является периодической, ее ДВПФ расходится на одной или нескольких гармонических частотах и ноль на всех остальных частотах. Это часто представлено использованием амплитудно- вариативной дельты Дирака.функции на частотах гармоник. Из-за периодичности существует только конечное число уникальных амплитуд, которые легко вычисляются с помощью гораздо более простого дискретного преобразования Фурье (ДПФ). (См. DTFT § Периодические данные .)
Связь с преобразованием Лапласа [ править ]
Билинейное преобразование [ править ]
Основная статья: Билинейное преобразование
Билинейное преобразование может быть использованы для преобразования непрерывного времени фильтров (представленные в области Лапласа) в дискретное время фильтров (представленных в Z-домене), и наоборот. Используется следующая замена:
преобразовать некоторую функцию в области Лапласа в функцию в Z-области ( преобразование Тастина ), или
из Z-области в область Лапласа. Посредством билинейного преобразования комплексная s-плоскость (преобразования Лапласа) отображается в комплексную z-плоскость (z-преобразования). Хотя это отображение (обязательно) является нелинейным, оно полезно тем, что отображает всю ось s-плоскости на единичный круг в z-плоскости. Таким образом, преобразование Фурье (которое является преобразованием Лапласа, вычисляемым на оси) становится преобразованием Фурье с дискретным временем. Это предполагает, что преобразование Фурье существует; то есть ось находится в области сходимости преобразования Лапласа.
Помеченное преобразование [ править ]
Основная статья: Преобразование с пометкой
Учитывая одностороннее Z-преобразование, X (z), функции с временной дискретизацией, соответствующее преобразование, помеченное звездочкой, производит преобразование Лапласа и восстанавливает зависимость от параметра дискретизации, T :
Обратное преобразование Лапласа - это математическая абстракция, известная как функция с импульсной выборкой .
Линейное уравнение разности с постоянным коэффициентом [ править ]
Уравнение линейной разности постоянных коэффициентов (LCCD) представляет собой представление линейной системы, основанной на уравнении авторегрессионного скользящего среднего .
Обе части приведенного выше уравнения можно разделить на α 0 , если оно не равно нулю, нормализуя α 0 = 1, и уравнение LCCD можно записать
Эта форма уравнения LCCD позволяет сделать более явным, что «текущий» выход y [n] является функцией прошлых выходов y [n − p] , текущего входа x [n] и предыдущих входов x [n− q] .
Передаточная функция [ править ]
Принятие Z-преобразования приведенного выше уравнения (с использованием законов линейности и сдвига во времени) дает
и переставляя результаты в
Нули и полюсы [ править ]
Из основной теоремы алгебры числителя имеет M корни ( что соответствует нулям H) , а знаменатель имеет N корни ( что соответствует полюсам). Переписываем передаточную функцию в терминах нулей и полюсов
где q k - k -й нуль, а p k - k-й полюс. Нули и полюсы обычно являются сложными, и когда они нанесены на комплексную плоскость (z-плоскость), это называется графиком полюс – ноль .
Кроме того, могут существовать нули и полюсы в точках z = 0 и z = ∞. Если мы примем во внимание эти полюса и нули, а также нули и полюсы нескольких порядков, количество нулей и полюсов всегда будет одинаковым.
Разложив знаменатель на множители, можно использовать частичное дробное разложение, которое затем можно преобразовать обратно во временную область. Это приведет к появлению импульсной характеристики и линейного уравнения разности постоянных коэффициентов системы.
Выходной ответ [ править ]
Если такая система H (z) управляется сигналом X (z), то на выходе будет Y (z) = H (z) X (z) . Выполняя частичное дробное разложение на Y (z) и затем выполняя обратное Z-преобразование, можно найти выход y [n] . На практике часто бывает полезно дробно разложить перед умножением этой величины на z, чтобы сгенерировать форму Y (z), которая имеет члены с легко вычисляемыми обратными Z-преобразованиями.
См. Также [ править ]
Расширенное Z-преобразование
Билинейное преобразование
Разностное уравнение (рекуррентное соотношение)
Дискретная свертка
Дискретное преобразование Фурье
Конечный импульсный отклик
Формальный степенной ряд
Производящая функция
Преобразование производящей функции
Преобразование Лапласа
Серия Laurent
Вероятностно-производящая функция
Звездное преобразование
Зак преобразовать
Регуляризация дзета-функции
Ссылки [ править ]
^ ER Kanasewich (1981). Анализ временных последовательностей в геофизике . Университет Альберты. С. 186, 249. ISBN 978-0-88864-074-1.
^ ER Kanasewich (1981). Анализ временной последовательности в геофизике (3-е изд.). Университет Альберты. С. 185–186. ISBN 978-0-88864-074-1.
^ Рагаццини, младший; Заде, Л.А. (1952). «Анализ систем выборки данных». Труды Американского института инженеров-электриков, Часть II: Приложения и промышленность . 71 (5): 225–234. DOI : 10,1109 / TAI.1952.6371274 . S2CID 51674188 .
^ Корнелиус Т. Леондес (1996). Реализация цифровых систем управления и вычислительная техника . Академическая пресса. п. 123. ISBN 978-0-12-012779-5.
^ Элиягу Ibrahim Юрий (1958). Системы контроля выборочных данных . Джон Вили и сыновья.
^ Элиягу Ibrahim Юрий (1973). Теория и применение метода Z-преобразования . Krieger Pub Co. ISBN 0-88275-122-0.
^ Элиягу Ибрагим Юрий (1964). Теория и применение метода Z-преобразования . Джон Вили и сыновья. п. 1.
^ Точно так же, как у нас есть одностороннее преобразование Лапласа и двустороннее преобразование Лапласа .
^ Больцерн, Паоло; Скаттолини, Риккардо; Скьявони, Никола (2015). Fondamenti di Controlli Automatici (на итальянском языке). MC Graw Hill Education. ISBN 978-88-386-6882-1.
^ а б в А. Р. Форузан (2016). «Область сходимости производной Z преобразования». Письма об электронике . 52 (8): 617–619. DOI : 10.1049 / el.2016.0189 .
Дальнейшее чтение [ править ]
Рефаат Эль Аттар, Конспекты лекций по Z-преобразованию , Lulu Press, Моррисвилл, Северная Каролина, 2005. ISBN 1-4116-1979-X .
Огата, Кацухико, Системы управления дискретным временем, 2-е изд. , Prentice-Hall Inc, 1995, 1987. ISBN 0-13-034281-5 .
Алан В. Оппенгейм и Рональд В. Шафер (1999). Обработка сигналов в дискретном времени, 2-е издание, Серия Prentice Hall Signal Processing. ISBN 0-13-754920-2 .
Внешние ссылки [ править ]
"Z-преобразование" , Энциклопедия математики , EMS Press , 2001 [1994]
Численное обращение Z-преобразования
Таблица Z-преобразований некоторых распространенных преобразований Лапласа
Запись Mathworld о Z-преобразовании
Потоки Z-Transform в Comp.DSP
График взаимосвязи между s-плоскостью преобразования Лапласа и Z-плоскостью преобразования Z