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

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

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

 

Заполним массив s: .

Проводим прямое преобразование Фурье: .

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

Размер массива - .

Для анализа вклада отдельных гармоник в исходный сигнал изобразим на графике модули и аргументы гармоник

 

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

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

 

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

 

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