華電潮流上機(jī)課程設(shè)計(jì)報(bào)告程序

1、電力系統(tǒng)潮流上機(jī)課程設(shè)計(jì)報(bào)告院系：電氣與電子工程學(xué)院班級(jí)：學(xué)號(hào)：學(xué)生姓名：指導(dǎo)教師：劉寶柱設(shè)計(jì)周數(shù)：成績(jī)：日期： 2012年 1月 5日、課程設(shè)計(jì)的目的與要求培養(yǎng)學(xué)生的電力系統(tǒng)潮流計(jì)算機(jī)編程能力，掌握計(jì)算機(jī)潮流計(jì)算的相關(guān)知識(shí)二、設(shè)計(jì)正文（詳細(xì)內(nèi)容見(jiàn)附錄）1 手算2 計(jì)算機(jī)計(jì)算3 思考題3.1 潮流計(jì)算的方法有哪些？各有何特點(diǎn)？答： 潮流計(jì)算的方法主要有高斯 -賽德?tīng)柕?、牛頓 - 拉夫遜迭代法和 P-Q分解法。 它們各自的特點(diǎn)如下：（1）高斯 - 賽德?tīng)柕ǚ譃橐怨?jié)點(diǎn)導(dǎo)納矩陣為基礎(chǔ)的高斯-賽德?tīng)柕ê鸵砸怨?jié)點(diǎn)阻抗矩陣為基礎(chǔ)的高斯 - 賽德?tīng)柕ǖ脑肀容^簡(jiǎn)單，要求的數(shù)字計(jì)算機(jī)的內(nèi)存量也

2、比較小，但收斂性差，當(dāng)系統(tǒng)規(guī) 模變大時(shí)，迭代次數(shù)急劇上升，往往出現(xiàn)迭代不收斂的情況；而阻抗法改善了電力系統(tǒng)潮流計(jì)算導(dǎo)納法德 收斂性問(wèn)題，在當(dāng)時(shí)獲得了廣泛的應(yīng)用，但是，阻抗法的主要缺點(diǎn)是占用計(jì)算機(jī)的內(nèi)存很大，每次迭代的 計(jì)算量很大。當(dāng)系統(tǒng)不斷擴(kuò)大時(shí)，這些缺點(diǎn)就更加突出。（2）牛頓- 拉夫遜法是數(shù)學(xué)中求解非線性方程式的典型方法，有較好的收斂性。只要在迭代過(guò)程中盡 可能保持方程式系數(shù)矩陣的稀疏性， 就可以大大提高牛頓潮流計(jì)算程序的計(jì)算效率， 牛頓法在收斂性、 內(nèi) 存要求、計(jì)算速度方面都超過(guò)了阻抗法，成為知道目前仍被廣泛采用的方法。（3）P-Q 分解法潮流計(jì)算派生于以極坐標(biāo)表示時(shí)的牛頓-拉夫遜法，它根

3、據(jù)電力系統(tǒng)的特點(diǎn)，抓住主要矛盾，對(duì)純數(shù)學(xué)的牛頓法進(jìn)行了改造。與牛頓法相比，P-Q 分解法的修正方程的系數(shù)矩陣 B和 B”分別是（n-1 ）和（ m-1）的方陣，替代了原有的（ n+m-2）階系數(shù)矩陣 J ；B、 B”在迭代過(guò)程中可以保持不變 且為對(duì)稱的系數(shù)矩陣，提高了計(jì)算速度，降低了對(duì)存儲(chǔ)容量的要求。 P-Q 分解法在計(jì)算速度方面有顯著地 提高，迅速得到了推廣。3.2 如果交給你一個(gè)任務(wù)， 請(qǐng)你用已有的潮流計(jì)算軟件計(jì)算北京城市電網(wǎng)的潮流， 你應(yīng)該做哪些工作？ （收集哪些數(shù)據(jù)，如何整理，計(jì)算結(jié)果如何分析）答： （1） 在進(jìn)行北京城市電網(wǎng)的潮流計(jì)算之前需要了解北京城市電網(wǎng)中所有的節(jié)點(diǎn)支路的相關(guān)數(shù)據(jù)

