flanging.m

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

function y=flanging(frase,D,a)      %tengo q meterle como se馻l de entrada un seno.
[x,fs,nbits]=wavread(frase);    %si es un seno esto no ace falta 
x=x';                           %y esto tampoco
y=x(1:D);
axv=0;
for i=1:length(x)-D
    %x1=zeros(1,D);  %Creo un vector de ceros de 1 a D(atraso de la se馻l)   
    %x2=[x1,x];      %ceros + se馻l
    ax=round((D/2)*(1+cos(2*pi*i/(length(x)-D))));       %multiplico la constante por cada valor del vator
    
    %x3=[x1,ax];     %concateno el vectro de ceros con el de la atenuacion 
    %y=x2+x3;
    %x=y;
    y=[y,x(i+D)+a*x(i+D-ax)];
end
sound(x)