【正文】 xx 理工學院畢業(yè)設計論文 36 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 xx 理工學院畢業(yè)設計論文 35 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 xx 理工學院畢業(yè)設計論文 34 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 xx 理工學院畢業(yè)設計論文 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)))。 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 xx 理工學院畢業(yè)設計論文 32 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)+... xx 理工學院畢業(yè)設計論文 29 j*TopoStructureAndBranchPara(CircleNumber,4))+j*TopoStructureAndBranchPara(CircleNumber,5)。 NAM(TopoStructureAndBranchPara(CircleNumber,2),TopoStruxx 理工學院畢業(yè)設計論文 28 ctureAndBranchPara(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 xx 理工學院畢業(yè)設計論文 26 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。謝謝 ! xx 理工學院畢業(yè)設計論文 23 參考文獻 [1] 楊以涵 .電力系統(tǒng)基礎 .水利電力出版社， 1986. [2] 于永源 . 電力系統(tǒng)分析 .中國電力出版社， 1996 年 . [3] 諸俊偉等 .電力系統(tǒng)分析 [M].北京：中國電力出版社， 1995. [4] 楊少兵 ,駱平 . 電力系統(tǒng)分析的教學軟件開發(fā)，電力系統(tǒng)潮流分析 .華北電力技術(shù)， 2020(10). [5] 韓禎祥，吳國炎等，電力系統(tǒng)分析，浙江大學出版社， 1993 年 . [6] 何仰贊等 .電力系統(tǒng)分析上冊 [M]．武漢：華中理工大學出版社 . [7] 李久勝等編 .電氣自動化英語 .哈爾濱：哈爾濱工業(yè)大學出版社，1999. [8] 尹克寧 .電力工程 .中國水利電力出版社， 1989. [9] 韋剛等 .電力系統(tǒng)分析要點與習題 .中國電力出版社， 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)絡分析 [M].北京：清華大學出版社1996. [15] 周全仁等 .電網(wǎng)計算與程序設計 [M].長沙：湖南科學技術(shù)出版社，1983. [16] 紀建偉等 .電力系統(tǒng)分析 .中國水利電力出版社， 1989. xx 理工學院畢業(yè)設計論文 24 附 錄 輸 入 原 始 數(shù) 據(jù)形 成 節(jié) 點 導 納 矩 陣設 電 壓 初 值設 迭 代 次 數(shù) k ＝ 0( 0 )e( 0 )f計 算 誤 差 向 量()kp? ()kQ?2( )kU?收 斂 否否求 雅 可 比 矩 陣 元 素K = k + 1修 正 節(jié) 點 電 壓解 修 正 方 程 ，求 解()ke?()kf?( 1 ) ( ) ( )k k kf f f?? ? ?( 1 ) ( ) ( )k k ke e e? ? ? ?K m a xk不 收 斂 停 機是是求 p u 節(jié) 點 無 功 功 率 ， 求 平衡 節(jié) 點 功 率停 機求 支 路 功 率 分 布 和 損 耗停 機否 牛頓－拉夫遜法計算潮流程序框圖