%% prog038 lectura de archivos de audiodevinfo

clearvars;
close all;
clc;

%% lectura de archivo monofonico
%usaremos el comando audioread
%[nombre_arreglo,frec_muestro] = audioread('nombre_archivo.wav');
%el archivo de audio debe estar en la misma carpeta que el programa
%en caso que no este deben colocar el camino entre los directorios
%[nombre_arreglo,frec_muestro] = audioread('C:\Users\Eddie\Desktop\Escritorio 2021\Cursos 2021 Semestre 01\Programación 2021 Sem01\Programas2021Sem 01\nombre_archivo.wav');

%lectura desde la misma carpeta
[x1,fs1] = audioread('violin.wav');

%ejemplo de lectura de otra carpeta
[x2,fs2] = audioread('C:\Users\Eddie\Downloads\violinLALA.wav');

%% generar el archivo de tiempo
Nt = length(x1);      %determinamos la longitud de arreglo que contiene  el sonido
t1 = (0:Nt-1)/fs1;    %tiempo con prelocalizacion

%ejes
tmin = t1(1);
tmax = t1(Nt);
amin = -1;
amax = 1;
%tiempo sin prelocalizacion
%for n = 1:Nt
%    t1(n) = (n-1)/fs;
%end;

%grafico
figure(1)
plot(t1,x1);
title('violin');
xlabel('tiempo t (s)');
ylabel('presion sonora p(t) (Pa)');
legend('p(t)');
axis([tmin, tmax, amin, amax]);
grid on;
box on;

%sonido
sound(x1,fs1);

%% lectura de sonido stereo
[x3,fs3] = audioread('MetalFatigue.wav');

%% generar el archivo de tiempo
Nt = length(x3);      %determinamos la longitud de arreglo que contiene  el sonido
t3 = (0:Nt-1)/fs3;    %tiempo con prelocalizacion

%grafico
figure(2)
plot(t3,x3);
title('metal fatigue');
xlabel('tiempo t (s)');
ylabel('presion sonora p(t) (Pa)');
legend('p(t)');
axis([tmin, tmax, amin, amax]);
grid on;
box on;

%sonido
sound(x3,fs3);

%separamos senhales con prelocalizacion
x3izq = x3(:,1);   %el 1 se refiere a la primera columna
x3der = x3(:,2);   %el 2 se refiere a la segunda columna

%separamos son prelocalizacion
%for n1 = 1:Nt
%   x3izq(n) = x3(n,1); %el 1 se refiere a la primera columna
%   x3der(n) = x3(n,2); %el 2 se refiere a la segunda columna
%end;

%graficamos por separado
figure(3)

subplot(2,1,1)
plot(t3,x3izq);
title('metal fatigue canal izquierdo');
xlabel('tiempo t (s)');
ylabel('presion sonora p(t) (Pa)');
legend('p_{izq}(t)');
axis([tmin, tmax, amin, amax]);
grid on;
box on;

subplot(2,1,2)
plot(t3,x3der);
title('metal fatigue canal derecho');
xlabel('tiempo t (s)');
ylabel('presion sonora p(t) (Pa)');
legend('p_{der}(t)');
axis([tmin, tmax, amin, amax]);
grid on;
box on;