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

hdxjm

潮流计算的：[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
    VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
% }, }6 x6 G/ F0 o# x+ q% n' C[F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, ...
4 X; i$ S+ r: e5 y! \    RATE_C, TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST] = idx_brch;
8 c6 ]! p$ G5 `/ }: [  \[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
1 o, u( X" r4 U7 P, k    GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;

; S: r9 x/ j" }( V' Y4 `0 Y( Z. K3 S%% default arguments
if nargin < 4
, G  I8 n( R4 `  K0 a) L    solvedcase = '';                %% don't save solved case
    if nargin < 3
        fname = '';                 %% don't print results to a file
) C6 m2 ^# J5 u# r        if nargin < 2
            mpopt = mpoption;       %% use default options
; |# [5 G$ ]; X1 c- y            if nargin < 1
                casename = 'case9'; %% default data file is 'case9.m'
            end
0 T* b: o2 `& h. e, [: K        end
    end
end
2 A, }5 b; D; h
. z' ^/ |7 V& j  ?2 l7 G6 j状态估计的：[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
7 p3 m8 x) ?; Q( [5 W( A8 j        VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
7 }9 T. P" _# ~. ]% ][GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
        GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;
  p( @* c. I  w[CP, CV, CC, CA, CONV_BUS, CONV_TYPE, CONV_TAP, CONV_XC, ...
        V_DC, I_DC, P_DC, Q_DC, ANGLE, MIMT, CRTL_TYPE, SETVALUE] = idx_dc;
( I; j. n7 N4 A: [$ v; M% u& L3 L
3 _/ [4 B, `0 ^$ n/ a6 E%% default arguments
if nargin < 4
    acmode = 4;
    dcmode = 1;
: v2 b% }- _& w        if nargin < 3
                                    %% for state estimation, form measurement system
& ?& t- @! i2 j& |- R8 V' Y+ o        fname = 'result.dat';
1 h  k5 C$ x, T: h% o" i4 F                if nargin < 2
4 X/ x/ Z, K% ^            measuresys = 'case14_result';
            if nargin < 1
' R- o. ~5 o, Y) y; Y              
                casename = 'case14_2';
            end
                end
. X$ Y0 T" {; ^% a5 [2 @/ j6 x        end
else
    flag = 1;
- J+ [9 ?; i2 D0 b, Jend

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

任我行2001

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

yuguanghai

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

任我行2001

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

yaoyaohsm

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

dayinfei

看不懂说的都是些什么。

electrons

默默打个酱油