%% animacion

clear all;
close all;
clc;

%este programa animara la grafica de funciones simples
%para ello necesitaresmos arreglos unidimensionales (N X 1)
%x,y,z que inicialmente podran identificarse como espacio
%ademas necesitamos una variable escalar (1 X 1) que es el tiempo t

%% inicializacion

%prealocacion - servira para acelerar los calculo
N = 1001;       %numero de puntos
x = zeros(N,1); %prealocamos la variable x (N X 1) para el espacio
p = zeros(N,1); %prealocamos la variable p (N X 1) para la presion sonora


%% repaso

%antes del calculo principal recordemos que el comando for
%tiene esta estructura general

%for var = valor_inicial : incremento : valor_final
%       instrucion 1
%       instrucion 2
%       instrucion 3
%       ...........
%       instrucion N
%end;

%% valores iniciales

x = 10*(0:N-1)/N;       %espacio
f = 1;                  %frecuencia
w = 2*pi*f;             %frecuencia angular
c = 1;                  %velocidad del sonido
k = w/c;                %número de inda
lamda = 2*pi/k;         %longitud de onda
A = 1.5;                %amplitud


%% calculo principal

%inicializamos el tiempo

for t = 0: 0.1 : 4*pi %inicio de ciclo del tiempo (t)
    

    %variable presion sonora p = (x,t)
    p = A*sin( w*t + k*x );  %onda plana simple
    
    
    %graficamos a cada iteración del tiempo se cambia
    %la figura generando la animación
    
    figure(1)
    plot(x,p);
    title('onda plana 1')
    xlabel('x (m)');
    ylabel('p(x,t) (m)');
    grid on;
    box on;
    axis([0 10 -2 2]);
    
    pause(0.01);
    
end;                  %fin del ciclo del tiempo