Graph3RealLife.m

close all
clear all

kValues = [4 9 16 36];
alpha = 0.2:0.05:0.8;
beta = 0.75;
numDev = 1;

%Switch the comments to switch between stopping point methods
%stoppingEq =  @(x) ceil(x/exp(1));
stoppingEq =  @(x) round(sqrt(x));
%stoppingEq =  @(x) 0.7 * round(sqrt(x));

%Number of synthetic or real highways you want/have
numSim = 30;

%          k                 lam                sim#
rates(1:length(kValues),1:length(alpha),1:numSim) = 10;
starts(1:length(kValues),1:length(alpha),1:numSim) = 10;
stops(1:length(kValues),1:length(alpha),1:numSim) = 10;
for simNum = 1:numSim
    for i = 1:length(alpha)
        %You can contruct using synthetic data or real data
        % if you use real data, set numSim to the number of 'Trip's you
        % have
        %highway = construct(densities(i), 1000);
        highway = reallife("Trip" + simNum + ".csv");
        
        for k = 1:length(kValues)
            %fprintf("k: %d, density: %d", kValues(k), densities(i));
            [rates(k,i, simNum), starts(k, i, simNum), stops(k,i, simNum)] ...
                = SGAS3(highway, kValues(k), stoppingEq, numDev, alpha(i));
        end
    end
end

avgPrice(1:length(kValues), 1:length(alpha)) = 50;
for i = 1:length(kValues)
    for j = 1:length(alpha)
        total = sum(rates(i, j, rates(i, j, :) >= 0));
        avgPrice(i, j) = total / length(rates(i, j, rates(i, j, :) >= 0));
    end
end

outOfGasRate(1:length(kValues), 1:length(alpha)) = 50;
for i = 1:length(kValues)
    for j = 1:length(alpha)
        total = sum(rates(i, j, rates(i, j, :) < 0));
        outOfGasRate(i, j) = -1 .* total / length(rates(i, j,:));
    end
end

avgStart(1:length(kValues), 1:length(alpha)) = 50;
for i = 1:length(kValues)
    for j = 1:length(alpha)
        total = sum(starts(i, j, stops(i, j, :) >= 0));
        avgStart(i, j) = total / length(starts(i, j, stops(i, j, :) >= 0));
    end
end
avgStop(1:length(kValues), 1:length(alpha)) = 50;
for i = 1:length(kValues)
    for j = 1:length(alpha)
        total = sum(stops(i, j, stops(i, j, :) >= 0));
        avgStop(i, j) = total / length(stops(i, j, stops(i, j, :) >= 0));
    end
end

axisColor = 'black';

for i = 1:length(kValues)
    subplot(length(kValues), 3, i*3-2);
    stem(alpha, avgPrice(i, :), 'Color', [1 0 0]);
    ylim([1.9 max([3 max(avgPrice(i, :)+0.2)])]);
    %xticks(alpha);
    ylabel("Avg. Gas Price");
    xlabel("alpha");
    title(sprintf("\\color{" + axisColor + "}k = %d", kValues(i)));
    set(gca,'color','none')
    set(gca,'XColor',axisColor,'YColor',axisColor)
    
    subplot(length(kValues), 3, i*3-1);
    stem(alpha, outOfGasRate(i, :), 'Color', [1 0 0]);
    ylim([0 max(outOfGasRate(i, :))+0.1]);
    %xticks(alpha);
    ylabel("% ran out of gas");
    xlabel("alpha");
    title(sprintf("\\color{" + axisColor + "}k = %d", kValues(i)));
    set(gca,'color','none')
    set(gca,'XColor',axisColor,'YColor',axisColor)
    
    subplot(length(kValues), 3, i*3);
    stem(alpha, avgStart(i, :), 'Color', [0 0 1]);
    hold on
    stem(alpha, avgStop(i, :), 'Color', [1 0 0]);
    ylim([0 max(avgStop(i, :))+0.1]);
    %xticks(alpha);
    ylabel("% highway before stop");
    xlabel("alpha");
    title(sprintf("\\color{" + axisColor + "}k = %d", kValues(i)));
    set(gca,'color','none')
    set(gca,'XColor',axisColor,'YColor',axisColor)
end
set(gcf, 'Position',  [100, 100, 2000, 700])
addpath('altmany-export_fig-b1a7288');
%export_fig fig.png -transparent