牛頓拉夫遜潮流計(jì)算[整理版].doc

1、牛頓拉夫遜潮流計(jì)算整理版牛頓拉夫遜潮流計(jì)算/ 整個(gè)程序?yàn)?:/ 牛拉法解潮流程序#include#include#define N 4 /節(jié)點(diǎn)數(shù)#define n_PQ 2 /PQ節(jié)點(diǎn)數(shù)#define n_PV 1 /PV節(jié)點(diǎn)數(shù)#define n_br 5 /void main()串聯(lián)支路數(shù)void disp_matrix(float *disp_p,int disp_m,int disp_n); /矩陣顯示函數(shù)float Us2*N=1.0,0,1.0,0,1.05,0,1.05,0; /電壓初值float PsN=0,-0.5,0.2; /有功初值float QsN=0,-0.3; /無功

2、初值float GNN,BNN; / struct各幾點(diǎn)電導(dǎo)電納/ 阻抗參數(shù)int nl; /左節(jié)點(diǎn)int nr; /右節(jié)點(diǎn)float R; /串聯(lián)電阻值float X; /串聯(lián)電抗值float Bl; /左節(jié)點(diǎn)并聯(lián)電導(dǎo)float Br; /右節(jié)點(diǎn)并聯(lián)電納ydatan_br=1,2,0,0.1880,-0.6815,0.6040,1,3,0.1302,0.2479,0.0129,0.0129,1,4,0.1736,0.3306,0.0172,0.0172,3,4,0.2603,0.4959,0.0259,0.0259,2,2,0,0.05,0,0;float Z2; /Z2=R2+X2各串聯(lián)阻

3、抗值的平方float eN,fN,dfe2*(N-1); /e,f存儲(chǔ)電壓的 x 軸分量和 y 軸分量 ,dfe存儲(chǔ)電壓修正值float mid1N,mid2N,dS2*(N-1); /mid1、 mid2 存儲(chǔ)計(jì)算雅克比行列式對角線元素的中間值 ,dS 存儲(chǔ)PQU的不平衡量float Jacob2*(N-1)2*(N-1),inv_J2*(N-1)2*(N-1);/ 雅克比行列式float dPQU=1.0; /PQU不平衡量最大值int kk=0; /迭代次數(shù)int i,j,k;float t;float Pijn_br; /存儲(chǔ)線路 i-j的有功float Qijn_br; /存儲(chǔ)線路

4、i-j的無功float Pjin_br; /存儲(chǔ)線路 j-i的有功float Qjin_br; /存儲(chǔ)線路 j-i的無功float dPijn_br; /存儲(chǔ)線路 i-j的有功損耗float dQijn_br; /存儲(chǔ)線路 i-j的無功損耗float AA,BB,CC,DD; /存儲(chǔ)線路潮流計(jì)算時(shí)的中間值/ 形成導(dǎo)納矩陣-for(i=0;iN;i+)for(j=0;jN;j+)Gij=0;Bij=0;for(i=0;in_br;i+)if(ydatai.nl!=ydatai.nr)Z2=(ydatai.R)*(ydatai.R)+(ydatai.X)*(ydatai.X);/ 串聯(lián)阻抗等效導(dǎo)納

5、值/ 非對角元素Gydatai.nl-1ydatai.nr-1=(-ydatai.R)/Z2;Bydatai.nl-1ydatai.nr-1=ydatai.X/Z2;Gydatai.nr-1ydatai.nl-1=(-ydatai.R)/Z2;Bydatai.nr-1ydatai.nl-1=ydatai.X/Z2;/ 對角元素Gydatai.nl-1ydatai.nl-1+=ydatai.R/Z2;Gydatai.nr-1ydatai.nr-1+=ydatai.R/Z2;Bydatai.nl-1ydatai.nl-1+=(-ydatai.X/Z2);Bydatai.nr-1ydatai.nr-

6、1+=(-ydatai.X/Z2);/ 并聯(lián)導(dǎo)納等效導(dǎo)納值Bydatai.nl-1ydatai.nl-1+=ydatai.Bl;Bydatai.nr-1ydatai.nr-1+=ydatai.Br;elseGydatai.nl-1ydatai.nr-1+=ydatai.R;Bydatai.nl-1ydatai.nr-1+=ydatai.X;printf(G=n);disp_matrix(*G,N,N);printf(B=n);disp_matrix(*B,N,N);/ 分離 e,ffor(i=0;i0.00001)/ 計(jì)算功率不平衡量for(i=0;iN-1;i+)mid1i=0;mid2i=

