|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
3 N: O8 o. A1 f7 _$ G, r' `! I9 R/ i: G+ m6 T4 N
function mpc=case5_01
$ e2 |4 Z1 M" z5 N% MATPOWER Case Format:Version 2
: P& F9 f+ o& s# Umpc.version='2';
+ C/ D$ Y; [1 O%%-----Power Flow Date-----%%
/ i: ^6 b2 d# ~8 Y%% system MVA base% M' _# [2 J4 u7 |( Q! C6 }7 j* o
mpc.baseMVA=100;7 \+ m+ M# d" j0 D2 c2 Z5 R8 @
%% bus date( q) ?; i9 p# y
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin' K% }( T. G) b1 T; Q
mpc.bus=[5 E& G* {. z3 x3 R! [" S
1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
: f+ b' U4 Y" I' P6 E 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;
9 I+ j# d9 H/ }. q+ L, ^3 u8 T; X3 t 3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;. P( C& k# }. z
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;2 y2 B8 ^* s+ }- d9 B/ r( o
5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;' Q& }+ o* ]& g( q+ x) s
];+ z( s( \* k+ W! u7 g# q
%% generator date! h2 }. D# G1 p1 X
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin 5 p; T6 i& I1 c* X$ _& ^
mpc.gen =[0 o% N/ @) E2 x6 c6 Q3 w5 Y
4 500 0 99990 -9999 1.050 100 1 600 0;
5 i) f' `! S, x4 W+ ] 5 0 0 99990 -9999 1.050 100 1 600 0;
( t2 ~7 m X* u- y4 W2 N$ s6 A ];
9 i9 T, H* A# p8 Q+ W( d+ U%% branch date; V- W3 p! A. a f- @4 s
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 5 m' T4 }* V" b1 o, g5 i" Y' H
mpc.branch=[( A S0 E# y: l3 ~* |
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
: ~# @3 q9 m1 l S 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
- t) V) b% F6 }1 B- Y/ X& K: K4 v 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;8 C3 w' l3 b! `* }) a6 C1 j
3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;% x, @6 j+ ^/ {4 |6 W) k8 R
2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;; l) b% B }$ W7 W
];
& @* q7 @; \3 W! [* [7 Vreturn; |
-
运行结果
|