%Ce fichier contient le squelette du code à écrire dans le TP. Vous
%n'avez qu'à décommenter (CTRL + T) les lignes à compléter au fur et à
%mesure de votre avancée dans l'exercice.
%
%Auteur: Patrick Marmaroli / Gilles Courtois
%Date de création: 16/04/2013
%Date de modification: 14/04/2015

clear all, close all, clc

%% chargement du signal audio
% [y,fs] = wavread('');

%% Filtrage du signal d'acquisition
% ---> filtre passe haut, fcoupure: 200 Hz

%% centrer, normaliser
% ---> ici enlever l'offset, et normaliser

%% Découpage des signaux
%--> découpez ici les signaux par fréquences

%% Algorithme de delay and sum beamforming
d = 0.03;   % distance inter-capteur
c = 343;    % vitesse du son
M = 8;      % nb de microphones

vang = linspace(0,180,90); %espace de recherche

y_d = zeros(size(y)); %variable intermediaire pouvant contenir les versions décalées des signaux d'entrées
RMS = zeros(length(vang),1); %root mean square vector

for numang = 1:length(vang)
    ang = vang(numang);
    
    % ---> ici manipulez les signaux (possibilité d'utiliser LTS2_Delay
    % )
    
    
    %calcul de l'énergie
    %     RMS(numang) = ...;
end

%% Affichage du résultat
figure;
plot(vang,RMS./max(abs(RMS)));
[val,ind] = max(RMS);
hold on
plot([vang(ind),vang(ind)],[0 1],'r');

