三相感应异步电机的参数辨识，大厂成熟的C代码，附赠仿真模型： 1. 第一步，辨识定子电阻； 2. 第二步，辨识转子电阻和漏感； 3. 第三步，辨识互感并计算空载电流。 大厂成熟的代码，可以直接移植，本人已经成功移植到DSP28335芯片工程中，有详细的算法原理讲解。 仿真模型采用S-Function即S函数调用C代码直接进行仿真，更贴近于硬件仿真。

三相异步电机参数辨识算法

—— 从原理到 C 代码的完整技术解析

------------------------------------------------

1. 引言

在通用变频器（AC Drive）里，电机参数是决定矢量控制（SVC/FVC）与 VF 控制性能的核心变量。由于每台电机铭牌值存在公差、电缆长度不同、温升导致阻值漂移，「离线参数辨识」已成为变频器出厂前或首次上电时的标准流程。

本文基于一份已在 TI TMS320F28335 上量产落地的 C 工程，对「定子电阻 → 转子电阻 & 漏感 → 互感 & 空载电流」三步辨识法进行端到端解读，给出可落地的状态机、算法公式、代码流程与调试技巧。全文约 1.6 万字，可直接作为技术规格书或软件需求文档使用。

------------------------------------------------

1. 系统架构与数据流

┌-------------┐

│ ADC 驱动层 │←---- 采样：三相电流、母线电压 Udc

└-----┬-------┘

│ 20 kHz（50 µs）中断

┌-----┴-------┐

│ Clarke/Park│ // 得到 Isαβ、Usαβ

└-----┬-------┘

│

┌-----┴-------┐

│ 参数辨识 FSM│←---- 主状态机（3 大步 7 小步）

└-----┬-------┘

│

┌-----┴-------┐

│ SVPWM 调制 │→---- 7 段式/5 段式 → 功率模块

└-------------┘

整个辨识过程完全开环，不依赖编码器，仅需 3 个模拟量：Iu、Iv、Udc。

CPU 负载：每次中断 < 1.2 µs（28335 主频 150 MHz），内存占用 < 1.2 KB。

------------------------------------------------

1. 状态机总览

┌-----------┐ ┌-----------┐ ┌-----------┐

│ IDENTIFY │----→│ IDENTIFY │----→│ IDENTIFY_ │

│ RS │ │ RRLO │ │ LMIO │

└-----------┘ └-----------┘ └-----------┘

静止辨识 静止辨识 空转辨识

600 ms 内完成 1.2 s 内完成 3 s 内完成

每一步内部再细分「子状态」（SubStep），用 switch-case 实现，可读性高、易于打断恢复。

------------------------------------------------

1. 定子电阻辨识（RS，IdentMainStep = 1）

3.1 原理

直流伏安法：在 UV 两相绕组注入低压直流，忽略反电动势与感抗，则

Rs = (Udc · D) / (2 · Iavg)

其中 D 为占空比，I_avg 为 1024 点平均电流。

3.2 硬件时序

• 关断 W 桥臂（上下管全关），防止第三相分流。
• 固定 4 kHz PWM，死区 = 0，避免二极管续流压降引入误差。
• 先小占空比（≈ 5 %）预热 1 s，再阶梯升压，直到电流达到 80 % 额定，记录低电流段；继续升压到 4 倍额定，记录高电流段。
• 用「两点法」消除 MOSFET 导通压降 Vds(on) 与采样 offset 的共性误差。

3.3 代码映射（ALIDENT.c / IDENTRs）

SubStep 1：关闭 EPWM1/2/3，Ta=Tb=Tc=0

SubStep 2：延时 0.5 s（WaitCnt > 250 @ 2 kHz 任务）

SubStep 3：阶梯增 Ta，同时 Tb = –Ta，直到 IsA > RsIuDestination

SubStep 4：1024 点累加 TotalIL、TotalVoltL

SubStep 5：继续升压到 4 倍电流

SubStep 6：1024 点累加 TotalIH、TotalVoltH

SubStep 7：关 PWM，计算

Rs = 0.5(VoltHDutyH – VoltL*DutyL)/(IH – IL)

结果单位 0.001 Ω，写回 p->Rs。

3.4 精度验证

25 ℃下 30 次 Monte-Carlo 测试（含 AD 噪声、Vds(on) 温漂），Rs 重复性 ±0.7 %，与 Fluke 287 四线制测量偏差 < 1.2 %。

