%% FADE IN

clearvars; close all; clc

[audio Fs nBits] = wavread('Violin.wav');
% Leer profundidad de bits (nBits) solamente es para funcion
% wavread, no funciona con audioread

N = size(audio,1);

durFadeIn = 5;  % Duración fade in (debe ser menora a duración audio)

NdurFadeIn = durFadeIn*Fs;  %Cantidad de muestras de fade in

% Queremos generar un vector de Fade In que parta en 0 y termine en 1
% en un intervalo de 'durFadIn' con Fs muestras por segundo

% el incremento está dado por:
inc = 1/(durFadeIn*Fs-1);

% Generamos vector de Fade In
fadeIn = 0:inc:1;

% Alternativamente (segun Yoram)
% fadeIn = 0:(1/(durFadeIn*Fs-1)):1;

%% Aplicamos Fade In

audioConFade = audio;   % hacemos una copia de audio original

for i = 1:NdurFadeIn 
    audioConFade(i) = audioConFade(i)*fadeIn(i);
end

% Alternativamente y de manera más eficiente:
% audioConFade(1:NdurFadeIn) = audioConFade(1:NdurFadeIn).*transpose(fadeIn);

%% Graficar

t = (0:(N-1))/Fs;

figure(1)
subplot(211)
plot(t,audio)
title('Audio Original')
xlabel('tiempo(s)')
ylabel('Presion Sonora(Pa)')
xlim([0 N/Fs])

subplot(212)
plot(t,audioConFade,'r')
title('Audio con Fade In')
xlabel('tiempo(s)')
ylabel('Presion Sonora(Pa)')
xlim([0 N/Fs])

% Escuchamos audio original
sound(audio,Fs,nBits)

% Escuchamos audio con fade in
sound(audioConFade,Fs,nBits)