clc
clear
load('final_table.mat')
load('a_ppm_pattern')
[f_table_size,~] = size(final_table);
ppm_range = max(final_table(:,1)) - min(final_table(:,1));
factor_ppm = f_table_size/ppm_range;
% peak_insten_2 = ones(size(a_ppm_pattern,1), size(final_table,2)-1);
for i = 1:size(a_ppm_pattern,1)
range_pattern = (max(final_table(:,1)) - a_ppm_pattern(i,:))*factor_ppm;
range_pattern = round(range_pattern);
peak_insten_1 = sum(final_table(range_pattern(1):range_pattern(2),2:end));
peak_insten_2(i,:) = peak_insten_1;
end
peak_insten_3 = peak_insten_2 -repmat(min(peak_insten_2),size(a_ppm_pattern,1),1);
peak_normal_t = peak_insten_3./repmat(sum(peak_insten_2),size(a_ppm_pattern,1),1);
ppm_aver = mean(a_ppm_pattern, 2);
Table_normal = [ppm_aver,peak_normal_t];
save('Table_normal','Table_normal') %Here is the table after normalization to the total intensity
bar (ppm_aver,peak_normal_t)
set(gca,'XDir','reverse')