Henrik Klaus Martin

close all, clear;
%% Histogram of measurements:

%load m100_10.txt; measurement = m100_10;
%load m1000_2.txt; measurement = m1000_2;
load m1000_2+100.txt; measurement = m1000_2_100;

%Testdata:
%measurement = [1 1; 2 2; 2 2; 3 3; 4 4; 2 2; 1 1; 4 1; 4 1]; % 3 * 2,4 + 2 * 1 + 1 * 3.

containers = 50;

%Histogram
delta = (max(measurement(:,1)) - min(measurement(:,1)))/containers;

histVal = zeros(1,containers);
for i = 1:size(measurement,1),
    if measurement(i,1) == min(measurement(:,1)),
        histVal(1) = histVal(1) + 1;
    else
        histVal(ceil((measurement(i,1)- min(measurement(:,1)))/delta)) ...
            = histVal(ceil((measurement(i,1)- min(measurement(:,1)))/delta)) + 1;
    end
end

% Weighted histogram
value = zeros(1,containers);
for i = 1:size(measurement,1),
    if measurement(i,1) == min(measurement(:,1)),
        value(1) = value(1) + measurement(i,2);
    else
        value(ceil((measurement(i,1)- min(measurement(:,1)))/delta))...
            = value(ceil((measurement(i,1)- min(measurement(:,1)))/delta))...
            + measurement(i,2);
    end
end

figure()
bar(min(measurement(:,1))+delta/2+(0:containers-1)*delta,histVal,1);
title('Histogram');
xlabel('x-value'), ylabel('Number');

figure()
bar(min(measurement(:,1))+delta/2+(0:containers-1)*delta,value,1);
title('Weighted histogram');
xlabel('x-value'), ylabel('sum of importance');

figure()
plot(measurement(:,1),measurement(:,2),'bx');
title('');
xlabel('x-value'), ylabel('importance');