function Z=kfb(v)
N=1024;

f3db=40;
f=-f3db:2*f3db/N:f3db;
S=exp(-(1.665*f/f3db).^2);
s=ifft(S(1:N));

s1=0:0.1:1;
y=[0    0.0034    0.0134    0.0302    0.0538    0.0844... 
  0.1224    0.1682    0.2278    0.3536    1.0000];
p1=polyfit(y,s1,4);
r=polyval(p1,s);
pp=r.^10;


R=fft(r);
RF=abs(R);
a0=RF(1,1);
RF=RF./a0;
HW1=sqrt(RF);

PP=fft(pp);
PPF=abs(PP);
a1=PPF(1,1);
PPF=PPF./a1;
HW2=sqrt(PPF);

 
v1=randn(1,1024);
v2=randn(1,1024);
v3=randn(1,1024);
v4=randn(1,1024);
v5=randn(1,1024);
v6=randn(1,1024);
v7=randn(1,1024);
v8=randn(1,1024);
v9=randn(1,1024);
v10=randn(1,1024);

v11=randn(1,1024);
v12=randn(1,1024);

V1=fft(v1);
V2=fft(v2);
V3=fft(v3);
V4=fft(v4);
V5=fft(v5);
V6=fft(v6);
V7=fft(v7);
V8=fft(v8);
V9=fft(v9);
V10=fft(v10);

V11=fft(v11);
V12=fft(v12);

W1=HW1.*V1;W2=HW1.*V2;W3=HW1.*V3;W4=HW1.*V4;W5=HW1.*V5;
W6=HW1.*V6;W7=HW1.*V7;W8=HW1.*V8;W9=HW1.*V9;W10=HW1.*V10;

W11=HW2.*V11;W12=HW2.*V12;

w1=ifft(W1);w2=ifft(W2);w3=ifft(W3);w4=ifft(W4);w5=ifft(W5);
w6=ifft(W6);w7=ifft(W7);w8=ifft(W8);w9=ifft(W9);w10=ifft(W10);
w11=ifft(W11);w12=ifft(W12);

w1=w1.^2;w2=w2.^2;w3=w3.^2;w4=w4.^2;w5=w5.^2;
w6=w6.^2;w7=w7.^2;w8=w8.^2;w9=w9.^2;w10=w10.^2;
w11=w11.^2;w12=w12.^2;

Y=w1+w2+w3+w4+w5+w6+w7+w8+w9+w10;
X=w11+w12;

z=X+Y;
z=sqrt(z);

y=xcorr(z);

Z=fft(y(1:N+1));

Z=10*log10(abs(Z(1:N/2+1))/sum(abs(Z)));
plot(Z);