4、， 并對(duì)節(jié)點(diǎn)和支路分類。處理 PQ節(jié)點(diǎn)時(shí)要了解節(jié)點(diǎn)的注入有功和無(wú)功功率；PV 節(jié)點(diǎn)要了解節(jié)點(diǎn)電壓大小，注入有功功率及節(jié)點(diǎn)所能提供的最大和最小無(wú)功功率；對(duì)于平衡節(jié)點(diǎn)要了解節(jié)點(diǎn)的電壓大小、相位及節(jié)點(diǎn) 所能提供的最大和最小有功無(wú)功功率，此外還需的節(jié)點(diǎn)數(shù)據(jù)有：支路的電阻電抗電納支路變壓器的變比及 線路所能傳輸?shù)淖畲笕萘康?。?） 計(jì)算結(jié)果分析：應(yīng)考慮 PQ節(jié)點(diǎn)的電壓是否過(guò)高或過(guò)低，分析 PV 節(jié)點(diǎn)的電壓幅值是否正常及無(wú)功 功率是否超出范圍，分析平衡節(jié)點(diǎn)有功無(wú)功功率是否在節(jié)點(diǎn)所能提供的范圍之內(nèi)，分析各支路的功率，看是否超出線路傳輸?shù)淖畲笕萘?，分析整個(gè)系統(tǒng)的網(wǎng)損是否達(dá)到標(biāo)準(zhǔn)3.3 設(shè)計(jì)中遇到的問(wèn)題和解決的

5、辦法。開(kāi)始時(shí)我們進(jìn)行原始數(shù)據(jù)的讀入和讀出，這一部分的程序比較簡(jiǎn)單，所以我們都很順利都寫出程 序。接下來(lái)就是節(jié)點(diǎn)導(dǎo)納矩陣的形成。編寫這一部分程序時(shí)我剛開(kāi)始沒(méi)有考慮到變壓器支路，所以當(dāng)原始 數(shù)據(jù)改為四節(jié)點(diǎn)四支路的數(shù)據(jù)時(shí)，節(jié)點(diǎn)導(dǎo)納矩陣的結(jié)果就出錯(cuò)了。經(jīng)過(guò)改正之后終于得出了正確的節(jié)點(diǎn)導(dǎo) 納矩陣。然后就是計(jì)算不平衡功率和雅可比矩陣。計(jì)算不平衡功率時(shí)候，為了方便計(jì)算，我將各節(jié)點(diǎn)的不 平衡功率放到一個(gè)一維數(shù)組里面。雅可比矩陣的形成這一部分上我用了很長(zhǎng)時(shí)間才將程序?qū)懗鰜?lái)。首先是 對(duì)角線上的 H,N,J,L 的計(jì)算公式寫的很復(fù)雜，與同學(xué)討論之后我發(fā)現(xiàn)可以先將每個(gè)節(jié)點(diǎn)的注入電流放到一 個(gè)數(shù)組里面，這樣后面用來(lái)求

6、H ,N ,J,L 時(shí)就比較方便。在這一模塊的過(guò)程中，我沒(méi)有注意到雅克比矩陣 中 PQ節(jié)點(diǎn)和 PV 節(jié)點(diǎn)的對(duì)應(yīng)關(guān)系， 而且我沒(méi)有將節(jié)點(diǎn)號(hào)取出來(lái)， 以致出現(xiàn)了錯(cuò)誤的編程思路和錯(cuò)誤的結(jié)果， 后來(lái)在同組同學(xué)的幫助下我又進(jìn)行了第二次的編程，思路是正確的，但是運(yùn)行出來(lái)的結(jié)果卻是錯(cuò)誤的，我 想應(yīng)該是我哪個(gè)公式出了錯(cuò)， 最終經(jīng)過(guò)反復(fù)檢查后終于得到了正確結(jié)果， 在編寫雅克比矩陣還需要注意 H、 J、N、 L 四個(gè)矩陣的形成以及循環(huán)的順序，然后就是如何將雅可比矩陣的各個(gè)元素放到一個(gè)二維數(shù)組里， 這一部分看起來(lái)很簡(jiǎn)單，但是卻很容易出錯(cuò)，經(jīng)過(guò)和同組同學(xué)的認(rèn)真思考和檢查，終于得出來(lái)了正確的雅 可比矩陣。三、課程設(shè)計(jì)總

