|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:% W3 D) m8 k& u% b$ j7 x7 V; U
' R" R; Z8 |+ ^$ \8 n" vfunction mpc=case5_01
7 m) H9 B" l; h' [+ e( a: A% MATPOWER Case Format:Version 21 H8 j$ c* k' P: L2 t) K' V
mpc.version='2';
- f2 Y" y# I4 ^. |$ @7 F: s4 G x" s+ K%%-----Power Flow Date-----%%' Q+ O& v$ J5 G. j
%% system MVA base; {& Y# g. r3 F
mpc.baseMVA=100;
4 P! s7 m0 {- I" R$ A4 l. K%% bus date
8 Z7 d2 N- q; S. ]! n, m% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin" t$ v& e# C, x# N P- E
mpc.bus=[
4 _! J9 k1 V |9 I" A" W+ ] 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;( w0 x% s6 ~% n O' e3 a
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
5 M# b: L5 ~: b$ R. t ~" f 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
7 [6 h. G$ m& L* g/ V5 h& C5 a 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
1 A( E- L: V- V# D 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;2 k! V- n Q, N+ o& a% f6 M
];
1 S/ O- y: G. H, j; k%% generator date
9 F$ ~6 ?+ ]/ |+ ?% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 1 ~" }: ?' _6 K/ X+ T. {
mpc.gen =[7 A* U( r6 ]* A8 A+ ?" w/ z
4 500 0 99990 -9999 1.050 100 1 600 0;: Q K& h$ V4 s# M! ?$ x2 S+ D( ~
5 0 0 99990 -9999 1.050 100 1 600 0;
$ H9 ^. W+ Z/ y9 t ];2 }6 K8 w, P, I ^2 @6 c0 c
%% branch date
- }! \. ?* Q: a, E2 F0 a' U% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
2 u( X9 S N* N, u9 A9 nmpc.branch=[4 U1 @3 r* h/ c; h. E T: h
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;; _8 b( W) G* G- m" o- P8 z
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;% Z$ F2 K1 t7 V+ g
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;5 ~2 j; V/ ?4 I3 E6 b2 @/ M: a7 B
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;0 k$ U" l) d& e4 y# n/ C/ H0 N
2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
5 r h0 P8 u% G% E ];# o' Z8 V9 o7 V. z3 q8 B
return; |
-
运行结果
|