如何求z变换

本文旨在为读者系统性地阐述如何求解z变换。文章将从基本定义出发，详细解析直接计算法、幂级数展开法以及利用常见变换对与性质进行求解的核心方法。内容涵盖从简单序列到复杂有理分式函数的处理技巧，包括部分分式分解与留数定理的应用，并对系统函数与收敛域的确定进行深入探讨。通过结合理论与实例，力求为工程与科研人员提供一份清晰、实用且具备深度的z变换求解指南。

       在数字信号处理与离散时间系统分析领域，z变换扮演着极为关键的角色。它不仅是分析线性时不变系统（线性时不变系统， Linear Time-Invariant System）稳定性与频率响应的有力工具，更是连接离散序列与复频域表示的核心桥梁。掌握如何求解一个给定序列的z变换，是深入理解并应用这一理论的基础。本文将避开过于抽象的纯理论推导，聚焦于一系列实用且系统的求解方法与技巧，引导读者逐步构建起求解z变换的完整知识框架。

一、 理解z变换的核心定义与收敛域

       求解z变换的第一步，是准确理解其数学定义。对于一个离散时间序列x[n]，其z变换X(z)定义为双边无穷级数：X(z) = Σ_n=-∞^∞ x[n] z^-n。这里的z是一个复变量。这个定义本身即提供了一种最基础的求解方法——直接求和法。但更重要的是，我们必须意识到，并非所有z值都能使该级数收敛。使级数收敛的所有z值的集合，称为z变换的收敛域（收敛域， Region of Convergence）。收敛域通常是一个以原点为中心的圆环区域，其形状由序列x[n]的性质决定。在求解过程中，明确标注收敛域与求得变换式本身同等重要，因为相同的z变换表达式配以不同的收敛域，可能对应着完全不同的时域序列。

二、 从简单序列入手：直接计算法

       对于形式简单的有限长序列或典型的无限长序列，直接根据定义进行求和是最直观的方法。例如，求解单位脉冲序列δ[n]的z变换。根据定义，X(z) = Σ_n=-∞^∞ δ[n] z^-n。由于δ[n]仅在n=0时为1，其余时刻为0，因此求和结果即为1，且收敛域为整个z平面。再比如，求解右边指数序列a^n u[n]（其中u[n]是单位阶跃序列）的z变换。将其代入定义式，得到X(z) = Σ_n=0^∞ (a^n) z^-n = Σ_n=0^∞ (a z^-1)^n。这是一个几何级数，当|a z^-1| < 1，即|z| > |a|时级数收敛，其和为1 / (1 - a z^-1)。因此，变换对为：a^n u[n] ←→ 1/(1 - a z^-1)，收敛域为|z| > |a|。通过直接计算这些基本序列，我们可以积累最基础的变换对库。

三、 利用已知变换对与线性性质

       z变换拥有一系列重要的性质，其中线性性质是最常用且强大的工具之一。该性质表明，若序列x1[n]和x2[n]的z变换分别为X1(z)和X2(z)，且收敛域分别为R1和R2，则对于任意常数a和b，序列a x1[n] + b x2[n]的z变换为a X1(z) + b X2(z)，其收敛域至少是R1与R2的交集。这意味着，对于一个复杂序列，若能将其分解为若干简单序列的线性组合，而这些简单序列的z变换是已知的，那么我们就可以直接写出其变换式。例如，求解序列x[n] = 2 δ[n] - 3 (0.5)^n u[n]的z变换。利用线性性质，分别求δ[n]和(0.5)^n u[n]的变换（后者可利用上一节的结果，将a替换为0.5），然后进行线性组合即可得到结果。

四、 时移性质的应用技巧

       时移性质是处理序列在时间轴上平移的关键。对于右边序列x[n] u[n]（其变换为X(z)），则有x[n-k] u[n-k] ←→ z^-k X(z)，其中k为正整数，代表向右延迟k个单位。这一性质在求解差分方程或分析系统响应时尤为重要。例如，已知单位阶跃序列u[n]的z变换为1/(1 - z^-1) (|z|>1)，那么延迟后的序列u[n-3]的z变换即为z^-3 / (1 - z^-1)。需要注意的是，时移性质的应用必须谨慎考虑序列的起始点，对于非因果序列或更一般的情况，需根据定义重新审视。

五、 指数加权序列的变换求解

       指数加权性质，也称为z域尺度变换性质，表述为：若x[n] ←→ X(z)，收敛域为R，则对于任意常数a（实数或复数），有a^n x[n] ←→ X(z/a)，且收敛域变为|a|R。这个性质极大地扩展了变换对库。例如，我们已经知道u[n] ←→ 1/(1 - z^-1)。那么，根据此性质，序列a^n u[n]的z变换即为将原变换X(z)中的z替换为z/a，得到1 / (1 - (z/a)^-1) = 1 / (1 - a z^-1)，这与我们通过直接计算得到的结果完全一致。该性质使得我们可以从一组基础变换，快速推导出更多相关序列的变换。

