数字控制系统中的Z变换与PID实现详解

1. 数字控制系统中的Z变换基础

1.1 从差分方程到Z变换

在数字控制系统中，我们处理的不是连续的模拟信号，而是离散时间序列。传统控制理论中的微分方程在这里演变为差分方程。以一个典型的二阶系统为例，其差分方程可以表示为：

xₙ = a₁xₙ₋₁ + a₂xₙ₋₂ + b₁uₙ₋₁ + b₂uₙ₋₂

这种方程直接求解相当繁琐，特别是当系统加入反馈后，计算复杂度会呈指数级增长。Z变换的出现完美解决了这个问题——它将差分方程转换为关于变量z的代数方程，使系统分析变得直观简单。

Z变换的定义式为：
X(z) = Σ xₙz⁻ⁿ (n从0到∞)

这个变换有几个关键特性值得注意：

• 线性性质：满足叠加原理
• 时移特性：xₙ₋ₖ ↔ z⁻ᵏX(z)
• 卷积定理：时域卷积对应z域乘积

1.2 传递函数与系统特性

通过Z变换，我们可以得到系统的传递函数H(z)=Y(z)/U(z)。以电机系统为例，其传递函数可能呈现如下形式：

H(z) = b₀ + b₁z⁻¹ / (1 + a₁z⁻¹ + a₂z⁻²)

传递函数的分母多项式被称为特征多项式，它的根决定了系统的固有特性：

• 极点位置决定系统稳定性（单位圆内稳定）
• 极点距离单位圆的远近影响响应速度
• 复数极点会产生振荡响应

在实际工程中，我们常用以下经验法则：

1. 主导极点应位于0.3-0.8半径范围内
2. 阻尼系数最好控制在0.5-0.7之间
3. 采样频率至少是系统带宽的10倍

2. PID控制器的数字实现

2.1 离散化PID算法

传统PID控制器的连续形式为：
u(t) = Kₚe(t) + Kᵢ∫e(t)dt + Kₚde(t)/dt

在数字系统中，我们需要对其进行离散化处理：

• 积分项采用累加近似：∫e(t)dt ≈ TₛΣeₖ
• 微分项采用差分近似：de(t)/dt ≈ (eₖ - eₖ₋₁)/Tₛ

这样得到的离散PID算法为：
uₙ = Kₚeₙ + KᵢTₛΣeₖ + Kₚ(eₙ - eₙ₋₁)/Tₛ

2.2 实用PID代码实现

以下是经过工程验证的PID控制器C++实现：

cpp复制class DigitalPID {
public:
    DigitalPID(float Kp, float Ki, fl