電力系統(tǒng)分析課程設(shè)計 潮流計算

1、武漢理工大學(xué)電力系統(tǒng)分析課程設(shè)計說明書目錄摘要11.任務(wù)及題目要求22.設(shè)計原理2 2.1設(shè)計思路 .2 2.2方法原理.33計算步驟64.結(jié)果分析8小結(jié)10參考文獻11附錄：源程序12本科生課程設(shè)計成績評定表24摘要潮流計算，就是根據(jù)給定電力系統(tǒng)的接線方式和參數(shù)等運行條件，確定電力系統(tǒng)穩(wěn)態(tài)運行各母線電壓、個支路電流與功率及網(wǎng)損。其是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運算。本設(shè)計即根據(jù)給定3條支路的運行條件，運用牛頓-拉夫遜的潮流計算方法，求出各線路的功率分布及功率損耗，未知節(jié)點的電壓。同時進行編程，在matlab運行得到結(jié)果。 關(guān)鍵詞：潮流計算 電力系統(tǒng) 牛頓-拉夫遜法1.任務(wù)及

2、題目要求對如下所訴系統(tǒng)編程進行潮流計算：節(jié)點數(shù)：3 支路數(shù)：3 計算精度：0.00010支路1： 0.0300+j0.0900 12支路2： 0.0200+j0.0900 23支路3： 0.0300+j0.0900 31節(jié)點1：PQ節(jié)點，S（1）=-0.5000-j0.2000節(jié)點2：PQ節(jié)點，S（2）=-0.6000-j0.2500節(jié)點3：平衡節(jié)點，U（3）=1.00000.0000經(jīng)分析可知題目所給的系統(tǒng)為三節(jié)點組成的環(huán)形回路，且均為線路沒有變壓器。有兩個為PQ節(jié)點，這類節(jié)點的有功功率P和無功功率Q是給定的，節(jié)點電壓是待求量。節(jié)點3為平衡節(jié)點，在潮流分布算出來以前，網(wǎng)絡(luò)中的功率損失是未知的

3、，因此，至少要有一個節(jié)點的有功功率不能給定，這個節(jié)點承擔(dān)了系統(tǒng)的有功功率平衡，故稱之為平衡節(jié)點。對系統(tǒng)進行潮流計算，要求各線路的功率分布及功率損耗，未知節(jié)點的電壓等，并進行編程運行得到結(jié)果。2. 設(shè)計原理2.1 設(shè)計思路本設(shè)計給出三條支路阻抗，和相關(guān)節(jié)點數(shù)據(jù) 節(jié)點1：PQ節(jié)點，S（1）=-0.5000-j0.2000 節(jié)點2：PQ節(jié)點，S（2）=-0.6000-j0.2500 節(jié)點3：平衡節(jié)點，U（3）=1.00000.0000進行潮流計算，求各線路的功率分布及功率損耗，未知節(jié)點的電壓。常用的計算潮流分析的方法有牛頓拉夫遜法，PQ分解法等。本次設(shè)計采用牛頓拉夫遜法進行計算。牛頓拉夫遜法是一種求

4、解非線性方程的數(shù)值解法 ，其要點是把非線性方程式的求解過程變成反復(fù)地對相應(yīng)的線性方程式進行求解的過程，即通常所稱的逐次線性化過程。由于便于編寫程序用計算機求解，應(yīng)用較廣。 2.1使用方法簡介牛頓拉夫遜法（NewtonRaphson法）是求解非線性方程代數(shù)方程組的有效迭代計算方法，同時是一種由泰勒展開式只取線性項所得到的線性近似。在牛頓拉夫遜法的每一次迭代過程中，對非線性方程通過線性化處理逐步近似。下面以單變量加以說明。 設(shè)欲求解的非線性代數(shù)方程為 (2-1) 設(shè)方程的真實解為x*，則必有f(x*)=0。用牛頓-拉夫遜法求方程真實解x*的步驟如下： 先給出解的近似值它與真解的誤差為，取x

