Skip to content
Snippets Groups Projects
Select Git revision
  • 6ead5d7bc6a5a8dc3bafd57c45c320a4e4ce5651
  • master default protected
  • developement_1 protected
  • Version_1.2.4
  • Version_1.2.3
  • Version_1.2.2
  • Version_1.2.1
  • Version_1.2.0
  • Version_1.0.1
  • Version_1.0.0
  • Version_0.1.0
  • Version_0.0.6
  • Version_0.0.5
  • Version_0.0.4
  • Version_0.0.3
  • Version_0.0.2
  • Version_0.0.1
17 results

GeometryEngine.cpp

Blame
  • Code owners
    Assign users and groups as approvers for specific file changes. Learn more.
    generate_DSSEConfData.m 6.24 KiB
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    % Code to generate the parameters for the testing of the state estimator,
    % the structure with the information of the measurements available and the
    % structure with the accuracies of those measurements
    %
    % @author Andrea Angioni <aangioni@eonerc.rwth-aachen.de>
    % @copyright 2018, Institute for Automation of Complex Power Systems, EONERC
    % @license GNU General Public License (version 3)
    %
    % dsse
    %
    % This program is free software: you can redistribute it and/or modify
    % it under the terms of the GNU General Public License as published by
    % the Free Software Foundation, either version 3 of the License, or
    % any later version.
    %
    % This program is distributed in the hope that it will be useful,
    % but WITHOUT ANY WARRANTY; without even the implied warranty of
    % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    % GNU General Public License for more details.
    %
    % You should have received a copy of the GNU General Public License
    % along with this program.  If not, see <http://www.gnu.org/licenses/>.
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    function [Test_SetUp,Combination_devices,Accuracy] = generate_DSSEConfData(GridData)
    % structure with parameters for the testing of the state estimator
    Test_SetUp = struct;
    Test_SetUp.N_MC = 10000; %number of MC simulation, for reliable results set between 1000 and 10 000
    Test_SetUp.limit1 = 0.000001; %treshold accuracy to interrupt newton rapson
    Test_SetUp.limit2 = 50; %maximum number of iterations
    Test_SetUp.time_steps = 1; %number of time step in each MC simulation
    % here the nodes/lines where the devices are installed should be indicated
    Combination_P = []; 
    Combination_Q = []; 
    Combination_Pseudo =  [2:GridData.Nodes_num];
    Combination_Vmagn =   [1];% 1:GridData.Nodes_num
    Combination_Vph =     [1]; 
    Combination_Imagn =   [1,5,15,35,105]; %1:GridData.Lines_num
    Combination_Iph =     [1,5,15,35,105]; 
    Combination_Pflow =   []; 
    Combination_Qflow =   [];
    % here the uncertainty of devices and pseudo measurements is indicated
    unc_dev = 0.01/3;
    unc_pseudo = 0.5/3;
    %the structure with the information of the measurements available
    Combination_devices=struct;
    Combination_devices.Combination_P = Combination_P;
    Combination_devices.Combination_Q = Combination_Q;
    Combination_devices.Combination_Vmagn = Combination_Vmagn;
    Combination_devices.Combination_Vph = Combination_Vph;
    Combination_devices.Combination_Imagn = Combination_Imagn;
    Combination_devices.Combination_Iph = Combination_Iph;
    Combination_devices.Combination_Pflow = Combination_Pflow;
    Combination_devices.Combination_Qflow = Combination_Qflow;
    Combination_devices.Combination_Pseudo = Combination_Pseudo;
    
    %the accuracy for each class of device is assigned
    Accuracy_P = sqrt(2*unc_dev^2); %accuracy of active power injection measurement
    Accuracy_Q = sqrt(2*unc_dev^2); %accuracy of active power injection measurement
    Accuracy_Vmagn = unc_dev; %accuracy of active power injection measurement
    Accuracy_Vph = unc_dev; %accuracy of active power injection measurement)
    Accuracy_Imagn = unc_dev; %branches with voltage phase angle measurement
    Accuracy_Iph = unc_dev; %branches with voltage phase angle measurement
    Accuracy_Pflow = sqrt(2*unc_dev^2); %branches with voltage phase angle measurement
    Accuracy_Qflow = sqrt(2*unc_dev^2); %branches with voltage phase angle measurement
    Accuracy_pseudo = unc_pseudo;
    %structure with the accuracies of the measurements
    Accuracy=struct;
    Accuracy.Accuracy_P=Accuracy_P;
    Accuracy.Accuracy_Q=Accuracy_Q;