电力系统实验南昌大学讲义.doc

实 验 报 告 实验课程： 电力系统分析 学生姓名： 李瑞欣 学 号： 6101113078 专业班级： 电气工程及其自动化132 指导老师： 徐敏 2015年 12月 日 南昌大学实验报告 学生姓名： 李瑞欣 学 号： 610113078 专业班级： 电气132 实验类型：□ 验证 □ 综合 □ 设计 □ 创新 实验日期： 实验成绩： 一、实验项目名称 电力网数学模型模拟实验 二、实验目的与要求： 本实验通过对电力网数学模型形成的计算机程序的编制与调试，获得形成电力网数学模型：节点导纳矩阵的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩阵。通过实验教学加深学生对电力网数学模型概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。 三、主要仪器设备及耗材 计算机、软件（已安装，包括各类编程软件C语言、C++、VB、VC等、应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等） 四、实验步骤 将事先编制好的形成电力网数学模型的计算程序原代码由自备移动存储设备导入计算机。 在相应的编程环境下对程序进行组织调试。 应用计算例题验证程序的计算效果。 对调试正确的计算程序进行存储、打印。 完成本次实验的实验报告。 五、实验数据及处理结果 运行自行设计的程序，把结果与手工计算结果相比较，验证所采用方法及所编制程序运行的正确性。 实验数据 见《电力系统分析》（上册）72页 例4-1 程序调试记录，逐条记录下程序调试过程中发现的问题及解决的方法。 1、程序调试过程出了几次错，例如在做复数矩阵的转置使错误使用了“ ’ ”命令(共轭转置，应改正为“ ’ ”(普通转置)。在编程时尽量不使用 上图演示了 2、考虑变压器支路的影响，修改原矩阵时要先修改Y（节点对地的导纳矩阵）再修改Z（节点之间的阻抗矩阵）否则会出错。因为在本程序中在修改Y时要用到原来的Z。 3、程序的设计方法是： 首先不考虑变压器支路的影响，手动化简变压器参数，然后编写程序，计算节点导纳矩阵，成功后再考虑变压器影响，修改源程序，最后把计算节点导纳矩阵部分做成函数文件，以便后续的实验调用。 本程序的算法见程序框图。在数据结构方面，我采用的是用三个矩阵分别输入网络中阻抗、导纳和变压器情况。 程序框图。 所编程序清单。 函数文件：JDDN.m function Yout = JDDN( Z,Y,T,N )% 计算节点导纳矩阵的函数% 有三个输入参数% 1.矩阵Z. 它代表了网络的节点间互阻抗矩阵，Y（1，2）表示1，2节点之间阻抗% 2.矩阵Y. 它代表了网络中各节点对地导纳，Y(1,2)表示1，2支路等效电路中1节点的对地导纳% 3.矩阵K. 它是有关变压器的信息,正k是变压器原边，-k是变压器副边，0是非变压器支路% 4.N表示总节点数，也是矩阵的阶数% 输出参数即为网络的节点导纳矩阵% 检查变压器支路，修改矩阵for m=1:N for n=1:N if T(m,n)0 k=T(m,n);%变比 Y(m,n)=Y(m,n)+1/( k*Z(m,n)/(k-1) ); Z(m,n)=Z(m,n)*k; elseif T(m,n)0 k=abs(T(m,n)); Y(m,n)=Y(m,n)+1/( k^2*Z(m,n)/(1-k) ); Z(m,n)=Z(m,n)*k; end endend% 计算节点导纳矩阵Yout=zeros(N);Zd=1./Z; %把网络互阻抗矩阵所有非零元素做成倒数以便求计算Zd(Zd==inf)=0;for m=1:N for n=1:N if m==n Yout(m,