|
Особенности алгоритма Герцеля, соотношение величнины выборки и частоты дискретизации |
|
|
|
Apr 24 2008, 05:32
|
Группа: Участник
Сообщений: 12
Регистрация: 30-10-07
Пользователь №: 31 878

|
Добрый день! Есть задача - определение в сигнале составлющей 22 Гц. АЦП достаточно мощный. Мне требуется в секунду делать порядка 10 проверок того, что сигнал 22Гц есть. Как оптимизировать частоту дискретизации и величину выборки на обработку для данной задачи?(при величине выборки 50 на высоких частотах дискретизации у меня получается полная фигня) В инете копал, но не нашёл прямого ответа на этот вопрос. Книга Блейхута же так написана, что мне, несмотря на физическое образование, не хватает математической подготовки, чтобы её расшифровать.
|
|
|
|
|
 |
Ответов
|
Apr 24 2008, 06:00
|

Эксперт
    
Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183

|
Цитата(YDM86 @ Apr 24 2008, 09:32)  Добрый день! Есть задача - определение в сигнале составлющей 22 Гц. АЦП достаточно мощный. Мне требуется в секунду делать порядка 10 проверок того, что сигнал 22Гц есть. Как оптимизировать частоту дискретизации и величину выборки на обработку для данной задачи?(при величине выборки 50 на высоких частотах дискретизации у меня получается полная фигня) В инете копал, но не нашёл прямого ответа на этот вопрос. Книга Блейхута же так написана, что мне, несмотря на физическое образование, не хватает математической подготовки, чтобы её расшифровать. Дело не в Герцеле. А в принципе. Герцель - это просто способ вычислять отсчёт ДПФ с помощью простейшего фильтра. Существует принцип непределённости - dF*dT ~ 1. Неопределённость частоты, время измерения... Поэтому, если измерять в течентии секунды - то все частоты находящиеся примерно в окрестности +-1гц от 22 гц будут схвачены. Или другими словами, если частота дискретизации АЦП 8кгц, частота Найквиста 4 кгц то необходимо использовать 4кгц/1гц=4000 точек, чтобы иметь непределённость частоты 1 гц. Где здесь 50 точек? 50 точек дают полосу измерения (неопределённость) 4000/50 = 80 гц. 80 гц много больше 22 гц и прихватывают много мусора, в том числе и огромный дрейф обычно присущий в АЦП на постоянной составляющей Цитата(eugen_pcad_ru @ Apr 24 2008, 09:56)  Алгоритм Герцеля есть по своей сути частный случай БПФ. Только ДПФ. БПФ из другой, быстрой, оперы ))
|
|
|
|
|
Apr 24 2008, 06:21
|
Группа: Участник
Сообщений: 12
Регистрация: 30-10-07
Пользователь №: 31 878

