Помогите разобраться с БПФ. («Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)

[an error occurred while processing this directive] Помогите разобраться с БПФ.
(«Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)

Отправлено TheMatrix 15 ноября 2006 г. 22:53

Понимаю, что надо учиьт матчасть, но происходит это долго.
Задача: провести измерения сигнала и разложить в спектр 0-10 Гц, с интервалом 0,5 Гц. Сильно большая точность не нужна.
Вопрос: сколько отсчётов делать? В течение какого времени? Использовать алгоритм rad 4?
Процессор ADSP-2105.

Составить ответ ||| Конференция ||| Архив

Ответы

C алгоритмом разобрался, теперь по интерпретации результатов. — TheMatrix (20.11.2006 17:40 83.167.82.221, 303 байт)
- Если у Вас 64 действительных отсчета, то достаточно иметь 32 (или 33 - зависит от реализации) комплексных отсчета. Это дело реализуется при помощи функции rfft_xx (для Блекфина она, например, называется rfft_fr16() ) — homekvn (20.11.2006 21:11 212.185.161.237, пустое)
  - Алгоритм уже написал. Частоты определяет, но не могу понять, в каком порядке он их выдаёт. То ли я не понимаю "формат" спектра, то ли я таки неправильно написал алгоритм. — TheMatrix (20.11.2006 21:39 83.167.82.186, пустое)
    - Нет, похоже всё работает, единичные частоты хорошо определяет. — TheMatrix (20.11.2006 22:58 83.167.82.162, пустое)
    - А небось в битреверсном порядке все выдается. Посмотрите какая у Вас функция приводит все в нормальный порядок (что-нибудь вроде bitrev() ) — homekvn (20.11.2006 21:44 212.185.161.237, пустое)
      - Всё, уже сам понял :) — TheMatrix (20.11.2006 23:19 83.167.82.207, 84 байт)
      - Нет, как я писал, выход нормальный, т.е. на последнем этапе расшифровка проводится. Тогда такой вопрос: — TheMatrix (20.11.2006 23:01 83.167.82.162, 110 байт)
        
        Тригонометрию советую вспомнить (+) — SM (21.11.2006 00:26 213.141.159.26, 108 байт)
  - Такая функция, как правило, и работает быстрее или за то же время выполняет сразу два БПФа. — homekvn (20.11.2006 21:26 212.185.161.237, пустое)
- Это не амплитуда и фаза, а действительная и мнимая часть. Амплитуда будет модуль, а фаза - аргумент. А вообще странноватый у Вас вопрос - реализовывать БПФ, не зная что оно выдает :) — SM (20.11.2006 19:46 85.21.237.237, пустое)
  - Да, да, ступил. Действительно модуль. Всё-таки, каков должен быть результат при преобразовании n периодов косинуса? — TheMatrix (20.11.2006 20:11 83.167.82.210, пустое)
