| |
 Цифровой полосовой фильтр, Реализация |
|
|
|
|
 Jun 27 2007, 10:21
|
Участник
Группа: Свой
Сообщений: 44
Регистрация: 30-05-06
Пользователь №: 17 578
|
Задача такова. Требуется реализовать цифровой полосовой фильтр на 2000Гц с полосой 10Гц и максимально высоким затуханием и что интересно нужна фаза сигнала (но про это потом). Проц не сигнальный процесс фильтрации необязательно непрерывный. Раньше ЦОС не занимался на считая простеньких ФНЧ 2-го порядка. Надыбал алгоритм БПФ. Алгоритм БПФ позволяет использовать целочисленную математику что хорошо и при этом выделяет фазу сигнала. Господа специалисты а может это то что мне надо. Если да то может быть тогда для одной частоты есть упрощение алгоритма.
Не подскажите заодно сколько отсчетов надо выбрать и какую частоту дискретизации. Плз. Фазу то я по любому думаю вытащу. Но с обычным фильтром траблы с непрерывной фильтрацией или с переходными процессами после подключения. Может кто уже нечто подобное реализовывал. Спасибо.
|
|
|
|
|
|
|
 |
Ответов
|
|
Jun 28 2007, 09:18
|
Участник
Группа: Свой
Сообщений: 44
Регистрация: 30-05-06
Пользователь №: 17 578
|
DRUID3 правильно пишет по теореме Котельникова *2 но надо с запасом Крякнутую Filter Solutions не нашел есть же mathlab. А вот по матлабу вопросик может ли он расчитывать и моделировать фильтр БПФ и как.
|
|
|
|
|
|
|
|
Jun 28 2007, 10:46
|
山伏
Группа: Свой
Сообщений: 1 827
Регистрация: 3-08-06
Из: Kyyiv
Пользователь №: 19 294
|
Цитата(e-tarasov @ Jun 28 2007, 07:29)  Не получится ли что добротность фильтра уменьшится из-за использования не целочисленной математики ведь ДПФ хорош в теории но на практике вопрос.  Нет, ДПФ хорош сам по себе и совершенно не имеет никакого отношения к целочисленной арифметике. Цитата(e-tarasov @ Jun 28 2007, 07:29) Что такое окно ханна и оверсемпленг. Окно - функция, на которую множат интервал накопления, это задает АЧХ каждого частотного отсчета. Если он вообще ни на что не умножен, то окно прямоугольное и АЧХ sinc(f). Оверсемпленг это когда отсчетов берут в несколько раз (а лучше на несколько порядков) больше чем нужно по Котельникову. Делают это для 1) уменьшения влияния случайных помех. 2)уменьшения апертурой погрешности АЦП. 3)при работе не с прямоугольными окнами для исключения их влияния на сигнал. Цитата(e-tarasov @ Jun 28 2007, 07:29) И вопрос про затухание к примеру насколько сравнимо с обычным цифровым. Если честно, первый раз такое слышу . Там нет никакого затухания, т.к. все можно домножить , но есть шумы (включая дискретизацию или округление) ну и такие понятия как коэффициент передачи, например при умножениях двух сигналов в формате int это важно. Цитата(e-tarasov @ Jun 28 2007, 12:18) DRUID3 правильно пишет по теореме Котельникова *2 но надо с запасом Крякнутую Filter Solutions не нашел есть же mathlab. А вот по матлабу вопросик может ли он расчитывать и моделировать фильтр БПФ и как. Может, но я, например его не использую принципиально  . Это связано с моими специфическими убеждениями, но советовать всем поступать так же я просто не в праве. А Filter Solutions насколько я помню фильтров, указывающих фазу, не строит, или я туплю  ? Цитата(e-tarasov @ Jun 28 2007, 07:29) Алгоритм Герцеля чем он хорош. Кстати, алгоритм Герцеля в Вашем случае оптимален. Благодарите el34. Я почему-то о нем и не подумал, но ведь у Вас все условия, так сказать. Узкий участок (хоть он работает для любого, но даже на таком он даст фору ДПФ) спектра, и действительный сигнал на входе. Помимо меньших вычислительных затрат Герцель имеет еще одно серьезное преимущество - у него нет шага сетки частот, т.е. настройка производится плавно. Да и времени накопления нет - т.е. время вычисления "наличия" частоты ~ее периоду и просто суперовские свойства для параллельной организации вычислений (например PLIS). Кстати на WiKi есть исходник.
--------------------
Нас помнят пока мы мешаем другим...
//--------------------------------------------------------
Хороший блатной - мертвый...
//--------------------------------------------------------
Нет старик, это те дроиды которых я ищу...
|
|
|
|
|
|
|
|
Jun 28 2007, 11:51
|
Местный
Группа: Новичок
Сообщений: 210
Регистрация: 3-11-06
Пользователь №: 21 936
|
Цитата(DRUID3 @ Jun 28 2007, 14:46) Нет, ДПФ хорош сам по себе и совершенно не имеет никакого отношения к целочисленной арифметике. Окно - функция, на которую множат интервал накопления, это задает АЧХ каждого частотного отсчета. Если он вообще ни на что не умножен, то окно прямоугольное и АЧХ sinc(f). Оверсемпленг это когда отсчетов берут в несколько раз (а лучше на несколько порядков) больше чем нужно по Котельникову. Делают это для 1) уменьшения влияния случайных помех. 2)уменьшения апертурой погрешности АЦП. 3)при работе не с прямоугольными окнами для исключения их влияния на сигнал. Если честно, первый раз такое слышу . Там нет никакого затухания, т.к. все можно домножить , но есть шумы (включая дискретизацию или округление) ну и такие понятия как коэффициент передачи, например при умножениях двух сигналов в формате int это важно. Может, но я, например его не использую принципиально . Это связано с моими специфическими убеждениями, но советовать всем поступать так же я просто не в праве. А Filter Solutions насколько я помню фильтров, указывающих фазу, не строит, или я туплю ? Кстати, алгоритм Герцеля в Вашем случае оптимален. Благодарите el34. Я почему-то о нем и не подумал, но ведь у Вас все условия, так сказать. Узкий участок (хоть он работает для любого, но даже на таком он даст фору ДПФ) спектра, и действительный сигнал на входе. Помимо меньших вычислительных затрат Герцель имеет еще одно серьезное преимущество - у него нет шага сетки частот, т.е. настройка производится плавно. Да и времени накопления нет - т.е. время вычисления "наличия" частоты ~ее периоду и просто суперовские свойства для параллельной организации вычислений (например PLIS). Кстати на WiKi есть исходник. Жуткая смесь "Оверсампленг" - это не увеличение числа отсчетов, это повышение ч-ты дискретизации. Алгоритм Герцеля - алгоритм вычисления ДПФ, со всеми вытекающими. Преимущества - экономия двойная по сравнению с лобовым вычислением. Других преимуществ нет. Вообще вопрос был о фильтрации. Если автору нужен отфильтрованный сигнал, то ДПФ и даже с окном его не спасут.
|
|
|
|
|
|
|
|
 Jun 28 2007, 12:10
