【正文】 。 B(m)=NodeData(i,2)NodeData(i,4)*sum。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 47 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。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 46 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 洛陽理工學(xué)院畢業(yè)設(shè)計論文 44 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)。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 43 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 洛陽理工學(xué)院畢業(yè)設(shè)計論文 42 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。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 41 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))=... 洛陽理工學(xué)院畢業(yè)設(shè)計論文 38 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))) 。 洛陽理工學(xué)院畢業(yè)設(shè)計論文 35 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 J for CircleNumber1=1:NumberOfBranch for CircleNumber2=1:NumberOfBranch NAM(CircleNumber1,CircleNumber2)=0。 if Temporary1 Temporary2 NumberOfNode=Temporary1。 Temporary1=max(TopoStructureAndBranchPara(:,1))。 TopoStructureAndBranchPara= csvread(Branchname)。)。,39。謝謝 ! 洛陽理工學(xué)院畢業(yè)設(shè)計論文 30 洛陽理工學(xué)院畢業(yè)設(shè)計論文 31 參考文獻 [1] 楊以涵 .電力系統(tǒng)基礎(chǔ) .水利電力出版社， 1986. [2] 于永源 . 電力系統(tǒng)分析 .中國電力出版社， 1996 年 . [3] 諸俊偉等 .電力系統(tǒng)分析 [M].北京：中國電力出版社， 1995. [4] 楊少兵 ,駱平 . 電力系統(tǒng)分析的教學(xué)軟件開發(fā)，電力系統(tǒng)潮流分析 .華北電力技術(shù) ， 2020(10). [5] 韓禎祥，吳國炎等，電力系統(tǒng)分析，浙江大學(xué)出版社， 1993 年 . [6] 何仰贊等 .電力系統(tǒng)分析上冊 [M]．武漢：華中理工大學(xué)出版社 . [7] 李久勝等編 .電氣自動化英語 .哈爾濱：哈爾濱工業(yè)大學(xué)出版社，1999. [8] 尹克寧 .電力工程 .中國水利電力出版社， 1989. [9] 韋剛等 .電力系統(tǒng)分析要點與習(xí)題 .中國電力出版社 ， 2020. [10] 陳珩 .電力系統(tǒng)穩(wěn)態(tài)分析 . 水利電力出版社 ， 1994. [11] 陸敏政 .電力工程 .北京：中國電力出版社， 1996. [12] Steven Holzner(美 )著 . 詳實翻譯組譯 . Visual Basic 6 技術(shù)內(nèi)幕 . 機械工業(yè)出版社 ， 2020. [13] 諸俊偉等 .電力系統(tǒng)分析 [M].北京：中國電力出版社， 1995. [14] 張伯明 ,陳壽孫 .高等電力網(wǎng)絡(luò)分析 [M].北京：清華大學(xué)出版社1996. 洛陽理工學(xué)院畢業(yè)設(shè)計論文 32 [15] 周全仁等 .電網(wǎng)計算與程序設(shè)計 [M].長沙：湖南科學(xué)技術(shù)出版社 ， 1983. [16] 紀(jì)建偉等 .電力系統(tǒng)分析 .中國水利電力出版社， 1989. 洛陽理工學(xué)院畢業(yè)設(shè)計論文 33 附 錄 輸 入 原 始 數(shù) 據(jù)形 成 節(jié) 點 導(dǎo) 納 矩 陣設(shè) 電 壓 初 值設(shè) 迭 代 次 數(shù) k ＝ 0( 0 )e( 0 )