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

 
 
> Быстрое преобразование Фурье не для степени 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
Ответов
almost
сообщение Sep 20 2011, 09:33
Сообщение #2


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

Группа: Свой
Сообщений: 199
Регистрация: 27-05-09
Из: Москва
Пользователь №: 49 648



Мне нравится как дискуссия по поводу спектра (в большинстве случаев) плавно уходит в метафизику =)
Немного отсебятенки к определению спектра:
Основополагающим стержнем держащим все здание под названием "спектр, спектральный анализ" является понятие о гармонической природе всех сигналов, т.е. "вера" (именно вера) в то, что у "чистого" синуса спектр бесконечно узкий (одна палка) и только эта "вера" позволяет нам утверждать что любой сигнал можно разложить в ряды Фурье (как правило оговаривается то, что сигнал должен быть "гармоническим", чтобы все было чисто с юридической точки зрения), а соответственно и вычислить спектр.
Ну а теперь если подумать, существуют ли гармонические сигналы? Конечно же нет, это сугубо математическая абстракция, ведь невозможно эмпирически доказать обратное (АЦП с частотой дискретизации бесконечность никогда не будет создан). Тем самым осуществляя преобразование Фурье мы самолично принимаем решение о входе в мир абстракций и уходе от "реальности". Вроде, пока что, ныне существующего понятия о сигналах вполне достаточно для реальных задач.
Могу ошибаться конечно же=)
Go to the top of the page
 
+Quote Post
ViKo
сообщение Sep 20 2011, 11:17
Сообщение #3


Универсальный солдатик
******

Группа: Модераторы
Сообщений: 8 634
Регистрация: 1-11-05
Из: Минск
Пользователь №: 10 362



Цитата(almost @ Sep 20 2011, 12:33) *
Мне нравится как дискуссия по поводу спектра (в большинстве случаев) плавно уходит в метафизику =)
...
Ну а теперь если подумать, существуют ли гармонические сигналы?
...
Могу ошибаться конечно же=)

Ничего страшного, если мы пытаемся дойти собственными размышлениями до понимания. Даже, наоборот.
Но в данном случае спектр - абсолютно точное математическое преобразование. Естественно, в реальности это преобразование можно осуществить с какой-то заданной точностью (хотя бы, разрядностью чисел ограничено). Но только и всего. Для сигнала, преобразованного в цифровую форму, появляются неточности, связанные с частотой дискретизации, количеством отсчетов в реализации, разрядностью АЦП, линейностью АЦП, джиттером частоты дискретизации, шумом.
Про гармонические сигналы. Если сигнал гармонический, т.е. содержит всего одну синусоидальную составляющую, то и спектр его - одна палка. Это, так сказать, разными словами говорится об одном. А если взять, например, прямоугольный сигнал, то в спектре у него - куча гармоник. Но вычислить его спектр тоже можно со сколь угодно большой точностью.
Синус - это не от балды взятый сигнал. Возьмите колесо, поставьте его в центр координатных осей X и Y, задайте на ободе колеса точку, и крутите это колесо с равномерной скоростью. И посмотрите, как изменяются координаты этой точки по Х и Y со временем.
И много в природе движется так, равномерно вращаясь. С определенной долей допущения, конечно.
Go to the top of the page
 
+Quote Post
almost
сообщение Sep 20 2011, 11:48
Сообщение #4


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

Группа: Свой
Сообщений: 199
Регистрация: 27-05-09
Из: Москва
Пользователь №: 49 648



Цитата(ViKo @ Sep 20 2011, 15:17) *
Ничего страшного, если мы пытаемся дойти собственными размышлениями до понимания. Даже, наоборот.

Я не в коем случае никого не обвиняю в не компетентности, дискуссию выше прочитал с удовольствием. Я просто в очередной раз подмечаю, что основы радиотехники имеют под собой достаточно хлипкий фундамент (и не я один). С другой стороны из этого хлипкого вырисовываются достаточно крепкие математические модели. Мое замечание, что все клонится к метафизике это не насмешка, я и в правду думаю что тут надо копать в сторону "А как же математика завязана с реальностью?", чтобы хоть чуть чуть взбодрить теорию о спектре, рядах фурье и прочем.

Цитата(ViKo @ Sep 20 2011, 15:17) *
Но в данном случае спектр - абсолютно точное математическое преобразование. Естественно, в реальности это преобразование можно осуществить с какой-то заданной точностью (хотя бы, разрядностью чисел ограничено). Но только и всего. Для сигнала, преобразованного в цифровую форму, появляются неточности, связанные с частотой дискретизации, количеством отсчетов в реализации, разрядностью АЦП, линейностью АЦП, джиттером частоты дискретизации, шумом.


