电流环阶跃响应参数详细分析：PI 参数、带宽与整定

1. 为什么看阶跃响应：电流环是系统的天花板

FOC 三闭环像俄罗斯套娃：电流环在里，速度环包电流环，位置环包速度环。任何外环的命令最终都要通过电流环落地——所以：

电流环带宽决定速度环上限：工程惯例速度环带宽取电流环的 1/5~1/10（见 06 篇三闭环整定）；
电流环阶跃响应决定"出力速度"：加速多快、抗扰动多快，都由它决定；
电流环参数最不能拍脑袋：Kp/Ki 与电机 R/L 直接挂钩，公式算出来就八九不离十——它是三环里唯一"可以解析设计"的环。

阶跃响应就是电流环的"体检报告"：给一个方波电流命令，看 iq 怎么追——多快（tr）、过不过冲（Mp）、稳不稳（ts），全在一条曲线上。

2. 模型：R/L 一阶对象 + PI 零极点对消

忽略反电动势与耦合项（低速/静止整定时的标准简化），d 轴或 q 轴电流通道就是一个一阶 RL 对象：

对象（一阶惯性）：     G(s) = 1 / (L·s + R)

        i(s)     Kp·s + Ki     1
闭环:   ──── = ────────── · ────────
        i*(s)      s        L·s + R

PI 的零点 (s = -Ki/Kp) 若能对消对象的极点 (s = -R/L)：
        Ki/Kp = R/L   ⟹   Ki = Kp · R/L   （零极点对消条件）

对消后闭环退化成干净的一阶低通：
        i(s)     Kp/L            Kp
        ──── = ────────── ,   ωc = ───  (rad/s)
        i*(s)   s + Kp/L          L

        —— 这就是"电流环 PI 可以解析设计"的全部秘密：
           ① 按对消条件定 Ki = Kp·R/L
           ② 按目标带宽定 Kp = ωc·L
            

两个式子就是电流环整定的标准答案——这也是电机参数计算器（tools/motor-param.html）里 PI 公式的来源。前提：R/L 估得准、反电动势被前馈或低速忽略。

3. 参数与闭环带宽：Kp 定带宽，Ki 定对消

参数	决定什么	解析关系	调整后果
Kp	闭环带宽（响应快慢）	f_bw ≈ Kp / (2π·L)	↑Kp 响应快；过大 → 受采样/延迟限制振荡
Ki	低频增益 / 稳态误差消除	对消点 Ki = Kp·R/L	↑Ki 偏离对消 → 引入零点、出现超调；过小 → 稳态误差
L、R	对象的极点位置	极点 = R/L	参数估错 → 对消失效，响应变形
采样/载频 fs	可实现带宽上限	f_bw_max ≈ fs / 10	载频越高，能压的带宽越宽（受延迟制约）

工程估算：L=1mH 的电机想要 500Hz 电流环带宽，Kp = 2π·500·0.001 ≈ 3.1（V/A）；Ki = Kp·R/L（R=0.5Ω → Ki ≈ 1570）。算完先上板，再按第 9 节实测微调。

4. 阶跃响应指标：tr / Mp / ts / 稳态误差

看阶跃曲线先认四个指标（示意：对消理想的一阶系统）：

指标	定义	一阶系统理论值	实际意义
上升时间 tr	10% → 90% 目标值	tr ≈ 2.2 / ωc	出力速度
超调量 Mp	峰值超过稳态值的百分比	0%（一阶无超调）	超调 = Ki 偏大或延迟
调节时间 ts	进入 ±2% 误差带	ts ≈ 4 / ωc	稳定速度
稳态误差	\|i(∞) − i*\|	0（PI 积分消除）	Ki=0 时才有

注意口径：tr 和 ts 都与带宽成反比（ωc 越大越快），超调则完全来自"非对消项"——看曲线时先问"有没有超调、tr 合不合理"，再决定动 Kp 还是 Ki。

5. Kp 扫描：从迟钝到振荡

固定 Ki 按对消条件（Ki = Kp·R/L），只拉 Kp，看响应形态（示意，fs=20kHz、L=1mH、R=0.5Ω）：

Kp	带宽 f_bw	响应形态	问题
0.5	≈ 80 Hz	慢、无超调、tr ≈ 4.4ms	太迟钝，速度环被拖累
2.0	≈ 320 Hz	理想一阶，tr ≈ 1.1ms	常规工作点
4.0	≈ 640 Hz	仍一阶但逼近采样/延迟极限	边缘，轻微过冲
8.0	≈ 1.3 kHz（>fs/10）	振荡/超调 20%+	撞上奈奎斯特/延迟，不能这么用

