Практическая часть. Определим функцию, задающую так называемый пилообразный сигнал
Определим функцию, задающую так называемый пилообразный сигнал

и изобразим ее на графике


Заполним массив s:
.
Проводим прямое преобразование Фурье:
.
Внимание! В том случае, когда в массиве s содержится
элементов, причем все числа действительные, следует использовать функцию fft. Во всех остальных случаях – функцию cfft. Массив g содержит комплексные коэффициенты дискретного преобразования Фурье.
Размер массива
-
.
Для анализа вклада отдельных гармоник в исходный сигнал изобразим на графике модули и аргументы гармоник


Проводим обратное преобразование Фурье, исключив гармоники с малым вкладом. Будем учитывать только гармоники с амплитудой не менее 0.3. Для отсечения слагаемых с малым вкладом воспользуемся функцией единичного скачка – функцией Хевисайда Ф

Для обратного преобразования Фурье используется функция ifft, если прямое преобразование осуществлялось с помощью fft, и cifft, если прямое преобразование осуществлялось с помощью cfft


Повторим преобразование Фурье, учтя слагаемые с амплитудой до 0.1.


Учет дополнительных гармоник существенно улучшил результат синтеза сигналов.