paulo.m

%% Constants
% Number of coefficients for LPC
number_of_coefficients = 4;

% Number of audio files per audio source
number_of_audio_files = 5;


%% Memory allocation
% Initialize audio sources' matrixes
A = zeros(number_of_audio_files, number_of_coefficients+1);
B = zeros(number_of_audio_files, number_of_coefficients+1);
D = zeros(number_of_audio_files, number_of_coefficients+1);
E = zeros(number_of_audio_files, number_of_coefficients+1);
Ez = zeros(number_of_audio_files, number_of_coefficients+1);
G = zeros(number_of_audio_files, number_of_coefficients+1);

audio_sources = {A, B, D, E, Ez, G};
audio_sources_str = {'A', 'B', 'D', 'E', 'Ez', 'G'};

number_of_audio_sources = size(audio_sources, 2);


%% Load data into memory
for i = 1:1:number_of_audio_sources
    % Load tmp vars with data
    tmp_source = audio_sources{i};
    tmp_source_str = audio_sources_str{i};

    % Fill audio sources' matrixes
    for j = 1:1:number_of_audio_files
        % Load file from disk
        data_file = strcat('./', tmp_source_str, num2str(j-1), '.wav');
        [data, Fs] = audioread(data_file);

        % Get coefficients
        [tmp_source(j, :), g] = lpc(data(:, 1), number_of_coefficients);
    end

    % Save data to audio_sources cell
    audio_sources{i} = tmp_source;
end


%% Data processing
% TODO


%% Data visualization
figure;
hold on;
grid on;

title('Scatter3');
xlabel('x');
ylabel('y');
zlabel('z');
for i = 1:1:number_of_audio_sources
    tmp_source = audio_sources{i};
    scatter3(tmp_source(1, 2), tmp_source(1, 3), tmp_source(1, 4));
end

hold off;