别再复制粘贴了！手把手教你从零编写MatPower的case文件（以6节点电网为例）
从电网拓扑图到MatPower案例：6节点系统建模实战指南
电力系统仿真研究中，MatPower作为一款开源工具包，其case文件的编写质量直接影响潮流计算结果的可靠性。许多初学者习惯直接修改现成案例，却对底层数据结构一知半解。本文将带您从一张6节点电网拓扑图出发，像搭积木一样构建完整的case文件，重点解析那些容易出错的参数设置。
1. 准备工作：理解电网结构与MatPower数据模型
拿到一张6节点电网拓扑图时，首先需要明确几个核心要素：
- 节点属性：哪些是负荷节点(PQ)、发电机节点(PV)和平衡节点(Slack)
- 支路参数：线路阻抗(r+jx)和变压器变比
- 设备容量：发电机出力范围与负荷大小
以一个典型6节点系统为例，其拓扑特征如下：
节点1(Slack)───┬───节点2(PV) │ 节点3(PQ)───节点4(PQ) │ 节点5(PV)───节点6(PQ)提示：实际建模前建议用铅笔标注每个节点的类型和关键参数，避免后续混淆
MatPower通过五个核心矩阵描述电网：
| 矩阵名称 | 功能描述 | 必填字段示例 |
|---|---|---|
| bus | 节点参数 | type, Pd, Vm |
| gen | 发电机参数 | Pg, Qmax, Vg |
| branch | 支路参数 | r, x, rateA |
| gencost | 发电成本 | (可选) |
| baseMVA | 基准容量 | 通常设为100 |
2. 构建bus矩阵：电网节点的数字化表达
bus矩阵是电网的骨架，每行对应一个物理节点。以下6节点系统的典型配置：
% 列说明：[bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin] bus = [ 1 3 0 0 0 0 1 1.05 0 230 1 1.05 0.95; 2 2 0 0 0 0 1 1.05 0 230 1 1.05 0.95; 3 1 90 30 0 0 1 1.0 0 230 1 1.05 0.95; 4 1 85 20 0 0 1 1.0 0 230 1 1.05 0.95; 5 2 0 0 0 0 1 1.05 0 230 1 1.05 0.95; 6 1 70 25 0 0 1 1.0 0 230 1 1.05 0.95; ];关键参数避坑指南：
type字段：
1：PQ节点（负荷）2：PV节点（发电机）3：平衡节点（参考节点）
标幺值转换：
- Pd/Qd使用有名值(MW/MVar)
- Vm/Va使用标幺值
常见错误：
- 混淆PQ/PV节点类型
- 忘记设置平衡节点
- 电压幅值超出合理范围(0.95-1.05 p.u.)
3. 配置gen矩阵：发电机参数详解
gen矩阵描述所有可控电源，注意与bus矩阵的对应关系：
% 列说明：[bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin] gen = [ 1 0 0 100 -100 1.05 100 1 250 10; % Slack发电机 2 80 0 80 -80 1.05 100 1 100 20; 5 70 0 60 -60 1.05 100 1 80 15; ];特殊参数说明：
- mBase：建议保持100与baseMVA一致
- Vg：PV节点电压控制目标值
- Pmax/Pmin：根据发电机实际容量设置
注意：Slack发电机虽然Pg设为0，但必须设置足够的调节容量(Pmax/Pmin)
4. 定义branch矩阵：连接节点的电网血脉
支路参数直接影响潮流分布，需要特别注意标幺值计算：
% 列说明：[fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax] branch = [ 1 2 0.02 0.06 0.030 250 250 250 0 0 1 -360 360; 1 4 0.08 0.24 0.025 250 250 250 0 0 1 -360 360; 2 3 0.06 0.18 0.020 250 250 250 0 0 1 -360 360; 3 4 0.01 0.04 0.015 250 250 250 0 0 1 -360 360; 4 5 0.08 0.24 0.025 250 250 250 0 0 1 -360 360; 5 6 0.04 0.12 0.010 250 250 250 0 0 1 -360 360; ];阻抗计算要点：
标幺值公式： [ Z_{pu} = \frac{Z_{actual}}{baseZ} = \frac{Z_{actual}}{(baseKV)^2/baseMVA} ]
典型错误：
- 直接使用欧姆值未转换
- 电阻/电抗比例不合理（正常x/r≈3-10）
5. 完整案例与验证技巧
整合所有矩阵的完整case文件结构：
function mpc = case6 mpc.version = '2'; mpc.baseMVA = 100; mpc.bus = [...]; % 填入前述bus矩阵 mpc.gen = [...]; % 填入前述gen矩阵 mpc.branch = [...]; % 填入前述branch矩阵验证模型正确性的三个步骤：
语法检查：
runpf('case6') % 基础潮流计算结果合理性验证：
- 平衡节点出力是否在合理范围
- PV节点电压是否维持设定值
- 线路负载率是否超限
可视化校验：
plot(mpc) % 需要安装MatPower的绘图工具
在首次调试时，我习惯先简化模型（如设所有负荷为0），逐步增加复杂度。曾经有个项目因为节点类型设置错误，导致迭代50次仍不收敛，最后发现是一个PV节点被误设为PQ节点。