matlab潮流计算 50节点,matlab潮流计算程序14节点.doc

matlab潮流计算程序14节点

JD 14;

ZL 20;

WC 0.01;

ZP [1,2,0.01335,0.04211,0.0,0.0;

1,3,0.0,0.20912,0.0,0.0;

1,4,0.0,0.55618,0.0,0.0;

1,10,0.05811,0.17632,0.0,0.0340;

1,11,0.06701,0.17103,0.0,0.0128;

2,10,0.05695,0.17388,0.0,0.0346;

2,12,0.0,0.25202,0.0,0.0;

2,14,0.05403,0.22304,0.0,0.0492;

3,4,0.0,0.11001,0.0,0.0;

3,13,0.0,0.17615,0.0,0.0;

4,5,0.03181,0.08450,0.0,0.0;

4,9,0.12711,0.27038,0.0,0.0;

5,6,0.08205,0.19207,0.0,0.0;

6,12,0.09498,0.19890,0.0,0.0;

7,8,0.22092,0.19988,0.0,0.0;

7,12,0.12291,0.25581,0.0,0.0;

8,9,0.17093,0.34802,0.0,0.0;

8,12,0.06615,0.13027,0.0,0.0;

10,11,0.04699,0.19797,0.0,0.0438;

10,14,0.01938,0.05917,0.0,0.0528];

JP [1,1,-0.478,0.039,1.0,0.0;

2,1,-0.076,-0.016,1.0,0.0;

3,1,0.0,0.0,1.0,0.0;

4,1,-0.295,-0.166,1.0,0.0;

5,1,-0.09,-0.058,1.0,0.0;

6,1,-0.035,-0.018,1.0,0.0;

7,1,-0.061,-0.016,1.0,0.0;

8,1,-0.135,-0.058,1.0,0.0;

9,1,-0.149,-0.05,1.0,0.0;

10,2,0.183,0.0,1.045,0.0;

11,2,-0.942,0.0,1.010,0.0;

12,2,-0.112,0.047,1.70,0.0;

13,2,0.0,0.174,1.90,0.0;

14,0,0.0,0.0,1.060,0.0];

[N1,N2] size ZP ;p max ZP :,1 ;q max ZP :,2 ;max1 max p,q ;

Y zeros max1 ;

for n 1:N1 i1 ZP n,1 ;j1 ZP n,2 ;r ZP n,3 ;x ZP n,4 ;g ZP n,5 ;b ZP n,6 ; m r*r+x*x; if i1 j1 Y i1,i1 Y i1,i1 +r/m-j*x/m+j*b; else Y i1,j1 Y i1,j1 -r/m+j*x/m; Y j1,i1 Y i1,j1 ; Y i1,i1 Y i1,i1 +r/m-j*x/m+j*b; Y j1,j1 Y j1,j1 +r/m-j*x/m+j*b; end

end%导纳矩阵生成

V zeros 1,JD ;XJ zeros 1,JD ;P zeros 1,JD ;Q zeros 1,JD ;V2 zeros 1,JD ;

G real Y ;B imag Y ;

for i 1:JD V i JP i,5 ; XJ i JP i,6 ; P i JP i,3 ; Q i JP i,4 ; V2 i V i *V i ;

end %赋值

DD 0;PD 1;

for i 1:JD e i V i *cos XJ i ; f i V i *sin XJ i ; end

while PD~ 0 PD 0;DD DD+1; n JD-1; H zeros n ;N zeros n ;J zeros n ;L zeros n ;R zeros n ;S zeros n ; for i 1:n for k 1:n if i~ k H i,k -B i,k *e i +G i,k *f i ; N i,k G i,k *e i +B i,k *f i ; J i,k -N i,k ; L i,k H i,k ; R i,k 0; S i,k 0; else a1 0;b1 0; for i1 1:n a1 a1+G i,i1 *e i1 -B i,i1 *f i1 ; b1 1+G i,i1 *f i1 -B i,i1 *e i1 ; end H i,k -B i,k *e i +G i,k *f i +b1; N i,k G i,k *e i +B i,k *f i +a1