реклама на сайте
подробности

 
 
> Алгоритм ускоренного преобразования фурье на 2**N - 1 точек
Костян
сообщение Jan 15 2013, 13:23
Сообщение #1


Знающий
****

Группа: Свой
Сообщений: 740
Регистрация: 24-07-06
Из: Minsk
Пользователь №: 19 059



Подскажите, существует ли ускоренный алгоритм преобразования фурье на 2**N - 1 точек (например 255 или 511 точек) ?
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
thermit
сообщение Jan 15 2013, 15:39
Сообщение #2


Знающий
****

Группа: Участник
Сообщений: 781
Регистрация: 3-08-09
Пользователь №: 51 730



Цитата
Костян:
Так не получиться. Для линейной свертки бы подошло, а для циклической - нет.


Чудеса, да и только... И чего это подойдет для линейной, что не подойдет для циклической?

Код
  x=[-1 -5 4 9 -8];% длина 5
  y=[5 -9 8 4 -1 4 -6 -2 9 1 -5]; %длина 11
%Результат свертки x и y имеет длину 15
%Через FFT-15
  z=ifft(fft([x zeros(1,10)]).*fft([y zeros(1,4)]));
%Через FFT-16
zz=ifft(fft([x zeros(1,11)]).*fft([y zeros(1,5)]));


z=
1.0e+002 *

-0.050000000000000
-0.160000000000000
0.570000000000000
-0.350000000000000
-1.080000000000000
1.610000000000000
-0.460000000000000
0.070000000000000
0.210000000000000
-1.400000000000000
0.660000000000000
1.260000000000000
-0.830000000000000
-0.530000000000000
0.400000000000000

zz=
1.0e+002 *

-0.050000000000000
-0.160000000000000
0.570000000000000
-0.350000000000000
-1.080000000000000
1.610000000000000
-0.460000000000000
0.070000000000000
0.210000000000000
-1.400000000000000
0.660000000000000
1.260000000000000
-0.830000000000000
-0.530000000000000
0.400000000000000
-0.000000000000000 %Этот 0 выкидывается, потому как за пределами результата

Какие проблемы?

ps
Впрочем, если охота попреодолевать трудности - флаг в руки...

Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- Костян   Алгоритм ускоренного преобразования фурье на 2**N - 1 точек   Jan 15 2013, 13:23
- - Serg76   Неужели так критично по скорости, что нельзя добав...   Jan 15 2013, 13:40
|- - Костян   QUOTE (Serg76 @ Jan 15 2013, 12:40) Неуже...   Jan 15 2013, 13:51
|- - Serg76   Цитата(Костян @ Jan 15 2013, 17:51) не в ...   Jan 15 2013, 14:00
- - thermit   ЦитатаКостян: не в скорости дело, нужна циклическа...   Jan 15 2013, 14:12
|- - Костян   QUOTE (thermit @ Jan 15 2013, 13:12) Сдел...   Jan 15 2013, 14:28
|- - Serg76   Цитата(Костян @ Jan 15 2013, 18:28) Но в ...   Jan 15 2013, 14:52
|- - Костян   QUOTE (Serg76 @ Jan 15 2013, 13:52) В том...   Jan 15 2013, 14:54
|- - Serg76   Цитата(Костян @ Jan 15 2013, 17:54) Спаси...   Jan 15 2013, 15:06
|- - Костян   QUOTE (Serg76 @ Jan 15 2013, 14:06) на са...   Jan 15 2013, 15:30
|- - Костян   QUOTE (thermit @ Jan 15 2013, 14:39) Чуде...   Jan 16 2013, 12:05
- - eugen_pcad_ru   Когда-то давно читал про это. По моему называется ...   Jan 16 2013, 05:11
- - TSerg   Можно тут взглянуть http://exfile.ru/396643   Jan 16 2013, 07:47
- - thermit   ЦитатаКостян: Возмите длину сигнала x и y равной 1...   Jan 17 2013, 08:37
- - alex_os   Цитата(thermit @ Jan 17 2013, 12:37) 2-й ...   Jan 18 2013, 04:54


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 3rd August 2025 - 09:25
Рейтинг@Mail.ru


Страница сгенерированна за 0.01386 секунд с 7
ELECTRONIX ©2004-2016