Отправлено andy_P 03 июля 2006 г. 12:15 В ответ на: Ну блоки у меня и так будут короче, чем ИХ фильтра, причем в несколько раз. Только это не решает проблемы. Может, я не понял чего (в частности, для меня непонятным является термин "паддинг данных") отправлено <font color=gray>homekvn</font> 03 июля 2006 г. 12:06

Делаем фильтрацию так:

1. H = FFT(ИХ) - получаем N точек импульсной характеристики. N - это длина ИХ без Вашего паддинга.
2. блоки данных берем длиной N/2, добиваем 0 (паддинг) Получаем блоки длиной N
3. FFT данных c паддингом (Т точек), умножаем на H
4. IFFT. Дальше overlap-save или overlap-add

При перестройке фильтра в частотной области получим те же N точек ИХ (нет нужды рассчитывать ЧХ фильтра таким образом, чтобы последние N отсчетов ИХ получались 0 как в Вашем случае)

Вроде бы так. Если наврал, то поправьте :-)

• Не могу солгаситься в одной детали, которая является существенной. — homekvn   (03.07.2006 14:29 212.185.161.237, 3559 байт)
  • Ответ+ — andy_P   (03.07.2006 18:06 80.82.63.185, 945 байт)
    • Хм-м! Про четные и нечетные отсчеты я в основном посте написал. Про свертку там же тоже упомянул, но проблема была в том, с чем сворачивать (+) — homekvn   (03.07.2006 20:02 212.185.161.237, 1907 байт)
      • Ответ+ — andy_P   (03.07.2006 20:43 80.82.63.185, 199 байт)
        • Но из этого можно что-то выжать только для определенного вида ИХ фильтра. Вобщем зря я это... Вы правы, ничего хорошего из этого не получится. — andy_P   (03.07.2006 21:09 80.82.63.185, пустое)
          • Жалко, что так. А как заманчиво выглядит обертка! Синтезировал фильтр в частотной области (прямо на лету) и тут же пропустил через него сигнал все в той же частотной области... Ан нет! Накаси-выкуси! Что ж буду брать иффт, брать буду ее действительную часть, затем снабжать все это добро нулями и обять ффтить... :-( — homekvn   (03.07.2006 21:27 212.185.161.237, пустое)
            • Либо свертка с синком в частотной области, чтобы окно на ИХ во временной наложить. Синк подусечь можно, но тогда алиасить будет. Кстати, если синк недлинный взять, то вычислений немного потребуется. Хотя о чем я :-) Ведь FFT c немеренной разрядностью у Вас в этом проекте? — andy_P   (03.07.2006 21:56 80.82.63.185, пустое)
              • А почему с синком сворачивать надо? Там же хрен знает сворачивать надо, вот с ней: fft(exp(-j*2pi*n/(2*N))) при n=0..N-1. А она ж вроде как не очень синк?... Да и правы Вы, будь она даже синком, укорачивать бы ее я не стал - опасно бо зело. — homekvn   (03.07.2006 22:05 212.185.161.237, пустое)
                • С синком прореженную нулями ЧХ сворачивать - эквивалентно наложению прямоугольного окна во временной области. — andy_P   (04.07.2006 10:47 80.82.63.185, пустое)
                  • Вроде так. Но тогда получается, что это fft(exp(-j*2*pi*n/(2*N))) n=0..N-1 синк и есть, а я в матлабе поглядел - не очень-то он синк напоминает, скорее совсем не напоминает. — homekvn   (04.07.2006 11:49 212.185.161.237, пустое)
                    • Ответ:+ — andy_P   (04.07.2006 12:31 80.82.63.185, 522 байт)
                      • Там немножко по-другому получается (+) — homekvn   (04.07.2006 13:04 212.185.161.237, 1027 байт, ссылка)
                        • Ответ: Можно и так. То что я написал нужно делать для разбавленной нулями длинной выборки. У Вас нечетные (или четные??? запутался) коэффициентя получаются сверткой с H0. Результат я думаю будет аналогичен.+ — andy_P   (04.07.2006 13:37 80.82.63.185, 193 байт)
                          • Если просто добавить нулей сзади, то четные отсчеты БПФ от такого сигнала совпадают с отсчетами БПФ исходного сигнала, в который мы нули не добавляли. Нечетные отсчеты и есть моя головная боль. Но то, что Вы написали это интересно (+) — homekvn   (04.07.2006 13:50 212.185.161.237, 502 байт)
                            • Я просто предлагаю проинтерполировать ЧХ H0 чтобы получить H. Если выбрать интерполирующий фильтра таким, как я написал, то отсчеты ИХ = ifft(H) будут равны 0 для n>N/2. А четные отсчеты H будут совпадать с отсчетами H0 — andy_P   (04.07.2006 14:15 80.82.63.185, 1 байт)
                              • Чего-то не получается (+) — homekvn   (04.07.2006 15:16 212.185.161.237, 1039 байт)
                                • Ответ:+ — andy_P   (04.07.2006 15:29 80.82.63.185, 182 байт)
                                  • Спасибо. Понял. Действительно так работает. Только это верно Вы ранее заметили, что сворачивать с синком в частотной области накладнее будет, чем в лоб сделать иффт, затем добавить нули, затем опять ффт. Это выгоднее будет даже несмотря на прореженность ЧХ, которую мы сворачиваем. (Только что сказанное справедливо для больших длин фильтров) — homekvn   (04.07.2006 15:51 212.185.161.237, пустое)
                                    • Тут единственный момент в том,что то с чем сворачиваете усекать можно и проверить, может рипплы в ИХ от N/2 до N будут приемлемыми+ — andy_P   (04.07.2006 16:01 80.82.63.185, 276 байт)
                                      • + — andy_P   (04.07.2006 18:51 80.82.63.185, 109 байт)
                                        • Еще не проверял, но, чувствую, в аудиообработке не пройдет - слышно будет. Проще уж фильтр в раза в полтора-два по длине сократить. — homekvn   (04.07.2006 18:56 212.185.161.237, пустое)
                                          • Hi-end у Вас немеренный получается :-) В коммуникашке все проще и сложнее - таких требований по разрядности нет, зато полосы широкие шибко :-) — andy_P   (04.07.2006 19:28 80.82.63.185, пустое)
                                            • А я-то причем? Мне сказали пилить - я пилю. По мне так вообще КИХ-фильтры для такого приложения, как у меня, не нужны. До селе прекрасный звук и при помощи БИХ-ов в кабине машины получался. А заказчик - увы - непереубеждаем. — homekvn   (04.07.2006 20:32 212.185.161.237, пустое)
                                              • Надеюсь, влезете в DSPшку со всем этим хозяйством. Успехов. — andy_P   (04.07.2006 20:57 80.82.63.185, пустое)
                    • Думаю, что так будет выглядеть "периодический" синк. — homekvn   (04.07.2006 12:13 212.185.161.237, пустое)
                • Ну а fft у меня на плавучке будет. — homekvn   (03.07.2006 22:06 212.185.161.237, пустое)
            • Тк ИХ у Вас действительная, то FFT половинной длины можно использовать. Все легче. — andy_P   (03.07.2006 21:48 80.82.63.185, пустое)
              • Да, эт-точно, можно. А можно и два канала в одном флаконе обработать, сделав один комплексный сигнал из двух действительных каналов. — homekvn   (03.07.2006 22:00 212.185.161.237, пустое)
        • тьфу, real(fft(exp(-j*2*pi/N*(0:N/2-1))))) = ones(1,N/2). — andy_P   (03.07.2006 20:46 80.82.63.185, пустое)
          • Да, результат действительно интересный. Раньше такого не замечал. Просто не приходило в голову рассматривать действительную часть от fft. Обычно меня только модуль интересует. Только вот поможет ли это?... — homekvn   (03.07.2006 21:06 212.185.161.237, пустое)

E-mail: info@telesys.ru