7、結(jié)或結(jié)論 兩周的潮流上機(jī)計(jì)算已經(jīng)結(jié)束了，雖然時(shí)間很短，可是我卻學(xué)到了很多。潮流上機(jī)計(jì)算除了需要一定 的電力系統(tǒng)穩(wěn)態(tài)知識(shí)以外， 還要有良好的 C語(yǔ)言基礎(chǔ)。 潮流上機(jī)的過(guò)程讓我們鞏固了潮流計(jì)算過(guò)程和方法， 以及穩(wěn)態(tài)的基礎(chǔ)知識(shí)。此外，從潮流數(shù)據(jù)的讀入和寫出，節(jié)點(diǎn)導(dǎo)納矩陣的形成，功率不平衡量的計(jì)算，到 形成雅克比矩陣和功率損耗的計(jì)算，我不斷地進(jìn)行編程和學(xué)習(xí)，在兩周的潮流上機(jī)過(guò)程中，我完成了程序 的編寫，掌握了計(jì)算機(jī)潮流計(jì)算的基本方法與步驟，了解了計(jì)算機(jī)潮流計(jì)算的優(yōu)點(diǎn)與重要性。潮流計(jì)算最關(guān)鍵的是進(jìn)行 C 語(yǔ)言的編程，沒(méi)有良好的 C語(yǔ)言基礎(chǔ)很難快速地編出潮流計(jì)算程序，編寫 時(shí)需要注意循環(huán)和嵌套的使用，我在

8、編寫的過(guò)程中由于循環(huán)沒(méi)有用好以及粗心等原因，導(dǎo)致多次編出來(lái)的 程序出現(xiàn)錯(cuò)誤，這在潮流計(jì)算中都是不應(yīng)該出現(xiàn)的。四、參考文獻(xiàn)1. 電力系統(tǒng)穩(wěn)態(tài)分析 ，陳珩，中國(guó)電力出版社， 2007 年，第三版；2. C 語(yǔ)言程序設(shè)計(jì) ，鄭玲，中國(guó)電力出版社， 2009 年，第一版； 附錄（設(shè)計(jì)流程圖、程序、表格、數(shù)據(jù)等）啟動(dòng)輸入原始數(shù)據(jù)1、設(shè)計(jì)流程圖2、程序 數(shù)據(jù)輸入模塊 struct Line int Num,NumI,NumJ; float R,X,B,K;struct Busint Num ;float Volt,Phase,GenP,GenQ,LoadP,LoadQ; int Type;struct S

9、huntint Num,NumI;float G,B;#include stdio.h#include string.h#include math.h#include stdlib.h#define NBUS 4#define NLINE 4/* Global variables */int nL,nSH,nB,nVA; float XNBUS;int L;void main() FILE *fp; int i,j,k,l,h,c;int i1,i2,i3;float d1,d2,d3,d4,d5,d6,r,x,g,b,e;struct Line sLNLINE;struct Bus sBNB

10、US;struct Shunt sSHNBUS;float YGNBUSNBUS,YBNBUSNBUS;i1=i2=i3=0;d1=d2=d3=d4=d5=d6=0.0;e=0.00001;for (i=0;iNBUS;i+)/*Read the input data*/if (fp=fopen( data.txt ,r )=NULL) printf(Can not open the file named data.txt n);exit(0);fscanf(fp, %d,%d,%d,&nB,&nL,&nSH);for (i=0;inB;i+) sBi.Num=sBi.Type=0;sBi.V

11、olt=1.0; sBi.Phase=sBi.GenP=sBi.GenQ=sBi.LoadP=sBi.LoadQ=0.0;fscanf(fp, %d,%f,%f,%f,%f,%f,%f,%d ,&i1,&d1,&d2,&d3,&d4,&d5,&d6,&i2);sBi.Num=i1;sBi.Volt=d1;sBi.Phase=d2;sBi.GenP=d3;sBi.GenQ=d4;sBi.LoadP=d5,sBi. LoadQ=d6;sBi.Type=i2;for (i=0;inL;i+)sLi.Num=sLi.NumI=sLi.NumJ=0; sLi.R=sLi.X=sLi.B=0.0;sLi.

