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

Спасибо! понял, поиграюсь. Благо, реальные данные уже есть в матлабе, есть материал для экспериментов.

Про FFT без децимации сразу на весь массив- сам бы так хотел. Но это у меня просто еще одна дополнительная опция в приборе, желательно бы обойтись без сотрясения основ уже сделанного. А CMSIS-овская уже используемая библиотека только FFT4096 умеет, вот и приходиться децимировать.
Мне бы вообще-то очень бы подошло что-то типа FFT длиной 16К или даже 32К точек float32, но это уже свой код применять вместо встроенных библиотек, про это тоже думал. Снижение быстродействия даже в 5 раз еще приемлемо, по сравнению с CMSIS. Тоже интересно подумать.

Попробуйте, я не смотрел, как это будет оптимизироваться на АРМ'е, но на DSP ложится прилично.
 FFT32k.zip ( 1010 байт ) Кол-во скачиваний: 47

Вы будете таки смеяться, но именно этот код я использую уже лет 5, про него и говорил.
Работает с небольшими вариациями в разных системах, начиная от мелкого 8-битного Майкрочипа, где данные храню во внешней микросхеме сериальной RAM. Но именно в АРМе я его тоже не пробовал.
Причем, действительно, работает всегда и везде. Нашел в Интернете.
Ну раз и Вы его тоже рекомендуете, значит он действительно неплох, не буду искать что-то еще.
Спасибо!

Интерполяцию окном Hanna смотри http://calculator2006.narod.ru/articles/haan.htm

Сообщение отредактировал ivan219 - Jun 16 2017, 14:06

отлично, спасибо!!!!
Скопировал себе на диск и пдф-ку сделал, а то в интернете бывает так- сегодня есть, а завтра -"404"

То есть расчет идет по трем точкам.
И странно там в тексте как-то окно задано, формула 1-cos(2pi*n/N) правильная? Она же двойку дает в максимуме, даже с рисунком в том же тексте не совпадает.
Думаю, просто опечатка, должно быть все-таки 0.5*(1-cos(2pi*n/N)), собственно так окно в первоисточнике и задается.

Будет спектр чуть выше.
И да скорей всего опечатка я у же и не помню почему так вышло.

Боюсь, что для Вашей жуткой точности не хватит интнрполяции по трем точкам. Лучше бы 5-7 точек или даже больше.