Заранее извиняюсь перед модератором за создание темы по столь мелкому и возможно уже задававшемуся вопросу. Господа Коты, поделитесь пожалуйста кодом БФП на Си без ассамблерных приблуд хотелось бы,если есть конечно В высшей математике не силен,потому и сам написать не в состоянии подобную процедуру. Прошу помощи. Заранее спасибо.
2 ARV: Пока что просто для изучения и попытки постичь данный алгоритм. Просто с ассамблером не очень дружу. Если пример для АВР студии но там все самое интересное вынесено в ASM. Может быть есть ещё какие то алгоритмы для получения частотного спектра сигнала кроме БФП и ДФП ?
алгоритмы-то есть, но для понимания они ничуть не проще, чем FFT или DFT. кстати, DFT - самое простое для понимания. если цель - освоить FFT вообще, то лучше делать это на компе, а не МК: и проще, и инстументарий удобнее, и готовых исходников море.
на мощные МК типа ARM, mega32 или dsPIC сишные тексты от компа переносятся практически без правки, а для простого AVR без ассемблерных вставок быстродействие такое, что FFT превращается в SFT - медленное преобразование Фурье
_________________ если рассматривать человека снизу, покажется, что мозг у него глубоко в жопе при взгляде на многих сверху ничего не меняется...
Качественное и безопасное устройство, работающее от аккумулятора, должно учитывать его физические и химические свойства, профили заряда и разряда, их изменение во времени и под влиянием различных условий, таких как температура и ток нагрузки. Мы расскажем о литий-ионных аккумуляторных батареях EVE и нескольких решениях от различных китайских компаний, рекомендуемых для разработок приложений с использованием этих АКБ. Представленные в статье китайские аналоги помогут заменить продукцию западных брендов с оптимизацией цены без потери качества.
Компания EVE выпустила новый аккумулятор серии PLM, сочетающий в себе высокую безопасность, длительный срок службы, широкий температурный диапазон и высокую токоотдачу даже при отрицательной температуре.
Эти аккумуляторы поддерживают заряд при температуре от -40/-20°С (сниженным значением тока), безопасны (не воспламеняются и не взрываются) при механическом повреждении (протыкание и сдавливание), устойчивы к вибрации. Они могут применяться как для автотранспорта (трекеры, маячки, сигнализация), так и для промышленных устройств мониторинга, IoT-устройств.
Здравствуйте коты разобравшиеся в преобразованиях Фурье!
Сейчас разбираюсь с ДПФ и БПФ... Поправите пожалуйста если не прав... Я так понимаю что мы берем N - количество выборок АЦП( например 1024), с частотой дискретизации например 44100 кГц... Потом после расчета действительной части, получаем N - амплитуд гармоник... Правильно ли я понял, что мы получаем 1024 частот гармонических колебаний, если да то какие это частоты (какие это 1024 частоты)?
Амплитуда гармоники вычисляется по формуле 4, если смотреть здесь. Это корень квадратный из суммы квадратов действительной и мнимой частей. Считая ДПФ по 1024 точкам, вы получаете 1024 отчета в частотной области. Всё верно.
При частоте дискретизации 44100 кГц (наверное имелось ввиду 44,100 кГц) частота гармоники считается f=m*Fдискр/N. N=1024, m изменяется от 0 до 1024. Таким образом, например 1-я гармоника равна 1*44100/1024=0,043 кГц, вторая = 0,086 кГц. Следовательно, шаг получается 43 Гц.
Советую почитать Richard G. Lyons Understanding Digital Signal Processing. В этой книге на пальцах разбирается ДПФ, БПФ, всякие БИХ фильтры и прочие ДСПшные методы.
_________________ Когда уже ничего не помогает - прочтите, наконец, инструкцию. Лучший оптимизатор находится у вас между ушей. (Майкл Абраш, программист Quake и QuakeII) Избыток информации ведёт к оскудению души - Леонтьев А. (сказано в 1965 г.)
Спасибо GRAF, за формулу... Kavka за сайтик тоже спасибо, я правда его видел но не стал рыть так как он на английском, по английски понимаю, но по русски все равно лучше, если сайт стоящий поковыряем... Буду разбираться если что обращусь за помощью
А вообще есть желание замутить анализатор спектра аудио, в качестве практики и тренировки... конечно уже есть куча проектов, но хочеться разобраться и реализовать самому
По ссылке который дал TRIOD Переделав код под IAR AVR, получил 11 мсек на одно преобразование. Мк Атмега8 тактовая 16 МГц. 128 точек. (можно использовать 32,64 или 128 точек) исходник на СИ
Вариант БПФ http://elm-chan.org/works/akilcd/report_e.html портированного в IAR_AVR. с примером в протеусе. Код чуть весомее оригинала. Время одного преобразования для 128 точек (7,18 мсек). МК Атмега 8, такт 16 МГц.
Вариант БПФ http://elm-chan.org/works/akilcd/report_e.html портированного в IAR_AVR. с примером в протеусе. Код чуть весомее оригинала. Время одного преобразования для 128 точек (7,18 мсек). МК Атмега 8, такт 16 МГц.
Вложение:
IAR_AVR_FFT_C_ASSEMB.rar
128 точек какой разрядности, 8 бит? 16? 32? и какой разрядности результаты?
Вариант БПФ http://elm-chan.org/works/akilcd/report_e.html портированного в IAR_AVR. с примером в протеусе. Код чуть весомее оригинала. Время одного преобразования для 128 точек (7,18 мсек). МК Атмега 8, такт 16 МГц.
Продолжается изучение ДПФ и БПФ... написал простое приложение на C# ДПФ на 1024 отсчета, в качестве исследуемого сигнала взял сумму 3 гармоник с разными амплитудами... Алгоритм прост и отлично работает:
Код:
for (int k = 0; k < 1024; k++) { for (int n = 0; n < 1024; n++) { x_real[k] += x_n[n] * Math.Cos( (-2) * Math.PI * k * n / 1024 ); x_imag[k] += x_n[n] * Math.Sin( (-2) * Math.PI * k * n / 1024); }
вопрос следующий: стоит ли бороться с зеркальным эффектом, или просто не учитывать эти частоты так как они находятся за частотой найквиста и все равно не могут быть вычислены верно ???
Вложения:
Комментарий к файлу: зеркальный эффект dft.JPG [72.39 KiB]
Скачиваний: 1347
Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 30
Вы не можете начинать темы Вы не можете отвечать на сообщения Вы не можете редактировать свои сообщения Вы не можете удалять свои сообщения Вы не можете добавлять вложения
Вход
Регистрация
Правила
FAQ
Поиск
