Вещественное значение частоты Опции

[nameless]

Прошу сильно не пинать.
Есть вектор отсчетов, полученных на интервале, в который укладывается пусть 0.1 - 10 периодов гармонического сигнала. То есть некратная (вещественная) частота (вернее оценка частоты) нас интересует. После дискретного преобразования Фурье ищем центр тяжести:
Ц = (SUMM(i=0 to [N/2]-1) of i*W(i)^2] / SUMM(i=0 to [N/2]-1) of W(i)^2] (у Гоноровского подсмотрел такой способ нахождения центра тяжести). Здесь W = sqrt(A^2 + B^2), A = SUMM(i=0 to N/2 - 1) of X(i)*sin(2*pi*i*K/N) B = ... понятное дело с косинусом. К = 0, 1, 2 ... N/2 - 1 целочисленный(!) Так вот ошибка при моделировании в районе Ц меньше 2 достигает 100%. Дальше - меньше, но не устраивает никак. Ткните носом пожалуйста.

[ivan219]

Да Y2 и X2 это коорденаты максимального бина. X1, Y1 и X3, Y3 это боковые бины. Соответственно (-1: 0: 1)

Упрощенный вареант.

Код

Result := (Y2 - Y1) / (2 * Y2 - Y1 - Y3) + X2 - 0.5;

X2, Y2 коорденаты максимального бина Y1, Y2 значения боковых бинов.

Сообщение отредактировал ivan219 - Mar 27 2011, 08:20