5、5;=x°+x°為第一次的修正估值 ，即 (2-2)其中x°為初始估值的增量，即x°=x¹-x°。設(shè)函數(shù)f(x)具有任意階導(dǎo)數(shù)，即可將上式在x°的鄰域展開為泰勒級數(shù)，即： （2-3）式中,分別為函數(shù)在處的一階導(dǎo)數(shù)，.，n階導(dǎo)數(shù)。如果差值很小，2-3式右端的一切高階項均可略去。于是，2-3便簡化為 0 (2-4) 這是對于變量的修正量的方程式，。解此方程可得修正量 (2-5)用所求的去修正近似解，變得 (2-6)可見，只要f'(x°)0，即可根據(jù)上式求出第一次的修正估值x¹，若恰巧有f(x¹

6、;)=0，則方程的真實解即為x*=x¹。若f(x¹)0，則用上述方法由x¹再確定第二次的修正估值x²， 這樣的迭代計算可以反復(fù)進行下去，迭代計算的通式是 (2-7)迭代過程的收斂判據(jù)為 (2-8)或 (2-9) 式中，為預(yù)先給定的小正數(shù)。 應(yīng)用牛頓法求解多變量非線性方程組2-1時，假定已給出各變量的初值，. ，令，. 分別為各變量的修正量，使其滿足方程 (2-10) 將上式中的n個多元函數(shù)在初始值附近展成泰勒級數(shù),并略去含有，的高階項，便得到. (2-11) 方程式2-11是對于修正量， 的線性方程組,稱為牛頓法的修正方程式.利用高斯消去法或三角分解法可

7、以解出修正量，。然后對初始近似值進行修正 (i=1,2,.,n) (2-13)如此反復(fù)迭代，在進行k1次迭代時，從求解修正方程式 (2-14)得到修正量，并對各變量進行修正得到： (2-16) 和 (2-17) 式中的X和分別是由n個變量和修正量組成的n維列向量；F(X)是由n個多元函數(shù)組成的n維列項量；J是n階方陣，稱為雅可比矩陣，它的第i、j個元素是第n個函數(shù)對第j個變量的偏導(dǎo)數(shù)；上角標(biāo)(k)表示陣的每一個元素都在點處取值。迭代過程一直到滿足收斂判據(jù) (2-18) 或 (2-19)為止。和為預(yù)先給定的小正數(shù)。需要指出的是，用牛頓-拉夫遜法求解非線性代數(shù)方程時，巧妙地設(shè)定一個合適的初始值是十

8、分重要的，若初始值設(shè)定的合適，不但可以保證迭代收斂（即向真實值逼近），而且可使迭代次數(shù)減少。相反，若初始值設(shè)定的不合適，則不但要增加迭代次數(shù)，而且還可能迭代發(fā)散（即遠(yuǎn)離真實值）或者循環(huán)（即在真實值附近往復(fù)變化）。為解決此問題，可在疊代前先粗略畫出函數(shù)f(x)是曲線，在該曲線與x軸的焦點附近設(shè)定合適的初始估值會有裨益。 3計算步驟用牛頓拉夫遜法計算電力系統(tǒng)潮流的基本步驟：(1)形成節(jié)點導(dǎo)納矩陣；(2)將各節(jié)點電壓設(shè)初值 ；(3)將節(jié)點初值代入相關(guān)求式，求出修正方程式的常數(shù)項向量；(4)將節(jié)點電壓初值代入求式，求出雅可比矩陣元素；(5)求解修正方程，求修正向量；(6)求取節(jié)點電壓的新值；(7)檢查

