%%prog 20230531 lectura de audios
clearvars;
clear all;
close all;
clc;

[y,Fs] = audioread('violin.wav');
N = length(y);
duracionSeg = N/Fs; %calculamos la duracion en segundos basado en la frecuencia de muestreo

t = linspace(0,duracionSeg,N); %genera un arreglo de N elementos que va desde 0 a duracionSeg

figure(1)
plot(t,y)
title('violin.wav')
xlabel('tiempo (s)')
ylabel('amplitud')
axis([0 t(end) -1 1]) % minimo_x maximo_x minimo_y maximo_y

%sound(y,Fs)

y_modificado = y * 40; %amplificamos la señal por 10

figure(2)
plot(t,y_modificado)
title('violin.wav modificado')
xlabel('tiempo (s)')
ylabel('amplitud')
axis([0 t(end) -10 10]) % minimo_x maximo_x minimo_y maximo_y

%ciclo for que verifique que todas las muestras se encuentran entre -1 y 1
%en caso que no esté dentro de esos valores, asignarle 1 o -1 según el signo del elementos
audiowrite('violin_amplificado.wav',y_modificado,Fs)

%%lectura de archivo modificado
[y2,Fs] = audioread('violin_amplificado.wav');
N2 = length(y);
duracionSeg2 = N2/Fs; %calculamos la duracion en segundos basado en la frecuencia de muestreo

t2 = linspace(0,duracionSeg2,N2); %genera un arreglo de N elementos que va desde 0 a duracionSeg

figure(3)
plot(t2,y2)
title('violin.wav')
xlabel('tiempo (s)')
ylabel('amplitud')
axis([0 t2(end) -1 1]) % minimo_x maximo_x minimo_y maximo_y

%sound(y2,Fs)