12、K=1.0;fscanf(fp, %2d %3d %3d %f %f %f %f ,&i1,&i2,&i3,&d1,&d2,&d3,&d4); sLi.Num=i1;sLi.NumI=i2;sLi.NumJ=i3;sLi.R=d1;sLi.X=d2;sLi.B=d3;sLi.K=d4;for (i=0;inSH;i+) sSHi.Num=sSHi.NumI=0;sSHi.G=sSHi.B=0.0; fscanf(fp, %2d %3d %f ,&i1,&i2,&d1); sSHi.Num=i1;sSHi.NumI=i2;sSHi.B=d1;if (fp!=NULL) fclose(fp);/*

13、Make Y Matrix*/for (i=1;inB+1;i+) for (j=1;jnB+1;j+)YGij=0.0;YBij=0.0;for (l=0; lnL; l+)i=sLl.NumI;j=sLl.NumJ;r=sLl.R;x=sLl.X;d1=r*r+x*x;g=r/d1;b=-x/d1;if (fabs(sLl.K-1.0)0.000001)YGii=YGii+g/sLl.B/sLl.B;YGjj=YGjj+g;YBii=YBii+b/sLl.B/sLl.B;YBjj=YBjj+b;YGij=YGij-g/sLl.B;YGji=YGji-g/sLl.B;YBij=YBij-b/

14、sLl.B;YBji=YBji-b/sLl.B;/* Check the Y matrix */if (fp=fopen( GGBB.txt ,w )=NULL)printf( Can not open the file named GGBB.txt n);exit(0);fprintf(fp,-Y Matrix-n);for (i=1;inB+1;i+) for (j=1;j0.000001)fprintf(fp, Y(%3d,%-3d)=(%10.5f,%10.5f)n ,i,j,YGij,YBij);if (fp!=NULL) fclose(fp);/* 設(shè)定電壓初值 */ float

15、VNBUSNBUS=0; for (i=1;iNBUS+1;i+) if (sBi-1.Type=0) Vi0=1.0;Vi1=0.0;for (i=1;iNBUS+1;i+) if (sBi-1.Type=1) Vi0=sBi-1.Volt;Vi1=0.0;for (i=1;iNBUS+1;i+)if (sBi-1.Type=2) Vi0=(sBi-1.Volt)*cos(sBi-1.Phase);/esVi1=(sBi-1.Volt)*sin(sBi-1.Phase);/fs/ 輸出電壓初值if (fp=fopen( 電壓初值 .txt , w )=NULL)printf(Can not

16、open the file named 電壓初值 .txt n );exit(0);fprintf(fp,- 電壓初值 -n );for (i=1;iNBUS+1;i+)for (j=1;j2;j+)fprintf(fp, Y(%2d )=(%10.5f,%10.5f)n ,i,Vi0,Vi1);if (fp!=NULL) fclose(fp);for (c=1;c+)/* 計(jì)算偏移量 */ float PQ1002=0,dV100=0,df100=0;for (i=1;iNBUS;i+)if (sBi-1.Type=0) / 計(jì)算 PQ節(jié)點(diǎn)for (j=1;jNBUS+1;j+)PQi0-=

17、Vi0*(YGij*Vj0-YBij*Vj1)+Vi1*(YGij*Vj1+YBij *Vj0);PQi1-=Vi1*(YGij*Vj0-YBij*Vj1)-Vi0*(YGij*Vj1+YBij *Vj0);if (sBi-1.Type=1) / 計(jì)算 PV節(jié)點(diǎn)for (j=1;jNBUS;j+)PQi0-=Vi0*(YGij*Vj0-YBij*Vj1)+Vi1*(YGij*Vj1+YBij *Vj0); dVi=Vi1*Vi1+Vi0*Vi0;h=1;for (i=1;iNBUS;i+)if (sBi-1.Type=0)dfh=PQi0+sBi-1.GenP-sBi-1.LoadP;h+;d

