Реализация алгоритмов дискретного преобразования Фурье на языке программирования Python | Статья в журнале «Молодой ученый»

Отправьте статью сегодня! Журнал выйдет 21 декабря, печатный экземпляр отправим 25 декабря.

Опубликовать статью в журнале

Автор:

Рубрика: Технические науки

Опубликовано в Молодой учёный №20 (154) май 2017 г.

Дата публикации: 20.05.2017

Статья просмотрена: 2265 раз

Библиографическое описание:

Кузин Д. А. Реализация алгоритмов дискретного преобразования Фурье на языке программирования Python // Молодой ученый. — 2017. — №20. — С. 39-43. — URL https://moluch.ru/archive/154/43570/ (дата обращения: 08.12.2019).



В настоящее время обработкой звукового сигнала занимаются множество специалистов: инженеры, звукорежиссеры, звукооператоры, музыканты. Все они преследуют цель синтеза звукового сигнала и (или) изменения параметров ЗС для выполнения инженерных (например, применение кодеков для сжатия звуковых сигналов) и художественных (например, улучшение звучания композиции путем тщательной настройки отдельных частот) задач.

Существует множество аппаратных и программных средств для синтеза и обработки звукового сигнала. Существует как платное ПО (Adobe Audition), так и бесплатное ПО (Audacity, Sonic Visualizer) для работы со звуком. Ядром (математической основой) данных программных средств является дискретное преобразование Фурье, а также его частный случай — быстрое преобразование Фурье. Кроме этого, используются различные особенности и свойства преобразования Фурье, например, симметричность, заполнение нулями, смещение.

Цель.

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

Результат работы.

Была осуществлена программная реализация свойств преобразования Фурье в Python.

Будем использовать «ipython», который представляет собой оболочку, которая очень удобна для интерактивного вызова команд python. Также используем «matplotpib» (для импорта при старте ipython пишем в терминале: «ipython-pylab»). Скрипт всё же удобно набирать в текстовом редакторе. В данном случае мы используем уже знакомый нам пакет «numpy», а также пару новых: треугольная функция (triang из scipy.signal) и реализацию дискретного преобразования Фурье в python (fft из scipy.fftpack). Если длина последовательности fft не будет равна 2n то вместо быстрого преобразования Фурье (fft) будет выполняться дискретное преобразование фурье (dft). Именно таким образом работает эта команда. Выполнив преобразование Фурье последними двумя строками кода получаем частотную и фазовую характеристики. Данная программа не выводит сразу никаких графиков, она считает все переменные. Построить графики можно вводя, например:

plot(x) — выводит входную функцию (треугольник)

plot(X) — выводит график преобразования Фурье

plot (mX) — выводит график частотного спектра

plot (pX) — выводит график фазового спектра

Рассмотрим график частотного спектра треугольного сигнала (рисунок 1). Он симметричен вокруг нулевого отсчета. Первая часть — это положительные значения (до 8 отсчета), а вторая часть — отрицательные. Это значит, что четырнадцатый отсчет является по факту также отсчетом минус первым. Данный график очень ясно дает понять ранее описанное свойство ДПФ — симметричность.

Рис. 1. Частотный спектр треугольного сигнала

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

− fftbuffer=np.zeros(15) — создаем буфер, заполненный нулям;

− fftbuffer [:8]=x [7:] — перемещаем конец треугольного сигнала в первую часть буфера;

− fftbuffer [8:]=x [:7] — перемещаем начало треугольного сигнала в последнюю часть буфера;

− X=fft(fftbuffer) — считаем преобразование Фурье буфера, заполненного значениями функции треугольника.

Построив график фазового спектра треугольника (который теперь расположен от -7 до 7) мы получаем симметричный график, который и должен был получиться (значения на рисунке в -15 степени, что очень близко к нулю) (рисунок 2):

Рис. 2. Частотный спектр треугольной функции

Теперь, используя практически тот же скрипт, но с добавлением функции чтения звуковых файлов, посчитаем преобразования Фурье реального звукового сигнала.

В файл скрипта заносим изменение в виде добавления функции «wavread», которая рассматривалась ранее в других скриптах, а также добавляем код, позволяющий центрировать звуковой файл. Таким образом получим скрипт, который выполняет ДПФ и БПФ:

Рис. 3. Скрипт, позволяющий осуществлять ДПФ с целью получения частотной и фазовой характеристики реального сигнала

Команды библиотеки «matplotlib» (plot(x) — построить форму сигнала, plot(x1) — построить выбранную для анализа часть сигнала) позволяют получить графики, которые необходимы нам для анализа. Например, после выполнения скрипта, запросим фазовую и частотную характеристики, введя plot(mX) и plot(pX). Также, чтобы визуализировать свойство заполнения нулями мы можем построить ДПФ буфер, используя plot(fftbuffer). На рисунке 4 приведены построенный с помощью скрипта python графики, демонстрирующие сразу несколько свойств преобразования Фурье — симметричность (на частотной характеристике и фазовой) и заполнения нулями (на графике буфера ДПФ).

Рис. 4. а) фрагмент ЗС, б) частотный спектр (включая положительную и отрицательную части), в) буфер ДПФ, г) фазовый спектр