9、是否收斂，如不收斂，則以各節(jié)點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉(zhuǎn)入下一步；(8)計算支路功率分布，節(jié)點無功功率和平衡節(jié)點功率。具體流程圖見下頁圖二。啟動輸入原始數(shù)據(jù)形成節(jié)點導(dǎo)納陣給定電壓初值e（0）、f（0）k=0根據(jù)公式計算用公式計算雅克比矩陣各元素高斯法解修正方程，求是計算平衡節(jié)點功率及全部線路功率輸出k=k+1否 圖2計算流程圖4.結(jié)果分析 經(jīng)過用C語言編程運算得到第一次迭代結(jié)果為： 節(jié)點導(dǎo)納矩陣: 6.66667+j-20.00000-3.33333+j10.00000-3.33333+j10.00000-3.33333+j10.00000 5.68627+j-2

10、0.58824-2.35294+j10.58824-3.33333+j10.00000-2.35294+j10.58824 5.68627+j-20.58824 Gauss_Seidel迭代結(jié)果： 第1次迭代PQ節(jié)點e及f的值：0.983500 -0.019500 第1次迭代PQ節(jié)點e及f的值：0.973610 -0.033997 第2次迭代PQ節(jié)點e及f的值：0.969642 -0.036485 第2次迭代PQ節(jié)點e及f的值：0.965875 -0.042673 Jacobi 矩陣: 19.31452 6.80237-9.57480-3.59699-7.5856018.98468 3.596

11、99-9.57480-9.51651-3.6463119.92097 5.76424 3.64631-9.51651-6.9773519.36506 不平衡值：-0.114256-0.054375-0.002280-0.007412 輸出 df,de 節(jié)點為 1 的電壓不平衡值： df=-0.00508 de=-0.00267 節(jié)點為 2 的電壓不平衡值： df=-0.00006 de=-0.00037 輸出迭代過程中的電壓值 節(jié)點為 1 的電壓修正新值： f=-0.04156 e= 0.96698 節(jié)點為 2 的電壓修正新值： f=-0.04273 e= 0.96551 還不能達到題目中的精

12、度要求，繼續(xù)進行迭代，經(jīng)6次迭代后得到結(jié)果。 各節(jié)點電壓為: U1=0.963186 + j-0.042323 U2=0.962333 + j-0.045732 U3=1.000000 + j0.000000 平衡節(jié)點功率為: S3=1.118788+j0.518288 輸電線路功率如下: 線路1-2的功率: 0.035453 + j-0.004289 線路2-1的功率: -0.035412 + j0.004412 線路2-3的功率: -0.564587 + j-0.254058 線路3-2的功率: 0.572847 + j0.291225 線路3-1的功率: 0.545942 + j0.22

13、7063 線路1-3的功率: -0.535453 + j-0.195598 輸電線路上損耗的功率為: 線路1-2上的功率損耗: 0.000041 + j 0.000123 線路2-3上的功率損耗: 0.008259 + j 0.037167 線路3-1上的功率損耗: 0.010488 + j 0.031465 網(wǎng)絡(luò)總損耗: 0.018789 + j 0.068755 滿足精度要求，計算出了潮流計算要求的各節(jié)點電壓，網(wǎng)絡(luò)中的功率分布及功率損耗。完成和達到了設(shè)計的要求和目的。小結(jié)在本次電力系統(tǒng)分析課程實習(xí)中，我感覺到考試僅僅體現(xiàn)了學(xué)習(xí)能力的一部分，而課程設(shè)計更能體現(xiàn)了綜合的學(xué)習(xí)實踐能力。原本模糊不

14、清的基本概念和實用方法，在操作本次課設(shè)的過程中都得到了明確和加強。課程設(shè)計讓我感受到比考試更大的壓力，有了更高的效率和收獲，深層次加深了知識的深度層次。望進行更多一些這樣的實踐性自主行的學(xué)習(xí)活動，受益良多。在本次課程設(shè)計得到了很多老師和同學(xué)的支持和幫助，特在此表示感謝，不再一一署名。參考文獻：1 何仰贊 溫增銀 電力系統(tǒng)分析 華中科技大學(xué)出版社 20062 于永源 楊綺雯 電力系統(tǒng)分析（第二版）中國電力出版社 20073 大電網(wǎng)最優(yōu)潮流計算 劉明波、謝敏、 趙維興 科學(xué)出版社 20104 電力系統(tǒng)分析的計算機算法 邱曉燕、 劉天琪 中國電力出版社 20095 薛定宇 控制系統(tǒng)計算機輔助設(shè)計 清

