reverberacion.asv

来自「eco and reverberation matlab codes」· ASV 代码 · 共 26 行

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


%sound(y)

    x4=zeros(1,necos*D-1);  %Creo un vector de ceros de 1 a D(atraso de la se馻l)*necos -1
                            %porque ya tiene un retardo anterior.    
    x5=[y,x4];      %la salida anterior la concateno con x4.
    x6=[x4,y];     %concateno el vectro de ceros con el de la atenuacion 
    y2=x5+x6;
    sound(y2)