精选优质文档-倾情为你奉上%电力系统极坐标下的牛顿-拉夫逊法潮流计算disp(电力系统极坐标下的牛顿-拉夫逊法潮流计算:);clearn=input(请输入结点数：n=);n1=input(请输入PV结点数：n1=);n2=input(请输入PQ结点数：n2=);isb=input(请输入平衡结点：isb=);pr=input(请输入精确度：pr=);K=input(请输入变比矩阵看:K=);C=input(请输入支路阻抗矩阵：C=);y=input(请输入支路导纳矩阵：y=);U=input(请输入结点电压矩阵：U=);S=input(请输入各结点的功率：S=);Z=zeros(1,n);N=zeros(n1+n2,n2);L=zeros(n2,n2);QT1=zeros(1,n1+n2);for m=1:n for R=1:n C(m,m)=C(m,m)+y(m,R); if K(m,R)=0 C(
