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

　　

【正文】 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStructureAndBranchPara(CircleNumber,2))+... +1/(TopoStructureAndBranchPara(CircleNumber,3)+... 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 洛陽理工學(xué)院畢業(yè)設(shè)計論文 39 end % [fname,pname] = uigetfile(39。*.dat39。,39。Select Node datafile39。)。 Nodename=strcat(pname,fname)。 NodeData = csvread(Nodename)。 floatmin=。 datamax=0。 m=1。 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)NodeDa洛陽理工學(xué)院畢業(yè)設(shè)計論文 40 ta(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。 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)))。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 41 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 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 洛陽理工學(xué)院畢業(yè)設(shè)計論文 42 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 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。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 43 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 if NodeData(l,1)==1 洛陽理工學(xué)院畢業(yè)設(shè)計論文 44 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)。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 45 m=m+1。 if NodeData(i,1)==1 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)))。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 46 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。 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)))。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 47 end A(n)=sum。 n=n+1。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 datamax=max(datamax,B(m))。