%% visualizacion sonido

clearvars;
clear all;
close all;
clc;

%% sonido

%hacemos sonar la frecuencia de 1000 Hz durante 2 segundos
tini = 0;
tfin = 2;
fs = 44100;
dt = 1/fs;
t = tini:dt:tfin;
x = 0.75*sin(2*pi*1000*t);
sound(x,fs);

figure(1)
plot(t,x);
title('sonido de 1000 Hz durante 2 segundos')
xlabel('tiempo t(s)');
ylabel('presion sonora p (Pa)');
legend('p(t)');
axis([min(t), max(t), -1.1*max(x), 1.1*max(x)]);
grid on;
box on;

%% modelo de propagacion sonora en una dimension
%la presion sonora sera una funcion del tiempo y del espacio
%es decir p = p(x,t) modelo de onda plana
%sin embargo para que este modelo tenga sentido debemos hacerlo
%en camara lenta

%variables principales
c = 1.0;        %velocidad del sonido en verdad c = 344 (m/s)
f = 2.0;        %frecuencia de 2 (Hz)
w = 2*pi*f;     %frecuencia angular
T = 1/f;        %periodo
l = c/f;        %longitud de onda
k = 2*pi/l;     %numero de onda
A = 1.0;

%espacio lineal
x = 0:0.01:5;     %espacio entre 0 a 5 metros
t = 1;            %tiempo fijo 1 segundo


for t = 0:0.01:2
    %presion sonora
    pxt = A*cos(w*t - k*x);

    %presion sonora x0 = 2(m)
    x0 = 2.0;
    px0t = A*cos(w*t - k*x0);

    figure(2)
    plot(x,pxt,'b',x0,px0t,'r o');
    title('presion sonora en el espacio');
    xlabel('espacio x(m)');
    ylabel('presion sonora p(Pa)');
    legend('p(x,t)');
    axis([0,5,-1.25,1.25]);
    grid on;
    box on;

    figure(3)
    hold on;
    plot(t,px0t,'r o');
    hold off;
    title('presion sonora en el tiempo');
    xlabel('tiempo x(s)');
    ylabel('presion sonora p(Pa)');
    legend('p(x,t)');
    axis([0,2,-1.25,1.25]);
    grid on;
    box on;


    pause(0.01);    %pausa por 0.01

end;