电力系统分析实验指导书.doc

PAGE 35 电力系统分析 实验指导书 郑州科技学院电气工程及其自动化教研室编 目 录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc365450297 实验一 生成节点导纳矩阵的matlab程序设计 PAGEREF _Toc365450297 \h 3 HYPERLINK \l _Toc365450298 实验二 支路追加法生成节点阻抗矩阵的matlab程序设计 PAGEREF _Toc365450298 \h 7 HYPERLINK \l _Toc365450299 实验三 牛顿-拉夫逊法电力系统潮流计算的matlab程序设计 PAGEREF _Toc365450299 \h 15 HYPERLINK \l _Toc365450300 实验四 简化模型的静态稳定计算的matlab程序设计 PAGEREF _Toc365450300 \h 26 HYPERLINK \l _Toc365450301 实验五 短路电流计算程序的实现 PAGEREF _Toc365450301 \h 33 PAGE 21 实验一 生成节点导纳矩阵的matlab程序设计 一、实验目的： 本实验通过对生成节点导纳矩阵的编程与调试，获得复杂电力系统的节点导纳矩阵，为电力系统的潮流计算程序编制打下基础。通过实验教学加深学生对节点导纳矩阵计算方法的理解，掌握变压器为非标准变比时的修正方法，熟悉MATLAB软件的使用方法，提高编制调试计算机程序的能力。 二、实验器材： 计算机、软件（已安装应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等）。 三、实验内容： 1.节点导纳矩阵的形成 节点导纳矩阵的对角元称自导纳，自导纳是节点以外的所有节点都接地，由节点向整个网络看到而得到的导纳。在节点加上单位大小的电压（的电压）时，由节点流向网络的电流就等于节点的自导纳 通过计算不难发现，就等于与节点连接的所有支路导纳的和。 节点导纳矩阵的非对角元称互导纳，即把节点j以外的节点全接地，而在节点j加以单位电压时，由节点流向j的电流加上负号的值。 同样，其实就是连接节点i和节点j支路的导纳值再加上负号。显然，。而且，如节点i，j之间没有直接联系，也不计两支路之间，例如两相邻电力线路之间的互感时，。 互导纳的这些性质决定了节点导纳矩阵是一个对称的稀疏矩阵。而且，由于每个节点所连接的支路数总有一定限度，随着网络中节点数的增加，非零元素数相对愈来愈少，节点导纳矩阵的稀疏度，即零元素数与总元素数的比值也就愈来愈高。 导纳矩阵的计算归结如下： 导纳矩阵的阶数等于电力系统网络的节点数。 导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数。 导纳矩阵的对角元素，即各节点的自导纳等于相应节点所连支路的导纳之和。 式中，为节点i与节点j间支路阻抗倒数，符号表示属于或与连接的，即内只包括与节点i直接相连的节点j。当节点i有接地支路时，还应包括的情况。 (4)导纳矩阵非对角元素等于节点i与节点j之间的导纳的负数 当i、j之间有多条并联支路时，求时应求所有并联支路导纳的代数和的负数 式中，L表示i、j之间并联支路数的条数。 按照以上计算式，对于实际网络均可根据给定的支路参数和连接情况，直观而简单地求出导纳矩阵。可以看出，用以上计算公式求得的导纳矩阵与根据定义得到的导纳矩阵是完全一致的。 2.变压器为非标准变比时的修正 无论采用有名制或标么制，凡涉及多电压级网络的计算，在精确计算时都必须将网络中所有参数和变量按实际变比归算到同一电压等级。实际上，在电力系统计算中总是有些变压器的实际变比不等于变压器两侧所选电压基准值之比，也就是不等于标准变比，而且变压器的变比在运行中是可以改变的。这将使每改变一次变比都要重新计算元件参数，很不方便。下面将介绍另一种可等值地体现变压器电压变换功能的模型。 等值变压器模型 图1 非标准变比时的修正电路 图2 以变压器阻抗表示 增加非零非对角元素为 节点i的自导纳，增加一个改变量为 节点j的自导纳，也增加一个改变量为 3. 程序代码及说明 n=input(请输入节点数:n=);nl=input(请输入支路数:nl=);%isb=input(请输入平衡母线节点号:isb=);%pr=input(请输入误差精度:pr=);B1=input(请输入由各支路参数形成的矩阵:B1=);%B2=input(请输入各节点参数形成的矩阵:B2=);X=input(请输入由节点号和接地支路参数形成的矩阵:X=); Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n); O=