《电力系统分析》课程.doc

电力系统分析电力系统分析 课程报告课程报告 题目 基于 Matlab 电力系统潮流计算 系系 别别 电气工程系 专业班级专业班级 学生姓名学生姓名 学生学号学生学号 指导教师指导教师 提交日期提交日期 2011 年 12 月 18 日 华南理工大学广州学院电气工程系课程报告 I 目 录 一 电力系统潮流计算机计算的意义和目的一 电力系统潮流计算机计算的意义和目的 1 1 1 11 1 潮流计算机计算的意义潮流计算机计算的意义 1 1 1 21 2 潮流计算机计算的目的潮流计算机计算的目的 2 2 1 31 3 设计内容设计内容 2 2 二 潮流计算的基本原理二 潮流计算的基本原理 2 2 2 12 1 潮流计算简介潮流计算简介 2 2 2 22 2 潮流计算方法潮流计算方法 3 3 2 32 3 M MATLABATLAB简介简介 7 7 三 潮流计算机计算的流程图三 潮流计算机计算的流程图 8 8 3 13 1 潮流计算流程图潮流计算流程图 8 8 3 23 2 潮流计算源程序图潮流计算源程序图 9 9 3 2 13 2 1 三机九节点系统三机九节点系统 9 9 3 33 3 运行计算结果及分析运行计算结果及分析 1515 四 总结四 总结 2323 五 参考文献五 参考文献 2525 附录 算例原始参数 25 华南理工大学广州学院电气工程系课程报告 1 一 电力系统潮流计算机计算的意义和目的 1 11 1 潮流计算机计算的意义潮流计算机计算的意义 潮流计算是电力系统的一项重要分析功能 是进行故障计算 继电保护整定 安 全分析的必要工具 电力系统已经与我们的生活息息相关 不可分割 进行电力系统潮流计算是保证 电力系统正常运行的必要计算 具体来讲电力系统潮流计算具有以下意义 1 在电网规划阶段 通过潮流计算 合理规划电源容量及接入点 合理规划网架 选择 无功补偿方案 满足规划水平的大 小方式下潮流交换控制 调峰 调相 调压的要求 2 在编制年运行方式时 在预计负荷增长及新设备投运基础上 选择典型方式进行潮 流计算 发现电网中薄弱环节 供调度员日常调度控制参考 并对规划 基建部门提出改进 网架结构 加快基建进度的建议 3 正常检修及特殊运行方式下的潮流计算 用于日运行方式的编制 指导发电厂开机 方式 有功 无功调整方案及负荷调整方案 满足线路 变压器热稳定要求及电压质量要 求 4 预想事故 设备退出运行对静态安全的影响分析及作出预想的运行方式调整方 案 总结为在电力系统运行方式和规划方案的研究中 都需要进行潮流计算以比较运 行方式或规划供电方案的可行性 可靠性和经济性 同时 为了实时监控电力系统的 运行状态 也需要进行大量而快速的潮流计算 因此 潮流计算是电力系统中应用最广泛 最基本和最重要的一种电气运算 基于电力系统计算对保证电力系统正常运行具有如此正要的意义 这就要求我们 能够快速准确的进行潮流计算 计算机技术的发展使电力系统机辅分析成为可能 各 种潮流计算软件也相继出现 MATLAB 使用方便 有着其他高级语言无法比拟的强大 的矩阵处理功能 MATLAB 拥有 600 多个工程数学运算函数 可实现潮流计算的矩阵 求积 求逆 稀疏矩阵形成 复数运算以及初等数学运算 同时 MATLAB 语言允许用 户以数学形式的语言编写程序 这样编程的工作量就大为减少 要达到较高的计算精 度 且兼顾矩阵程序设计的难易程度 使 MATLAB 成为首选潮流计算的计算机语言 因此本次设计提出了基于 MATLAB 潮流计算软件的分析与设计 该软件能快速准确 华南理工大学广州学院电气工程系课程报告 2 的对电力系统潮流进行计算 并具有一定的辅助分析功能 通过电力系统潮流计算课程报告该环节 使学生熟悉电气工程中主要电力设备的 特性 数学模型 相互关系及计算方法 为进一步掌握和研究电气工程规划 设计和 运行等问题打下良好的基础 1 21 2 潮流计算机计算的目的潮流计算机计算的目的 电力系统潮流计算机计算的目的 1 掌握电力系统潮流计算的基本原理 2 掌握并能熟练运用一门计算机语言 MATLAB 语言或 C 语言或 C 语言 3 采用计算机语言对潮流计算进行计算机编程 1 31 3 设计内容设计内容 1 根据电力系统网络推导电力网络数学模型 写出节点导纳矩阵 2 赋予各节点电压变量 直角坐标系形式 初值后 求解不平衡量 3 形成雅可比矩阵 4 求解修正量后 重新修改初值 从 2 开始重新循环计算 5 求解的电压变量达到所要求的精度时 再计算各支路功率分布 功率损耗和 平衡节点功率 6 上机编程调试 7 书写课程报告 二 潮流计算的基本原理 2 12 1 潮流计算简介潮流计算简介 利用电子计算机进行潮流计算从 20 世纪 50 年代中期就已经开始 此后 潮流计 算曾采用了各种不同的方法 这些方法的发展主要是围绕着对潮流计算的一些基本要 求进行的 电力系统潮流计算属于稳态分析范畴 不涉及系统元件的动态特性和过渡 过程 因此其数学模型不包含微分方程 是一组高阶非线性方程 非线性代数方程组 的解法离不开迭代 因此 潮流计算方法首先要求它是能可靠的收敛 并给出正确答 案 随着电力系统规模的不断扩大 潮流问题的方程式阶数越来越高 目前已达到几 千阶甚至上万阶 对这样规模的方程式并不是采用任何数学方法都能保证给出正确答 华南理工大学广州学院电气工程系课程报告 3 案的 这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法 知道现在 潮流算法的研究仍然非常活跃 但是大多数研究都是围绕改进牛顿法和 P Q 分解法进 行的 此外 随着人工智能理论的发展 遗传算法 人工神经网络 模糊算法也逐渐 被引入潮流计算 但是 到目前为止这些新的模型和算法还不能取代牛顿法和 P Q 分 解法的地位 由于电力系统规模的不断扩大 对计算速度的要求不断提高 计算机的 并行计算技术也将在潮流计算中得到广泛的应用 成为重要的研究领域 通过几十年的发展 潮流算法日趋成熟 近几年 对潮流算法的研究仍然是如何 改善传统的潮流算法 即高斯 塞德尔法 牛顿法和快速解耦法 牛顿法 由于其在求 解非线性潮流方程时采用的是逐次线性化的方法 为了进一步提高算法的收敛性和计 算速度 人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来 于是产生了二阶 潮流算法 后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点 提出了采用直角坐标的保留非线性快速潮流算法 岩本伸一等提出了一种保留非线性 的快速潮流计算法 但用的是指教坐标系 因而没法利用 P Q 解耦 为了更有利于大电 网的潮流计算 将此原理推广用于 P Q 解耦 这样 既利用了保留非线性的快速算法 在 迭代中使用常数雅克比矩阵 又保留了 P Q 解耦的优点 对于一些病态系统 应用非 线性潮流计算方法往往会造成计算过程的振荡或者不收敛 从数学上讲 非线性的潮 流计算方程组本来就是无解的 这样 人们提出来了将潮流方程构造成一个函数 求 此函数的最小值问题 称之为非线性规划潮流的计算方法 优点是原理上保证了计算 过程永远不会发散 如果将数学规划原理和牛顿潮流算法有机结合一起就是最优乘子 法 另外 为了优化系统的运行 从所有以上的可行潮流解中挑选出满足一定指标要 求的一个最佳方案就是最优潮流问题 最优潮流是一种同时考虑经济性和安全性的电 力网络分析优化问题 OPF 在电力系统的安全运行 经济调度 可靠性分析 能量管 理以及电力定价等方面得到了广泛的应用 可信域和线性搜索方法是保证最优化算法 全局收敛性能的两类技术 将内点法和可信域 线性搜索方法有机结合 构造新的优 化算法 是数学规划领域的研究热点 对于一些特殊性质的潮流计算问题有直流潮流 计算方法 随机潮流计算方法和三相潮流计算方法 2 22 2 潮流计算方法潮流计算方法 2 2 1牛顿 拉夫逊法概述 电力系统潮流计算是电力系统分析中的一种最基本的计算 是对复杂电力系统正 常和故障条件下稳态运行状态的计算 潮流计算的目标是求取电力系统在给定运行状 华南理工大学广州学院电气工程系课程报告 4 态的计算 即节点电压和功率分布 用以检查系统各元件是否过负荷 各点电压是否 满足要求 功率的分布和分配是否合理以及功率损耗等 对现有电力系统的运行和扩 建 对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮 流计算为基础 潮流计算结果可用如电力系统稳态研究 安全估计或最优潮流等对潮 流计算的模型和方法有直接影响 实际电力系统的潮流技术那主要采用牛顿 拉夫逊法 牛顿 拉夫逊法 简称牛顿法 在数学上是求解非线性代数方程式的有效方法 其要 点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程 即 通常所称的逐次线性化过程 2 2 2 高斯 赛德尔法 高斯 塞德尔法原理比较简单 主要以节点导纳矩阵为基础 下面简单介绍下其原 理和潮流计算过程 1 高斯 塞德尔法的基本原理 设有 n 个联立的非线性方程 0 0 0 21 212 211 nn n n xxxf xxxf xxxf 2 1 解此方程组可得 nnn n n xxxgx xxxgx xxxgx 21 2122 2111 2 2 若已经求得各变量的第 k 此迭代值 2 1 k n kk xxx 则第 k 1 次迭代值为 1 1 2 1 11 1 2 1 11 1 2 2 11 1 1 k n kkk n k n kkk k n kkk xxxgx xxxgx xxxgx 2 3 只要给定变量的初值 0 0 2 0 1 n xxx 就可以按式 2 10 迭代计算 一直进行到所 有变量都满足收敛条件 1 k i k i xx 即可 华南理工大学广州学院电气工程系课程报告 5 2 高斯 塞德尔潮流计算过程 假设有 n 个节点的电力系统 没有 PV 节点 平衡节点编号为 s 功率方程可写成下列 复数方程式 n ij j jij i ii ii i UY U jQP Y U 1 1 2 4 对每一个 PQ 节点都可列出一个方程式 因而有 n 1 个方程式 在这些方程式中 注入功率 i P 和 i Q 都是给定的 平衡节点电压也是已知的 因而只有 n 1 个节点的电压 为未知量 从而有可能求得唯一解 将上式写成高斯 塞德尔法的迭代形式 1 11 1 1 1 i j n ij k jij k jij k i ii ii k i UYUY U jQP Y U 2 5 如系统内存在 PV 节点 假设节点 p 为 PV 节点 设定的节点电压为 Up0 假定高斯 塞德尔迭代法已完成第 k 次迭代 接着要做第 k 1 次迭代前 先按下式求出节点 p 的 注入无功功率 Im 1 1 k j n j pj K P k p UYUQ 2 6 然后代入下式 求出 p 点电压 n pj j k jpj k p k pp pp k p UY U jQP Y U i 1 1 1 1 2 7 在迭代过程中 按上式求得的节点 p 的电压大小不一定等于设定的节点电压 Up0 所有在下一次的迭代中 应以设定的 Up0 对电压进行修正 但其相角仍保持上式所求 得的值 使得 1 0 1 k pp k p UU 2 8 如果所求得 PV 节点的无功功率越限 则无功功率在限 该 PV 节点转化为 PQ 节点 归纳起来 高斯 塞德尔迭代法计算潮流的步骤为 1 设定各节点电压的初值 并给定迭代误差判据 2 对每一个 PQ 节点 以前一次迭代的节点电压值代入功率迭代方程式求出新值 华南理工大学广州学院电气工程系课程报告 6 3 对于 PV 节点 求出其无功功率 并判断是否越限 如越限则将 PV 节点转化为 PQ 节点 4 判别各节点电压前后二次迭代值相量差的模是否小于给定误差 如不小于 则回 到第 2 步 继续进行计算 否则转到第 5 步 5 根据功率方程求出平衡节点注入功率 6 求支路功率分布和支路功率损耗 2 2 3 PQ 分解法 PQ 分解法是牛顿法的一种简化方法 它利用了电力系统特有的运行特性 改进和 提高了运行速度 由牛顿法的修正方程进行展开可得 U U LKQ U U NP 2 9 根据电力系统的运行特性进行简化 1 考虑到电力系统中有功功率分布主要受节点电压相角的影响 无功功率分布主要 受节点电压幅值的影响 所以可以近似的忽略电压幅值变化对有功功率和电压相 位变化对无功功率分布的影响 即 U U LQHP KN 0 0 2 10 2 根据电力系统的正常运行条件还可作下列假设 1 电力系统正常运行时线路两端的电压相位角一般变化不大 不超过 10 20 度 2 电力系统中一般架空线路的电抗远大于电阻 3 节点无功功率相应的导纳 Q U U 远小于该节点的自导纳的虚部 用算式表示如下 iiii ijijij ij BUQ BG 2 sin 1cos 2 11 由以上假设 可得到雅克比矩阵的表达式 iiiiiii ijjiijij BUL BUUL 2 2 12 华南理工大学广州学院电气工程系课程报告 7 修正方程式为 UUB U U UBUQ UBUP 2 13 U 为节点电压有效值的对角矩阵 B 为电纳矩阵 由节点导纳矩阵中各元素的虚部 构成 根据不同的节点还要做一些改变 1 在有功功率部分 要除去与有功功率和电压相位关系较小的因素 如不包含各 输电线路和变压器支路等值 型电路的对地电纳 2 在无功功率部分 PV 节点要做相应的处理 则修正方程表示为 UBQU UBPU 1 1 2 14 一般 由于以上原因 B 和 B 是不相同的 但都是对称的常数矩阵 PQ 分解法的特点 1 以一个 n 1 阶和一个 n m 1 阶线性方程组代替原有的 2n m 1 阶线性方程组 2 修正方程的系数矩阵 B 和 B 为对称常数矩阵 且在迭代过程中保持不变 3 P Q 分解法具有线性收敛特性 与牛顿 拉夫逊法相比 当收敛到同样的精度时 需要的迭代次数较多 4 P Q 分解法一般只适用于 110KV 及以上电网的计算 因为 35KV 及以下电压等级 的线路 r x 比值很大 不满足上述简化条件 可能出现迭代计算不收敛的情况 2 4 4 拟牛顿算法 拟牛顿法是从牛顿法派生出来的新的算法 它一出现就引起广泛的重视 近年 来 拟牛顿法的研究十分活跃 它成为解非线性方程组及优化问题的重要方法 它能 在计算电力系统的潮流分布中 成功地减少每步迭代的计算量 并保持着超线性收敛 速度 2 32 3 MatlabMatlab 简介简介 2 3 1 Matlab 概述 MATLAB Matrix Laboratory 为美国 Mathworks 公司 1983 年首次推出的一套高性 能的数值分析和计算软件 其功能不断扩充 版本不断升级 MATLAB 将矩阵运算 数 华南理工大学广州学院电气工程系课程报告 8 值分析 图形处理 编程技术结合在一起 为用户提供了一个强有力的科学及工程问 题的分析计算和程序设计工具 它还提供了专业水平的符号计算 文字处理 可视化 建模仿真和实时控制等功能 是具有全部语言功能和特征的新一代软件开发平台 MATLAB 具有编程效率高 用户使用方便 扩充能力强 语句简单 内涵丰富 高效方 便的矩阵和数组运算 方便的绘图功能等特点 给用户带来了极大的方便 2 3 1 matlab GUI 简介 图形用户界面 GUI 是用户与计算机程序之间的交互方式 是用户与计算机进行 信息交流的方式 计算机在屏幕显示图形和文本 若有扬声器还可产生 声音 用户通 过输入设备 如 键盘 鼠标 跟踪球 绘制板或麦克风 与计算机通讯 用户界面 设定了如何观看和如何感知计算机 操作系统或应用程序 通常 多是根据悦目的结 构和用户界面功能的有效性来选择计算机或程序 图形用户界面或 GUI 是包含图形对 象 如 窗口 图标 菜单和文本的用户界面 以某种方式 选择或激活这些对象 通 常引起动作或发生变化 最常见的激活方法是用鼠标或其它点击设备去控制屏幕上的 鼠标指针的运动 按下鼠标按钮 标志着对象的选择或 其它动作 Matlab 作为强大的数学计算软件 同样也提供了图像用户界面设计的功能 在 matlab 中 基本的图形用户界面对象包含 3 类 用户控件对象 uicontrol 下拉式 菜单对象 uimenu 和快捷菜单对象 uicontexmenu 根据这些对象可以设计出界 面友好 操作方便的图形用户界面 三 潮流计算机计算的流程图 3 13 1 潮流计算流程图潮流计算流程图 华南理工大学广州学院电气工程系课程报告 9 图 3 1 潮流计算机计算流程图 3 23 2 潮流计算源程序图潮流计算源程序图 3 2 13 2 1 三机九节点系统三机九节点系统 3 2 1 1 主函数 Sbase MVA 100 fid fopen Nodedata txt N textscan fid s u d f f f f f f fclose fid busnumber size N 1 1 for i 1 busnumber Bus i name N 1 i Bus i type N 2 i Bus i no i Bus i Base KV N 3 i Bus i PG N 4 i Bus i QG N 5 i 华南理工大学广州学院电气工程系课程报告 10 Bus i PL N 6 i Bus i QL N 7 i Bus i pb N 8 i Bus i V 1 0 Bus i angle 0 end fid fopen Aclinedata txt A textscan fid s s f f f f fclose fid aclinenumber size A 1 1 for i 1 aclinenumber Acline i fbname A 1 i Acline i tbname A 2 i Acline i Base KV A 3 i Acline i R A 4 i Acline i X A 5 i Acline i hB A 6 i for k 1 busnumber if strcmp Acline i fbname Bus k name Acline i fbno Bus k no end if strcmp Acline i tbname Bus k name Acline i tbno Bus k no end end end fid fopen Transdata txt T textscan fid s f f s f f f f fclose fid tansnumber size T 1 1 for i 1 tansnumber Trans i fbname T 1 i Trans i fbBase KV T 2 i Trans i fbrated KV T 3 i Trans i tbname T 4 i Trans i tbBase KV T 5 i Trans i tbrated KV T 6 i Trans i R T 7 i Trans i X T 8 i for k 1 busnumber if strcmp Trans i fbname Bus k name Trans i fbno Bus k no end if strcmp Trans i tbname Bus k name 华南理工大学广州学院电气工程系课程报告 11 Trans i tbno Bus k no end end Trans i k Trans i tbrated KV Trans i fbBase KV Trans i fbrated KV Trans i tbBase KV tempx Trans i fbrated KV 2 Trans i fbBase KV 2 Trans i X tempx Trans i X Trans i R tempx Trans i R end N 0 Trans 1 Trans 2 for Y G jB matrix G B B2 FormYmatrix Bus busnumber Acline aclinenumber Trans tansnumber B B B2 B dlmwrite Gmatrix txt G delimiter t precision 6 dlmwrite Bmatrix txt B delimiter t precision 6 G B B2 pause JP JQ FormJPQmatrix Bus B B2 busnumber JP iJP inv JP JQ iJQ inv JQ pause maxiteration 0 for i 1 busnumber NodeV i Bus i V Nodea i Bus i angle VX i Bus i V cos Bus i angle VY i Bus i V sin Bus i angle dQGQL i Bus i QG Bus i QL dPGPL i Bus i PG Bus i PL end NodeV NodeV Nodea Nodea VX VX VY VY dQGQL dQGQL dPGPL dPGPL pause for nointer 1 10 maxdP 1 华南理工大学广州学院电气工程系课程报告 12 maxdQ 1 epsilon 0 000001 noiteration 0 while maxdP epsilon deltaP deltaQ maxdP maxdQ da iJP deltaP dV iJQ deltaQ Nodea Nodea da NodeV NodeV dV noiteration noiteration 1 if noiteration 20 break end end for i 1 busnumber Bus i V NodeV i NodeV i NodeV i Bus i Base KV Bus i angle Nodea i Nodea i Nodea i 180 pi end noiteration Nodea Nodea NodeV NodeV Clear 3 2 1 2 子函数 生成 G B 矩阵 function G B X FormYmatrix Bus busnumber Acline aclinenumber Trans tansnumber Y zeros busnumber X zeros busnumber for i 1 busnumber Y i i Y i i Bus i pb j end for i 1 aclinenumber f Acline i fbno t Acline i tbno Y f f Y f f Acline i hB j 1 Acline i R Acline i X j Y t t Y t t Acline i hB j 1 Acline i R Acline i X j Y f t Y f t 1 Acline i R Acline i X j Y t f Y t f 1 Acline i R Acline i X j 华南理工大学广州学院电气工程系课程报告 13 X f f X f f 1 Acline i X X t t X t t 1 Acline i X X f t 1 Acline i X X t f 1 Acline i X end for i 1 tansnumber f Trans i fbno t Trans i tbno Y f f Y f f 1 Trans i R Trans i X j Y t t Y t t 1 Trans i R Trans i X j Trans i k 2 Y f t Y f t 1 Trans i R Trans i X j Trans i k Y t f Y t f 1 Trans i R Trans i X j Trans i k X f f X f f 1 Trans i X X t t X t t 1 Trans i X X f t 1 Trans i X X t f 1 Trans i X end G real Y B imag Y end 生成 JP JQ 矩阵 function JP JQ FormJPQmatrix Bus B B2 busnumber JP B JQ B2 for i 1 busnumber if Bus i type 1 for k 1 busnumber JQ i k 0 JQ k i 0 JP i k 0 JP k i 0 end JQ i i 1 JP i i 1 end if Bus i type 3 for k 1 busnumber JQ i k 0 JQ k i 0 end JQ i i 1 end end 华南理工大学广州学院电气工程系课程报告 14 end 计算偏节点 PQ 差量 function deltaP deltaQ maxdP maxdQ FormdPQvector Bus NodeV Nodea dQGQL dPGPL B G busnumber deltaQ dQGQL deltaP dPGPL maxdP 0 maxdQ 0 for i 1 busnumber if Bus i type 1 deltaQ i 0 deltaP i 0 end if Bus i type 3 deltaQ i 0 y1 0 y2 0 y3 0 for k 1 busnumber if B i k 0 G i k 0 y1 y1 G i k VX k B i k VY k y2 y2 G i k VY k B i k VX k y3 y3 NodeV k G i k cos Nodea i Nodea k B i k sin Nodea i Nodea k end end deltaP i deltaP i y3 NodeV i deltaP2 i deltaP2 i y1 VX i y2 VY i Bus i V end if Bus i type 2 y1 0 y2 0 y3 0 y4 0 for k 1 busnumber if B i k 0 G i k 0 y1 y1 G i k VX k B i k VY k y2 y2 G i k VY k B i k VX k y3 y3 NodeV k G i k cos Nodea i Nodea k B i k sin Nodea i Nodea k y4 y4 NodeV k G i k sin Nodea i Nodea k B i k cos Nodea i Nodea k end end deltaP i deltaP i y3 NodeV i deltaP2 i deltaP2 i y1 VX i y2 VY i Bus i V 华南理工大学广州学院电气工程系课程报告 15 deltaQ i deltaQ i y4 NodeV i deltaQ2 i deltaQ2 i y1 VY i y2 VX i Bus i V end if maxdP abs deltaP i maxdP abs deltaP i end if maxdQ abs deltaQ i maxdQ abs deltaQ i end deltaP i deltaP i NodeV i deltaQ i deltaQ i NodeV i end end 3 33 3 运行计算结果及分析运行计算结果及分析 3 3 1 三机九节点系统 3 3 1 1 输入数据 节点数据 Nodedata txt bus1 118 0 0 0 0 0 bus2318 1 63 0 0 0 0 bus33 18 0 85 0 0 0 0 bus42 230 0 0 0 0 0 bus52 230 0 0 1 25 0 5 0 bus62 230 0 0 0 9 0 30 bus72 230 0 0 0 0 0 bus82 230 0 0 1 0 0 35 0 bus92 230 0 0 0 0 0 支路数据 Aclinedata txt bus4 bus5 230 0 01 0 085 0 088 bus4 bus6 230 0 017 0 092 0 079 bus5 bus7 230 0 032 0 161 0 153 bus6 bus9 230 0 039 0 17 0 179 bus7 bus8 230 0 0085 0 072 0 0745 bus8 bus9 230 0 0119 0 1008 0 1045 变压器数据 Transdata txt bus1 18 0 18 0 bus4 230 230 0 0 0 0576 bus2 18 0 18 0 bus7 230 230 0 0 0 0625 bus3 18 0 18 0 bus9 230 230 0 0 0 0586 3 3 1 2 输出数据 Sbase MVA 100 N 9x1 cell 9x1 uint32 9x1 int32 9x1 double 9x1 double 9x1 double 华南理工大学广州学院电气工程系课程报告 16 9x1 double 9x1 double 8x1 double busnumber 9 A 6x1 cell 6x1 cell 6x1 double 6x1 double 6x1 double 6x1 double aclinenumber 6 T 3x1 cell 3x1 double 3x1 double 3x1 cell 3x1 double 3x1 double 3x1 double 3x1 double tansnumber 3 华南理工大学广州学院电气工程系课程报告 17 NodeV 1 1 1 1 1 1 1 1 1 Nodea 0 0 0 0 0 0 华南理工大学广州学院电气工程系课程报告 18 0 0 0 dQGQL 0 0 0 0 0 5000 0 3000 0 0 3500 0 dPGPL 0 1 6300 0 8500 0 1 2500 0 9000 0 1 0000 0 noiteration 9 Nodea 0 9 6687 4 7711 2 4066 4 3499 4 0173 3 7991 0 6215 1 9256 NodeV 18 18 18 227 220 224 229 227 231 3 3 2 习题 11 3 3 3 2 1 输出结果 maxd 1 k 1 PQ2 10 5000 0 0 10 0000 0 6000 0 4500 10 0000 0 4000 0 3000 10 0000 0 4000 0 2800 10 0000 0 6000 0 4000 10 0000 0 4000 0 3000 10 0000 0 5000 0 3500 华南理工大学广州学院电气工程系课程报告 19 10 0000 0 5000 0 4000 x 0 0024 x 0 0025 x 0 0037 x 0 0098 x 0 0337 x 0 0041 x 0 1817 maxd 0 1553 V2 10 5000 10 1553 9 9206 9 9661 9 8091 9 9432 9 9550 9 9346 maxd 0 1909 k 2 V 10 5000 10 1553 9 9206 9 9661 9 8091 9 9432 9 9550 9 9346 PQ2 10 5000 0 0 10 0000 0 6000 0 4500 10 0000 0 4000 0 3000 10 0000 0 4000 0 2800 10 0000 0 6000 0 4000 10 0000 0 4000 0 3000 10 0000 0 5000 0 3500 10 0000 0 5000 0 4000 x 0 0024 x 0 0025 x 0 0038 x 0 0100 x 0 0351 华南理工大学广州学院电气工程系课程报告 20 x 0 0042 x 0 1763 maxd 2 3602e 004 V2 10 5000 10 1556 10 0772 9 8865 9 9672 9 7512 9 7633 10 0909 maxd 0 1565 maxd 0 1581 maxd 0 1920 k 3 V 10 5000 10 1556 10 0772 9 8865 9 9672 9 7512 9 7633 10 0909 PQ2 10 5000 0 0 10 0000 0 6000 0 4500 10 0000 0 4000 0 3000 10 0000 0 4000 0 2800 10 0000 0 6000 0 4000 10 0000 0 4000 0 3000 10 0000 0 5000 0 3500 10 0000 0 5000 0 4000 x 0 0024 x 0 0026 x 0 0039 x 0 0096 x 0 0340 x 0 0040 x 0 1762 maxd 1 0915e 004 V2 华南理工大学广州学院电气工程系课程报告 21 10 5000 10 1557 10 0774 10 0435 9 9675 9 9102 9 9221 10 0911 maxd 2 5488e 004 maxd 0 1570 maxd 0 1590 k 4 V 10 5000 10 1557 10 0774 10 0435 9 9675 9 9102 9 9221 10 0911 PQ2 10 5000 0 0 10 0000 0 6000 0 4500 10 0000 0 4000 0 3000 10 0000 0 4000 0 2800 10 0000 0 6000 0 4000 10 0000 0 4000 0 3000 10 0000 0 5000 0 3500 10 0000 0 5000 0 4000 x 0 0024 x 0 0025 x 0 0038 x 0 0096 x 0 0340 x 0 0040 x 0 1761 maxd 2 0435e 005 V2 10 5000 10 1557 10 0775 10 0438 9 9677 9 9105 9 9224 10 0912 maxd 1 1416e 004 maxd 2 5935e 004 maxd 3 4787e 004 华南理工大学广州学院电气工程系课程报告 22 k 5 V 10 5000 10 1557 10 0775 10 0438 9 9677 9 9105 9 9224 10 0912 PQ2 10 5000 0 0 10 0000 0 6000 0 4500 10 0000 0 4000 0 3000 10 0000 0 4000 0 2800 10 0000 0 6000 0 4000 10 0000 0 4000 0 3000 10 0000 0 5000 0 3500 10 0000 0 5000 0 4000 x 0 0024 x 0 0025 x 0 0038 x 0 0096 x 0 0340 x 0 0040 x 0 1761 maxd 1 8943e 007 V2 10 5000 10 1557 10 0776 10 0439 9 9677 9 9107 9 9226 10 0913 maxd 2 0613e 005 maxd 1 1454e 004 maxd 1 3115e 004 k 6 V 10 5000 10 1557 10 0776 10 0439 9 9677 9 9107 9 9226 10 0913 PQ2 10 5000 0 0 10 0000 0 6000 0 4500 华南理工大学广州学院电气工程系课程报告 23 10 0000 0 4000 0 3000 10 0000 0 4000 0 2800 10 0000 0 6000 0 4000 10 0000 0 4000 0 3000 10 0000 0 5000 0 3500 10 0000 0 5000 0 4000 x 0 0024 x 0 0025 x 0 0038 x 0 0096 x 0 0340 x 0 0040 x 0 1761 maxd 3 2818e 008 V2 10 5000 10 1557 10 0776 10 0439 9 9677 9 9107 9 9226 10 0913 maxd 1 9630e 007 maxd 2 0684e 005 maxd 2 1077e 005 k 7 V 10 5000 10 1557 10 0776 10 0439 9 9677 9 9107 9 9226 10 0913 四 总结 4 14 1 三机九节点系统结果分析三机九节点系统结果分析 busnumber 9 aclinenumber 6 tansnumber 3 表示此次计算为 9 节点 6 支路 3 变压器的网络 B B2 分别为 B 和 B JP 为除平衡节点外的矩