вопрос по БПФ

[SEreGA]

Помогите разобраться с проблемой.
Делаю БПФ некого *.wav файла. В нем записаны несколько частот, длительностью по 3 секунды 1000, 1400, 2100 и тд, между частотами тишина. Частота дискретизации звука записанного файла 8кГ . Проблема в том что когда использую расчет по 8192 точкам получаю АЧХ с зеркальной составляющей относительно 2кГ. (т.е. на спектре я имею 2 пика (1000Гц и 3000Гц), (1400 и 2600), (2100 и 1900). При расчете по 64, 128, 256, 512, 1024, 2048 и 4096 точках таких проблем нет спектр имеет один основный пик и гармоники.
Подскажите с чем это связано и как можно решить эту проблему? есть ли алгоритм пересчета? Частоту дискретизации увеличивать до 16кГ нет возможности.
Заранее благодарен!

[MBo]

странный эффект
Вот, ввожу 8192 точки, три серии с частотами 1000, 1400, 2100 Гц, получаю ожидаемую картинку:

var
  Re, Im: array of Single;
  i, N: Integer;
  Amp: Double;
begin
  N := 8192;
  SetLength(Re, N);
  SetLength(Im, N);
  for I := 0 to N div 4 - 1 do begin
    Re[i] := Sin(1000 * i * 2 * Pi / N);
    Re[i + N div 4 + 500] := Sin(1400 * i * 2 * Pi / N);
    Re[i + 2 * (N div 4) + 1000] := Sin(2100 * i * 2 * Pi / N);
  end;
  for I := 0 to N - 1 do begin
    Im[i] := 0;
    Series1.AddXY(i, Re[i]);
  end;
  FFT1D(Re, Im, N, 1);
  for I := 0 to N - 1 do
    Series2.AddXY(i, Hypot(Re[i], Im[i]));