三相异步电机dq模型Simulink仿真：从方程到曲线-育师

简介：这份三相异步电机本体模型面向电机原理学习、本体设计分析与Simulink仿真入门者，核心是用数学公式搭建的成品模型，可用于观察启动、稳态运行及不同工况下的电流、电压、转矩、转速等波形，解决缺少可修改参考模型、难以直观验证电机特性的问题。压缩包共14个文件，约547KB，包含6个txt说明、6张jpg波形或结构图、1个html页面和1个doc文档，便于先读说明再对照图形理解模型结构、参数含义与各文件用途。已有249人学习下载，说明模型在同类资源中具备一定参考价值与复用基础。读者可据此熟悉三相异步电机的建模思路与参数设置，对照电流、电压、转矩、转速等波形检查运行性能，并在此基础上调整负载、修改参数或扩展故障工况，作为课程学习、课题预研和电机控制仿真的起点，也为后续课程设计与科研仿真提供可直接修改的基础框架。

1. 从一台转不起来的电机说起：这套模型要解决的三个现实问题

做三相异步电机本体模型，绕不开的一个动作是把定子、转子的电压方程、磁链方程、转矩方程和运动方程一行行敲进 Matlab 的 Simulink，再用积分器把它们串成一台会转的电机。这件事听起来像电气工程课设的标准流程，真动手才知道坑在哪里：同一组参数，有人搭出来 0.3 秒就爬到 1470 r/min，有人搭出来第一个采样点就是 NaN。差别往往不在公式抄错，而在坐标系选型、积分初值和求解器步长这三件小事上。

这套成品模型把四阶 dq 电压方程、坐标变换、电磁转矩与机械运动拆成彼此独立的子系统，参数集中在同一个脚本里改，电气量和机械量分开观测。它适合三类人：正在做电机拖动课设、需要一份能跑出曲线的东西交差的学生；做传动方向毕设、要给矢量控制算法准备一个被控对象的同学；以及手上有个变频器控制策略、想先拿理想电机验证逻辑的从业者。下面按「数学骨架 → 总体架构 → 关键实现 → 排查 → 验收」的顺序，把这份模型从公式到曲线的每一步说清楚。

2. 数学骨架：从三相时变方程压到四阶 dq 模型

2.1 三相静止坐标系下的电压与磁链方程长什么样

先从最原始的三相 abc 坐标系写起。定子三相绕组各自满足电压方程，转子如果是鼠笼型，三相端电压恒为零：

u_a = R_s·i_a + dψ_a/dt u_b = R_s·i_b + dψ_b/dt u_c = R_s·i_c + dψ_c/dt 0 = R_r·i_ar + dψ_ar/dt 0 = R_r·i_br + dψ_br/dt 0 = R_r·i_cr + dψ_cr/dt

磁链方程才是麻烦的来源。定子自感是常数，L_aa = L_ls + L_ms；定子互感在理想三相对称条件下近似为-L_ms/2；真正随位置变化的是定转子之间的互感，它和转子电角度 θ_r 呈余弦关系，写成L_ar = L_m·cos(θ_r)，另外两相则分别是cos(θ_r - 120°)和cos(θ_r + 120°)。转子绕组对定子的互感同理，只是角度反过来。把这六组方程展开，得到的是一个 6×6 的电感矩阵，矩阵里至少有一半元素是 θ_r 的三角函数。

如果直接在 Simulink 里照这个形式搭，模型会变成什么样子？六个积分器、十几个 Trigonometric Function 模块、几十个乘加环节，每一拍都要重新计算电感矩阵的每个元素，然后再做一次 6×6 的矩阵与向量相乘。连线一多，查找信号源就成了体力活；更麻烦的是这种时变系数结构天生容易引出代数环——电感值依赖角度，角度由转速积分而来，而转速又依赖电磁转矩，转矩又依赖电流，电流的导数里还带着电感。仿真跑起来慢、改一个参数要重算一堆系数，都是这套结构的副作用。

2.2 Clarke 与 Park 变换：把时变电感变成常数矩阵

工程上的处理办法是把三相静止坐标先压到两相静止坐标，再旋转到同步坐标系。第一步 Clarke 变换采用等幅值形式，三相三线制下零序分量为零，可以略去：

i_α = (2/3)·(i_a - 0.5·i_b - 0.5·i_c) i_β = (2/3)·(√3/2·i_b - √3/2·i_c)

