最新潮流上机课程设计报告 联系客服

发布时间 : 星期六 文章最新潮流上机课程设计报告更新完毕开始阅读f0a7dff77c1cfad6195fa754

课程设计报告

( 2009-- 20010年度第二学期)

名 称:电力系统潮流上机 院 系: 班 级: 学 号: 学生姓名: 指导教师: 设计周数: 两周

成 绩:

日期:

课程设计报告

一、课程设计的目的与要求

培养学生的电力系统潮流计算机编程能力,掌握计算机潮流计算的相关知识

二、设计正文(详细内容见附录)

1. 手算

2. 计算机计算 3. 思考题

三、课程设计总结或结论

1 由于计算机潮流是上学期学过的内容,c语言的学习也是几学期前的事情了,所以此次试验前不得不好好去复习,在开始的实验过程中我有些盲目,不知从何下手,于是在参考了别人的节点导纳矩阵后,终于有了些思路,学过的内容大部分也想起来了。我发现此次的潮流程序大部分都用到了循环和嵌套,我也在这些方面出了很多的错,以至于不得不花很长时间去修改。

2 我们做的实验室九个节点的,其实节点的个数并不是很重要,只要程序编好了,只要在数组的规模上进行变动,不论是几个节点的都可以进行潮流计算。只是在实验中我又多次在循环的初值上我又多次犯了错误,错将i=0或i=1混淆,导致数据结果出错。另外我还更加意识到c语言格式书写的重要性。记得再算雅克比矩阵式由于一开始循环下的大括号太过混乱,语法并没有出错,导致数据结果出错,然我就接了很长时间,最后在同学的帮忙下才终于发现并改正了错误。 四、参考文献

1. 《电力系统计算:电子数字计算机的应用》,西安交通大学等合编。北京:水利电力出版社; 2. 《现代电力系统分析》,王锡凡主编,科学出版社; 3. 《电力系统稳态分析》,陈珩,中国电力出版社,2007年,第三版;

附录(包括:1.手算过程及结果;2.计算机计算流程图、表格、数据;3.思考题答案) 2、程序流程图 启动 输入原始数据

形成节点导纳矩阵

设非平衡节点电压初值ei(0)fi(0); 令迭代次数k=0 对PQ节点计算?Pi(k).?Qi(k)对PV节点计算?Pi(k),?Ui2(k)

令节点号i=1 1 课程设计报告

形成,i>n?

迭代是否收敛 ?e(k)max雅客比矩阵是否计算雅客比矩阵各元素(k)(k)(k)(k)(k) Hij,Nij,Jij,L(ijk),Rij,Sij增加迭代次数k=k+1 增加节点号i=i+1 解修正方程,由?Pi(k).?Qi(k)?Ui2(k)及雅客比矩阵用高斯法求各节点的电压增量?e求(k)i计算节点的新电压 ,?fi (k)ei(k?1)?ei(k)??ei(k) fi(k?1)?fi(k)??fi(k)出max?e(k),?f(k)max ,?f(k)max??计算平衡节点的功率及线路功率 停止 附录:程序中定义的变量名函数名的含义及功能 config.h文件中已经给出的变量名 变量名 gBus gLine gY_G,gY_B gDelta_P gDelta_Q gDelta_PQ gJaccobi gDelta_f,gDelta_e,gDelta_fe gf,ge flow.cpp文件中已经给出的模板函数名,

2

含义 存储节点相关信息 存储线路相关信息 存储节点导纳数据 存储有功不平衡量 PQ节点存储无功不平衡量 PV节点存储电压不平衡量 存储所有不平衡量 存储雅可比矩阵 存储修正量 存储电压值 课程设计报告

模板函数名 void test() void GetData() void GetYMatrix() void SetInitial() void GetUnbalance() void GetJaccobi() void GetRevised() void GetNewValue() 功能 测试NEquation.h文件的解方程功能 获得data.txt文件中的数据 计算节点导纳矩阵 设初值 计算不平衡量 计算雅可比矩阵 计算修正量 计算新的迭代值 自定义变量含义说明

自定义变量 P[Bus_Num][Bus_Num] Q[Bus_Num][Bus_Num] PS QS PW QW H[Bus_Num][Bus_Num] J[Bus_Num][Bus_Num] N[Bus_Num][Bus_Num] L[Bus_Num][Bus_Num] R[Bus_Num][Bus_Num] S[Bus_Num][Bus_Num] 含义 存储线路有功功率 存储线路无功功率 存储平衡节点有功功率 存储平衡节点无功功率 存储网损中有功功率 存储网损中无功功率 雅克比矩阵元素 雅克比矩阵元素 雅克比矩阵元素 雅克比矩阵元素 雅克比矩阵元素 雅克比矩阵元素 ---节点注入功率---

P(0)= 0.72212 Q(0)= 0.26902 P(1)= 1.63000 Q(1)= 0.04395 P(2)= 0.85000 Q(2)= -0.08345 P(3)= -0.00000 Q(3)= 0.00000 P(4)= -1.25000 Q(4)= -0.50000 P(5)= -0.90000 Q(5)= -0.30000 P(6)= 0.00000 Q(6)= 0.00000 P(7)= -1.00000 Q(7)= -0.35000 P(8)= -0.00000 Q(8)= -0.00000

---线路传输功率---

P_mn[0]= 0.72212,Q_mn[0]= 0.26902,P_nm[0]= -0.72212,Q_nm[0]= -0.23718 P_mn[1]= 1.63000,Q_mn[1]= 0.04395,P_nm[1]= -1.63000,Q_nm[1]= 0.11549 P_mn[2]= 0.85000,Q_mn[2]= -0.08345,P_nm[2]= -0.85000,Q_nm[2]= 0.12442 P_mn[3]= 0.76535,Q_mn[3]= -0.03690,P_nm[3]= -0.75477,Q_nm[3]= -0.07887 P_mn[4]= -0.24523,Q_mn[4]= -0.27113,P_nm[4]= 0.24625,Q_nm[4]= 0.06069 P_mn[5]= -0.84177,Q_mn[5]= -0.11956,P_nm[5]= 0.86466,Q_nm[5]= -0.07858 P_mn[6]= -0.59037,Q_mn[6]= -0.13023,P_nm[6]= 0.60375,Q_nm[6]= -0.18511 P_mn[7]= 0.41077,Q_mn[7]= 0.22215,P_nm[7]= -0.40823,Q_nm[7]= -0.38044 P_mn[8]= 0.31135,Q_mn[8]= 0.01503,P_nm[8]= -0.30963,Q_nm[8]= -0.16977

3