Elavläsning 2-3-dubblas med min Actaris-elmätare som inte kan läsa elförbrukning korrekt. I mitt fall ett nätaggregat.Löst

close all; clear all; clc; clf;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plots color-coded energy use from time-stamped and 
% grouped kWh meter readings
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% data row format:
% [YYYY MM DD HH MM ss, meter reading, group code}]
data = [2013 12 30 12 00 00, 28887, 2
    2013 12 30 17 40 00, 28891, 2
    2013 12 31 06 25 00, 28927, 2
    2013 12 31 13 00 00, 28933, 2
    2013 12 31 19 10 00, 28938, 2
    2013 12 31 22 30 00, 28941, 2
    2014 1 1 11 15 00, 28978, 2
    2014 1 2 00 10 00, 28991, 2
    2014 1 2 04 25 00, 29005, 1
    2014 1 2 05 15 00, 29005, 2
    2014 1 2 07 40 00, 29012, 2
    2014 1 2 12 55 00, 29030, 2
    2014 1 2 15 40 00, 29033, 2
    2014 1 2 23 00 00, 29039, 2
    2014 1 3 06 10 00, 29065, 2
    2014 1 3 18 00 00, 29075, 2
    2014 1 3 23 55 00, 29081, 2
    2014 1 4 09 00 00, 29111, 2
    2014 1 4 18 00 00, 29119, 2
    2014 1 4 23 55 00, 29125, 2
    2014 1 5 09 00 00, 29155, 2
    2014 1 5 15 55 00, 29162, 2
    2014 1 5 17 35 00, 29164, 2
    2014 1 5 23 50 00, 29169, 1
    2014 1 6 08 50 00, 29170, 1
    2014 1 6 11 00 00, 29171, 1
    2014 1 6 14 10 00, 29172, 1
    2014 1 6 17 10 00, 29173, 2
    2014 1 6 18 10 00, 29175, 2
    2014 1 6 19 10 00, 29176, 2
    2014 1 6 20 10 00, 29177, 2
    2014 1 6 21 10 00, 29178, 2
    2014 1 6 23 10 00, 29179, 2
    2014 1 6 23 55 00, 29180, 2
    2014 1 7 02 00 00, 29187, 2
    2014 1 7 03 05 00, 29190, 1
    2014 1 7 07 30 00, 29191, 1
    2014 1 7 11 00 00, 29192, 1
    2014 1 7 12 30 00, 29193, 1
    2014 1 7 17 00 00, 29195, 1
    2014 1 7 19 00 00, 29196, 1
    2014 1 7 21 00 00, 29197, 1
    2014 1 7 23 55 00, 29199, 1
    2014 1 8 10 00 00, 29200, 1
    2014 1 8 12 05 00, 29201, 1
    2014 1 8 13 50 00, 29202, 1
    2014 1 8 16 05 00, 29203, 2
    2014 1 8 18 05 00, 29205, 2
    2014 1 8 20 05 00, 29207, 1
    2014 1 8 22 05 00, 29208, 1
    2014 1 8 23 55 00, 29209, 1
    2014 1 9 06 45 00, 29210, 1
    2014 1 9 10 05 00, 29211, 1
    2014 1 9 11 45 00, 29212, 3
    2014 1 9 14 15 00, 29219, 3
    2014 1 9 16 20 00, 29224, 3
    2014 1 9 17 20 00, 29227, 1
    2014 1 9 22 00 00, 29228, 3
    2014 1 9 23 59 00, 29234, 3
    2014 1 10 01 05 00, 29237, 3
    2014 1 10 08 25 00, 29256, 1
    2014 1 10 10 35 00, 29258, 1
    2014 1 10 12 35 00, 29259, 1
    2014 1 10 16 35 00, 29260, 1
    2014 1 10 21 30 00, 29264, 1
    2014 1 10 23 30 00, 29265, 4
    2014 1 11 08 35 00, 29267, 4
    2014 1 11 11 35 00, 29268, 3
    2014 1 11 15 35 00, 29280, 4
    2014 1 11 19 35 00, 29280, 4];
