File audio

clear
close all

[x, Fc] = audioread('test2.wav');

sound(x, Fc);   % Suona nella periferica audio

X = fft(x);         % Calcolo della fft
[Ncampioni Ncanali] = size(X);
X = 1/Ncampioni*fftshift(X);    
P = abs(X).^2;      % Periodogramma
P = P./max(max(P)); % Normalizzazione
P = 10*log10(P);    % Periodigramma in dB

dF = Fc/Ncampioni; 
f = -Fc/2:dF:Fc/2-dF;  % Asse delle frequenze

figure;
subplot(2,1,1);
plot(f,P(:,1));
ylabel('Periodogramma canale SX');
axis([0 15000 -100 0]);
grid on;
subplot(2,1,2);
plot(f,P(:,2));
xlabel('Frequenza [Hz]');
ylabel('Periodogramma canale DX');
axis([0 15000 -100 0]);
grid on;

pause();

% Sovrappone il segnale all'originale con un ritardo di 2 secondi
ritardo_secondi = 2;
rit_campioni = ritardo_secondi*Fc;
y = zeros(Ncampioni+rit_campioni,2);
y(1:Ncampioni,:) = x;
y(rit_campioni+1:Ncampioni+rit_campioni,:) = ... 
          y(rit_campioni+1:Ncampioni+rit_campioni,:) + x;

sound(y, Fc);   % Suona nella periferica audio