毕业论文--基于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);endfor i=1:n??? S(i)=B2(i,1)-B2(i,2);??? B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;while IT2~=0??? IT2=0;a=a+1;??? for i=1:n??????? if i~=isb ??????????? C(i)=0;??????????? D(i)=0;??????????? for j1=1:n???????????????? C(i)= C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);???????????????? D(i)= D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);??????????? end??????????? P1=C(i)*e(i)+f(i)*D(i);??????????? Q1=f(i)*C(i)-D(i)*e(i);??????????? V2=e(i)^2+f(i)^2;%117页malihong打if B2(i,6)~=3??? DP=P(i)-P1;??? DQ=Q(i)-Q1;??? for j1=1:n??????? if j1~=isbj1~=i??????????? X1=-G(i,j1)*e(i)-B(i,j1)*f(i);??????????? X2=B(i,j1)*e(i)-G(i,j1)*f(i);??????????? X3=X2;??????????? X4=-X1;??????????? p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;??????????? m=p+1;??????????? J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;???? elseif j1==ij1~=isb??????????? X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);??????????? X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);??????????? X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);??????????? X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);??????????? p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;??????????? J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;?????? end??? endelse??? DP=P(i)-P1;??? DV=V(i)^2-V2;??? for j1=1:n??????? if j1~=isbj1~=i???????????? X1=-G(i,j1)*e(i)-B(i,j1)*f(i);??????????? X2=B(i,j1)*e(i)-G(i,j1)*f(i);??????????? X5=0;???????