Добро пожаловать, гость
:: алгоритмы  и методы :: :: олимпиадные задачи :: :: связь :: :: о сайте :: :: форум ::

Форум работает в режиме архива, только для чтения и поиска.
Архив 2004 Архив 2007 Архив 2013

 
 
Опции темы Поиск в этой теме Опции просмотра
  #1  
Старый 03.12.2007, 09:20
Новичок

Отправить личное сообщение для SEreGA Посмотреть профиль Найти все сообщения от SEreGA
 
Регистрация: 03.12.2007
Сообщений: 1

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

Отправить личное сообщение для MBo Посмотреть профиль Найти все сообщения от MBo
 
Регистрация: 21.09.2006
Адрес: Новосибирск
Сообщений: 1,374

странный эффект
Вот, ввожу 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]));
Изображения:
Тип файла: jpg fft.jpg (17.8 Кб, 240 просмотров)

Последний раз редактировалось MBo, 03.12.2007 в 11:04.
 


Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск
Опции просмотра


Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
Вопрос по ДСТ NepsteR Математические алгоритмы (другое) 1 21.07.2007 18:19
Вопрос по сплайну (интерполяция) Jek Вычислительная геометрия 3 21.06.2007 08:13
вопрос по обработке изображения tumanovalex Обработка изображений, звук, графика 0 10.04.2007 13:32
вопрос по обработке. Роман Обработка изображений, звук, графика 1 17.03.2007 17:16
вопрос по фурье преобразованию. Kate Ovechkina Математические алгоритмы 22 07.12.2006 09:00