简介：一份围绕MATLAB变压器仿真的系统性分析文档，面向电气工程专业学生、电力系统设计人员及MATLAB仿真入门者。内容以电磁感应原理和变压器基本结构为起点，详细梳理了数学模型、等效电路、空载损耗、负载损耗及动态暂态过程，并逐一给出空载、负载、暂态特性曲线分析与故障诊断实例。目录同时包含研究背景与意义、国内外研究现状、MATLAB环境搭建、模型输入参数确定和仿真模型验证修正等模块，逻辑完整，可作为课程设计、毕业设计或工程报告的撰写参考。资源为单篇docx文档，压缩包约60KB，无额外代码或模型文件，胜在知识点集中、结构清晰，适合快速理解变压器仿真建模的整体思路与关键参数计算方法。目前已有54人学习，适合需要系统梳理变压器特性分析框架的读者。
1. 用MATLAB搭一台能算准的变压器，比你想的更依赖模型选择
变压器仿真在电力电子和电力系统里被问得最多的问题不是“怎么搭模型”，而是“为什么我搭的模型算出来的空载电流和铭牌对不上”。问题通常不出在 MATLAB 操作上，而是出在模型层级选错：拿一个只适合稳态分析的简化等效电路去算励磁涌流，或者拿一个只适合工频的线性电感模型去看 DC-DC 变换器里的高频励磁电流，结果自然离谱。这篇博文直接讲我平时做变压器仿真会用的一套路径：先从铭牌数据反推等效电路参数，再用 M 脚本做稳态特性分析，最后落到 Simulink 的 Simscape 物理模型去做暂态校验。整个过程用到的工具只有 MATLAB 基础模块和 Simulink，不需要额外工具箱，适合做课程设计、科研预研，也适合工程上需要快速估算变压器运行特性的场景。读者最好有一点电路理论和 MATLAB 基础，即使没有，按步骤也能跑通。
2. 变压器仿真的第一层：从铭牌反推等效电路参数
2.1 变压器等效电路为什么是仿真的地基
MATLAB 里画变压器模型，本质上是把磁路耦合关系折算成电路。最常见的 T 型等效电路里有励磁支路（并联的激磁电阻 Rc 和激磁电感 Lm）、原副边漏抗（折算后的 R1、R2、L1、L2），四个参数就足以描述一台工频变压器的稳态行为了。关键是这四个值没法直接测，铭牌上只有容量 Sn、电压 U1n/U2n、空载损耗 P0、空载电流 I0%、短路损耗 Pk、短路电压 Uk%，把这些铭牌值换算成等效电路参数就是仿真的第一道工序。
常踩的坑是直接把短路损耗当成绕组铜耗去反推电阻，却忘了短路试验时励磁支路也有损耗。工程上宋任穷忽略励磁支路的影响，但如果要精度，短路试验数据算出的铜耗应该做温度折算，MATLAB 脚本里一般会把 75 摄氏度的电阻值作为输出标准。另一个常见坑是空载电流的百分数换算：铭牌写 I0% = 1.2%，指的是空载电流占额定电流的 1.2%，实际空载电流有效值要乘上额定电流，而不是直接把百分比当成安培数。
我倾向先把标幺值算明白，再用基准值反推有名值。这样做的理由是：标幺值下的励磁阻抗和漏抗一般在固定范围内（励磁阻抗标幺值通常在 20~200 之间，漏抗标幺值通常在 0.04~0.15 之间），一旦参数反推结果超出这个区间，基本可以断定铭牌数据读错或换算出错，这比对着有名值空想靠谱得多。
2.2 MATLAB 脚本实现参数反推的最小代码
以下脚本把铭牌数据转成 MATLAB 变量，再计算出折算到原边的等效电路参数。代码故意精简成线性流程，方便看懂每一步在算什么：
%% 变压器铭牌数据（示例：100kVA，10kV/0.4kV，Yyn0） Sn = 100e3; % 额定容量 VA U1n = 10e3; % 原边额定电压 V U2n = 400; % 副边额定电压 V f = 50; % 额定频率 Hz P0 = 320; % 空载损耗 W I0_pct = 1.8; % 空载电流百分数 % Pk = 1450; % 短路损耗 W（已折算到75℃） Uk_pct = 4.0; % 短路电压百分数 % %% 由铭牌数据计算基值 I1n = Sn / U1n; % 原边额定电流 A Zbase = U1n^2 / Sn; % 原边阻抗基值 Ω I0 = I0_pct / 100 * I1n; % 空载电流有效值 A %% 励磁支路参数（由空载试验折算） Rc = U1n^2 / P0; % 激磁电阻 Ω % 励磁阻抗模值：Zmag = U1n / I0 Zmag = U1n / I0; % 励磁阻抗 = Rc // jXm，解出 Xm（注意不是串联关系） Xm = 1 / sqrt(1/Zmag^2 - 1/Rc^2); Lm = Xm / (2*pi*f); % 励磁电感 H %% 短路阻抗参数（由短路试验折算） Zk_pu = Uk_pct / 100; % 短路阻抗标幺值 Zk = Zk_pu * Zbase; % 短路阻抗有名值 Ω Rk_pu = Pk / Sn; % 短路电阻标幺值（铜耗/容量） Rk = Rk_pu * Zbase; % 短路电阻 Ω % 漏抗：Zk^2 = Rk^2 + Xk^2 Xk = sqrt(Zk^2 - Rk^2); Lk = Xk / (2*pi*f); % 短路电感 H %% 输出折算到原边的 T 型等效电路参数 fprintf('激磁电阻 Rc = %.2f Ω\n', Rc); fprintf('激磁电感 Lm = %.2f H\n', Lm); fprintf('绕组等效电阻 Rk = %.2f Ω\n', Rk); fprintf('绕组等效漏感 Lk = %.2f H\n', Lk);参数换算逻辑分三段：第一段先算额定电流和阻抗基准值，后面所有换算都基于这两个值；第二段由空载试验反推励磁支路，注意 Rc 和 Xm 是并联关系，不能直接用Zmag = Rc + jXm去解方程；第三段做短路试验换算，短路阻抗、短路电阻、短路电抗构成阻抗三角形，用勾股定理拆出漏抗。脚本输出里 Rk 和 Lk 是折算到原边的总值，如果后续搭建 Simulink 模型需要把绕组参数分开到原副边，一般按各占一半处理，副边折算值再按变比折算回去。
2.3 参数合理性校验表
反推出来的参数不要直接用，先拿经验区间做一次体检：
| 参数 | 经验范围（标幺值） | 越界可能原因 |
|---|---|---|
| 激磁电阻 Rc_pu | 50~500 | 空载损耗 P0 单位写错（W/kW 混淆） |
| 激磁电抗 Xm_pu | 20~200 | I0% 与 I0 实际值换算错误 |
| 短路电阻 Rk_pu | 0.005~0.02 | Pk 温度未折算到 75℃ |
| 短路电抗 Xk_pu | 0.04~0.15 | Uk% 取成了百分数分子但没除以 100 |
Rc_pu 和 Xm_pu 计算时记得各自除以 Zbase 再比对。这个环节最大的价值是挡掉大部分“仿真跑出来波形看着很怪”的问题——参数不对，后面 Simulink 里怎么调都白搭。
3. 稳态特性分析：电压调整率、效率曲线与负载工况仿真
3.1 为什么用 MATLAB 脚本而不是直接拉 Simulink
变压器稳态特性分析的核心是看外部特性（输出电压随负载变化的曲线）和效率特性（效率随负载系数变化的曲线）。这两条曲线用解析公式就能算，Simulink 的仿真步进反而带来不必要的数值开销。我通常的做法是先写 M 脚本把特性曲线全部算出来，画出图，确认趋势正确后再进 Simulink 做带暂态的验证。这样做的另一个好处是脚本本身就是活的文档：边界条件改一个数，整组曲线自动更新，比对着示波器截屏解释半天效率高得多。
电压调整率用简化公式ΔU% = β(Rk_pu*cosφ + Xk_pu*sinφ) * 100%就能算，其中 β 是负载系数（实际负载/额定负载），φ 是负载阻抗角。这个公式忽略了纵向分量的二阶项，工程上足够，但如果想精确可以按复数运算直接算。效率公式则是η = βSncosφ / (βSncosφ + P0 + β²Pk)，P0 是铁耗近似不变，β²Pk 是铜耗随负载平方变化，两条曲线的交点在 β = sqrt(P0/Pk) 处。
3.2 用 M 脚本生成负载特性与效率曲线
以下脚本基于第 2 节算出的参数，直接得到负载调整率和效率曲线：
%% 基于等效电路参数的稳态特性分析 % 参数沿用上节变量：Sn, P0, Pk, Rk_pu, Xk_pu beta = 0:0.05:1.5; % 负载系数 0~1.5 倍额定负载 cosphi = 0.8; % 负载功率因数 sinphi = sqrt(1 - cosphi^2); %% 电压调整率（感性负载） deltaU_pct = beta * (Rk_pu*cosphi + Xk_pu*sinphi) * 100; %% 效率曲线 % 负载输出功率 = β*Sn*cosφ，损耗 = P0 + β²*Pk Pout = beta * Sn * cosphi; % 输出有功功率 W Ploss = P0 + beta.^2 * Pk; % 总损耗 W eta = Pout ./ (Pout + Ploss) * 100; % 效率 % %% 最大效率点 beta_max_eta = sqrt(P0 / Pk); eta_max = (beta_max_eta*Sn*cosphi) / ... (beta_max_eta*Sn*cosphi + P0 + beta_max_eta^2*Pk) * 100; fprintf('最大效率点负载系数 β = %.2f，最大效率 = %.2f%%\n', ... beta_max_eta, eta_max); %% 绘图 figure('Color','w'); subplot(2,1,1); plot(beta, deltaU_pct, 'b-o', 'LineWidth', 1.5); grid on; box on; xlabel('负载系数 \beta'); ylabel('电压调整率 (%)'); title('变压器负载特性：电压调整率曲线'); % 标记额定工况点 hold on; plot(1, deltaU_pct(beta==1), 'ro', 'MarkerSize', 10, 'LineWidth', 2); legend('电压调整率', '额定点', 'Location', 'northwest'); subplot(2,1,2); plot(beta, eta, 'r-s', 'LineWidth', 1.5); grid on; box on; xlabel('负载系数 \beta'); ylabel('效率 (%)'); title('变压器效率特性曲线'); hold on; % 标记最大效率点 plot(beta_max_eta, eta_max, 'ko', 'MarkerSize', 10, 'LineWidth', 2); legend('效率', '最大效率点', 'Location', 'southeast');画图代码里用了subplot(2,1,1)把电压调整率和效率分上下两个子图展示，这样可以直接观察过载区的特性变化。函数fprintf把最大效率点输出在命令行。注意beta.^2用了点乘，是因为beta是向量，不加点乘 MATLAB 会尝试做矩阵乘法报错或算出完全错误的结果。legend里的Location参数建议按曲线实际位置调整，如果曲线在左上角相交，就改成'southwest'。
3.3 读取数据系列导入工作区的技巧
实际项目中铭牌数据经常会以表格形式提供（CSV、Excel 或 .docx 里的表格），手动敲进脚本既慢又容易抄错。MATLAB 里我一般用readtable直接读：
% 读取铭牌数据表（假设已有 transformer_data.xlsx） data = readtable('transformer_data.xlsx'); Sn = data.Sn; % 单位 VA P0 = data.P0; % 单位 W Pk = data.Pk; % 单位 W这种做法的好处是脚本与数据分离，多台变压器只需要在同一张表里加行，脚本重跑即可横向对比。如果表头是中文或者含特殊字符，readtable会自动生成变量名，取字段时拿data.(genvarname('空载损耗(Ｗ)'))这种写法即可，不过建议做表的时候直接用英文变量名，省得后面引用时写一长串转义。数据文件放工作目录下，用pwd确认当前路径正确，readtable报文件不存在时八成是路径问题，用trace或者直接disp(pwd)排查。
4. Simulink 与 Simscape 物理建模：从稳态走向暂态
4.1 Simscape 变压器模型和 SimPowerSystems 旧模型的选型差异
MATLAB 里搭变压器模型有两条路：一是 Simscape Electrical 的 Transformer 模块（含线性与饱和两种）、二是用受控电压源/电流源手搭理想变压器加外接漏抗励磁支路。Simscape 物理模型的好处是端口方向自动处理，可以直接用 Three-Phase Transformer 做三相连接，坏处是模型封装了内部细节，出问题时不好查是哪一段磁路饱和。受控源手搭的好处是每一条支路都暴露在外面，很适合教学演示和学生理解结构，坏处是端口正负方向和参考地稍不留神就接反。
我一般这么定夺：三相工频变压器用 Simscape 的 Three-Phase Transformer (Two Windings) 省时间，单相变压器研究励磁涌流或 DC 偏磁时用手搭模型，因为饱和特性的自定义更直接。MATLAB 新版本里 SimPowerSystems 已经并入 Simscape Electrical，模块库搜Transformer出来的结果和旧版差别很多，如果搜不到先确认是否在 Simscape 而不是 Simulink 库。
4.2 Simulink 最小可跑模型搭建步骤
以下步骤基于 MATLAB R2020a 以后版本界面，操作位置差异不会太大：
- 在 MATLAB 命令窗口输入
simulink打开 Simulink 起始页，新建 Blank Model。 - 打开库浏览器（Library Browser），在搜索框输入
Simscape > Electrical > Specialized Power Systems > Fundamental Blocks > Elements，拖入Two-Winding Transformer（单相双绕组变压器）或Three-Phase Transformer (Two Windings)（三相双绕组变压器）。 - 双击变压器模块，把第 2 节算出的参数填入 Units 选择
SI，勾选Saturation可打开饱和特性，不勾就是线性模型。 - 原边接
AC Voltage Source，副边接Series RLC Load，用Powergui模块设置仿真类型为Continuous。 - 用
Voltage Measurement和Current Measurement把电压电流信号引出到Scope，用于观察波形。
参数填的时候最容易错的是并联支路（magnetizing branch）和串联支路（series branch）的物理量分布：Simscape 变压器的对话框里一般有Winding 1和Winding 2两个绕组各自的电阻漏抗输入框，外加Magnetizing resistance Rm和Magnetizing inductance Lm，整个等效电路和第 2 节算出的 Rc、Lm、Rk、Lk 完全对应。注意规范的是，Simscape 里Rk要拆成两半分别填进原副边的绕组电阻框，如果直接整值填入某一侧绕组，短路试验仿真出来的阻抗会偏大一倍。
4.3 空载合闸和绕组短路的暂态特性分析
稳态特性仿真跑通后，真正的技术点落在暂态分析。空载合闸是变压器最典型的暂态过程，磁通从 0 冲到约 2 倍稳态磁通峰值，在剩磁影响下可能冲到 2.5 倍甚至更高，导致励磁涌流高达额定电流的 6~8 倍。Simulink 里模拟这个场景只需要把AC Voltage Source的Initial Phase改成 0 度（电压过零点合闸是最严酷工况），把仿真停止时间设为 0.2 秒，用Scope看电流波形，可以看到明显偏置的尖顶波形。
短路暂态则把负载换成Three-Phase Fault，设置故障时间从 0.05 秒开始，持续 0.1 秒，观察副边短路电流。此时电流往往达到额定电流的 4~10 倍（取决于 Uk% 大小），如果模型里没有饱和特性，短路电流峰值会按线性模型线性放大，实际饱和会让电流峰值略低。但我一般会把饱和特性勾上，“宁可慢一点也要把磁路饱和的影响算进去”，当需要提交报告时这个对比图很有说服力。用powergui里的Steady-State工具可以快速确认稳态工作点，避免启动阶段的超调影响对稳态值判断。
5. 仿真与实测误差的偏差溯源校核方法
5.1 把两组模型结果叠在同一张图上验证
同一台变压器，M 脚本算稳态，Simulink 算暂态，两套结果的衔接点应该落在负载特性曲线上。我通常在 Simulink 模型里让变压器带额定负载运行，等暂态消退后（大约 0.1 秒以后），用mean函数取电压电流平均值，然后和脚本算出的电压调整率对比，误差超过 5% 就要查参数量纲或接线。
%% 从 Simulink 输出 out 结构取稳态数据 V2_measured = mean(out.V2.Data(end-1000:end)); % 取末端稳态段电压 I2_measured = mean(out.I2.Data(end-1000:end)); % 取末端稳态段电流 V2_sim = U2n * (1 - deltaU_pct(21)/100); % 脚本预测值（β=1.0 处） err = abs(V2_measured - V2_sim) / V2_sim * 100; fprintf('Simulink 与 M 脚本预测电压偏差：%.2f%%\n', err);out.V2.Data的长度取决于仿真步数和保存设置，直接取最后 1000 个点做平均可以滤掉尾部毛刺，拿到稳定的数值。如果波形还没完全进入稳态就把数据截断，均值会偏高或偏低，这种情况可以把仿真时间加长。在 Simulink 配置参数里用Data Import/Export勾选输出tout、yout也能达到同样效果，不过out结构体方式在新版本中更统一。
5.2 变压器仿真没法复现时的分层排查顺序
如果仿真波形和预期完全对不上，按下面顺序逐级排查：
- 第一步：看 Powergui 是不是
Continuous，Phasor 模式只适合工频稳态，看不了励磁涌流的暂态尖峰。 - 第二步：看变压器模块的
Measurement端口有没有接阻抗，Simscape 里测量端口是无穷大内阻的电压表，不接东西会导致虚拟谣言，仿真直接报错变形。 - 第三步：看饱和特性的
Saturation characteristic表（电流-磁链折线点），默认参数是适应小变压器的，如果只改容量不改这条折线，大变压器跑到高磁密区时波形会出现异常断裂。 - 第四步：看仿真步长，默认
ode23tb可以跑，但遇到极短暂态过程建议改成ode15s并设最大步长 1e-4 秒，否则漏掉电流尖峰。
表格形式给出典型异常和处置方向：
| 故障现象 | 可能原因 | 处置措施 |
|---|---|---|
| 空载电流波形完全对称无尖峰 | 饱和特性未启用 | 勾选 Saturation 并填入磁链-电流折线点 |
| 稳态电压比铭牌低 10% 以上 | 漏抗参数填错位置 | 检查原副边绕组的漏感分配，原副边各填一半 |
| 励磁涌流峰值明显偏小 | 合闸角恰好 90° | 把电压源初始相位改为 0° 重测最恶劣合闸条件 |
| 仿真速度极慢且报数值奇异 | 励磁电感量级错误 | 检查 Lm 是否按 H 为单位填入，mH 直接填会导致过大的 dψ/dt |
这套排查流程能覆盖我在工程中遇到的大部分问题。实际仿真时如果发现 Simulink 结果和脚本值有系统性偏差，先不要急着调参，重新确认两套模型采用的是同一组基准值（容量、频率、电压基准）。很多团队里脚本用线电压、Simulink 模块里填相电压，一个sqrt(3)的因子差就会让所有结果对不上。写文档时把基准值表格贴在模型说明里，比事后解释高效得多。
%% 最后：保存工作区参数供下次仿真直接加载 save('transformer_params.mat', 'Sn', 'U1n', 'U2n', 'f', ... 'P0', 'I0_pct', 'Pk', 'Uk_pct', 'Rc', 'Lm', 'Rk', 'Lk');在这个脚本里，save命令把铭牌参数和换算结果全部存在 mat 文件里，下次开新工程直接load('transformer_params.mat')就能恢复整个工作区，这在归档模型版本或做多台变压器对比时非常实用。经常做变压器仿真的话，建议把这段参数换算和特性分析封装成函数文件transformer_analysis.m，入参是铭牌数据，出参是结构体，后续所有脚本和 Simulink 模型都从这个函数取数，避免每个脚本里维护一套变量名，减少复制粘贴造成的低级错误。
本文还有配套的精品资源，点击获取
