|
|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:
& i, N. |4 q9 ~( j: y
% n: }% E* U4 k2 T+ a" Dfunction mpc=case5_019 Q: z3 ?! W) u! i
% MATPOWER Case Format:Version 2$ X; O V* O" f- Z% x( y
mpc.version='2';
8 ] |5 ?2 t6 h( H% G%%-----Power Flow Date-----%%: z+ ~" k9 |" H0 A- m: M
%% system MVA base$ B& m8 r+ r: q1 m
mpc.baseMVA=100;4 k8 A+ }( d: j, Z1 T7 Z. ~9 f
%% bus date( E8 v# u3 h( ]
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
+ g- t9 K4 a$ F- q( R7 u/ Y0 y# dmpc.bus=[
. V0 B, [& {6 t( D' ^% k9 T( r 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;' f, d2 r% @: M/ O; e: `
2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;# s+ D3 `5 ~+ b: T, j# h; @+ s9 _
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;( ^1 [. Y# V% s
4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
2 M6 p! A7 S" Y" T+ u5 Z& i2 G 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;# A' Z4 o, Q- h5 ]4 G+ W* F
];7 @8 [ [9 d! R% |) g
%% generator date( W* S; r- w& }3 T7 |7 L* Z
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
% X5 u" g2 N8 U# x, m0 I% ^mpc.gen =[8 G$ ?2 o- z2 n3 F* j% ]2 V; B
4 500 0 99990 -9999 1.050 100 1 600 0; I) t7 L7 K( o
5 0 0 99990 -9999 1.050 100 1 600 0;- p8 ?; ~. c9 T8 v! p' j
];. q# a i1 R! D* t9 @4 j
%% branch date- m- I& F& _* o
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax - ^2 ]# J: t2 J$ p" `8 p0 C) |$ U* H
mpc.branch=[9 R' C' E: c" e. G, |
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;9 _1 {, H0 ~3 P
3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;9 @2 Q) `% y: }2 ?, {8 ^: I- W
3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
; F: u5 S- ]7 s6 f" w" B. s 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
: I1 Q) L- b5 m+ k2 N 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
9 t; S6 a$ r- A- V' e$ j9 _ ];
4 N, u' L% A- lreturn; |
-
运行结果
|