潮流计算matlab实例计算.docx

《潮流计算matlab实例计算.docx》由会员分享，可在线阅读，更多相关《潮流计算matlab实例计算.docx（10页珍藏版）》请在优知文库上搜索。

1、潮流计算（malab）实例计算潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料：学习相关软件（软件自选:本设计选择MatIab进行设计）.O.45-jO.I50.4j().05O.O8+jO240.0l+j0.0390drl0.006OQrl+80dO.O4+jO.I2-0.2+j0,2）0.6+j0.l系统接线困中节点1为平衡节点.节点2+心4、.5为PQ节点口2 .在给定的电力网络上画出等值电路图.3 .运用计算机进行潮流计算.4 .编写设计说明书.一、设计原理1 .牛顿-拉夫逊原理牛顿迭代法是取x之后,在这个根底上,找到比x更接近的方程的跟,一步一步迭代,从而找到更接近方程

2、根的近似跟.牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程fS）=0的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根.电力系统潮流计算,一般来说,各个母线所供负荷的功率是的,各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵,然后由节点导纳矩阵列写功率方程，由于功率方程里功率是的，电压的幅值和相角是未知的,这样潮流计算的问题就转化为求解非线性方程组的问题了.为了便于用迭代法解方程组,需要将上述功率方程改写成功率平衡方程,并对功率平衡方程求偏导,得出对应的雅可比矩阵,给未知节点赋电压初值,一般为额定电压,将初值带入功率平衡方程,得到功率不平衡量,这样由功率不平

3、衡量、雅可比矩阵、节点电压不潮流计算（matlab）实例计算平衡量（未知的）构成了误差方程,解误差方程,得到节点电压不平衡量,节点电压加上节点电压不平衡量构成新的节点电压初值,将新的初值带入原来的功率平衡方程,并重新形成雅可比矩阵,然后计算新的电压不平衡量,这样不断迭代,不断修正,一般迭代三到五次就能收敛.牛顿一拉夫逊迭代法的一般步骤：（I）形成各节点导纳矩阵Y.（2）设个节点电压的初始值U和相角初始值e还有迭代次数初值为0.（3）计算各个节点的功率不平衡量.（4）根据收敛条件判断是否满足,假设不满足那么向下进行.（5）计算雅可比矩阵中的各元素.（6）修正方程式个节点电压（7）利用新值自第步开

4、始进入下一次迭代,直至到达精度退出循环.（8）计算平衡节点输出功率和各线路功率2 .网络节点的优化D静态地按最少出线支路数编号这种方法由称为静态优化法.在编号以前.首先统计电力网络个节点的出线支路数,然后,按出线支路数有少到多的节点顺序编号.当由n个节点的出线支路相同时,那么可以按任意次序对这n个节点进行编号.这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也2）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的,在编号过程中认为固定不变的，事实上,在节点消去过程中，每消去一个节点以后,与该节点相连的各节点的出线支路数将发生变化（增加

5、,减少或保持不变）.因此,如果每消去一个节点后,立即修正尚未编号节点的出线支路数,然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果,动态按最少出线支路数编号方法的特点就是按出线最少原那么编号时考虑了消去过程中各节点出线支路数目的变动情况.3 .MATLAB编程应用MaHab是MatrixLaboratory的缩写,主要包括：一般数值分析,矩阵运算、数字信号处理、建模、系统限制、优化和图形显示等应用程序.由于使用MaUab编程运算与人进行科学计算的思路和表达方式完全一致,所以不像学习高级语言那样难于掌握,而且编程效率和计算效率极高,还可在计算机上直接输出结果和精美的图形拷贝,所以

6、它确实为一高效的科研助手.二、设计内容1.设计流程图启城/9输入原始数据/计算节点的新电压e（*+i）=e*）+/*）iiiUle-5fori=l:4forj=l:4H(i,j)=O;N(i,j)=O;M(i,j)=OL(izj)=0;oP(i)=0;oQ(i)=0;endendfori=l:4潮流计算(matlab)实例计算forj=l:5oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);0Q(i)=0Q(i)-U(i)*U(j)*(G(izj)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j

7、);end0P(i)=0P(i)+P(i);0Q(i)=0Q(i)+Q(i);endfx=0P,0Q,;为求雅克比矩阵3当1=：时候求H,N,M,L如下：fori=l:4forj=l:4ifi=jH(i,j)=-U(i)*Uj)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);N(i,j)=-Ui)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);1.(i,j)=H(i,j);M(i,j)=-N(i,j);endendendH,NzMzL彩当1=:时H,N,M,L如下：fori=l:4forj=l:5ifi=

8、jH(i,i)=H(i,i)+U(i)*U(j)*(G(irj)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);N(izi)=N(izi)-U(i)*U(j)*(G(izj)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);L(i,i)=L(izi)-U(i)*U(j)*(G(izj)*sin(e(i)-e(j)-B(izj)*cos(e(i)-e(j);endendN(i,i)=N(i,i)-2*(U(i)J2*G

9、(i,i)；1.(i,i)=L(i,i)+2*(U(i)J2*B(i,i)；endJ=HzN;MZL%J为雅克比矩阵ox=-(inv(J)*fx);fori=l:4oe(i)=ox(i);oU(i)=ox(i+4)*U(i);endfori=l:4e(i)=e(i)+oe(i);U(i)=U(i)+oU(i);endcount=count+l;endoxzUzezcount名求节点注入的净功率i=5;forj=l:5潮流计算(matlab)实例计算P(i)=U(i)*U(j)*(G(irj)*cos(e(i)-e(j)+B(irj)*sin(e(i)-e(j)+P(i);Q(i)=U(i)*U

10、(j)*(G(izj)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j)+Q(i);endS(5)=P(5)+Q(5)*sqrt(-l);S考求节点注入电流I=Y*U3.运行结果eY值：ComnandWindowi1N忙MATbCB?WatChthisWidEC.SeeCenncaorreadRettina3忸rtEd.10.8333-32,5000i-1.6667+5,000011.6667+5.uiiui-2.5000+7.SOOOi-1.6667+5.0()00112.9167-38.75001-10.0000+30.00001-L.25LIL1-1.6667+b.u

11、uui-iu.uuu+30.uuu0i12.9167-38.75001-1.25uu+3.7500i-2.5000+7.5(XX)I-1.2500+3.750013.7500-11.25004-5.uuuu+15.uuu0i-1.25uu4-3.750016.25D0迭代过程：count-+15.000013.7500i-I8.75001U05. OOOO5.005. OOOO30. OOOO7. 50005, OOOO30.003. 75007. 50003. 7500N1.66671-56672. 5000L 66671.666710.00001.25001QOOOO1.2500-1. 6

12、6672.5000Q-1.6667T.6QS7-2.5000-10.OOOO-1.6667-10.0000O1.25-1. 2500-2.505.OQOO5.IJUUU7.50005.00003u.IJUUU5.000030.003.75007.50a3.75000J三-33.10006.00005.uOOj7.500j-10.53331.66776672.600u5.00-33.975030.00a1.666712.S41710.co5.OOOO30.uuOO-38.75005.75001.666710.COOQ-12.91671.25007.5000Q3.7500-11.250?2.50001.2500-3.750011.1333-1.6667-1.6667-2.50uJ-31.60005.iJUUiJ5.JUUJ7.5Uuu-1.6667-1.666712.9917-10.Duuu