%% Frecuencia de Muestreo - Parte 3

clearvars; close all; clc

% Modifiquemos el programa anterior para poder incorporar
% una frecuencia de oscilación como variable

Fs = 44100;     % Frecuencia de muestreo de audio (calidad CD)
T = 1/Fs;       % Intervalo
tMax = 10;       % Duración temporal de arreglo
N = tMax*Fs;    % Número de muestras en arreglo
f = 20;         % Frecuencia de oscilación en [Hz]
Amp = 0.8;      % Amplitud de onda

t5 = (0:(N-1))/Fs;  % Arreglo de tiempo

y5 = Amp*sin(2*pi*f*t5);    % Ahora incluimos frecuencia y amplitud

% Graficamos
figure(1)
plot(t5,y5)
xlabel('Tiempo(s)')
ylabel('Nivel (Pa)')
title('Onda Senoidal')
grid on
box on
xlim([0 tMax])
ylim([-1.1 1.1])

%% Probar el teorema nyquist-shannon

% Demostración de ALIASING

for f = 5:5:100
    
    fprintf('Frec = %d [Hz]\n',f)
    
    Fs = 100;     % Frecuencia de muestreo
    T = 1/Fs;       % Intervalo
    tMax = 1;       % Duración temporal de arreglo
    N = tMax*Fs;    % Número de muestras en arreglo
    % f = 20;         % Frecuencia de oscilación en [Hz]
    Amp = 0.8;      % Amplitud de onda
    
    t5 = (0:(N-1))/Fs;  % Arreglo de tiempo
    
    y5 = Amp*sin(2*pi*f*t5);    % Ahora incluimos frecuencia y amplitud
    
    clf     % Limpia el gráfico en la figura
    
    % Graficamos
    figure(2)
    plot(t5,y5)
    xlabel('Tiempo(s)')
    ylabel('Nivel (Pa)')
    title('Onda Senoidal')
    grid on
    box on
    xlim([0 1])
    ylim([-1.1 1.1])
    
    pause       % Pausa la ejecución del programa hasta orden de usuario
    
end