Отправлено andy_P 04 июля 2006 г. 15:29 В ответ на: Чего-то не получается (+) отправлено <font color=gray>homekvn</font> 04 июля 2006 г. 15:16

Этап 5 должен выгладеть как
new_X1 = 1/16*ifft(fft(X1).*fft(sinc));

и все у Вас получится :-)

Реализовывать это в лоб не надо. Проще будет посчитать именно циклическую свертку

• Спасибо. Понял. Действительно так работает. Только это верно Вы ранее заметили, что сворачивать с синком в частотной области накладнее будет, чем в лоб сделать иффт, затем добавить нули, затем опять ффт. Это выгоднее будет даже несмотря на прореженность ЧХ, которую мы сворачиваем. (Только что сказанное справедливо для больших длин фильтров) — 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, пустое)

E-mail: info@telesys.ru