15、華大學(xué)出版社 1996附錄：源程序#include <stdio.h>#include <math.h>#define M 40 /*矩陣階數(shù)*/#define N 10 /*最大迭代次數(shù)*/ int i,j,k,a,b,l,z; /* 循環(huán)變量 */ int n, /* 節(jié)點數(shù) */ m, /* 支路數(shù) */ dd, /*對地支路數(shù)*/ pq, /* PQ節(jié)點數(shù) */ pv, /* PV節(jié)點數(shù) */ byq; /* 變壓器數(shù) */ float eps, /* 精度 */ bb, /* 變比 */ max; /*指標(biāo)標(biāo)志*/ float yyM; /*中間變量*/ st

16、atic float GMM,BMM,B1MM,B2MM,B3MM; struct jiedian /* 定義節(jié)點結(jié)構(gòu)體 */ int num,s; /* num為節(jié)點號，s為節(jié)點類型*/ float p,q,u,f,v; float dp,dq,df,du; jiedianM; struct zhilu /* 定義支路結(jié)構(gòu)體 */ int num; int p1,p2; /*支路的兩個節(jié)點*/ float r,x; /*支路的電阻與電抗*/ zhiluM; FILE *fp1,*fp2; /*文件指針*/ void data1() /* 讀取數(shù)據(jù)函數(shù) */ int h,numb,a,b; a

