潮流计算和状态估计程序请教

hdxjm

潮流计算的：[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
/ w# t* L3 K; r- S    VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
. ?% l3 F/ x6 N# f[F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, ...
    RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST] = idx_brch;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
* F/ b/ J, h1 q4 r5 B8 [, `    GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;
) E7 z2 C; Q# B' N+ r
%% default arguments
2 @/ f  R  ]0 L5 zif nargin < 4
    solvedcase = '';                %% don't save solved case
8 [2 N% k, K" a. i+ R0 y    if nargin < 3
7 s$ t5 T1 l2 c, v        fname = '';                 %% don't print results to a file
        if nargin < 2
            mpopt = mpoption;       %% use default options
$ J5 ?: H( ]1 G8 j" x            if nargin < 1
                casename = 'case9'; %% default data file is 'case9.m'
# N0 \) Y- L1 P/ u3 r' {  B            end
        end
- s3 ]3 g% p5 h. P; N    end
end
# \  ?$ x! P. W
状态估计的：[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
3 H0 h& q' j, E  r5 f) d        VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
' v, M! {  t6 ^5 A/ Q0 e        GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;
2 M" H3 M: K6 R1 K# R- O[CP, CV, CC, CA, CONV_BUS, CONV_TYPE, CONV_TAP, CONV_XC, ...
! V( |2 x! [+ l; n% ~7 g        V_DC, I_DC, P_DC, Q_DC, ANGLE, MIMT, CRTL_TYPE, SETVALUE] = idx_dc;
2 d: q" ^  Z" F# Q4 n& p
. W2 g& [3 s8 Y- ^, J# _! _8 A%% default arguments
if nargin < 4
    acmode = 4;
& e1 S; @/ M: t/ X    dcmode = 1;
        if nargin < 3
0 u+ {. ?2 J% O- A4 ^* L- S: ~7 `                                    %% for state estimation, form measurement system
        fname = 'result.dat';
, u- e" \8 t) Z! y8 a3 T1 l                if nargin < 2
. j6 G* K7 O# |$ ]# m            measuresys = 'case14_result';
* F7 \4 [% C" L! s            if nargin < 1
              
$ F  A1 f. k5 X) T: B/ v- _( B                casename = 'case14_2';
            end
, }+ U7 L4 a# h3 A0 x                end
        end
else
    flag = 1;
& B9 w5 Y7 X, ^% Q4 A4 i8 m) J' vend

上面这段程序什么意思啊？

任我行2001

根据输入量个数不同，进而进行不同的操作！

yuguanghai

怎么还是matlab啊,我晕!!!!!!!!!!!!!!!!!

任我行2001

这段话是MATPOWER中的语句。前面几句是将三个输入的矩阵对每一列给定了一个INDEX，方便到时候引用。
后面nargin是根据输入变量的多少，进行不同的具体操作。

yaoyaohsm

看了好几个星期的这样程序 还是看不懂啊

dayinfei

看不懂说的都是些什么。

electrons

默默打个酱油