Хелп. Нужен цифровой фильтр: сишный исходник + расчет коэффициентов. («Телесистемы»: Конференция «Микроконтроллеры и их применение»)

Разработка, производство и продажа радиоэлектронной аппаратуры

| Карта сайта | Пишите нам | В избранное

Требуется программист в Зеленограде
- обработка данных с датчиков; ColdFire; 40 тыс.
e-mail:jobsmp@pochta.ru

Телесистемы | Электроника | Конференция «Микроконтроллеры и их применение»

Хелп. Нужен цифровой фильтр: сишный исходник + расчет коэффициентов.

Отправлено Tow 01 сентября 2008 г. 11:51

Помогите собрать цифровой фильтр НЧ.
Почитал я конференцию и понял что не так быстро это делается как мне хотелось бы :)
Как я понял ничего универсально нет.

Что я имею как исходные данные:
Микроконтроллер с встроенным АЦП (Дельта-Сигма)
Скорость оцифровки - 1.5 кГц
Нужно фильтровать сигнал в полосе частот от 3 до 300 Гц.

Что посоветуете? От чего отталкиваться?
Если можно укажите на простое доходчивое чтиво на эту тему что бы меньше формул, а то куда не глянь, то все ими аж пестрит!

Составить ответ | Вернуться на конференцию

Ответы

Ну вот рекурсивный 2-го порядка, он же биквад - Глупый_пИнгвин (01.09.2008 18:43:29 62.117.81.4, 881 байт)
- Ну ты что имел ввиду? Пишеш BW2_A,BW2_В а где они применены? Что это вообще...{+} - VasilyS (01.09.2008 23:27:9 80.92.96.25, 206 байт)
Какой конкретно МК и какие параметры фильтра нужны (крутизна спадов АЧХ в дБ/окт, линейность/допустимые пульсации АЧХ в полосе пропускания, разрядность данных)? И как понять "Нужно фильтровать сигнал в полосе частот от 3 до 300 Гц" - это перестраиваемый в этих пределах ФНЧ или фиксированная полоса после фильтрации ВЧ и НЧ? - MBedder (01.09.2008 12:01:17 80.92.96.19, пустое)
- Ответ: - Tow (01.09.2008 12:14:6 77.123.60.36, 397 байт)
  - Это значит что у вас полосовой фильтр с частотой среза 3 и 300гц. Да надо е...{+} - MARS (01.09.2008 13:39:16 212.40.41.26, 395 байт)
  - Перестраиваемый, значит надо его коэффициенты менять, а это значит что надо их где то хранить и объем зависит от типа фильтра, его порядка и т.д. - GDI (01.09.2008 13:12:27 195.131.97.15, пустое)
    - Интересно, а зачем менять коэффициенты, если перестраивается только частота среза ФНЧ реального времени? - MBedder (01.09.2008 13:42:16 80.92.96.19, пустое)
      - Не, ну можно, конечно и частотой дискретизации менять частоту среза при постоянных коэффициентах., мы же не знаем что вопрошающему больше подходит. - GDI   (01.09.2008 14:48:2 195.131.97.15, пустое)
        
        Он описал АЦП и МК - все уже прекрасно знаем, поэтому перестройка коэффициентов в данном случае идет лесом - MBedder   (01.09.2008 15:00:5 80.92.96.19, пустое)
        
        Автор не указал с каким дискретом надо полосу регулировать, а я никогда не работал с силабсами и потому не знаю возможностей перестройки частоты дискретизации встроенного АЦП. С внешним АЦП было бы, безусловно легче это сделать. - GDI   (01.09.2008 15:12:40 195.131.97.15, пустое)
        
        Какаянах перестройка частоты дискретизации?? АЦП молотит себе в нужном фиксированном темпе и плюхает данные в кольцевой буфер, а фильтр оттуда выбирает N-ные отсчеты - MBedder   (01.09.2008 15:24:28 80.92.96.19, пустое)
        
        Если автора устроит получать частоты среза кратные основной(300-150-100гц и т.д.), а если ему надо с равномерным шагом в 10Гц, то что тогда делать, децимацию применять, но она тоже с фильтрацией связана и мы рекурсивно приходим опять к изменяемым коэффициентам? Кстати, брать N-ые отсчеты, тоже желательно используя интерполирующие фильтры. - GDI   (01.09.2008 15:57:11 195.131.97.15, пустое)
        
        Читай корень - АЦП на 1500 в/с, верхняя граница фильтра - 100 ГЦ. Даже без интерполяции в лоб можно получить срез на 100, 94, 88, 83 и т.д. ГЦ (чем ниже, тем мельче шаг) - MBedder   (01.09.2008 16:06:3 80.92.96.19, пустое)
        
        А можно по подробнее, может я чего то недопонял? Если мы считаем фильтр на 100Гц, при АЦП на 1500в/с, то беря каждый второй отсчет мы получим 750в/с и тот же фильтр у нас по этим выборкам даст 50Гц, как получить промежуточный результат 94, 88 и т.д? - GDI   (01.09.2008 17:10:28 195.131.97.15, пустое)
        
        На фильтр достаточно десятка-другого выборок, а то и меньше - вот и прикинь, что типа для 100 ГЦ нужна будет каждая 20-я, для 95 - каждая 21-я и т.д. - MBedder   (01.09.2008 19:29:43 80.92.96.19, пустое)
        
        Я понял мысль, да так будет работать, но очень плохо, надо будет отфильтровать все частоты выше половины полученной частоты выборок иначе они будут биения давать, а это еще один ФВЧ аппаратный или программный. Ладно, завершим дискуссию, домой уже пора :) - GDI   (01.09.2008 19:46:47 195.131.97.15, пустое)
      - я менял коффициенты. Скажем перестроить с 4 на 8 гц. Надо частоту с 500 под...{+} - MARS   (01.09.2008 13:46:38 212.40.41.26, 276 байт)
        
        Какаянах еще ОСь? У аффтара есть АЦП и МК, АЦП пусть гонит данные в длинный буфер, а фильтр просто берет в обработку каждую N-ную выборку в зависимости от выбранной частоты среза - все коэффициенты и структура фильтра остаются неизменными - MBedder   (01.09.2008 13:55:38 80.92.96.19, пустое)
        
        Да понятно. я про свою задачу. Менял коффициенты. - MARS   (01.09.2008 14:32:32 212.40.41.26, пустое)
  - Это не высшая математика, а алфавит - позорно браться за решение задачи, не удосужившись разобраться даже в терминологии. Для сведения - простой фильтр с пологой АЧХ делается в несколько строк элементарного кода, а сложный с крутыми срезами - гораздо сложнее. Коэффициенты и структура фильтра при перестройке частоты среза могут оставаться неизменными, а изменяться может либо частота подачи отсчетов на фильтр, либо длина задержки отсчетов - MBedder (01.09.2008 13:08:42 80.92.96.19, пустое)
