реклама на сайте
подробности
 
 Wiki Wiki  Photo Photo  Forum Forum  Reviews Reviews  Help (!) Help (!)
 
> Быстрое преобразование Фурье не для степени 2, как довести число отсчетов до стпени двойки
Sergey Makarov
сообщение Sep 14 2011, 12:45
Сообщение #1


Частый гость
**

Группа: Свой
Сообщений: 147
Регистрация: 19-10-05
Из: Россия, г.Чебоксары
Пользователь №: 9 845
Добрый день!
Стоит задача осуществить быстрое преобразовние Фурье для 500 отсчетов. Аппаратно DSP поддерживает быстрое преобразование по модулю 2. Вопрос, как довести кол-во точек до 512? Достаточно ли заполнить недостающие значения нулями? Сильно ли исказится результат преобразования после такого допущения? Если можно ссылкой на математическое обоснование поделитесь плз. Или все таки придется до 512 точек дополнять интерполировав существующие значения? Кто чем может помогите плз.
Go to the top of the page
 
+Quote Post
 
 Start new topic
Ответов
GetSmart
сообщение Oct 18 2011, 08:22
Сообщение #2


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753
Господа/Товарищи!
А дробный ДПФ/БПФ никто не видел? Мне вот понадобилось (особено для малых N [16..256]). В частности хотелось бы определять промежуточные дробные частоты с веществ. значениями 0.5, 0.25, 0.75 и др. от 1/T.
Знаю точно, что дробный ДПФ работает не хуже целого по части нахождения пиков гармоник. Испытывал на медленном алгоритме. Но вдруг есть побыстрее.

Сообщение отредактировал GetSmart - Oct 18 2011, 08:25


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
SPACUM
сообщение Oct 18 2011, 10:49
Сообщение #3


Частый гость
**

Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531
Цитата(GetSmart @ Oct 18 2011, 12:22) *
Знаю точно, что дробный ДПФ работает не хуже целого по части нахождения пиков гармоник. Испытывал на медленном алгоритме. Но вдруг есть побыстрее.

Согласен.
Похоже выше 7й гармоники дополнение нулями работает удовлетворительно. А ниже нужны особые методы.
У меня еще и шумы, так что менее 2000 точек мне неинтересно, а больше долго.
А для 30 точек без шума может метод барона Прони подойдет?
Ничего другого по вопросу в литературе не встречал.
Еще мне кажется для частот в начале спектра Фурье и в конце спектра Фурье предел Крамера-Рао недостижим с запасом.
И для Вашего случая особенно. А у Вас получилось близко?


--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Oct 18 2011, 10:58
Сообщение #4


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753
Цитата(SPACUM @ Oct 18 2011, 15:49) *
А для 30 точек без шума может метод барона Прони подойдет?

Цитата
Барон де Прони(Гаспар Рише) был другом Наполеона и одним из основоположников метрической системы мер.
Имя барона Гаспара де Прони входит в список 72 имён на Эйфелевой башне. Он действительно разработал
свой метод спектрального анализа при котором вычислялись истинные амплитуды и частоты и не было
никакой спектральной утечки. Соответствующая программа на ФОРТРАНе приведена в книге
С.Л. Марпл-мл. "Цифровой спектральный анализ и его приложения" (гл. 11). Сигнал приближался
затухающими синусоидами. Этот метод применяется не везде из-за большого объема расчетов.
Однако результат восхитительный.

Никогда о нём не слышал. Будем почитать.

Цитата(SPACUM @ Oct 18 2011, 15:49) *
И для Вашего случая особенно. А у Вас получилось близко?

Скорее нет. Хотя в 32 точках всего 16 гармоник и я затрудняюсь правильно ответить.


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
SPACUM
сообщение Oct 18 2011, 11:01
Сообщение #5


Частый гость
**

Группа: Участник
Сообщений: 161
Регистрация: 22-06-09
Из: Москва
Пользователь №: 50 531
Цитата(GetSmart @ Oct 18 2011, 14:58) *
Никогда о нём не слышал. Будем почитать.

http://www.padabum.com/d.php?id=3411
стр. 365.

Сообщение отредактировал SPACUM - Oct 18 2011, 11:05


--------------------
Ты можешь знать все что угодно, но пока ты не доказал это на практике, ты не знаешь ничего!© Ричард Бах
Go to the top of the page
 
+Quote Post
GetSmart
сообщение Oct 18 2011, 11:21
Сообщение #6


.
******

Группа: Участник
Сообщений: 4 005
Регистрация: 3-05-06
Из: Россия
Пользователь №: 16 753
Цитата(SPACUM @ Oct 18 2011, 16:01) *
http://www.padabum.com/d.php?id=3411
стр. 365.

