0001 function test_se
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018 idx.idx_zPF = [1;2];
0019 idx.idx_zPT = [3];
0020 idx.idx_zPG = [1;2;3];
0021 idx.idx_zVa = [];
0022 idx.idx_zQF = [];
0023 idx.idx_zQT = [];
0024 idx.idx_zQG = [];
0025 idx.idx_zVm = [2;3];
0026
0027
0028 measure.PF = [0.12;0.10];
0029 measure.PT = [-0.04];
0030 measure.PG = [0.58;0.30;0.14];
0031 measure.Va = [];
0032 measure.QF = [];
0033 measure.QT = [];
0034 measure.QG = [];
0035 measure.Vm = [1.04;0.98];
0036
0037
0038 sigma.sigma_PF = 0.02;
0039 sigma.sigma_PT = 0.02;
0040 sigma.sigma_PG = 0.015;
0041 sigma.sigma_Va = [];
0042 sigma.sigma_QF = [];
0043 sigma.sigma_QT = [];
0044 sigma.sigma_QG = [];
0045 sigma.sigma_Vm = 0.01;
0046
0047
0048 nbus = 3;
0049 [success, measure, idx, sigma] = checkDataIntegrity(measure, idx, sigma, nbus);
0050 if ~success
0051 error('State Estimation input data are not complete or sufficient!');
0052 end
0053
0054
0055 casename = 'case3bus_P6_6.m';
0056 type_initialguess = 2;
0057 [baseMVA, bus, gen, branch, success, et, z, z_est, error_sqrsum] = run_se(casename, measure, idx, sigma, type_initialguess);