plot_octave_bands.m

% Data
rooms = {'MR', 'KT', 'LR', 'ST', 'UL'};
frequencies = [500, 1000, 2000, 4000, 8000];

% Ground truth (approximated from the image)
ground_truth = [
    0.55, 0.60, 0.65, 0.65, 0.60;  % MR
    0.45, 0.50, 0.55, 0.60, 0.65;  % KT
    0.15, 0.15, 0.15, 0.15, 0.15;  % LR
    1.80, 1.85, 1.90, 1.95, 1.90;  % ST
    0.45, 0.45, 0.50, 0.55, 0.50   % UL
];

% Baked data (corrected)
baked = [
    0.19, 0.21, 0.20, 0.21, 0.21;  % MR
    0.18, 0.22, 0.22, 0.22, 0.23;  % KT
    0.50, 0.50, 0.49, 0.48, 0.45;  % LR
    0.72, 0.62, 0.64, 0.60, 0.62;  % ST
    0.44, 0.43, 0.43, 0.38, 0.47   % UL
];

% Realtime data
realtime = [
    0.33, 0.42, 0.40, 0.42, 0.40;  % MR
    0.33, 0.44, 0.40, 0.44, 0.39;  % KT
    0.29, 0.36, 0.34, 0.35, 0.33; % LR
    0.64, 0.73, 0.70, 0.71, 0.70;  % ST
    0.32, 0.41, 0.39, 0.42, 0.41   % UL
];

% Create figure with subplots
figure('Position', [100, 100, 1200, 400]);

for i = 1:length(rooms)
    subplot(1, 5, i);
    hold on;
    
    % Plot ground truth
    plot(frequencies, ground_truth(i,:), 'k-', 'LineWidth', 2);
    
    % Plot JND lines (20%)
    plot(frequencies, ground_truth(i,:)*1.2, 'k--', 'LineWidth', 1);
    plot(frequencies, ground_truth(i,:)*0.8, 'k--', 'LineWidth', 1);
    
    % Plot data
    plot(frequencies, baked(i,:), 'gs-', 'LineWidth', 1.5, 'MarkerSize', 6);  % Changed from 'y' to 'g' for green
    plot(frequencies, realtime(i,:), 'rd-', 'LineWidth', 1.5, 'MarkerSize', 6);
    
    % Customize plot
    title(rooms{i});
    xlabel('Frequency (kHz)');
    if i == 1
        ylabel('RT60 (s)');
    end
    xlim([400 9000]);
    set(gca, 'XScale', 'log');
    set(gca, 'XTick', frequencies);
    set(gca, 'XTickLabel', {'0.5', '1', '2', '4', '8'});
    grid on;
    
    % Adjust y-axis limits based on data
    if i == 1  % MR
        ylim([0 1.4]);
    elseif i == 2  % KT
        ylim([0 2.1]);
    elseif i == 3  % ST
        ylim([0 2.5]);
    else  % UL
        ylim([0 1.5]);
    end
    
    % Add legend only to the first subplot
    if i == 1
        legend('Ground Truth', 'JND', '', 'Baked', 'Realtime', 'Location', 'best');
    end
    
    hold off;
end

% Add overall title
sgtitle('RT60s over different frequency bands for the rooms');