function bh92lobe = genbh92lobe()
% ==> generation of the main lobe of the spectrum of a Blackman-Harris 92dB window
bh92size = 4096;
bh92lobe = zeros(bh92size,1);
bh92N = 512;
bh92const = [.35875, .48829, .14128, .01168];
bh92Theta = -4*2*pi/bh92N;
bh92ThetaIncr = 8*2*pi/bh92N/bh92size;
for i=1:bh92size
    for m=0:3
        bh92lobe(i)=bh92lobe(i)-bh92const(m+1)/2*...
            (sine2sine(bh92Theta-m*2*pi/bh92N,bh92N)+...
            sine2sine(bh92Theta+m*2*pi/bh92N,bh92N));
    end;
    bh92Theta = bh92Theta + bh92ThetaIncr;
end;
bh92lobe = bh92lobe/bh92lobe(bh92size/2+1);

function x = sine2sine(x, N)
% sine2sine function !!!
x = sin((N/2)*x) / sin(x/2);