简介这份资源面向电力系统分析的学习者与研究人员提供基于MATLAB实现的牛顿-拉夫森潮流计算程序覆盖6节点与9节点IEEE标准测试系统可用于理解稳态下电压、功率分布与线路潮流的求解过程适合课程设计、算法验证与科研入门等场景。压缩包共2个文件均为.m脚本整体约2KB分别对应6节点与9节点系统的潮流计算实现结构精简、便于直接运行与二次修改。目前已有426人学习下载说明其在教学与练习中具有一定参考价值。读者可从中获得完整的节点建模、初值设定、迭代收敛判断与结果输出流程借助脚本理解PQ、PV与平衡节点的处理方式并对照IEEE标准系统验证算法正确性为后续扩展至更复杂网络打下基础。1. 六节点潮流计算从标题里的“9_bus”说起它到底在算什么很多人第一次看到power_flow_9_bus_6节点潮流这种命名会犯迷糊标题里既有 9 又有 6到底几个节点我一开始也踩过这个坑。常见情况是项目目录或数据集编号沿用了某个 9 节点算例的框架但实际参与潮流求解的母线只有 6 条剩下的是被等值掉的中间节点或纯标注节点。所以真正要跑的是6 节点潮流计算属于电力系统稳态分析里最小可复现的一类算例。它能解决什么问题给定发电机出力、负荷功率和网络阻抗求出每条母线的电压幅值与相角、每条支路的功率分布和网损。适合谁适合刚接触牛顿-拉夫逊法、想手推一遍雅可比矩阵的人也适合做配网重构、微电网调度前先拿小算例验证代码的工程师。6 节点规模小到可以手算校验又足够暴露收敛性、PV/PQ 节点转换、平衡节点选取这些真实问题比 3 节点有代表性比 30 节点好调试。2. 六节点算例怎么建模节点类型、导纳矩阵与数据表2.1 先定节点类型再谈方程潮流计算的本质是解一组非线性代数方程。6 个节点里必须有一个平衡节点Slack负责兜住全网功率不平衡量其余节点分为 PQ 节点注入有功、无功已知和 PV 节点注入有功、电压幅值已知。节点类型定错雅可比矩阵就奇异迭代直接发散。我一般这样分配节点 1 设为平衡节点电压幅值 1.05 p.u.、相角 0节点 2、3 设为 PV 节点模拟两台带励磁的发电机节点 4、5、6 设为 PQ 节点接负荷。这样 6 节点系统共有 2 个 PV、3 个 PQ、1 个 Slack未知量个数为 2×3PQ 的电压幅值和相角 2×1PV 的相角 8 个方程数正好匹配。提示平衡节点不要选在重负荷中心否则迭代初期功率不平衡量过大容易震荡。2.2 支路参数与节点导纳矩阵6 节点算例的支路通常给成 R、X、B/2 的形式。下面是我常用的数据表结构单位统一为标幺值基准功率 100 MVA。支路首端末端R (p.u.)X (p.u.)B/2 (p.u.)1120.0100.0800.0102140.0150.1200.0123230.0080.0600.0084250.0120.1000.0105360.0100.0900.0096450.0200.1500.0157560.0180.1300.013节点导纳矩阵 Ybus 是 6×6 复数矩阵。对角元是节点自导纳等于连在该节点上所有支路导纳之和加上对地充电电容非对角元是互导纳等于两节点间支路导纳的负值。没有直接相连的节点互导纳为 0。2.3 用 Python 组装 Ybus 的最小代码import numpy as np # 支路数据: [首端, 末端, R, X, B/2]节点编号从 1 开始 branches [ [1, 2, 0.010, 0.080, 0.010], [1, 4, 0.015, 0.120, 0.012], [2, 3, 0.008, 0.060, 0.008], [2, 5, 0.012, 0.100, 0.010], [3, 6, 0.010, 0.090, 0.009], [4, 5, 0.020, 0.150, 0.015], [5, 6, 0.018, 0.130, 0.013], ] n 6 Ybus np.zeros((n, n), dtypecomplex) for f, t, R, X, B2 in branches: z R 1j * X y 1.0 / z i, j f - 1, t - 1 # 对角元累加支路导纳和充电电容 Ybus[i, i] y 1j * B2 Ybus[j, j] y 1j * B2 # 互导纳为负的支路导纳 Ybus[i, j] - y Ybus[j, i] - y print(np.round(Ybus, 4))这段代码的逻辑很直白遍历每条支路把串联导纳y和对地充电电容jB2累加到两端节点的自导纳上再把-y写到互导纳位置。参数说明R、X是串联电阻和电抗B2是半充电电纳单位都是标幺值。如果算例给的是变压器支路还要引入非标准变比那时互导纳不再对称需要额外处理变比列。3. 牛顿-拉夫逊法求解雅可比矩阵与迭代收敛3.1 功率不平衡量怎么算牛顿-拉夫逊法的核心是每次迭代求解修正方程[ΔP; ΔQ] J · [Δθ; ΔV/V]。先根据当前电压估计值算出各节点注入功率再和给定值相减得到不平衡量。对 PQ 节点有功和无功不平衡量都要算对 PV 节点只算有功不平衡量无功由电压幅值约束反推。def calc_power(Ybus, V, theta): # V 为幅值向量, theta 为相角向量(弧度) Vc V * np.exp(1j * theta) S Vc * np.conj(Ybus Vc) return S.real, S.imag这里S V · conj(Ybus · V)是节点注入功率的标准公式。注意conj只作用在电流上别写错位置否则无功符号会整体反掉这是新手最常见的翻车点之一。3.2 雅可比矩阵的四个分块雅可比矩阵按[H N; M L]分块H 是 ΔP 对 Δθ 的偏导N 是 ΔP 对 ΔV/V 的偏导M 是 ΔQ 对 Δθ 的偏导L 是 ΔQ 对 ΔV/V 的偏导。6 节点规模下我建议先手推一遍 H 和 N 的公式再写代码否则调试时根本不知道哪块出错。对 i≠j 的情况H 和 N 的非对角元有简洁形式对角元则要加上自身注入功率项。具体公式在任何电力系统分析教材里都有这里不重复抄重点是你写完后要用数值差分校验一次。# 数值差分校验雅可比矩阵的示例 def numerical_jacobian(func, x, eps1e-6): n len(x) J np.zeros((n, n)) f0 func(x) for k in range(n): xp x.copy() xp[k] eps J[:, k] (func(xp) - f0) / eps return J把解析雅可比和数值差分结果对比误差应在 1e-4 量级以内。如果差很多多半是偏导公式里漏了项或者符号写反。3.3 迭代终止条件与收敛判据我一般用最大功率不平衡量作为收敛判据阈值取 1e-6 p.u.。6 节点算例通常 4 到 6 次迭代收敛。如果超过 15 次还在震荡不要硬跑先检查三件事平衡节点选得对不对、PV 节点无功是否越限、支路参数有没有把 X 写成负的。注意收敛判据用无穷范数还是二范数会影响迭代次数但不影响最终解。工程上习惯用无穷范数因为能直接看出哪个节点最不满足。4. 避坑与排查六节点潮流计算里最容易翻车的五件事4.1 现象迭代发散电压幅值越算越大原因通常是雅可比矩阵对角元符号错误或者平衡节点功率没有正确排除在方程之外。解决方法是先用平启动所有 PQ 节点电压设 1.0∠0°跑一次如果平启动都发散基本可以确定是矩阵组装问题而不是初值问题。4.2 现象PV 节点无功出力远超合理范围原因是 PV 节点在迭代中无功不受约束如果电压设定值不合理无功会跑到几千 Mvar。解决方法是给 PV 节点加无功上下限一旦越限就转换为 PQ 节点把无功固定在限值上重新迭代。这个逻辑在 6 节点算例里很容易验证。4.3 现象网损算出来是负的原因是支路功率方向定义混乱或者充电电容符号搞反。解决方法是统一约定支路功率从首端流向末端为正充电电容作为容性无功发出。算完总网损应该等于所有发电机出力减去所有负荷差值为正。4.4 现象Ybus 矩阵不对称理想情况下 Ybus 应该对称。如果不对称检查是不是有变压器变比没处理或者互导纳赋值时只写了一侧。6 节点纯线路算例里Ybus 必须严格对称不对称就是代码 bug。4.5 现象收敛后电压相角差超过 30 度这说明支路电抗相对有功传输过大或者负荷水平设得太高。6 节点算例如果相角差超过 30 度实际系统已经接近静态稳定极限潮流解虽然数学上存在但工程上不可接受。解决方法是降低负荷或增加并联补偿。5. 从能跑到好用结果校验、网损分解与参数扫描5.1 用功率平衡做第一道校验收敛后第一件事不是看电压而是做功率平衡校验。把所有发电机有功相加减去所有负荷有功再减去支路有功损耗结果应该接近零。无功同理但要记得加上充电电容发出的无功。这个校验能抓出大部分符号错误。# 功率平衡校验 P_gen np.array([0, 1.5, 1.0, 0, 0, 0]) # 各节点发电机有功 P_load np.array([0, 0.5, 0.3, 0.8, 0.6, 0.4]) P_loss 0.0 for f, t, R, X, B2 in branches: i, j f - 1, t - 1 # 用收敛后的电压计算支路电流和损耗 # 此处省略电压取值重点看校验逻辑 pass # 理论上 P_gen.sum() - P_load.sum() - P_loss ≈ 0参数说明P_gen和P_load都是标幺值平衡节点的发电机出力由潮流结果反推不参与给定。校验阈值取 1e-4 即可太小没意义太大说明没收敛好。5.2 网损分解到支路6 节点算例的网损可以逐条支路算公式是I²R。把每条支路电流算出来乘以电阻累加就是总有功损耗。这样做的好处是能看出哪条支路损耗占比最大为后续扩容或无功补偿提供依据。我一般会输出一张支路损耗表按损耗从大到小排序。5.3 参数扫描看收敛边界想真正理解 6 节点算例可以做一组参数扫描固定网络参数逐步增加总负荷观察迭代次数和最低电压的变化。当负荷增加到某个临界点潮流方程无解迭代开始发散。这个临界点就是电压崩溃点对理解 PV 曲线很有帮助。# 负荷倍增扫描 for scale in np.arange(1.0, 2.5, 0.1): P_load_scaled P_load * scale # 调用潮流求解函数记录迭代次数和最低电压 # 当迭代次数超过 20 或电压低于 0.7 时停止这个扫描不需要改网络结构只改负荷比例跑一遍就能看到 6 节点系统的承载极限。我习惯把结果画成曲线横轴是负荷倍数纵轴是最低电压幅值拐点位置一目了然。5.4 一个容易被忽略的技巧初值选择平启动不是唯一选择。如果算例接近重载平启动可能不收敛这时可以用直流潮流结果作为初值相角由直流潮流给出电压幅值全设 1.0。直流潮流是线性的求解极快给出的相角已经接近真实解能显著改善牛顿法的收敛性。这个技巧在 6 节点算例上效果很明显尤其是负荷较重的时候。我自己写潮流程序有个习惯每次改完雅可比矩阵先用数值差分校验一遍再跑正式算例。这个习惯帮我省下了大量排查符号错误的时间。6 节点算例虽小但它是验证代码正确性的最佳试验台别急着上大系统先把这 6 个节点跑稳、跑透。希望帮到你。本文还有配套的精品资源点击获取 六节点潮流计算实战：从Ybus组装到牛顿-拉夫逊法收敛避坑
简介这份资源面向电力系统分析的学习者与研究人员提供基于MATLAB实现的牛顿-拉夫森潮流计算程序覆盖6节点与9节点IEEE标准测试系统可用于理解稳态下电压、功率分布与线路潮流的求解过程适合课程设计、算法验证与科研入门等场景。压缩包共2个文件均为.m脚本整体约2KB分别对应6节点与9节点系统的潮流计算实现结构精简、便于直接运行与二次修改。目前已有426人学习下载说明其在教学与练习中具有一定参考价值。读者可从中获得完整的节点建模、初值设定、迭代收敛判断与结果输出流程借助脚本理解PQ、PV与平衡节点的处理方式并对照IEEE标准系统验证算法正确性为后续扩展至更复杂网络打下基础。1. 六节点潮流计算从标题里的“9_bus”说起它到底在算什么很多人第一次看到power_flow_9_bus_6节点潮流这种命名会犯迷糊标题里既有 9 又有 6到底几个节点我一开始也踩过这个坑。常见情况是项目目录或数据集编号沿用了某个 9 节点算例的框架但实际参与潮流求解的母线只有 6 条剩下的是被等值掉的中间节点或纯标注节点。所以真正要跑的是6 节点潮流计算属于电力系统稳态分析里最小可复现的一类算例。它能解决什么问题给定发电机出力、负荷功率和网络阻抗求出每条母线的电压幅值与相角、每条支路的功率分布和网损。适合谁适合刚接触牛顿-拉夫逊法、想手推一遍雅可比矩阵的人也适合做配网重构、微电网调度前先拿小算例验证代码的工程师。6 节点规模小到可以手算校验又足够暴露收敛性、PV/PQ 节点转换、平衡节点选取这些真实问题比 3 节点有代表性比 30 节点好调试。2. 六节点算例怎么建模节点类型、导纳矩阵与数据表2.1 先定节点类型再谈方程潮流计算的本质是解一组非线性代数方程。6 个节点里必须有一个平衡节点Slack负责兜住全网功率不平衡量其余节点分为 PQ 节点注入有功、无功已知和 PV 节点注入有功、电压幅值已知。节点类型定错雅可比矩阵就奇异迭代直接发散。我一般这样分配节点 1 设为平衡节点电压幅值 1.05 p.u.、相角 0节点 2、3 设为 PV 节点模拟两台带励磁的发电机节点 4、5、6 设为 PQ 节点接负荷。这样 6 节点系统共有 2 个 PV、3 个 PQ、1 个 Slack未知量个数为 2×3PQ 的电压幅值和相角 2×1PV 的相角 8 个方程数正好匹配。提示平衡节点不要选在重负荷中心否则迭代初期功率不平衡量过大容易震荡。2.2 支路参数与节点导纳矩阵6 节点算例的支路通常给成 R、X、B/2 的形式。下面是我常用的数据表结构单位统一为标幺值基准功率 100 MVA。支路首端末端R (p.u.)X (p.u.)B/2 (p.u.)1120.0100.0800.0102140.0150.1200.0123230.0080.0600.0084250.0120.1000.0105360.0100.0900.0096450.0200.1500.0157560.0180.1300.013节点导纳矩阵 Ybus 是 6×6 复数矩阵。对角元是节点自导纳等于连在该节点上所有支路导纳之和加上对地充电电容非对角元是互导纳等于两节点间支路导纳的负值。没有直接相连的节点互导纳为 0。2.3 用 Python 组装 Ybus 的最小代码import numpy as np # 支路数据: [首端, 末端, R, X, B/2]节点编号从 1 开始 branches [ [1, 2, 0.010, 0.080, 0.010], [1, 4, 0.015, 0.120, 0.012], [2, 3, 0.008, 0.060, 0.008], [2, 5, 0.012, 0.100, 0.010], [3, 6, 0.010, 0.090, 0.009], [4, 5, 0.020, 0.150, 0.015], [5, 6, 0.018, 0.130, 0.013], ] n 6 Ybus np.zeros((n, n), dtypecomplex) for f, t, R, X, B2 in branches: z R 1j * X y 1.0 / z i, j f - 1, t - 1 # 对角元累加支路导纳和充电电容 Ybus[i, i] y 1j * B2 Ybus[j, j] y 1j * B2 # 互导纳为负的支路导纳 Ybus[i, j] - y Ybus[j, i] - y print(np.round(Ybus, 4))这段代码的逻辑很直白遍历每条支路把串联导纳y和对地充电电容jB2累加到两端节点的自导纳上再把-y写到互导纳位置。参数说明R、X是串联电阻和电抗B2是半充电电纳单位都是标幺值。如果算例给的是变压器支路还要引入非标准变比那时互导纳不再对称需要额外处理变比列。3. 牛顿-拉夫逊法求解雅可比矩阵与迭代收敛3.1 功率不平衡量怎么算牛顿-拉夫逊法的核心是每次迭代求解修正方程[ΔP; ΔQ] J · [Δθ; ΔV/V]。先根据当前电压估计值算出各节点注入功率再和给定值相减得到不平衡量。对 PQ 节点有功和无功不平衡量都要算对 PV 节点只算有功不平衡量无功由电压幅值约束反推。def calc_power(Ybus, V, theta): # V 为幅值向量, theta 为相角向量(弧度) Vc V * np.exp(1j * theta) S Vc * np.conj(Ybus Vc) return S.real, S.imag这里S V · conj(Ybus · V)是节点注入功率的标准公式。注意conj只作用在电流上别写错位置否则无功符号会整体反掉这是新手最常见的翻车点之一。3.2 雅可比矩阵的四个分块雅可比矩阵按[H N; M L]分块H 是 ΔP 对 Δθ 的偏导N 是 ΔP 对 ΔV/V 的偏导M 是 ΔQ 对 Δθ 的偏导L 是 ΔQ 对 ΔV/V 的偏导。6 节点规模下我建议先手推一遍 H 和 N 的公式再写代码否则调试时根本不知道哪块出错。对 i≠j 的情况H 和 N 的非对角元有简洁形式对角元则要加上自身注入功率项。具体公式在任何电力系统分析教材里都有这里不重复抄重点是你写完后要用数值差分校验一次。# 数值差分校验雅可比矩阵的示例 def numerical_jacobian(func, x, eps1e-6): n len(x) J np.zeros((n, n)) f0 func(x) for k in range(n): xp x.copy() xp[k] eps J[:, k] (func(xp) - f0) / eps return J把解析雅可比和数值差分结果对比误差应在 1e-4 量级以内。如果差很多多半是偏导公式里漏了项或者符号写反。3.3 迭代终止条件与收敛判据我一般用最大功率不平衡量作为收敛判据阈值取 1e-6 p.u.。6 节点算例通常 4 到 6 次迭代收敛。如果超过 15 次还在震荡不要硬跑先检查三件事平衡节点选得对不对、PV 节点无功是否越限、支路参数有没有把 X 写成负的。注意收敛判据用无穷范数还是二范数会影响迭代次数但不影响最终解。工程上习惯用无穷范数因为能直接看出哪个节点最不满足。4. 避坑与排查六节点潮流计算里最容易翻车的五件事4.1 现象迭代发散电压幅值越算越大原因通常是雅可比矩阵对角元符号错误或者平衡节点功率没有正确排除在方程之外。解决方法是先用平启动所有 PQ 节点电压设 1.0∠0°跑一次如果平启动都发散基本可以确定是矩阵组装问题而不是初值问题。4.2 现象PV 节点无功出力远超合理范围原因是 PV 节点在迭代中无功不受约束如果电压设定值不合理无功会跑到几千 Mvar。解决方法是给 PV 节点加无功上下限一旦越限就转换为 PQ 节点把无功固定在限值上重新迭代。这个逻辑在 6 节点算例里很容易验证。4.3 现象网损算出来是负的原因是支路功率方向定义混乱或者充电电容符号搞反。解决方法是统一约定支路功率从首端流向末端为正充电电容作为容性无功发出。算完总网损应该等于所有发电机出力减去所有负荷差值为正。4.4 现象Ybus 矩阵不对称理想情况下 Ybus 应该对称。如果不对称检查是不是有变压器变比没处理或者互导纳赋值时只写了一侧。6 节点纯线路算例里Ybus 必须严格对称不对称就是代码 bug。4.5 现象收敛后电压相角差超过 30 度这说明支路电抗相对有功传输过大或者负荷水平设得太高。6 节点算例如果相角差超过 30 度实际系统已经接近静态稳定极限潮流解虽然数学上存在但工程上不可接受。解决方法是降低负荷或增加并联补偿。5. 从能跑到好用结果校验、网损分解与参数扫描5.1 用功率平衡做第一道校验收敛后第一件事不是看电压而是做功率平衡校验。把所有发电机有功相加减去所有负荷有功再减去支路有功损耗结果应该接近零。无功同理但要记得加上充电电容发出的无功。这个校验能抓出大部分符号错误。# 功率平衡校验 P_gen np.array([0, 1.5, 1.0, 0, 0, 0]) # 各节点发电机有功 P_load np.array([0, 0.5, 0.3, 0.8, 0.6, 0.4]) P_loss 0.0 for f, t, R, X, B2 in branches: i, j f - 1, t - 1 # 用收敛后的电压计算支路电流和损耗 # 此处省略电压取值重点看校验逻辑 pass # 理论上 P_gen.sum() - P_load.sum() - P_loss ≈ 0参数说明P_gen和P_load都是标幺值平衡节点的发电机出力由潮流结果反推不参与给定。校验阈值取 1e-4 即可太小没意义太大说明没收敛好。5.2 网损分解到支路6 节点算例的网损可以逐条支路算公式是I²R。把每条支路电流算出来乘以电阻累加就是总有功损耗。这样做的好处是能看出哪条支路损耗占比最大为后续扩容或无功补偿提供依据。我一般会输出一张支路损耗表按损耗从大到小排序。5.3 参数扫描看收敛边界想真正理解 6 节点算例可以做一组参数扫描固定网络参数逐步增加总负荷观察迭代次数和最低电压的变化。当负荷增加到某个临界点潮流方程无解迭代开始发散。这个临界点就是电压崩溃点对理解 PV 曲线很有帮助。# 负荷倍增扫描 for scale in np.arange(1.0, 2.5, 0.1): P_load_scaled P_load * scale # 调用潮流求解函数记录迭代次数和最低电压 # 当迭代次数超过 20 或电压低于 0.7 时停止这个扫描不需要改网络结构只改负荷比例跑一遍就能看到 6 节点系统的承载极限。我习惯把结果画成曲线横轴是负荷倍数纵轴是最低电压幅值拐点位置一目了然。5.4 一个容易被忽略的技巧初值选择平启动不是唯一选择。如果算例接近重载平启动可能不收敛这时可以用直流潮流结果作为初值相角由直流潮流给出电压幅值全设 1.0。直流潮流是线性的求解极快给出的相角已经接近真实解能显著改善牛顿法的收敛性。这个技巧在 6 节点算例上效果很明显尤其是负荷较重的时候。我自己写潮流程序有个习惯每次改完雅可比矩阵先用数值差分校验一遍再跑正式算例。这个习惯帮我省下了大量排查符号错误的时间。6 节点算例虽小但它是验证代码正确性的最佳试验台别急着上大系统先把这 6 个节点跑稳、跑透。希望帮到你。本文还有配套的精品资源点击获取