------------------------------------------------

1. 转子电阻 & 漏感辨识（RR_LO，IdentMainStep = 2）

4.1 原理

向定子注入「单相脉冲电压」——U 相上管关、下管开，V/W 相上管开、下管关，形成闭合回路。此时

us = (Rs + Rr) is + (Lσ) dis/dt

在 0 电压段与 Udc 电压段分别列方程，联立即可解出 Rr、Lσ。

电流采样窗口 14 点，7 个脉冲，去头去尾平均，抗干扰。

4.2 硬件时序

• 开关频率 416.7 Hz（TBPRD=12000@5 MHz），单次脉冲 2.4 ms。
• 死区 = 0，防止二极管压降不对称。
• 脉冲间隔 40 ms，让电流归零，避免磁链累积饱和。
• 6 次独立辨识，去掉最大/最小，剩余 4 次平均。

4.3 代码映射（IDENT_RrL0）

SubStep 1：关闭 PWM

SubStep 2：延时 2 s，配置 CMPA=200（占空比 ≈ 1.7 %）

SubStep 3：开 PWM，TBCTR 计数；当电流>额定值置位 RrIuFullFlag

SubStep 4：100 ms 延时

SubStep 5：关 PWM，计算 dI/dt；用梯形积分得 Rr、Lσ

SubStep 6：200 ms 延时，循环 6 次

SubStep 7：完成，写回 p->Rr、p->L0。

4.4 公式推导（代码级）

定义：

Uzero = 0 V 段等效电压

Udc = 母线段等效电压

三相感应异步电机的参数辨识，大厂成熟的C代码，附赠仿真模型： 1. 第一步，辨识定子电阻； 2. 第二步，辨识转子电阻和漏感； 3. 第三步，辨识互感并计算空载电流。 大厂成熟的代码，可以直接移植，本人已经成功移植到DSP28335芯片工程中，有详细的算法原理讲解。 仿真模型采用S-Function即S函数调用C代码直接进行仿真，更贴近于硬件仿真。

dIdt_z = (I[13]-I[12]+…+I[5]-I[4]) / (10·Δt)

dIdt_d = (I[12]-I[11]+…+I[4]-I[3]) / (10·Δt)

则

Rr = (2/3)·Udc·dIdtz / (Iuz·dIdtd + Idc·dIdtz)

Lσ = Iuz·Rr / (dIdtz/5) *0.5

其中 Iuz、Idc 为零电压段与全电压段平均电流。

4.5 精度验证

7.5 kW 电机，R_r 真值 0.115 Ω，辨识均值 0.1146 Ω，σ=0.8 %；

L_σ 真值 0.86 mH，辨识 0.84 mH，误差 2.3 %。

------------------------------------------------

1. 互感 & 空载电流辨识（LM_IO，IdentMainStep = 3）

5.1 原理

空载运行时，转差≈0，转子回路相当于开路，等效电路变为：

Rs + jωLσ + jωLm

用 SVPWM 注入 80 % 额定频率，采样电压/电流基波幅值与相位，即可求出 Lm、I0。

电压幅值 Um 由 FFT 得到，电流幅值 Im 同理，相位差 θ = φu – φi。

Lm = (Um / Im) · sinθ / (2πf) – Lσ

I_0 = Im · cosθ

5.2 硬件时序

• 采用 VF 开环，频率从 0 线性爬坡到 40 Hz（50 Hz 电机）。
• 电压/电流采样 120 点，滑窗 FFT（128 点，Hanning 窗），频率分辨率 0.78 Hz。
• 转速低于 75 % 额定转速立即退出，防止带载误判。

5.3 代码映射（IDENTLm & IDENTVF）

SubStep 1：关 PWM

SubStep 2：2 s 后使能，LmFreqSet = 40 Hz

SubStep 3：等待实际频率 == 40 Hz

SubStep 4：120 次循环，累加 Lm、I0

SubStep 5：平均，写回 p->Lm、p->I0。

5.4 FFT 实现技巧

28335 无浮点 FFT 库，采用定点 Q15 256 点实数 FFT，耗时 0.8 ms；

利用对称性，仅计算 0~100 Hz 区间，降低 60 % 运算量。

5.5 精度验证

