《电力系统实验.docx》由会员分享，可在线阅读，更多相关《电力系统实验.docx（3页珍藏版）》请在三一办公上搜索。
1、实验一生成节点导纳矩阵的Inatlab程序设计一、实验目的：本实验通过对生成节点导纳矩阵的编程与调试，获得复杂电力系统的节点导纳矩阵，为电力系统的潮流计算程序编制打下根底。通过实验教学加深学生对节点导纳矩阵计算方法的理解，掌握变压器为非标准变比时的修正方法，熟悉MAT1.AB软件的使用方法，提高编制调试计算机程序的能力。二、实验器材：计算机、软件（己安装应用软件MAT1.AB等）、移动存储设备（学生自备，软盘、U盘等）。三、实验内容：1.节点导纳矩阵的形成节点导纳矩阵的对角元YH=1,2,./2）称自导纳，自导纳Yii是节点i以外的所有节点都接地，由节点i向整个网络看到而得到的导纳。在节点i加
2、上单位大小的电压（=1的电压）时，由节点i流向网络的电流就等于节点，的自导纳通过计算不难发现，就等于与节点i连接的所有支路导纳的和。节点导纳矩阵的非对角元%（，=1,2,./；J=i）称互导纳，即把节点j以外的节点全接地，而在节点,加以单位电压时，由节点i流向）的电流加上负号的值。同样，/其实就是连接节点，和节点/支路的导纳值再加上负号。显然，Yij=Yji.而且，如节点i,J之间没有直接联系，也不计两支路之间，例如两相邻电力线路之间的互感时，4=0。互导纳的这些性质决定了节点导纳矩阵是一个对称的稀疏矩阵。而且，由于每个节点所连接的支路数总有一定限度，随着网络中节点数的增加，非零元素数相对愈来
3、愈少，节点导纳矩阵的稀疏度，即零元素数与总元素数的比值也就愈来愈高。导纳矩阵的计算归结如下：（1）导纳矩阵的阶数等于电力系统网络的节点数。（2）导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数。（3）导纳矩阵的对角元素，即各节点的自导纳等于相应节点所连支路的导纳之和。式中，为节点i与节点/间支路阻抗马倒数，符号i表示/属于i或与i连接的即内只包括与节点i直接相连的节点人当节点，有接地支路时，还应包括/=0的情况。（4）导纳矩阵非对角元素%等于节点i与节点J之间的导纳的负数当人/之间有多条并联支路时，求修时应求所有并联支路导纳的代数和的负数式中，1.表示A）之间并联支路数的条
4、数。按照以上计算式，对于实际网络均可根据给定的支路参数和连接情况，直观而简单地求出导纳矩阵。可以看出，用以上计算公式求得的导纳矩阵与根据定义得到的导纳矩阵是完全一致的。2 .变压器为非标准变比时的修正无论采用有名制或标么制，凡涉及多电压级网络的计算，在精确计算时都必须将网络中所有参数和变量按市价变比归算到同一电压等级。实际上，在电力系统计算中总是有些变压器的实际变比不等于变压器两侧所选电压基准值之比，也就是不等于标准变比，而且变压器的变比在运行中是可以改变的。这将使每改变一次变比都要从新计算元件参数，很不方便。下面将介绍另一种可等值地表达变压器电压变换功能的模型。等值变压器模型图1非标准变比时
5、的修正电路图2以变压器阻抗表示增加非零非对角元素为节点i的自导纳，增加一个改变量为节点，的自导纳，也增加一个改变量为3 .程序代码及说明n=input(,请输入节点数:e);nl=inputC请输入支路数:nl=);%isb=inputC请输入平衡母线节点号:isb=);%pr=input(请输入误差精度:pr=);Bl=inputC请输入由各支路参数形成的矩阵:Bl=);%B2=inpu(请输入各节点参数形成的矩阵:B2=);X=inpuR请输入由节点号和接地支路参数形成的矩阵:X=%Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=ze
6、ros(1,n);S1=zeros(nl);%zeros(m,n)产生mXn的零矩阵,ZerOS(n)产生nXn的全0方阵。fori=l:n%节点数ifX(i,2)*=0;%matlab关系符号：(大于),=(大于等于)，(小于),=(小于等于)，=(等于)二(不等于)p=X(i,l);Y(p,p)=lX(i,2);%接地支路J.尸点除代表矩阵对应元素相除endendfori=l:nl%支路数181&6)=0%折算到哪一侧的标志,()非标准变比在4侧非标准变比在在P侧P=Bl(i,l);q=Bl(i,2);%Bl(i,l),B1(i,2)为支路编号elsep=Bl(i,2);q=Bl(i,l)
7、;endY(p,q)=Y(p,q)-l./(Bl(i,3)*Bl(i,5);%计算节点P和q间的互导纳,Bl(i,3)为支路阻抗,Bl(i,5)为变压器支路变比七Bl(i,5)=l时不是变压器支路Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1.(B1(i,3)*Bl(i,5)2)+B1(i,4)./2;%对角元,节点q的自导纳,Bl(i,4)为支路对地容抗Y(p,p)=Y(p,p)+1.Bl(i,3)+Bl(i,4)./2;%对角元,节点P的自导纳end%输出导纳矩阵dispC导纳矩阵Y=);disp(Y);四、实验数据：电力系统接线图如图3所示，根据图中数据可得Bl和X矩阵：图3
8、电力系统接线图Bl=l20.03i01.050;230.08+0.3i0.5i10;240.1+0.35i010;340.04+0.25i0.5i10;350.015i01.051;X=10;20;30;40;50;五、实验过程1 .手算图3所示网络的导纳矩阵。2 .认真阅读源程序，建立matlab的M程序，输入节点导纳矩阵计算程序。3 .输入实验数据,用Bl和X矩阵代替Bl和X输入语句并调试程序。Bl矩阵Bl(i,l)为支路始端编号p；81d2)为支路末端编号4，且pq；Bl(i,3)为支路阻抗；Bl(i,4)为支路对地容抗，无填0；Bl(i,5)为变压器支路变比，非变压器支路填1；Bl(i,6)为变压器变比折算到哪一侧的标志,0非标准变比在q侧,1非标准变比在在p侧，非变压器支路填0；X矩阵X(i,l)为节点编号；X(i,2)为该节点的接地支路阻抗，如果节点没有接地支路填“0”；4 .正确输入节点数和支路数进行计算，将计算结果与手算结果进行比拟。六、要求认真书写实验报告并上交实验数据编写:许琅，饶美丽2013-12
