Отправлено andy_P 05 сентября 2005 г. 14:24 В ответ на: Помогите советом... /+/ отправлено <font color=gray>Николай</font> 02 сентября 2005 г. 09:31

Кушает ресурсы в основном только первая секция (та что крутится на наибольшей частоте дискретизации) - звеньев там много не надо - достаточно подавить после Fd_next/2 - Fd_next - частота дискретизации следующей ступени. Так как фильтр широкий, то по фазе он не сильно нагадит в полосе 0..55 Hz. Остальные ступени можно делать что ких, что бих - много ресурсов не съедят, если сначала скажем раза в 4 частоту дисретизации уменьшите.

• Ответ: Минус один - каждая децимация будет мало-мало алиасить в полосу. Те чем меньше стадий децимации, тем лучше — andy_P   (05.09.2005 16:54 80.82.63.185, пустое)
  • В принципе, да. Но алиазинга можно избежать, если сделать децимацию только после первого биквада, а после второго и третьего оставить все как есть. Эти второй и третий биквады весь алиазинг предыдущего шага и "сожрут". Возможны и иные схемы, когда децимация делается не по "полной программе" — homekvn   (05.09.2005 17:06 212.185.161.237, пустое)
    • Да вроде не сожрут что в полосу 0..55 Hz попало? — andy_P   (05.09.2005 17:42 80.82.63.185, пустое)
      • Так в полосу 0..55 и не попадут, поскольку я предлагаю децимацию сделать только после первого каскада и, скажем, не в 58 раз, как это в принципе можно было бы сделать, доведя частоту дискретизации до 110 Гц, а только в 10-12 раз. Тогда алиазинг в полосе 0..55 будет ничтожно малым, поскольку за первым каскадом будут следовать два каскада без децимации вообще. И вот эти два каскада и "сожрут" алиазинг от первого каскада почти полностью. — homekvn   (06.09.2005 11:07 212.185.161.237, пустое)
• Можно и ких. Только все-таки MIPS-ов и памяти он больше скушает (+) — homekvn   (05.09.2005 16:10 212.185.161.237, 761 байт)
  • Ответ: дак там полоса 55 hz на частоте дискретизации > 6kHz. Делаем срез первой ступени 100. Имеем ~4 октавы или 24 dB с секции (если децимировать на 2 первый раз). Те 2 секций по моему тупому разумению на первой ступени должно хватить за глаза. Наверное, можно уложиться и в одну секцию, если как-о извратиться — andy_P   (05.09.2005 16:42 80.82.63.185, пустое)
    • Да КИХ без прореживания через БПФ на такой-то частоте на калькуляторе можно успеть посчитать... Ловкость рук, так сказать... — -=ВН=-   (05.09.2005 17:09 194.190.181.231, пустое)
      • Ни-ни (+) — homekvn   (05.09.2005 17:43 212.185.161.237, 1335 байт)
        • Ну считайте. — -=ВН=-   (05.09.2005 18:03 194.190.181.231, 561 байт)
          • Согласен, что упорядочивать после Фурье не нужно. Но на мой взгляд задача не стоит того, чтобы ее через Фурье решать. Кстати, Фурье, на мой взгляд, надо будет 32-разрядное взять. Иначе после FFT и IFFT шумы округления будут достаточно большими. — homekvn   (06.09.2005 10:59 212.185.161.237, пустое)
            • Для 16-ти разрядного проц. его с бл. пл. запятой можно сделать, 16-ти разрядным. — -=ВН=-   (06.09.2005 12:04 194.190.181.231, пустое)
  • Про память да, памяти КИХ сожрет больше. — -=ВН=-   (05.09.2005 16:22 194.190.181.231, 39 байт)
    • Ну, согласен. MIPSы с учетом прореживания, а также того, что при вычислении КИХ достаточно 16-бит, будут такими, что то на то и выйдет. Но задержка... Сам ких имеет неслабую задержку (около 300 семплов на тактовой частоте), да еще и многокаскадность внесет свой вклад. А если еще после этого апсемплить надо?.. — homekvn   (05.09.2005 16:45 212.185.161.237, пустое)
      • Какую задержку внесет многокаскадность? Есть фильтр с какой-то пер. ф-ей. И все. Никого не интересует как он реализован. Многкаскадно или немногокаскадно. — -=ВН=-   (05.09.2005 17:05 194.190.181.231, 137 байт)
        • Не совсем так. Это несложно показать (+) — homekvn   (05.09.2005 17:25 212.185.161.237, 1169 байт)
          • Ну и посчитайте теперь, какой фильтр, но однокаскадный, будет полностью эквивалентен Вашим 2 каскадам. Будет ли его длина 59 отчетов. — -=ВН=-   (05.09.2005 17:35 194.190.181.231, пустое)
            • В том-то все и дело, что прежде, чем писать, посчитал. Все так, как написал, и есть. Иначе бы и писать не стал (+) — homekvn   (05.09.2005 17:49 212.185.161.237, 394 байт)
              • Да что Вы посчитали, если эквив. длина Ваших 2 секций=92 отсчетам. В пересчете на исх. частоту дискр. — -=ВН=-   (05.09.2005 17:53 194.190.181.231, пустое)
                • Да, Ваша правда. Прошу прощения за бузеш. Если сделать фильтр с таким же подавлением, то получается примерно одинаковая задержка (+-1 семпл). — homekvn   (06.09.2005 10:52 212.185.161.237, пустое)
                • Прошу прощения. Завтра отвечу. Сейчас ухожу. — homekvn   (05.09.2005 18:05 212.185.161.237, пустое)
  • Сорри (я подумал, что ких, а написано бих). Беру пред. слова назад — homekvn   (05.09.2005 16:12 212.185.161.237, пустое)