极坐标法潮流计算.docVIP

课程设计（论文）任务书 专业 学生姓名 学 号 题名称 课程名称 潮流计算课程设计 课程编号 121202306 设计地点 综合仿真实验室 课程设计（论文）目的 1.掌握电力系统极坐标下的牛顿-拉夫逊计算的基本原理； 2.掌握并能熟练运用一门计算机语言（MATLAB语言或FORTRAN或C语言或C++语言）； 3.采用计算机语言对极坐标下的牛顿-拉夫逊计算进行计算机编程计算。 通过课程设计, 使学生巩固电力系统潮流计算的基本原理与方法，掌握潮流计算的数值求解方法(节点导纳矩阵，修正方程)，开发系统潮流计算的计算程序。让学生掌握用计算机仿真分析电力系统的方法。同时，通过软件开发，也有助于计算机操作能力和软件开发能力的提高。 已知技术参数和条件 在图所示的简单电力系统中，系统中节点1、2为节点，节点3为PV节点，节点4为平衡节点，已给定，，，，，，网络各元件参数的标幺值如表2所示，给定电压的初始值如表2所示，收敛系数。试求： 采用极坐标下的牛顿-拉夫逊计算图1网络的潮流分布。 任务和要求 任务： 熟练掌握计算机语言，并采用计算机编程进行下列计算： 根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；掌握潮流计算的数值求解方法(节点导纳矩阵，修正方程)，开发系统潮流计算的计算程序。 要求： 1.手工计算，手写，采用A4纸，得出计算结果。2.编写程序：它包括程序源代码；程序说明；部分程序的流程图；程序运行结果，电子版。 注：1此表由指导教师填写，经教研室审批指导教师学生生效； 2此表1式3份，学生、指导教师、教研室各1份。 [1] 何仰赞等.电力系统分析[M]. 武汉：华中理工大学出版社，2002.3 [2] 西安交通大学等.电力系统计算[M].北京：水利电力出版社，1993.12 五、进度安排 2010年12月20日：下达课程设计的计划书，任务书，设计题目及分组情况。 2010年12月21日-23日：学生完成潮流计算的手工计算。 2010年12月24日：讲述课程设计编程的思路、要求。举例：用MATLAB软件编写的部分程序。 2010年12月25日-30日：学生编写程序。 2011年1月1日-3日：上机调试程序，得出正确结果。 2011年1月4日-5日：整理课程设计报告。 2011年1月6日：学生答辩。 六、教研室审批意见 教研室主任（签字）： 年 月 日 七|、主管教学主任意见 主管主任（签字）： 年 月 日 八、备注 指导教师（签字）： 学生（签字）： 设计主题 题目一：在下图所示的简单电力系统中，系统中节点1、2为节点，节点3为PV节点，节点4为平衡节点，已给定，，，，，，网络各元件参数的标幺值如表1所示，给定电压的初始值如表2所示，收敛系数。试求： 采用极坐标下的牛顿-拉夫逊计算图示网络的潮流计算。 表1 网络各元件参数的标幺值 支路 电阻 电抗 输电线路 变压器变比k 1—2 0.02 0.06 0.01 — 1—3 0.01 0.03 0.01 — 2—3 0.03 0.07 — — 2—4 0.0 0.05 — 0.9625 3—4 0.02 0.05 — — 表2各节点电压（初值）标幺值参数 节点i 1 2 3 4 1.00+j0.0 1.0+j0.0 1.0+j0.0 1.05+j0.0 3 潮流计算流程图 本次课程设计采用极坐标下的牛顿-拉夫逊计算网络的潮流计算。其牛顿-拉夫逊潮流计算程序框图如下所示。 图3.1 极坐标下的牛顿-拉夫逊潮流计算程序框图 4 手工计算 插入手写的潮流计算过程 5 MATLAB程序设计 5.1 程序 %电力系统极坐标下的牛顿-拉夫逊法潮流计算 disp(电力系统极坐标下的牛顿-拉夫逊法潮流计算:); clear n=input(请输入结点数：n=); n1=input(请输入PV结点数：n1=); n2=input(请输入PQ结点数：n2=); isb=input(请输入平衡结点：isb=); pr=input(请输入精确度：pr=); K=input(请输入变比矩阵:K=); C=input(请输入支路阻抗矩阵：C=); y=input(请输入支路导纳矩阵：y=); U=input(请输入结点电压矩阵：U=); S=input(请输入各结点的功率：S=); Z=zeros(1,n