17、=1; fp1=fopen("d:input.txt","r"); if(fp1=NULL) printf(" can not open file !n"); exit(0); fscanf(fp1,"%d,%d,%d,%d,%d,%d,%fn",&n,&m,&dd,&pq,&pv,&byq,&eps); /*輸入節(jié)點數(shù),支路數(shù),對地支路數(shù),PQ節(jié)點數(shù)，PV節(jié)點數(shù)，變壓器數(shù)和精度*/ j=1;k=pq+1; for(i=1;i<=n;i+) /*輸入節(jié)點

18、類型的輸入功率和節(jié)電電壓初值*/ fscanf(fp1,"%d,%d",&numb,&h);if(h=1) /*類型h=1是PQ節(jié)點*/ fscanf(fp1,",%f,%f,%f,%fn",&jiedianj.p,&jiedianj.q,&jiedianj.u,&jiedianj.f); /*輸入PQ節(jié)點的數(shù)據(jù)*/ jiedianj.num=numb; jiedianj.s=h; j+; if(h=2) /*類型h=2是pv節(jié)點*/ fscanf(fp1,",%f,%f,%f,%fn"

19、,&jiediank.p,&jiediank.v,&jiediank.u,&jiediank.f); /*輸入PV節(jié)點的數(shù)據(jù)*/ jiediank.num=numb; jiediank.s=h; k+; if(h=3) /*類型h=3是平衡節(jié)點*/ fscanf(fp1,",%f,%fn",&jiediann.u,&jiediann.f); /*輸入平衡節(jié)點的數(shù)據(jù)*/ jiediann.num=numb; jiediann.s=h; for(i=1;i<=m;i+) /*輸入支路阻抗*/ fscanf(fp1,"

20、;%d,%d,",&zhilui.num,&h);if(h=0) /*輸入非變壓器支路阻抗*/fscanf(fp1,"%d,%d,%f,%fn", &zhilui.p1,&zhilui.p2,&zhilui.r,&zhilui.x); if(h=1) /*輸入變壓器支路阻抗*/ fscanf(fp1,"%f,%d,%d,%f,%fn",&bb,&zhilui.p1,&zhilui.p2,&zhilui.r, &zhilui.x); zhilui.r=zhil

21、ui.r*bb; zhilui.x=zhilui.x*bb; b=m+a; zhilub.num=b; zhilub.p1=zhilui.p1; zhilub.p2=0; zhilub.r=zhilui.r*bb/(1-bb); zhilub.x=zhilui.x*bb/(1-bb); zhilub+1.num=b zhilub+1.p1=zhilui.p2; zhilub+1.p2=0; zhilub+1.r=zhilui.r/(bb-1); zhilub+1.x=zhilui.x/(bb-1); a=a+2; fclose(fp1); if(fp2=fopen("d:output

22、.txt","w")=NULL) printf(" can not open file!n"); exit(0); fprintf(fp2,"nn * 原始數(shù)據(jù) *n"); fprintf(fp2," =n") fprintf(fp2," 節(jié)點數(shù):%2d 支路數(shù):%2d 對地支路數(shù):%2d 變壓器數(shù):%2d n PQ節(jié)點數(shù):%2d PV節(jié)點數(shù):%2d 精度:%fn",n,m+2*byq,dd,byq,pq,pv,eps); fprintf(fp2," -n"); f

23、or(i=1;i<=pq;i+) fprintf(fp2," PQ節(jié)點 節(jié)點%2d P%d=%+f Q%d=%+fn", jiediani.num,jiediani.num,jiediani.p,jiediani.num,jiediani.q); for(i=pq+1;i<=pq+pv;i+) fprintf(fp2," PV節(jié)點 節(jié)點%2d P%d=%+f V%d=%+fn", jiediani.num,jiediani.num,jiediani.p,jiediani.num,jiediani.v); fprintf(fp2,"

24、平衡節(jié)點 節(jié)點%2d u%d=%+f f%d=%+fn", jiediann.num,jiediann.num,jiediann.u,jiediann.num,jiediann.f); fprintf(fp2," -n"); for(i=1;i<=m+2*byq;i+) fprintf(fp2," 支路%2d 相關(guān)節(jié)點:%2d,%2d R=%+f X=%+fn", i,zhilui.p1,zhilui.p2,zhilui.r,zhilui.x); fprintf(fp2," =n"); void form_y() /*

25、 形成節(jié)點導(dǎo)納矩陣函數(shù) */ float S; for(i=1;i<=m+2*byq;i+) /*節(jié)點導(dǎo)納矩陣的主對角線上的導(dǎo)納*/ for(j=0;j<=n;j+) if(zhilui.p1=j)|(zhilui.p2=j) S=zhilui.r*zhilui.r+zhilui.x*zhilui.x; if(S=0) continue; Gjj+=zhilui.r/S; Bjj+=-zhilui.x/S; for(i=1;i<=m+2*byq;i+) /*對地導(dǎo)納*/ j=zhilui.p1; k=zhilui.p2; if(j=0)|(k=0) S=zhilui.r*zh

26、ilui.r+zhilui.x*zhilui.x; if(S=0) G0i=B0i=Gi0=Bi0=0;continue; if(j=0) G0k+=zhilui.r/S; B0k+=-zhilui.x/S; Gk0=G0k; Bk0=B0k; if(k=0) G0j+=zhilui.r/S; B0j+=-zhilui.x/S; Gj0=G0j; Bj0=B0j; for(k=1;k<=m-dd;k+) /*節(jié)點導(dǎo)納矩陣非主對角線上的導(dǎo)納*/ i=zhiluk.p1; j=zhiluk.p2;S=zhiluk.r*zhiluk.r+zhiluk.x*zhiluk.x;if(S=0) co

27、ntinue;Gij+=-zhiluk.r/S;Bij+=zhiluk.x/S;Gji=Gij;Bji=Bij; for(i=1;i<=(pq+pv);i+) /*形成B'矩陣*/ for(j=1;j<=(pq+pv);j+) B1ij=Bjiediani.numjiedianj.num; for(i=1;i<=pq;i+) /*形成B''矩陣*/ for(j=1;j<=pq;j+) B2ij=Bjiediani.numjiedianj.num; /*輸出節(jié)點導(dǎo)納矩陣、系數(shù)矩陣B'、系數(shù)矩陣B''*/ fprintf(fp

28、2,"nn * 計算結(jié)果 *n"); fprintf(fp2,"=n") fprintf(fp2,"n 節(jié)點導(dǎo)納矩陣為:"); for(i=1;i<=n;i+) fprintf(fp2,"n "); for(j=1;j<=n;j+) fprintf(fp2,"%+8.6f+j%+8.6f ",Gij,Bij); fprintf(fp2,"n -n"); fprintf(fp2,"n 節(jié)點導(dǎo)納矩陣實部G為:"); for(i=1;i<=n;

29、i+) fprintf(fp2,"n "); for(j=1;j<=n;j+) fprintf(fp2,"%+8.6f ",Gij); fprintf(fp2,"n -n"); fprintf(fp2,"n 節(jié)點導(dǎo)納矩陣虛部B為:"); for(i=1;i<=n;i+) fprintf(fp2,"n "); for(j=1;j<=n;j+) fprintf(fp2,"%+8.6f ",Bij); fprintf(fp2,"n -n");

30、fprintf(fp2,"n 系數(shù)矩陣B'為:"); for(i=1;i<=pv+pq;i+) fprintf(fp2,"n "); for(j=1;j<=pv+pq;j+) fprintf(fp2,"%+8.6f ",B1ij); fprintf(fp2,"n -n"); fprintf(fp2,"n 系數(shù)矩陣B''為:"); for(i=1;i<=pq;i+) fprintf(fp2,"n "); for(j=1;j<=pq

31、;j+) fprintf(fp2,"%+8.6f ",B2ij); fprintf(fp2,"n =n"); void gaoshi(int A,float aaMM,float bbM) /*高斯計算函數(shù)*/int i,j,z; int k=0; float a; for(i=k+1;i<A;i+) for(z=0;(z+i)<A;z+) aai+zk=aai+zk/aakk; /*校驗過程*/ bbi+z=bbi+z-aai+zk*bbk; for(j=k+1;j<A;j+) aai+zj=aai+zj-aai+zk*aakj; k

32、=k+1; bbA-1=bbA-1/aaA-1A-1; for(i=A-2,k=A-2;i>=0;i-) a=0; for(j=i+1;j<A;j+) a=aaij*bbj+a; bbi=(bbi-a)/aaii; float solve(float max) /*迭代函數(shù)*/ float aaM,bbM;/*求節(jié)點有功功率不平衡量dp*/ for(i=1;i<=(pq+pv);i+) /*求dp*/ jiediani.dp=jiediani.p; for(j=1;j<=n;j+) jiediani.dp=jiediani.dp-jiediani.u*jiedianj.

33、u*(Gjiediani.numjiedianj.num*cos(jiediani.f-jiedianj.f)+Bjiediani.numjiedianj.num*sin(jiediani.f-jiedianj.f); for(i=1;i<=(pq+pv);i+) /*定義各節(jié)點的dP/U*/ if(jiediani.s!=3) yyi=fabs(jiediani.dp/jiediani.u); max=yy1; /*求dP/U的最大值*/ for(i=1;i<=(pq+pv);i+) if(yyi>max)max=yyi; for(i=1;i<=(pq+pv);i+)

34、 fprintf(fp2," 有功功率不平衡量dP%d=%+fn",jiediani.num,jiediani.dp); fprintf(fp2,"n -n"); /*求各節(jié)點的相位角f，df*/ a=1; /*裝換為高斯函數(shù)的標(biāo)準(zhǔn)形式*/ for(i=1;i<=n;i+) if(jiediani.s=3)a=a+1;continue; aai-a=jiediani.dp/jiediani.u; b=1; for(j=1;j<=n;j+) if(jiedianj.s=3)b=b+1;continue; B3i-aj-b=0.0-B1i+1-a

35、j+1-b*jiediani.u; gaoshi(pv+pq),B3,aa); /*代入高斯函數(shù)求解方程*/ a=1; for(i=1;i<=n;i+) if(jiediani.s=3)a+;continue; jiediani.df=aai-a; jiediani.f=jiediani.f+jiediani.df; for(i=1;i<=(pq+pv);i+) fprintf(fp2," df%d=%f 電壓的相位角f%d=%+fn",jiediani.num,jiediani.df,jiediani.num,jiediani.f); fprintf(fp2,

36、"n -n"); /*求節(jié)點的無功功率不平衡量dq*/ for(i=1;i<=n;i+) /*求dq*/ jiediani.dq=jiediani.q; for(j=1;j<=n;j+) jiediani.dq=jiediani.dq-jiediani.u*jiedianj.u*(Gjiediani.numjiedianj.num*sin(jiediani.f-jiedianj.f)-Bjiediani.numjiedianj.num*cos(jiediani.f-jiedianj.f); for(i=1;i<=pq;i+) /*定義各節(jié)點的dQ/U*/

37、yyi=fabs(jiediani.dq/jiediani.u); for(i=1;i<=pq;i+) /*求所有dP/U和dQ/U中的最大值*/ if(yyi>max)max=yyi; for(i=1;i<=pq;i+) fprintf(fp2," 無功功率不平衡量dQ%d=%+fn",jiediani.num,jiediani.dq); fprintf(fp2,"n -n"); /*計算各節(jié)點電壓*/ for(i=1;i<=n;i+) /*裝換為高斯函數(shù)的標(biāo)準(zhǔn)形式*/ for(j=1;j<=n;j+) B3i-1j-1=

38、0-B2ij; a=1; for(i=1;i<=n;i+) if(jiediani.s!=1)a=a+1;continue; bbi-a=jiediani.dq/jiediani.u; gaoshi(pq,B3,bb); /*代入高斯函數(shù)求解方程*/ a=1; for(i=1;i<=n;i+) if(jiediani.s=3)a+;continue; jiediani.du=bbi-a; jiediani.u=jiediani.u+jiediani.du; for(i=1;i<=n;i+) if(jiediani.s=1) fprintf(fp2," dU%d=%f

39、 電壓的大小U%d=%+fn",jiediani.num,jiediani.du,jiediani.num,jiediani.u); fprintf(fp2,"n =n"); return(max); /*返回循環(huán)標(biāo)志值*/float mul_Re(x1,y1,x2,y2) /*求實部*/float x1,x2,y1,y2;float x; x=x1*x2-y1*y2; return(x); float mul_Im(x1,y1,x2,y2) /*求虛部*/float x1,x2,y1,y2; float y; y=x1*y2+x2*y1; return(y);

40、void data2() /* 潮流計算結(jié)果 */ float sp=0,sq=0,mo1,mo2,mo3,mo4,mo5,mo6,ff,ttM; int i1,j1,i11=0,j11=0; static float PMM,QMM; fprintf(fp2,"nn 各節(jié)點電壓為:"); for(i=1;i<=n;i+) ttM=jiediani.u;jiediani.u=(jiediani.u)*(cos(jiediani.f); /*由極坐標(biāo)轉(zhuǎn)化為直角坐標(biāo)*/ jiediani.f=(ttM)*(sin(jiediani.f); fprintf(fp2,"n U%d=%+8.6f + j %+8.6f", jiediani.num,jiediani.u,jiediani.f); /*輸出直角坐標(biāo)節(jié)點電壓*/ fprintf(fp2,"n =n"); fprintf(fp2,"nn 平衡節(jié)點功率為:"); for(i=1;i<=n;i+) i1=jiediann.num; j1=jiediani.num; sp+=mul_Re(Gi1j1,-Bi1j1,jiediani.u,-jie