六、 处理初值不为零的右边序列

       在实际问题中，序列并非总是从n=0开始。对于形如x[n] u[n+n0]（n0为正整数）的序列，即序列在n=-n0时已开始，直接应用标准右边序列的变换公式或性质可能会出错。此时，可靠的方法是回归定义进行求和。例如，求解序列x[n] = b^n u[n+2]的z变换。根据定义，X(z) = Σ_n=-2^∞ b^n z^-n。我们可以令m = n+2，则n = m-2，代入后求和指标m从0到∞，从而将求和式转化为已知的几何级数形式进行计算。这种方法虽然步骤稍多，但能确保结果的准确性。

七、 求解左边序列的z变换

       左边序列，即当n大于某个值时序列值全为零的序列，其z变换的求解需要特别注意收敛域。考虑一个典型的左边指数序列：-a^n u[-n-1]。根据定义，其z变换为X(z) = Σ_n=-∞^-1 (-a^n) z^-n。令m = -n，则求和变为Σ_m=1^∞ (-a^-m) z^m = - Σ_m=1^∞ (a^-1 z)^m。这同样是一个几何级数，当|a^-1 z| < 1，即|z| < |a|时收敛，其和为 - (a^-1 z) / (1 - a^-1 z) = 1 / (1 - a z^-1)。有趣的是，其表达式与右边序列a^n u[n]的变换式完全相同，但收敛域截然不同（|z| < |a|）。这再次强调了收敛域在z变换中的决定性作用。

八、 幂级数展开法及其适用场景

       对于某些非有理函数形式的z变换，或者当我们需要从已知的X(z)反求其对应的序列x[n]时，幂级数展开法（也称为长除法）是一种实用工具。该方法的核心是将X(z)在特定的收敛域内展开为z的正幂或负幂的级数形式，级数中z^-n项的系数即为序列值x[n]。例如，若X(z) = 1 / (1 - 1.5z^-1 + 0.5z^-2)，且已知收敛域为|z|>1，表明对应一个右边序列。我们可以通过多项式长除法，将X(z)展开为z^-1的幂级数：1 + 1.5z^-1 + 1.75z^-2 + ...，从而直接读出x[0]=1, x[1]=1.5, x[2]=1.75, ...。这种方法直观，尤其适用于求取序列的前若干项。

九、 部分分式分解法：处理有理分式的利器

       在工程实践中，绝大多数系统的传递函数，即序列z变换的表达式X(z)，都是关于z^-1的有理分式。此时，部分分式分解法成为求解其逆变换（或理解其构成）的最系统方法。其基本思路是，将复杂的真分式X(z)分解为若干个简单分式之和，这些简单分式对应着我们熟知的典型序列的变换式。分解的关键在于求出分母多项式的根（极点），并根据极点的类型（单极点、重极点）确定分解形式。例如，对于单极点情况，X(z)可分解为Σ A_k / (1 - p_k z^-1)的形式，其中每个分式都对应一个指数序列p_k^n u[n]（若收敛域为|z|>|p_k|）。通过求解各系数A_k，即可轻松写出时域序列。

十、 留数定理在z反变换中的应用

       对于更一般的情况，或者当部分分式分解较为繁琐时，基于复变函数理论的留数定理提供了求解z反变换的普适性公式：x[n] = (1/(2πj)) ∮_C X(z) z^n-1 dz，其中积分路径C是收敛域内的一条逆时针闭合曲线。该公式表明，序列值x[n]等于X(z)z^n-1在其所有被C包围的极点上的留数之和。尽管计算涉及复积分，但在实际应用中，对于有理函数形式的X(z)，求留数通常转化为求导或代数运算，比直接积分更为简便。留数定理是验证其他方法结果正确性的理论基石，尤其在处理非有理函数或复杂收敛域时不可或缺。

十一、 系统函数与零极点分析

       在离散系统分析中，我们经常直接面对系统函数H(z)，即系统单位脉冲响应h[n]的z变换。求解H(z)本身往往不是最终目的，而是要通过其表达式分析系统特性。H(z)通常表示为两个z^-1多项式之比。令分子多项式为零得到的根称为零点，令分母多项式为零得到的根称为极点。零极点图（在z平面上标注零极点的位置）是分析系统频率响应、稳定性与因果性的强大可视化工具。例如，系统的稳定性要求所有极点都位于单位圆内。因此，求解一个系统描述（如差分方程）的z变换得到H(z)后，首要步骤之一就是确定其极点位置并判断稳定性。

十二、 由差分方程求解系统函数

       线性常系数差分方程是描述离散时间系统的常见方式。对差分方程两边同时取z变换，并充分利用时移性质，是求解该系统系统函数H(z)的标准流程。考虑一个一般形式的差分方程：Σ_k=0^N a_k y[n-k] = Σ_r=0^M b_r x[n-r]。假设系统初始静止（即零初始条件），对每一项应用时移性质：y[n-k] ←→ z^-k Y(z)， x[n-r] ←→ z^-r X(z)。代入方程并整理，可得Y(z) / X(z) = H(z) = (Σ b_r z^-r) / (Σ a_k z^-k)。通过这一步骤，我们成功地将时域的差分方程转化为了复频域（z域）的代数方程，并直接得到了系统函数的有理分式表达式。

