电力系统分析潮流计算课程序设计及其MATLAB程序设计

学号 班内序号 4推后代潮流计算的原理 5后代潮流计算迭代过程 5 五、确定前推回代支路次序 7六、前推回代计算输入文件 8 七、配电网潮流计算的要点 11 (2)在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方1945K(Ω)R节点编号节点类型1根节点5叶节点7叶节点节点初始电(kV)P(MVA)i0Q(MVA)i099叶节点10叶节点12叶节点13叶节点看懂前推回代法计算程序；2.报告叙述计算原理及计算流程；3.绘制计算流程框图；4.确定前推回代支路次序(广度优先，或深度优先)并编写前推回代计算输入文5.整理专利计算结果；6.总结配电网潮流计算的要点；2.赋予各节点电压变量(直角坐标系形式)初值后，求解不平衡量；5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和节点6.上机编程调试；连调；果作比较分析。8.准备计算机演示答辩，书写该课程设计说明书(必须计算机打印)方程牛顿迭代法(Newton'smethod)又称为牛顿-拉夫逊(拉弗森)方法(Newton-Raphsonmethod世纪提出的一种在实数域和复数域上近似求解方程的方法。化的一种近似方法。把f(x)在x0点附近展开成泰勒级数f(x)=f(x0)+(x－f(x(n))/f'(x(n。))值，则系统母线电压幅值的微小变化U对母线有功功率的改变P影响很小。同样，功率方程在用极坐标形式表示时，它的修正方程式可简化为的变化而且又都是不对称矩阵。对牛顿法的进一步简化，即把式(2-1)中的系数矩阵9在一般情况下，线路两端的电压相角ij是不大的(不超过10°~20°)，因此，LDiU2iiiiiii ijijijL=UUBijijij「UBU|1111...|H|...|n-1n-1,1n-1n-1,111122UBU2222UBUn-1n-1,22...UBU]11,n-1n-1|...U2B2,n-1Un-1||...Un-1Bn-1,n-1Un-1」=|=|||LU2||11||B21...||...BB1,n=1||1||......||||U=UB'UD1D1「UBU|1111|LUBUmm11UBU...UBU]UBU...UBU]||||U」n-1||||||.||=1U2B...B]「UB12...B1B...BU2]||||U」=UB"UD2D2(2-5)将式(2-4)和(2-5)代人式|lU.i=Uiej9i=Ui(cos9i+jsin9i)中，得到D1D1D2D1D1和D1D111BBBB…Bn-1,2B……… iisijijijijj=1iisijijijijj=1(2-10)式(2-8)、(2-9)和(2-10)构成了PQ分解法迭代过程的基本方程式。理节省了大量的计算机内存，对于任何种类的配电网只要有合理的R/X值，此方法均可代过程NBNB-1Pi)=∑PL(j)+∑LP(j)Qi=∑QL(j)+∑LQ(j)即认为迭代收敛，则迭代过程结束。在开始计算的时候先把各节点上的LP(i)、LQ(i)分别设为零，利用式(1)计算出是否小于设定阈值判断是否收敛，如果满足，则输出结果，如果不满足,必须重复上述。第二步，利用第一步所得的支路首端功率和本步骤刚算出的本支路始节点的电压(对电源点为已知电压)，从电源点开始逐条支路进行计算，求得各支路终节点的电压。第二种是逐条追加支路的方法。首先从根节点(电源点)开始接出第一条支路，引出一个新节点，以后每次追加的支路都必须从已出现的节点接出遵循这个原则逐条追加支输入原始数据形成导纳矩阵Y值e0，f0置K＝0雅可比矩阵各元素是计算线路功率及平衡节点功率输出修正各节点电压kekfkfkK=K+1[PQ,FT,RX]=case115();%调用数据文件NN=size(PQ,1);%节点数NB=size(FT,1);%支路数数V=PQ(:,1);%V初始电压相量maxd=1k0;whilemaxd>PQ2=PQ;%每一次迭代各节点的注入有功和无功相同PL=;fori=1:NBkf=FT(i,1);kt=FT(i,2);%前推始节点号%前推终节点号x=(PQ2(kf,2)^2+PQ2(kf,3)^2)/(V(kf)^2);%计算沿线电流平方APQ1(i,1)=RX(i,1)*x+PQ2(kf,2);%计算支路首端有功/MWRX(i,1)~RPQ1(i,2)=RX(i,2)*x+PQ2(kf,3);%计算沿支路的无功损耗/MvarRX(i,2)~XPL=PL+RX(i,1)*x;%累积计算整个网络的功率损耗angle(1)=;fori=NB:-1:1kf=FT(i,2);kt=FT(i,1);%回代始节点号%回代终节点号00dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2))/V(kf);dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1))/V(kf);V2(kt)=sqrt((V(kf)-dv1)^2+(dv2)^2);angle(kt)=atand(dv2/(V(kf)-dv1))+angle(kf);maxd=abs(V2(2)-V(2));V2(1)=V(1);fori=3:1:NNifabs(V2(i)-V(i))>maxd;maxd=abs(V2(i)-V(i));maxdk=k+1PQ1V=V2anglePQ1V=V2anglePLend%节点电压角度计算结果单位度function[PQ,FT,RX]=case115()PQ=[%节点电压有功无功dv1dv2%计算支路末端电压/kV%计算支路电压相角(单位度)21FT=[%首端末端5411437698326282XXRX=[%Rk=8PQ1=V=angle=0PL=1.需要注意系统中所有母线的电压是否在允许的范围以内。2.系统中各种元件(线路、变压器等)是否会出现过负荷。3.在进行电压和功率分布计算以前，先要对网络等值电路作些简化处理。4.确立前推回代的计算次序。乏不[1]何仰赞温增银.电力系统分析.华中科技大学出版社，2002