Алгоритм Герцеля

[ivan219]

Интересует возможность динамического расчета, но не на каждом такте дискретизации, а блоками по 6000 семплов.
Частота дискретизации 48кГц размер блок Герцеля 48000 семплов. И хотелось бы 8 раз в секунду обновлять результат добовляя по 6000 семплов. Вычислять 8 раз в секунду полностью весь алгоритм с новыми данными не хочется.

[Бахурин Сергей]

ну как вариант рассчитывать на каждом такте а выводить каждый 6000 -й такт. Второй вариант поупражняться в математике и попробовать для вашего случая вывести алгоритм чтобы рассчитывал нужными вам блоками. Мне помнится были вопросы о несоответствии приведенной в статье схемы, и схемы показанной в литературе. Хотелось узнать вопросы сняты или они еще актуальны?

[ivan219]

Идея выводить каждый 6000 мне то же в голову приходила, но думал, есть более шустрый алгоритм.
О не соответствии вопрос снят в двух книгах одинаковый ошибочный алгоритм он хорошо считает амплитуду, а вот с фазой накладки. Ваш алгоритм работает идеально.

Ещё при динамическом расчёте у вас на вход подаётся комплексный сигнал, как тут быть, если у меня сигнал вещественный?

[Бахурин Сергей]

подавать мнимую часть равную 0. По поводу книг думаю в статье стоит сделать замечание о несоответствии приведенных в них структурных схем алгоритма Герцеля.

[ivan219]

Сейчас посмотрел статьи в Интернете и другие книги такое ощущение что всё содрано с одного источника
Ваш алгоритм лучше. Интересно он был получен из за случайной ошибки или вы его придумали? Если вы то можете назвать его в свою честь Алгоритм Герцель - Бахурин

http://sonder.ru/content/view/473/32/

[Бахурин Сергей]

я просто аккуратно проделал вывод и получил результат отличный от того что было нарисовано в книжке. Проверил экспериментально алгоритм в книжке явно налажал мой вывод дал верный результат. Получив моральное удовлетворение написал статью на сайт.

[ivan219]

Всё отличие заключается в дополнительном коэффициенте, который присутствует в вашем алгоритме http://calculator2006.narod.ru/goertzelmod1.gif
Мне кажется, что этот факт нужно отразить в вашей статье, так как в сети очень много информации по алгоритму Герцеля. И вашь могут счесть ошибочным.

Интересная статейка http://www.dsprelated.com/showarticle/102.php только я так и не понял что тут подрузумевалось под не стабильностью.

[ivan219]

Немного оптимизировал алгоритм динамического перерасчёта для вещественного сигнала.

goertzelmod_html_m1746c635.gif (31.45 КБ) 18172 просмотра

На одно умножение меньше.

goertzelmod_html_m1746c6351.gif (8.73 КБ) 18172 просмотра

[ivan219]

Ещё одна оптимизация алгоритма.
Всего два умножения.

[IgorV]

Немного оптимизировал алгоритм динамического перерасчёта для вещественного сигнала.



goertzelmod_html_m1746c6351.gif

Тут все верно? У меня почему-то ощущение, что I и Q местами поменяны