资源描述：
《PQ分解法进行潮流计算程序.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、%本程序的功能是利用PQ分解法进行潮流计算clearclcn=input('请输入节点数：n=');nl=input('请输入支路数：nl=');isb=input('请输入平衡母线节点号：isb=');pr=input('请输入误差精度：pr=');disp('请输入由支路参数形成的矩阵：B1=');%输入B1[fname,pname]=uigetfile('*.xls');file=[pnamefname];x=xlsread(file);B1=[x(:,1),x(:,2),(x(:,3)+x(:,4)*i),x(:,5)*i,x
2、(:,6),x(:,7)];disp('请输入由节点参数形成的矩阵：B2=');%输入B2[fname,pname]=uigetfile('*.xls');file=[pnamefname];y=xlsread(file);B2=[(y(:,1)+y(:,2)*i),(y(:,3)+y(:,4)*i),y(:,5),y(:,6),y(:,7),y(:,8)];disp('请输入由节点号及其对地阻抗形成的矩阵：X=');%输入X[fname,pname]=uigetfile('*.xls');file=[pnamefname];z=xl
3、sread(file);X=z;na=input('请输入PQ节点数na=');Y=zeros(n);YI=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);fori=1:nifX(i,2)~=0;p=X(i,1);Y(p,p)=1./X(i,2);endendfori=1:nlifB1(i,6)==0p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
4、YI(p,q)=YI(p,q)-1./B1(i,3);Y(q,p)=Y(p,q);YI(q,p)=YI(p,q);Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;YI(q,q)=YI(q,q)+1./B1(i,3);Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;YI(p,p)=YI(p,p)+1./B1(i,3);end%求导纳矩阵G=real(Y);B=imag(YI);BI=imag(Y);fori=1:nS(i)=B2(i,1)-B2(i,2);BI(i,
5、i)=BI(i,i)+B2(i,5);endP=real(S);Q=imag(S);fori=1:ne(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4);endfori=1:nifB2(i,6)==2V(i)=sqrt(e(i)^2+f(i)^2);O(i)=atan(f(i)./e(i));endendfori=2:nifi==nB(i,i)=1./B(i,i);elseIC1=i+1;forj1=IC1:nB(i,j1)=B(i,j1)./B(i,i);endB(i,i)=1./B(
6、i,i);fork=i+1:nforj1=i+1:nB(k,j1)=B(k,j1)-B(k,i)*B(i,j1);endendendendp=0;q=0;fori=1:nifB2(i,6)==2p=p+1;k=0;forj1=1:nifB2(j1,6)==2k=k+1;A(p,k)=BI(i,j1);endendendendfori=1:naifi==naA(i,i)=1./A(i,i);elsek=i+1;forj1=k:naA(i,j1)=A(i,j1)./A(i,i);endA(i,i)=1./A(i,i);fork=i+1:n
7、aforj1=i+1:naA(k,j1)=A(k,j1)-A(k,i)*A(i,j1);endendendendICT2=1;ICT1=0;kp=1;kq=1;K=1;DET=0;ICT3=1;whileICT2~=0
8、ICT3~=0ICT2=0;ICT3=0;fori=1:nifi~=isbC(i)=0;fork=1:nC(i)=C(i)+V(k)*(G(i,k)*cos(O(i)-O(k))+BI(i,k)*sin(O(i)-O(k)));endDP1(i)=P(i)-V(i)*C(i);DP(i)=DP1(i)./V(i);DE
9、T=abs(DP1(i));ifDET>=prICT2=ICT2+1;endendendNp(K)=ICT2;ifICT2~=0fori=2:nDP(i)=B(i,i)*DP(i);ifi~=nIC1=i+1;fork=
