[covsh 0]

есть прямоугольный импульс который попадаеться на RC цепь с частотный коэффициент передачи k(jw) , нужно посмотреть что будет ны выходе

 

в книжке нашел, что для RC равен k(jw)=1/(1+j*w*tau), (w=омега)

 

FFT рабочее, проверял

так вот после fft имеем массив дейтвительных и мнимых частей коэффициентов.

далее, k(jw) умножаю на (1-j*w*tau) и получаю

k(jw)=1/(1+w*w*tau*tau)+j(-w*tau)/(1+w*w*tau*tau)

т.е. действительныю и мнимую чать k(jw), которые потом перемножаю с коэффициентами Фурье и делаю обратное преобразование.

 

Так вот всесто нужного пилообразного сигнала получаю совсем другое.

 

Где я ошибся? исходник нужен?

[Jools 0]

Как перемножаешь с коэффициентами Фурье? А вообще, лучше, конечно, исходник.

[-=Space=- 0]

Кстати насчет умножений.

Коэффициенты , они же все меньше 1.

Т.е. можно в типе short int поставить точку после 15 разряда

получим числа от -1 до ~1

Но для сишника они так и остануться -32768 +32767

 

Вообщем получается что то вроде этого

 

short int a,b,r;

long int c;

с = a*b;

r = c>>16;

 

a и b - множители , r - произведение

 

так вот , оптимально ли это по произодительности?

 

-----добавлено------

 

Сегодня пропарился с 16 битной целочисленной арифметикой

Получается не очень.

Или переполнения или слишком много шумов в спектре.

Можно конечно попробовать сдвигать точку вправо, на каждом шаге

 

Стоит ли вообще заморачиваться?

Реализуемо ли FFT в целых числах?

[Eugeno 0]

> т.е. действительныю и мнимую чать k(jw), которые потом перемножаю

> с коэффициентами Фурье и делаю обратное преобразование.

 

А перемножение с комплексно-сопряжёнными коэффициентами делал?

 

Сегодня пропарился с 16 битной целочисленной арифметикой

Получается не очень.

Или переполнения или слишком много шумов в спектре.

Можно конечно попробовать сдвигать точку вправо, на каждом шаге

Стоит ли вообще заморачиваться?

Реализуемо ли FFT в целых числах?

 

Конечно реализуемо.

Коэффициенты-то конечно в диапазоне [-1; 1], но по ходу вычисления FFT происходит ещё и суммирование. Поэтому надо сигнал брать с запасом в log2(N) бит. 16 бит для сигнала в этом случае может оказаться мало. Если тема тебе ещё интересна, пиши, расскажу про некоторые ухищрения.

[-=Space=- 0]

2Eugeno:

 

идея "сдвигать точку вправо, на каждом шаге" удалась, алгоритм функционирует.

После некоторых оптимизаций удалось довести скорость FFT1024 до ~300 мс на проце msp430f161 с 3,6864МГц кварцем. Подсчет модуля считается столькоже , но я его почти неоптимизировал. Пока на этом остановился.

Сейчас переключился на разрюхивание TMS320v5509a

[Eugeno 0]

идея "сдвигать точку вправо, на каждом шаге" удалась, алгоритм функционирует.

<{POST_SNAPBACK}>

В зависимости от необходимой точности можно уменьшить битность коэффициентов, что позволит ускорить fft раза в полтора.

 

А для 55-й серии должны быть готовые варианты fft.

Для c6x-серии есть оптимизированные версии от TI и для версии с пошаговым сдвигом, и без онного, и с упрощёнными коэффициентами.