九節(jié)點(diǎn)系統(tǒng)潮流計(jì)算編程牛N-R法

1、如圖所示系統(tǒng)，試計(jì)算潮流分布，相關(guān)數(shù)據(jù)見版潮流計(jì)算用戶手冊P121。#includesrxA-zvi發(fā)電1<>#include<>floatdivRe(floatb1,floatb2,floatb3,floatb4)(floata1r;a1r=(b1*b3+b2*b4)/(b3*b3+b4*b4);return(a1r);floatdivIm(floatb1,floatb2,floatb3,floatb4)(floata1i;a1i=(b2*b3-b1*b4)/(b3*b3+b4*b4);return(a1i);floatmulRe(floatb1,floatb2,fl

2、oatb3,floatb4)(floata2r;a2r=b1*b3-b2*b4;return(a2r);floatmulIm(floatb1,floatb2,floatb3,floatb4)(floata2i;a2i=b2*b3+b1*b4;return(a2i);floatMax(floata,intn)inti;floatmax;max=fabs(a0);for(i=1;i<n;i+)if(fabs(ai)>max)max=fabs(ai);return(max);voidmain()inti,j,k,h,km;intT=16;floateps,sumpi1,sumpi2,su

3、mqi1,sumqi2,max,sumir,sumii,I1r,I1i,t,xx,xxx;floatpi08,qi08,detpi8,detqi8,Iir08,Iii08,J01616,detsi16,detui16,delta_p99,delta_q99,a1632,ni1616,H88,N88,J88,L88,ei19fi19,sp99,sq99;staticfloatybr99=,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;staticfloatybi99=,0,0,0,0,0

4、,0,0,0,0,0,0,0,0,0,16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,16,0,0,0,-16,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,;staticfloatyd99=0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),(0,0,0,0,0,0,0),0,0,0,0,0,0,0),0,0,0;floatei09=,;floatfi09=;floatpi9=0,0,0,0;floatqi9=0,0,0,0,0,0;h=0;km=15;eps=;doh+

5、=1;printf("nNowThe%dthtimes'n",h);for(i=0;i<8;i+)printf("ei0%d=%付”,i,ei0i);printf("fi0%d=%ft",i,fi0i);for(i=0;i<8;i+)printf("pi%d=%ft",i,pii);printf("qi%d=%ft",i,qii);sumpi2=0;sumqi2=0;for(i=0;i<8;i+)for(j=0;j<9;j+)sumpi1=(ei0i*(ybrij*ei0吐

6、ybiij*fi0j)+fi0i*(ybrij*fi0j+ybiij*ei0j);sumpi2+=sumpi1;pi0i=sumpi2;printf("pi0%d=%ft",i,pi0i);sumpi2=0;for(i=0;i<8;i+)for(j=0；j<9；j+)(sumqi1=(fi0i*(ybrij*ei0j-ybiij*fi0j)-ei0i*(ybrij*fi0j+ybiij*ei0j)；sumqi2+=sumqi1;qi0i=sumqi2;printf("qi0%d=%ft",i,qi0i);sumqi2=0;for(i=0;i&

7、lt;8;i+)detpii=pii-pi0i;detqii=qii-qi0i;if(i=6|i=7)qi0i=ei0i*ei0i+fi0i*fi0i;detqii=i;printf("detpi%d=%ft",i,detpii);printf("detqi%d=%ft",i,detqii);/*/節(jié)點(diǎn)的注入電流表達(dá)式for(i=0;i<8;i+)Iii0i=0;Iir0i=0;for(i=0;i<8;i+)for(j=0;j<9;j+)Iir0i+=ybrij*ei0j-ybiij*fi0j;Iii0i+=ybrij*fi0j+ybi

8、ij*ei0j;)*/求解NHJL矩陣for(i=0;i<8;i+)for(j=0;j<8;j+)if(i=j)if(i=6|i=7)Hij=-ybiij*ei0j+ybrij*fi0j+Iii0i;Nij=ybrij*ei0j+ybiij*fi0j+Iir0i;Jij=2*fi0i;Lij=2*ei0i;)elseHij=-ybiij*ei0j+ybrij*fi0j+Iii0i;Nij=ybrij*ei0j+ybiij*fi0j+Iir0i;Jij=-ybrij*ei0j-ybiii*fi0j+Iir0i;Lij=-ybiij*ei0j+ybrij*fi0j-Iii0i;)els

