Доброе время суток.
Непонятно какими критериями пользоваться при выборе разрядности переменных D0, D1, D2

Допустим для x выбираем 10 бит, для k - 16 бит, разрядность D0?
Интуиция подсказывает, что надо учитывать размер блока N, где бы про это почитать?

сформулируйте вопрос свой. А то совершенно не понятно о чем речь идет?

Да всё, алгоритм уже закодирован и работает.




Применяли на практике алгоритм Герцеля? Выражение в первом сообщении реализует вычисление мощности одной частоты (подробности об алгоритме есть в гугле), проблема была в определении необходимой разрядности переменных сдвигового регистра (d0, d1, d2), так сказать методика.

A16 + B16 -> C17 ?

N = N1 + N2; // N - итоговая разрядность суммы, N1 - разрядность исходных операторов
N2 = Log2(S)+1; // S - число сложений

А как учесть множитель k, который по модулю может быть до 2? Просто добавить ещё один бит? Или "+1" во второй формуле (N2 = Log2(S)+1) это и учитывает?

Как раз это и учитывает: Герцель = вычисление одонго бина Фурье. Увеличение разрядности в фурье как раз log2(s), где S - длина фурье. соответственно на коэффициент k закладываем ещё 1 бит и получаем увеличение разрядности log2(s)+1. Исходная разрядность N1, тогда N=N1+log2(S)+1, всё как изложено выше.

Алгоритм Герцеля - модификация рекурсивного ДПФ.
А в нем ошибки накапливаются неслабо даже с плавающей запятой.
Поэтому Герцеля нельзя использовать для измерения спектра, а только в таких
триггерных схемах, как ДТМФ.
И значит, формулу ошибки с логарифмом надо умножить надвое.

Откуда такая информация?