简介：电机设计MATLAB源代码合集，由作者基于电机学典型章节整理而成，共17个m文件，压缩包整体仅11KB，主要面向电机方向学生与初步接触电机设计的工程师。代码覆盖直流/交流电机模型、磁路与电路分析、转矩功率计算、控制策略设计与参数辨识等内容，各文件对应具体习题或实验场景，每个脚本都保留清晰变量与赋值结构，便于结合教材逐项对照学习并直接修改复用。文件小巧，可直接在MATLAB中运行，快速理解从电机建模到仿真分析的完整思路，还能替换参数、扩展算法做更深入的验证。已有265人浏览学习，适合想在电机设计实践中提升MATLAB应用能力、弥补理论到仿真鸿沟的初学者，是一份轻量实用的入门参考。
1. 电机设计MATLAB源代码：为什么方案阶段先跑脚本
很多做电机的人一提到电磁设计，第一反应是开有限元仿真。但实际产品开发里，第一版电磁方案绝大多数不是仿真出来的，而是一段几百行的 MATLAB 脚本先算出来的。原因很直接：一次瞬态场仿真少说几十分钟，脚本几十毫秒出一版。
额定功率、电压、转速这几个数填进去，主尺寸、磁路、绕组、效率一整条线就出来了。方案可行了，再进有限元细校，能省掉大量无效的仿真排队。这篇文章把一套实际的电机设计 MATLAB 源代码拆开讲，从计算模型到每个模块的参数，再到最容易翻车的几个地方。新手能照着跑通一条完整设计流程，熟手可以拿去做方案排比和参数扫描。
2. 先看脚本背后的物理模型：主尺寸、磁路、绕组一条线怎么走
2.1 从额定指标到主要尺寸：先定转矩再反推 D²L
设计入口永远是那几个额定值：额定功率 Pn、额定相电压 Un、额定转速 nn、极对数 p。第一步不是直接套功率公式，而是先算电磁转矩。转矩才决定电机有效体积。同样功率下，3000 转和 30000 转的电机体积完全不同，用功率推尺寸必然出错。
算转矩之前要先估一个效率和功率因数，把额定功率放大成电磁功率。电磁功率除以机械角速度得到电磁转矩，然后利用转矩与尺寸的核心关系：
T_em = C_T · D² · L · A · B_g
这套源代码里，D²L 才是真正要拿下的目标值。A 是线负荷，B_g 是气隙磁密平均值，C_T 是综合系数。C_T 和绕组分布、磁密波形都有关系，工程上取 0.45~0.65 都可以，最好用自家图纸反推标定。
拿到 D²L 之后，再给一个长径比 λ = L / D。小电机取 0.8~1.2，高速电机可以拉到 2 以上。定子内径和铁芯长度就从 D²L 和 λ 唯一解出来。这一步在脚本里只有三行，但系数 C_T 和 λ 的选择决定了整个方案的性格。
2.2 磁路校核：气隙、齿、轭三处磁密不能只看一个
尺寸定下来，紧接着就是磁路。磁路校核的核心对象是每极磁通：
Φ = α_i · τ · L · B_g
τ 是极距，α_i 是极弧系数，按磁钢极弧换算，一般取 0.6~0.75。拿到每极磁通后，磁通走的三条路要分别校核：气隙、定子齿、定子轭。
气隙磁密是设计初值，直接给。齿部和轭部则需要通过几何尺寸来验算。齿磁密 B_t = Φ / A_t，A_t 是每极下齿部总截面积；轭磁密 B_c = Φ / (2A_c)，磁通沿轭部走一圈，左右两半各承担一半。
齿磁密一般按不超过 1.5~1.7T 控制，轭磁密按 1.3~1.5T 控制。超过硅钢饱和点，励磁电流会飙升，功率因数迅速变差。脚本的做法是先按限值反推最小齿宽和最小轭高，而不是给一个初值再盲调。反推出来的尺寸如果和冲片图对不上，就改槽形或者改磁负荷，这是方案迭代的正常循环。
BH 曲线在这个阶段就要介入了。齿磁密和轭磁密确定后，查硅钢 BH 曲线得到对应磁场强度 H，在 MATLAB 里用 interp1 完成。这里有个口径问题：脚本用平均值磁通核算，齿轭磁密后续和峰值磁密对表时要乘一个波形系数，注释里不写清楚迟早会被自己坑到。
2.3 绕组与槽满率：决定这版方案能不能“绕得下”
主尺寸和磁密过了，绕组才是工艺可行性的关键。脚本里要连续算四个数：每极每相槽数 q、每相串联匝数 N_ph、每槽导体数 N_s、槽满率 SF。
q = Z_s / (2·p·m)，m 是相数，三相取 3。整数槽常见 q=3、4，现在分数槽也很普遍。每相串联匝数由反电动势公式决定：
N_ph = U_n / (4.44 · f · Φ · K_dp)
4.44 就是 2π/√2 的简化写法，f 是电频率。K_dp 是绕组系数初值，整数槽双层绕组可以给 0.92~0.96，分数槽必须重算，这是后话。
每槽导体数再代入 N_s = 2·m·a·N_ph / Z_s，a 是并联支路数。双层绕组 N_s 必须取偶数，因为每槽上下两个线圈边各分一半导体。线径按电流密度定，槽满率则是绕组能不能塞进槽的硬指标。散绕组槽满率做到 0.38~0.45 是正常区间，超过 0.5 基本没法批量生产。
这三步串起来，就是一套电机设计脚本的骨架：输入额定值 → 转矩反推尺寸 → 磁路校核 → 绕组和槽满率 → 性能估算。下面直接看可运行的代码。
3. 可运行的电机设计源代码：输入参数与主尺寸、磁路模块
3.1 输入参数区怎么组织：额定值、材料数据、初值分开放
我一般把脚本头部的输入区按四段组织：额定性能、迭代初值、电磁负荷、材料与限值。全部用 SI 单位，长度用米，磁通用韦伯。文件名随意，下面这段可以直接存成 .m 文件跑：
%% 电机电磁方案初算脚本 % 适用：永磁同步/异步电机初步尺寸估算 % 口径：Bg 为气隙磁密平均值，与每极磁通计算保持一致 clear; clc; % ---- 额定性能输入 ---- Pn = 750; % 额定输出功率, W Un = 220; % 额定相电压有效值, V nn = 1500; % 额定转速, rpm p = 2; % 极对数 % ---- 迭代初值 ---- eta = 0.85; % 预估效率, 参考同类电机 pf = 0.90; % 预估功率因数, 永磁同步取 0.9~0.95 % ---- 电磁负荷初值 ---- Bg = 0.65; % 气隙磁密平均值, T A = 18000; % 线负荷, A/m alpha_i = 0.68; % 极弧系数 K_dp = 0.95; % 绕组系数初值, 整数槽双层绕组 j = 4.5; % 电流密度, A/mm^2 % ---- 材料与结构限值 ---- delta = 0.5e-3; % 气隙长度, m B_ds = 1.6; % 定子齿磁密限值, T B_cs = 1.4; % 定子轭磁密限值, T k_Fe = 0.95; % 叠压系数 m = 3; % 相数 a = 1; % 并联支路数 rho_cu = 2.2e-8; % 铜电阻率, 按 100 摄氏度附近取, Ohm*m % 硅钢 BH 曲线, 实测数据优先 B_fe = [0.00 0.50 1.00 1.20 1.40 1.50 1.60 1.70 1.80 1.90 2.00]; H_fe = [0 60 100 180 400 900 1800 3500 7000 14000 28000];分段的作用是让方案排比时尽量少动材料参数。我只改第一段的额定值，电磁负荷和材料限值按项目实际水平固定。eta 和 pf 是迭代初值，第一次跑完把效率计算值回代，一般两三轮就收敛。
Bg 的口径容易引起团队内争吵。这里是平均值口径，和后面每极磁通计算一致；如果直接拿齿磁密峰值去对表，会发现脚本算出来偏低。注释里写清楚口径，比三个月后翻代码时靠猜要省事得多。
BH 曲线建议用供应商实测表替代手册值。手册数据偏理想，饱和点滞后，插值到 1.6T 附近时误差会直接传导到铁耗和功率因数估算。
3.2 主要尺寸计算模块：从电磁转矩反推定子内径和铁芯长度
输入参数区下面就是第一段核心计算。这里用转矩反推 D²L，再用长径比分解成内径和长度：
%% 主要尺寸计算 f = nn * p / 60; % 电频率, Hz n_rps = nn / 60; % 机械转速, 转/秒 P_em = Pn / eta; % 电磁功率, W T_ref = P_em / (2 * pi * n_rps); % 电磁转矩, N*m % 转矩方程反推 D^2*L % T_em = C_T * D^2 * L * A * Bg C_T = 0.5; % 综合系数, 按电机类型取 0.45~0.65 DL2 = T_ref / (C_T * A * Bg); % D^2*L, m^3 % 长径比分解 lambda = 0.9; % 小型永磁电机取 0.8~1.2 Di1 = (DL2 / lambda)^(1/3); % 定子内径, m Lef = lambda * Di1; % 铁芯长度, mC_T 是这套脚本里最“玄学”的系数。它对结果影响很大，但无法从铭牌直接读到。我习惯先用一台相似结构的量产电机反推一份 C_T 存着，设计新电机时优先用反推值，实在没有参考再按 0.5 起步。
lambda 直接决定电机“胖还是瘦”。转轴细、转动惯量敏感的场景取大值，盘式或扁平结构取小值。脚本里先给 0.9 跑通流程，再按实际空间约束修正。
3.3 磁路校核模块：BH 曲线插值与齿轭磁密计算
尺寸出来以后，立刻做磁路校核。这部分的目标不是算一个结果，而是判断给定尺寸和限值下结构件能不能放得下：
%% 磁路校核: 每极磁通与齿轭最小尺寸 tau = pi * Di1 / (2 * p); % 极距, m Phi = alpha_i * tau * Lef * Bg; % 每极磁通, Wb % 先按整数槽选槽数 q = 3; % 每极每相槽数 Zs = 2 * p * m * q; % 定子槽数 L_Fe = Lef * k_Fe; % 铁芯净长, m % 按齿磁密限值反推最小齿宽 bt_min = Phi / (B_ds * L_Fe * Zs / (2 * p)); % 按轭磁密限值反推最小轭高 hc_min = Phi / (2 * B_cs * L_Fe); % BH 曲线插值, 外插需要控制风险 H_t = interp1(B_fe, H_fe, B_ds, 'linear', 'extrap'); fprintf('每极磁通: %.6f Wb\n', Phi); fprintf('最小齿宽: %.1f mm\n', bt_min * 1e3); fprintf('最小轭高: %.1f mm\n', hc_min * 1e3); fprintf('齿部 1.6T 对应 H: %.0f A/m\n', H_t);核心逻辑是先给磁密限值，反推几何尺寸下限。齿宽和轭高出来后，和实际冲片图做比对：如果槽数定了、槽形画出来齿宽比下限小，磁密就会超标，需要加外径或减小每极磁通。
interp1 加了 'extrap'，但这是双刃剑。BH 表横轴只到 2.0T，外插到 2.2T 也不会报错，但算出来的 H 已经没有物理意义。后面第 5 章会专门说这个坑。磁路校核在脚本阶段做到这一步就够用，齿磁密、轭磁密都压在限值附近，方案才是安全的。
4. 绕组与性能估算模块：让脚本从能算变成能用
4.1 每槽导体数、线径与槽满率校核
磁路过后就是绕组。磁路决定方案能不能成立，绕组决定方案能不能绕得出来。这一段代码把每槽导体数、线径和槽满率一次算完：
%% 绕组计算 % 由相电压反推每相串联匝数 N_ph = Un / (4.44 * f * Phi * K_dp); % 换算每槽导体数 Ns = N_ph * 2 * m * a / Zs; Ns = ceil(Ns); % 向上取整 if mod(Ns, 2) ~= 0 % 双层绕组必须取偶数 Ns = Ns + 1; end % 用圆整后的 Ns 修正每相串联匝数 N_ph = Ns * Zs / (2 * m * a); % 线径估算 I_ph = Pn / (3 * eta * Un * pf); % 相电流, A, 注意是 3 相 S_cu = I_ph / (a * j); % 每根导体截面积, mm^2 d_cu = sqrt(4 * S_cu / pi); % 裸铜线径, mm d_out = d_cu * 1.06; % 漆包线外径, mm % 槽满率校核 S_slot = 60; % 槽毛面积, mm^2, 冲片图量取 S_eff_in = 0.75; % 槽有效面积系数 SF = Ns * pi * (d_out / 2)^2 / (S_slot * S_eff_in); fprintf('每相串联匝数: %.1f\n', N_ph); fprintf('每槽导体数: %d\n', Ns); fprintf('裸铜线径: %.2f mm, 漆包外径: %.2f mm\n', d_cu, d_out); fprintf('槽满率: %.1f%%\n', SF * 100);每相串联匝数公式里的 4.44 是 2π/√2，反电动势要和额定电压匹配。这里按额定相电压设计，意味着额定点反电动势接近电压值，电流和转矩才能对齐。
相电流公式特别容易少乘一个 3。Un 是相电压时，三相总功率 Pn = 3·U_ph·I·eta·pf，所以 I_ph 必须除以 3；如果用的是线电压，则用根号 3。两种口径差出三倍，效率估算直接崩掉。
线径那里的 1.06 是按漆膜厚度粗略放大。小线径漆膜占比更高，0.2mm 以下的线最好查线规表精确算。槽面积 S_slot 这里是示例值，实际从冲片图量取；槽有效面积系数 0.75 已经在扣绝缘、槽楔和制造公差，不要再额外打折。
4.2 铜耗、铁耗与效率估算
绕组算完，最后一步是性能初估。损耗分铜耗、铁耗、机械损耗和杂散损耗，依次叠加：
%% 损耗估算 % 每相电阻 l_turn = 2 * (Lef + 0.03); % 平均匝长, m, 端部按 30mm 粗估 R_ph = rho_cu * N_ph * l_turn / S_cu; % 每相电阻, Ohm P_cu = 3 * I_ph^2 * R_ph; % 铜耗, W % 铁耗: 比损耗粗估, 系数按实际硅钢牌号拟合 G_fe = 2.5; % 硅钢质量初值, kg k_fe = 2.2; % 1T/50Hz 比损耗系数 P_fe = k_fe * (f / 50)^1.3 * Bg^2 * G_fe; % 机械与杂散损耗按经验比例 P_mec = 0.01 * Pn; P_stray = 0.005 * Pn; P_loss = P_cu + P_fe + P_mec + P_stray; eta_calc = Pn / (Pn + P_loss); fprintf('铜耗: %.1f W\n', P_cu); fprintf('铁耗: %.1f W\n', P_fe); fprintf('估算效率: %.1f%%\n', eta_calc * 100);每相电阻的计算逻辑是：N_ph 匝串联，每匝平均长度 l_turn，截面积 S_cu（单位已经是 mm²，MATLAB 里 R_ph = rho * length / area 直接算）。rho_cu 取 100 摄氏度附近的值，更接近实际运行工况。
铁耗是整段估算误差最大的地方。脚本用比损耗 × 铁重粗估，k_fe 需要按具体硅钢牌号拟合，供应商给出的 1T/50Hz 比损耗直接替换即可。铁重 G_fe 在几何不完整时按同类电机估，等冲片图出来再回代，误差能控制在 15% 以内。
效率计算值如果有 5 个百分点的误差很正常，脚本的价值在于趋势判断：改一档铁芯长度，效率往哪个方向走、槽满率是否还能接受，这个趋势是可信的。
4.3 结果输出与绘图：设计评审看什么
算完一批参数后，我习惯统一打印成一张方案卡，方便横向比较。上面的 fprintf 已经输出主要结果，这里再补一个简单的方案判定：
%% 方案合理性检查 check_pass = true; if SF > 0.50 fprintf('[警告] 槽满率 %.1f%% 超过 50%%, 绕线有难度\n', SF * 100); check_pass = false; end if B_ds > 1.7 fprintf('[警告] 齿磁密 %.2fT 偏高, 硅钢可能饱和\n', B_ds); check_pass = false; end if eta_calc < 0.7 fprintf('[警告] 效率低于 70%%, 检查输入初值\n', eta_calc); check_pass = false; end这组输出里我最看重槽满率和磁密超限警告。效率初估可以后期标定，但槽满率超 0.5 是工艺硬伤，磁密超饱和是原理硬伤，两个都不该拖到有限元阶段才发现。
5. 电机设计脚本的常见坑与排查：五条血泪经验
5.1 BH 曲线插值返回 NaN，整条脚本崩掉
现象：改小一档功率后，脚本运行到磁路校核直接报错，变量区出现 NaN，后面的铁耗和效率全部失去意义。
原因：齿磁密或轭磁密算到 1.7T 以上，超出了 BH 表横轴范围。interp1 默认不处理外插，查不到值就返回 NaN。
解决：所有 BH 插值统一走一个封装函数，内部做边界判断。更稳妥的做法是在磁密超限时直接报错，提醒你调整尺寸或磁负荷：
function H = getHfromBH(B, B_fe, H_fe) if B > max(B_fe) error('磁密 %.2fT 超出 BH 表范围, 请检查磁路设计\n', B); end H = interp1(B_fe, H_fe, B, 'linear'); end5.2 槽满率一算就是 100% 以上
现象：第一版方案算出来槽满率 1.2，线径粗得根本塞不进槽，但绕组匝数看起来又合理。
原因：槽面积用了毛面积，没乘有效面积系数；线径用了裸铜径，没算漆包线外径；线规选的太粗也没有换并联支路数去调整。
解决：槽满率必须用漆包线外径和有效槽面积计算，SF = N_s·d_out²/(S_slot·S_eff_in)。有效面积系数取 0.7~0.8，先扣绝缘和工艺空隙。如果 SF 还是超，优先增加并联支路数而不是盲目缩小线径。
5.3 效率算出来是负的，或者超过 100%
现象：脚本跑完打印估算效率 -20%，检查公式没发现明显错误。
原因：相电流公式口径错了。Un 是相电压时，I_ph = Pn / (3·eta·Un·pf)，少乘一个 3 会让电流大三倍，铜耗算出来比输入功率还大，效率自然成负值。
解决：统一电压口径并写进注释。相电压配 3，线电压配根号 3，两者不要混用。效率出现超过 100% 或负值，先查这一行，九个 case 里八个都是这里的问题。
5.4 只改一个参数，结果全乱
现象：改完铁芯长度，磁密变了但槽满率没变；改回原来的值，结果却和第一版对不上。
原因：脚本按分区执行，工作区残留了旧变量。某个模块没从头跑，Zs 改过但绕组模块没重新执行，旧 N_s 和新 Z_s 混在一起，输出自然错乱。
解决：脚本阶段每次修改后从顶部 clear 全部重跑。要做参数扫描，必须把三个模块封装成函数，输入输出明确，不依赖工作区状态。这也是后面第 6 章的基础。
5.5 分数槽电机绕组系数还按整数槽填，反电动势差很多
现象：极槽配合改成 30 槽 4 极后，K_dp 还填 0.95，有限元一算反电动势差 8%。
原因：q = 30/(4×3) = 2.5，不是整数槽，分布系数需要按分数槽重新推导，短距系数也要按实际线圈跨距算。直接把 0.95 填进去，磁路和绕组完全对不上。
解决：把绕组系数当设计变量，在脚本里按 q 计算分布系数，再乘短距系数：
q = Zs / (2 * p * m); alpha = pi / (m * q); % 槽距电角 Kd = sin(q * alpha / 2) / (q * sin(alpha / 2)); % 短距系数 Ky 按线圈跨距 y 计算 % K_dp = Kd * Ky分数槽的绕组系数值不能凭经验填，每次极槽配合变化都要从公式推一遍。
6. 让脚本从单点计算变成方案工具：参数扫描与有限元对账
6.1 把计算模块封装成函数，一页循环扫完一版方案
单点算完只是开始。真正让这套源代码产生价值的，是把前几章的模块封装成函数，然后做参数扫描。我习惯先定义这样一个入口函数：
function [Di1, Phi, Ns, SF, eta_calc] = em_design(Pn, Un, nn, p, A, Bg, Lef_target) % 输入: 额定功率/电压/转速/极对数/线负荷/磁密/铁芯长度 % 输出: 内径/每极磁通/每槽导体数/槽满率/效率 % 内部复用第 3、4 章的计算模块 f = nn * p / 60; n_rps = nn / 60; P_em = Pn / eta; T_ref = P_em / (2 * pi * n_rps); DL2 = T_ref / (0.5 * A * Bg); Di1 = (DL2 / Lef_target)^(1/3); % ... 后续计算照搬前面模块 end封装好之后，一页循环就能扫完一组方案：
Lef_list = linspace(0.06, 0.12, 7); res = []; for i = 1:length(Lef_list) [Di1_i, Phi_i, Ns_i, SF_i, eta_i] = ... em_design(750, 220, 1500, 2, 18000, 0.65, Lef_list(i)); res(i,:) = [Di1_i*1e3, Phi_i, Ns_i, SF_i, eta_i]; end plot(Lef_list*1e3, res(:,5)*100); xlabel('铁芯长度/mm'); ylabel('估算效率/%');扫描的重点不是找什么“最优解”，而是看拐点。效率随铁芯长度增加趋缓、槽满率随匝数变化的方向、磁密是否在某一段长度下突破限值，这些都是比单点数值更可靠的设计信息。
6.2 空载反电动势对账：脚本和有限元能不能对上
脚本出方案后，一定要和有限元对一次账。最值得对齐的量是空载反电动势，它是磁路的直接体现，比转矩更好对标：
E_0 = 4.44 · f · N_ph · K_dp · Φ
操作方法是把额定转速、空载工况放进有限元，读出相反电动势波形基波幅值，再和脚本的 E_0 对比。偏差在 5% 以内，说明极弧系数、每极磁通、绕组系数这套模型是自洽的。偏差超过 10%，不要急着调有限元，回头查极弧系数取值和磁钢工作点。
我在实际项目里养成的习惯是：脚本里每个常数都写上来源和适用口径，不然三个月后回看代码，谁也说不清那个 0.5 是怎么来的。脚本的强项是快，但再快也替代不了样机标定；用一次有限元或一次实验数据把脚本标定过，它才真正变成你自己的设计工具。希望帮到你。
本文还有配套的精品资源，点击获取