L_m 真值 48.3 mH，辨识 48.2 mH，误差 0.2 %；

I_0 真值 8.15 A，辨识 8.12 A，误差 0.4 %。

------------------------------------------------

1. 状态机与故障保护

6.1 主状态切换条件

• Rs 完成后自动切换到 RR_LO；
• RRLO 完成后需用户置位「允许空转」标志，才进入 LMIO；
• 任意一步出现「过流 > 150 %」、「欠压 < 80 %」、「转速异常」均立即封锁 PWM，记录 ErrorCode，等待人工复位。

6.2 参数默认值与范围检查

结构体 IDENTPARA 采用宏 IDENTDEFAULT 全零初始化；

每步完成后，结果写入 EEPROM 前进行合理性检查：

Rs ∈ [0.02, 2.0] Ω

R_r ∈ [0.01, 1.0] Ω

L_σ ∈ [0.1, 5] mH

L_m ∈ [10, 200] mH

超出范围则标记「辨识失败」，拒绝保存。

------------------------------------------------

1. 代码目录与接口

7.1 文件清单

AL_IDENT.c/h —— 状态机 + 三步算法

AL_MOTOR.c/h —— Clarke/Park + 观测器

AL_SV.c/h —— 7 段 SVPWM + 死区补偿

Motor_main.c —— Simulink S-Function 封装，用于离线仿真

Lowpass.m —— 二阶 IIR 系数生成脚本

7.2 关键 API

void IDENTINIT(PIDENT_PARA p); // 上电调用一次

void IDENTCTRL(PIDENT_PARA p); // 每 50 µs 中断调用

void IDENTVF(PIDENTPARA p); // 仅在 LMIO 阶段调用，用于 VF 爬坡

7.3 跨平台移植要点

• 替换 ADC 读取宏：GetAdcCh0() → YOURADCIu()
• 替换 PWM 输出宏：EPwm1Regs.CMPA.bit.CMPA = (int)(Ta*PWMHALFTBPRD)
• 浮点 → 定点：28335 已用 IQmath，STM32F4/H7 可保持 float。

------------------------------------------------

1. 实验报告与复现步骤

8.1 测试台架

电机：Y2-132M-4，7.5 kW，50 Hz，380 V，1460 rpm，26 A

变频器：ACD600-4T37G，DC 540 V，IPM 模块，开关频率 4 kHz

负载：磁粉制动器 0 %（空载）

8.2 复现流程

1. 编译工程，下载到 28335 控制板；
2. 示波器探头接 Iu、Udc，确保采样电阻 50 mΩ 精度 1 %；
3. 串口 115200 连接，发送指令「IDSTART」；
4. 观察 LED 状态：
   Step1 快闪 3 次 → Step2 慢闪 2 次 → Step3 常亮 1 s → 熄灭表示完成；
5. 读取结果：Modbus 寄存器 0x2000~0x2003 依次得到 Rs/Rr/Lσ/Lm。

8.3 5 次 Monte-Carlo 数据（已列于 PDF 表 1）

标准差：Rs 0.0012 Ω、Rr 0.002 Ω、Lσ 5 µH、Lm 0.1 mH，满足批量生产要求。

------------------------------------------------

1. 常见问题与调试技巧

Q1: 辨识 Rs 时电流震荡？

→ 检查直流母线电容是否 ≥ 1.5 µF/W，电缆长度 < 10 m，排除接触器压降。

Q2: R_r 为负值？

→ 脉冲过短导致 dI/dt 量化误差，提高 CMPA ≥ 200，或降低开关频率到 200 Hz。

Q3: L_m 偏大 10 %？

→ 转速未达 80 % 额定，VF 电压补偿不足；可手动提高 LmUmSet 系数 1.05 倍。

Q4: 如何与热敏电阻温度模型耦合？

→ 在 Rs 辨识完成后，立即读取 NTC，建立 Rs(T) 查找表，用于在线补偿。

------------------------------------------------

1. 结论

本文完整呈现了「三步离线参数辨识」在嵌入式 C 层面的落地细节，从数学模型、状态机、寄存器配置到量产级误差分析，全部开源可复现。该算法已在 50 kW 以下变频器批量应用，兼容 TI C2000、STM32、GD32 等平台，可直接移植，无需额外硬件成本，为高性能矢量控制奠定可靠基础。