реклама на сайте
подробности

 
 
> Цифровой полосовой фильтр, Реализация
e-tarasov
сообщение Jun 27 2007, 10:21
Сообщение #1


Участник
*

Группа: Свой
Сообщений: 44
Регистрация: 30-05-06
Пользователь №: 17 578



Задача такова. Требуется реализовать цифровой полосовой фильтр на 2000Гц с полосой 10Гц и максимально высоким затуханием и что интересно нужна фаза сигнала (но про это потом). Проц не сигнальный процесс фильтрации необязательно непрерывный. Раньше ЦОС не занимался на считая простеньких ФНЧ 2-го порядка. Надыбал алгоритм БПФ. Алгоритм БПФ позволяет использовать целочисленную математику что хорошо и при этом выделяет фазу сигнала. Господа специалисты а может это то что мне надо. Если да то может быть тогда для одной частоты есть упрощение алгоритма.
Не подскажите заодно сколько отсчетов надо выбрать и какую частоту дискретизации. Плз. Фазу то я по любому думаю вытащу. Но с обычным фильтром траблы с непрерывной фильтрацией или с переходными процессами после подключения. Может кто уже нечто подобное реализовывал. Спасибо.
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
e-tarasov
сообщение Jun 28 2007, 09:18
Сообщение #2


Участник
*

Группа: Свой
Сообщений: 44
Регистрация: 30-05-06
Пользователь №: 17 578



DRUID3 правильно пишет по теореме Котельникова *2 но надо с запасом Крякнутую Filter Solutions не нашел есть же mathlab. А вот по матлабу вопросик может ли он расчитывать и моделировать фильтр БПФ и как.
Go to the top of the page
 
+Quote Post
DRUID3
сообщение Jun 28 2007, 10:46
Сообщение #3


山伏
*****

Группа: Свой
Сообщений: 1 827
Регистрация: 3-08-06
Из: Kyyiv
Пользователь №: 19 294



Цитата(e-tarasov @ Jun 28 2007, 07:29) *
Не получится ли что добротность фильтра уменьшится из-за использования не целочисленной математики ведь ДПФ хорош в теории но на практике вопрос.

biggrin.gif Нет, ДПФ хорош сам по себе и совершенно не имеет никакого отношения к целочисленной арифметике.
Цитата(e-tarasov @ Jun 28 2007, 07:29) *
Что такое окно ханна и оверсемпленг.

Окно - функция, на которую множат интервал накопления, это задает АЧХ каждого частотного отсчета. Если он вообще ни на что не умножен, то окно прямоугольное и АЧХ sinc(f).
Оверсемпленг это когда отсчетов берут в несколько раз (а лучше на несколько порядков) больше чем нужно по Котельникову. Делают это для
1) уменьшения влияния случайных помех.
2)уменьшения апертурой погрешности АЦП.
3)при работе не с прямоугольными окнами для исключения их влияния на сигнал.


Цитата(e-tarasov @ Jun 28 2007, 07:29) *
И вопрос про затухание к примеру насколько сравнимо с обычным цифровым.

Если честно, первый раз такое слышу biggrin.gif . Там нет никакого затухания, т.к. все можно домножить biggrin.gif , но есть шумы (включая дискретизацию или округление) ну и такие понятия как коэффициент передачи, например при умножениях двух сигналов в формате int это важно.


Цитата(e-tarasov @ Jun 28 2007, 12:18) *
DRUID3 правильно пишет по теореме Котельникова *2 но надо с запасом Крякнутую Filter Solutions не нашел есть же mathlab. А вот по матлабу вопросик может ли он расчитывать и моделировать фильтр БПФ и как.

Может, но я, например его не использую принципиально smile.gif . Это связано с моими специфическими убеждениями, но советовать всем поступать так же я просто не в праве. А Filter Solutions насколько я помню фильтров, указывающих фазу, не строит, или я туплю 07.gif ?

Цитата(e-tarasov @ Jun 28 2007, 07:29) *
Алгоритм Герцеля чем он хорош.

Кстати, алгоритм Герцеля в Вашем случае оптимален. Благодарите el34. Я почему-то о нем и не подумал, но ведь у Вас все условия, так сказать. Узкий участок (хоть он работает для любого, но даже на таком он даст фору ДПФ) спектра, и действительный сигнал на входе. Помимо меньших вычислительных затрат Герцель имеет еще одно серьезное преимущество - у него нет шага сетки частот, т.е. настройка производится плавно. Да и времени накопления нет - т.е. время вычисления "наличия" частоты ~ее периоду и просто суперовские свойства для параллельной организации вычислений (например PLIS). Кстати на WiKi есть исходник.


--------------------
Нас помнят пока мы мешаем другим...
//--------------------------------------------------------
Хороший блатной - мертвый...
//--------------------------------------------------------
Нет старик, это те дроиды которых я ищу...
Go to the top of the page
 
+Quote Post
-=ВН=-
сообщение Jun 28 2007, 11:51
Сообщение #4


