caxapa.ru :: Норот, каг-то не парился, а тут, вдруг, заметил, что BPF IIR странно реагирует на постоянку на входе. Он её не давит. Остаётся какой-то постоянный уровень, пропорциональный постоянке, хотя обещается, что "давим фпыль". И дажэ плывучка не помогает, т.е. не ошибка от ограниченной разрядности при фикспойнте. Шозах?

mse homjak (25.04.2026 23:24, просмотров: 5419)

Норот, каг-то не парился, а тут, вдруг, заметил, что BPF IIR странно реагирует на постоянку на входе. Он её не давит. Остаётся какой-то постоянный уровень, пропорциональный постоянке, хотя обещается, что "давим фпыль". И дажэ плывучка не помогает, т.е. не ошибка от ограниченной разрядности при фикспойнте. Шозах?

Ответить

- Для полосового и/или ФВЧ сумма коэффициентов в числителе передаточной характеристики цифрового фильтра должна строго равнятся НУЛЮ - argus98(27.04.2026 11:43, +1)
  - Так оно и есть. B0=-B2 B1=0. - mse homjak(27.04.2026 19:53)
    - а можно рекурентное выражение глянуть - IBAH(27.04.2026 20:24)
      - Дык, стандартное выражэние для каноник форм2: Z₀=N_in*A0*Gain+Z_-1*A1+Z_-2*A2. Out=Z₀*B0+Z_-1*B1+Z_-2*B2. - mse homjak(27.04.2026 20:39)
        
        Если я правильно понял эту запись, чтобы проверять на равенство 0, надо сначала раскрыть скобки - IBAH(27.04.2026 20:49)
        
        Не. Ничо не надо раскрывать. Каноник форм2 очень хорошо разделяет вычислительные ветки вокруг линии задержки. Сперва выщитывается Z0. Потом Out. А потом все Z сдвигается по линии задержки для следующей итерации. Qedit выдаёт Gain, A0,A1,A2,B0,B1,B2. Nin, входные данные. mse homjak(1 знак., 27.04.2026 21:00, картинка)
  - Я тоже про эту фишку знаю. А доказательства придумать не могу. - IBAH(1 знак., 27.04.2026 16:56, ссылка)
    - Доказательство элементарно: argus98(297 знак., 27.04.2026 17:12)
      - Спасибо. - Nikolay_Po(27.04.2026 17:13)
- Как раз зависимость от добротности это нормально. Все проблемы растут пропорционально, кажись, квадрату добротности. Чтобы снизить требования к чувствительности коэффициентов надо работать над топологией. А как синтезировал фильтр? Когда-то писал полосовик во флоате, с расчетом коэффициентов налету. Синтезировал методом билинейного отображения. Если IBAH(49 знак., 27.04.2026 10:51)
  - Синтезировал Qedit. - mse homjak(27.04.2026 20:01)
  - Вот IBAH(2053 знак., 27.04.2026 11:31)
    - переменная х1 нигде не используется. - Kpoк(27.04.2026 11:49)
      - а вот dX1, действительно нигде не используется, эта переменная была объявлена глобальной, зачем , не помню, може для тестов - IBAH(27.04.2026 13:30)
      - Все учтено могучим ураганом! x2=x1; x1=x0; - IBAH(27.04.2026 12:26)
- Должен работать и давить постоянку в 0. framer(584 знак., 26.04.2026 20:47, картинка, картинка)
  - ХЗ, постоянку не давит. Несколько фильтров сгенерил, сделал 12 десятичных разрядов после запятой и всё такое. Остаётся. И чем больше добротность, тем больше. Можэт это Qedit глючит, ХЗ. Но как ППФ работает. Частоту туда-сюда гонял, давит. Шум вводил, тожэ давит. - mse homjak(26.04.2026 21:10)
    - Нужна диверсификация инструментов. Проверить коэфициенты. framer(212 знак., 26.04.2026 21:28)
      - Нет там всё нормально. Импульсная харакктеристика падающая. В симуляции есть затухающий переходной процэсс и выход на линию. mse homjak(1 знак., 26.04.2026 21:35, картинка)
        
        Минимальное изменеие коэффициента привело к появлению постоянки на выходе framer(26 знак., 26.04.2026 21:46, картинка)
        
        Походу, точности Qedit не хватает. Хотя у него есть и double. Или сам алгоритм рассчёта багованный. - mse homjak(26.04.2026 21:59)
        
        Относительно диверсификации инструментов . Можно попробовать так. Но мопед не мой :) framer(2339 знак., 26.04.2026 22:19)
        
        Спасибо, попробую. - mse homjak(26.04.2026 22:25)
- Фильтр IIR, по моему, в ноль не обязан выходить. - BlackMorda(1 знак., 26.04.2026 20:32, ссылка)
  - Походу, да. - mse homjak(26.04.2026 21:35)
    - Стремиться должен. - Kpoк(27.04.2026 07:25)
- Повысьте разрядность расчетов. И ноль вернётся к нулю. - Kpoк(26.04.2026 18:47)
  - В Экселе 12 после запятой. Куда ещо точнее? - mse homjak(26.04.2026 20:45)
    - Фильтры в Экселе.... Есть ли пределы бесстыдству? - Kpoк(27.04.2026 07:28)
