matlab讨论

ssunny1234

   大家有什么matlab方面的问题，心得欢迎在此讨论

ssunny1234

上传一个MATLAB教程PPT课件

继续心灵之旅

讨论都没有人啊，上传文件啊

electrons

如题：在运用matpower求解潮流时，显示索引超出矩阵维数    该如何解决问题？？？
& D+ w! \9 G1 P! k以下是数据文件：
function mpc=case4_2
  K& Y$ u. d- @8 D+ T: M%CASE power flow data for modified 4 buses and 2 generators
# m) R: Y8 L' T5 h( I) x3 Y: j%matpower case format:version 2
. W& A" c8 B2 N$ ]" f/ {mpc.version='2';
%%-----POWER FLOW DATA-----%%
%% system MVA base
, C" q/ Q( s" u8 ?mpc.baseMVA=100;

+ B# l" j; c+ J! `  _+ a5 Q%% bus data
9 t* v/ g# D/ `7 p% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
mpc.bus=[
    1     1     0     0     0      0      1     1      0     100      1      1.1     0.94
+ t4 L  c& M  p7 j  I3 @    2     1    50    30     0    0.05     1     1      0     100      1      1.1     0.94
    3     2     0     0     0      0      1     1      0     100     1.05    1.1     0.94
2 v: T2 W) S2 }1 u+ A$ K' A    4     3    15    10     0      0      1     1      0     100     1.05    1.1     0.94
4 p& T/ `; C% m8 s    ];

%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
mpc.gen=[
    4     0     0      99990      -9999      1.05      100     1      30    0
( N  y1 J' m3 Z+ g6 O4 ]    5    20     0      99990      -9999      1.05      100     1      30    0
    ];
7 P& }1 O8 U% o/ N1 v' I
$ }( @, y; m" t, a' t) b9 d%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status angmin amgmax
mpc.branch=[
/ B/ ?2 S4 f! @1 C% B1 L, Y' W    1   2        0       0.188     -0.0774      0   0   0   10   0   1   -360   360
# W: K0 p, G5 q/ k# Z. A6 V    1   3    0.130165   0.247934   0.025864     0   0   0    0   0   1   -360   360
: j* m+ q3 ^2 N* H    1   4    0.173554   0.330579   0.189673     0   0   0    0   0   1   -360   360
+ h0 [  k/ d+ U    4   3    0.260331   0.495868   0.051728     0   0   0    0   0   1   -360   360
    ];

* A# r, d9 s+ \# ~return;