Наверно, последний вопрос по алгоритму. — TheMatrix (20.11.2006 12:29 83.167.93.188, 1003 байт)
Наверно, последний вопрос по алгоритму. — TheMatrix (20.11.2006 12:28 83.167.93.188, 1003 байт)
Т.е. перед применением БПФ надо 2 секунды собирать данные? Можно быстрее? — TheMatrix (17.11.2006 18:35 83.167.93.225, пустое)
- Нет. Сами задали разрешение 0.5 Гц. — andy_P (17.11.2006 18:59 unknown, 80.82.63.185, пустое)
  - Это так, но разве нельзя по небольшому количеству выборок судить о частоте сигнала? Например половина синусоиды - это менее периода. Я уже понял, как писали - можно частоту выборок повышать кратно, т.е. кратно уменьшать время. Я правильно понял? — TheMatrix (17.11.2006 19:56 83.167.92.249, 329 байт)
    - Можно, но если обладаете достаточной априорной информацией о характере сигнала. Например, что это один синус. Или еще что-то заранее известное. Если Вы о сигнале ничего кроме его полосы не знаете, то низзя. — SM (18.11.2006 01:08 213.141.159.26, пустое)
    - Так Вам что надо? О частоту сигнала измерить (!!!одна синусоида!!!) или спектр сигнала получить с шагом 0.5 Hz (исходный пост) - читай например разрешить две сигусоиды отличающиеся на 0.5 Hz? — andy_P (17.11.2006 20:46 unknown, 80.82.63.185, 381 байт)
      - Эти посты я написал исходя из поста homekvn. — TheMatrix (17.11.2006 21:16 83.167.93.184, пустое)
      - Надо спектр. — TheMatrix   (17.11.2006 21:14 83.167.93.184, 230 байт)
        
        По затратам это все мелочевка (частота дискретизации очень низкая (ее Вам SM с homekvn писали )). Но спектр с заданными Вами разрешением быстрее не измерить не получится. К тому же чтобы оценку спектра хорошую получить нужно усреднение, а это потребует еще времени. — andy_P   (18.11.2006 10:42 80.82.36.126, пустое)
        
        Мне всё никак не ответят на вопрос: homekvn писал о том, что можно КРАТНО повышать частоту измерений. — TheMatrix   (18.11.2006 10:54 83.167.93.227, 127 байт)
        
        Нет не так. Можно понизить чатоту дискретизации и делать меньше вычислений, если у Вас она (частота дискретизации) сильно больше требуемой полосы. У вас полоса 10 Hz, минимальная частота дисертизации должна быть больше 20. Если у Вас она выше, то можгно децимировать. На времени реакции (времени за которое данные собираются) это никак не скажется. Просто обсчитывать будете меньше отсчетов. — andy_P   (18.11.2006 11:32 80.82.36.45, пустое)
        
        Усреднение - это ладно, главное чтобы время реакции системы не было таким большим. — TheMatrix   (18.11.2006 10:50 83.167.93.227, пустое)
        
        Да не вызовет это уже никаких затрат. Затраты разрешением определяются. Впрочем (+) — homekvn   (17.11.2006 23:42 84.146.24.133, 386 байт)
        
        А задача такая. — TheMatrix   (18.11.2006 09:38 83.167.92.229, 473 байт)
        
        В Вашем случае я бы отказался от БПФ вообще, а сделал бы так, как опишу здесь. (+) — homekvn   (18.11.2006 12:58 84.146.49.20, 3716 байт)
        
        Про банк фильтров - прям как у меня описали... ;-) А warped-FFT - я очень хочу! Поделитесь пожалуйста материалами... (+) — ASergej_R19   (18.11.2006 21:46 80.250.160.170, 60 байт)
        
        Отправил кое-что. Проверьте почту. Если не дошло, напишите, повторю. — homekvn   (18.11.2006 22:36 84.146.49.20, пустое)
        
        Все дошло... Большое спасибо... (-) — ASergej_R19   (18.11.2006 22:45 80.250.160.170, пустое)
        
        Матлабовский код также отправил сегодня — homekvn   (20.11.2006 16:15 212.185.161.237, пустое)
        
        Да, я получил - большое спасибо... (-) — ASergej_R19   (20.11.2006 16:21 80.250.160.170, пустое)
        
        Спасибо, но есть ограничения на задание. — TheMatrix   (18.11.2006 13:38 83.167.93.229, 508 байт)
        
        Если память данных только полкило, то кроме как БПФ на 40 (а лучше 64) отсчетов (при частоте дискретизации 20 Гц) в голову ничего не лезет. Тем более, что частоту дискретизации входного сигнала выбираете Вы сами. Это большое преимущество. (+) — homekvn   (18.11.2006 20:00 84.146.49.20, 1217 байт)
        
        Памяти 0.5 k 16-разрядных слов... — TheMatrix   (18.11.2006 22:30 83.167.93.167, 218 байт)
        
        Храните таблицу в памяти программ. Размер таблицы - длина преобразования N (в этом случае должна быть кратна 4) деленная на 4 плюс 1. Т.е. храним синус (или косинус) на интервале от 0 до пи/2 с шагом 2*пи/N. Впрочем (+) — homekvn   (18.11.2006 22:41 84.146.49.20, 234 байт)
        
        Добавлю (+) — homekvn   (18.11.2006 13:12 84.146.49.20, 914 байт)
        
        А у нас еще суббота:) — Бяка   (18.11.2006 13:10 83.237.18.245, пустое)
        
        Да и у нас на два часа позже столицы нашей Родины все наступает :-) А вообще.. — homekvn   (18.11.2006 13:15 84.146.49.20, 75 байт)
        
        Если у нас суббота, то у вас тем более воскресенья еще нет:) — Бяка   (18.11.2006 14:21 83.237.18.245, пустое)
        
        Столица нашей Родины - Москва, а не Ванкувер ;-) — homekvn   (18.11.2006 15:59 84.146.49.20, пустое)
