Отправлено ASergej_R19 27 сентября 2006 г. 19:13

Ничего не получается...
В продолжение http://www.telesys.ru/wwwboards/dsp/262/messages/38179.shtml

Действительно обнаружил предельный цикл, все попытки что-то сделать - ни к чему не привели... Я безнадежен?

Простейший IIR HP Баттерворта 2-го порядка на частоте дискретизации 1500 Гц имеет срез -3 дБ на 0.4 Гц. Коэффициенты получил 1.31, на вход данные 24.8 - арифметика 32-х разрядной точности фиксированная точка и не получается - шумы прут (в основном - в виде постоянной составляющей) в зависимости от частоты естесственно, но выходят "на полку" в районе -60 дБ, что не приемлемо...
Его вообще возможно такой построить используя 32-х разрядную арифметику и получить цифровой SNR за -120 дБ? Тут есть какая-то профессиональная хитрость?

Конечно советы использовать 64 битную арифметику и/или плавующую точку - хорошие, но загрузку процессора жалко...
Строить этот HP аналоговым способом - тоже не интересно - хочется разобраться с цифровой обработкой...

Заранее всем спасибо за советы и терпение...

• Я уж писал об этом. Попробую еще раз кинуть эту идейку (+) — homekvn   (27.09.2006 21:17 212.185.161.237, 338 байт)
  • У меня на входе постоянной составляющей нету - подавал идеальный синус... (-) — ASergej_R19   (28.09.2006 11:19 80.250.160.170, пустое)
    • А... пардон. Я Ваш предыдущий вопрос перепутал с другим как две капли воды похожим на Ваш. Но тогда (+) — homekvn   (28.09.2006 11:54 212.185.161.237, 685 байт)
      • Я так и делал... Данные 24 бит, привел их к 24.8... Видимо косяк в реализации, ибо BH и SM написали, что с самим фильтром все ок... Будем искать (с)... (-) — ASergej_R19   (28.09.2006 12:02 80.250.160.170, пустое)
        • Не, то, о чем Вы написали, это не совсем то, о чем писал я. Я не предлагал увеличивать разрядность, если Ваш сигнал удовлетворяет некоторым условиям (его амплитуда не превышает определенного значения). — homekvn   (28.09.2006 12:54 212.185.161.237, пустое)
          • Допустим, если макс. амплитуда не превышает 2^18, то можно все сделать в имеющейся разрядности. Для этого перед фильтрацией надо увеличить амплитуду в 2^5 раз, далее профильтровать, а потом выход сдвинуть вправо на 5 разрядов. Постоянная составляющая тогда ликвидируется. — homekvn   (28.09.2006 12:58 212.185.161.237, пустое)
            • А, понял... Нет, к сожалению амплитуда сигналов может быть любой в пределах этих 24 бит... Но - спасибо за идею... (-) — ASergej_R19   (28.09.2006 13:21 80.250.160.170, пустое)
              • Ответ: Если у Вас данные занимают все 24 разряда, то где запас по динамическому диапазону для переменных состояния?(+) — Vadim Kudryavtsev   (29.09.2006 13:35 59.11.26.53, 661 байт)
                • Ответ (+) — homekvn   (29.09.2006 14:01 212.185.161.237, 116 байт)
                  • Но вообще, да (+) — homekvn   (29.09.2006 14:05 212.185.161.237, 211 байт)
• Интересно, откуда люди берут на входе данный 24.8, при этом подспудно рассматривая возможность их аналоговой фильтрации?:-) — -=ВН=-   (27.09.2006 19:44 193.125.71.140, пустое)
  • Я просто сумбурно немного написал... Извините... (+) — ASergej_R19   (27.09.2006 19:53 80.250.160.170, 189 байт)
    • И сколько же разрядность Вашего АЦП? — SM   (27.09.2006 19:59 213.208.180.51, пустое)
      • АЦП Cirrus CS5381, 24-бит дельта-сигма... (-) — ASergej_R19   (27.09.2006 20:02 80.250.160.170, пустое)
    • Да они все простые, фильтры. А единственный путь уменьшить "шум" - увеличивать разрядность. Хотите Вы этого или нет. Если конечно нет потусторонних граблей в реализации. — -=ВН=-   (27.09.2006 19:59 193.125.71.140, пустое)
      • Ну тогда бы у меня вообще не работал бы ни один фильтр... (+) — ASergej_R19   (27.09.2006 20:09 80.250.160.170, 146 байт)
        • Ну раз в косяке сильно сомневаетесь, то вперед, к повышенной разрядности. — -=ВН=-   (27.09.2006 20:13 193.125.71.140, пустое)
  • Да и получить "цифровой SNR" в -120дб аналоговым фильтром не видится возможным. Он даст "аналоговый SNR" :) — SM   (27.09.2006 19:48 213.208.180.51, пустое)
    • А под "цифровым SNR" имелось ввиду шумы вычислений фильтра... (-) — ASergej_R19   (27.09.2006 19:55 80.250.160.170, пустое)
• А какая структура фильтра? Тут такое дело (+) — SM   (27.09.2006 19:41 213.208.180.51, 349 байт)
  • Да... Извините, забыл написать... (+) — ASergej_R19   (27.09.2006 19:50 80.250.160.170, 290 байт)
    • И что, никак не получается понять, на каком этапе вычислений эта гадость образуется? Пустите модель пошагово, да найдите того, кто переполняется или наоборот, слишком обрезает результат. — SM   (27.09.2006 19:58 213.208.180.51, пустое)
      • Проблема в том, что коэффициенты у фильтра видимо... (+) — ASergej_R19   (27.09.2006 20:14 80.250.160.170, 232 байт)
        • Да косяк у Вас. Фильтр нормальный. Не поленился проверил, правда на плав. зап. На целых лень. — -=ВН=-   (27.09.2006 21:09 193.125.71.140, 281 байт)
          • Да вот и с целыми проверил. Те же 2 синуса, но 24 разрядных, формат 1.23. Промеж. и вых. данные в формате 2.30, т.е. 32 разряда. "Шумы" в районе -120 дб. — -=ВН=-   (27.09.2006 21:44 193.125.71.140, пустое)
            • М-да... Странно... Пошел искать косяк... Огромное спасибо... (-) — ASergej_R19   (28.09.2006 11:22 80.250.160.170, пустое)
        • А как Вы их к 1.31 умудрились привести? Если они явные 2.30? — SM   (27.09.2006 20:29 213.208.180.51, пустое)
          • вынес 2 в scale фактор... (+) — ASergej_R19   (27.09.2006 20:34 80.250.160.170, 128 байт)
            • Если завтра будет время, поглядю поподробнее на сей фильтр — SM   (27.09.2006 20:38 213.208.180.51, пустое)
              • Бью челом! Хоть одним глазком... (+) — ASergej_R19   (27.09.2006 20:43 80.250.160.170, 280 байт)
                • Фильтрует, только в путь (+) — SM   (27.09.2006 21:36 213.208.180.51, 1035 байт)
                  • Большое спасибо... Странное дело... Пошел искать ошибку у себя... (-) — ASergej_R19   (28.09.2006 11:23 80.250.160.170, пустое)