⚙️ 位移定理的双面性：时移与频移

Z变换的位移定理主要分为两大类：时域位移定理和频域位移定理。理解这两者的区别与联系是掌握Z变换应用的关键。

1. 时域位移定理

时域位移定理描述了序列x[n]在时间轴上平移k个单位后，其Z变换的变化规律。这是求解差分方程最核心的工具。

? 双边Z变换的时移性质

若 ，收敛域为R，则对于整数k：

Z{x[n-k]} = z^{-k}X(z)

这说明时域的延迟（k>0）对应于Z域的乘法因子z^(-k)，即引入一个极点在原点（如果k>0）或零点在原点（如果k<0）。

? 单边Z变换的时移性质（重点）

在工程实际中，我们常处理因果信号或具有初始条件的系统，此时需使用单边Z变换。单边Z变换的时移定理更为复杂，因为它必须考虑平移后“丢失”或“产生”的初始值。

右移（延迟）k步 (k>0)：

Z⁺{x[n-k]} = z^{-k}X⁺(z) - Σ_{m=0}^{k-1} x[m]z^{-(k-m)}

左移（超前）k步 (k>0)：

Z⁺{x[n+k]} = z^{k}X⁺(z) - Σ_{m=0}^{k-1} x[m]z^{k-m}

注：单边Z变换的右移定理中，减去的项代表了被移出时间窗口的初始条件信息，这是求解非零初始条件差分方程的关键。

2. 频域位移定理（调制性质）

频域位移定理，也称为Z域的尺度变换性质或调制性质。它描述了时域序列乘以指数序列后，其Z变换在Z域中的变化。

? 定理表述

若 ，收敛域为R，则对于复常数a：

Z{a^n x[n]} = X(z/a)

这意味着时域的指数加权对应于Z域变量的尺度变换。收敛域变为 |a|R。

? 典型应用：复指数调制

当 a = e^(jω₀) 时，该定理变为：

Z{e^{jω₀n} x[n]} = X(z e^{-jω₀})

这一性质在通信系统中极为重要，用于将基带信号的频谱搬移到载波频率处。在数字滤波器设计中，它也用于频率响应的平移。

? 工程应用：从理论到实践

Z变换的位移定理不仅是数学推导的工具，更是解决实际工程问题的利器。以下是两个主要应用场景：

① 求解线性常系数差分方程

离散系统的动态行为通常由差分方程描述。例如，一个简单的一阶系统：

y[n] - 0.5y[n-1] = x[n]

利用时移定理，我们可以将方程两边取Z变换，将差分方程转化为代数方程，从而快速求解系统响应Y(z)，再通过逆Z变换得到时域解y[n]。这对于分析电路暂态过程、经济模型预测等非常有效。

② 数字滤波器设计

在设计IIR（无限脉冲响应）滤波器时，位移定理帮助我们将滤波器的传递函数H(z)与时域脉冲响应h[n]联系起来。通过部分分式展开和逆Z变换，工程师可以确定滤波器的冲激响应，进而评估其稳定性和相位特性。

? 发展历程与关键节点

1947年

W. R. Evans 首次引入Z变换概念，用于分析离散数据控制系统。

1950s

H. B. Curry 和 J. R. Ragazzini 等人完善了Z变换的数学理论，确立了位移定理在求解差分方程中的地位。

1965年

Cooley-Tukey FFT算法的提出，虽然主要基于DFT，但极大地推动了基于Z变换的数字信号处理算法的实时应用。

现代

Z变换理论被广泛应用于5G通信、音频处理、医学图像重建等领域，位移定理依然是DSP教材的核心章节。

? 典型例题与代码实现

为了帮助读者更好地理解，我们提供一个经典的求解差分方程的例子，并附上MATLAB代码验证。

例题：求解差分方程

已知系统差分方程为 y[n] - 0.8y[n-1] = x[n]，初始条件 y[-1] = 1，输入 x[n] = u[n]（单位阶跃）。求y[n]。

解：

  1. 两边取单边Z变换：
    Y(z) - 0.8[z⁻¹Y(z) + y[-1]] = X(z)
    （注意：这里使用了单边Z变换的右移定理，并代入了初始条件y[-1]=1）
  2. 代入数值：
    Y(z) - 0.8z⁻¹Y(z) - 0.8(1) = X(z)
    Y(z)(1 - 0.8z⁻¹) = X(z) + 0.8
  3. 求解Y(z)：
    Y(z) = X(z)/(1 - 0.8z⁻¹) + 0.8/(1 - 0.8z⁻¹)
  4. 逆Z变换：
    y[n] = (0.8)ⁿ u[n] + (0.8)(0.8)ⁿ u[n] ... (具体计算略，需结合X(z)=z/(z-1))

MATLAB 验证代码

% 定义系统系数
b = 1;      % 分子系数
a = [1, -0.8]; % 分母系数
% 定义输入信号 u[n]
n = 0:20;
x = ones(1, length(n));
% 使用 filter 函数求解差分方程 (内部隐含了位移定理的应用)
y = filter(b, a, x);
% 考虑初始条件 y[-1]=1，需要手动修正或设置 initial conditions
zi = filter(designfilt('lowpassiir', ...)); % 此处仅为示意，实际应使用 filtic
% 更准确的包含初始条件的求解：
zi = filtic(a, b, 1); % 根据y[-1]=1构造初始状态
y_ic = filter(b, a, x, zi);
plot(n, y_ic, 'o-');
title('System Response with Initial Condition');
xlabel('n');
ylabel('y[n]');