|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
. h4 L* f8 Y+ q% [! c2 C+ ?& h" L5 z( y8 l5 J
function mpc=case5_01: `# \* {" O" v2 _5 w% M X. r! J
% MATPOWER Case Format:Version 2
, s1 p s9 H5 X4 J. U; @% Z H/ cmpc.version='2';5 C2 v9 S1 Q) y6 _/ A0 |
%%-----Power Flow Date-----%%" _# K8 Z" F$ Z& s5 ^" N
%% system MVA base4 }- H7 J% }$ V0 H
mpc.baseMVA=100;4 d& H* K9 ]$ G F1 I
%% bus date
- M2 s0 p- I5 }9 q- h% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
2 @9 C! E' T/ P% v* `; [" Fmpc.bus=[8 J9 n: ^7 ]6 A* Z
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
% [3 Y. ]( r1 U5 {' v 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
( b) H# m; n: ~' p) | 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;3 K2 S7 ^6 [) q) ?& O0 k
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94; t' Q5 B! C! X# ]. B
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
; r2 F* f' e" ]! R ];- o M. O a9 w' y @3 u
%% generator date" e1 L, q; ~4 I- M& l5 h
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
; z/ S# Y9 w; bmpc.gen =[
) V5 u1 C0 \, A+ m! } 4 500 0 99990 -9999 1.050 100 1 600 0;
! C1 h( p/ }3 ~! T' |7 {3 T 5 0 0 99990 -9999 1.050 100 1 600 0;
4 T3 s+ {* X5 H4 c6 W( N3 b ];
* b- D8 Q0 V8 o4 O$ \! S: L%% branch date
& M+ w1 z& L) m5 q+ X, ^1 j% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax ; s, A5 e* U- C6 e! O2 [5 m1 w
mpc.branch=[) r# r8 c+ S' `8 w8 w
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;& k7 h% W3 v: B& e
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;' {) o1 n( ~! S; c4 q! j
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;& Q" a5 P, [; C- j3 L" I: ]
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
1 ^/ a& ~# X4 L2 P! a- g6 G2 q 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;' t+ P5 _1 [6 Q
];! C8 P6 Z& I) [
return; |
-
运行结果
|