规律：Kp 是响应速度的旋钮，但存在"载频墙"——把带宽压到 fs/10 以上，离散化延迟开始主导，曲线从"一阶"变"二阶振荡"。这就是为什么载频 20kHz 的系统很少把电流环带宽调过 2kHz。

6. Ki 偏离对消：超调从哪来

对消条件是 Ki = Kp·R/L。把 Ki 调大（Ki > Kp·R/L），PI 零点跑向高频，闭环出现额外零点 + 极点对，超调随之而来：

示例：Kp=2（f_bw≈320Hz），R=0.5Ω，L=1mH → 对消 Ki = 1000

    Ki = 500 （偏小）  → 零点太近，响应慢 + 小幅欠冲，稳态仍无差
    Ki = 1000（对消）  → 干净一阶，0% 超调，tr ≈ 1.1ms
    Ki = 2000（偏大）  → 零点引入，超调 ≈ 8%，tr 略快但振荡尾巴
    Ki = 4000（过大）  → 超调 ≈ 20%，ts 明显变长（积分的"回摆"）

结论：Ki 的唯一"正确值"在 R/L 附近。
      超调出现时先怀疑 Ki 偏大，而不是 Kp——很多人调反了。
            

为什么偏大的 Ki 会加超调：PI 零点 (Ki/Kp) 与对象极点不对消后，零点在闭环里表现为"微分前馈"效应——命令一到就推过头，然后积分再慢慢拉回来，形成过冲。这也是为什么"先按公式对消、再微调 Kp"是最稳的整定路线。

7. 数字实现的两大杀手：采样频率与一拍延迟

理论是连续域，板上是离散域。两个效应会把"算好的阶跃响应"变差：

采样频率 fs（=PWM 载频）：电流采样与 PWM 同步，一个 PWM 周期算一次 PI。带宽不能超过 fs/10 附近（工程经验），否则离散化误差让相位滞后累积、响应振荡；
一拍延迟（计算 + 占空比更新）：本周期采的电流，算完要下一个周期才生效——等效在环路里加 e^(−sTs) 延迟。延迟让"临界 Kp"变小：同样带宽下，延迟越大超调越多。

仿真里最直观（第 10/11 节）：同一组 Kp/Ki，关掉延迟是干净一阶；打开一拍延迟，超调马上冒出来——板上看到的"多出来的超调"，一大半是延迟给的，不是参数给的。对策：采样时刻对准 PWM 中点、占空比下一周期更新前先算完（双缓冲）、Kp 按"延迟折减"取保守值。

8. 电压限幅与积分饱和：anti-windup 不能省

阶跃命令很大时（比如 10A 阶跃），PI 输出会撞上母线电压限幅 ±V_bus：

现象：积分项继续累积 → 命令撤销后积分"撑着"输出 → 电流过冲/振荡
      这就是积分饱和（windup）。

对策（anti-windup，三选一即可）：
 ① 钳位积分：v_sat 时冻结积分累加（最简单，最常用）
 ② 反算补偿：v_sat 时把多出的量从积分里扣回去（背计算法）
 ③ 条件积分：|e| 大时清零积分（适合大阶跃场景）

C 里的钳位实现（第 10 节代码同款）：
    v_raw = Kp*e + acc;
    if (v_raw >  V_MAX) { v_out =  V_MAX; }   /* 限幅 */
    else if (v_raw < -V_MAX) { v_out = -V_MAX; }
    else { v_out = v_raw; acc += Ki*e*Ts; }   /* 仅未饱和时积分 */
    （注意：e 符号向 — 限幅时积分为零即"冻结"）*/
            

整定大阶跃响应时，先确认曲线是不是撞限幅了——撞限幅的响应跟 PI 参数没关系，是"出力上限"决定的（母线电压 / 电机反电动势的物理墙）。第 11 节仿真器里把阶跃幅值调大就能看到这个分界线。

9. 实测方法：示波器看电流阶跃

板上验证三步（配合 11 篇的工具链思路）：

① 给阶跃：上位机/CANopen 发电流环阶跃命令（比如 iq: 0 → 1A → 0）
            —— 用示波器电流探头或 DAC 输出 iq 观测点

② 看曲线：记录 iq 波形，量 tr / Mp / ts（第 4 节口径）

③ 判读调整：
    - 无超调、tr 偏慢          → Kp 调大（按 ωc=Kp/L 反推）
    - 有超调且振荡            → 先查是否撞限幅；再查延迟；最后 Ki 偏大
    - 稳态有差                → Ki 太小（或对消点 R/L 估错）
    - 高转速下响应变形        → 反电动势/耦合项进入，需前馈（见 06 篇）

推荐顺序：先空载、低速、小阶跃整定；再带载、高速验证——每步都记曲线。
            

