|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:3 G; e5 w1 d+ B$ V' f
! W1 B. S9 j# Y( l# n0 f
function mpc=case5_014 N2 s+ ?1 t% e+ \
% MATPOWER Case Format:Version 2
& a4 l7 o; p, ]mpc.version='2';
9 I" w; O7 v. E" E' [' _%%-----Power Flow Date-----%%5 _# `/ i- A1 s
%% system MVA base, O+ X8 O( l: O5 h8 l) L
mpc.baseMVA=100;
: L4 p, E! q! n' K%% bus date6 Y! c: V' g) k
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
: l% ?/ a5 [3 Wmpc.bus=[/ v6 r) H/ `1 x+ {4 Z' y( S2 h
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
, N, a9 s! k: B% a 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;; r2 V( s) r6 j
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
: S- P8 U! w3 X" b8 A/ ? 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;/ M7 ^( E' ?% s2 u3 S
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;9 u$ g3 J; x' q3 s" C
];7 l7 l; x4 { B9 a* X
%% generator date
4 s. L5 F3 E# r8 P( Z& G% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 9 T! ]. W3 a$ M8 f6 R* f
mpc.gen =[, S; p' W* c1 U1 p2 h7 i) {
4 500 0 99990 -9999 1.050 100 1 600 0;/ j5 w% W0 C" c& G- N3 E
5 0 0 99990 -9999 1.050 100 1 600 0;
' t$ z$ K7 i( g1 @ ];
6 x) `; `& ~" T) C%% branch date* X$ q& A5 l# @+ B t
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 4 ?2 p; T8 z' ^ l
mpc.branch=[6 A3 v- z# [7 t, C5 n
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;* G, a* |' @0 R( o2 H* x: u
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
+ q# r9 J7 m5 |9 ?! U2 f 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;$ V6 x, q; L. P0 C! S
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
* P3 m2 o h1 w3 S 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
; ]& W) |, o/ L ];
( X. f8 n" i; @return; |
-
运行结果
|