|
Так всё же - мне отправляться от величины выборки для Герцеля и подбирать под неё частоту дискретизации или наоборот? Требуется просто обнаружить 22 Гц и всё. Я пробовал 50 точек в выборке с частотой дискретизации 200 Гц(именно по алгоритму Герцеля). При этом был хороший результат - при сигнале 22,5 Гц мощность(22Гц-овой компоненты), рассчитанная по Герцелю была меньше на порядок, чем если подать именно 22Гц. При 500Гц дискр-ии уже хуже - что 22Гц подавай, что 23, а результат вычисления 22 герцовой компоненты почти один и тот же. Предложите, пожауйста рецепт - я раньше не работал ни сБПФ, ни с ДПФ, так что в предыдущий постах понял не всё, простите меня,теоретика.
|
|
|
|
Сообщений в этой теме
YDM86 Особенности алгоритма Герцеля Apr 24 2008, 05:32 eugen_pcad_ru Цитата(YDM86 @ Apr 24 2008, 08:32) Как оп... Apr 24 2008, 05:56  DRUID3 Цитата(eugen_pcad_ru @ Apr 24 2008, 08:56... Apr 24 2008, 06:42   fontp Цитата(DRUID3 @ Apr 24 2008, 10:42) Это н... Apr 24 2008, 07:22   TSerg Цитата(DRUID3 @ Apr 24 2008, 10:42) Это н... Apr 24 2008, 08:10    DRUID3 Цитата(fontp @ Apr 24 2008, 10:22) Это в ... Apr 24 2008, 10:18     TSerg Цитата(DRUID3 @ Apr 24 2008, 14:18) Все-т... Apr 24 2008, 11:50      DRUID3 Цитата(Stanislav @ Apr 24 2008, 13:31) Ва... Apr 24 2008, 19:02       Stanislav Цитата(DRUID3 @ Apr 24 2008, 23:02) Спаси... Apr 24 2008, 19:18        DRUID3 Цитата(Stanislav @ Apr 24 2008, 22:18) На... Apr 24 2008, 19:53         TSerg Цитата(DRUID3 @ Apr 24 2008, 23:53) А тут... Apr 25 2008, 06:11 YDM86 ЦитатаНемного странно что у Вас с повышением часто... Apr 24 2008, 08:53 fontp Цитата(YDM86 @ Apr 24 2008, 12:53) Т.е. е... Apr 24 2008, 09:08  Stanislav Цитата(fontp @ Apr 24 2008, 13:08) Принци... Apr 24 2008, 10:23 Stanislav Цитата(YDM86 @ Apr 24 2008, 12:53) Т.е. е... Apr 24 2008, 10:04 GetSmart Цитата(Stanislav)Совершенно неверное утверждение. ... Apr 24 2008, 10:16 fontp Цитата(GetSmart @ Apr 24 2008, 14:16) Ген... Apr 24 2008, 10:23  Stanislav Цитата(fontp @ Apr 24 2008, 14:23) Если н... Apr 24 2008, 10:42   fontp Цитата(Stanislav @ Apr 24 2008, 14:42) Ин... Apr 24 2008, 10:49    Stanislav Цитата(fontp @ Apr 24 2008, 14:49) По бло... Apr 24 2008, 10:54   GetSmart Цитата(Stanislav @ Apr 24 2008, 16:42) Во... Apr 24 2008, 10:55    fontp Цитата(GetSmart @ Apr 24 2008, 14:55) Как... Apr 24 2008, 11:27     Stanislav Цитата(DRUID3 @ Apr 24 2008, 23:53) спас... Apr 24 2008, 20:23      fontp Цитата(Stanislav @ Apr 25 2008, 00:23) ОК... Apr 25 2008, 07:28       Stanislav Цитата(fontp @ Apr 25 2008, 11:28) Класси... Apr 25 2008, 12:27        fontp Цитата(Stanislav @ Apr 25 2008, 16:27) Из... Apr 25 2008, 12:35         Stanislav Цитата(fontp @ Apr 25 2008, 16:35) У Марп... Apr 25 2008, 12:43          eugen_pcad_ru Отступление от темы: виноват, алгоритм Герцеля прн... Apr 25 2008, 13:16 Stanislav Цитата(GetSmart @ Apr 24 2008, 14:16) Ген... Apr 24 2008, 10:31 GetSmart Цитата(fontp)Он сказал примерно то же самое, тольк... Apr 24 2008, 10:33 GetSmart Если автору нужно регистрировать появление частоты... Apr 24 2008, 11:35 GetSmart Цитата(TSerg)Он не говорил об уровне ничего.
...
О... Apr 24 2008, 12:01 fontp Цитата(GetSmart @ Apr 24 2008, 16:01) Авт... Apr 24 2008, 12:04 TSerg Цитата(GetSmart @ Apr 24 2008, 16:01) Авт... Apr 24 2008, 12:25 GetSmart Цитата(fontp)Энергию отклика фильтра в полосе прим... Apr 24 2008, 12:11 fontp Цитата(GetSmart @ Apr 24 2008, 16:11) Ну ... Apr 24 2008, 12:13 GetSmart Цитата(TSerg)И ни слова о какой-либо динамике ее и... Apr 24 2008, 12:32 fontp Цитата(GetSmart @ Apr 24 2008, 16:32) Тяж... Apr 24 2008, 12:38 TSerg Цитата(GetSmart @ Apr 24 2008, 16:32) Тяж... Apr 24 2008, 13:00 GetSmart Цитата(fontp)А ширина полосы пофиг, сигнал достато... Apr 24 2008, 13:07 fontp Цитата(GetSmart @ Apr 24 2008, 17:07) При... Apr 24 2008, 13:20 GetSmart Цитата(fontp)Стабильный, значит - идеальный цуг си... Apr 24 2008, 13:29 fontp Цитата(GetSmart @ Apr 24 2008, 17:29) Это... Apr 24 2008, 13:49 GetSmart Цитата(fontp)Вы просто очень упёртый. Вы не нашли ... Apr 24 2008, 15:02 GetSmart Цитата(YDM86)Добрый день! Есть задача - опреде... Apr 24 2008, 16:56 Doka граждане, к порядку!
модер... Apr 24 2008, 18:00 Stanislav Цитата(Doka @ Apr 24 2008, 22:00) граждан... Apr 24 2008, 18:16 GetSmart Цитата(Stanislav)Готов показать, что GetSmart не н... Apr 24 2008, 18:36 GetSmart Цитата(Stanislav)Кроме того, если о полезном сигна... Apr 24 2008, 20:39 YDM86 Не было меня несколько дней - к сожалению. Отвечаю... Apr 29 2008, 03:47 Stanislav Цитата(YDM86 @ Apr 29 2008, 07:47) Есть н... Apr 29 2008, 05:07 TSerg Цитата(YDM86 @ Apr 29 2008, 07:47) Объект... Apr 29 2008, 11:15 ZiB А есть параметры передатчика ?
Было бы проще предл... May 7 2008, 04:40
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|