Какой-то он избыточный. Мне нужен довольно простой результат. Даже не "оценивание". А просто спектр с определённой степенью точности и с промежуточными частотами. При разрешении 1/8 бина БПФ точность амплитуды уже будет не менее 95% (точно не измерял, но не ниже 90%). Мне нужна примерно такая точность амплитуды и частОты 1/8 целого. И всё. Если амплитуда синуса ниже шума, то вообще не важен результат. Я пока что планировал итеративно вычислять ДПФ дробных частот в районе пиков на БПФ до достижения нужной точности частоты. Но с БПФ плохо то, что он амплитуду (пики) выдаёт неверную для дробных частот. Придётся анализировать окресности всех бинов до 1/2 от самого высокого. Хорошо хоть точек мало.

Сообщение отредактировал GetSmart - Oct 18 2011, 11:37


--------------------
Заблуждаться - Ваше законное право :-)
Go to the top of the page
 
+Quote Post
_Anatoliy
сообщение Oct 18 2011, 11:49
Сообщение #7


Утомлённый солнцем
******

Группа: Свой
Сообщений: 2 646
Регистрация: 15-07-06
Из: г.Донецк ДНР
Пользователь №: 18 832
Цитата(GetSmart @ Oct 18 2011, 12:21) *
Какой-то он избыточный. Мне нужен довольно простой результат. Даже не "оценивание". А просто спектр с определённой степенью точности и с промежуточными частотами. При разрешении 1/8 бина БПФ точность амплитуды уже будет не менее 95% (точно не измерял, но не ниже 90%). Мне нужна примерно такая точность амплитуды и частОты 1/8 целого. И всё. Если амплитуда синуса ниже шума, то вообще не важен результат. Я пока что планировал итеративно вычислять ДПФ дробных частот в районе пиков на БПФ до достижения нужной точности частоты. Но с БПФ плохо то, что он амплитуду (пики) выдаёт неверную для дробных частот. Придётся анализировать окресности всех бинов до 1/2 от самого высокого. Хорошо хоть точек мало.