第二步 Park 变换，旋转角 θ_e 取同步电角速度的积分：

i_d = i_α·cos(θ_e) + i_β·sin(θ_e) i_q = -i_α·sin(θ_e) + i_β·cos(θ_e)

反变换就是把矩阵转置回去。放到 Simulink 里，这段用 MATLAB Function 模块写最省事，代码和公式一一对应，评审的人一眼能看懂：

function [i_dq, i_abc] = abc_dq(i_abc, th_e) % 三相静止坐标与同步旋转坐标的双向变换（等幅值原则） % 输入：i_abc 三相电流/电压（A 或 V），th_e 同步电角度（rad） % 输出：i_dq = [d 轴分量; q 轴分量]，i_abc 为反变换结果 C = 2/3 * [1, -1/2, -1/2; 0, sqrt(3)/2, -sqrt(3)/2]; % Clarke，等幅值 ial = C * i_abc(:); Pk = [cos(th_e), sin(th_e); -sin(th_e), cos(th_e)]; % Park i_dq = Pk * ial; i_abc = (Pk \ ial) * ones(1,3) .* [1;1;1]'; % 反变换占位，实际按需展开 end

这里有一个高频踩坑点：Clarke 变换的系数和转矩公式的系数必须配套。用等幅值变换（系数 2/3），dq 轴电流的幅值等于三相电流的幅值，电磁转矩公式前面要乘 1.5；用等功率变换（系数 √(2/3)），转矩公式前面就不带 1.5。两套体系混用，转矩会差 1.5 倍，转速稳态值看着「差不多」但负载一加就露馅。我一般会在参数脚本里留一个P.ktr常量专门记这个系数，避免后面改变换时忘了同步改转矩。

变换带来的收益很直接：电感矩阵从随角度变化的 6×6 变成常数矩阵，磁链和电流之间是线性关系，整个电机被压缩成四个状态变量（i_ds、i_qs、i_dr、i_qr）或者四个磁链状态，模型规模砍掉一半以上，仿真速度也跟着上来。

2.3 转矩方程与运动方程：两条积分链怎么接上

压到 dq 坐标以后，电压方程和磁链方程都变成定常系数：

u_ds = R_s·i_ds + dψ_ds/dt - ω_e·ψ_qs u_qs = R_s·i_qs + dψ_qs/dt + ω_e·ψ_ds 0 = R_r·i_dr + dψ_dr/dt - (ω_e - ω_r)·ψ_qr 0 = R_r·i_qr + dψ_qr/dt + (ω_e - ω_r)·ψ_dr ψ_ds = L_s·i_ds + L_m·i_dr ψ_qs = L_s·i_qs + L_m·i_qr ψ_dr = L_r·i_dr + L_m·i_ds ψ_qr = L_r·i_qr + L_m·i_qs

其中L_s = L_ls + L_m、L_r = L_lr + L_m，ω_e 是同步电角速度，ω_r = p·ω_m 是转子电角速度（p 为极对数，ω_m 为机械角速度）。注意转子那两个方程右端是 0，这是鼠笼短接的物理事实，不是漏写——很多模型在这里给转子留了个电压输入端口，接上理想电压源，电机就悄悄变成了绕线转子电机，空载电流和起动转矩全对不上。

电磁转矩用磁链和电流的叉乘形式最省事，也最容易和前面的变换系数对齐：

T_e = 1.5·p·(ψ_ds·i_qs - ψ_qs·i_ds)

机械侧是标准的转动方程，加上粘性阻尼项更贴近实际：

J·dω_m/dt = T_e - T_L - B·ω_m dθ_m/dt = ω_m θ_e = p·θ_m

这四组方程连起来就是一条闭环积分链：定子电压 → 电流导数 → 电流积分 → 磁链 → 转矩 → 转速 → 角度 → 坐标变换角 → 回到电压方程里的旋转项。链上有两处积分器：一处是电流（或者磁链），一处是转速。仿真能不能收敛、波形干不干净，八成取决于这两个积分器的初值、限幅和求解器步长设得对不对。

3. 总体架构：五个子系统与一条不产生代数环的信号链

3.1 顶层模块划分与端口约定

顶层不追求「一个模块搞定一切」，而是按物理角色切成五块，每块单独封装、单独测。这样做的直接好处是：出问题时能快速定位到是电磁部分算错还是机械部分接反，而不是对着一坨连线猜。

