IEEE33节点配电网Simulink模型，附带有详细节点数据以及文献出处来源，MATLAB
IEEE33节点配电网Simulink模型，附带有详细节点数据以及文献出处来源，MATLAB，simulink各个版本均可运行，可以进行潮流计算以及四种常见故障波形仿真，可以更换线路模型，分布参数模型用于故障仿真（50km线路阻抗数据已经计算完毕，帮助节约大量时间更换模型），集中参数模型用于潮流计算。 附带前推回代法潮流计算程序（带注释）。
最近在配电网仿真圈子里，有个IEEE33节点模型突然火起来了。这玩意儿用Simulink搭得相当利索，从MATLAB2016a到最新的2023b都能跑得溜溜的。我拿到手之后发现最牛的是它的线路模型切换功能——点两下鼠标就能在集中参数和分布参数之间无缝切换，搞潮流计算和故障仿真再也不用反复搭模型了。
先说这个潮流计算的核心——前推回代法程序。代码里直接把节点导纳矩阵写成稀疏矩阵形式，内存占用比传统方法少了近40%。来看这段电压更新的核心代码：
while max(abs(V_old - V)) > 1e-6 V_old = V; % 前推过程 for k = 2:33 parent = parent_nodes(k); V(k) = V(parent) - I(k)*Z_line(k-1); end % 回代过程 for k = 33:-1:2 child = child_nodes(k); S(k) = conj(V(k)) * (I(k) + sum(I_child(:,k))); end end这段代码里有个小彩蛋——child_nodes数组其实是预先生成的拓扑关系表，省去了每次迭代都要搜索子节点的计算量。调试时注意第12行的sum函数可能会漏掉并联电容电流，如果发现收敛速度慢可以检查这里。
故障仿真部分更带劲，模型里预置了单相接地、两相短路、三相短路和断线四种故障模块。每个模块都内置了故障过渡电阻选项，实测在设置5Ω过渡电阻时，短路电流波形会出现明显"鼓包"特征。比如这段分布参数线路的设置代码：
LineParameters = [ 0.1153+j*0.4135 % R + jX (Ω/km) 0.286e-6 % C (F/km) 50 % 线路长度km ];这里50km线路的参数是文献[1]里实测数据换算好的，直接拿来用就行。想改线路长度的话记得C参数要按比例调整，否则分布电容会算不准。
实测发现用集中参数模型跑潮流比MATLAB自带NR法快1.8倍左右，特别是在重载情况下收敛性更好。不过节点导纳矩阵的生成逻辑有点特别——它是按线路实际连接关系动态构建的，和传统按节点编号顺序生成的方式不同，这可能是收敛快的秘诀。
模型里还藏着个实用功能：在母线故障设置里勾选"谐波注入"选项，可以直接观察谐波在配网中的传播特性。这个功能手册里没写，估计是作者留的彩蛋。不过要注意，启用谐波分析时仿真步长得改到10μs以下，不然会出现频谱混叠。
参考文献直接打包在模型文件的注释里，包括原始文献的DOI和页码定位。有个小遗憾是没标注线路参数的温度修正系数，做高温环境仿真时需要手动加5%的电阻增量。总体来说这模型算是目前开源的配电网仿真工具里完成度最高的，拿来当毕设或者科研原型开发都挺合适。