【正文】 ntf( ,)。 printf(有功功率增量 P[2]=%f,P[2])。 printf(\n)。 printf( ,)。 printf(電壓增量 V22=%f,V22)。當(dāng) ji? 時(shí) 22()0iiij i ij ijjiiij i ij ijjiijjPQG e B fefPQB e G ffeVVef??? ??? ? ? ? ? ??? ???? ???? ? ? ??? ???? ??????? ?? 當(dāng) ij? 時(shí) ????????????????????????????????????????????????????????????????????iiiiiinkiiiiiikikkikiinkiiiiiikikkikiinkiiiiiikikkikiinkiiiiiikikkikiiffVeeVfBeGfBeGfQfGeBeBfGeQfGeBeBfGfPfBeGfBeGeP22)()()()(221111 以下為程序 ： //****形成雅克比矩陣 ********************** for(j=1。j++) {if(1==j) {float c=0,d=0。 for(m=1。m++) {c+=G[1][m]*e[m]B[1][m]*f[m]。 } 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 18 J[1*N1][j*N1]=cG[1][j]*e[1]B[1][j]*f[1]。 J[1*N][j*N1]=d+B[1][j]*e[1]G[1][j]*f[1]。 } else {J[1*N1][j*N1]=G[1][j]*e[1]B[1][j]*f[1]。 J[1*N1][j*N]=B[1][j]*e[1]G[1][j]*f[1]。 } } for(j=1。j++) {if(2==j) {float c=0,d=0。 for(m=1。m++) {c+=G[2][m]*e[m]B[2][m]*f[m]。 } J[2*N1][j*N1]=cG[2][j]*e[2]B[2][j]*f[2]。 J[2*N][j*N1]=2*e[2]。 } else {J[2*N1][j*N1]=G[2][j]*e[2]B[2][j]*f[2]。 J[2*N1][j*N]=B[2][j]*e[2]G[2][j]*f[2]。 } } printf(雅克比矩陣是 :\n)。i5。j5。 printf( )。 } LU 分解法求修正方程 LU 分解 ,又稱 Gauss 消去法 ,可把任意方陣分解成下三角矩陣的基本變換形式 (行交換 )和上三角矩陣的乘積。 其中 L 為下三角矩陣的基本變換形式 ,U為上三角矩陣。原始的求法 x=A^（ 1） *b,某些情況下，如果矩陣A 中的數(shù)非常小，我認(rèn)為不是因?yàn)榇髷?shù)除以小數(shù)誤差大么， 1/A 算出的誤差會(huì)很大。 求修正方程的程序如下 //********計(jì)算修正方程 ************* for(i=1。i++) {L[i][i]=1。iM。 L[i][1]=J[i][1]/U[1][1]。nM。jM。 for(s=0。s++) sigma1+=L[n][s]*U[s][j]。 } for(i=n。i++) { sigma2=0。s=n1。 L[i][n]=(J[i][n]sigma2)/U[n][n]。 b[2]=Q[1]。 b[4]=V22。iM。 for(n=1。n++) sigma1+=L[i][n]*y[n]。 } for(i=M1。i) { sigma2=0。nM。 x[i]=(y[i]sigma2)/U[i][i]。xe[2]=x[3]。xf[2]=x[4]。 for(i=1。i++) {e[i]+=xe[i]。 } for(i=1。i++) {printf(e[%d]=,i)。 printf( ,)。i3。 printf(%f,f[i])。 } printf(\n) 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 22 計(jì)算網(wǎng)絡(luò)中功率分布 最后要計(jì)算出平衡節(jié)點(diǎn)的功率和網(wǎng)絡(luò)中的功率分布。 經(jīng)過查找相關(guān)的資料得到： “ 多年的實(shí)踐證明，牛頓法具有很好的二次收斂性，是求解多元非線性方程的經(jīng)典算法，至今仍是電力系統(tǒng)潮流計(jì)算的主流。 ” —— 中國自動(dòng)化網(wǎng) . 改善調(diào)度員潮流計(jì)算收斂性的措施 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 25 通過本次電力系統(tǒng)分析課程設(shè)計(jì)，使我了解了自己在哪些方面有缺陷。于是我認(rèn)真地看電力系統(tǒng)分析下冊(cè)有關(guān)潮流計(jì)算的牛頓 拉夫遜法！了解了何謂 PQ 節(jié)點(diǎn)， PV節(jié)點(diǎn)，平衡節(jié)點(diǎn)等。 1 是 PQ 節(jié)點(diǎn)，節(jié)點(diǎn) 2 是PV節(jié)點(diǎn)，節(jié)點(diǎn) 3 是平衡節(jié)點(diǎn)。后來在網(wǎng)上找到了有關(guān) C 語言的相關(guān)程序！經(jīng)過了自己的仔細(xì)研讀，了解了各個(gè)函數(shù)的作用。因此在今后的學(xué)習(xí)過程中要端正學(xué)習(xí)態(tài)度。為將來的學(xué)習(xí)和工作打下良好的基礎(chǔ)！ 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 26 參考文獻(xiàn) [1] 何仰贊等 .電力系統(tǒng)分析上冊(cè) [M]． 武漢：華中理工大學(xué)出版社 . [2] 何仰贊等 .電力系統(tǒng)分析下冊(cè) [M]． 武漢：華中理工大學(xué)出版社 . [3] 諸俊偉等 .電力系統(tǒng)分析 [M]. 北京：中國電力出版社， 1995. [4] 周全仁等 .電網(wǎng)計(jì)算與程序設(shè)計(jì) [M].長(zhǎng)沙：湖南科學(xué)技術(shù)出版社 ,1983. [5]丁化成 .單片機(jī)應(yīng)用技術(shù) [A]. 北京 :北京航空航天大學(xué)出版社 ,2020. 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 27 附錄： include include include define N 2 define M 5 main() {double G[4][4],B[4][4],J[5][5]。 float Qs[4]={0,1},V2s=1, float V22,max,P3,Q3。 int i,j,n,s,k1=0。 //****************計(jì)算導(dǎo)納矩陣 ******************* G[1][1]=。 G[2][2]=。 G[3][3]=。 G[1][2]=G[2][1]=。 G[1][3]=G[3][1]=。 G[2][3]=G[3][2]=。 for(i=1。i++) {for(j=1。j++) 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 28 {printf(%f+(%f)j,G[i][j],B[i][j])。 } printf(\n)。 //******************************************** //計(jì)算各節(jié)點(diǎn)不平衡量 loop1: printf(迭代次數(shù) k1=%d\n,k1)。i3。 for(j=1。j++) {a+=G[i][j]*e[j]B[i][j]*f[j]。 } P[i]=Ps[i](e[i]*a+f[i]*b)。//計(jì)算無功功率的增量 V22=V2s*V2se[2]*e[2]。 printf( ,)。 printf(有功功率增量 P[2]=%f,P[2])。 printf(\n)。 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 29 printf( ,)。 printf(電壓增量 V22=%f,V22)。 //************篩選出最大值 *********************** max=fabs(P[1])fabs(P[2])?fabs(P[1]):fabs(P[2])。 max=maxfabs(V22)?max:fabs(V22)。 //******************************************** while (k1=4) { //****形成雅克比矩陣 ********************** for(j=1。j++) {if(1==j) {float c=0,d=0。 for(m=1。m++) {c+=G[1][m]*e[m]B[1][m]*f[m]。 } J[1*N1][j*N1]=cG[1][j]*e[1]B[1][j]*f[1]。 J[1*N][j*N1]=d+B[1][j]*e[1]G[1][j]*f[1]。 } 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 30 else {J[1*N1][j*N1]=G[1][j]*e[1]B[1][j]*f[1]。 J[1*N1][j*N]=B[1][j]*e[1]G[1][j]*f[1]。 } } for(j=1。j++) {if(2==j) {float c=0,d=0。 for(m=1。m++) {c+=G[2][m]*e[m]B[2][m]*f[m]。 } J[2*N1][j*N1]=cG[2][j]*e[2]B[2][j]*f[2]。 J[2*N][j*N1]=2*e[2]。 } else {J[2*N1][j*N1]=G[2][j]*e[2]B[2][j]*f[2]。 J[2*N1][j*N]=B[2][j]*e[2]G[2][j]*f[2]。 } } printf(雅克比矩陣是 :\n)。i5。j5。 printf( )。 } //********計(jì)算修正方程 ************* for(i=1。i++) {L[i][i]=1。iM。 L[i][1]=J[i][1]/U[1][1]。nM。jM。 for(s=0。s++) sigma1+=L[n][s]*U[s][j]。 } for(i=n。i++) { 武 漢理工大學(xué)《 電力系統(tǒng)分析 》課程設(shè)計(jì)說明書 32 sigma2=0。s=n1。 L[i][n]=(J[i][n]sigma2)/U[n][n]。 b[2]=Q[1]。 b[4]=V22。iM。 for(n=1。n++) sigma1+=L[i][n]*y[n]。 } for(i=M1。i) { sigma2=0。nM。 x[i]=(y[i]sigma2)/U[i][i]。xe[2]=x[3]。xf[2]=x[4]。 for(i=1。i++) {e[i]+=xe[i]。 } for(i=1。i++) {printf(e[%d]=,i)。 printf( ,)。i3。 printf(%f,f[i])。 } printf(\n)。 goto loop1。j4。 b1+=G[3][j]*f[j]+B[3][j]*e[j]。 Q3=f[3]*a1e[3]*b1。 printf(\n)