信号方向的约定要提前定死：电磁本体接收的是定子 dq 电压，机械运动接收的是转矩和负载转矩，坐标变换接收的是电角度。端口名建议直接带上单位，比如uds_V、wm_radps、Te_Nm。看着啰嗦，但两个人对接模型的时候，能省掉一轮「你这个 ω 是机械角还是电角」的扯皮。转差率单独算一路：s = (ω_e - ω_r)/ω_e，接到示波器上和转速并排看，判断稳态点比看转速直观。

3.2 建模方式选型：基础模块、MATLAB Function 还是 S-Function

同一个方程有四种搭法，选哪种取决于你要交付什么。

| 方案 | 搭建成本 | 可读性 | 仿真速度 | 适用场景 | | 基础模块拼乘加 | 高 | 低 | 中 | 课堂演示、需要展示每个运算环节 | | MATLAB Function | 低 | 高 | 中 | 本体模型首选，公式即代码 | | S-Function（C） | 高 | 中 | 高 | 大批量参数扫描、要生成嵌入式代码 | | Simscape Electrical 现成模块 | 极低 | 低 | 慢 | 系统级联调，不关心内部方程 |

这份模型走的是 MATLAB Function 路线。理由很实在：四个状态方程写成矩阵形式只有几行，注释能直接标上物理量，导师或者评审翻代码时不用逆着连线找信号；改方程就是改一行代码，不用在画布上挪模块。代价是每次仿真步进都有一次函数调用开销，参数得通过结构体传进去，不能用工作区里的裸变量直接引用——这一点在下一节展开。

S-Function 不是不能用，如果后面要拿这套模型去做几百组参数的批量扫描，把方程固化成 C Mex 能省不少时间。但对课设和毕设这个量级，MATLAB Function 的调试便利性远比那点速度重要。

3.3 参数脚本与工作区变量：一次改参数，全模型生效

模型里最忌讳的事情是把 1.85、0.21 这种数字直接写进模块对话框。改一台电机就得翻遍整个画布，漏改一个漏感，转矩曲线就会莫名偏低。正确做法是所有参数集中在一个初始化脚本里：

% init_motor.m —— 三相异步电机参数集中初始化 clear P; P.fn = 50; % 额定频率 Hz P.p = 2; % 极对数 P.Un = 380; % 线电压有效值 V P.Upk = P.Un * sqrt(2) / sqrt(3); % 相电压峰值 V（约 310 V） P.Rs = 1.85; % 定子电阻 Ω P.Rr = 1.65; % 转子电阻折算值 Ω P.Lls = 0.01; % 定子漏感 H P.Llr = 0.01; % 转子漏感 H P.Lm = 0.20; % 励磁电感 H P.Ls = P.Lls + P.Lm; % 定子自感 H P.Lr = P.Llr + P.Lm; % 转子自感 H P.J = 0.03; % 转动惯量 kg·m^2 P.B = 1e-3; % 粘性阻尼 N·m·s P.we = 2*pi*P.fn; % 同步电角速度 rad/s P.Ts = 1e-4; % 求解器最大步长 s P.tend = 1.5; % 仿真时长 s assignin('base', 'P', P); % 写入基础工作区，供模型回调引用

这段脚本有几个细节值得说。P.Upk单独算出来，是因为很多人把 380 V 直接当相电压峰值塞进u_ds，结果磁链偏大、电流偏大、转速却看不出明显异常，属于那种「曲线看着像对的、数值全错」的隐蔽 bug。P.Ts和P.tend放进参数结构体，是为了让批量脚本能覆盖它们。

脚本的调用位置建议放在模型的InitFcn回调里（Model Properties → Callbacks → InitFcn 填init_motor），这样每次打开模型或者点仿真，参数都会自动刷新，不会出现「改了脚本没生效」的错觉。用assignin写到 base 工作区而不是 Model Workspace，是为了让 MATLAB Function 内部的P参数能正常解析，同时也方便在命令行里直接查P.Ls。

4. 关键实现：磁链、转矩、运动三段的代码与参数落地

4.1 电磁本体：四阶状态方程的矩阵写法与积分器初值

电磁本体是整个模型的心脏。把 2.3 节的四组方程写成矩阵形式，代码量能压到十几行，而且不容易抄错符号：

