|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Evgenij Masherov 2:5020/175.2 25 Dec 2002 10:55:29
To : Igor Rudym
Subject : Re: Инверсный фильтр
--------------------------------------------------------------------------------
Tue Dec 24 2002 15:42, Igor Rudym wrote to Evgenij Masherov:
IR>>> фомула следующая
IR>>> W(w) = G(w)/K(w)
IR>>> где G = FFT(g); g = k*y;
IR>>> y - начальный вектор, k - орератор преобразования (например
IR>>> смаза)
IR>>> Hо это не работает! Даже когда k берется следующим образом:
IR>>> k(x-s) = 1, если x-s=0 и 0, иначе.
EM>> Это известные штуки, связанные с тем, что БПФ считает данные
EM>> циклическими. Hужно делать сие через свертку с наложением.
Гольденьберг, Матюшкин и Поляк. Цифровая обработка сигналов.
Отнес и Эноксон.
Бендат и Пирсол.
Вкратце - умножение ПФ от вектора Х на ПФ от вектора У эквивалентно вычислению
свертки этих векторов и взятию от нее ПФ. Hо свертки циклической, т.е. за
последним элементом этого вектора оказывается первый. При работе "в лоб" это
порождает неприятности на границах. Поэтому вектора удлинябт вдвое, дополняя
нулями, так что за последним настоящим идут нули, а на первый при свертке
наползают они же. При этом свернутый вектор оказываетя вдвое длиннее
требуемого, и используется сумма первой половины на данном шаге и второй
половины на предыдущем (а вторая половина данного - запоминается для
следующего).
Евгений Машеров АКА СанитарЖеня
--- ifmail v.2.15dev5
* Origin: FidoNet Online - http://www.fido-online.com (2:5020/175.2)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
