Взять кусочек 4096-точечного DFT в окрестности максимума 128-точечного FFT?, Есть ли смысл... Опции

[tmtlib]

Как насчёт следующего алгоритма:
1) Берём FFT 128 точек, грубо находим максимум
2) В окрестности максимума запускаем DFT с шагом, как у FFT на 4096 точки, но только для маленького диапазона частот. Т.е. получаем кусок спектра из 50-60 точек с повышенной точностью.

Есть смысл? Как это будет соотноситься с поиском центра масс только по FFT-128?

[fontp]

Как насчёт следующего алгоритма:
1) Берём FFT 128 точек, грубо находим максимум
2) В окрестности максимума запускаем DFT с шагом, как у FFT на 4096 точки, но только для маленького диапазона частот. Т.е. получаем кусок спектра из 50-60 точек с повышенной точностью.

Есть смысл? Как это будет соотноситься с поиском центра масс только по FFT-128?

1.А как вы его возьмёте? У вас обычно нет данных для этого... Добавить нулей? Можно, уже много раз обсуждалось...
Или Вы просто не хотите делать большой БПФ4096, а хотите поиметь только часть частот около максимума?? А данные добавлять всё равно нулевые...Тоже возможно, сделать на нужных частотах Герцеля, но не БПФ, поскольку это одно и то же по результату. В любом случае потом по точкам вблизи максимума подгоняют параболу (или даже функцию спектрального окна, но это не для реального времени)


2.Вот Вам простой и самый точный из известных мне способ интерполяции вблизи максимума, если без дополнительных преобразований. Значительно точнее чем по центру масс. Но как и в случае с добавлением нулей будет содержать небольшие систематические ошибки. Т.е. при очень высоком SNR будет работать теоретически не оптимально

Y is a three-element complex vector with the
% DFT output magnitude maximizer as the center element Y(kmax-1), Y(kmax), Y(kmax+1)).
%
% Returns -0.5 < x < 0.5, which is the fraction of the sample
% spacing (i.e., bin width) about the center element where the
% peak is estimated to be.

% Ref: Macleod, M.D., "Fast Nearly ML Estimation of the Parameters
% of Real or Complex Single Tones or Resolved Multiple Tones,"
% IEEE Trans. Sig. Proc. Vol 46 No 1,
% January 1998, pp141-148.

ref = Y(2); % Isolate phase reference.
R = real(Y.*conj(ref)); % Generate phase corrected coefficient vector.
gamma = (R(1)-R(3))/((2*R(2))+R(1)+R(3)); % Calculate offset.
delta = (sqrt(1 + 8*gamma*gamma)-1)/(4*gamma); % Final estimate.
x = delta;

Можете найти полжение максимума по этой формуле, потом прогнать Герцель в окрестности. Герцель это фильтр, который в точности даёт результат БПФ. но только в одной частотной точке. Можно найти в первом приближении центральную частоту по Macleod, прогнать Герцель на этой частоте f0 и двух соседних f0-0.25df, f0+0.25df, df шаг БПФ на 128. Потом провести параболу и посмотреть на какой частоте её вершина. Называют это авторы ML-extension

[tmtlib]

1.А как вы его возьмёте? У вас обычно нет данных для этого... Добавить нулей? Можно, уже много раз обсуждалось...
Или Вы просто не хотите делать большой БПФ4096, а хотите поиметь только часть частот около максимума??

Да, я хотел делать именно так:
БПФ128 находит максимум X,
потом идёт цикл ДПФ4096, но i пробегает значения не 0..4095, а (X*32-10)..(x*32+10)
В обоих случаях ничего нулей не добавляю.

пс Вот смотрю приведённый вами интересный метод.