function [di, psi, Te] = motor_dq(u_dq, wm, i_dq, P) % 三相异步电机 dq 轴本体方程（电动机惯例，同步旋转坐标系） % 输入：u_dq = [u_ds; u_qs] 定子电压，V % wm 转子机械角速度，rad/s % i_dq = [i_ds; i_qs; i_dr; i_qr] 四阶电流状态，A % 输出：di 电流导数，A/s % psi 四阶磁链，Wb % Te 电磁转矩，N·m R = diag([P.Rs, P.Rs, P.Rr, P.Rr]); L = [P.Ls, 0, P.Lm, 0; 0, P.Ls, 0, P.Lm; P.Lm, 0, P.Lr, 0; 0, P.Lm, 0, P.Lr]; psi = L * i_dq; % 磁链 = 电感矩阵 × 电流 we = P.we; % 同步电角速度 wr = P.p * wm; % 转子电角速度 W = [0, -we, 0, 0; we, 0, 0, 0; 0, 0, 0, -(we - wr); 0, 0, (we - wr), 0]; % 旋转电动势的反对称矩阵 u = [u_dq(1); u_dq(2); 0; 0]; % 鼠笼转子端电压恒为 0 di = L \ (u - R*i_dq - W*psi); % 状态方程：di/dt = L⁻¹(u - Ri - Wψ) Te = 1.5 * P.p * (psi(1)*i_dq(2) - psi(2)*i_dq(1)); end

代码里的W矩阵对应电压方程中的旋转项，第一块[0 -we; we 0]乘上[ψ_ds; ψ_qs]正好得到-we·ψ_qs和we·ψ_ds；第二块的旋转频率换成转差角频率we - wr，对应转子方程。符号一旦搞反，最典型的症状是空载时电流一直往上爬、转速却起不来，或者电机一开始就工作在发电状态。

L \ (...)每步都做一次 4×4 求解。电感矩阵是定常的，可以把inv(L)提前算好放进参数结构体，运行时改成矩阵乘法，能省一点时间；但对这个规模的模型，差别肉眼可忽略，我更倾向保留\的写法，物理意义更清楚。

积分器的设置比代码本身更影响结果。四个电流状态各配一个 Integrator 模块（或者用 State-Space 的向量形式），初值一律给 0——转子静止、定子未通电，电流为零是唯一合理的初始条件。积分器上限建议给到 ±50 A 左右，防止数值发散时直接冲成 Inf 导致整个仿真卡死；下限对应给负值。这个限幅不是为了物理正确，纯粹是给调试留一颗后悔药：波形被削平了你立刻知道出事了，而不是等到仿真报错才回头找原因。

4.2 转矩与机械运动：负载阶跃怎么加、角度怎么防溢出

转矩计算已经在电磁本体的函数里输出了，单独拆一个 Gain 模块出来接 Scope 也合理，方便答辩时指着曲线讲。真正需要单独搭的是机械运动部分，因为它要接负载：

function [dwm, dth] = mech_dynamics(Te, TL, wm, P) % 机械运动方程：J·dωm/dt = Te - TL - B·ωm % 输入：Te 电磁转矩 N·m，TL 负载转矩 N·m，wm 机械角速度 rad/s % 输出：dwm 角加速度 rad/s²，dth 机械角速度（供积分得角度） dwm = (Te - TL - P.B * wm) / P.J; dth = wm; end

负载转矩TL用 Signal Builder 或者 Step 模块给：0 到 0.6 s 给 0，0.6 s 阶跃到 10 N·m，就能复现「空载启动—突加负载」这个最经典的工况。这里有个容易忽略的点：TL的符号必须和Te在同一套正方向约定下。电动机惯例下，电磁转矩为正驱动，负载转矩为正阻碍，两者相减。如果负载模块里手滑给了负号，电机会越转越快，稳态转速冲到同步转速以上，看起来像「永动机」，实际是符号错了。

角度积分那条链要单独处理。θ_m一直积分下去，几分钟的仿真时长就能累积到几千弧度，sin、cos对超大角度的计算精度会下降，严重时坐标变换会出现肉眼可见的抖动。处理办法是在角度积分器后面接一个mod运算，把机械角折算到[0, 2π)再乘极对数：

th_e = P.p * mod(th_m, 2*pi); % 电角度折算，避免大角度三角函数精度损失

mod模块放在 MATLAB Function 里或者直接用 Simulink 的 Math Function（选 mod）都行。这一步加上以后，长时间仿真（比如 10 s 以上）的波形会明显干净。

