Lab 1a raman

% clear
% close all
clc

%%

labels{1} = '30 sec, 20 ml/min';
labels{2} = '30 sec, 30 ml/min';
labels{3} = '30 sec, 40 ml/min';
labels{4} = '120 sec, 20 ml/min';

data{1} = load('0.5min_20ml.txt');
data{2} = load('0.5min_30ml.txt');
data{3} = load('0.5min_40ml.txt');
data{4} = load('2min_20ml.txt');

faint = 0.25;
clr_f{1} = [0.6350 0.0780 0.1840 faint];
clr_f{2} = [0.4660 0.6740 0.1880 faint];
clr_f{3} = [0 0.4470 0.7410 faint];
clr_f{4} = [0.4940 0.1840 0.5560 faint];

% data_smooth = cell(1,4);
% for i = 1:length(data)
%     x = data{i}(:,1);
%     y = data{i}(:,2);
%     yy = smooth(y,'rlowess')
%     data_smooth{i} = yy;
% end

clr{1} = [0.6350 0.0780 0.1840];
clr{2} = [0.4660 0.6740 0.1880];
clr{3} = [0 0.4470 0.7410];
clr{4} = [0.4940 0.1840 0.5560];


figure

val_range{1} = [100 240];
val_range{2} = [1250 1400];
val_range{3} = [1540 1640];
val_range{4} = [2500 2750];
span = [0.1 0.5 0.1 0.4];
ylims{1} = [40 140];
ylims{2} = [40 90];
ylims{3} = [0 400];
ylims{4} = [40 100];

for i = 1:length(val_range)
    A = find(data{i}(:,1) >= val_range{i}(1), 1);
    B = find(data{i}(:,1) >= val_range{i}(2), 1);
    ind_range{i} = [A:B];
end

%
% fit_types = ['gauss2' 'gauss1' 'gauss3' 'gauss1'];

for i = 1:length(data)

    subplot(2,length(val_range),[1:length(val_range)])
    hold on
    grid on
    plot(data{i}(:,1),data{i}(:,2),'-','color',clr{i})
    xlabel('wavenumber, cm^{-1}')
    ylabel('intensity, a.u.')
    legend(labels)
    ftype()

    for j = 1:length(val_range)
        subplot(2,length(val_range),length(val_range)+j)
        [x1,y1,yy1] = smoother(data{i},ind_range{j},span(j));
        hold on
        grid on
        plot(x1,y1,'color',clr_f{i})
        plot(x1,yy1,'color',clr{i},'Linewidth',2)
        xlim(val_range{j})
        xlabel('wavenumber, cm^{-1}')
        ylabel('intensity, a.u.')
        ylim(ylims{j})
        ftype()
    end

%     for j = 1:length(val_range)
%         subplot(2,length(val_range),length(val_range)+j)
%         hold on
%         grid on
%         [x1,y1,yy1] = fitter(data{i},ind_range{j},fit_types(j));
%         xlim(val_range{j})
%         xlabel('wavenumber, cm^{-1}')
%         ylabel('intensity, a.u.')
%         ftype()
%     end

%     subplot(2,4,5)
%     [x1,y1,yy1] = smoother(data{i},ind_range{1});
%     hold on
%     grid on
%     plot(x1,y1,'color',clr_f{i})
%     plot(x1,yy1,'color',clr{i})
%     xlim(val_range{1})
%     xlabel('wavenumber, cm^{-1}')
%     ylabel('intensity, a.u.')
%
%     subplot(2,4,6)
%     [x,y,yy] = smoother(data{i},ind_range{2});
%     hold on
%     grid on
%     plot(x,y,'color',clr_f{i})
%     plot(x,yy,'color',clr{i})
%     xlim(val_range{2})
%     xlabel('wavenumber, cm^{-1}')
%     ylabel('intensity, a.u.')
%
%     subplot(2,4,7)
%     [x,y,yy] = smoother(data{i},ind_range{3});
%     hold on
%     grid on
%     plot(x,y,'color',clr_f{i})
%     plot(x,yy,'color',clr{i})
%     xlim(val_range{3})
%     xlabel('wavenumber, cm^{-1}')
%     ylabel('intensity, a.u.')
%
%     subplot(2,4,8)
%     [x,y,yy] = smoother(data{i},ind_range{4});
%     hold on
%     grid on
%     plot(x,y,'color',clr_f{i})
%     plot(x,yy,'color',clr{i})
%     xlim(val_range{4})
%     xlabel('wavenumber, cm^{-1}')
%     ylabel('intensity, a.u.')
end


%%

function [x,y,yy] = smoother(data,range,span);
x = data(range,1);
y = data(range,2);
yy = smooth(y,span,'loess');
end

function [x,y,yy] = fitter(data,range,degree)

faint = 0.25;
clr_f{1} = [0.6350 0.0780 0.1840 faint];
clr_f{2} = [0.4660 0.6740 0.1880 faint];
clr_f{3} = [0 0.4470 0.7410 faint];
clr_f{4} = [0.4940 0.1840 0.5560 faint];

clr{1} = [0.6350 0.0780 0.1840];
clr{2} = [0.4660 0.6740 0.1880];
clr{3} = [0 0.4470 0.7410];
clr{4} = [0.4940 0.1840 0.5560];

x = data(range,1);
y = data(range,2);
F = fit(x,y,degree);

plot(x,y,'color',clr_f{i})
plot(F,x,y,'color',clr{i})
end

function ftype()
set(gca,'FontSize',13,'fontWeight','bold')
set(findall(gcf,'type','text'),'FontSize',13,'fontWeight','bold')
end