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

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;
[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, ...
0 H) t  w) r& h$ }8 F7 O    GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;

%% default arguments
if nargin < 4
    solvedcase = '';                %% don't save solved case
2 U; m8 E" T3 k; O    if nargin < 3
        fname = '';                 %% don't print results to a file
        if nargin < 2
            mpopt = mpoption;       %% use default options
4 p2 d8 |; K. R6 n, x+ ~            if nargin < 1
: d, K/ A7 h3 k$ N                casename = 'case9'; %% default data file is 'case9.m'
: A6 J$ m$ I# ?1 B6 p$ O8 C" a7 T            end
) \3 v: z, v/ W9 b5 ^) v; r        end
    end
$ s) t& h3 ~$ d% iend

( {! v1 t$ V  }/ O. b& K4 P状态估计的：[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
& L, {( a$ Z  F8 |# S7 `' Q8 n7 |, w        VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, ...
        GEN_STATUS, PMAX, PMIN, MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN] = idx_gen;
* j) {/ v( {9 |# s0 u[CP, CV, CC, CA, CONV_BUS, CONV_TYPE, CONV_TAP, CONV_XC, ...
: t" k8 ]8 ]9 q. p* }& {5 i        V_DC, I_DC, P_DC, Q_DC, ANGLE, MIMT, CRTL_TYPE, SETVALUE] = idx_dc;
* G  B5 _) _, N- Y/ K
# C0 C0 @6 h1 e%% default arguments
5 ?$ G, T+ H1 m+ gif nargin < 4
4 q  L$ b! ^" n1 x# \4 @( J    acmode = 4;
    dcmode = 1;
( z0 P2 b. o3 k: W        if nargin < 3
                                    %% for state estimation, form measurement system
1 m% d: v5 B+ [' v6 t        fname = 'result.dat';
8 C1 ]0 x" ]6 D" N7 K1 k8 i                if nargin < 2
            measuresys = 'case14_result';
! m1 C/ C; y6 j; g' A7 @            if nargin < 1
              
                casename = 'case14_2';
            end
                end
        end
2 J  f: z$ q  I9 e' b! M% relse
    flag = 1;
: z# K: @% v8 J8 R2 o% p' [0 d+ A! Yend

4 D+ Y7 F  f; n上面这段程序什么意思啊？

任我行2001

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

yuguanghai

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

任我行2001

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

yaoyaohsm

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

dayinfei

看不懂说的都是些什么。

electrons

默默打个酱油