在数学和信号处理中,Z变换(英语:Z-transform)把一连串离散的实数或复数信号,从时域转为复频域表示。
可以把它认为是拉普拉斯变换的离散时间等价。在时标微积分中会探索它们的相似性。
历史
现在所知的Z变换的基本思想,拉普拉斯就已了解,而1947年W. Hurewicz用作求解常系数差分方程的一种容易处理的方式。[1] 后来由1952年哥伦比亚大学的采样控制组的雷加基尼和查德称其为“Z变换”。[2][3]
E. I. Jury后来发展并推广了改进或高级Z变换。[4][5]
Z变换中包含的思想在数学里称作母函数方法,该方法可以追溯到1730年的时候,棣莫弗与概率论结合将其引入。[6]
从数学的角度,当把数字序列视为解析函数的(洛朗)展开时,Z变换也可以看成是洛朗级数。
定义
像很多积分变换一样,Z变换可以有单边和双边定义。
双边Z变换
双边Z变换把离散时域信号 x[n] 转为形式幂级数 X(Z)。
![X(z)={\mathcal {Z))\{x[n]\}=\sum _((n=-\infty ))^((\infty ))x[n]z^((-n))](https://wikimedia.org/api/rest_v1/media/math/render/svg/12f6e27003f8c3271124b8af3ea0092c2906ae3e)
当中
是整数,
是复数变量,其表示方式为

其中 A 为 z 的模,j 为虚数单位,而 ɸ 为辐角(也叫相位角),用弧度表示。
单边Z变换
另外,只对 n ≥ 0 定义的 x[n],单边Z变换定义为
![X(z)={\mathcal {Z))\{x[n]\}=\sum _((n=0))^((\infty ))x[n]z^((-n)).](https://wikimedia.org/api/rest_v1/media/math/render/svg/a3e560ddcffcbab6fa176f4d2dd8e3fe60905b55)
在信号处理中,这个定义可以用来计算离散时间因果系统的单位冲激响应。
单边Z变换的一个重要例子是概率母函数,其中 x[n] 部分是离散随机变量取 n 值时的概率,而函数 X(z) 通常写作 X(s),用 s = z−1 表示。Z变换的性质(在下面)在概率论背景下有很多有用的解释。
地球物理学定义
地球物理中的Z变换,通常的定义是 z 的幂级数而非 z−1 的。例如,Robinson、Treitel[7]和Kanasewich都使用这个惯例。[8] 地球物理定义为:
![X(z)={\mathcal {Z))\{x[n]\}=\sum _((n))x[n]z^((n)).](https://wikimedia.org/api/rest_v1/media/math/render/svg/af64bf848f2f92b8aab0469ae4c87827d8092916)
这两个定义是等价的;但差分结果会有一些不同。例如,零点和极点的位置移动在单位圆内使用一个定义,在单位圆外用另一个定义。[7][8]
因此,需要注意特定作者使用的定义。
逆Z变换
逆Z变换为
![x[n]={\mathcal {Z))^((-1))\{X(z)\}={\frac {1}{2\pi j))\oint _((C))X(z)z^((n-1))dz](https://wikimedia.org/api/rest_v1/media/math/render/svg/872e380a9d155a1ee7a3cb5e2ee0e4f033927995)
其中 C 是完全处于收敛域(ROC)内的包围原点的一个逆时针闭合路径。在 ROC 是因果的情况下(参见例2),这意味着路径 C 必须包围 X(z) 的所有极点。
这个曲线积分的一个特殊情形出现在 C 是单位圆的时候(可以在ROC包含单位圆的时候使用,总能保证 X(z) 是稳定的,即所有极点都在单位圆内)。逆Z变换可以化简为逆离散傅里叶变换:
![x[n]={\frac {1}{2\pi ))\int _((-\pi ))^((+\pi ))X(e^((j\omega )))e^((j\omega n))d\omega .](https://wikimedia.org/api/rest_v1/media/math/render/svg/957cad6f61b3feec604ba454617acaea7beae9be)
有限范围 n 和有限数量的均匀间隔的 z 值的Z变换可以用Bluestein的FFT算法方便地计算。离散时间傅里叶变换 (DTFT)—不要与离散傅里叶变换(DFT)混淆—是通过将 z 限制在位于单位圆上而得到的一种Z变换的特殊情况。
收敛域
收敛域(ROC)是指Z变换的求和收敛的复平面上的点集。
![ROC=\left\{z:\left|\sum _((n=-\infty ))^((\infty ))x[n]z^((-n))\right|<\infty \right\}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bad42c66ac900377f174a50adb30ad596be17f9a)
例1(无ROC)
令 x[n] = (0.5)n。在区间 (−∞, ∞) 上展开 x[n] 成为
![x[n]=\left\{\cdots ,0.5^((-3)),0.5^((-2)),0.5^((-1)),1,0.5,0.5^{2},0.5^{3},\cdots \right\}=\left\{\cdots ,2^{3},2^{2},2,1,0.5,0.5^{2},0.5^{3},\cdots \right\}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/5d0a16581c6c01c2dbd61d6e345d0c1daf45a4ef)
观察上面的和
![\sum _((n=-\infty ))^((\infty ))x[n]z^((-n))\to \infty .](https://wikimedia.org/api/rest_v1/media/math/render/svg/e39c1521fe62d231dfb0fae8a8583d4fad0882b0)
因此,没有一个 z 值可以满足这个条件。
例2(因果ROC)
ROC用蓝色表示,单位圆用灰色虚点圆表示(外圈者,而 |z| = 0.5 这个圆用虚线圆表示(内圈者)
令
(其中 u 是单位阶跃函数)。在区间 (−∞, ∞) 上展开 x[n] 得到
![x[n]=\left\{\cdots ,0,0,0,1,0.5,0.5^{2},0.5^{3},\cdots \right\}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/6a7beddbdd74691d956130f78850030ad4d8877e)
观察这个和
![\sum _((n=-\infty ))^((\infty ))x[n]z^((-n))=\sum _((n=0))^((\infty ))0.5^{n}z^((-n))=\sum _((n=0))^((\infty ))\left({\frac {0.5}{z))\right)^{n}={\frac {1}{1-0.5z^((-1)))).](https://wikimedia.org/api/rest_v1/media/math/render/svg/4bcec1d978fd88c533b13fc81a3a7b1dcb784bba)
最后一个等式来自无穷几何级数,而等式仅在 |0.5z−1| < 1 时成立,可以以 z 为变量写成 |z| > 0.5。因此,收敛域为 |z| > 0.5。在这种情况下,收敛域为复平面“挖掉”原点为中心的半径为 0.5 的圆盘。
例3(非因果ROC)
ROC用蓝色表示,单位圆用灰色虚点圆表示(用眼睛看会呈红色),而 |z| = 0.5 这个圆用虚线圆表示
令
(其中 u 是单位阶跃函数)。在区间 (−∞, ∞) 上展开 x[n] 得到
![x[n]=\left\{\cdots ,-(0.5)^((-3)),-(0.5)^((-2)),-(0.5)^((-1)),0,0,0,0,\cdots \right\}.](https://wikimedia.org/api/rest_v1/media/math/render/svg/7bf2b69a4dafb9fc491500100fe2de6c44bfcf08)
观察这个和
![\sum _((n=-\infty ))^((\infty ))x[n]z^((-n))=-\sum _((n=-\infty ))^((-1))0.5^{n}z^((-n))=-\sum _((m=1))^((\infty ))\left({\frac {z}{0.5))\right)^((m))=1-{\frac {1}{1-0.5^((-1))z))={\frac {1}{1-0.5z^((-1))))](https://wikimedia.org/api/rest_v1/media/math/render/svg/a57fe1a23dddbf87776696e5e0073d09629cead3)
再次使用无穷几何级数,此等式只在 |0.5−1z| < 1 时成立,可以用 z 为变量写成 |z| < 0.5。因此,收敛域为 |z| < 0.5。在这种情况下,收敛域为中心在原点的半径为 0.5 的圆盘。
本例与上例的不同之处仅在收敛域上。这是意图展示只有变换结果是不够的。
实例结论
实例2和3清楚地表明,当且仅当指定收敛域时,x[n] 的Z变换 X(z) 才是唯一的。画因果和非因果情形的零极点图表明,在这两种情况下收敛域都不包含极点位于 0.5 的情形。这可以拓展到多个极点的情形:收敛域永远不会包含极点。
在例2中,因果系统产生一个包含 |z| = ∞ 的收敛域,而例3中的非因果系统产生包含 |z| = 0 的收敛域。
ROC表示为蓝色圆环 0.5 < |z| < 0.75
在有多个极点的系统中,收敛域可以既不包含 |z| = ∞ 也不包含 |z| = 0。画出的收敛域与一个圆形带。例如,
![x[n]=0.5^{n}u[n]-0.75^{n}u[-n-1]](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a35aa48e0a80015443d04e9c2af649eb8979eab)
的极点为 0.5 与 0.75。收敛域会是 0.5 < |z| < 0.75,不包含原点和无穷大。这样的系统称为混合因果系统,因为它包含一个因果项 (0.5)nu[n] 和一个非因果项 −(0.75)nu[−n−1]。
一个系统的稳定性可以只通过了解收敛域来确定。如果收敛域包含单位圆(即 |z| = 1),那么系统是稳定的。在上述系统中因果系统(例2)是稳定的,因为 |z| > 0.5 包含单位圆。
如果给定一个没有收敛域的Z变换(即模糊的 x[n]),可以确定一个唯一的 x[n] 满足下列:
如果你要稳定性,收敛域必须包含单位圆;如果你需要一个因果系统,收敛域必须包含无穷大,并且系统函数应为一个右边序列。如果你需要一个非因果系统,那么收敛域必须包含原点,且系统函数为左边序列。如果你既要稳定性,也要因果性,系统函数的所有极点都必须在单位圆内。
可以找到唯一的 x[n]。
性质
Z变换性质
|
时域
|
Z域
|
证明
|
收敛域
|
记法
|
|
|
|
|
线性
|
|
|
|
包含 ROC1 ∩ ROC2
|
时间膨胀
|
r: 整数
|
|
|
|
降采样
|
|
|
ohio-state.edu 或 ee.ic.ac.uk
|
|
时移
|
|
|
|
ROC,除了 k > 0 时 z = 0 和 k < 0 时 z = ∞
|
Z域的
尺度性质
|
|
|
|
|
时间反转
|
|
|
|
|
共轭复数
|
|
|
|
|
实部
|
|
|
|
|
虚部
|
|
|
|
|
微分
|
|
|
|
|
卷积
|
|
|
|
包含 ROC1 ∩ ROC2
|
互相关
|
|
|
|
包含 与 的ROC的交集
|
一阶差分
|
|
|
|
包含 X1(z) 与 z ≠ 0 的ROC的交集
|
累积
|
|
|
|
|
乘法
|
|
|
|
-
|
帕塞瓦尔定理
![\sum _((n=-\infty ))^((\infty ))x_{1}[n]x_{2}^{*}[n]\quad =\quad {\frac {1}{j2\pi ))\oint _{C}X_{1}(v)X_{2}^{*}({\tfrac {1}{v^{*))})v^((-1)){\mathrm {d))v](https://wikimedia.org/api/rest_v1/media/math/render/svg/ff45b737972bd6dc88fc06588ae6e08910d74e8b)
初值定理:如果 x[n] 为因果的,那么
![x[0]=\lim _((z\to \infty ))X(z).](https://wikimedia.org/api/rest_v1/media/math/render/svg/815494d8476445adef605f74b5b5a6765fb203c4)
终值定理:如果 (z−1)X(z) 的极点在单位圆内,则
![x[\infty ]=\lim _((z\to 1))(z-1)X(z).](https://wikimedia.org/api/rest_v1/media/math/render/svg/1462160ef0f0d8de8000f78372f2a1b21c2a6031)
常见的Z变换对表
这里:
![u:n\mapsto u[n]={\begin{cases}1,&n\geq 0\\0,&n<0\end{cases))](https://wikimedia.org/api/rest_v1/media/math/render/svg/08c15373dbd58410ab17d3c2c7ebe2123e276298)
是单位阶跃函数而
![\delta :n\mapsto \delta [n]={\begin{cases}1,&n=0\\0,&n\neq 0\end{cases))](https://wikimedia.org/api/rest_v1/media/math/render/svg/37c07b1365b8d2566271e4752535678f92aedafa)
是离散时间单位冲激函数。两者通常都不认为是真正的函数,但由于它们的不连续性把它们看成是分布(它们在 n = 0 处的值通常无关紧要,除非在处理离散时间的时候,它们会变成衰减离散级数;在本章节中对连续和离散时间域,都在 n = 0 处取 1,否则不能使用下表中收敛域一栏的内容)。同时列出两个“函数”,使得(在连续时间域)单位阶跃函数是单位冲激函数的积分,或(在离散时间域)单位阶跃函数是单位冲激函数的求和,因此要令他们的值在 n = 0 处为 1。
|
信号,![x[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/864cbbefbdcb55af4d9390911de1bf70167c4a3d) |
Z变换, |
ROC
|
1 |
![\delta [n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/f2a6caf535cb44fa3526b2f320330a805edfdfaa) |
1 |
所有 z
|
2 |
![\delta [n-n_{0}]](https://wikimedia.org/api/rest_v1/media/math/render/svg/2bdb0265027e056f16fce87ab282b57cb03c4f8c) |
 |
|
3 |
![u[n] \,](https://wikimedia.org/api/rest_v1/media/math/render/svg/e693a2911b29e6c8d440d97e46d27760559af7c5) |
 |
|
4 |
![e^((-\alpha n))u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/4074d2265c2c901f6c15c3fbcbbf54d3a420d2eb) |
 |
|
5 |
![-u[-n-1]](https://wikimedia.org/api/rest_v1/media/math/render/svg/34dd7dba0f892e5bcad792136d96cd5f5548a327) |
 |
|
6 |
![nu[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8a28e84b105a96db578fb6e6b047465468b77ec) |
 |
|
7 |
![- n u[-n-1] \,](https://wikimedia.org/api/rest_v1/media/math/render/svg/b9c7bfd00539cf805ba91e15a60b73576194dbd1) |
 |
|
8 |
![n^{2}u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/ea874c2bd6b83f29b93caf0cbe50ee9131eaebc2) |
 |
|
9 |
![- n^2 u[-n - 1] \,](https://wikimedia.org/api/rest_v1/media/math/render/svg/bcc9d247970a92c7c6a69da9b5a272190dadcd24) |
 |
|
10 |
![n^{3}u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/e9e2a53a00fc122eed75716c0c58cf9e58a0f38d) |
 |
|
11 |
![-n^{3}u[-n-1]](https://wikimedia.org/api/rest_v1/media/math/render/svg/7bd111768ef860fc18a2c93e5dc2fb4b03dfab8c) |
 |
|
12 |
![a^{n}u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/bef62e50254aa3175939a01611766c01f9bf7b39) |
 |
|
13 |
![-a^{n}u[-n-1]](https://wikimedia.org/api/rest_v1/media/math/render/svg/4718b1c4477718ebeb49ac1fc41415cadeadf1e7) |
 |
|
14 |
![na^{n}u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5885cf352282908bc931ed56ad572fa84f6235c) |
 |
|
15 |
![-na^{n}u[-n-1]](https://wikimedia.org/api/rest_v1/media/math/render/svg/34d2fae6bc70beb0ec9d5881b38a29d427823fad) |
 |
|
16 |
![n^{2}a^{n}u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/8b68eda406d1e088553723c0395d4ce2cdeff46e) |
 |
|
17 |
![-n^{2}a^{n}u[-n-1]](https://wikimedia.org/api/rest_v1/media/math/render/svg/97bf891898f79a3cd0cf05030244592b6aaad421) |
 |
|
18 |
![\cos(\omega _{0}n)u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/6911a3c468c99d1dc042b3b5015b48108d9476aa) |
 |
|
19 |
![\sin(\omega _{0}n)u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/f59eccb10aa95ef5ba0a1ed904aee27526fe377d) |
 |
|
20 |
![a^{n}\cos(\omega _{0}n)u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b90c4e4b46e7725d99960e3f99a846c65a5d5da) |
 |
|
21 |
![a^{n}\sin(\omega _{0}n)u[n]](https://wikimedia.org/api/rest_v1/media/math/render/svg/4663af1f68929f4e26833381893076c001dfbebb) |
 |
|
与傅里叶级数和傅里叶变换的关系
对于区域 |z|=1(称为单位圆)内的 z 值,我们可以通过定义 z=ejω 来用单一实变量的函数来表示该变换。于是双边变换就简化为了傅里叶级数:
-
![\sum _((n=-\infty ))^((\infty ))x[n]\ z^((-n))=\sum _((n=-\infty ))^((\infty ))x[n]\ e^((-j\omega n)),](https://wikimedia.org/api/rest_v1/media/math/render/svg/f14ea7fbaab1d911450572d089188f9ed49bd4d1)
|
|
(Eq.1)
|
也被称作 x[n] 序列的离散时间傅里叶变换(DTFT)。这个以 2π 为周期的函数是傅里叶变换的周期性求和,这使得它成为广泛使用的分析工具。要理解这一点,令 X(f) 为任意函数 x(t) 的傅里叶变换,该函数以某个间隔 T 采样就与 x[n] 序列相等。于是 x[n] 序列的DTFT可以写作:
![\underbrace {\sum _((n=-\infty ))^((\infty ))\overbrace {x(nT)}^((x[n]))\ e^((-j2\pi fnT))}_(({\text{DTFT))))={\frac {1}{T))\sum _((k=-\infty ))^((\infty ))X(f-k/T).](https://wikimedia.org/api/rest_v1/media/math/render/svg/a954eea32605d9f5457cf746677b7d04bc788ef6)
若T的单位是秒,
的单位即为赫兹。比较两个数列可得
为标准化频率,单位是radians per sample。数值ω=2π对应
Hz. ,而且在替换
后, Eq.1可以表示为傅里叶变换X(•):
![\sum _((n=-\infty ))^((\infty ))x[n]\ e^((-j\omega n))={\frac {1}{T))\sum _((k=-\infty ))^((\infty ))\underbrace {X\left({\tfrac {\omega }{2\pi T))-{\tfrac {k}{T))\right)}_((X\left({\frac {\omega -2\pi k}{2\pi T))\right))).](https://wikimedia.org/api/rest_v1/media/math/render/svg/3539536a288102ecbb30aaf9c3aedb72241a0df2)
若数列x(nT)表示线性时不变系统的冲激响应,这些函数也称为频率响应,当x(nT)是周期性数列,其DTFT在一或多个共振频率发散,在其他频率均为零。这一般会用在共振频率,振幅可变的狄拉克δ函数表示。因为其周期性,只会有有限个振幅,可以用较简单许多的离散傅里叶变换来计算。(参照离散傅里叶变换#周期性)
和拉氏变换的关系
双线性变换
双线性变换可以用在连续时间滤波器(用拉氏域表示)和离散时间滤波器(用Z域表示)之间的变换,其变换关系如下:

将一个拉氏域的函数
变换为Z域下的
,或是

从Z域变换到拉氏域。借由双线性变换,复数的s平面(拉氏变换)可以映射到复数的z平面(Z变换)。这个变换是非线性的,可以将S平面的整个jΩ轴映射到Z平面的单位圆内。因此,傅里叶变换(在jΩ axis计算的拉氏变换)变成离散时间傅里叶变换,前提是假设其傅里叶变换存在,也就是拉氏变换的收敛区域包括jΩ轴。
线性常系数差分方程
线性常系数差分(LCCD)方程是基于自回归滑动平均的线性系统表达形式。
![\sum _((p=0))^((N))y[n-p]\alpha _((p))=\sum _((q=0))^((M))x[n-q]\beta _((q))](https://wikimedia.org/api/rest_v1/media/math/render/svg/f003f262cbe68dc97f84a5dac9e927a3f2c52858)
上面等式两边可以同时除以 α0,如果非零,正规化 α0 = 1,LCCD方程可以写成
![y[n]=\sum _((q=0))^((M))x[n-q]\beta _((q))-\sum _((p=1))^((N))y[n-p]\alpha _((p)).](https://wikimedia.org/api/rest_v1/media/math/render/svg/fea05d580f1b1fecb9f9aa3ec3d80a9e0d8b02b1)
LCCD方程的这种形式有利于更加明确“当前”输出 y[n] 是过去输出 y[n−p]、当前输入 x[n] 与之前输入 x[n−q] 的一个函数。
传递函数
对上述方程去Z变换(使用线性和时移法则)得到

整理结果

零点和极点
由代数基本定理得知分子有 M 个根(对应于 H 的零点)和分母有 N 个根(对应于极点)。用极点和零点重新整理传递函数为

其中 qk 为 k 阶零点,pk 为 k 阶极点。零点和极点通常是复数,当在复平面(z平面)作图时称为零极点图。
此外,在 z = 0 和 z = ∞ 也可能存在零点和极点。如果我们把这些极点和零点以及高阶零点和极点考虑在内的话,零点和极点的数目总会相等。
通过对分母因式分解,可以使用部分分式分解可以转换回时域。这样做会导出系统的冲激响应和线性常系数差分方程。
输出响应
如果一个系统 H(z) 由信号 X(z) 驱动,那么输出为 Y(z) = H(z)X(z)。通过对 Y(z) 部分分式分解并取逆Z变换可以得到输出 y[n]。在实际运用中,在分式分解
之后再乘 z 产生 Y(z) 的一个形式(含有很容易计算逆Z变换的项)往往很有用。
延伸阅读
- Refaat El Attar, Lecture notes on Z-Transform, Lulu Press, Morrisville NC, 2005. ISBN 978-1-4116-1979-1.
- Ogata, Katsuhiko, Discrete Time Control Systems 2nd Ed, Prentice-Hall Inc, 1995, 1987. ISBN 978-0-13-034281-2.
- Alan V. Oppenheim and Ronald W. Schafer (1999). Discrete-Time Signal Processing, 2nd Edition, Prentice Hall Signal Processing Series. ISBN 978-0-13-754920-7.
外部链接
|
---|
| 理论 | |
---|
| 子领域 | |
---|
| 技术 | |
---|
| 采样 |
- 混叠
- 奈奎斯特率 / 频率
- 升采样
- 降采样
- 过采样
- 欠采样
- 采样率
- 量化
|
---|
|