А интерполировать по трём бинам не пробовали?Я когда-то игрался с этим,получались неплохие результаты но дальше модели не пошло.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- Sergey Makarov   Быстрое преобразование Фурье не для степени 2   Sep 14 2011, 12:45
- - fontp   Нулей добавить - нормально. Фактически это означае...   Sep 14 2011, 13:58
- - SPACUM   Цитата(Sergey Makarov @ Sep 14 2011, 16:4...   Sep 14 2011, 14:32
|- - getch   Листинг программы на C++ (N - любое)   Sep 14 2011, 17:51
- - Alexey Lukin   Чем больше добавляется нулей, тем точнее получаетс...   Sep 14 2011, 19:34
|- - GetSmart   Цитата(Alexey Lukin @ Sep 15 2011, 00:34)...   Sep 14 2011, 19:39
|- - thermit   Цитата(Alexey Lukin @ Sep 14 2011, 22:34)...   Sep 14 2011, 20:29
- - Alexey Lukin   thermit, вы что-то не в кассу... Результат ДПФ ...   Sep 14 2011, 21:26
|- - GetSmart   Цитата(Alexey Lukin @ Sep 15 2011, 02:26)...   Sep 14 2011, 21:57
- - Alexey Lukin   Ничего не искажается! Наоборот, становится вид...   Sep 14 2011, 22:55
|- - GetSmart   Цитата(Alexey Lukin @ Sep 15 2011, 03:55)...   Sep 16 2011, 20:59
|- - Alexey Lukin   Цитата(GetSmart @ Sep 17 2011, 00:59) ДПФ...   Sep 17 2011, 02:39
|- - GetSmart   Цитата(Alexey Lukin @ Sep 17 2011, 07:39)...   Sep 17 2011, 08:37
|- - Alexey Lukin   Цитата(GetSmart @ Sep 17 2011, 12:37) ДПФ...   Sep 17 2011, 15:26
|- - thermit   Цитата(Alexey Lukin @ Sep 17 2011, 18:26)...   Sep 17 2011, 16:04
|- - GetSmart   Цитата(Alexey Lukin @ Sep 17 2011, 20:26)...   Sep 17 2011, 16:18
|- - Alexey Lukin   Цитата(GetSmart @ Sep 17 2011, 20:18) Вам...   Sep 17 2011, 16:58
|- - GetSmart   Цитата(Alexey Lukin @ Sep 17 2011, 21:58)...   Sep 17 2011, 17:16
- - Sergey Makarov   Спасибо большое за ответы! Кое что уже началос...   Sep 15 2011, 05:54
|- - SPACUM   Цитата(Sergey Makarov @ Sep 15 2011, 09:5...   Sep 15 2011, 06:22
|- - Sergey Makarov   Цитата(SPACUM @ Sep 15 2011, 10:22) Это з...   Sep 15 2011, 06:43
|- - SPACUM   Цитата(Sergey Makarov @ Sep 15 2011, 10:4...   Sep 15 2011, 07:05
- - thermit   ЦитатаAlexey Lukin: Результат ДПФ — это и есть спе...   Sep 15 2011, 06:45
- - bookevg   В книге Гольденберга 1985г на стр.220, если не вып...   Sep 15 2011, 07:02
|- - fontp   QUOTE (bookevg @ Sep 15 2011, 11:02) В кн...   Sep 15 2011, 07:18
|- - SPACUM   Цитата(bookevg @ Sep 15 2011, 11:02) дать...   Sep 15 2011, 08:22
- - Sergey Makarov   как все сложно. Многое надо обновить в памяти, инф...   Sep 15 2011, 08:08
|- - bahurin   Почитайте для начала тута и здеся.   Sep 15 2011, 08:28
|- - SPACUM   Цитата(Sergey Makarov @ Sep 15 2011, 12:0...   Sep 15 2011, 09:04
- - thermit   ЦитатаAlexey Lukin: ДПФ это не свёртка. ДПФ это ум...   Sep 17 2011, 10:30
- - GetSmart   Alexey Lukin, я кажется понял ваше непонимание про...   Sep 17 2011, 11:49
- - Alexey Lukin   Слово "спектр" имеет множество значений,...   Sep 17 2011, 16:11
|- - thermit   Цитата(Alexey Lukin @ Sep 17 2011, 19:11)...   Sep 17 2011, 19:35
- - Alexey Lukin   Дальнейший спор мне представляется бессмысленным. ...   Sep 17 2011, 18:49
|- - GetSmart   Цитата(Alexey Lukin @ Sep 17 2011, 23:49)...   Sep 17 2011, 19:29
- - Alexey Lukin   thermit, вы не томите, подскажите ваше определение...   Sep 17 2011, 19:53
- - thermit   inf S(w)=SUM( x(n) * exp(-j*w*n) ) n=-inf...   Sep 17 2011, 20:43
|- - Alexey Lukin   О! Наконец-то Гуру-киллер сказал умную вещь...   Sep 17 2011, 21:12
- - GetSmart   Щас скажу умную вещь Спектры бывают разные. Поэто...   Sep 17 2011, 20:59
- - GetSmart   Тогда скажу ещё одну умную вещь. Можно не хлопать ...   Sep 17 2011, 21:29
- - ViKo   Сомневающимся, что дополнение нулями улучшает пред...   Sep 18 2011, 07:20
|- - Alexey Lukin   Цитата(ViKo @ Sep 18 2011, 11:20) для Ale...   Sep 18 2011, 15:36
- - Sergey Makarov   А можно поближе к земным материям? я чет из всего ...   Sep 20 2011, 05:00
|- - ViKo   Цитата(Sergey Makarov @ Sep 20 2011, 08:0...   Sep 20 2011, 06:38
||- - Sergey Makarov   Цитата(ViKo @ Sep 20 2011, 10:38) А вы кн...   Sep 20 2011, 07:06
|- - SPACUM   Цитата(Sergey Makarov @ Sep 20 2011, 09:0...   Sep 20 2011, 12:43
|- - Sergey Makarov   Цитата(SPACUM @ Sep 20 2011, 16:43) Не по...   Sep 20 2011, 13:42
- - almost   Мне нравится как дискуссия по поводу спектра (в бо...   Sep 20 2011, 09:33
|- - ViKo   Цитата(almost @ Sep 20 2011, 12:33) Мне н...   Sep 20 2011, 11:17
||- - almost   Цитата(ViKo @ Sep 20 2011, 15:17) Ничего ...   Sep 20 2011, 11:48
||- - ViKo   Цитата(almost @ Sep 20 2011, 14:48) Я вот...   Sep 20 2011, 11:59
|||- - almost   Цитата(ViKo @ Sep 20 2011, 15:59) Колесо....   Sep 20 2011, 12:23
||- - fontp   QUOTE (almost @ Sep 20 2011, 15:48) Я вот...   Sep 21 2011, 08:15
|||- - GetSmart   Цитата(fontp @ Sep 21 2011, 13:15) ... за...   Sep 23 2011, 20:51
||- - Ofei   Цитата(almost @ Sep 20 2011, 15:48) ... Я...   Sep 28 2011, 11:24
|- - GetSmart   Цитата(almost @ Sep 20 2011, 14:33) Мне н...   Sep 20 2011, 23:26
- - thermit   Чем спектр синуса отличается от спектра косинуса?   Sep 24 2011, 10:19
|- - GetSmart   Цитата(thermit @ Sep 24 2011, 15:19) Чем ...   Sep 24 2011, 19:17
- - Дмитрий_Б   Почему именно гармонические функции из всех ортого...   Sep 28 2011, 17:54
- - Alexey Lukin   Это хорошая причина. Есть и другие: например, наше...   Sep 29 2011, 08:43
- - GetSmart   Почитал все посты тредстартера. Не понял что надо ...   Oct 5 2011, 19:25
- - Kluwert   Цитата(Sergey Makarov @ Sep 14 2011, 16:4...   Oct 7 2011, 06:52
- - Alexey Lukin   Не забывайте, что свёртка с sinc происходит при лю...   Oct 7 2011, 07:48
- - andyp   ЦитатаДобрый день! Стоит задача осуществить бы...   Oct 8 2011, 11:49
- - Sergey Makarov   Цитата(GetSmart @ Oct 5 2011, 23:25) Почи...   Oct 10 2011, 13:54
|- - SPACUM   Цитата(Sergey Makarov @ Oct 10 2011, 17:5...   Oct 10 2011, 14:46
|- - andyp   Цитатавобщем путем моделирования стало ясно, что д...   Oct 10 2011, 17:50
|- - ViKo   Цитата(Sergey Makarov @ Oct 10 2011, 16:5...   Oct 11 2011, 07:56
- - Sergey Makarov   Цитата(SPACUM @ Oct 10 2011, 18:46) Тк. а...   Oct 11 2011, 11:44
|- - andyp   Цитата(Sergey Makarov @ Oct 11 2011, 15:4...   Oct 11 2011, 12:06
||- - GetSmart   Цитата(_Anatoliy @ Oct 18 2011, 16:49) А ...   Oct 18 2011, 11:57
||- - _Anatoliy   Цитата(GetSmart @ Oct 18 2011, 12:57) Над...   Oct 18 2011, 12:05
|- - SPACUM   Цитата(GetSmart @ Oct 18 2011, 15:21) Как...   Oct 18 2011, 14:59
- - Alexey Lukin   А что такое дробный ДПФ? Если надо определять коэф...   Oct 18 2011, 10:00
|- - GetSmart   Цитата(Alexey Lukin @ Oct 18 2011, 15:00)...   Oct 18 2011, 10:48
|- - GetSmart   Цитата(Alexey Lukin @ Oct 18 2011, 15:00)...   Oct 18 2011, 14:48
|- - Alexey Lukin   Цитата(GetSmart @ Oct 18 2011, 18:48) Пра...   Oct 18 2011, 15:41
|- - GetSmart   Цитата(Alexey Lukin @ Oct 18 2011, 20:41)...   Oct 18 2011, 16:34
|- - SPACUM   Цитата(GetSmart @ Oct 18 2011, 20:34) У м...   Oct 18 2011, 18:41
|- - Alexey Lukin   Цитата(GetSmart @ Oct 18 2011, 20:34) Тол...   Oct 18 2011, 19:02
|- - fontp   QUOTE (Alexey Lukin @ Oct 18 2011, 23:02)...   Oct 19 2011, 13:03
|- - GetSmart   Цитата(fontp @ Oct 19 2011, 18:03) Шум са...   Oct 19 2011, 15:48
|- - fontp   QUOTE (GetSmart @ Oct 19 2011, 19:48) В о...   Oct 20 2011, 07:20
- - ViKo   Так спектр и усреднить можно. По множеству реализа...   Oct 18 2011, 13:47
- - GetSmart   Цитата(SPACUM @ Oct 18 2011, 23:41) ... 1...   Oct 18 2011, 19:24
- - Alexey Lukin   Цитата(GetSmart @ Oct 18 2011, 23:24) ДПФ...   Oct 18 2011, 20:12
|- - GetSmart   Цитата(Alexey Lukin @ Oct 19 2011, 01:12)...   Oct 20 2011, 13:03
|- - SPACUM   Цитата(GetSmart @ Oct 20 2011, 17:03) Поч...   Oct 21 2011, 07:25
||- - GetSmart   Цитата(SPACUM @ Oct 21 2011, 12:25) Посмо...   Oct 21 2011, 12:45
||- - SPACUM   CRLB можно вычислить по этой формуле: RMS(df) / Fs...   Oct 21 2011, 13:11
||- - GetSmart   Цитата(SPACUM @ Oct 21 2011, 18:11) CRLB ...   Oct 21 2011, 13:53
||- - SPACUM   Цитата(GetSmart @ Oct 21 2011, 17:53) А н...   Oct 21 2011, 15:11
|- - Alexey Lukin   Цитата(GetSmart @ Oct 20 2011, 17:03) Сей...   Oct 21 2011, 08:15
- - SPACUM   Цитата(GetSmart @ Oct 18 2011, 23:24) 5. ...   Oct 19 2011, 06:06
2 страниц V   1 2 >

 

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

 

RSS Текстовая версия Сейчас: 22nd June 2025 - 01:12
Рейтинг@Mail.ru


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