|
楼主 |
发表于 2013-9-15 11:11:54
|
显示全部楼层
一下是代码:8 @1 ^) y! q+ ~
0 y, Q9 w+ P( \$ j) |
function mpc=case5_01
0 \: ]" r& ^3 t% C$ H% MATPOWER Case Format:Version 2
: _$ \6 D! t( E) B5 s# rmpc.version='2';, |7 M9 B( W& v. _# |
%%-----Power Flow Date-----%%, M, D2 ]1 Q5 b( f
%% system MVA base' E0 p! Q5 x7 {: [, w; Q
mpc.baseMVA=100;) ^# x& m# y! W: n0 X
%% bus date% E0 @! j4 k; U) { l. `
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin* M0 g2 g" U" b1 Z8 k3 n9 x+ o
mpc.bus=[
O2 q: r) [' |9 k5 t2 A 1 1 160 80 0 0 1 1 0 100 1 1.1 0.94;
% P: q+ x0 f4 } ` 2 1 200 100 0 0 1 1 0 100 1 1.1 0.94;3 C% h& L* a6 B# o) O# K I
3 1 370 130 0 0 1 1 0 100 1 1.1 0.94;
5 {& p9 W2 C8 n; Z 4 2 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
0 o5 l! j! l# M( }. U9 @ 5 3 0 0 0 0 1 1.050 0 100 1 1.1 0.94;
& n: _: I% t' B/ m# C ];
0 u& O' D9 U& }# {7 U7 Y2 f%% generator date" K" Q. G5 t& l, E! L! o$ N
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
* T4 O* c8 @6 q8 Zmpc.gen =[
5 j( i1 A2 }0 X5 T8 e 4 500 0 99990 -9999 1.050 100 1 600 0;) T, Q8 a0 a0 Z
5 0 0 99990 -9999 1.050 100 1 600 0;
5 n, W! c* P' a9 Q0 `. N# m ];
: n3 {2 u3 @4 ?5 \1 C% c%% branch date
" x" O- P, |# N; @' e% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax 1 L6 r( f& l6 n: _" ]8 y4 G( t
mpc.branch=[/ b1 N5 A% i+ o) E. n9 ?& l! l0 o
2 1 0.04 0.25 0.5 0 0 0 0 0 1 -360 360;
; Z4 A( u" o. R: J8 m) h 3 1 0.1 0.35 0 0 0 0 0 0 1 -360 360;
; a- \8 K1 x1 Y s1 e 3 2 0.08 0.3 0.5 0 0 0 0 0 1 -360 360;
% x9 [6 m6 \/ P4 K/ c. F 3 5 0 0.03 0 0 0 0 1.05 0 1 -360 360;
9 \% e" n9 D) X6 { 2 4 0 0.015 0 0 0 0 1.05 0 1 -360 360;
) `" M! }: W* Z% f! f7 z ];7 m) Q+ R$ h- O' U: @2 Z! _
return; |
-
运行结果
|