5节点电力系统牛顿-拉夫逊法潮流计算讲述.doc

（ 二 〇 一 四 年 十 二 月 程序设计 % 本程序的功能是用牛顿拉夫逊法进行潮流计算 n=input(请输入节点数:n=); nl=input(请输入支路数:nl=); 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=zeros(1,n); O=zeros(1,n);S1=zeros(nl); for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B