Ответ (+) — SM (16.11.2006 00:41 213.141.159.26, 651 байт)
Ответ (+) — homekvn (16.11.2006 00:27 84.146.52.150, 1975 байт)
- Про БПФ — TimMATtU (17.11.2006 11:31 217.71.134.23, 217.71.128.84, 552 байт)
  - Ответ (+) — homekvn (17.11.2006 11:58 212.185.161.237, 919 байт)
    - В догонку (+) — homekvn (17.11.2006 12:07 212.185.161.237, 926 байт)
      - Не забывайте еще, что работаем мы не с непрерывными сигналами, а с последовательностями (и последовательностями ограниченными по спектру) (+) — homekvn   (17.11.2006 12:34 212.185.161.237, 207 байт)
        
        Все наоборот. Последовательности благодаря дискретному характеру во времени бесконечны и периодичны по спектру, а значит неограничены. В отличие от непрерывных сигналов, которые могут быть ограничеными. — SM   (18.11.2006 14:58 213.141.159.26, пустое)
        
        Выразился, возможно, не совсем понятно. Имел в виду то, что рассматривая периодическое продолжение последовательности ограниченной длины N, мы тем самым автоматически ограничиваем ее спектр. И ограничиваем величиной N/2 (или если не рассматривать отрицательных отсчетов, то величиной N). В случае (+) — homekvn   (18.11.2006 16:50 84.146.49.20, 367 байт)
        
        Стормозил маленько. Написал бред (+) — homekvn   (19.11.2006 20:51 84.146.18.63, 367 байт)
        
        Это почему в физическом смысле не удовлетворяет (+) — SM   (20.11.2006 00:08 213.141.159.26, 137 байт)
        
        Угу. Ваша правда. Совсем торможу. — homekvn   (20.11.2006 01:58 84.146.18.63, пустое)
        
        Ответ: — SM   (18.11.2006 14:57 213.141.159.26, пустое)
        
        О! А бывают последовательности "по спектру неограниченные"? — st256   (18.11.2006 13:45 217.151.231.218, 307 байт)
        
        А по поводу ограниченности или неограниченности спектра последовательности предлагаю иное видение этого вопроса. Для этого прочитайте пост SM аккурат над Вашим. Впрочем, я далее уже по этому поводу тоже высказался. — homekvn   (18.11.2006 16:58 84.146.49.20, пустое)
        
        Не возможно отличить две частоты, попадающие в один интервал. Вот и все. :-) — homekvn   (18.11.2006 16:39 84.146.49.20, пустое)
        
        Так вопрос, если я правильно понял, как раз с этим и связан - Какое отношение имеет коэфициент ДПФ к "настоящему" отсчёту непрерывного преобразования? — fontp   (17.11.2006 13:21 192.168.10.211, 213.171.46.37, пустое)
        
        А я и написал прямо - никакого. ДПФ "похож" на ПФ только в тех смыслах, о которых я сказал. Точный результат (и то, с точностью до умножения на непрерывную дельта-функцию) можно получить лишь для периодических сигналов, причем только в том случае, если длину ДПФ мы взяли в точности кратную периоду сигнала. — homekvn   (17.11.2006 13:31 212.185.161.237, пустое)
        
        Мда... Это ничего если я Вас поправлю? А пишите Вы много, но не всегда правильно. — st256   (18.11.2006 13:51 217.151.231.218, 166 байт)
        
        Не согласен. (+) — homekvn   (18.11.2006 16:37 84.146.49.20, 563 байт)
        
        Если у Вас плохо с психикой, то предупредите, пожалуйста, заранее. И ни в коем случае, по примеру некоторых недоучек не пишите в Лигу Наций. Ее уже давно нету. — st256   (18.11.2006 18:45 217.151.231.218, 1016 байт)
        
        Для таких невинных забав, как общение на форуме моя психика вполне пригодна. Покорнейше благодарю за предусмотрительность и обеспокоенность. :-) Ну а по теме (+) — homekvn   (18.11.2006 19:39 84.146.49.20, 1379 байт)
        
        Давайте так, — st256   (18.11.2006 20:06 217.151.231.218, 714 байт)
        
        st256, в следующий раз, прежде чем вспоминать мой ник всуе, просто сразу пошлите себя куда-нибудь. Хоть на Тибет, что ли. Ладно? :-) А то мне не досуг. — -=ВН=-   (20.11.2006 12:49 193.125.71.140, пустое)
        
        Это не я вспоминаю. Это о Вас вспоминают другие люди. Вот я их сразу на Вас и посылаю. — st256   (20.11.2006 16:49 217.151.231.218, пустое)
        
        Не, мил человек, у других я ссылок имени меня этой ветке не наблюдаю. Так что пошлите лучше себя. На Тибет. Или в Тибет? В общем - в Гималаи. Это будет справедливо. И не ешьте на ночь сырых помидоров. — -=ВН=-   (20.11.2006 17:38 193.125.71.140, пустое)
        
        Я не эту ветку имел ввиду. Я имел ввиду прославление Вас, а так же постановку мне в пример Вашего светлого образа в письмах иногда до меня доходящих. Ну вот авторов таких писем я сразу по определенному и отправляю: т.е. прямо к Вам. — st256   (20.11.2006 18:01 217.151.231.218, пустое)
        
        Не нужно гротесков. В письме, которое я Вам направил более месяца назад, нет никакого "прославления". А что я Вам написал тогда, готов повторить еще не единожды. — homekvn   (20.11.2006 18:36 212.185.161.237, пустое)
        
        Странный Вы, Корея на Вас повлияла ужасно плохо... Мне чихать 3 раза на содержание писем Вами получаемых и отправляемых. Так же, как и чихать на прославление Вашего светлого образа в этой конфе, если таковое было. Вам было предложено послать себя в связи с содержанием Вашего поста. Что непонятного? — -=ВН=-   (20.11.2006 18:26 193.125.71.140, пустое)
        
        В смысле по определенному маршруту к Вам отправляю. — st256   (20.11.2006 18:03 217.151.231.218, пустое)
        
        Ясно. (+) — homekvn   (18.11.2006 21:39 84.146.49.20, 1964 байт)
        
        Я прикрываю дискуссию. — st256   (19.11.2006 09:06 217.151.231.218, 567 байт)
        
        А дискуссии и не было, я вам пытался ее предложить, Вы предпочли ей пассажи о моем психическом состоянии и моем низком уровне образования. (+) — homekvn   (19.11.2006 20:38 84.146.18.63, 1170 байт)
        
        Поправочка - с периодом T у Вас полный порядок. — st256   (18.11.2006 18:54 217.151.231.218, 35 байт)
        
        Ну вот давайте подискутируем, если у Вас еще не отпало желание, тогда и посмотрим, у кого с теорией лучше. — homekvn   (18.11.2006 19:41 84.146.49.20, пустое)
        
        Мне тяжело с Вами дискутировать. Вам показываешь, где Вы неправы, а Вы - в партком с петицией на мое аморальное поведение. Хоть партком уже давно отменили, но тем не менее оказаться в центре очередного скандала мне не хочется. Даже то, что Вам никто не ответит, а меня никто не накажет, мало успокаивает: представьте, что Вы на экскурсии в сумсшедшем доме - вроде Вас это не касается, но все равно боязно... — st256   (18.11.2006 19:53 217.151.231.218, пустое)
        
        Ну не хотите, не дискутируйте. Просто проверьте, например, в Матлабе. Код для проверки займет 2 строчки (с выводом графика - три) — homekvn   (18.11.2006 20:03 84.146.49.20, пустое)

Перейти к списку ответов ||| Конференция ||| Архив ||| Главная страница ||| Содержание