Цифровой фильтр Опции

[ASV]

На сколько мне известно порядок фильтра характеризует коэффициент затухания. В цифровых фильтрах, как я убедился он то же влияет на коэффициент затухания. А зависит ли Порядок фильтра от количества выборок? Могу ли я поступить следующим образом:
- есть сигнал 100Гц;
- я оцифровую его с частотой 1000Гц;
- получаеться десять выборок;
- фильтр 8-го порядка;
Вопрос: имею ли я право пропустить через фильтр все 10 отсчетов? И что получиться. Фильтр НЧ с КИХ.
Может это все глупо, но я столкнулся с цифровыми фильтрами первый раз.

[BVU]

Для выяснения причин данной проблемы нужно подать на вход системы (АЦП) синус 316гц, и контролировать, что получается на выходе фильтра 316гц и 430гц. Так же синус и для 430гц. Если будет 'пролаз' значит фильтры не работают, как полосовые или имеют не расчетную полосу, а на много больше. Так же характеристику фильтра можно снять при помощи анализатора спектра (Брюль и Кьер например или другой...), что еще комфортней для анализа и настройки фильра.

[skopus]

Для выяснения причин данной проблемы нужно подать на вход системы (АЦП) синус 316гц, и контролировать, что получается на выходе фильтра 316гц и 430гц. Так же синус и для 430гц. Если будет 'пролаз' значит фильтры не работают, как полосовые или имеют не расчетную полосу, а на много больше. Так же характеристику фильтра можно снять при помощи анализатора спектра (Брюль и Кьер например или другой...), что еще комфортней для анализа и настройки фильра.

http://webfile.ru/766585 в этом файле картинки с показаниями со входов и выходов фильтров.

Да фильтр 430гц пропускает и 316 гц. Да он не работает. Вопрос ПОЧЕМУ??? Рачетные коэффициенты верны, программная реализация верна. Почему же он не работает. Может дело в математике?

http://webfile.ru/769075 - это принципиальная схема. Там же на картинке 3 решаемых уравнения.
Но на самом деле там есть спорный момент. По уравнениям я не вижу на втором сумматоре прибавления сигнала задержанного на 2 отсчета (который называется D2) . А если этого нет тогда то, что у меня называется "Рекурсивной частью звена", на самом деле - цифровой резонатор...

Посмотрите пожалуйста!

Если я прав, то этим можно объяснить то, что амплитуда сигнала частотой 316гц на выходе фильтра, расчитанного на эту частоту, в несколько раз превышает входной сигнал... Но почему тогда фильтр на 430гц одинаково резонирует как на 430 так и на 316?

[BVU]

Для выяснения причин данной проблемы нужно подать на вход системы (АЦП) синус 316гц, и контролировать, что получается на выходе фильтра 316гц и 430гц. Так же синус и для 430гц. Если будет 'пролаз' значит фильтры не работают, как полосовые или имеют не расчетную полосу, а на много больше. Так же характеристику фильтра можно снять при помощи анализатора спектра (Брюль и Кьер например или другой...), что еще комфортней для анализа и настройки фильра.

http://webfile.ru/766585 в этом файле картинки с показаниями со входов и выходов фильтров.

Да фильтр 430гц пропускает и 316 гц. Да он не работает. Вопрос ПОЧЕМУ??? Рачетные коэффициенты верны, программная реализация верна. Почему же он не работает. Может дело в математике?

http://webfile.ru/769075 - это принципиальная схема. Там же на картинке 3 решаемых уравнения.
Но на самом деле там есть спорный момент. По уравнениям я не вижу на втором сумматоре прибавления сигнала задержанного на 2 отсчета (который называется D2) . А если этого нет тогда то, что у меня называется "Рекурсивной частью звена", на самом деле - цифровой резонатор...

Посмотрите пожалуйста!

Если я прав, то этим можно объяснить то, что амплитуда сигнала частотой 316гц на выходе фильтра, расчитанного на эту частоту, в несколько раз превышает входной сигнал... Но почему тогда фильтр на 430гц одинаково резонирует как на 430 так и на 316?

Судя по картинкам (http://webfile.ru/766585) у фильтра происходит возбуждение - отсюда и увеличение амплитуды выходного сигнала (это видно на эффекте модуляции), по всей видимости из за неустойчивости рекурсивной части звена. Что конкретно и почему сказать немогу. Попробуйте взять за основу уже работающую на практике модель фильтра (смена математической концепции)...

[Stanislav]

...Судя по картинкам (http://webfile.ru/766585) у фильтра происходит возбуждение - отсюда и увеличение амплитуды выходного сигнала (это видно на эффекте модуляции), по всей видимости из за неустойчивости рекурсивной части звена. Что конкретно и почему сказать немогу. Попробуйте взять за основу уже работающую на практике модель фильтра (смена математической концепции)...

Да, я тоже считаю, что происходит самовозбуждение из-за недостаточной точности вычислений, хотя, чтобы утверждать точно, нужно моделировать, а на это нет времени. Для уменьшения эффекта нужно отодвинуть полюса фильтра внутрь единичной окружности (в Z-плоскости). Для этого достаточно модифицировать к-ты фильтра: A1 умножить на число 0,9<b<0.98, а А2 - на b^2. Сначала выберите b близким к 0,9. Если устойчивость не будет достигнута - причина ошибки в другом. Возможно, эффекты округления промежуточных данных и коэффициентов фильтра. Так что рекомендую прислушаться к совету - использовать матлаб.

ЗЫ. Из-за постоянки такая беда тоже может быть. Попробуйте перед укладыванием числа в буфер, вычесть из него постоянку.

Сообщение отредактировал Stanislav - Jan 25 2006, 16:13