Как найти ядро ​​свертки в частотной области?


12

У меня есть два вектора пространственных данных (каждый около 2000 элементов в длину). Один - это свернутая версия другого. Я пытаюсь определить ядро, которое произвело бы такую ​​свертку. Я знаю, что могу сделать это путем нахождения обратного преобразования Фурье отношения преобразований Фурье выходного и входного векторов. Действительно, когда я делаю это, я получаю более или менее ту форму, которую я ожидал. Тем не менее, мой вектор ядра имеет ту же размерность, что и два входных вектора, когда на самом деле свертка использовала только одну пятую (~ 300-400) точек. Тот факт, что я получаю правильную форму, но неверное количество очков, заставляет меня думать, что я не совсем правильно использую функции ifft и fft. Кажется, если бы я действительно поступал правильно, это должно было произойти естественно. На данный момент я просто делаю;

FTInput = fft(in); 
FtOutput = fft(out);
kernel = ifft(FtOutput./FTInput).

Это правильно, и я должен правильно интерпретировать выходной вектор, или я упростил задачу? Я уверен, что это последнее, я просто не уверен, где.


2
Ваши входные данные дополняются нулями с обеих сторон до длины ядра свертки? Должно быть, иначе вы потеряете информацию, которая может быть причиной этих артефактов.
оставил около

Ответы:


5

Если в вашем сигнале присутствует шум, простое деление области Фурье вызовет множество ошибок в вашем результате. Некоторые способы избежать этого - использование так называемого двухканального БПФ ( часть 1 и часть 2 ). Я также могу предложить деконволюцию с помощью адаптивных фильтров, в частности, фильтры LMS или NLMS ([наименьшие] наименьшие средние квадраты), которые легко понять, не слишком дороги с точки зрения циклов ЦП в случае длинных сигналов. Адаптивные фильтры LMS очень устойчивы к шуму.


Это сработало, спасибо, и я обнаружил совершенно новый тип вещей, о которых я не знал.
Боулер

@Phonon эти ссылки работают? Что это за подпись точно?
Спейси

@Mohammad: Это для доступа к статьям и книгам технического обзора Brüel & Kjær. Регистрация бесплатна и дает доступ ко многим хорошим статьям.
Тор
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.