牛顿—拉夫逊法潮流计算MATLAB程序(共5页).doc

精选优质文档-倾情为你奉上牛顿拉夫逊法潮流计算程序By Yuluo%牛顿-拉夫逊法进行潮流计算n=input(请输入节点数：n=);n1=input(请输入支路数：n1=);isb=input(请输入平衡母线节点号：isb=);pr=input(请输入误差精度：pr=);B1=input(请输入由支路参数形成的矩阵：B1=);B2=input(请输入各节点参数形成的矩阵：B2=);X=input(请输入由节点参数形成的矩阵：X=);Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=seros(1,n);O=zeros(1,n);S1=zeros(n1);for i=1:n if X(i,2)=0; p=X(i,1); Y(p,p)=1./X(i,2); end end for i=1:n1 if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2)