%% programa solución de ecuacion de onda 02

clearvars;
close all;
clc;

%este programa grafica y anima la solucion de la ecuacion de onda
%mediante su solucion usando serie de fourier

%% parametros
%% datos inciales
L = 1;          %longitud de la cuerda
A = 1;          %amplitud
c = 1;          %velocidad de propagacion de ondas
b = 0.1;        %amortiguamiento

x = 0:L/100:L;  %vector de distancia
Nx = length(x); %longitud de vector x

%% ciclo principal de animacion

N = 10;          %numero de terminos de la sumatoria

p(1:Nx) = 0;   %incializacion de presion

for t = 0:0.01:5
    
    p(1:Nx) = 0;   %incializacion de presion
    
    for n = 1:N
      
        kn = (2*n -1)*pi/(2*L);
        bn = b/(2*c^2);
        wn = (2*n -1)*pi*c/(2*L);
        xin = bn/wn;
        wd = wn*sqrt(1 - xin);
        Bn = (8*L)/( ((2*n-1)*pi)^2*c*sqrt(1-xin) );
        Bn = Bn*( sin( (2*n+3)*pi/8 ) - sin( (6*n+1)*pi/8  ) );
               
        %p(1:Nx) = p(1:Nx) + Bn*exp(-bn*t)*sin(wd*t)*sin(kn*x);
        p = p + Bn*exp(-bn*t)*sin(wd*t)*sin(kn*x);
    end;
    
    figure(1)
    plot(x,p)
    title('onda sonora tubo');
    xlabel('x (m)');
    ylabel('p(x,t) (Pa)')
    legend(['t = ',num2str(t),' (s)']);
    axis([0, L, -1.1, 1.1])
    grid on;
    box on;
    
    pause(0.01);
    
end;