БПФ по основанию 4, Нужен алгоритм. Опции

[ivan219]

Уковонибудь есть алгоритм вычесления БПФ по основанию 4.

[ivan219]

Спасибо конечно но я в С++ не силён может на Пскале есть???

Сообщение отредактировал ivan219 - Dec 12 2008, 19:59

[vadkudr]

http://www.fftw.org/newsplit.pdf
Вот здесь есть алгоритм еще быстрее, чем сплит-радикс 4. Проверял на матлабе - работает.
И на сам сайт интересно посмотреть.

[diwil]

http://www.fftw.org/newsplit.pdf
Вот здесь есть алгоритм еще быстрее, чем сплит-радикс 4. Проверял на матлабе - работает.
И на сам сайт интересно посмотреть.

это сплит.... Теоретически он работает несколько быстреев смысле требуемого числа операций. Однако при его реализации оказывается, что помимо выполнения основных операций сложения-умножения, нужно еще крутить большое количество счетчиков. Это можно хорошо делать на жесткой логике (типа ПЛМ). Однако при реализации на процессоре (-рах) большое количество счетчиков не позволяет разместить все промежуточные данные в регистрах и приходится задействовать или стэк или память, что может привести к заметному замедлению вычисления.

Не могу говорить про все ДСП... да и вообще про ДСП
но, скажем, на АРМе радикс4 будет самым быстрым.

Хотя, на код матлаба из статьи было бы посмотреть оч интересно.

[vadkudr]

Хотя, на код матлаба из статьи было бы посмотреть оч интересно.

Посмотрите - писался просто для проверки теории в статье. Поэтому заведомо работает медленно.

Запуск:

>> a=randn(1,32)
a =
Columns 1 through 8
-0.4326 -1.6656 0.1253 0.2877 -1.1465 1.1909 1.1892 -0.0376
Columns 9 through 16
0.3273 0.1746 -0.1867 0.7258 -0.5883 2.1832 -0.1364 0.1139
Columns 17 through 24
1.0668 0.0593 -0.0956 -0.8323 0.2944 -1.3362 0.7143 1.6236
Columns 25 through 32
-0.6918 0.8580 1.2540 -1.5937 -1.4410 0.5711 -0.3999 0.6900
>> b=Alg_02_NewSplitRadixFFT(a)
b =
Columns 1 through 4
2.8652 -4.0386 - 2.6955i -3.3874 + 0.9221i -2.1831 - 3.6784i
Columns 5 through 8
3.5893 + 5.2095i 0.9802 + 5.1840i -0.8291 + 0.9154i -2.6340 + 2.6884i
Columns 9 through 12
-5.0758 - 1.0582i -4.9958 + 0.0682i 7.7442 + 0.5439i 2.8367 + 9.6884i
Columns 13 through 16
2.7128 + 4.6691i -1.1440 - 4.9140i 0.4670 + 5.2595i -0.8161 - 2.9031i
Columns 17 through 20
-3.1601 -0.8161 + 2.9031i 0.4670 - 5.2595i -1.1440 + 4.9140i
Columns 21 through 24
2.7128 - 4.6691i 2.8367 - 9.6884i 7.7442 - 0.5439i -4.9958 - 0.0682i
Columns 25 through 28
-5.0758 + 1.0582i -2.6340 - 2.6884i -0.8291 - 0.9154i 0.9802 - 5.1840i
Columns 29 through 32
3.5893 - 5.2095i -2.1831 + 3.6784i -3.3874 - 0.9221i -4.0386 + 2.6955i
>> c=fft(a)
c =
Columns 1 through 4
2.8652 -4.0386 - 2.6955i -3.3874 + 0.9221i -2.1831 - 3.6784i
Columns 5 through 8
3.5893 + 5.2095i 0.9802 + 5.1840i -0.8291 + 0.9154i -2.6340 + 2.6884i
Columns 9 through 12
-5.0758 - 1.0582i -4.9958 + 0.0682i 7.7442 + 0.5439i 2.8367 + 9.6884i
Columns 13 through 16
2.7128 + 4.6691i -1.1440 - 4.9140i 0.4670 + 5.2595i -0.8161 - 2.9031i
Columns 17 through 20
-3.1601 -0.8161 + 2.9031i 0.4670 - 5.2595i -1.1440 + 4.9140i
Columns 21 through 24
2.7128 - 4.6691i 2.8367 - 9.6884i 7.7442 - 0.5439i -4.9958 - 0.0682i
Columns 25 through 28
-5.0758 + 1.0582i -2.6340 - 2.6884i -0.8291 - 0.9154i 0.9802 - 5.1840i
Columns 29 through 32
3.5893 - 5.2095i -2.1831 + 3.6784i -3.3874 - 0.9221i -4.0386 + 2.6955i
>>

Прикрепленные файлы

 NewSplitRadix.zip ( 2.12 килобайт ) Кол-во скачиваний: 35