潮流计算matlab-牛顿拉夫逊法法.doc

《潮流计算matlab-牛顿拉夫逊法法.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、用matlab潮流计算(牛顿拉夫逊法)%主程序：[dfile,pathname]=uigetfile('ieee14.m','SelectDataFile');ifpathname==0error('youmustselectavaliddatafile')elselfile=length(dfile);%stripoff.meval(dfile(1:lfile-2));%打开数据文件endglobaln;globalm;[nb,mb]=size(bus);%节点重新编号[nl,ml]=size(line);nSW=0;%平衡节点数目n

2、PV=0;%PV节点数目nPQ=0;%PQ节点数目fori=1:nb,%nb为总节点数type=bus(i,6);iftype==3,nSW=nSW+1;%统计平衡节点数目SW(nSW,:)=bus(i,:);elseiftype==2,nPV=nPV+1;%统计PV节点数目PV(nPV,:)=bus(i,:);elsenPQ=nPQ+1;%统计PQ节点数目PQ(nPQ,:)=bus(i,:);endend;bus=[PQ;PV;SW];newbus=[1:nb]';f=bus(:,1);nodenum=[newbusbus(:,1)]

3、;bus(:,1)=newbus;fori=1:nlforj=1:2fork=1:nbifline(i,j)==nodenum(k,2)line(i,j)=nodenum(k,1);breakendendendendY=y(bus,line);%形成节点导纳矩阵K=0;%迭代次数初值Kmax=10;%最大迭代次数eps1=1.0e-10;eps2=1.0e-10;m=nPQ;n=nb;Um=eye(m,m);myf=fopen('output1.dat','w');forK=1:Kmaxfori=1:mforj=1:mifi==jUm(

4、i,j)=bus(i,2);endendendb=dPQ(Y,bus);C=jac(bus,Y);dX=Cb';dx=dX';[nx,mx]=size(dx);fori=1:n-1%计算相角bus(i,3)=bus(i,3)-dX(i,1);endB=dx(nx,n:mx)Um;%计算电压差bus(1:m,2)=bus(1:m,2)-B';%计算电压值dx(nx,n:mx)=B;fprintf(myf,'--第%d次迭代时雅可比矩阵--',K)fprintf(myf,'');fori=1:(n+m-1)forj=1:(n+m-1)

5、fprintf(myf,'%8.6f',C(i,j));fprintf(myf,'');endfprintf(myf,'');endfprintf(myf,'--第%d次迭代时dPQ的误差--',K)fprintf(myf,'');fori=1:(n+m-1)fprintf(myf,'%8.6e',b(1,i));fprintf(myf,'');endfprintf(myf,'');fprintf(myf,'--第%d次迭代时dx(误差)--',K)fprintf(myf,'');fori=1:(n+m-1)fpri

6、ntf(myf,'%8.6e',dX(i,1));fprintf(myf,'');endfprintf(myf,'');fprintf(myf,'第%d次迭代后节点电压(仅PQ节点)',K)fprintf(myf,'');fori=1:mfprintf(myf,'%8.6f',bus(i,2));fprintf(myf,'');endfprintf(myf,'');fprintf(myf,'第%d次迭代后相角(角度)',K)fprintf(myf,'');fori=1:nfprintf(myf,'%8.6f',b

7、us(i,3)180/pi);fprintf(myf,'');endfprintf(myf,'');if(max(abs(dx))

8、n%计算PV节点、平衡节点的无功Q1=0;T=0;forj=1:nT=bus(k,2)bus(j,2)(real(Y(k,j))sin(bus(k,3)-bus(j,3))-imag(Y(k,j))cos(