还剩37页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
课程设计任务书课程名称电力系统稳态分析课程设计设计题目电力系统的潮流分布计算—牛拉法指导教师时间1周
一、教学要求电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点,提高学生综合能力为目标,尽可能结合实际工程进行设计内容的安排要充分考虑学校现有的设备,设计时间及工程实际需要,并使学生初步学会运用所学知识解决工程中的实际问题
二、设计资料及参数
(一)设计原始资料
1、待设计电气设备系统图
2、电力系统网络各元件参数
3、电力系统电气元件的使用规范
4、电力工程电气设计手册
(二)设计参考资料
1、《电力系统稳态分析》,陈珩,中国电力出版社,2007,第三版
2、《电力系统分析》,韩祯祥,浙江大学出版社,2005,第三版
3、《电力系统分析课程实际设计与综合实验》,祝书萍,中国电力出版社,2007,第一版
三、设计要求及成果
1.根据给定的参数或工程具体要求,收集和查阅资料;学习相关软件(软件自选)
2.在给定的电力网络上画出等值电路图
3.运用计算机进行潮流计算
4.编写设计说明书基本要求
1.编写潮流计算程序;
2.在计算机上调试通过(?);
3.运行程序并计算出正确结果(?);
4.写出课程设计报告(包括以下内容)(1份)
(1)程序框图;
(2)源程序;
(3)符号说明表;
(4)算例及计算结果
5.编写计算说明书(1份)
四、进度安排根据给定的参数或工程具体要求,收集和查阅资料(半天)学习软件(MATLAB或C语言等)(一天半)编程计算复杂系统潮流计算(三天)编写计算设计书(一天)
五、评分标准课程设计成绩评定依据包括以下几点工作态度(占10%);基本技能的掌握程度(占20%);程序编写是否合理是否有运行结果40%;课程设计说明书编写水平占30%分为优、良、中、合格、不合格五个等级考核方式设计期间教师现场检查;评阅设计报告系统接线图其中节点1为平衡节点,节点
2、
3、
4、5为PQ节点前言一课程设计任务书…………………………………………………………2第1章电力系统潮流分布简述…………………………………………………6第1节潮流计算的发展趋势……………………………………6第2节潮流计算的发展史………………………………………6第3节潮流计算的意义…………………………………………8第2章设计程序……………………………………………………………………8第1节潮流计算的基本步骤……………………………………………8第2节潮流计算设计资料及参数……………………………………10第3节P-Q法潮流计算流程图…………………………………………10第4节潮流计算的matlab程序…………………………………………11第5节潮流计算程序运行结果显示……………………………………16附录感想……………………………………………………………………………39参考……………………………………………………………………………39第一章电力系统潮流计算简述电力系统的潮流分布计算(牛顿——拉夫逊直角坐标法)第一节潮流计算的意义 1在电网规划阶段通过潮流计算合理规划电源容量及接入点合理规划网架选择无功补偿方案满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求 2在编制年运行方式时在预计负荷增长及新设备投运基础上选择典型方式进行潮流计算发现电网中薄弱环节供调度员日常调度控制参考并对规划、基建部门提出改进网架结构加快基建进度的建议 3正常检修及特殊运行方式http://baike.baidu.com/view/
1269321.htm\t_blank下的潮流计算用于日运行方式的编制指导发电厂开机方式有功、无功调整方案及负荷调整方案满足线路、变压器热稳定要求及电压质量要求 4预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案 总结为在电力系统运行方式http://baike.baidu.com/view/
3087326.htm\t_blank和规划方案的研究中,都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性同时,为了实时监控http://baike.baidu.com/view/
647394.htm\t_blank电力系统的运行状态,也需要进行大量而快速的潮流计算因此,潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算在系统规划设计和安排系统的运行方式时,采用离线潮流计算;在电力系统运行状态http://baike.baidu.com/view/
2525114.htm\t_blank的实时监控中,则采用在线潮流计算第二节潮流计算的发展史 利用电子计算机进行潮流计算从20世纪50年代中期就已经开始此后,潮流计算曾采用了各种不同的方法,这些方法的发展主要是围绕着对潮流计算的一些基本要求进行的对潮流计算的要求可以归纳为下面几点
(1)算法的可靠性或收敛性
(2)计算速度和内存占用量
(3)计算的方便性和灵活性 电力系统潮流计算属于稳态分析范畴,不涉及系统元件的动态特性和过渡过程因此其数学模型不包含微分方程,是一组高阶非线性方程非线性代数方程组的解法离不开迭代,因此,潮流计算方法首先要求它是能可靠的收敛,并给出正确答案随着电力系统规模的不断扩大,潮流问题的方程式阶数越来越高,目前已达到几千阶甚至上万阶,对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法 在用数字计算机求解电力系统潮流问题的开始阶段,人们普遍采用以节点导纳矩阵http://baike.baidu.com/view/
10337.htm\t_blank为基础的高斯-赛德尔迭代法这个方法的原理比较简单,要求的数字计算机的内存量也比较小,适应当时的电子数字计算机http://baike.baidu.com/view/
907732.htm\t_blank制作水平和电力系统理论水平,于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法 20世纪60年代初,数字计算机已经发展到第二代,计算机的内存和计算速度发生了很大的飞跃,从而为阻抗法的采用创造了条件阻抗矩阵是满矩阵,阻抗法要求计算机储存表征系统接线和参数的阻抗矩阵这就需要较大的内存量而且阻抗法每迭代一次都要求顺次取阻抗矩阵中的每一个元素进行计算,因此,每次迭代的计算量很大 阻抗法改善了电力系统潮流计算问题的收敛性,解决了导纳法无法解决的一些系统的潮流计算,在当时获得了广泛的应用,曾为我国电力系统设计、运行和研究作出了很大的贡献但是,阻抗法的主要缺点就是占用计算机的内存很大,每次迭代的计算量很大当系统不断扩大时,这些缺点就更加突出为了克服阻抗法在内存和速度方面的缺点,后来发展了以阻抗矩阵为基础的分块阻抗法这个方法把一个大系统分割为几个小的地区系统,在计算机内只需存储各个地区系统的阻抗矩阵及它们之间的联络线的阻抗,这样不仅大幅度的节省了内存容量,同时也提高了计算速度 克服阻抗法缺点的另一途径是采用牛顿-拉夫逊法(以下简称牛顿法http://baike.baidu.com/view/
1504001.htm\t_blank)牛顿法是数学中求解非线性方程式的典型方法,有较好的收敛性解决电力系统潮流计算问题是以导纳矩阵为基础的,因此,只要在迭代过程中尽可能保持方程式系数矩阵http://baike.baidu.com/view/
3503330.htm\t_blank的稀疏性,就可以大大提高牛顿潮流程序的计算效率自从20世纪60年代中期采用了最佳顺序消去法以后,牛顿法在收敛性、内存要求、计算速度方面都超过了阻抗法,成为直到目前仍被广泛采用的方法 在牛顿法的基础上,根据电力系统的特点,抓住主要矛盾,对纯数学的牛顿法进行了改造,得到了P-Q分解法P-Q分解法在计算速度方面有显著的提高,迅速得到了推广 牛顿法的特点是将非线性方程线性化20世纪70年代后期,有人提出采用更精确的模型,即将泰勒级数的高阶项也包括进来,希望以此提高算法的性能,这便产生了保留非线性的潮流算法另外,为了解决病态潮流计算,出现了将潮流计算表示为一个无约束非线性规划问题的模型,即非线性规划潮流算法 近20多年来,潮流算法的研究仍然非常活跃,但是大多数研究都是围绕改进牛顿法和P-Q分解法进行的此外,随着人工智能理论的发展,遗传算法、人工神经网络http://baike.baidu.com/view/
19743.htm\t_blank、模糊算法也逐渐被引入潮流计算但是,到目前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位由于电力系统规模的不断扩大,对计算速度的要求不断提高,计算机的并行计算技术也将在潮流计算中得到广泛的应用,成为重要的研究领域第三节潮流计算的发展趋势 通过几十年的发展,潮流算法日趋成熟近几年,对潮流算法的研究仍然是如何改善传统的潮流算法,即高斯-塞德尔法、牛顿法和快速解耦法牛顿法,由于其在求解非线性潮流方程时采用的是逐次线性化的方法,为了进一步提高算法的收敛性和计算速度,人们考虑采用将泰勒级数的高阶项或非线性项也考虑进来,于是产生了二阶潮流算法后来又提出了根据直角坐标形式的潮流方程是一个二次代数方程的特点,提出了采用直角坐标的保留非线性快速潮流算法对于一些病态系统,应用非线性潮流计算方法往往会造成计算过程的振荡或者不收敛,从数学上讲,非线性的潮流计算方程组本来就是无解的这样,人们提出来了将潮流方程构造成一个函数,求此函数的最小值问题,称之为非线性规划潮流的计算方法优点是原理上保证了计算过程永远不会发散如果将数学规划原理和牛顿潮流算法有机结合一起就是最优乘子法另外,为了优化系统的运行,从所有以上的可行潮流解中挑选出满足一定指标要求的一个最佳方案就是最优潮流问题最优潮流是一种同时考虑经济性和安全性的电力网络http://baike.baidu.com/view/
429324.htm\t_blank分析优化问题OPF在电力系统的安全运行、经济调度、可靠性分析、能量管理以及电力定价等方面得到了广泛的应用第2章程序设计第一节潮流计算的基本步骤
1.形成节点导纳矩阵YB由上图可得相应的节点导纳矩阵Y=[
6.25-
18.75i-5+15i-
1.25+
3.75i00;-5+15i
10.834-
32.5i-
1.667+5i-
1.667+5i-
2.5+
7.5i;-
1.25+
3.75i-
1.667+5i
12.917-
38.75i-10+30i0;0-
1.667+5i-10+30i
12.917-
38.75i-
1.25+
3.75i;0-
2.5+
7.5i0-
1.25+
3.75i
3.75-
11.25i]由于节点1为平衡节点,其他的节点均为PQ节点,系数矩阵,阶数相同又应对该等值网络,不存在除去与有功功率和电压相位或无功功率和电压大小关系较小的因素的可能性,这两个矩阵,完全相同他们就由导纳矩阵的虚部部分中除第一行第一列外的各个元素所组成.2计算各节点有功功率不平衡量取=
1.06,=0;====
1.0;====0计算各节点有功功率不平衡量=-各节点电压的相位角弧度由矩阵方程式-=3计算各节点无功功率不平衡量按下式计算各节点无功功率不平衡量=-由矩阵方程式计算各节点电压大小-=4计算平衡节点功率和线路功率==第二节潮流计算设计资料及参数题目
(一)设计原始资料网络接线如图所示,各支路导纳均以标幺值标于图中节点注入功率分别为,,,,其中节点1连接的实际上相当于给定功率的发电厂设节点5电压保持定值,试运用以极坐标表示的高斯—塞德尔法、牛顿—拉夫逊法、P-Q分解法计算该系统的潮流分布计算精度要求各节点电压修正量不大于第三节P-Q分解法的程序流程图第四节P-Q计算设计题目潮流分布的完整的MATLAB程序clccleardisp节点总数为;N=5disp平衡节点为;1dispPQ节点为;JD=
[2345]Y=[
6.25-
18.75i-5+15i-
1.25+
3.75i00;-5+15i
10.834-
32.5i-
1.667+5i-
1.667+5i-
2.5+
7.5i;-
1.25+
3.75i-
1.667+5i
12.917-
38.75i-10+30i0;0-
1.667+5i-10+30i
12.917-
38.75i-
1.25+
3.75i;0-
2.5+
7.5i0-
1.25+
3.75i
3.75-
11.25i]G=realY;B=imagY;B1=B2:52:5;B2=B1;b1=invB1;b2=invB2;dt2=0;dt3=0;dt4=0;dt5=0;u2=
1.0;u3=
1.0;u4=
1.0;u5=
1.0;p2=
0.20;p3=-
0.45;p4=-
0.40;p5=-
0.60;q2=
0.20;q3=-
0.15;q4=-
0.05;q5=-
0.10;k=0;wucha=1;N1=4;whilewucha
10.00001;dt1=0;u1=
1.06;form=2:5forn=1:5ptn=um*un*Gmn*cosdtm-dtn+Bmn*sindtm-dtn;enddisp各节点有功率不平衡量为dpm=pm-sumptendform=2:5u1=diagu;u2=invu12:52:5;u3=flipudrot90u2:5;u4=u31:4:;dp1=flipudrot90dp2:5;dp2=dp11:4:;dtt=-b1*dp2/u4*u2;dtt=dtt:1:1;dtt=flipudrot90dttdisp各节点电压的相位角dt为dtm=dtm+dttm-1endform=2:5forn=1:5qtn=um*un*Gmn*sindtm-dtn-Bmn*cosdtm-dtn;enddisp各节点无功率不平衡量为qqm=qm-sumqtendform=2:5qq1=flipudrot90qq2:5;qq2=qq11:4:;ut=-b2*qq2/u4;disp各节点电压的修正量utdisp各节点电压的大小为um=um+utm-1endforn=1:5Un=un*cosdtn+j*sindtn;endforn=1:5In=Y1n*Un;enddisp各节点功率s为S1=U1*sumconjIform=1:5forn=1:5Smn=Um*conjUm-conjUn*conj-Ymn;endenddisp各节点功率sij为Swucha=maxabsutk=k+1End第五节程序运行结果显示节点总数为N=5平衡节点为ans=1PQ节点为JD=2345Y=
6.2500-
18.7500i-
5.0000+
15.0000i-
1.2500+
3.7500i00-
5.0000+
15.0000i
10.8340-
32.5000i-
1.6670+
5.0000i-
1.6670+
5.0000i-
2.5000+
7.5000i-
1.2500+
3.7500i-
1.6670+
5.0000i
12.9170-
38.7500i-
10.0000+
30.0000i00-
1.6670+
5.0000i-
10.0000+
30.0000i
12.9170-
38.7500i-
1.2500+
3.7500i0-
2.5000+
7.5000i0-
1.2500+
3.7500i
3.7500-
11.2500i各节点有功率不平衡量为dp=
00.5000各节点有功率不平衡量为dp=
00.5000-
0.3750各节点有功率不平衡量为dp=
00.5000-
0.3750-
0.4000各节点有功率不平衡量为dp=
00.5000-
0.3750-
0.4000-
0.6000dtt=-
0.0370-
0.0855-
0.0918-
0.1086各节点电压的相位角dt为dt=0-
0.0370000dtt=-
0.0370-
0.0855-
0.0918-
0.1086各节点电压的相位角dt为dt=0-
0.0370-
0.085500dtt=-
0.0370-
0.0855-
0.0918-
0.1086各节点电压的相位角dt为dt=0-
0.0370-
0.0855-
0.09180dtt=-
0.0370-
0.0855-
0.0918-
0.1086各节点电压的相位角dt为dt=0-
0.0370-
0.0855-
0.0918-
0.1086各节点无功率不平衡量为qq=
01.2119各节点无功率不平衡量为qq=
01.2119-
0.0773各节点无功率不平衡量为qq=
01.2119-
0.0773-
0.1919各节点无功率不平衡量为qq=
01.2119-
0.0773-
0.1919-
0.3196各节点电压的修正量ut=
0.
03950000.
00810000.
00640000.0001000各节点电压的大小为u=
1.
06001.
03951.
00001.
00001.0000各节点电压的修正量ut=
0.
03950000.
00810000.
00640000.0001000各节点电压的大小为u=
1.
06001.
03951.
00811.
00001.0000各节点电压的修正量ut=
0.
03950000.
00810000.
00640000.0001000各节点电压的大小为u=
1.
06001.
03951.
00811.
00641.0000各节点电压的修正量ut=
0.
03950000.
00810000.
00640000.0001000各节点电压的大小为u=
1.
06001.
03951.
00811.
00641.0001各节点功率s为S1=
1.1389+
0.2403i各节点功率sij为S=
00.7229+
0.1332i
0.4160+
0.1071i00-
0.7133-
0.1044i
00.3110+
0.0850i
0.3469+
0.0845i
0.6671+
0.1416i-
0.4029-
0.0676i-
0.3052-
0.0677i
00.2083-
0.0128i00-
0.3398-
0.0633i-
0.2078+
0.0141i
00.0714+
0.0033i0-
0.6499-
0.0900i0-
0.0710-
0.0021i0k=1各节点有功率不平衡量为dp=0-
0.4117-
0.3750-
0.4000-
0.6000各节点有功率不平衡量为dp=0-
0.
41170.0498-
0.4000-
0.6000各节点有功率不平衡量为dp=0-
0.
41170.
04980.0762-
0.6000各节点有功率不平衡量为dp=0-
0.
41170.
04980.
07620.1209dtt=-
0.
01050.
00120.
00170.0035各节点电压的相位角dt为dt=0-
0.0474-
0.0855-
0.0918-
0.1086dtt=-
0.
01050.
00120.
00170.0035各节点电压的相位角dt为dt=0-
0.0474-
0.0843-
0.0918-
0.1086dtt=-
0.
01050.
00120.
00170.0035各节点电压的相位角dt为dt=0-
0.0474-
0.0843-
0.0901-
0.1086dtt=-
0.
01050.
00120.
00170.0035各节点电压的相位角dt为dt=0-
0.0474-
0.0843-
0.0901-
0.1050各节点无功率不平衡量为qq=0-
0.1367-
0.0773-
0.1919-
0.3196各节点无功率不平衡量为qq=0-
0.
13670.0175-
0.1919-
0.3196各节点无功率不平衡量为qq=0-
0.
13670.
01750.0237-
0.3196各节点无功率不平衡量为qq=0-
0.
13670.
01750.
02370.0378各节点电压的修正量ut=-
0.
00380000.
00020000.
00040000.0009000各节点电压的大小为u=
1.
06001.
03581.
00811.
00641.0001各节点电压的修正量ut=-
0.
00380000.
00020000.
00040000.0009000各节点电压的大小为u=
1.
06001.
03581.
00831.
00641.0001各节点电压的修正量ut=-
0.
00380000.
00020000.
00040000.0009000各节点电压的大小为u=
1.
06001.
03581.
00831.
00681.0001各节点电压的修正量ut=-
0.
00380000.
00020000.
00040000.0009000各节点电压的大小为u=
1.
06001.
03581.
00831.
00681.0009各节点功率s为S1=
1.3260+
0.2508i各节点功率sij为S=
00.9151+
0.1435i
0.4109+
0.1072i00-
0.8998-
0.0977i
00.2415+
0.0816i
0.2741+
0.0809i
0.5423+
0.1344i-
0.3981-
0.0687i-
0.2378-
0.0707i
00.1903-
0.0114i00-
0.2695-
0.0672i-
0.1899+
0.0125i
00.0639+
0.0036i0-
0.5306-
0.0995i0-
0.0636-
0.0026i0k=2各节点有功率不平衡量为dp=
00.
04190.
04980.
07620.1209各节点有功率不平衡量为dp=
00.0419-
0.
00440.
07620.1209各节点有功率不平衡量为dp=
00.0419-
0.0044-
0.
00440.1209各节点有功率不平衡量为dp=
00.0419-
0.0044-
0.0044-
0.0058dtt=
0.
00160.
00060.
00060.0008各节点电压的相位角dt为dt=0-
0.0459-
0.0843-
0.0901-
0.1050dtt=
0.
00160.
00060.
00060.0008各节点电压的相位角dt为dt=0-
0.0459-
0.0838-
0.0901-
0.1050dtt=
0.
00160.
00060.
00060.0008各节点电压的相位角dt为dt=0-
0.0459-
0.0838-
0.0895-
0.1050dtt=
0.
00160.
00060.
00060.0008各节点电压的相位角dt为dt=0-
0.0459-
0.0838-
0.0895-
0.1043各节点无功率不平衡量为qq=
00.
01530.
01750.
02370.0378各节点无功率不平衡量为qq=
00.0153-
0.
00050.
02370.0378各节点无功率不平衡量为qq=
00.0153-
0.0005-
0.
00060.0378各节点无功率不平衡量为qq=
00.0153-
0.0005-
0.0006-
0.0001各节点电压的修正量ut=
1.0e-003*
0.
77760000.
51180000.
54910000.6921000各节点电压的大小为u=
1.
06001.
03661.
00831.
00681.0009各节点电压的修正量ut=
1.0e-003*
0.
77760000.
51180000.
54910000.6921000各节点电压的大小为u=
1.
06001.
03661.
00881.
00681.0009各节点电压的修正量ut=
1.0e-003*
0.
77760000.
51180000.
54910000.6921000各节点电压的大小为u=
1.
06001.
03661.
00881.
00731.0009各节点电压的修正量ut=
1.0e-003*
0.
77760000.
51180000.
54910000.6921000各节点电压的大小为u=
1.
06001.
03661.
00881.
00731.0016各节点功率s为S1=
1.2936+
0.2441i各节点功率sij为S=
00.8855+
0.1383i
0.4081+
0.1057i00-
0.8712-
0.0955i
00.2475+
0.0815i
0.2799+
0.0806i
0.5495+
0.1333i-
0.3954-
0.0678i-
0.2437-
0.0701i
00.1889-
0.0122i00-
0.2751-
0.0664i-
0.1885+
0.0132i
00.0632+
0.0033i0-
0.5376-
0.0976i0-
0.0629-
0.0023i0k=3各节点有功率不平衡量为dp=0-
0.0057-
0.0044-
0.0044-
0.0058各节点有功率不平衡量为dp=0-
0.
00570.0003-
0.0044-
0.0058各节点有功率不平衡量为dp=0-
0.
00570.
00030.0004-
0.0058各节点有功率不平衡量为dp=0-
0.
00570.
00030.
00040.0005dtt=
1.0e-003*-
0.2442-
0.1346-
0.1413-
0.1674各节点电压的相位角dt为dt=0-
0.0461-
0.0838-
0.0895-
0.1043dtt=
1.0e-003*-
0.2442-
0.1346-
0.1413-
0.1674各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0895-
0.1043dtt=
1.0e-003*-
0.2442-
0.1346-
0.1413-
0.1674各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1043dtt=
1.0e-003*-
0.2442-
0.1346-
0.1413-
0.1674各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044各节点无功率不平衡量为qq=0-
0.0020-
0.0005-
0.0006-
0.0001各节点无功率不平衡量为qq=0-
0.
00200.0001-
0.0006-
0.0001各节点无功率不平衡量为qq=0-
0.
00200.
00010.0001-
0.0001各节点无功率不平衡量为qq=0-
0.
00200.
00010.
00010.0001各节点电压的修正量ut=
1.0e-004*-
0.9362000-
0.5617000-
0.6002000-
0.7119000各节点电压的大小为u=
1.
06001.
03651.
00881.
00731.0016各节点电压的修正量ut=
1.0e-004*-
0.9362000-
0.5617000-
0.6002000-
0.7119000各节点电压的大小为u=
1.
06001.
03651.
00871.
00731.0016各节点电压的修正量ut=
1.0e-004*-
0.9362000-
0.5617000-
0.6002000-
0.7119000各节点电压的大小为u=
1.
06001.
03651.
00871.
00721.0016各节点电压的修正量ut=
1.0e-004*-
0.9362000-
0.5617000-
0.6002000-
0.7119000各节点电压的大小为u=
1.
06001.
03651.
00871.
00721.0015各节点功率s为S1=
1.2987+
0.2445i各节点功率sij为S=
00.8900+
0.1387i
0.4087+
0.1058i00-
0.8755-
0.0954i
00.2468+
0.0815i
0.2792+
0.0806i
0.5488+
0.1333i-
0.3960-
0.0677i-
0.2430-
0.0701i
00.1891-
0.0121i00-
0.2745-
0.0665i-
0.1887+
0.0132i
00.0633+
0.0033i0-
0.5369-
0.0977i0-
0.0630-
0.0023i0k=4各节点有功率不平衡量为dp=
1.0e-003*
00.
73710.
26180.
44250.5138各节点有功率不平衡量为dp=
1.0e-003*
00.7371-
0.
04820.
44250.5138各节点有功率不平衡量为dp=
1.0e-003*
00.7371-
0.0482-
0.
06390.5138各节点有功率不平衡量为dp=
1.0e-003*
00.7371-
0.0482-
0.0639-
0.1022dtt=
1.0e-004*
0.
29030.
13660.
14310.1567各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044dtt=
1.0e-004*
0.
29030.
13660.
14310.1567各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044dtt=
1.0e-004*
0.
29030.
13660.
14310.1567各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044dtt=
1.0e-004*
0.
29030.
13660.
14310.1567各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044各节点无功率不平衡量为qq=
1.0e-003*
00.
25640.
09530.
09800.1310各节点无功率不平衡量为qq=
1.0e-003*
00.2564-
0.
01640.
09800.1310各节点无功率不平衡量为qq=
1.0e-003*
00.2564-
0.0164-
0.
01490.1310各节点无功率不平衡量为qq=
1.0e-003*
00.2564-
0.0164-
0.0149-
0.0256各节点电压的修正量ut=
1.0e-004*
0.
11300000.
06130000.
06560000.0753000各节点电压的大小为u=
1.
06001.
03651.
00871.
00721.0015各节点电压的修正量ut=
1.0e-004*
0.
11300000.
06130000.
06560000.0753000各节点电压的大小为u=
1.
06001.
03651.
00881.
00721.0015各节点电压的修正量ut=
1.0e-004*
0.
11300000.
06130000.
06560000.0753000各节点电压的大小为u=
1.
06001.
03651.
00881.
00731.0015各节点电压的修正量ut=
1.0e-004*
0.
11300000.
06130000.
06560000.0753000各节点电压的大小为u=
1.
06001.
03651.
00881.
00731.0016各节点功率s为S1=
1.2981+
0.2445i各节点功率sij为S=
00.8894+
0.1387i
0.4087+
0.1058i00-
0.8750-
0.0954i
00.2469+
0.0815i
0.2793+
0.0806i
0.5489+
0.1333i-
0.3960-
0.0677i-
0.2431-
0.0701i
00.1891-
0.0121i00-
0.2746-
0.0664i-
0.1887+
0.0132i
00.0633+
0.0033i0-
0.5370-
0.0977i0-
0.0630-
0.0023i0k=5各节点有功率不平衡量为dp=
1.0e-003*0-
0.0924-
0.0482-
0.0639-
0.1022各节点有功率不平衡量为dp=
1.0e-003*0-
0.
09240.0067-
0.0639-
0.1022各节点有功率不平衡量为dp=
1.0e-003*0-
0.
09240.
00670.0080-
0.1022各节点有功率不平衡量为dp=
1.0e-004*0-
0.
92370.
06670.
07970.1255dtt=
1.0e-005*-
0.3627-
0.1670-
0.1760-
0.1966各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044dtt=
1.0e-005*-
0.3627-
0.1670-
0.1760-
0.1966各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044dtt=
1.0e-005*-
0.3627-
0.1670-
0.1760-
0.1966各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044dtt=
1.0e-005*-
0.3627-
0.1670-
0.1760-
0.1966各节点电压的相位角dt为dt=0-
0.0461-
0.0839-
0.0896-
0.1044各节点无功率不平衡量为qq=
1.0e-004*0-
0.3235-
0.1636-
0.1492-
0.2558各节点无功率不平衡量为qq=
1.0e-004*0-
0.
32350.0199-
0.1492-
0.2558各节点无功率不平衡量为qq=
1.0e-004*0-
0.
32350.
01990.0180-
0.2558各节点无功率不平衡量为qq=
1.0e-004*0-
0.
32350.
01990.
01800.0276各节点电压的修正量ut=
1.0e-005*-
0.1456000-
0.0813000-
0.0872000-
0.1025000各节点电压的大小为u=
1.
06001.
03651.
00881.
00731.0016各节点电压的修正量ut=
1.0e-005*-
0.1456000-
0.0813000-
0.0872000-
0.1025000各节点电压的大小为u=
1.
06001.
03651.
00871.
00731.0016各节点电压的修正量ut=
1.0e-005*-
0.1456000-
0.0813000-
0.0872000-
0.1025000各节点电压的大小为u=
1.
06001.
03651.
00871.
00731.0016各节点电压的修正量ut=
1.0e-005*-
0.1456000-
0.0813000-
0.0872000-
0.1025000各节点电压的大小为u=
1.
06001.
03651.
00871.
00731.0016各节点功率s为S1=
1.2982+
0.2445i各节点功率sij为S=
00.8895+
0.1387i
0.4087+
0.1058i00-
0.8751-
0.0954i
00.2469+
0.0815i
0.2793+
0.0806i
0.5489+
0.1333i-
0.3960-
0.0677i-
0.2431-
0.0701i
00.1891-
0.0121i00-
0.2746-
0.0664i-
0.1887+
0.0132i
00.0633+
0.0033i0-
0.5370-
0.0977i0-
0.0630-
0.0023i0k=6附录感想电力系统的潮流分布计算非常重要,它涉及到系统运行的方方面面在电力系统运行方式和规划方案的研究中,都需要进行潮流分布计算以比较运行方式或规划供电方案的可行性,可靠性和经济性同时为了满足实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算因此,潮流计算是电力系统中应用最广泛,最基本和最中央的一种电气运算说实话看到这个设计时很是头疼,matlab没有接触过(这是第一次)c语言也忘得差不多了,只好从零开始学习matlab,但由于种种原因学的乱七八糟的,没怎么明白只好找了些参考,整完了这次的论文PAGE37。