|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:1 r- S: K- B2 b) @, j
) @6 F9 @& g4 cfunction mpc=case5_01( [# F" H( i O, g- M9 { w+ m1 N
% MATPOWER Case Format:Version 2
6 N" o I2 B/ |9 S3 a2 Qmpc.version='2';
1 a: f) g& |0 A! w' u%%-----Power Flow Date-----%%
' y9 y3 |7 I8 {( ^% B- H; U%% system MVA base9 U& `9 C9 H9 I. N) e
mpc.baseMVA=100;
; r( Q3 J! O. [ }%% bus date
9 a) h6 ^+ a- g- r% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
7 ^% p% I3 \ G$ smpc.bus=[
* T1 V, U% Z- Q$ S u 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;, V" ^8 M. O. ?' L. M
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;& l0 C) V7 y. `0 H& m8 f* d
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
: v0 s; r' T7 B$ ` 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
. G7 u" S, @3 V. ~ 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
2 `3 I4 ]0 N1 K$ } ];$ G; X+ a* a3 |3 u5 T7 G
%% generator date
2 W1 @8 e1 T+ J1 L- {% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
8 {) `- Y, n! _+ T( Nmpc.gen =[1 y+ v. j) n/ Q
4 500 0 99990 -9999 1.050 100 1 600 0;
) Y" D/ a' z. y$ S( d. r" k 5 0 0 99990 -9999 1.050 100 1 600 0;& o2 V* Y7 Y: m$ v7 G% }1 z
];
) ^4 x }4 l6 N- A%% branch date; a5 Y( ?6 J# e8 r* n
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 4 x: {/ a; X& M7 S5 c$ }( E
mpc.branch=[
$ l% F! ~& ]* B+ H5 } 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
{, { N; W2 N5 d 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;( S: z! R& V4 D" e: d( d$ O
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;0 { M2 g) V, S& ]
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;8 m. L6 P( z% ^6 Z* W
2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;# ^# m r6 m$ y' G$ J
];
} A; y7 p6 E: d: J! L6 creturn; |
-
运行结果
|