4.3 参数化封装与批量仿真：一组脚本跑完四种工况

模型调通以后，往往要跑好几组工况出图：空载启动、突加负载、额定负载、堵转。手动改一个参数点一次仿真，四组下来半小时就没了，还容易记混哪张图对应哪组设置。用Simulink.SimulationInput批处理，一次性跑完并归档：

% run_cases.m —— 一次跑完四种工况，结果存到 results 目录 cases = struct( ... 'name', {'no_load', 'step_load', 'rated', 'locked'}, ... 'tl', {0, 10, 19.1, 0}, ... 'tend', {1.0, 1.5, 1.5, 0.2}); for k = 1:numel(cases) in = Simulink.SimulationInput('asyn_motor_dq'); in = in.setVariable('TL_cfg', cases(k).tl); % 覆盖负载转矩 in = in.setModelParameter('StopTime', num2str(cases(k).tend)); in = in.setModelParameter('SaveState', 'on', ... 'StateSaveName', 'xFinal'); % 存末态，便于续算 out = sim(in); save(fullfile('results', [cases(k).name '.mat']), 'out'); end

SimulationInput相比老式的set_param有两个实惠：一是参数覆盖只作用于这一次仿真，不会污染模型文件，跑完打开模型还是原始值；二是天然支持parsim并行，参数组多的时候能把时间压下来。SaveState打开是为了某些需要从稳态续算的场景，比如先空载跑到稳态，再把末态作为堵转工况的初值，省掉一段过渡过程。

仿真配置本身也要固定下来，别每次靠手点。在模型回调或者脚本里统一设一遍：

set_param('asyn_motor_dq', 'Solver', 'ode23tb', ... 'MaxStep', '1e-4', 'RelTol', '1e-4', 'AbsTol', '1e-6', ... 'StopTime', '1.5', 'SaveOutput', 'on', 'SaveFormat', 'StructureWithTime');

ode23tb适合这种含刚性环节的方程（定子时间常数比机械时间常数小两个数量级），变步长配MaxStep = 1e-4，既能抓住开关细节又不至于慢得离谱。相对容差放到 1e-4 是经验值，再松波形开始起毛刺，再紧仿真时间翻倍，收益不明显。

5. 排查手册：模型发散、转矩反号、稳态偏差的五类现场问题

模型搭完第一次点仿真，十有八九跑不出教科书上那条漂亮的曲线。下面这五类问题出现频率最高，基本覆盖了从「一运行就崩」到「看着对但数值不对」的全部区间。

5.1 启动第一拍就报 NaN 或者直接发散

现象是仿真刚跑几十微秒就报积分失败，示波器上电流一瞬间冲到 1e6 量级然后变成 NaN，或者转速在 0.01 s 内飞到几万转。

根因通常有两个。一是代数环：转矩依赖电流，转速依赖转矩，而电流方程里的旋转项又依赖转速，如果这几段之间全部用直通连接、没有经过任何积分器或延时，Simulink 会在同一时刻反复迭代求解，数值一旦超出范围就炸。二是积分器初值和电压初值冲突：定子电压第一拍就给了 310 V 峰值，电流初值是 0，而电感又很小，di/dt在第一个步长里极大，步长稍大就会越界。

处理办法是把反馈链断开一拍。转矩计算使用上一时刻的电流状态（电流积分器的输出本身就滞后一个步长，只要不把di直接引到转矩计算里，一般不会成环）；如果确实存在直通路径，在中间插一个 Unit Delay 或者 Memory 模块，代价是引入一个采样周期的延迟，对这个模型完全可接受。另一手是给电流积分器加限幅（±50 A 量级），并且在仿真开始的 1 ms 内把定子电压做斜坡上升，从 0 爬到额定峰值，避开第一拍的数值冲击。斜坡用 Ramp 模块或者 Saturation 加 Gain 都能做，加完以后启动电流的波形反而更接近真实软起。

5.2 转矩方向反了，电机一启动就工作在发电状态

现象是转速积分出来是负的，或者转速为正但电磁转矩长期为负，转差率算出负数，定子电流相位和电压对不上。