%% Tick settings and time axis labels
ticksize = 1;                                 % (days)
ticklabel = 'dd/mm';                         % date format for x-ticks
timelabel = 'Datum';                             % x-axis label
hourlabel = 'HH:MM';
%% Group titles array (shown in legends)
grouptitles = {'Utrustning urkopplad' 
    'Utrustning inkopplad med gammal nätdel' 
    'Utrustning inkopplad med gammal nätdel (allt annat urdraget)'
    'Utrustning inkopplad med ny nätdel'};
%% group colors
colors= [...
    0   69  232
    255 11  17
    228 44  228
    17  239 12
    0   221 66 
    193 45  193]/255;
% groups and colors must have at least as many elements as the highest
% number group
%% Figure titles and Y axis labels
kWhPlotTitle = 'Mätarställning';
kWhAxisLabel = 'Ack. energiförbrukning/kWh';
PowerPlotTitle = 'Effekt';
PowerAxisLabel = 'Effektförbrukning/kW';
%% Split data to separate arrays
t = datenum(data(:, 1:6));
meter = data(:, 7)';
meter = meter - meter(1);                % "tare" values (optional)
group = data(:, 8);
%% Create ticks from time array
ticks = [t(1)-0.5, t(1)+0.5:ticksize:t(end)+ticksize];
%% Plot total energy use
figure(1)
% title(kWhPlotTitle);
hold on
%% Plot as single lines, with different groups mapped to different colors
handles = zeros(size(grouptitles));
hold on
for i = 1:size(t)-1,
    handle = plot([t(i), t(i+1)], [meter(i), meter(i+1)], ...
        'color', colors(group(i), :), 'linewidth', 3);
    
    % Save handles for legend
    if handles(group(i)) == 0
        handles(group(i)) = handle;
    end
end
%% Set labels and stuff
set(gca, 'XTick', ticks)                  % use "ticks" as ticks
datetick('x', ticklabel, 'keepticks')     % set readable date labels
ylabel(kWhAxisLabel)
xlabel(timelabel)
grid on
set(gca, 'XMinorGrid','on')
legend(handles(find(handles ~= 0)), grouptitles(find(handles ~= 0)), 'location', 'northwest')
%% Create second axis with hours display
ax1 = gca;
ax2 = axes('Position',get(ax1,'Position'),...
           'XAxisLocation','top',...
           'YAxisLocation','right',...
           'Color','none',...
           'XColor','k','YColor','k');
set(ax2, 'XTick', ticks(1):0.5:ticks(end));
set(ax2, 'YTick', []);
datetick('x', hourlabel, 'keepticks');
%% Plot deltaQ/deltaT (Power use)
figure(2)
% title(PowerPlotTitle);
dQ = diff(meter);        % find difference between adjacent data values
dt = diff(t)'.*24;   % --||-- adjacent time stamps ( .*24 for hours)
%% Insert dummy data since diff() lowers size
dQ = [dQ 0];  
dt = [dt dt(end)];
P = dQ./dt;
%% Plot power in filled rectangles
handles = zeros(size(grouptitles));
for i = 1:size(t)-1,
    handle = patch([t(i); t(i); t(i+1); t(i+1)], [0; P(i); P(i); 0], colors(group(i), :));
    
    % Save handles for legend
    if handles(group(i)) == 0
        handles(group(i)) = handle;
    end
end
%% Set labels and stuff
ylabel(PowerAxisLabel)
xlabel(timelabel)
set(gca, 'XTick', ticks)         % use "ticks" as ticks
datetick('x', ticklabel, 'keepticks')     % set readable date labels
grid on
set(gca, 'XMinorGrid','on')
legend(handles(find(handles ~= 0)), grouptitles(find(handles ~= 0)), 'location', 'northwest')
%% Create second axis with hour display
ax2 = gca;
ax2 = axes('Position',get(ax1,'Position'),...
           'XAxisLocation','top',...
           'YAxisLocation','right',...
           'Color','none',...
           'XColor','k','YColor','k');
set(ax2, 'XTick', ticks(1):0.5:ticks(end), 'YTick', []);
datetick('x', hourlabel, 'keepticks');