Разложение числа на простые множители:
1600 = 2 · 2 · 2 · 2 · 2 · 2 · 5 · 5
Микс БПФ по двойке и пятерке.

+1.
Подробнее можно посмотреть здесь: Блейхут Р. Быстрые алгоритмы цифровой обработки сигналов.

О вариантах БПФ
Parallel.ru

Еще:
Морфизприбор - БПФ для ГАС:
 NontraditionalFFT.pdf ( 186.41 килобайт ) Кол-во скачиваний: 119

К слову - нетрадиционные БПФ, это сейчас - нетрадиционные. Тогда, в 70-х традиционными были скорее нетрадиционные сегодня, т.к. цена быстродействия за лишнюю точку была высока.

P.S.
Тогда, когда еще не все были повально увлечены БПФ по степеням двойки, достаточно интересным и более быстрым был алгоритм Винограда. Представляя одномерное БПФ как многомерное, с разложением по простым множителям, удавалось значительно повысить скорость обработки.
 PNAS_1976_Winograd_1005_6.pdf ( 237.93 килобайт ) Кол-во скачиваний: 75

Годная статья но маленькая, сейчас "нетрадиционное ффт" называется mixed radix fft. Вот например


А вот такой вопрос топикстартеру, скажем первая гармоника (50 Гц) изменяется на 0.1 %, тогда 50я гармоника изменится примерно на пол бина FFT и чего Вы тогда намеряете?

Для ТС-ра, даже и это - слишком много, надеюсь, что - пока.

P.S.
Родной мой, ну ты хватил - давай еще CRAY и пр. сюда подтащим.
ТС-ру нужна реализация на микроконтроллере, какие уж тут ядра и параллельности.

В 4G, он же LTE, есть размерность FFT, равная 1536, так что там 3^1*2^9. И сейчас есть нетрадиционные размерности

Не пугайте юношу страшным словом "Метрология".

Ага для CRAY есть аналогичное . На самом деле ценность математики, в этих статьях - это возможность записать любое FFT в виде
одной формулы, а там CRAY или Intel I-7 или АRM это уже вторично.

Для тех, кто в такой теме давно - может и да.
Но, зачем ТС-ра пугать в его теме?

И вправду, если нужны значения только на определенных частотах (гармоники) нафига бпф делать? Считайте мощность только в нужных точках (сдвигайте спектр на столько на сколько надо герц и осредняйте).
И про разрешающую спотобность тоже правильно упомянули тут, хотя так же правильно сказали, что это не наше дело.
В целом общая идея комментариев, что интерполяция тут не нужна.

Задайте частоту дискретизации АЦП 6400 Гц. За 128 выборок наберете 20 мс. Вычислив БПФ, получите 64 гармоники 50 Гц сигнала, с полосой каждого бина 50 Гц.

Потом для большей точности измерения амплитуды усредняйте результаты БПФ.

Не можете уйти от 8 кВыб/с? Интерполируйте на 4, децимируйте на 5. Хотя я этого садомазо из-за упертости схемотехника (или кого?) не понимаю.

Да, АЧХ каждого бина не плоская, на краях загибается. Думаю, можно забить. Или взять 256 выборок, БПФ по 256, получить шаг 25 Гц. Лучше 512 или больше, чтобы знать, как относить палки спектра к нужным гармоникам сигнала.

10 раз усреднённый спектр по 20мс и один по 200мс совсем не одно и то же.

А то!

По-хорошему, я бы вообще делал десять замеров периодов, начиная каждый замер в момент прерывания по zero-detected.
Мало ли там, что - по ГОСТУ-у

Я бы вообще не считал (количество замеров)... Усреднял бы и усреднял непрерывно замеры.