deconvolution и нули

[ushahid]

Ако се опитаме да deconvolve квадратна вълна или отгоре шапка от всякакви функции, ние винаги са изправени пред проблема за разделение на нула.Има ли решение за предотвратяване на този проблем?
например
а = триъгълен импулс
B = Top Hat
H =?(система за отговор)

за deconvolution
А = FFT (а)
B = FFT (б)

H = реалната (ifftshift (ifft (A / B)));

сега, тъй като Б е sinc функция, така че има нули в него, която не позволява правилното функциониране deconvolution.Всички предложения за избягване на този проблем дели на нула и да получите по-разумна преценка на нашата система за отговор?

 

[vadkudr]

Можете да използвате deconvolution в момента домейн.В момента операторът домейн извивка е представена от toeplitz матрица.Deconvolution изисква invertion на тази матрица.Това може да стане по-бърз начин с помощта на Lewinson-Дърбин алгоритъм със сложност O (N ^ 2).
По начина, по който съществува по-бързо (макар и малко сложна) алгоритъм със сложност O (N * log2 (N) ^ 2)

виждам
R. Кумар.Бърз алгоритъм за решаване на toeplitz система от уравнения.IEEE Trans.на Acous., реч и Регистрация.Proc., 33:254 - 267, февруари, 1985

MATLAB и реализация на моята начална страница
http://vadkudr.boom.ru/Algorithms/vadkudr_Algorithms_eng.html