Programa

fs1=10000;

z1=wavrecord(2*fs1,fs1);

x=z1;

fs2=5000;

z2=wavrecord(2*fs2,fs2);

h=z2;

y=conv(x,h);

subplot (2,1,1),stem (y)

% Loop para evaluar la sumatoria

acumulador=0;

for n = 1:length(x)

      acumulador = acumulador + x(n);

end

acumulador

%Usando una matriz, para llenar

clear all

clc

%1. Definicion de las entradas y de la respuesta al impulso del filtro a implementar.

e=[4 5 2 3 -5 6];

h=[3 4 0 1 2 6 7 ];

%2. creando matriz resultado en ceros

m=length(h)+length(e)-1;

n=length(h);

w=zeros(m,n);

y=zeros(1,n);

W=zeros(1,length(e)+length(h)-1);

a=1;

%3. Realizando convolucion y llenado de matriz

for k=1:m

    for b=n:-1:2

        y(b)=y(b-1);

    end

    if k<=(m-n+1)

        y(1)=e(k);

    else

        y(1)=0;

    end

    for i=1:n

        w(k,i)= h(i)*y(i);

        a=a+1 ;

    end

    a=1;

end

for j=1:m

    xnuevo=0;

    xant=0;

    for c=1:n

        xant=w(j,c);

        xnuevo=xant+xnuevo;

    end

    W(j)=xnuevo;

end

%Señal de salida

Solucion_Codigo=W;

Solucion_Matlab=conv(e,h);

display(Solucion_Codigo)

display(Solucion_Matlab)