9、eif(i=6|i=7)Hij=ybrij*fi0j-ybiij*ei0j;Nij=ybrij*ei0j+ybiij*fi0j;Jij=0;Lij=0;)elseHij=ybrij*fi0j-ybiij*ei0j;Nij=ybrij*ei0j+ybiij*fi0j;Jij=-ybiij*fi0j-ybrij*ei0j;Lij=ybrij*fi0j-ybiij*ei0j;)/形成jacobian矩陣for(i=0;i<16;i+)for(j=0;j<16;j+)if(i%2=0&&j%2=0)J0ij=Hi/2j/2;elseif(i%2=0&&j%2

10、!=0)J0ij=Ni/2(j-1)/2;elseif(i%2!=0&&j%2=0)J0ij=J(i-1)/2j/2;elseJ0ij=L(i-1)/2(j-1)/2;)/for(i=0;i<16;i+)/for(j=0;j<16;j+)/printf("J0%d%d=%ft",i,j,J0ij);/*/求detui/*for(i=0;i<16;i+)if(i%2=0)detsii=detpii/2;elsedetsii=detqi(i-1)/2;/將detp和detq用一個(gè)數(shù)組表示for(i=0;i<T;i+)for(j=0;j&l

11、t;(2*T);j+)if(j<16)aij=J0ij;elseif(j=T+i)aij=;elseaij=;for(i=0;i<T;i+)t=aki/aii;for(j=0;j<(2*T);j+)xx=aij*t;akj=akj-xx;)for(i=0;i<T;i+)(t=aii;for(j=0;j<(2*T);j+)aij=aij/t;for(i=0;i<T;i+)for(j=0;j<T;j+)niij=aij+T;/*printf("逆矩陣為：n");for(i=0;i<T;i+)(for(j=0;j<T;j+)p

12、rintf("%10.3f",niij);printf("n");*/xxx=;for(i=0;i<T;i+)xxx=;for(j=0;j<T;j+)xxx=xxx+niij*detsij;)detuii=xxx;/檢測detui滿足要求與否max=Max(detui,16);printf("max=%fn",max);*/算第n次迭代后的ufor(i=0;i<T;i+)if(i%2=0)fi1i/2=fi0i/2+detuii;elseei1(i-1)/2=ei0(i-1)/2+detuii;*for(i=0;i&

13、lt;8;i+)/下一次迭代賦初值ei0i=ei1i;fi0i=fi1i;for(i=0;i<8;i+)printf("ei0%d=%ft",i,ei0i);printf("fi0%d=%ft",i,fi0i);for(i=0;i<8;i+)pii=detpii+pi0i;qii=detqii+qi0i;while(max>eps&&h<km);printf("Alldo%dtimesn",h);sumir=0;sumii=0;/*for(i=0;i<9;i+)/平衡節(jié)點(diǎn)功率計(jì)算I1r=m

14、ulRe(ybr8i,-ybi8i,ei0i,-fi0i);I1i=mulIm(ybr8i,-ybi8i,ei0i,-fi0i);sumir+=I1r;sumii+=I1i;pi8=mulRe(ei08,fi08,sumir,sumii);qi8=mulIm(ei08fi08,sumir,sumii);printf("S9=%f+j%fn",pi8,qi8);sumpi1=0;sumpi2=0;sumqi1=0;sumqi2=0;for(j=0;j<9;j+)if(i!=j&&ybiij!=0)sumpi1=mulRe(ei0i,fi0i,ydij);

15、sumqi1=mulIm(ei0i,fi0i,ydij);sumpi2=mulRe(ei0i-ei0j,fi0i-fi0j,-ybrij,-ybiij);sumqi2=mulIm(ei0i-ei0j,fi0i-fi0j,-ybrij,-ybiij);sumpi1+=sumpi2;sumqi1+=sumqi2;spij=mulRe(ei0i,fi0i,sumpi1,-sumqi1);sqij=mulIm(ei0i,fi0i,sumpi1,-sumqi1);printf("S%d=%f+j%fn",(i+1)*10+j+1,spij,sqij);for(i=0;i<9;i+)for(j=0;j<9;j+)if(j!=i&&ybiij!=0)delta_pij=spij+spji;