Чтобы фазовый спектр был более читаем, применим развертывание фазы. Для этого изменим строку с расчетом фазовой характеристики:

pX=np.unwrap(np.angle(X))

И построим фазовую характеристику. Рисунок 5 изображает фазовую характеристику нашего входного сигнала в удобном для анализа виде.

Рис. 5. Фазовая характеристика звукового сигнала

Таким образом удалось реализовать несколько свойств и теорем Фурье в виде Python скрипта, чтобы лучше понять концепции симметрии, заполнения нулями, развертывания фазы. Данный скрипт позволяет выполнять Фурье анализ фрагмента звукового сигнала.

Литература:

  1. Официальный сайт Соник Визуализатор для скачивания и получения документации. Режим доступа: http://www.sonicvisualiser.org/
  2. Официальный сайт проекта Audacity для скачивания и получения документации. Режим доступа: http://www.audacityteam.org/
  3. Официальный сайт пакета iPython. Режим доступа: http://www.ipython.org/
  4. Электроакустика и звуковое вещание: Учебное пособие для вузов / Э. И. Вологдин, А. П. Ефимов, А. А. Фадеев, И. А. Алдошина, Г. П. Катунин, Л. Н. Кацнельсон, Ю. А. Ковалгин — 2007. – 872c.
  5. Цифровая обработка сигналов в трактах звукового вещания: Учебное пособие для вузов / О. Б. Попов, С. Г. Рихтер — 2007. – 341 с.
  6. Акустика: Справочник / А. П. Ефимов, А. В. Никонов, М. А. Сапожников, В. И. Шоров — 1989. – 336 с.
  7. Сайт Wikipedia. Discrete Fourier Transform. Режим доступа: https://en.wikipedia.org/wiki/Discrete_Fourier_transform
  8. Сайт Wikipedia. Fast Fourier Transform. Режим доступа: https://en.wikipedia.org/wiki/Fast_Fourier_transform
  9. Официальный сайт компании National Instruments предоставляющей платформы для проведения исследований, а также программный продукт LabVIEW. Режим доступа: http://www.ni.com/ru-ru.html
Основные термины (генерируются автоматически): фазовая характеристика, частотный спектр, звуковой сигнал, свойство преобразования, треугольный сигнал, дискретное преобразование, быстрое преобразование, треугольная функция, фазовый спектр, заполнение нулями.


Похожие статьи

Алгоритмы преобразования Фурье и их применение при анализе...

Дискретное преобразование Фурье (ДПФ) (DiscreteFourierTransform, DFT) имеет вид: (2). ДПФ ставит в соответствие N отсчетам сигнала x(n), где n = 0…N — 1, N отсчетов дискретного спектра X(k), при этом предполагается, что и сигнал, и спектр сигнала являются...

Исследование процесса цифровой обработки сигнала при работе...

Ключевые слова: дискретное Преобразование Фурье, спектр сигнала, белый шум, импульс.

Быстрое преобразование Фурье (БПФ, FFT) — алгоритм быстрого вычисления дискретного преобразования Фурье (ДПФ).

Преобразование Фурье и преобразование Хартли

Кроме того, используя собственные свойства, преобразование введенное Хартли позволяет произвести расчёты косвенных вычислений спектра мощности Фурье данной функции, пользуясь только арифметикой.

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

Дискретное преобразования Фурье (ДПФ) на конечной области для функции одной переменной , где определяется как [1]: , Данное преобразование является основным в частотной области улучшения изображений.

Контроль структуры магнитного поля МПФС ЛБВ методом...

Если исследуемый сигнал задан в виде дискретных значений, то используют дискретное преобразование Фурье (ДПФ). Для уменьшения времени расчета ДПФ используют алгоритм быстрого преобразования Фурье (БПФ) [2, 3].

Аналого-цифровое преобразование | Статья в журнале...

аналого-цифровое преобразование, входной сигнал, непрерывная функция, аналоговый сигнал, процедура преобразования, цифровой код, цифровая форма, входной аналоговый сигнал, аналого-цифровой...

Генератор звуковых сигналов аддитивного синтеза...

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

‒ Вынесены в качестве коэффициентов функции такие важные параметры сигнала

Преобразование Фурье: [Электронный ресурс]: [веб-сайт] — Режим доступа: https...

Предварительная обработка речевых сигналов для системы...

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

Алгоритмы преобразования Фурье и их применение при анализе...

Дискретное преобразование Фурье (ДПФ) (DiscreteFourierTransform, DFT) имеет вид: (2). ДПФ ставит в соответствие N отсчетам сигнала x(n), где n = 0…N — 1, N отсчетов дискретного спектра X(k), при этом предполагается, что и сигнал, и спектр сигнала являются...

Исследование процесса цифровой обработки сигнала при работе...

Ключевые слова: дискретное Преобразование Фурье, спектр сигнала, белый шум, импульс.

Быстрое преобразование Фурье (БПФ, FFT) — алгоритм быстрого вычисления дискретного преобразования Фурье (ДПФ).

Преобразование Фурье и преобразование Хартли