Местный
***

Группа: Новичок
Сообщений: 210
Регистрация: 3-11-06
Пользователь №: 21 936



Цитата(DRUID3 @ Jun 28 2007, 14:46) *
biggrin.gif Нет, ДПФ хорош сам по себе и совершенно не имеет никакого отношения к целочисленной арифметике.

Окно - функция, на которую множат интервал накопления, это задает АЧХ каждого частотного отсчета. Если он вообще ни на что не умножен, то окно прямоугольное и АЧХ sinc(f).
Оверсемпленг это когда отсчетов берут в несколько раз (а лучше на несколько порядков) больше чем нужно по Котельникову. Делают это для
1) уменьшения влияния случайных помех.
2)уменьшения апертурой погрешности АЦП.
3)при работе не с прямоугольными окнами для исключения их влияния на сигнал.
Если честно, первый раз такое слышу biggrin.gif . Там нет никакого затухания, т.к. все можно домножить biggrin.gif , но есть шумы (включая дискретизацию или округление) ну и такие понятия как коэффициент передачи, например при умножениях двух сигналов в формате int это важно.
Может, но я, например его не использую принципиально smile.gif . Это связано с моими специфическими убеждениями, но советовать всем поступать так же я просто не в праве. А Filter Solutions насколько я помню фильтров, указывающих фазу, не строит, или я туплю 07.gif ?
Кстати, алгоритм Герцеля в Вашем случае оптимален. Благодарите el34. Я почему-то о нем и не подумал, но ведь у Вас все условия, так сказать. Узкий участок (хоть он работает для любого, но даже на таком он даст фору ДПФ) спектра, и действительный сигнал на входе. Помимо меньших вычислительных затрат Герцель имеет еще одно серьезное преимущество - у него нет шага сетки частот, т.е. настройка производится плавно. Да и времени накопления нет - т.е. время вычисления "наличия" частоты ~ее периоду и просто суперовские свойства для параллельной организации вычислений (например PLIS). Кстати на WiKi есть исходник.

Жуткая смесь smile.gif
"Оверсампленг" - это не увеличение числа отсчетов, это повышение ч-ты дискретизации.
Алгоритм Герцеля - алгоритм вычисления ДПФ, со всеми вытекающими. Преимущества - экономия двойная по сравнению с лобовым вычислением. Других преимуществ нет.
Вообще вопрос был о фильтрации. Если автору нужен отфильтрованный сигнал, то ДПФ и даже с окном его не спасут. smile.gif
Go to the top of the page
 
+Quote Post
DRUID3
сообщение Jun 28 2007, 12:10
Сообщение #5


山伏
*****

Группа: Свой
Сообщений: 1 827
Регистрация: 3-08-06
Из: Kyyiv
Пользователь №: 19 294



Цитата(-=ВН=- @ Jun 28 2007, 14:51) *
Жуткая смесь smile.gif
"Оверсампленг" - это не увеличение числа отсчетов, это повышение ч-ты дискретизации.

Да, именно это я и хотел сказать... Прошу прощения у тех, кто не может или ленится выхватывать из контекста... blush.gif

Цитата(-=ВН=- @ Jun 28 2007, 14:51) *
Алгоритм Герцеля - алгоритм вычисления ДПФ, со всеми вытекающими. Преимущества - экономия двойная по сравнению с лобовым вычислением. Других преимуществ нет.

Что Вы говорите... biggrin.gif Хотя, может в Вашей реализации это так и есть ... lol.gif


--------------------
Нас помнят пока мы мешаем другим...
//--------------------------------------------------------
Хороший блатной - мертвый...
//--------------------------------------------------------
Нет старик, это те дроиды которых я ищу...
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- 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
|- - -=ВН=-   Цитата(DRUID3 @ Jun 28 2007, 16:10) Да, и...   Jun 28 2007, 12:30
|- - DRUID3   Цитата(-=ВН=- @ Jun 28 2007, 15...   Jun 28 2007, 12:39
|- - -=ВН=-   Цитата(DRUID3 @ Jun 28 2007, 16:39) Оу-у-...   Jun 28 2007, 12:47
|- - DRUID3   Цитата(-=ВН=- @ Jun 28 2007, 15...   Aug 14 2007, 12:45
|- - -=ВН=-   Цитата(DRUID3 @ Aug 14 2007, 16:45) Про...   Aug 14 2007, 15:36
||- - DRUID3   Цитата(-=ВН=- @ Aug 14 2007, 18...   Aug 15 2007, 10:20
||- - blackfin   Цитата(DRUID3 @ Aug 15 2007, 14:20) Ну чт...   Aug 15 2007, 10:56
||- - -=ВН=-   Цитата(DRUID3 @ Aug 15 2007, 14:20) Даже...   Aug 15 2007, 12:26
|- - 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


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 20th July 2025 - 21:05
Рейтинг@Mail.ru


Страница сгенерированна за 0.0146 секунд с 7
ELECTRONIX ©2004-2016