18、fh=PQi1+sBi-1.GenQ-sBi-1.LoadQ;h+;if (sBi-1.Type=1)dfh=PQi0+sBi-1.GenP-sBi-1.LoadP;h+;dfh=sBi-1.Volt*sBi-1.Volt-dVi;h+; /* 輸出偏移量 */ if (c=1)if (fp=fopen( 第一次迭代偏移量 .txt , w )=NULL)printf( Can not open the file named 第一次迭代偏移量 .txt n );exit(0);fprintf(fp,- 第一次迭代偏移量 -n );for (i=1;ih;i+)fprintf(fp, df(%2

19、d )=(%10.5f)n ,i,dfi);if (fp!=NULL) fclose(fp);/* 計(jì)算雅克比矩陣 */float ga2*(NBUS-1)=0,gb2*(NBUS-1)=0;floatH2*(NBUS-1)2*(NBUS-1)=0,J2*(NBUS-1)2*(NBUS-1)=0,L2*(NBUS-1)2*(NBUS-1)=0,N 2*(NBUS-1)2*(NBUS-1)=0,JJ2*(NBUS-1)2*(NBUS-1)=0,R2*(NBUS-1)2*(NBUS-1)=0,S 2*(NBUS-1)2*(NBUS-1)=0;for (i=1;iNBUS;i+)for (j=1;j

20、NBUS;j+)if (i=j)gai=-(PQi0*Vi0+PQi1*Vi1)/(Vi0*Vi0+Vi1*Vi1);/ 計(jì)算 agbi=-(PQi0*Vi1-PQi1*Vi0)/(Vi0*Vi0+Vi1*Vi1);/ 計(jì)算 bHii=-Vi0*YBii+YGii*Vi1+gbi;Nii=Vi0*YGii+YBii*Vi1+gai;Jii=-Vi0*YGii-YBii*Vi1+gai;Lii=Vi1*YGii-YBii*Vi0-gbi;Rii=2*Vi1;Sii=2*Vi0;elseHij=-Vi0*YBij+YGij*Vi1; / 非對(duì)角線元素正確。Nij=Vi0*YGij+YBij*Vi1

21、;Jij=-Vi0*YGij-YBij*Vi1;Lij=-Vi0*YBij+YGij*Vi1;Rij=0.0;Sij=0.0;h=1;for (i=1;iNBUS;i+)l=1;if (sBi-1.Type=0)for (j=1;jNBUS;j+)JJhl=Hij; l+;JJhl=Nij; l+;h+;l=1;for (j=1;jNBUS;j+)JJhl=Jij;l+;JJhl=Lij;l+;h+;l=1;if (sBi-1.Type=1)for (j=1;jNBUS;j+)JJhl=Hij; l+;JJhl=Nij; l+;h+;l=1;for (j=1;jNBUS;j+)JJhl=Rij

