【導(dǎo)讀】節(jié)點的已知量求解未知量，潮流計算在電力系統(tǒng)中有著獨特的作用。擴(kuò)建中各種計算提供必要的依據(jù)。計算潮流分布的方法很多，本設(shè)計主要用的是基于極坐標(biāo)的牛頓－拉夫遜法。我們可以畫出計算框圖，用MATLAB編寫出程序，來代替?zhèn)鹘y(tǒng)。復(fù)雜電力系統(tǒng)是一個包括大量母線、支路的龐大系統(tǒng)。統(tǒng)進(jìn)行潮流分析時，采用人工計算的方法已經(jīng)不再適用。計算機(jī)計算已逐漸成為。來驗證本設(shè)計在實際運行中的優(yōu)越性。

　　

【正文】 tureAndBranchPara(CircleNumber,4))) 。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStruxx 理 工學(xué)院畢業(yè)設(shè)計論文 28 ctureAndBranchPara(CircleNumber,1))=NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))。 else NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1))=... NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,1))+... +1/(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))+j*TopoStructureAndBranchPara(CircleNumber,5)。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))=... NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))+... +1/(TopoStructureAndBranchPara(CircleNumber,3)+... xx 理 工學(xué)院畢業(yè)設(shè)計論文 29 j*TopoStructureAndBranchPara(CircleNumber,4))+j*TopoStructureAndBranchPara(CircleNumber,5)。 NAM(TopoStructureAndBranchPara( CircleNumber,1),TopoStructureAndBranchPara( CircleNumber,2))=... 1/(TopoStructureAndBranchPara(CircleNumber,3)+... j*TopoStructureAndBranchPara(CircleNumber,4))。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,1))=NAM(TopoStructureAndBranchPara(CircleNumber,1),TopoStructureAndBranchPara(CircleNumber,2))。 end end % [fname,pname] = uigetfile(39。*.dat39。,39。Select Node datafile39。)。 Nodename=strcat(pname,fname)。 NodeData = csvread(Nodename)。 floatmin=。 datamax=0。 m=1。 xx 理 工學(xué)院畢業(yè)設(shè)計論文 30 n=1。 for i=1:NumberOfNode1 if NodeData(i,1)==1 sum=0。 sum1=0。 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)))。 sum1=sum1+NodeData(w,4)*(real(NAM(i,w))*sind(NodeData(i,5)NodeData(w,5))imag(NAM(i,w))*cosd(NodeData(i,5)NodeData(w,5)))。 end A(n)=sum。 n=n+1。 A(n)=sum1。 n=n+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 datamax=max(datamax,B(m))。 m=m+1。 xx 理 工學(xué)院畢業(yè)設(shè)計論文 31 B(m)=NodeData(i,3)NodeData(i,4)*sum1。 datamax=max(datamax,B(m))。 m=m+1。 else sum=0。 for r=1:NumberOfNode sum=sum+NodeData(r,4)*(real(NAM(i,r))*cosd(NodeData(i,5)NodeData(r,5))+imag(NAM(i,r))*sind(NodeData(i,5)NodeData(r,5)))。 end A(n)=sum。 n=n+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 datamax=max(datamax,B(m))。 m=m+1。 end end while datamaxfloatmin m=1。 for i=1:NumberOfNode1 xx 理 工學(xué)院畢業(yè)設(shè)計論文 32 n=1。 for k=1:NumberOfNode1 if i==k YMatrix(m,n)=NodeData(i,4)^2*imag(NAM(i,i))。 if NodeData(i,1)==1 YMatrix(m,n)=YMatrix(m,n)+NodeData(i,4)*A(m+1)。 end n=n+1。 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)))。 n=n+1。 end if NodeData(k,1)==1 if i==k YMatrix(m,n)=NodeData(i,4)*real(NAM(i,i))A(m)。 n=n+1。 else xx 理 工學(xué)院畢業(yè)設(shè)計論文 33 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)))。 n=n+1。 end end end m=m+1。 if NodeData(i,1)==1 n=1。 for l=1:NumberOfNode1 if i==l YMatrix(m,n)=NodeData(i,4)^2*real(NAM(i,i))+NodeData(i,4)*A(m1)。 n=n+1。 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)))。 n=n+1。 end xx 理 工學(xué)院畢業(yè)設(shè)計論文 34 if NodeData(l,1)==1 if i==l YMatrix(m,n)=NodeData(i,4)*imag(NAM(i,i))A(m)。 n=n+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)))。 n=n+1。 end end end m=m+1。 end end X=inv(YMatrix)*B39。 datamax=0。 m=1。 for i=1:NumberOfNode1 NodeData(i,5)=NodeData(i,5)+X(m)。 m=m+1。 if NodeData(i,1)==1 xx 理 工學(xué)院畢業(yè)設(shè)計論文 35 NodeData(i,4)=NodeData(i,4)+X(m)。 m=m+1。 end end m=1。 n=1。 for i=1:NumberOfNode1 if NodeData(i,1)==1 sum=0。 sum1=0。 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)))。 sum1=sum1+NodeData(v,4)*(real(NAM(i,v))*sind(NodeData(i,5)NodeData(v,5))imag(NAM(i,v))*cosd(NodeData(i,5)NodeData(v,5)))。 end A(n)=sum