十三、 收敛域的确定原则与系统特性关联

       如前所述，收敛域的确定至关重要。对于有理函数形式的X(z)，其收敛域由极点位置决定，并且必须满足以下原则：收敛域是一个以原点为中心的圆环，且不包含任何极点；对于右边序列，收敛域延伸到最外层极点之外；对于左边序列，收敛域延伸到最内层非零极点之内；对于双边序列，收敛域是一个环状区域。结合系统分析，若H(z)代表一个因果系统的系统函数，则其收敛域必为最外层极点之外的区域。若系统稳定，则该收敛域必须包含单位圆。因此，求解z变换后，根据物理背景或问题约束（如因果性、稳定性）来确定正确的收敛域，是完整解答的必要环节。

十四、 利用z变换求解序列的卷积

       时域卷积定理是z变换最优雅且实用的性质之一。它指出，两个序列在时域的卷积，对应于它们z变换在z域的乘积。即，若x[n] ←→ X(z)， h[n] ←→ H(z)，则x[n] h[n] ←→ X(z) H(z)。这为计算复杂卷积提供了一条捷径：先分别求出两个序列的z变换，然后相乘得到乘积的z变换，最后再通过z反变换（如部分分式分解法）求出卷积结果序列。这种方法尤其适用于序列长度较长或直接进行时域卷积计算量巨大的情况，将复杂的卷积运算转化为相对简单的代数乘法和标准反变换过程。

十五、 帕塞瓦尔定理与能量计算

       z变换还与序列的能量计算密切相关。帕塞瓦尔定理（帕塞瓦尔定理， Parseval's Theorem）的一种形式表明，一个序列的总能量可以通过其z变换在单位圆上的积分来计算。具体地，Σ_n=-∞^∞ |x[n]|^2 = (1/(2πj)) ∮_C X(z) X(1/z) z^-1 dz，其中积分路径C为单位圆。这个定理在信号处理中用于分析信号能量、设计滤波器等。虽然它不直接用于“求解”一个序列的z变换表达式，但它深刻揭示了时域能量与复频域表示之间的内在联系，是z变换理论体系完整性的重要体现。

十六、 数字滤波器设计中的z变换应用

       在数字滤波器设计中，z变换是核心理论基础。无论是无限长脉冲响应滤波器（无限长脉冲响应滤波器， Infinite Impulse Response Filter）还是有限长脉冲响应滤波器（有限长脉冲响应滤波器， Finite Impulse Response Filter），其设计过程都紧密围绕系统函数H(z)展开。例如，在无限长脉冲响应滤波器设计中，通常先根据模拟滤波器原型，通过双线性变换等方法得到数字滤波器的H(z)。这个H(z)的表达式直接决定了滤波器的差分方程系数。因此，“求解z变换”在此语境下，更多地意味着根据性能指标推导或综合出所需的H(z)函数，并确保其零极点分布满足稳定性和频率响应要求。

十七、 结合实例的综合求解演练

       现在，我们通过一个综合实例来串联多个求解技巧。设一个因果系统由差分方程y[n] - 0.7y[n-1] = x[n]描述。首先，对其取z变换（利用线性及时移性质），在零初始条件下得到Y(z) - 0.7z^-1Y(z) = X(z)。整理得系统函数H(z) = Y(z)/X(z) = 1 / (1 - 0.7z^-1)，极点位于z=0.7。由于系统因果，收敛域为|z|>0.7。若要求单位脉冲响应h[n]，即x[n]=δ[n]时的响应，此时X(z)=1，故H(z)=Y(z)。利用已知变换对，直接得到h[n] = (0.7)^n u[n]。若输入为x[n]=u[n]，其变换为X(z)=1/(1-z^-1) (|z|>1)，则输出变换Y(z)=H(z)X(z)=1/[(1-0.7z^-1)(1-z^-1)]。通过部分分式分解，可求得输出序列y[n]。此例展示了从系统描述到响应求解的全过程。

十八、 总结与资源延伸

       求解z变换是一项融合了定义理解、性质运用、数学技巧与工程直觉的技能。从最基础的定义求和，到灵活运用线性、时移、指数加权等性质；从处理右边、左边序列，到掌握幂级数、部分分式、留数定理等核心方法；再到将其应用于系统分析、差分方程求解与滤波器设计，构成了一套层次分明的方法论体系。建议学习者在掌握本文所述方法的基础上，进一步参考权威教材，如奥本海姆的《离散时间信号处理》（Discrete-Time Signal Processing），并通过大量练习巩固。最终，将z变换视为一种思维语言，便能游刃有余地在时域与复频域之间切换，深刻洞察离散时间系统的本质。