Угумс, спектр это математическая "фича", позволяющая объединить эмпирический опыт тысячи исследователей с учетом джиттера/разрядности/ и прочего. Т.е. в моем понимание понятие спектра это инструмент, а не объяснение.

Цитата(ViKo @ Sep 20 2011, 15:17) *
Про гармонические сигналы. Если сигнал гармонический, т.е. содержит всего одну синусоидальную составляющую, то и спектр его - одна палка. Это, так сказать, разными словами говорится об одном. А если взять, например, прямоугольный сигнал, то в спектре у него - куча гармоник. Но вычислить его спектр тоже можно со сколь угодно большой точностью.
Синус - это не от балды взятый сигнал. Возьмите колесо, поставьте его в центр координатных осей X и Y, задайте на ободе колеса точку, и крутите это колесо с равномерной скоростью. И посмотрите, как изменяются координаты этой точки по Х и Y со временем.
И много в природе движется так, равномерно вращаясь. С определенной долей допущения, конечно.


Я вот не могу догнать почему именно так, почему именно синус? Моск разрывается просто =)
Есть у Юнга достаточно интересное произведение - Синхроничность называется, то же моск взрывает.
Go to the top of the page
 
+Quote Post
fontp
сообщение Sep 21 2011, 08:15
Сообщение #5


Эксперт
*****

Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183



QUOTE (almost @ Sep 20 2011, 15:48) *
Я вот не могу догнать почему именно так, почему именно синус? Моск разрывается просто =)
Есть у Юнга достаточно интересное произведение - Синхроничность называется, то же моск взрывает.


Спектральное разложение в математике, обобщенно - это собственные функции оператора.
А комплексные синусы - это собственные функции конкретного оператора сдвига. Именно поэтому синусы занимают центральное место в анализе процессов в однородном пространстве или времени. Ничего удивительного

QUOTE (GetSmart @ Sep 21 2011, 03:26) *
Более подробно о спектре. Допустим того же синуса. Спектр вообще это не одна палка. Спектр любого сигнала это целое множество "палок". Но у синуса это одна ярковыраженная палка (ненулевой длины) и все соседние нулевой длины. Вот это и есть спектр синуса при правильной мат.обработке. То есть спектр у синуса прерывный. У других не ЧМ сигналов тоже спектр в идеале прерывный, но из нескольких палок, например АМ модуляция. Мне например смешно, когда мне показывают какую-то "жуткую" картинку с пологими сторонами "палки" и говорят что это самый правильный спектр синуса, на котором видно множество правильных деталей. То есть "детали" приписывают именно сигналу, а не ограниченности данных и кривой обработке. С одной стороны, если нужно отловить на спектре только самую ярковыраженную палку или несколько их с большей точностью значения частоты, то это да, смысл есть в дополнении нулями. Но надо прекрасно понимать, что при этом возникают искажения более слабых гармоник, да и главные пики имеют амплитуду ниже оригинальной сигнала когда его период не кратен. То есть точность определения частоты возрастает, точность определения амплитуды уменьшается для пиков спектра, для слабых сигналов всё намного хуже.


Не вступая в бесполезную дискуссию, и соглашаясь с Вами в отношении "деталей", замечу только, что напротив - спектр "прерывный" бывает только в идеале, в абстракции. Физические спектры всегда непрерывны, хотя бы из-за нестабильности несущей или ограниченности времени наблюдения. "Прерывные" (линейчатые) спектры всего лишь удобная абстракция, когда от некоторых деталей можно отвлечься
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
||- - ViKo   Цитата(almost @ Sep 20 2011, 14:48) Я вот...   Sep 20 2011, 11:59
|||- - almost   Цитата(ViKo @ Sep 20 2011, 15:59) Колесо....   Sep 20 2011, 12:23
|||- - 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   Господа/Товарищи! А дробный ДПФ/БПФ никто не в...   Oct 18 2011, 08:22
|- - SPACUM   Цитата(GetSmart @ Oct 18 2011, 12:22) Зна...   Oct 18 2011, 10:49
|- - GetSmart   Цитата(SPACUM @ Oct 18 2011, 15:49) А для...   Oct 18 2011, 10:58
|- - SPACUM   Цитата(GetSmart @ Oct 18 2011, 14:58) Ник...   Oct 18 2011, 11:01
|- - GetSmart   Цитата(SPACUM @ Oct 18 2011, 16:01) http:...   Oct 18 2011, 11:21
|- - _Anatoliy   Цитата(GetSmart @ Oct 18 2011, 12:21) Как...   Oct 18 2011, 11:49
||- - 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 Текстовая версия Сейчас: 7th August 2025 - 14:07
Рейтинг@Mail.ru


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