22、;l+;JJhl=Sij;l+;h+;/*if輸出雅克比矩陣(c=1) if (fp=fopen( printf( Can not open the file named exit(0);fprintf(fp, fprintf(fp, for (i=1;ih;i+) fprintf(fp, fprintf(fp, n for (i=1;ih;i+)*/第一次迭代雅克比矩陣 .txt , w )=NULL)雅克比矩陣 .txt n );第一次迭代雅克比矩陣 );n ,c); %dtt,i);); ,i);%2dfprintf(fp,for (j=1;jh;j+) fprintf(fp, %10.

23、5ft ,JJij); fprintf(fp, n );if (fp!=NULL) fclose(fp);/* 高斯法求解方程組 */ float U100=; int LH100=; float t,d;l=1;for (i=1;ih;i+)LHi=0;for (k=1;kh;k+)d=0.0;for (j=k;jd)d=fabs(JJkj);/ 在一行中找到一個(gè)最大值賦值 d,并用JSK 記住這個(gè)最大值所在的列號(hào)LHk=j;if(fabs(d)0.000001)/ 如果 d的數(shù)值太小，做為被除數(shù)將帶來(lái)很大的誤差l=0; elseif (LHk!=k)for (i=1;ih;i+)t=JJi

24、k;JJik=JJiLHk; JJiLHk=t;/ 進(jìn)行列交換，讓最大值始終在對(duì)角元上if(l=0)break ; for(j=k+1;jh;j+)JJkj=JJkj/JJkk;/ 對(duì)角元上的元素消為dfk=dfk/JJkk;for (i=k+1;ih;i+)for (j=k+1;j0.00001)/ 用追趕法，解方程組，求未知數(shù) xUh-1=dfh-1;for (i=h-2;i=0;i-) t=0.0;for (j=i+1;jh;j+)t=t+JJij*Uj;Ui=(dfi-t);/* 輸出高斯結(jié)果 */if (c=1)if (fp=fopen( 第一次迭代高斯 .txt ,w )=NULL

25、) printf( err );exit(0);for (i=1;ih;i+)fprintf(fp,%f ,Ui);fprintf(fp,n ); fclose(fp);/* 得到電壓值 */ h=1;for (i=1;iNBUS;i+)if (sBi-1.Type=0) Vi1+=Uh; h+; Vi0+=Uh; h+;if (sBi-1.Type=1) Vi1+=Uh; h+; Vi0+=Uh; h+; /* 輸出電壓值 */if (c=1)if (fp=fopen( 第一次迭代電壓值 .txt , w )=NULL) 電壓值 .txt n );printf( Can not open t

26、he file namedexit(0);fprintf(fp,- 第一次迭代后的電壓初值 -n );for (i=1;iNBUS+1;i+)for (j=1;j2;j+)fprintf(fp, Y(%2d )=(%10.5f,%10.5f)n ,i,Vi0,Vi1); if (fp!=NULL) fclose(fp);/* 求最大變化值 */float w;w=0;for (i=1;ifabs(w)w=Ui;if (w0.00001)break ;/* 輸出電壓終值 */if (fp=fopen( 7電壓終值 .txt ,w )=NULL)printf(Can not open the fi

27、le named 電壓終值 .txt n );exit(0);fprintf(fp,- 電壓終值 -n );fprintf(fp, 循環(huán) %d次 n ,c);for (i=1;iNBUS+1;i+)for (j=1;j2;j+)fprintf(fp, Y(%2d )=(%10.5f,%10.5f)n ,i,Vi0,Vi1);if (fp!=NULL) fclose(fp);/ 計(jì)算平衡節(jié)點(diǎn)的功率float dp,ps,qs,dsp,dsq,sumgen;ps=0;qs=0;for (i=1;iNBUS+1;i+)if (sBi-1.Type=2)for (j=1;jNBUS+1;j+)ps+=

28、Vi0*(YGij*Vj0-YBij*Vj1)+Vi1*(YGij*Vj1+YBij*Vj 0); / 有功qs+=Vi1*(YGij*Vj0-YBij*Vj1)-Vi0*(YGij*Vj1+YBij*Vj0); / 無(wú)功/ 各節(jié)點(diǎn)功率float g1,g2,b1,b2;float sP100100=,sQ100100=;for (l=0; lnL; l+)i=sLl.NumI;j=sLl.NumJ;if (fabs(sLl.K-1.0)0.000001) /*Normal lines or transformers*/sPij=Vi1*(-Vi1*sLl.B+(Vi1-Vj1)*g-(Vi0

29、-Vj0)*b)+Vi0*(V i0*sLl.B+(Vi1-Vj1)*b+(Vi0-Vj0)*g);sQij=Vi0*(-Vi1*sLl.B+(Vi1-Vj1)*g-(Vi0-Vj0)*b)-Vi1*(Vi0 *sLl.B+(Vi1-Vj1)*b+(Vi0-Vj0)*g);sPji=Vj1*(-Vj1*sLl.B+(Vj1-Vi1)*g-(Vj0-Vi0)*b)+Vj0*(Vi0 *sLl.B+(Vj1-Vi1)*b+(Vj0-Vi0)*g);sQji=Vj0*(-Vj1*sLl.B+(Vj1-Vi1)*g-(Vj0-Vi0)*b)-Vj1*(Vi0 *sLl.B+(Vj1-Vi1)*b+(V

30、j0-Vi0)*g);else /*abnormal transformer ratio*/k=sLl.B;/* 簡(jiǎn)化公式所設(shè)中間變量 */g1=g*(1-k)/k/k; b1=b*(1-k)/k/k; g2=g*(k-1)/k; b2=b*(k-1)/k; g/=k;b/=k;sPij=Vi0*(g*(Vi0-Vj0)-b*(Vi1-Vj1)+g1*Vi0-b1*Vi1)+Vi1*(g*(Vi1-Vj1)+b*(Vi0-Vj0)+g1*Vi1+b1*Vi0);sQij=Vi1*(g*(Vi0-Vj0)-b*(Vi1-Vj1)+g1*Vi0-b1*Vi1)-Vi0*(g*(Vi1-Vj1)+b

31、*(Vi0-Vj0)+g1*Vi1+b1*Vi0);sPji=Vj0*(g*(Vj0-Vi0)-b*(Vj1-Vi1)+g2*Vj0-b2*Vj1)+Vj1*(g*(Vj1-Vi1)+b*(Vj0-Vi0)+g1*Vj1+b1*Vj0);sQji=Vj1*(g*(Vj0-Vi0)-b*(Vj1-Vi1)+g2*Vj0-b2*Vj1)+Vj0*(g*(Vj1-Vi1)+b*(Vj0-Vi0)+g1*Vj1+b1*Vj0);/* 網(wǎng)絡(luò)總損 */ dsp=ps;dsq=qs; sumgen=ps;for (i=1;iNBUS+1;i+)dsp+=sBi-1.GenP-sBi-1.LoadP; dsq

32、+=sBi-1.GenQ-sBi-1.LoadQ; sumgen+=sBi-1.GenP;dp=dsp/sumgen*100; / 輸電效率 /* 輸出功率情況 */if (fp=fopen( 8功率情況 .txt ,w )=NULL)printf(Can not open the file named 功率情況 .txt n );exit(0);fprintf(fp,- 功率情況 -n );fprintf(fp, 平衡節(jié)點(diǎn)功率 S=%10.5f+j%10.5fn ,ps,qs);for (i=1;iNBUS+1;i+)for (j=1;j0.000001)fprintf(fp, S(%3d,

33、%-3d)=(%10.5f,j%10.5f)n ,i,j,sPij,sQij); fprintf(fp,網(wǎng)損為 %10.5f+j%10.3f, 輸電效率為 %10.3fn ,dsp,dsq,100-dp);if (fp!=NULL) fclose(fp);/* 輸出最終結(jié)果 */if (fp=fopen( 9最終結(jié)果 .txt ,w )=NULL)printf(Can not open the file named 最終結(jié)果 .txt n );exit(0);fprintf(fp,- 最終結(jié)果 -n );fprintf(fp,nnn );fprintf(fp, 一共循環(huán)迭代計(jì)算 %d次n ,c

34、);fprintf(fp,nnn );fprintf(fp,- 電壓終值 -n );for (i=1;iNBUS+1;i+)for (j=1;j2;j+)fprintf(fp, U(%2d )=(%8.5f,%8.5f)n ,i,Vi0,Vi1);fprintf(fp,nnn );fprintf(fp, 平衡節(jié)點(diǎn)功率 S=%8.5f+j%7.5fn ,ps,qs);fprintf(fp,nnn );fprintf(fp,- 各節(jié)點(diǎn)間功率 -n );for (i=1;iNBUS+1;i+)for (j=1;j0.000001)fprintf(fp, S(%2d,%-2d)=(%8.5f,j%8.

35、5f)n ,i,j,sPij,sQij); fprintf(fp, nnn );fprintf(fp, 網(wǎng)損為： %7.5f+j%7.5f 。 nn 網(wǎng)損率為：百分之 %7.5fn ,dsp,dsq,dp);if (fp!=NULL) fclose(fp);3、實(shí)驗(yàn)數(shù)據(jù)與表格(1)、節(jié)點(diǎn)導(dǎo)納陣： -Y Matrix-Y( 1,1 )=( 1.01534, -7.49200)Y( 1,2 )=( -0.56148, 2.30208)Y( 1,3 )=( 0.00000, 3.03030)Y( 1,4 )=( -0.45386, 1.89107)Y( 2,1 )=( -0.56148, 2.302

36、08)Y( 2,2 )=( 1.04225, -4.67651)Y( 2,4 )=( -0.48077, 2.40385)Y( 3,1 )=( 0.00000, 3.03030)Y( 3,3 )=( 0.00000, -2.75482)Y( 4,1 )=( -0.45386, 1.89107)Y( 4,2 )=( -0.48077, 2.40385)Y( 4,4 )=( 0.93463, -4.26159)2）第一次迭代雅克比矩陣12317.621040.992652-1.038037.362973-2.30208-0.5614840.56148-2.302085-3.333330.00000

37、60.000000.00000- 第二次迭代雅克比矩陣 12317.611080.767022-1.365347.248863-2.18607-0.7977740.79777-2.186075-3.333330.4741260.000000.00000- 第三次迭代雅克比矩陣 12317.534620.745742-1.354027.163693-2.14287-0.7884940.78849-2.142875-3.298080.4849560.000000.00000- 第四次迭代雅克比矩陣 12317.533680.744532-1.353077.162683-2.14208-0.7881940.78819-2.142085-3.297720.4859560.000000.00000- 第五次迭代雅克比矩陣 12317.533690.744462-1.353007.162703-2.14208-0.7881840.78818-2.142085-3.297710.4860160.000000.00000456-2.30208-0.56148-3.030300.000000.56148-2.30208-0.00000-3.030304.826121.01821