7、0;for(j=0;jN;j+)mid1i=mid1i+Gij*ej-Bij*fj;mid2i=mid2i+Gij*fj+Bij*ej;dS2*i=Psi-(ei*mid1i+fi*mid2i);if(in_PQ)dS2*i+1=Qsi-(fi*mid1i-ei*mid2i);elsedS2*i+1=Us2*i*Us2*i-(ei*ei+fi*fi);dPQU=0;for(i=0;i2*(N-1);i+)if(dSi0&dPQU0&dPQU0.00001)kk+;/ 形成雅克比行列式-for(i=0;i2*(N-1);i+)for(j=0;j2*(N-1);j+)Jacobij=0;for(j

8、=0;jN-1;j+)/ 求 H,N for(i=0;iN-1;i+)if(i!=j)Jacob2*i2*j=Bij*ei-Gij*fi;Jacob2*i2*j+1=-Gij*ei-Bij*fi;elseJacob2*i2*i=Bii*ei-Gii*fi-mid2i;Jacob2*i2*i+1=-Gij*ei-Bij*fi-mid1i;/ 求 J,L for(i=0;in_PQ;i+)if(i!=j)Jacob2*i+12*j=Gij*ei+Bij*fi; Jacob2*i+12*j+1=Bij*ei-Gij*fi;elseJacob2*i+12*i=Gij*ei+Bij*fi-mid1i;

9、Jacob2*i+12*i+1=Bij*ei-Gij*fi+mid2i;/ 求 R,Sfor(i=n_PQ;iN-1;i+)if(i=j)Jacob2*i+12*i=-2*fi;Jacob2*i+12*i+1=-2*ei;/ 雅克比行列式求逆-/ 初始化 inv_JNN for(i=0;i2*(N-1);i+) for(j=0;j2*(N-1);j+)if(i!=j)inv_Jij=0; else inv_Jij=1;/ 將原矩陣化簡為對角陣for(i=0;i2*(N-1);i+)for(j=0;j2*(N-1);j+)if(i!=j)t=Jacobji/Jacobii;for(k=0;k2*

10、(N-1);k+)Jacobjk-=Jacobik*t;inv_Jjk-=inv_Jik*t;/ 原矩陣各對角元素化為 1，畫出逆矩陣for(i=0;i2*(N-1);i+)if(Jacobii!=1)t=Jacobii;for(j=0;j2*(N-1);j+)inv_Jij=inv_Jij/t;/ 求電壓修正值-for(i=0;i2*(N-1);i+)dfei=0;for(j=0;j2*(N-1);j+)dfei-=inv_Jij*dSj;for(i=0;iN-1;i+)ei+=dfe2*i+1;fi+=dfe2*i;elsebreak;/ 循環(huán)結(jié)束-/ 求平衡節(jié)點(diǎn)功率-mid1N-1=0;

11、mid2N-1=0;for(j=0;jN;j+)mid1N-1=mid1N-1+GN-1j*ej-BN-1j*fj;mid2N-1=mid2N-1+GN-1j*fj+BN-1j*ej;PsN-1=eN-1*mid1N-1+fN-1*mid2N-1;QsN-1=fN-1*mid1N-1-eN-1*mid2N-1;for(i=n_PQ;iN-1;i+)Qsi=fi*mid1i-ei*mid2i;/-/ 顯示輸出結(jié)果printf(kk=%dn,kk);printf(P=);for(i=0;iN;i+)printf(%9.4f,Psi);printf(nQ=);for(i=0;iN;i+)printf

12、(%9.4f,Qsi);printf(ne=);for(i=0;iN;i+)printf(%9.4f,ei);printf(nf=);for(i=0;iN;i+)printf(%9.4f,fi);printf(n);/ 求線路上的潮流/ 計(jì)算 Sij for(i=0;in_br;i+)if(ydatai.nl!=ydatai.nr)Z2=(ydatai.R)*(ydatai.R)+(ydatai.X)*(ydatai.X);AA=-fydatai.nl-1*ydatai.Bl+(eydatai.nl-1-eydatai.nr-1)*ydatai.R/Z2+(fydatai.nl-1-fydat

13、ai.nr-1)*ydatai.X/Z2;BB=-eydatai.nl-1*ydatai.Bl-(fydatai.nl-1-fydatai.nr-1)*ydatai.R/Z2+(eydatai.nl-1-eydatai.nr-1)*ydatai.X/Z2;Piji=eydatai.nl-1*AA-fydatai.nl-1*BB; Qiji=eydatai.nl-1*BB+fydatai.nl-1*AA; printf(S%d%d=%9.4f+j%9.4fn,ydatai.nl,ydatai.nr,Piji,Qiji);dPiji=Piji+Pjii;dQiji=Qiji+Qjii;printf

14、(n);/ 計(jì)算 Sji for(i=0;in_br;i+)if(ydatai.nl!=ydatai.nr)Z2=(ydatai.R)*(ydatai.R)+(ydatai.X)*(ydatai.X);CC=-fydatai.nr-1*ydatai.Br+(eydatai.nr-1-eydatai.nl-1)*ydatai.R/Z2+(fydatai.nr-1-fydatai.nl-1)*ydatai.X/Z2; DD=-eydatai.nr-1*ydatai.Br-(fydatai.nr-1-fydatai.nl-1)*ydatai.R/Z2+(eydatai.nr-1-eydatai.nl-1)*ydatai.X/Z2;Pjii=eydatai.nr-1*CC-fydatai.nr-1*DD;Qjii=eydatai.nr-1*DD+fydatai.nr-1*CC; printf(S%d%d=%9.4f+j%9.4fn,ydatai.nr,ydatai.nl,Pjii,Qjii);printf(n);/ 計(jì)算 dSij for(i=0;in_br;i+)if(ydatai.nl!=ydatai.nr)dPiji=Piji+Pjii;dQiji=Qiji+Qjii;printf(