Re: Инверсный фильтр

From

Evgenij Masherov (2:5020/175.2)

To

Igor Rudym

Date

2002-12-27T12:56:47Z

Area

RU.ALGORITHMS

From: "Evgenij Masherov" <EMasherow@nsi.ru>

Fri Dec 27 2002 11:28, Igor Rudym wrote to Evgenij Masherov:

 
 IR>>>>>>> фомула следующая
 IR>>>>>>> W(w) = G(w)/K(w)

 IR>>>>>>> где  G = FFT(g); g = k*y;
 IR>>>>>>> y - начальный вектор, k - орератор преобразования (например  смаза)

 EM>>>> Вкратце - умножение ПФ от вектора Х на ПФ от вектора У
 EM>>>> эквивалентно вычислению свертки этих векторов и взятию от нее
 EM>>>> ПФ. Но свертки циклической, т.е. за последним элементом этого
 EM>>>> вектора оказывается первый.

 IR>>> И как это применительно к моему вопросу? А есть хоть како-нить 
 IR>>> примерчик в инете? Как брать дискретную циклическую
 IR>>> свертку? Свертка - это итеграл от умножения X*Y?

 EM>> Свертка - это сумма вида
 EM>> z[j]=SUM x[i]*y[j-i]
 EM>> Связь ее с Фурье в том, что Фурье от свертки есть произведение Фурье
 EM>> от свертываемых векторов.

 IR> Понятно, но у меня не произведение, а деление G/K. тогда свертка будет:
 IR> z[j] = SUM g[i]/k[j-i]
 IR> если так, но дело в том, что k - содержит нули.

 Ну, деление есть умножение на обратное...
Но что до нулей - наличие их означает, что при прямом преобразовании из
сигнала что-то полностью ушло, и восстановить невозможно. Если есть точный
нуль - то просто на него не делят. Намного интереснее ситуация, когда вместо
точного нуля у нас ненулевое значение, полностью обусловленное ошибкой. Тогда
деля на это, как правило, малое число, мы усиливаем ошибку чрезвычайно.
Здесь нужно регуляризацию применить...

Евгений Машеров АКА СанитарЖеня

--- ifmail v.2.15dev5
 * Origin: FidoNet Online - http://www.fido-online.com (2:5020/175.2)