Целочисленные алгоритмы ЦОС, надо сформулировать правила Опции

[Zelepuk]

Существует рабочий алгоритм(несколько фильтров и БПФ). Код работает на плавучке на обычном десктопе.
Стоит задача портировать алгоритм на процессор с фиксированой точкой.
Отсчёты, над которыми производятся преобразования, поступают в формате 1.15.
Понимаю что тупое преобразование float в int ни к чему хорошему не приведёт.
Тогда стоит вопрос: какие требования к оперрациям нужно предъявить?

Коэффиценты фильтров и БПФ должны быть тоже в 1.15?
Как производить умножение, деление, суммирование и вычитание, чтобы избегать переполнений?

(лучше если примеры коротенькие на С будут)

Используется 32-битный процессор на ARM9E ядре

Заранее спасибо.



Сообщение отредактировал Zelepuk - Mar 6 2011, 07:48

[Ruslan1]

еще может интересно будет, насчет целочисленной реализации БПФ посмотрите вот это

Вроде нареканий не было и чистый си, хоть куда лепится И очень экономно подходит к используемой памяти, как RAM так и ROM.
Используются 16-битные знаковые числа. Очень удобно напрямую данные с АЦП подсовывать, ну или после нормирования, чтобы максимально использовать все биты.

[MSP430F]

еще может интересно будет, насчет целочисленной реализации БПФ посмотрите вот это

Вроде нареканий не было и чистый си, хоть куда лепится И очень экономно подходит к используемой памяти, как RAM так и ROM.
Используются 16-битные знаковые числа. Очень удобно напрямую данные с АЦП подсовывать, ну или после нормирования, чтобы максимально использовать все биты.

Скажите, а можно этот код переделать на 4096 точек (там сейчас 1024) ? Ответьте, пожалуйста. Очень надо.

[Zelepuk]

Скажите, а можно этот код переделать на 4096 точек (там сейчас 1024) ? Ответьте, пожалуйста. Очень надо.

конечно можно! За отдельную плату готов предоставить готовую реализацию на 4096 точек

P.S. Протестировано на MSP430F47197 в проекте счётчика электроэнергии.

Сообщение отредактировал Zelepuk - Aug 23 2013, 07:57

[MSP430F]

конечно можно! За отдельную плату готов предоставить готовую реализацию на 4096 точек

P.S. Протестировано на MSP430F47197 в проекте счётчика электроэнергии.

Вопрос такой. Какова точность целочисленного 32-разрядного БПФ (вся арифметика не должна выходить за пределы 32 разрядов) ?
Если у меня данные честные 16-битные и я хочу получить на выходе спектр с диапазоном тоже в 90 дБ, то целочисленных вычислений в 32 разрядов для этого не достаточно, не так ли ?

[Zelepuk]

Вопрос такой. Какова точность целочисленного 32-разрядного БПФ (вся арифметика не должна выходить за пределы 32 разрядов) ?
Если у меня данные честные 16-битные и я хочу получить на выходе спектр с диапазоном тоже в 90 дБ, то целочисленных вычислений в 32 разрядов для этого не достаточно, не так ли ?

почему недостатачно? как вы считаете? у меня всё было достатачно точно, единственное что: нужно данные масштабировать в нужных местах.

Сообщение отредактировал Zelepuk - Aug 23 2013, 09:03

[MSP430F]

почему недостатачно? как вы считаете? у меня всё было достатачно точно, единственное что: нужно данные масштабировать в нужных местах.

Вот над этим я сейчас и бьюсь. Вот здесь человек пишет:
" Перед автором стояла задача написать преобразование Фурье на 2048 точек при разрядности исходных данных 16 бит. Из-за отсутствия арифметического сопроцессора пришлось делать целочисленное преобразование, что создало некоторые трудности. При разрядности исходных данных 16 бит разрядность коэффициентов должна быть не менее 16, чтобы не происходило потери точности. Их произведение содержит 32 разряда. 2048 точек дают еще 11 дополнительных разрядов, а это значит, что в 32-разрядное процессорное слово промежуточные данные не помещаются. Вычисление каждой “бабочки” ведется с точностью 64 разряда, а результат округляется до 32 разрядов. "

И как же можно тогда все посчитать в 32 разрядах без потери точности ?

[alex_os]

Вот над этим я сейчас и бьюсь. Вот здесь человек пишет:
" Перед автором стояла задача написать преобразование Фурье на 2048 точек при разрядности исходных данных 16 бит. Из-за отсутствия арифметического сопроцессора пришлось делать целочисленное преобразование, что создало некоторые трудности. При разрядности исходных данных 16 бит разрядность коэффициентов должна быть не менее 16, чтобы не происходило потери точности. Их произведение содержит 32 разряда. 2048 точек дают еще 11 дополнительных разрядов, а это значит, что в 32-разрядное процессорное слово промежуточные данные не помещаются. Вычисление каждой “бабочки” ведется с точностью 64 разряда, а результат округляется до 32 разрядов. "

И как же можно тогда все посчитать в 32 разрядах без потери точности ?

Ужас какой ! При 16 бит входных данных сплошной long long.
Вы определитесь для себя что такое "без потери точности". Можно сделать ФФТ с 16ти битными данными и 16ти битными
поворачивающими множителями, можно сделать с 32 битными данными и 16 ти битными поворачивающими множителями,
кажется АRMы умеют делать умножения 16 * 32 бит c сохранением старших 32 бит.

[DASM]

Ужас какой ! При 16 бит входных данных сплошной long long.
Вы определитесь для себя что такое "без потери точности". Можно сделать ФФТ с 16ти битными данными и 16ти битными
поворачивающими множителями, можно сделать с 32 битными данными и 16 ти битными поворачивающими множителями,
кажется АRMы умеют делать умножения 16 * 32 бит c сохранением старших 32 бит.

UMULL 32x32 = 64 полных во всех с индексом М армах