简介：一份面向电气类专业课程设计/毕业设计场景的Matlab变压器仿真建模与特性分析完整设计报告。报告以Matlab电力系统模块为平台，给出三相双绕组变压器（500kV/230kV、450MVA）仿真模型搭建思路，涵盖变压器饱和特性、磁滞与剩磁影响，以及空载、负载、空载合闸、副边突然短路等典型工况的仿真过程与结果分析，适合需要掌握Simulink变压器建模方法或撰写电气方向设计文档的学习者参考。资源为单个docx文档，共1个文件，压缩包仅321KB，轻量便于下载，文档内包含设计任务书、学生日志、摘要、方案论证、系统设计与图表等完整结构，可直接作为格式与内容参照。目前已有154人学习，对需要快速完成变压器仿真方向课程报告或设计报告的同学有较高参考价值。
1. 为什么电气专业设计报告里的变压器特性分析，第一版验证要放在 matlab 仿真上
电气专业的变压器设计报告，最难交差的不是计算书，而是“算出来的特性值对不对”。电压调整率算到 4.7%，空载电流按 2.5% 取，这些数只能靠试验验证；但实验室不一定有对应容量的三相变压器，更不可能让你短路顶额定电流。matlab 的 Simulink 恰好把这条验证链补上：把设计参数标幺化之后填进变压器仿真模型，空载、短路、负载三种试验在同一张模型图上切换，特性分析数据直接落回报告表格。
这篇按课程设计最常见的路径讲：模型怎么选、参数怎么填、试验怎么做、批量扫描怎么跑，最后把效率曲线整理成能放进设计报告的高清图。适合正在写变压器设计报告，或者手里有一堆仿真波形却不知道下一步怎么处理和印证的人。整条链路不依赖真实设备，也不需要额外的硬件在环环境。
2. 变压器仿真建模：Simscape 变压器模块怎么选、参数怎么填
2.1 三种建模路径，设计报告里选哪种更稳
在 matlab 环境里做变压器仿真建模，常见做法有三条路，差别很大。
第一种是直接使用 Simscape Electrical（旧称 SimPowerSystems）里的变压器模块，比如 Three-Phase Transformer (Two Windings)。这类模块内部已经处理好了标幺值换算、磁路饱和、铁损电阻和漏感分相，填入铭牌参数就能跑，是绝大多数设计报告的主路径。
第二种是用 Simulink 基础模块搭微分方程，按 u₁ = R₁i₁ + L₁ di₁/dt + M di₂’/dt 的关系逐相建立。这条路适合解释气隙、匝数、漏感之间的解析关系，但饱和特性要自己查表插值，三相耦合绕组的处理非常容易出错。报告里用来推导公式可以，作为特性分析依据就不够可靠。
第三种是写成状态方程之后封装成 Level-2 MATLAB S-Function，灵活性最高，但调试成本也高，一般只在研究磁滞模型、损耗模型这类课题里用。
设计报告选第一种最合理。注意一点：Simscape Electrical 在 matlab 中作为独立工具箱存在，安装后在命令行敲ver('SimscapeElectrical')能查到版本和许可证信息再继续，否则建模型到一半会报模块找不到。另外，如果你做的是反激、推挽这类开关电源变压器，Simscape 的工频变压器模型并不适用，磁性元件要自己按磁芯参数建非线性模型，那是另一个方向。
2.2 用 Two Windings 模型搭一个最小可运行试验台
以一台 100kVA、10kV/0.4kV、Dyn11 连接组别的配电变压器为例，把最小试验台搭起来。这个模型同时用于后面的空载、短路和负载试验，一次建好，反复切换。
搭建步骤：
- MATLAB 命令行运行
simulink，在库浏览器里定位 Simscape > Electrical > Specialized Power Systems > Fundamental Blocks > Elements。 - 拖入 Three-Phase Transformer (Two Windings)、Three-Phase Source、Three-Phase RLC Load 和两个 Three-Phase V-I Measurement。
- 在电源和变压器之间放一个断路器（Three-Phase Breaker），负载侧再放一个，用来切换试验工况。
- 从 Specialized Power Systems 里补一个 Powergui 模块到模型顶层，这一步不能省，Simscape Electrical 模型没有 Powergui 会直接拒绝仿真。
- 把求解器设为 ode23tb，停时先给 0.2s。
变压器模块的对话框参数是标幺值体系，直接填设计值即可。以这台变压器为参考，参数表如下。
| 参数项 | 填入值 | 对应含义 |
|---|---|---|
| Nominal power and frequency | 100e3 VA, 50 Hz | 额定容量与电网频率 |
| Winding 1 (V1 Ph-Ph, R1, L1) | 10e3 V, 0.002 pu, 0.04 pu | 一次侧线电压、短路电阻与漏感标幺值 |
| Winding 2 (V2 Ph-Ph, R2, L2) | 400 V, 0.002 pu, 0.04 pu | 二次侧线电压、归算到本侧的短路参数 |
| Magnetization resistance Rm | 200 pu | 对应额定电压下的铁耗，量级由空载损耗决定 |
| Magnetization inductance Lm | 50 pu | 对应空载电流约 2% 额定电流的励磁需求 |
| Saturation characteristic | [0 0.02 0.6 1.3; 0 1.0 1.1 1.2] | 第一行电流 pu，第二行磁通 pu，膝点放在 1.2 pu |
参数填好后，用一段最简脚本做冒烟测试，确认模型能正常跑。
% transformer_min_run.m % 最小验证脚本：加载模型、设置停时、运行并取出二次侧电压 model = 'tf_oc_test'; load_system(model); set_param(model, 'StopTime', '0.2'); set_param(model, 'Solver', 'ode23tb'); simOut = sim(model, 'ReturnWorkspaceOutputs', 'on'); % 从日志里取二次侧电压瞬时值和对应时间轴 v2 = simOut.logsout.getElement('V2').Values.Data; t = simOut.logsout.getElement('V2').Values.Time; % 只取 0.1s 之后的稳态段计算 RMS V2_rms = sqrt(mean(v2(t > 0.1).^2)); fprintf('二次侧电压 RMS = %.1f V\n', V2_rms);脚本的前三段是模型路径和仿真配置，核心在最后三行。用t > 0.1做逻辑索引，是为了躲开合闸瞬间的励磁涌流和直流偏移分量，变压器刚接上电网的前几个周波波形根本不是稳态正弦，直接对全波形做 RMS 会把电压算高。这里用sqrt(mean(v2(mask).^2))而不是rms函数，是为了不依赖 Signal Processing Toolbox，纯 matlab 环境也能运行。
2.3 空载电流偏大偏小，先查励磁支路参数而不是变比
许多仿真模型第一次跑空载试验，二次侧电压正常，但空载电流和设计值对不上。这不是变比算错，而是励磁支路参数没有和目标特性挂钩。
Simscape 变压器模块里有两个励磁参数：Rm 对应铁耗，Lm 对应空载电流。两者是并联关系，量级上没有必然联系。很多模板把 Rm 和 Lm 都留成默认的 500 pu，跑出来空载电流只有 0.2% 额定电流，空载损耗也只有额定容量的 0.2%。对小容量变压器来说，这个数字明显偏低，报告里没法交代。
估算原则很简单：空载电流 I₀ 取额定电流的 2%~4%，那么励磁感抗标幺值 Xm ≈ 1 / I₀_pu，也就是 25~50 pu 左右；空载损耗 P₀ 约占额定容量的 0.5%~1.5%，铁损电阻标幺值 Rm ≈ 1 / P₀_pu，大约 70~200 pu。两个参数先按这个量级填，跑完空载试验后再用实测仿真波形反推修正。
另一个常见问题是饱和特性曲线。模块里的 Saturation characteristic 用两行向量描述磁化曲线，默认的线性区只能走到 1.2 pu。如果额定工作点设置得离膝点太近，变压器在负载投入或电压抬升时会出现异常的谐波电流，波形失真明显。膝点建议放在 1.2 pu~1.3 pu，把额定工作点留在线性区内。励磁涌流的本质就是阶跃电压作用在非线性励磁支路上的磁链响应，膝点位置直接决定涌流峰值和衰减时间，这一步不调好，后面短路试验和负载试验的波形都会带着莫名的畸变。
3. 变压器特性分析怎么做：空载、短路与负载试验的仿真步骤
3.1 空载试验：一次侧加额定电压，二次侧开路
空载试验的目标是测出励磁参数、变比和空载损耗。在 Simulink 里操作很简单：断开二次侧断路器，让负载侧完全开路，一次侧加额定线电压 10kV，仿真时长 0.2s。
电压源直接给额定值就一定对吗？要确认源的内阻设置。Three-Phase Source 里如果不勾选 Specify internal resistance, 默认是理想电压源，对空载试验没问题；但如果是把某个带内阻的发电机模型当作电源，空载试验里源阻抗和励磁阻抗会分压，10kV 的源电压到变压器端子上只剩 9.2kV 左右，空载电流会系统性偏小。设计报告里最好在变压器一次侧加一个 Three-Phase V-I Measurement 直接量端子电压，不要拿源电压当变压器电压。
数据处理时关注三个量：一次侧电流 I₀、三相总功率 P₀、二次侧电压 U₂₀。P₀ 在低压配电变压器里约占额定容量的 0.5%~1.5%，I₀ 约占额定电流的 2%~4%，U₂₀ 应该在 400V 上下，偏差超过 1% 就说明变比或连接组别设置有问题。
空载损耗基本等于铁耗，因为此时一次侧铜耗 I₀²R₁ 只有几十瓦量级，相对于千瓦级的铁耗可以忽略。把 P₀ 和 I₀ 代回公式，就能得到 Rm 和 Lm 的仿真辨识值，和 2.2 节里的估算量级对照。波形上注意看 I₀ 的包络：如果合闸后电流峰值达到额定电流的 5~8 倍且衰减很慢，说明膝点设置偏低，励磁支路饱和过深。
3.2 短路试验：二次侧短接，一次侧升到短路电压
短路试验在仿真里比真实试验简单，但数据陷阱更多。先把二次侧断路器合上，并把负载换成 1e-6 欧姆的三相电阻模拟短路，然后调节一次侧电压源的幅值，让一次侧电流刚好等于额定电流，记录此时的电压 U₁k 和三相总功率 Pₖ。
这里有一个真实试验里能避免、仿真里反而容易犯的错：直接拿 10kV 额定电压去顶短路模型。仿真不会冒烟，但电流会飙到额定电流的几十倍，变成了一场短路故障仿真，而不是特性分析。正确做法是先估算短路阻抗标幺值蕶加上铜耗损耗给出的短路电压百分比 um，依此电压调试源幅值。默认 R=0.002 pu、X=0.04 pu 时，短路电压大约 4%，所以一次侧电压调到 400V 附近，电流就到额定值了。
短路试验的测量数据用来分离铜耗和漏抗。三相总功率 Pₖ 就是额定电流下的铜耗，有功分量 Rk 由 Pₖ = 3·I₁n²·Rk 得到，漏抗 Xk 由短路阻抗模和电阻分量勾股相减得到。断开磁路饱和的影响：额定短路电流下铁芯磁通远低于额定值，励磁支路基本不参与分流，所以短路参数直接代表一次侧和二次侧归算电阻、漏感的串联值。
仿真时长给 0.1s 足够了，短路电流经过最初一两个周波的直流偏移之后迅速进入稳态。取 RMS 时同样只取后半段波形，不然残存的直流分量会让电流有效值算大。
3.3 负载特性与电压调整率的仿真读数
负载特性是变压器特性分析里和设计报告结论直接相关的一节。做法是不改变一次侧额定电压，把二次侧负载分别设成纯阻性、阻感性、阻容性，在额定容量的一定比例下记录 U₂ 和 I₂，画出外特性曲线，并计算电压调整率。
电压调整率定义为 ΔU% = (U₂₀ − U₂) / U₂₀ × 100%。负载电流、负载功率因数一致时，感性负载的电压调整率最大，阻性次之，容性负载可能出现负调整率，也就是负载越重电压反而越高。这是因为容性电流对主磁通起助磁作用，在仿真波形上表现为 U₂ 随电流上升的“上翘”。
在一组典型参数下扫描得到的数量级如下表，用于核对你的模型是否合理。
| 负载类型 | 功率因数 | 负载系数 1.0 时 ΔU% | 负载系数 1.2 时 ΔU% |
|---|---|---|---|
| 纯阻性 | 1.0 | 3.5 | 4.2 |
| 阻感性 | 0.8 滞后 | 5.8 | 7.1 |
| 阻容性 | 0.8 超前 | -1.6 | -2.0 |
这几个值的绝对大小由短路阻抗标幺值决定，短路阻抗越大，电压调整率越大。如果你的模型跑出来阻性负载电压调整率超过 8%，回去查 R₁、L₁ 的标幺值是不是填得太大了。
% extract_char.m % 从仿真输出中提取二次侧线电压、电流有效值 function [U2, I2] = extract_char(simOut, t_start) v2 = simOut.logsout.getElement('V2').Values.Data; i2 = simOut.logsout.getElement('I2').Values.Data; t = simOut.logsout.getElement('V2').Values.Time; mask = t > t_start; U2 = sqrt(mean(v2(mask).^2)); I2 = sqrt(mean(i2(mask).^2)); end这段代码把电压电流的波形数据压成两个有效值，后面批量扫负载时每次调用一次即可。电流 I₂ 默认是相电流还是线电流，取决于接线的是三相 V-I Measurement 的哪个输出端口；如果之后要把仿真值和铭牌电流对比，线电流和相电流之间的 √3 关系要先明确。
4. 用 m 脚本批量扫负载系数，把特性分析变成可复现的自动化流程
4.1 用 set_param 改负载，用 sim 批量跑
手动改负载再点运行，做 12 个工况能花掉一下午，而且每轮都在手动搬数据。电气专业设计报告要的是可复现的曲线，不是几张孤立波形。把参数扫描写成 m 脚本是这条链路的正解。
这里用 Three-Phase RLC Load 的 Active Power 参数来控制负载轻重，比直接改阻抗值更贴近变压器铭牌习惯。每轮循环只改一个参数，然后调用sim获取输出。
% batch_load_scan.m % 扫负载系数 0.1~1.2，完成阻性负载特性试验 model = 'tf_load_test'; load_system(model); Pn = 100e3; % 额定容量 beta_list = []; % 负载系数 U2_list = []; % 二次侧线电压 RMS I2_list = []; % 二次侧电流 RMS P2_list = []; % 三相总输出有功 for beta = 0.1:0.1:1.2 % 负载有功按额定容量的 beta 倍设置，保持功率因数为 1 set_param([model '/Load'], 'ActivePower', num2str(beta * Pn)); simOut = sim(model, 'StopTime', '0.2'); logsout = simOut.logsout; % 波形数据在 logsout 里，对应模块输出端右键 Log Selected Signals t = logsout.getElement('V2').Values.Time; v2 = logsout.getElement('V2').Values.Data; i2 = logsout.getElement('I2').Values.Data; p2 = logsout.getElement('P2').Values.Data; mask = t > 0.1; % 跳过合闸暂态 beta_list(end+1) = beta; %#ok<SAGROW> U2_list(end+1) = sqrt(mean(v2(mask).^2)); I2_list(end+1) = sqrt(mean(i2(mask).^2)); P2_list(end+1) = mean(p2(mask)); end这个脚本有三个关键点。第一，set_param改参数必须发生在每次sim调用之前，而且参数名要用模块对话框里看到的字段名，写错不会报错而是静默不生效，这是 set_param 机制最坑的地方。第二，Simulink 的输出日志不会自动记录所有信号，必须在模块输出端子上右键选择 Log Selected Signals，或者加 To Workspace 模块指定变量名，否则logsout.getElement会空指针。第三，sim函数在 MATLAB R2010a 之后默认以 SimulationOutput 对象返回数据，ReturnWorkspaceOutputs参数控制是否同时写回基础工作区，保留默认的'on'可以让调试时随时查看中间变量。
4.2 循环里清理数据，而不是循环后补处理
批量扫描最容易出的问题不是脚本报错，而是数据里混着暂态。每一轮仿真都是合闸启动，0.1s 之后才算进入稳态。上面脚本用t > 0.1统一截断，对大多数配电变压器模型够用；但如果饱和膝点调得很低，励磁涌流可能持续 10 个周波以上，0.2s 停时就不够了。判断标准以一次侧电流包络为准：波形包络不再起伏之后，再留 2~3 个周波作为采样窗口。
这一步如果做得不干净，后面画的效率曲线和电压调整率曲线在轻载段会带毛刺，因为轻载时暂态电流占比相对更大，截断误差更明显。宁可把停时调长到 0.3s，也不要为了省时间在数据里烧杂质。
4.3 批量仿真太慢，切加速模式而不是换电脑
12 组工况跑下来，如果每组都要 5 秒以上，检查两个地方：仿真模式和解算器。默认 Normal 模式下，Simulink 每轮都重新编译模型；同一模型重复仿真时可以在第一轮之后切换加速模式。
set_param(model, 'SimulationMode', 'accelerator'); % 或直接传参数 simOut = sim(model, 'StopTime', '0.2', 'SimulationMode', 'rapid');accelerator会为首轮生成仿真目标代码，后续轮次直接执行编译后的代码，对非线性变压器模型通常能快 3~10 倍。rapid模式更快，但要求模型完全支持快速加速目标，且每次sim调用都要重新把参数传给外部模式，在某些模型里反而慢。我的建议是先用accelerator，速度不够再评估rapid。注意两种模式下数值结果和 Normal 模式会有微小差异，设计报告里如果要对多个工况做横向对比，所有工况必须用同一种仿真模式，混用会让曲线出现本不该有的阶跃。
解算器方面，ode23tb是专为刚性电路问题设计的，变压器这种带储能元件和快速开关动作的模型用它最稳。ode45在空载试验里能跑，短路试验里经常因为步长过小卡住。
5. 效率曲线与最大效率点校验：把仿真数据整理成报告图表
5.1 取 P1 和 P2，效率计算的输入别用错功率
效率特性分析需要一次侧输入有功 P₁ 和二次侧输出有功 P₂，两个量要同时来自仿真，不能用 P₂ 除以 P₂ 加损耗来凑。在模型里给一次侧、二次侧各加一个 Active & Reactive Power 模块（位于 Simscape > Electrical > Specialized Power Systems > Fundamental Blocks > Measurements），输出端接 To Workspace，批量扫描时就会自动记录每个工况下三相总有功。
eta = P2_list ./ P1_list;注意 P1 是三相总输入有功，不是单相功率，Active & Reactive Power 模块输出已经是三相合计功率。如果自己用瞬时电压电流相乘再求均值，记得要写成mean(v_a.*i_a + v_b.*i_b + v_c.*i_c)，漏掉相序对应关系是这类手写功率计算最常用的出错点。
5.2 最大效率点：仿真峰值和理论公式对照
变压器效率曲线有一个反直觉的结论：最大效率点不在额定工况，而在负载系数等于铁耗与铜耗之比平方根的位置。理论公式 β_opt = sqrt(P₀ / Pₖ)，P₀ 来自空载试验，Pₖ 来自短路试验，这两个值在前面已经拿到，直接代入。
P0 = 1200; % 空载损耗，来自空载试验 Pk = 2000; % 短路损耗，来自短路试验 beta_opt = sqrt(P0 / Pk); eta = P2_list ./ (P1_list); % 直接用测量数据，不套公式 [eta_max, idx] = max(eta); fprintf('仿真最大效率 %.2f%% @ beta=%.2f\n', eta_max*100, beta_list(idx)); fprintf('理论最优负载系数 beta_opt = %.2f\n', beta_opt);把仿真扫到的峰值点和理论公式算出的 β_opt 对照，偏差在 0.05 以内说明损耗模型设置合理。偏差过大时，先查空载损耗 P₀ 是不是被 Rm 的标幺值带偏了，再查短路损耗 Pₖ 是不是在电流没到额定值时就采了数。
5.3 曲线图导出：300dpi 起步，尺寸和字体先定好
设计报告里的效率曲线一般插在“特性分析”章节，图片清晰度直接影响评审观感。R2020a 之后的 matlab 用exportgraphics导图，之前的老代码用print，导出效果差别很大。
figure('Units', 'centimeters', 'Position', [2 2 14 9]); plot(beta_list, eta * 100, 'o-', 'LineWidth', 1.2); xlabel('负载系数 \beta'); ylabel('效率 / %'); grid on; title('变压器效率特性曲线'); exportgraphics(gcf, 'efficiency_curve.png', 'Resolution', 300);exportgraphics按打印分辨率渲染，矢量线型和文字在 300dpi 下不会出现边缘锯齿，直接插进 Word 或者 LaTeX 都不需要二次处理。14cm×9cm 的画布尺寸适配文档半栏排版，字号用 matlab 默认的 10pt 就够，放大到全页宽也不会糊。
本文还有配套的精品资源，点击获取