根因基本是 dq 轴的定义和转矩公式不配套。转矩公式T_e = 1.5·p·(ψ_ds·i_qs - ψ_qs·i_ds)对应的是 q 轴超前 d 轴 90° 的旋转方向；如果 Park 变换里把i_q = i_α·sinθ + i_β·cosθ的符号写反，等于把 q 轴翻了个方向，转矩符号跟着反。另一种情况是极对数 p 漏乘：转子电角速度写成wr = wm，转矩公式里又乘了 p，两边不一致，低速时看不出来，接近同步速时转差算错，转矩突变。

排查办法不复杂：给一组直流量手动验算。令i_ds = 5 A、i_qs = 5 A、i_dr = i_qr = 0，手算ψ_ds = L_s·5、ψ_qs = L_s·5，代进转矩公式得到T_e = 1.5·p·(L_s·25 - L_s·25) = 0。这个结果说明纯 dq 电流不产生转矩（因为磁链和电流同相位），再把i_qs改成 5、i_ds改成 0，转矩应该为正。用两组手算值去卡代码输出，符号对不对立刻见分晓。

5.3 稳态转速停在同步转速以上，或者转差率明显偏大

现象是空载时转速显示 1500 r/min 甚至更高，加载后转速掉得离谱，转差率算出来 10% 以上。

前者多半是负载转矩符号写反，把T_L当成了驱动力，或者阻尼项B·ω_m前面多了个负号，机械方程变成J·dω_m/dt = T_e + T_L，电机自然一路加速。检查方法很直接：把T_L强制置零，看空载稳态转速是否落在 1470~1490 r/min 区间，落不进去就是机械侧的问题。后者通常是转子电阻或者漏感填错，特别是把折算到定子侧的转子电阻当成实际值再折一次，阻值偏大，转差自然偏大。参数表里Rr一定要注明「折算值」，避免二次折算。

还有一种隐蔽情况：励磁电感Lm填得比Ls还大。物理上Lm必须小于L_s和L_r，因为自感等于漏感加励磁电感。填反了以后漏感系数 σ 变成负数，方程的性质都变了，波形会呈现一种「说不上哪里不对但就是不对」的状态。

5.4 波形上有高频毛刺，负载突变时数值抖动

现象是稳态电流上叠着一层细密的高频振荡，负载阶跃的那一瞬间曲线出现尖刺甚至短时溢出。

这几乎都是求解器步长的问题。定子时间常数L_s/R_s大约 0.1 s，但电流环里还有更快的动态，固定步长设成 1e-3 就会把细节全部略过。改用变步长ode23tb，MaxStep限制到 1e-4，RelTol收到 1e-4，毛刺基本消失。如果还不行，把AbsTol从默认的 auto 改成 1e-6，让求解器在小信号段也不放松精度。

另一个来源是负载阶跃本身。Step 模块在跳变点产生一个数学上的不连续，求解器为了跨过这个点会自动缩小步长，表现为那一拍的计算时间变长、波形出现一个尖峰。这不是模型错，是数值方法的正常反应。真要消掉，把阶跃改成 1 ms 的斜坡，物理上更接近接触器吸合的过程，波形也干净。

5.5 换一台电机参数后结果全变，但公式没动过

现象是把参数脚本里的数值换成另一台电机的铭牌数据，仿真要么跑飞，要么稳态点完全不对。

先查单位。铭牌上给的是功率、电压、电流、转速，模型要的是电阻、电感、惯量。定子电阻 1.85 Ω 和 0.00185 kΩ 差别巨大；电感 0.21 H 和 0.21 mH 差三个数量级。再查电压口径：380 V 是线电压有效值，进 dq 方程前必须转成相电压峰值（约 310 V），直接把 380 塞进去，磁链会偏大 22%，空载电流跟着偏大，稳态转速却可能只差几转，非常容易蒙混过关。

然后是惯量。转动惯量对启动时间的影响是线性的，J填大十倍，启动过程慢十倍，看着像「电机没劲」，实际是惯量给错了。最后查极对数，2 极对和 4 极对的电机，同步转速差一倍，任何转速相关的判据都要跟着改。我一般会在参数脚本末尾加一段自检：

assert(P.Lm < P.Ls && P.Lm < P.Lr, '励磁电感必须小于自感'); assert(P.Rs > 0 && P.Rr > 0, '电阻必须为正'); fprintf('同步转速 %.0f r/min，额定转矩约 %.1f N·m\n', ... 60*P.fn/P.p, 3000/(2*pi*P.fn/P.p));

跑一次脚本看到这两行输出，比跑完仿真再回头找问题省事得多。