牛顿-拉夫逊法潮流计算matlab程序

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

1、牛顿-拉夫逊法潮流计算matlab程序主程序"PowerFlow_NR.m"function[bus_res,S_res]=PowerFlow_NR_2%牛顿-拉夫逊法解潮流方程的主程序[bus,line]=OpDF_;%打开数据文件的子程序，返回bus（节点数据）和line（线路数据）回主程序[nb,mb]=size(bus);[nl,ml]=size(line);%计算bus和line矩阵的行数和列数[bus,line,nPQ,nPV,nodenum]=Num_(bus,line);%对节点重新

2、排序的子程序Y=y_(bus,line);%计算节点导纳矩阵的子程序myf=fopen('Result.m','w');fprintf(myf,'---------------bylongdinhohe---------------http://blog.sina.com.cn/longdinhohe---------------nn');fclose(myf);%在当前目录下生成“Result.m”文件，写入节点导纳矩阵formatlongEPS=1.0e-10;%设定误差精度fort=1:100

3、%开始迭代计算，设定最大迭代次数为100，以便不收敛情况下及时跳出[dP,dQ]=dPQ_(Y,bus,nPQ,nPV);%计算功率偏差dP和dQ的子程序J=Jac_(bus,Y,nPQ);%计算雅克比矩阵的子程序UD=zeros(nPQ,nPQ);fori=1:nPQUD(i,i)=bus(i,2);%生成电压对角矩阵enddAngU=J[dP;dQ];dAng=dAngU(1:nb-1,1);%计算相角修正量dU=UD*(dAngU(nb:nb+nPQ-1,1));%计算电压修正量bus(1:n

4、PQ,2)=bus(1:nPQ,2)-dU;%修正电压bus(1:nb-1,3)=bus(1:nb-1,3)-dAng;%修正相角if(max(abs(dU))

5、计算线路的等效Yt和Ym的子程序，以计算线路潮流bus_res=bus_res_(bus);%计算节点数据结果的子程序S_res=S_res_(bus,line,YtYm);%计算线路潮流的子程序myf=fopen('Result.m','a');fprintf(myf,'---------------牛顿－拉夫逊法潮流计算结果----------nn节点计算结果：n节点节点电压节点相角（角度）节点注入功率n');fori=1:nbfprintf(myf,'%2.0f',bus_res(i,1));

6、fprintf(myf,'%10.6f',bus_res(i,2));fprintf(myf,'%10.6f',bus_res(i,3));fprintf(myf,'%10.6f+j%10.6fn',real(bus_res(i,4)),imag(bus_res(i,4)));endfprintf(myf,'n线路计算结果：n节点I节点J线路功率S(I,J)线路功率S(J,I)线路损耗dS(I,J)n');fori=1:nlfprintf(myf,'%2.0f',S_res(i,1));fprint

7、f(myf,'%2.0f',S_res(i,2));fprintf(myf,'%10.6f+j%10.6f',real(S_res(i,3)),imag(S_res(i,3)));fprintf(myf,'%10.6f+j%10.6f',real(S_res(i,4)),imag(S_res(i,4)));fprintf(myf,'%10.6f+j%10.6fn',real(S_res(i,5)),imag(S_res(i,5)));endfclose(myf);%迭代结束后继续在“Result.m

8、”写入节点计算结果和线路计算结果程序结束子程序1"OpDF_.m"作用为打开数据文件function[bus,line]=OpDF_[dfile,pathname]=uigetfile('*.m','SelectDataFile');%数据文件类型为m文件，窗口标题为“SelectDataFile”ifpathname==0error('youmustselectavaliddatafile')%如果没有选择有效文件，则出现错误提示elselfile=length(d