Кодирование деконволюции с использованием Python

Question

Прежде чем начать, я должен сказать вам, что у меня нулевые знания о DSP в Python.
Я хочу деконволюцию двух звуковых сигналов, используя python, чтобы я мог извлечь импульсную характеристику помещения, входной сигнал был синусоидальным, а выходной записал его.
Я написал кусок кода, но он не работал, я пытался слишком долго и действительно безрезультатно.
Может кто-нибудь, пожалуйста, помогите мне с кодом, который вычисляет БПФ ввода и вывода, а затем вычислить h iFFT их доли и построить его.

Answer 1

Деконволюция является некорректной сложной проблемой при наличии шума и пространственно-вариативного размытия. Я предполагаю, что у вас нет пространственно-вариативной проблемы, поскольку вы используете FFT, поэтому вы можете использовать модуль восстановления из пакета python skimage (вместо программирования алгоритма на низком уровне с FFT).

Здесь вы можете изучить пример кода с одним из реализованных методов в модуле восстановления.

Я рекомендую вам прочитать O'Leary et al. Книга, если вы хотите узнать больше. У всех авторов этой книги есть более продвинутые книги на эту великую тему.