Встала задача реализовать алгоритм дискретного косинусного преобразования на контроллере на базе процессора PowerPC. Решающую роль играет производительность. В связи с этим возникло несколько вопросов.

1. Известны алгоритмы быстрого преобразования фурье, работающие с выборками размера 2^n. Существуют ли подобные алгоритмы для DCT?

2. Существуют ли целочисленные алгоритмы? Плавающая арифметика сильно бы замедлила вычисления.

3. Существуют ли инкрементальные алгоритмы? Через равные промежутки времени на контроллер поступают значения извне; и каждый раз нужно пересчитывать DCT. При этом самое старое значение выбрасывается из выборки, а с другой стороны к ней добавляется новое. Есть ли возможность как-то сократить вычисления в этом случае?

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

Буду благодарен тем, кто поделится опытом, приведёт ссылки на конкретные алгоритмы

http://www.hpl.hp.com/techreports/96/HPL-96-03.html
еще на ftp.analog.com есть примеры реализации быстрого 2d dct 8х8, но на асме для blackfin'a.
3. imho каждый раз придется все пересчитывать.
4. отрицательных частот после dct вроде не должно получаться.