实测与仿真对照是最高效的：先用第 11 节仿真器把参数调到"理论曲线好看"，再上板复现——两者对不上，先怀疑采样/延迟/限幅，而不是参数。

10. 可编译 C 代码：离散 PI + 一阶对象阶跃仿真

完整可编译：离散电流环（PI + anti-windup + 一拍延迟开关）+ R/L 对象递推，跑阶跃并打印指标：

/* iq_step.c —— 电流环阶跃响应离散仿真（零极点对消 + anti-windup）*/
#include <stdio.h>
#include <math.h>

/* 电机与环路参数（示意）*/
static const double L  = 1e-3;     /* 相电感 H */
static const double R  = 0.5;      /* 相电阻 Ω */
static const double FS = 20000.0;  /* 采样/载频 Hz */
static const double TS = 1.0/FS;
static const double KP = 2.0;     /* 按 f_bw≈320Hz */
static const double KI = KP*R/L;   /* 零极点对消 Ki = Kp·R/L */
static const double VMAX = 24.0;   /* 母线限幅 V */
#define DELAY_ONE  1    /* 1=开一拍延迟, 0=关 */

int main(void) {
    double i = 0, acc = 0, v_out = 0, v_dly = 0;
    double iref = 1.0;        /* 1A 阶跃 */
    double i_peak = 0, t10 = -1, t90 = -1, t_settle = -1;
    int n, settled = 0;
    const int N = 4000;      /* 200ms */

    for (n = 0; n < N; n++) {
        double e = iref - i;
        double v_raw = KP*e + acc;
        /* 限幅 + anti-windup：未饱和才积分 */
        if (v_raw >  VMAX)      { v_out =  VMAX; }
        else if (v_raw < -VMAX) { v_out = -VMAX; }
        else                    { v_out = v_raw; acc += KI*e*TS; }
                /* 一拍延迟：关=直接用本拍输出；开=用上一拍输出(v_dly)，本拍排队下一周期生效 */
        double v_app;
        if (DELAY_ONE) { v_app = v_dly; v_dly = v_out; }
        else           { v_app = v_out; }
        /* 一阶对象递推: i[k+1] = i[k] + Ts/L·(v - R·i) */
        i += TS/L*(v_app - R*i);
        /* 指标采集 */
        if (i > i_peak) i_peak = i;
        if (t10 < 0 && i >= 0.1*iref) t10 = n*TS;
        if (t90 < 0 && i >= 0.9*iref) t90 = n*TS;
        if (!settled && n > 200 && fabs(i-iref) < 0.02*iref) { settled = 1; t_settle = n*TS; }
    }
    printf("Kp=%.1f Ki=%.1f fs=%.0fHz 延迟=%s\n", KP, KI, FS, DELAY_ONE ? "开" : "关");
    printf("tr = %.3f ms (%.3f→%.3f)\n", (t90-t10)*1e3, t10*1e3, t90*1e3);
    printf("ts = %.3f ms\n", t_settle*1e3);
    printf("超调 Mp = %.1f%%  稳态 i = %.4f A\n", (i_peak-iref)/iref*100, i);
    printf("理论带宽 f_bw = %.0f Hz (Kp/2πL)\n", KP/(2*M_PI*L));
    return 0;
}
/* 编译: gcc iq_step.c -lm -o iq_step && ./iq_step
   把 DELAY_ONE 改 0 再看一次——超调就是延迟给的 */
            

这个仿真就是第 11 节演示器的命令行版——改参数、重编译、看指标，适合批量扫参数；交互版适合直观拖滑块。

11. 交互演示：电流环阶跃响应仿真器

拖 Kp / Ki / L / R / 载频五个滑块，看 iq 阶跃曲线与四项指标实时变化——先试"零极点对消"（Ki 自动 = Kp·R/L 附近），再故意拉大 Ki 看超调，最后打开一拍延迟对比（示意仿真）：

电流环阶跃响应仿真器（iq: 0 → 1A）

Kp (0.5~8)

Ki (100~6000)

L(mH) (0.3~3)

R(Ω) (0.1~2)

载频 fs(kHz) (5~100)

母线 Vbus(V) (5~48)

一拍延迟（仿真 PWM 更新）

默认即对消点（Kp=2、Ki=1000、L=1mH、R=0.5Ω、fs=20kHz、延迟关）：干净一阶，0% 超调、tr≈1.1ms——对照第 3 节公式 f_bw≈320Hz。
把 Ki 拉到 4000：超调约 20% 冒出来（第 6 节）；勾上"一拍延迟"再拉高 Kp：超调与振荡尾巴更早出现（第 7 节）。
把 Kp 拉到 8：tr 变快但 ts 反而拖长（离散零阶保持的尾巴）；再勾上"一拍延迟"：超调约 8% 冒出来——这就是第 7 节的"载频墙"。