【正文】 S(TopoStructureAndBranchPara(i,2),TopoStructureAndBranchPara(i,1))=(abs(V(TopoStructureAndBranchPara(i,2))))^2*(j)*TopoStructureAndBranchPara(i,5)+V(TopoStructureAndBranchPara(i,2))*... conj((V(TopoStructureAndBranchPara(i,2))V(TopoStructureAndBranchPara(i,1)))*1/(TopoStructureAndBranchPara(i,3)+j*TopoStructureAndBranchPara(i,4)))。 S(TopoStructureAndBranchPara(i,2),TopoStructureAndBranchPara(i,1))=(abs(V(TopoStructureAndBranchPara(i,2))))^2*(conj((1TopoStructureAndBranchPara(i,5))/(TopoStructureAndBranchPara(i,5)^2*(TopoStructureAndBranchPara(i,3)+j*TopoStructureAndBranchPara(i,4)))))+... V(TopoStructureAndBranchPara(i,2))*conj((V(TopoStructureAndBranchPara(i,2))V(TopoStructureAndBranchPara(i,1)))*(1/(TopoStructureAndBranchPara(i,5)*(TopoStructureAndBranchPara(i,3)+... j*TopoStructureAndBranchPara(i,4)))))。 end Balance=(NodeData(NumberOfNode,4)*sumNodeData(NumberOfNode,5)*sum1)+j*(NodeData(NumberOfNode,5)*sum+NodeData(NumberOfNode,4)*sum1)。 for u=1:NumberOfNode sum=sum+real(NAM(NumberOfNode,u))*real(V(u))imag(NAM(NumberOfNode,u))*imag(V(u))。 end sum=0。 m=m+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 end A(n)=sum。 else sum=0。 datamax=max(datamax,B(m))。 m=m+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 A(n)=sum1。 end A(n)=sum。 for v=1:NumberOfNode sum=sum+NodeData(v,4)*(real(NAM(i,v))*cosd(NodeData(i,5)NodeData(v,5))+imag(NAM(i,v))*sind(NodeData(i,5)NodeData(v,5)))。 for i=1:NumberOfNode1 if NodeData(i,1)==1 sum=0。 end end m=1。 if NodeData(i,1)==1 NodeData(i,4)=NodeData(i,4)+X(m)。 for i=1:NumberOfNode1 NodeData(i,5)=NodeData(i,5)+X(m)。 datamax=0。 end end end m=m+1。 else YMatrix(m,n)=NodeData(i,4)*(real(NAM(i,l))*sind(NodeData(i,5)NodeData(l,5))imag(NAM(i,l))*cosd(NodeData(i,5)NodeData(l,5)))。 end if NodeData(l,1)==1 if i==l YMatrix(m,n)=NodeData(i,4)*imag(NAM(i,i))A(m)。 else YMatrix(m,n)=NodeData(i,4)*NodeData(l,4)*(real(NAM(i,l))*cosd(NodeData(i,5)NodeData(l,5))+imag(NAM(i,l))*sind(NodeData(i,5)NodeData(l,5)))。 for l=1:NumberOfNode1 if i==l YMatrix(m,n)=NodeData(i,4)^2*real(NAM(i,i))+NodeData(i,4)*A(m1)。 end end end m=m+1。 else YMatrix(m,n)=NodeData(i,4)*(real(NAM(i,k))*cosd(NodeData(i,5)NodeData(k,5))+imag(NAM(i,k))*sind(NodeData(i,5)NodeData(k,5)))。 end if NodeData(k,1)==1 if i==k YMatrix(m,n)=NodeData(i,4)*real(NAM(i,i))A(m)。 else YMatrix(m,n)=NodeData(i,4)*NodeData(k,4)*(real(NAM(i,k))*sind(NodeData(i,5)NodeData(k,5))imag(NAM(i,k))*cosd(NodeData(i,5)NodeData(k,5)))。 if NodeData(i,1)==1 YMatrix(m,n)=YMatrix(m,n)+NodeData(i,4)*A(m+1)。 for i=1:NumberOfNode1 n=1。 m=m+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 end A(n)=sum。 else sum=0。 datamax=max(datamax,B(m))。 m=m+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 A(n)=sum1。 end A(n)=sum。 for w=1:NumberOfNode sum=sum+NodeData(w,4)*(real(NAM(i,w))*cosd(NodeData(i,5)NodeData(w,5))+imag(NAM(i,w))*sind(NodeData(i,5)NodeData(w,5)))。for i=1:NumberOfNode1 if NodeData(i,1)==1 sum=0。m=1。 floatmin=。Nodename=strcat(pname,fname)。Select Node datafile39。*.dat39。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,1))=NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))=... NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))+... +1/(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))+j*TopoStructureAndBranchPara(CircleNumber,5)。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,1))=NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))=... NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))+... 1/... (TopoStructureAndBranchPara(CircleNumber,5)^2*(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))) 。end%The following program is to form the Nodal Admittance Matrix % and the Topologic structure and Branch Parametres are arranged% I,J,R,X,C/K, and pay attention to the inpedence is in the side % Node I and the ratio of transformer is in the side of Node Jfor CircleNumber1=1:NumberOfBranch for CircleNumber2=1:NumberOfBranch NAM(CircleNumber1,CircleNumber2)=0。if Temporary1 Temporary2 NumberOfNode=Temporary1。Temporary1=max(TopoStructureAndBranchPara(:,1))。TopoStructureAndBranchPara= csvread(Branchname)。)。,39。謝謝!22參考文獻(xiàn)[1] ，1986.[2] 于永源. ，1996年.[3] [M].北京：中國(guó)電力出版社