|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
# }/ B4 E( B' j8 n+ B
( u4 p- V* a: z: x9 `8 _/ ]function mpc=case5_01
2 s5 m- U- ~: h6 A% MATPOWER Case Format:Version 2" G" R5 A* `+ I }
mpc.version='2';9 S% H1 L% e& P, h" J5 ]
%%-----Power Flow Date-----%%* B* R8 @9 `7 Y, x }
%% system MVA base c- w' [0 t3 c+ e
mpc.baseMVA=100;
; N* Q3 B0 T# a( g+ ]%% bus date
& k3 P3 D5 h& E. i% y/ r. o* X, y% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
, c/ q6 m+ d1 Z3 w# t2 } rmpc.bus=[% v U6 S' O. f# `
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
4 S6 ?& S! o; b3 q: u. b2 z 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;% h2 l) Y4 f# @3 }
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;6 c' ~, d$ g/ v5 q, X0 r1 J! [
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
! Z3 T6 Q( c% B* w7 b# N$ o 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;9 J, A4 F9 l0 I2 b" b3 D
];. a% R) a$ t. l; |! i9 K$ K
%% generator date/ S9 S; }3 J" G1 E! _/ I' R2 _2 Y
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
( q) D% }3 E* k; R- b- Zmpc.gen =[
J: e% e2 C+ [' m6 J# [ 4 500 0 99990 -9999 1.050 100 1 600 0;0 B2 i0 @9 w/ F4 V2 o/ g
5 0 0 99990 -9999 1.050 100 1 600 0;
5 |% Y4 z) @6 v! e ];5 f- t+ S* ~! h+ q9 n9 I
%% branch date
" O. e1 y& d* P% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax
5 J2 Q0 b1 q! h" G# u7 Y% pmpc.branch=[
; G/ O6 { l8 [" ]# L$ J* E 2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
' U' s( V. _( {& V" t 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;2 F9 k5 M8 L: z& Y
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
7 X9 i! I8 h; C) H. x- o7 u 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
2 Q8 f0 ?/ Y7 J, H8 _. Q0 u 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;+ Y5 I( s8 j* `. a
];% y' o! s% O- b0 i$ N* S* ]
return; |
-
运行结果
|