Python: Частотный анализ звуковых файлов

Я создаю несколько звуковых файлов, которые воспроизводят тона на разных частотах с определенным количеством гармоник.
В конечном итоге эти звуки будут воспроизводиться на устройстве с небольшим динамиком.

У меня есть кривая частотной характеристики динамика, и я хочу сделать на Python следующее:

1. Построить частотный спектр звукового файла. Мне нужно взять БПФ файла и построить его с помощью gnuplot
2. Применить нелинейную передаточную функцию на основе кривой частотной характеристики в таблице данных.
3. Постройте результат после применения функции.

] Кто-нибудь знает:

• Какой самый простой способ сделать это?
• или Приложение (основанное на GNU / Linux), которое могло бы сделать это за меня?