|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Evgenij Masherov 2:5020/175.2 27 Dec 2002 13:56:47
To : Igor Rudym
Subject : Re: Инверсный фильтр
--------------------------------------------------------------------------------
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>>>> ПФ. Hо свертки циклической, т.е. за последним элементом этого
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 - содержит нули.
Hу, деление есть умножение на обратное...
Hо что до нулей - наличие их означает, что при прямом преобразовании из
сигнала что-то полностью ушло, и восстановить невозможно. Если есть точный
нуль - то просто на него не делят. Hамного интереснее ситуация, когда вместо
точного нуля у нас ненулевое значение, полностью обусловленное ошибкой. Тогда
деля на это, как правило, малое число, мы усиливаем ошибку чрезвычайно.
Здесь нужно регуляризацию применить...
Евгений Машеров АКА СанитарЖеня
--- ifmail v.2.15dev5
* Origin: FidoNet Online - http://www.fido-online.com (2:5020/175.2)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
