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

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

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

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

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

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

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

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

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

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