发布时间 : 星期一 文章电力系统短路计算的程序实现更新完毕开始阅读044c7f06bed5b9f3f90f1c59
电力系统短路计算的程序实现
致谢
大学四年的学习生活即将结束了，本次毕业设计也已经接近尾声，作为一个本科生的毕业设计，由于经验的匮乏，难免有许多考虑不周全的地方，但在各位老师的督促指导下，以及同学们的关心和帮助下，我的各方面都有所提高，整体素质也有了很大的进步。这一切都是和我系的老师对我的指引和教育是分不开的。
首先我要衷心感谢我的指导老师赵熹老师，这篇论文是在赵老师的悉心指导下，结合相关的书籍、资料完成的，更是倾注了她大量的心血。赵老师在毕业设计过程中给予了我细心指导和无私帮助，从初稿到定稿，赵老师不厌其烦，一审再审，大到篇章布局的偏颇，小到语句格式的瑕疵，都一一予以指出，帮助我解决了论文中出现的许多难题。她的教诲，令我终身难忘。在论文完成之际，特向赵老师致以衷心的感谢！
同时感谢学校各位领导和老师，是他们在我学业和生活中给予了极大的关心和帮助。是他们的默默工作，为我们所有的学生创建了一个安定和有序的学习环境。
感谢我的室友们，从遥远的家来到这个陌生的城市里，是你们和我共同维系着彼此之间兄弟般的感情，维系着寝室那份家的融洽。四年的时光，仿佛就在昨天。
谨向我的父母和家人表示诚挚的谢意。他们是我生命中永远的依靠和支持，他们无微不至的关怀，是我前进的动力；他们的殷殷希望，激发我不断前行。没有他们就没有我，我的点滴成就都来自他们。
感谢老师们百忙中审阅我的论文！
在论文即将完成之际，我的心情非常激动，从开始进入课题到论文的顺利完成，有多少可敬的师长、同学、朋友给了我无言的帮助，在这里请接受我诚挚的谢意！
-44-
电力系统短路计算的程序实现
参考文献
[1] 陈跃.电力系统分册(第二版). 中国水利水电出版社, 2008. [2] 温步瀛. 电力工程基础. 中国电力出版社, 2006.
[3] 刘天琪 邱晓燕. 电力系统分析概论. 科学出版社, 2005.
[4] 陈亚明. 电力系统计算程序及其应用. 水利电力出版社, 2000.
[5] 曹绳敏. 电力系统课程设计及毕业设计参考资料. 水利电力出版社, 1995.
[6] 同济大学应用数学系.线性代数.高等教育出版社, 2005.
[7] 李维波.MATLAB在电气工程中的应用. 中国电力出版社,2006. [8] 张文勤.电力系统清楚.中国电力出版社, 1998.
[9] 李光琦. 电力系统暂态分析. 中国电力出版社, 1994. [10] 刘万顺. 电力系统故障分析. 中国电力出版社, 1998. [11] 周光荣. 电力系统故障分析. 清华大学出版社, 1998. [12] 马大强. 电力系统机电暂态过程. 水利电力出版社, 1998. [13] 杨冠强. 电力系统自动装置原理. 水利电力出版社, 1995. [14] 夏道止. 电力系统分析. 浙江大学出版社. 1993.
-45-
电力系统短路计算的程序实现
附件1对称短路源程序
clear
ZZ(1,2)=j*0.1; ZZ(1,3)=j*0.1;ZZ(2,3)=j*0.1;
j*10j*10???j*26.6266??YB??j*10?j*33.2933j*10?;
?j*10j*10?j*19.96???n=3; k=3; for i=1:n if i= =k Ⅱ(1)=1; else
Ⅱ(1)=0; end end
Z(:,k)= YB\\Ⅱˊ; Zk=Z(:,k) k,Ik=1/Z(k,k) for i=1:n
U(i)=1-Z(i,k)*Ik; end Un=Uˊ for i=1:n for j=1:n if i I(i,j)=(U(i)-U(j))/ZZ(i,j); ij(1)=i;ij(2)=j; ij,Iij=I(i,j) end end end -46- 电力系统短路计算的程序实现 附件2不对称短路源程序 clear ZZ1(1,2)=j*0.1; ZZ1(1,3)=j*0.1; ZZ1(2,3)=j*0.1; ZZ2(1,2)=j*0.1; ZZ2(1,3)=j*0.1; ZZ2(2,3)=j*0.1; ZZ0(1,2)=j*0.2; ZZ0(1,3)=j*0.2; ZZ0(2,3)=j*0.2; ??j*26.6266j*10j*10?Y1= ??j*10?j*33.2933j*10?; ??j*10j*10?j*19.96????j*10j*10?Y2= ??j*26.6266?j*10?j*33.2933j*10?; ??j*10j*10?j*19.96????Y0= ??j*30j*5j*5??j*5?j*50j*5?; ??j*5j*5?j*10??????j*39.96j*10j*10j*20?j*10?j*59.96j*100YY1= ??j*10j*10?j*19.960?j*2000?j*30??0j*4000YY2=YY1; N1=3; N2=5; k=3; fault=1; for p=1:N1 if p= =k I(p)=1; else I(p)=0; -47- 0?j*40?0??0?; ??j*60??