Nuherz Filter Solution. Делает именно то что Вам надо и без зауми. В осле есть. Могу на фтп Электроникс выложить или на Рапидшару. Архив RAR 7 Мб . - Doppler (01.09.2008 11:59:31 85.223.109.75, пустое, ссылка)
- Выложите куда Вам удобно. /* Только скажите как его оттудова стащить */ - Tow (01.09.2008 12:17:49 77.123.60.36, пустое)
  - Тяни с Рапиды. Просто кликни ссылку. - Doppler (01.09.2008 12:50:3 85.223.109.75, пустое, ссылка)
    - Стащил. - Tow (01.09.2008 13:08:47 77.123.60.36, пустое)
    - Правильная ссылка. - Doppler (01.09.2008 13:03:2 85.223.109.75, пустое, ссылка)
- Какой номер? - Codavr (01.09.2008 12:14:42 213.187.100.43,213.187.127.99, пустое)
  - 11 - Doppler (01.09.2008 12:50:37 85.223.109.75, пустое)
Уточнение. - Tow (01.09.2008 11:57:29 77.123.60.36, 176 байт)
- Нихрена себе, уточнение.... Ну так нагенерируйте коэффициентов через 10Гц. Плавно хочется? Тогда рассчитывайте их сами, я не помощник. - Doppler (01.09.2008 12:01:16 85.223.109.75, пустое)
Матлаб такое делает и коэффициенты и Си-шный код, только код не для всякого МК подойдет, естественно, код только для фильтра. - GDI (01.09.2008 11:57:3 195.131.97.15, пустое)
- Матлабом к сожалению не владею. Говорят есть специализированное ПО (и довол...{+} - Tow (01.09.2008 11:59:7 77.123.60.36, 219 байт)
  - Я себе фвч (чатсота среза 4 гц, дискритизация 500гц) считал в интернете. П...{+} - MARS (01.09.2008 12:12:54 212.40.41.26, 395 байт)
    - Заранее спасибо. - Tow (01.09.2008 12:15:29 77.123.60.36, пустое)
      - просматривай конфу вечером... - MARS   (01.09.2008 12:52:39 212.40.41.26, пустое)
        
        Хорошо. - Tow   (01.09.2008 12:59:43 77.123.60.36, пустое)
        
        Во или найди первоисточник в электронном варианте - MARS   (01.09.2008 13:25:6 212.40.41.26, пустое, ссылка)
        
        Пример расчета циф. фильтра. - MARS   (01.09.2008 13:26:11 212.40.41.26, пустое)
        
        Медленно качает. - Tow   (01.09.2008 13:43:53 77.123.60.36, пустое)
  - А можно и не владеть. Там есть готовые функции для расчёта коеффициентов и неплохой хелп к ним. Просто наберите в командной строке, например, rcosine(1,4,'fir/sqrt',0.125,7) и получите массив коеффициентов. - Сидоргек (01.09.2008 12:04:9 89.209.87.196, пустое)