reverberacion.m

来自「eco and reverberation matlab codes」· M 代码 · 共 19 行

function y=reverberacion(frase,D,a,necos)

[x,fs,nbits]=wavread(frase);
x=x';   %traspuesta de x porque wavread te da las muestras en columnas


for i=1:necos
 
x1=zeros(1,D);  %Creo un vector de ceros de 1 a D(atraso de la se馻l)
x2=[x,x1];      %concateno los dos vectores
ax=a*x;         %multiplico la constante por cada valor del vator
x3=[x1,ax];     %concateno el vectro de ceros con el de la atenuacion 
y=x2+x3;        %primer eco
x=y;
end

sound(x)