牛顿-拉夫逊潮流计算.docVIP

牛顿-拉夫逊潮流计算 学生：李敏 2007109522 指导老师：李咸善 潮流计算原理：节点分为三类，PQ节点，PV节点，平衡节点，每个节点涉及4个量：P(i),Q(i),U(i),Q(i)根据各节点的功率，电压列写平衡方程，组成多维非线性方程，用数学方法（牛顿拉夫逊）将方程组转化为线性的，可以解方程组。在牛顿法中通过节点给定值可求得PQ节点和pv节点的电压相位和PQ节点的电压有效值U(i),及PV节点的Q(i),便可求得所有节点的参量，实现整个系统参数的明了化。牛顿拉夫逊的基本原理数学表述如下： 牛顿-拉夫逊法计算潮流 由线路参数得出节点导纳矩阵Y 功率方程： 其中 雅可比矩阵各元素的表达式 用MATLAB语言编程调试 心得体会： 李老师布置这个课外题已经很长时间了，一直没有静下心来做，因为前阵子在忙着学习cad考试耽搁了一些时间，所以没有及时的完成，这个任务我很妥了一段时间，觉得只有一点C语言基础想自己动手编一个潮流程序还是很困难的，所以就一直没有开始。到学期后期感觉就不想动手了，但是您一直支持我们完成它，前些天在您的提醒下我还是下决心还是抓紧时间完成。以前就很想了解matlab这个软件的，只知道它的功能很强大，但一直没能学习。 我很佩服也很喜欢李老师的，博学多识，为人师表，是我见过做好的老师之一。我很能理解老师的良苦用心，希望我们学生不要拘泥于课本，还要在课外多扩展自己的知识，丰富提高自己的专业素质，努力使自己成为电力行业的业内人士。我不能丢失掉这次机会，于是我从图书馆借来两本matab的书，看看基础，再结合课本上的潮流计算的内容，试着编写程序，在做的过程中，遇到了很多问题： 比如基本语言语法不会，矩阵的表示和运算； 比如那5个节点，开始我用2345为PQ节点，1为平衡节点，这与书上编号有些差别，这个也可以做出来节点导纳矩阵，但到了雅各比出现了一些问题，我怎么也没改过来，后来向做过的同学请教，也参考了一些他们的意见，后来还是换成标准形式就可以了，还有公式中循环的使用范围的确定等等 三天我什么也没做，早上很早就起来，晚上12点睡觉，除了吃饭就是看书，编程调试，修改，查书找指令，看类似的课件，问同学、、、实在很累，但我坚持下来了，看到我最后运行的结果，我感到很欣慰，终于印证了那一句话：天下无难事，只要肯攀登。 总之花了大概三天，我基本了解了matlab语言的运用，也知道用它来处理矩阵问题的优越性，不愧为“矩阵实验室”，也会用它来做一些工程计算，算是入门吧，但是它的功能实在强大，在以后的学习中在逐渐加以掌握。 第一次编程，难免会出现一些问题，如果有错误，恳请老师指正。 程序：%???·?-ê?êy?Y %FROM TO R X zdata=[5 1 0.02 0.06 5 2 0.08 0.24 1 2 0.06 0.18 1 3 0.06 0.18 1 4 0.04 0.12 2 3 0.01 0.03 3 4 0.08 0.24]; nl=zdata(:,1);nr=zdata(:,2);R=zdata(:,3);X=zdata(:,4); nbr=length(zdata(:,1));nbus=max(max(nl),max(nr)); Z=R+sqrt(-1)*X; y=ones(nbr,1)./Z; %?§?·μ??é Y=zeros(nbus,nbus); %??Y3?ê??ˉ?a0 for k=1:nbr %·??????a??μ?êy?μ if nl(k)0nr(k)0 Y(nl(k),nr(k))=Y(nl(k),nr(k))-y(k); Y(nr(k),nl(k))=Y(nl(k),nr(k)); end end for n=1:nbus %?????a??êy?μ for k=1:nbr if nl(k)==n|nr(k)==n Y(n,n)=Y(n,n)+y(k); else , end end end %Y?aμ??é???ó G