Для анализа качества метод FFT-IFFT после обратного преобразования нужно брать прямое преобразование но с размером в 4 или да же 16 раз большем. Тогда будет видна истин.

Это все пропускающий фильтр с 1 коэффициентом усиления. Но не линейной ФЧХ

Теория http://vunivere.ru/work2336/page4

Vlad27 писал(а): Рассчитать их по билинейному преобразованию нельзя ввиду
нереальности такого аналогового прототипа с единичной
передачей в бесконечность.

Вот аналоговый прототип.http://calculator2006.narod.ru/FileInet/fv.gif
Ищи по фразе фазовращатели.

Если известны частоты то вы можете воспользоваться Герцелем в этом случае вы получите максимум точности. Правда с потерей производительности. Но для ПК и такого приложения это не критично.

Если вы хотите найти длину периода то рекомендую все же искать не через вершины а через 0.

Сделайте децимацию.
Тогда длинна периода станет меньше.
Но до децимации нужно сделать фильтрацию что бы не было наложения.

С аналоговыми фильтрами и частотой дискретизации расчёт простой. Чем выше частота дискретизации тем меньший порядок аналогового фильтра нужен для заданной полосы. Самый плохой вариант это когда крайняя частота вашего полезного сигнала находится рядом с частотой Найквиста. В этом случае порядок анало...

Вы сделали не фильтр!
Вы сделали -> генератор сигнала и умножили его на что то с начало подумал оконная функция но и на неё да же не похоже.

Все делают обычно так:
int i = 0; i < N - 1; i++

Интересные у вас результаты.
Не задумывались написать статью по всему что у вас наработано?
И опубликовать в журнале.

А вы по отдельности всё посчитайте! (Cos(a*t) * Cos(w*t)) - (Cos(a*t) * Sin(w*t))
При умножении вещественного сигнала на комплексный вы получите комплексный сигнал разностной и суммарной частоты.