当前位置：首页 > news >正文

【MATLAB】基于 MATLAB/Simulink 的无刷直流电机（BLDC）转速控制模糊 PID 算法

news 2026/5/25 10:12:22

基于 MATLAB/Simulink 的无刷直流电机（BLDC）转速控制模糊 PID 算法

摘要：针对永磁无刷直流电机（Brushless DC Motor, BLDC）在常规 PI 转速控制下存在的启动超调大、稳态转矩波动剧烈的问题，本文提出一种由三个并行模糊逻辑控制器（Fuzzy Logic Controller, FLC）在线整定 PID 三个增益的调速控制策略。以 48 V/1 kW/700 rpm 等级 BLDC 为控制对象，在 MATLAB/Simulink 环境下以信号流形式搭建电机数学模型，结合六步换相逻辑构成被控对象；进而以转速误差 e 与误差变化率 de 为模糊输入，构建 5×5 共 25 条 Mamdani 规则的模糊推理系统，输出 ΔKp、ΔKi、ΔKd 实时叠加到 PID 基线增益上。仿真结果表明：在 650 rpm 转速指令下，所提方法的转速超调由固定增益 PI 的 28.2% 降至 3.7%，±2% 进入时间由 97 ms 缩短至 19.4 ms；稳态转速误差小于 0.02 rpm，控制效果显著优于固定增益 PI。

关键词：无刷直流电机；模糊 PID；转速控制；Simulink 仿真；anti-windup

1 引言

无刷直流电机以效率高、转矩密度大、可靠性好等优点，在电动汽车、伺服系统、工业自动化等领域应用广泛。其转速控制最常用的方法是经典的 PI 或 PID 控制，结构简单、易于工程实现。然而 BLDC 本身具有强非线性（换相纹波、相电流方波等）、参数不确定（温度引起的 Rs 漂移、负载变化）以及电气-机械时间常数差异大等特点，固定增益 PID 在大范围转速指令下难以兼顾快速性与平稳性：增益高则启动超调与稳态振荡严重，增益低则响应迟缓。

模糊控制具有不依赖精确模型、能够嵌入工程经验、对参数摄动鲁棒的优点，将其与 PID 结合形成"模糊整定 PID"是一类有效折中。本文工作的主要贡献是：

在 Simulink 中以信号流形式自建一套 5 维状态的 BLDC 数学模型，配合 Stateflow 实现的六步换相逻辑，避免对 Simscape 物理建模工具箱版本差异的依赖，模型可由脚本程序化生成与重建；
构建三个独立的 Mamdani 模糊推理系统，分别对 Kp、Ki、Kd 进行在线整定，并对模糊规则的方向、论域、PID 抗积分饱和模式进行系统性调整以适配本电机参数；
在统一的电机对象与统一的工况下，定量比较固定增益 PI 与模糊整定 PID 的瞬态与稳态指标，给出完整的对比证据链。

2 BLDC 数学模型

2.1 状态空间方程

考虑三相星形连接、梯形反电动势的 BLDC，忽略涡流损耗与齿槽转矩，相变量下的电气方程为：

{diadt=1Ls(va−Rsia−ea)dibdt=1Ls(vb−Rsib−eb)dicdt=1Ls(vc−Rsic−ec) \begin{cases} \displaystyle\frac{di_a}{dt}=\frac{1}{L_s}(v_a-R_s i_a - e_a) \\ \displaystyle\frac{di_b}{dt}=\frac{1}{L_s}(v_b-R_s i_b - e_b) \\ \displaystyle\frac{di_c}{dt}=\frac{1}{L_s}(v_c-R_s i_c - e_c) \end{cases}⎩⎨⎧dtdia=Ls1(va−Rsia−ea)dtdib=Ls1(vb−Rsib−eb)dtdic=Ls1(vc−Rsic−ec)

机械方程为：

Jdωdt=Te−Bω−TL,dθmdt=ω J\frac{d\omega}{dt}=T_e-B\omega-T_L,\quad \frac{d\theta_m}{dt}=\omegaJdtdω=Te−Bω−TL,dtdθm=ω

电磁转矩由三相瞬时电功率除以机械角速度得到：

Te=eaia+ebib+ecicω T_e=\frac{e_a i_a + e_b i_b + e_c i_c}{\omega}Te=ωeaia+ebib+ecic

三相反电动势 e_a、e_b、e_c 取决于电角度 θ_e = Pp·θ_m，呈 120° 相差的梯形波形，幅值与 ω 成正比、比例系数为反电动势常数 ke。状态向量取 x = [i_a, i_b, i_c, ω, θ_m]^T，由 Stateflow EMChart 实现状态导数与副输出（i_abc、e_abc、T_e、ω）的同步计算，外部通过 Integrator 块完成 5 维状态积分。

2.2 六步换相

由机械角度 θ_m 经极对数 Pp 映射至电角度 θ_e ∈ [0, 2π)，将一个电周期划分为六个 60° 扇区。每个扇区中两相导通、一相悬空，导通极性由扇区号决定，第三相电压取 0。占空比 d ∈ [0, 1] 与母线电压 Vdc 直接相乘得到导通相的端电压幅值，构成最简化的等效驱动方式；该简化省略了 PWM 高频开关纹波，但保留了真实换相纹波，足以体现 BLDC 在控制层面的关键特征。

2.3 电机参数

电机额定参数与电气-机械参数列于表 1。Rs、Ls、ke、Pp 取同等级 48 V/1 kW 商用 BLDC 的经验值；转动惯量 J 与摩擦系数 B 按"机电时间常数 ≈ 2 s"匹配启动响应的时间尺度，并保持 J/B 比例使稳态摩擦转矩量级合理。

项	符号	取值	单位
母线电压	Vdc	48	V
额定功率	P_N	1	kW
最高转速	n_max	700	rpm
参考转速	n_ref	650	rpm
定子电阻	Rs	0.36	Ω
定子电感	Ls	0.6	mH
反电动势常数	ke	0.05	V·s/rad
极对数	Pp	4	—
转动惯量	J	1×10⁻²	kg·m²
摩擦系数	B	5×10⁻³	N·m·s

表 1 BLDC 电机仿真参数

3 控制策略

3.1 固定增益 PI 基线

基线控制器采用单转速闭环：

u(t)=Kpe(t)+Ki ⁣∫0t ⁣e(τ)dτ,e(t)=nref−n(t) u(t)=K_p e(t)+K_i\!\int_0^t\! e(\tau)d\tau,\quad e(t)=n_{\text{ref}}-n(t)u(t)=Kpe(t)+Ki∫0te(τ)dτ,e(t)=nref−n(t)

输出 u(t) 限幅到 [0, 1] 作为占空比指令送入换相模块。增益取 (Kp, Ki) = (0.15, 5)，由对所提电机模型在 7 组备选值的网格试凑后选定，作为后续模糊整定 PID 的对照基线。

3.2 模糊整定 PID 总体结构

所提模糊整定 PID 控制系统的数据流如图 1 所示。

3.3 模糊推理系统设计

3.3.1 输入隶属度

转速误差 e 论域取 [−700, 700] rpm，覆盖最高转速 ±n_max；误差变化率 de 论域取 [−60000, 60000] rpm/s，覆盖启动期 |de|_max ≈ 50000 rpm/s。两输入均采用 5 档三角形隶属度，标签为 {NB, NS, Z, PS, PB}，对称分布，相邻函数 50% 重叠，如图 2 所示。