(完整版)基于MATLAB牛顿拉夫逊法进行潮流计算.docVIP

%本程序的功能是用牛顿拉夫逊法进行潮流计算 n=input( 请输入节点数 :n=); nl=input( 请输入支路数 :nl=); isb=input( 请输入平衡母线节点号 :isb=); pr=input( 请输入误差精度 :pr=); B1=input( 请输入由各支路参数形成的矩阵 :B1=); B2=input( 请输入各节点参数形成的矩阵 :B2=); 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)+B1(i,4)./2; end %求导纳矩阵 disp(导纳矩阵 Y=); disp(Y); G=real(Y);B=imag(Y); for i=1:n e(i)=real(B2(i,3)); f(i)=imag(B2(i,3)); V(i)=B2(i,4); e