|
山伏
Группа: Свой
Сообщений: 1 827
Регистрация: 3-08-06
Из: Kyyiv
Пользователь №: 19 294
|
Цитата(-=ВН=- @ Jun 28 2007, 14:51) Жуткая смесь "Оверсампленг" - это не увеличение числа отсчетов, это повышение ч-ты дискретизации. Да, именно это я и хотел сказать... Прошу прощения у тех, кто не может или ленится выхватывать из контекста...  Цитата(-=ВН=- @ Jun 28 2007, 14:51) Алгоритм Герцеля - алгоритм вычисления ДПФ, со всеми вытекающими. Преимущества - экономия двойная по сравнению с лобовым вычислением. Других преимуществ нет. Что Вы говорите... Хотя, может в Вашей реализации это так и есть ... 
--------------------
Нас помнят пока мы мешаем другим...
//--------------------------------------------------------
Хороший блатной - мертвый...
//--------------------------------------------------------
Нет старик, это те дроиды которых я ищу...
|
|
|
|
|
|
|
|
Jun 28 2007, 12:30
|
Местный
Группа: Новичок
Сообщений: 210
Регистрация: 3-11-06
Пользователь №: 21 936
|
Цитата(DRUID3 @ Jun 28 2007, 16:10) Да, именно это я и хотел сказать... Прошу прощения у тех, кто не может или ленится выхватывать из контекста... Ну если Вы именно это и хотели сказать , то необходимость оверсампленга при применении окна вызывает большие сомнения. : Цитата(DRUID3 @ Jun 28 2007, 16:10) Что Вы говорите... Хотя, может в Вашей реализации это так и есть ... Да это у кого угодно так
|
|
|
|
|
|
|
|
Jun 28 2007, 12:47
|
Местный
Группа: Новичок
Сообщений: 210
Регистрация: 3-11-06
Пользователь №: 21 936
|
Цитата(DRUID3 @ Jun 28 2007, 16:39) Оу-у-у... Да у нас тут теоретики завелись... Ну ничего, еще пару лет посидите над практикой и поймете зачем при применении окна - оверсемплинг... Вобщем ответ такой же, как и на первую реплику... P.S.: вот за это я и ненавижу Матлаб... Ясно. Сказать, зачем нужен оверсамленг при применении окна, и как он уничтожает влияния окна на сигнал, Вы не можете. Аналогично и с мифическими преимуществми (кроме вычислительных) алгоритма Герцеля. Этого собственно следовало ожидать. Ваша ссылка на 2-х годичную практику выжала у меня слезу. Скупую мужскую слезу. Я плачу, мой друг. Я плачу.
|
|
|
|
|
|
|
|
Aug 15 2007, 12:26
|
Местный
Группа: Новичок
Сообщений: 210
Регистрация: 3-11-06
Пользователь №: 21 936
|
Цитата(DRUID3 @ Aug 15 2007, 14:20) Даже не пытаюсь, Ну Вам что, цитаты приводить Ваши, где Вы говорили о каких-то преимуществах, причем не вычислительных, Герцеля по сравнению с ДПФ? Теперь у Вас не преимущества, а особенности и не ДПФ, а БПФ. Вы уж о чем-нибудь одном валяйте. А про "оверсамленг" Вы что-то совершенно не то говорите. Вообще. Он никак не связан с окнами. И окна с ним. "Оверсамленг" - увеличение частоты дискретизации. Все. Ничего более. И спектральное разрешение и характеристики окон определяются длительностью анализируемого куска сигнала (окна) и никак с часттотой дискретизации не связаны. Что Вы хотели сказать Вашим примером - непонятно абсолютно. Вам наверное тоже.
|
|
|
|
|
|
Сообщений в этой теме
 e-tarasov Цифровой полосовой фильтр Jun 27 2007, 10:21 DRUID3 ...частота 2000 берем 4096 частоту дискретизации (... Jun 27 2007, 11:47  Самурай Цитата(DRUID3 @ Jun 27 2007, 15:47) ...ча... Jun 27 2007, 12:11 DRUID3 Цитата(Самурай @ Jun 27 2007, 15:11) Не н... Jun 27 2007, 12:56 anton Цитата(Самурай @ Jun 27 2007, 16:11) Не н... Jun 27 2007, 14:43 DRUID3 Цитата(anton @ Jun 27 2007, 17:43) Фильтр... Jun 27 2007, 15:25 el34 DRUID3>Кстати забыл сказать, что ДПФ для такого... Jun 27 2007, 13:30 el34 DRUID3>А ето вообще нИасилил. Сместить 2 kHz ко... Jun 27 2007, 18:24 e-tarasov Спасибо за ответы. Насчет ДПФ. Не получится ли что... Jun 28 2007, 04:29 Alex11 И еще интересный вопрос к автору: "цифровой п... Jun 28 2007, 06:27 blackfin Цитата(DRUID3 @ Aug 15 2007, 14:20) Ну чт... Aug 15 2007, 10:56 blackfin Цитата(DRUID3 @ Aug 14 2007, 16:45) Чтобы... Aug 14 2007, 17:29 BORIV Цитата(blackfin @ Aug 14 2007, 21:29) При... Sep 11 2009, 07:12 alex_os Цитата(e-tarasov @ Jun 27 2007, 14:2... Jul 1 2007, 07:05 el34 alex_os>Сделайте FIR с децимацией. Т.е примерно... Jul 1 2007, 08:41 alex_os Цитата(el34 @ Jul 1 2007, 12:41) alex_os... Jul 2 2007, 06:10 e-tarasov Не знаю пока что такое алгоритм Герцеля и в чем ег... Jul 2 2007, 10:26 el34 alex_os>Почему Вы думаете , что многоступенчата... Jul 2 2007, 10:53 fontp Цитата(el34 @ Jul 2 2007, 14:53) e-taraso... Jul 2 2007, 11:58 alex_os Цитата(el34 @ Jul 2 2007, 14:53) потому, ... Jul 3 2007, 05:02 blackswan мне необходимо было сделать полосовой БИХ-фильтр д... Jul 7 2007, 09:09 BratherLU Загляните в хелп матлаба - Help -> Full Product... Jul 7 2007, 11:14 blackswan а других источников нет?... в моей версии хелп как... Jul 8 2007, 17:18 BratherLU http://www.busim.ee.boun.edu.tr/~resources/fdq.pdf... Jul 9 2007, 05:48 blackswan я там нашел только в виде блок-схем, но в них не о... Jul 20 2007, 14:25 fontp Да, лучше не спорьте, тем более с =ВН=. Тем более ... Aug 15 2007, 09:07 Lalka А не мог бы кто-нибудь выложить алгоритм как можно... Sep 10 2009, 06:41 bahurin Мое мнение на данный вопрос. Сразу предупреждаю в ... Sep 10 2009, 12:53 sup-sup Цитата(Lalka @ Sep 10 2009, 09:41) А не м... Sep 10 2009, 19:53 bahurin Цитата(sup-sup @ Sep 10 2009, 23:53)... Sep 11 2009, 06:58 pokos Непонятно, насколько прямоугольным должен быть фил... Sep 15 2009, 11:25 анатолий Интересный полосовой фильтр получается из фазового... Sep 27 2009, 12:50
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|