- Сколько остается конкретно при нулевом значении на входе? У меня было что как раз конечная разрядность оставляла постоянку. - ГopeГope(26.04.2026 16:20, )
  - Ноль, в ноль. даю +/-32767 и на выходе, число, +/-835, при +\-16384, 417 с копейками. при чом, разницы с фикспойнтом и флотом, нет. Да и как BPF он работает, как положено. Но, чем уже полоса(добротность, походу), тем на постоянке число больше. А по переменке, душит внеполосное, как положэно. И частоту сигнала, и шум. - mse homjak(26.04.2026 18:26)
    - Может это не полосовой, а ФНЧ с высокой добротностью полюса? :) - IBAH(27.04.2026 10:51)
      - Нет. Поведение, как у полосового. Давит +/- от центральной. Давит шум. Всё понятно, кроме постоянки. - mse homjak(27.04.2026 20:05)
        
        Тогда АЧХ в студию. Т.е. давит от центральной, а на постоянке подъём? - IBAH(1 знак., 27.04.2026 20:31, картинка)
        
        нет. щас пробежал по частотам, низ давится всего раз в 20, верх, фпыль. Видимо, эллиптический. щас баттерворта посчитаю. - mse homjak(27.04.2026 22:09)
        
        баттерворт помог, но не принципиально. подавил раз в 100 на 1/4 от частоты пропускания фильтра. Видно, это баг Qedit. - mse homjak(27.04.2026 22:27)
        
        а какой может быть батерворт для полосового фильтра 1 порядка??? этот КюЭдит просто кю - IBAH(27.04.2026 22:38)
        
        Кароч, придушил Питона(это не то. что вы подумали). Он выдал коэффициенты, которые позволили при отстройке, примерно, 8% получить подавление под 40дб.Постоянку придушил под 50дб. Правда, полоса у него, получилась раз в 10 уже, чем у Qedit. Ну щас приведём, плус-мнус к одному знаменателю. - mse homjak(27.04.2026 23:39)
        
        нет. степень подавления от добротности не меняется - mse homjak(27.04.2026 23:46)
        
        Я бы вообще код запросил. - Kpoк(27.04.2026 21:40)
        
        чо его запрашивать, он в Экселе. Чисто по формуле, где-то сверху писал. - mse homjak(27.04.2026 22:11)
    - болван говорит, что "это нормально" ибо это проблема FP. POV(2 знак., 26.04.2026 18:30, ссылка, ссылка)
      - Я ужэ промодэлировал в Экселе. Уж там-то плывучка, моё почтение.12 знаков после запятой выставил. Пох. - mse homjak(26.04.2026 19:07)
- Есть такое мнение, что "постоянка" - это голая абстракция. Всегда есть какое то ее изменение, хотя бы при включении. То есть нужно понимать какой уровень подавления, например, на частоте 0,001 Гц. А это совсем не "постоянка". ))) - my504(26.04.2026 12:23)
  - Ну как, абстракция... Цифруем переменку. Подвешиваем её за пол-динамического диапазона и на АЦП. Но если пол-диапазона не совпадает с 0x8000, то имеем смещение. и, соотвецтвенно, если дальше RMS, то у Хьюстона проблема. Особенно, на малом сигнале. 50мВ погрешности квазинуля дали 3 lsb смещения. Повышение добротности даст больше. - mse homjak(26.04.2026 18:40)
    - Простите, сами мы не местные, из колхоза "Стальное вымя". А у вас там унутре не вычитается ли из сигнала эта постоянная составляющая в виде константы? - Бapбoc(26.04.2026 20:11)
      - Ессно, вычитается. Придётся вводить периодическую коррекцию нуля. - mse homjak(26.04.2026 21:05)
        
        Вы и в экселе вычитаете? Лично я ничего не понял. Изначально нужно тестировать фильтр с чистым синусом. Тогда ничего не нужно вычитать. А если требуется удалять постоянку, то постоянку считают отдельным очень НЧ фильтром, а не удаляют константу. - my504(26.04.2026 21:17)
        
        Естественно, первично тестировал в чистом виде, без смещения. Но в жэлезе увидел аномально высокий уровень RMS. Стал разбираться и вышел на влияние постоянки. В Экселе я её добавил и "оппаньки"! Периодическое приведение постоянки в "0" вычитанием из сигнала и будет таким очень НЧ. Можно дажэ назвать цифровой УВХ. - mse homjak(26.04.2026 21:26)
    - Прежде всего на фильтр нужно подать синтетический синус - сначала без смещения, а потом с заданным смещением. Иначе непонятно что вы измеряете. То ли то качество АЦП, то ли внешний антиалиасинг, то ли входной буфер... - my504(26.04.2026 20:05)
      - Я в Экселе смотрю. Увидел в жэлезе. - mse homjak(26.04.2026 20:46)
        
        У вас рекурсивный фильтр. Сколько вы в экселе сделали итераций? Может стоит в железе работать с синусом из флеша? - my504(26.04.2026 20:53)
        
        Моделька занимает 10000 строк. В жэлезе не поэкспериментировать с разрядностью, частотой, шумом. Можно найти минимум-максимум сигналов, чтобы не вылезти из разрядности. - mse homjak(26.04.2026 21:12)
        
        Стесняюсь спросить, модель фильтра 10000 строк? - Kpoк(27.04.2026 08:59)
        
        10000 итераццый. 1 строка - итераццыя. - mse homjak(27.04.2026 10:24)
        
        То есть за 10000 шагов с уровня 30000 опускается до 800, а дальше никак и намёка нет? Ну, значит точно, где-то это 800 превращаются в ноль. - Kpoк(27.04.2026 11:08)
        
        Не просто опускается. Тут картинко. - mse homjak(1 знак., 27.04.2026 19:52, ссылка)

Программируемая логика и ЦОС