Кроме того, используя собственные свойства, преобразование введенное Хартли позволяет произвести расчёты косвенных вычислений спектра мощности Фурье данной функции, пользуясь только арифметикой.

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

Дискретное преобразования Фурье (ДПФ) на конечной области для функции одной переменной , где определяется как [1]: , Данное преобразование является основным в частотной области улучшения изображений.

Контроль структуры магнитного поля МПФС ЛБВ методом...

Если исследуемый сигнал задан в виде дискретных значений, то используют дискретное преобразование Фурье (ДПФ). Для уменьшения времени расчета ДПФ используют алгоритм быстрого преобразования Фурье (БПФ) [2, 3].

Аналого-цифровое преобразование | Статья в журнале...

аналого-цифровое преобразование, входной сигнал, непрерывная функция, аналоговый сигнал, процедура преобразования, цифровой код, цифровая форма, входной аналоговый сигнал, аналого-цифровой...

Генератор звуковых сигналов аддитивного синтеза...

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

‒ Вынесены в качестве коэффициентов функции такие важные параметры сигнала

Преобразование Фурье: [Электронный ресурс]: [веб-сайт] — Режим доступа: https...

Предварительная обработка речевых сигналов для системы...

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

Похожие статьи

Алгоритмы преобразования Фурье и их применение при анализе...

Дискретное преобразование Фурье (ДПФ) (DiscreteFourierTransform, DFT) имеет вид: (2). ДПФ ставит в соответствие N отсчетам сигнала x(n), где n = 0…N — 1, N отсчетов дискретного спектра X(k), при этом предполагается, что и сигнал, и спектр сигнала являются...

Исследование процесса цифровой обработки сигнала при работе...

Ключевые слова: дискретное Преобразование Фурье, спектр сигнала, белый шум, импульс.

Быстрое преобразование Фурье (БПФ, FFT) — алгоритм быстрого вычисления дискретного преобразования Фурье (ДПФ).

Преобразование Фурье и преобразование Хартли

Кроме того, используя собственные свойства, преобразование введенное Хартли позволяет произвести расчёты косвенных вычислений спектра мощности Фурье данной функции, пользуясь только арифметикой.

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

Дискретное преобразования Фурье (ДПФ) на конечной области для функции одной переменной , где определяется как [1]: , Данное преобразование является основным в частотной области улучшения изображений.

Контроль структуры магнитного поля МПФС ЛБВ методом...

Если исследуемый сигнал задан в виде дискретных значений, то используют дискретное преобразование Фурье (ДПФ). Для уменьшения времени расчета ДПФ используют алгоритм быстрого преобразования Фурье (БПФ) [2, 3].

Аналого-цифровое преобразование | Статья в журнале...

аналого-цифровое преобразование, входной сигнал, непрерывная функция, аналоговый сигнал, процедура преобразования, цифровой код, цифровая форма, входной аналоговый сигнал, аналого-цифровой...

Генератор звуковых сигналов аддитивного синтеза...

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

‒ Вынесены в качестве коэффициентов функции такие важные параметры сигнала

Преобразование Фурье: [Электронный ресурс]: [веб-сайт] — Режим доступа: https...

Предварительная обработка речевых сигналов для системы...

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

Алгоритмы преобразования Фурье и их применение при анализе...

Дискретное преобразование Фурье (ДПФ) (DiscreteFourierTransform, DFT) имеет вид: (2). ДПФ ставит в соответствие N отсчетам сигнала x(n), где n = 0…N — 1, N отсчетов дискретного спектра X(k), при этом предполагается, что и сигнал, и спектр сигнала являются...

Исследование процесса цифровой обработки сигнала при работе...

Ключевые слова: дискретное Преобразование Фурье, спектр сигнала, белый шум, импульс.

Быстрое преобразование Фурье (БПФ, FFT) — алгоритм быстрого вычисления дискретного преобразования Фурье (ДПФ).

Преобразование Фурье и преобразование Хартли

Кроме того, используя собственные свойства, преобразование введенное Хартли позволяет произвести расчёты косвенных вычислений спектра мощности Фурье данной функции, пользуясь только арифметикой.

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

Дискретное преобразования Фурье (ДПФ) на конечной области для функции одной переменной , где определяется как [1]: , Данное преобразование является основным в частотной области улучшения изображений.

Контроль структуры магнитного поля МПФС ЛБВ методом...

Если исследуемый сигнал задан в виде дискретных значений, то используют дискретное преобразование Фурье (ДПФ). Для уменьшения времени расчета ДПФ используют алгоритм быстрого преобразования Фурье (БПФ) [2, 3].

Аналого-цифровое преобразование | Статья в журнале...

аналого-цифровое преобразование, входной сигнал, непрерывная функция, аналоговый сигнал, процедура преобразования, цифровой код, цифровая форма, входной аналоговый сигнал, аналого-цифровой...

Генератор звуковых сигналов аддитивного синтеза...

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

‒ Вынесены в качестве коэффициентов функции такие важные параметры сигнала

Преобразование Фурье: [Электронный ресурс]: [веб-сайт] — Режим доступа: https...

Предварительная обработка речевых сигналов для системы...

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

Задать вопрос