Перейти к содержанию
    

Спектральный анализ Фурье

3 минуты назад, IhorU сказал:

Я уже говорил что пишу с ошибками. Так я еще не говорил что в детстве часто дрался и особенно любил бить снобов. 

Ремесленникам привет

Видимо вам, молодой человек, на какой-то другой форум.

Здесь практически одни, насколько я понимаю, с вашей точки зрения, снобы -  и подозреваю - практически все как минимум в 2 раза старше вас.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

12 минут назад, a123-flex сказал:

Видимо вам молодой человек, на другой форум.

Здесь практически одни, насколько я понимаю, с вашей точки зрения, снобы -  и подозреваю - практически все как минимум в 2 раза старше вас.

Пожалуйтесь модератору за неуважение к сединам, и нарушение правил форума.

Пока, ВЫ мне надоелис сударь, толку с вас все равно никакого, кроме нравооучений времен очаковских и покорения Крыма

Изменено пользователем IhorU

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

1 час назад, IhorU сказал:

Хоть по прикалуюсь, если не получу поддержку

А какого рода вы хотите получить поддержку? Давайте я вас поддержу - женщины более жалостливые :).

 

Выбираете массив длиной в целую степень двойки, чтобы FFT на нем работал. Скажем, на 2048 точек. Заполняете его через АЦП измерениями. Вызываете процедуру FFT для преобразования в частотную область.  Лучше выбрать ту, что работает не с комплексными, а с реальными входными данными, строя частотный спектр на том же месте (т.к. у МК памяти мало).

Потом вырезаете из частного спектра зоны в тех областях, где обычно бывают специфические биоритмы. Вырезать лучше не резким отсечением на краях, а по возможности плавно, умножая спектр на колоколообразную функцию с вершиной над точкой выделяемой частоты. Лично я  предпочитаю "окно Ханна" (об этом здесь: https://ru.wikipedia.org/wiki/Оконное_преобразование_Фурье), поскольку оно легко вычисляется - верхняя половинка косинуса. Затем, то что получилось, возвращаете назад  в действительную область функцией обратного преобразования Фурье (IFFT), той самой, что идет в паре с той функций прямого преобразования FFT, которым пользовались в начале. А дальше глядите на то, что получилось, и наслаждаетесь :).

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Годах этак в 2013-2014 делал я проект кардиографа, с анализом снимаемых данных. Было там множество и фильтров и БПФ. И отдельно - анализаторы этих данных, выявляющие различные патологии в сердечном ритме. Всю теорию и математику делал не сам - мне алгоритмы (даже с исходниками) предоставили коллеги из Будапешта. Они эту теорию нарабатывали и клинически обкатывали до этого несколько лет (почти 10 лет!) на базе какой-то больницы и под руководством тамошнего профессора с соответствующей специализацией (кроме него в команде были и другие врачи). Обкатывали по базе накопленных данных в клинике, да и с реальными пациентами насколько помню имели дело. Обкатывали алгоритмы на PC (win).

Следующим этапом работ было: перенести эти обкатанные алгоритмы для выполнения непосредственно на устройстве, снимающем кардиограммы.

Сначала, по плану работ, я должен был только написать код, снимающий данные и управляющий всей периферией (драйвера) и предоставляющий поток входных сырых данных с АЦП коду мат.обработки и забирающий результаты обработки и передающий их в канал. Но, так как девайс был на базе DSP, а коллеги, писавшие и отлаживавшие это всё на PC, не имели никакого опыта разработки ни под DSP ни вообще под embedded, то у них не получилось: не хватало производительности DSP-ядра для обработки. Сильно не хватало - в несколько раз.

Поэтому адаптировать код под DSP было предложено мне. Я сам в предметной области (кардиологии) ничего не шарил (да и сейчас не особо). Но будапештсткие коллеги сделали для меня подробное описание своих алгоритмов, очень качественное. И предоставили полные исходники. Я это перенёс на девайс. Потом оптимизировал (для ускорения). Хватило си-оптимизации, но с активным использованием intrinsic-функций и на чистом си. Код хорошо оптимизировался компилятором (проверялось и допиливалось и допиливалось много итераций по ассемблерным листингам).

Итого: удалось добиться, что при включении всех опций обработки и анализа (макс. вычислительная загрузка), загрузка DSP-ядра была около 20%.

А теперь внимание! Делалось это всё на OMAP-L137 с тактовой ~295МГц (не максимальная конечно) для обоих ядер (ARM+DSP). Обработка на DSP-ядре велась в double-значениях (на этом ядре double - поддерживается системой команд, 2 такта на double-MAC-операцию если не изменяет склероз). DSP-ядро занималось исключительно математикой прикладной обработки данных. Все манипуляции с железом и потоками данных делало ARM-ядро (чтоб разгрузить DSP).

 

Итого: DSP-ядро C674x с аппаратным double при тактовой ~295МГц загружено на ~20%. А вы говорите абдурина.  :biggrin::biggrin::biggrin:

+теория и консультации в прикладной области от команды врачей. А не пиво-инет.

 

PS: ОЗУ там потребовалось немного для всей этой обработки: всё влезло во внутреннее ОЗУ + дополнительно 1...2 МБ SDRAM.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Надеюсь никто из жителей форума не будет помогать молодому и раннему любителю пива, что так любит дерзить старшим, строить полиграф для каких то явно весьма сомнительных целей.

Вы молодой человек не с Украины случаем ?

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

4 минуты назад, a123-flex сказал:

Надеюсь никто из жителей форума не будет помогать молодому и раннему любителю пива, что так любит дерзить старшим, строить полиграф для каких то явно весьма сомнительных целей.

Да ничего он не построит, не беспокойтесь. Разве что очередную свистоперделку-мигалку лампочками.  :pioneer_smoke:

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

1 час назад, jcxz сказал:

Да ничего он не построит, не беспокойтесь. Разве что очередную свистоперделку-мигалку лампочками.  :pioneer_smoke:

Я не беспокоюсь. За нас обо всем побеспокоился тот, кто создал этот прекрасный мир, и сила действия ВСЕГДА равна силе противодействия.

Но во избежание ненужных движений предлагаю неправильные движения не поддерживать. Ни на шаг.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

8 минут назад, a123-flex сказал:

Надеюсь никто из жителей форума не будет помогать молодому и раннему любителю пива, что так любит дерзить старшим, строить полиграф для каких то явно весьма сомнительных целей.

Вы молодой человек не с Украины случаем ?

Во порадовал, я из пиндоского Трампостана. Так будет лучше? Сейчас пожалуюсь на сегрегацию по национальному признаку в ООН

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

2 минуты назад, IhorU сказал:

Во порадовал, я из пиндоского Трампостана. Так будет лучше? Сейчас пожалуюсь на сегрегацию по национальному признаку в ООН

Вперед и с песней

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

53 минуты назад, Xenia сказал:

А какого рода вы хотите получить поддержку? Давайте я вас поддержу - женщины более жалостливые :).

 

Большое спасибо, вы уже секономили мне кучу времени. Меня всегда любили женщины преподаватели. Судьба ласкает молодых и ряннных. ляля ляля...

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

8 hours ago, IhorU said:

Хоть по прикалуюсь, если не получу поддержку

Только это... не сильно увлекайтесь. Иначе настанет некий возрастной порог, ну, например в 30 или 40 лет, когда поймёте, что всю жизнь потратили на "прикалюхи", не заработали ничего. А времени уже не так много, ибо семья, дом, ипотека и т.п.

А по теме: берите любую книгу по ЦОС, и изучайте. Без некой минимальной базы вы ничего не сделаете. За вас сделают другие.

Вот весьма приятная книга. Не смотрите, что написана в стиле комикса. Это у японцев подход такой. А информация там действительно полезная для первого знакомства.

8 hours ago, IhorU said:

Я уже говорил что пишу с ошибками. Так я еще не говорил что в детстве часто дрался и особенно любил бить снобов.

Сколько вам лет, милейший? Не очень много, или уже достаточно, но вы ещё в детстве счастливом? Нашли, чем гордиться.

8 hours ago, IhorU said:

А больше мне и не надо было, разберусь.

Мечты, мечты...

8 hours ago, jcxz said:

А вы говорите абдурина.

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

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

9 hours ago, Xenia said:

Выбираете массив длиной в целую степень двойки....

Заполняете его через АЦП измерениями.

Вызываете процедуру FFT для преобразования в частотную область.....

Потом вырезаете из частного спектра зоны .... умножая спектр на колоколообразную функцию с вершиной над точкой выделяемой частоты. Лично я  предпочитаю "окно Ханна" ....

 А дальше глядите на то, что получилось, и наслаждаетесь :).

Бог с ним, с любителем пива.

Имею я вопрос к Ксении.

Я не слишком большой специалист в ЦОС. Мои занятия этой частью математики ограничены моим классом решаемых задач.

Но мне показалось странным использование оконной функции в частотной области. Мне всегда представлялось логичным, что с оконной функцией сворачивают ИСХОДНЫЙ сигнал ВО ВРЕМЕННОЙ области. И делают это для того, чтобы уменьшить краевые эффекты прямоугольного окна, которые и вызывают боковые лепестки в АЧХ фильтров.

Так в чем тут прикол, Ксения?

PS. Собственно ранее я написал нашему пассионарию от пива, что нет никакой необходимости делать полное частотное преобразования сигнала, если нужно выделить лишь несколько частотных каналов. То есть достаточно определить длину буфера через требуемую полосу пропускания каналов, заполнять этот буфер в режиме скользящего окна, умножать это скользящее окно на оконную функция выбранную по критерию формы АЧХ, выполнять ДПФ (а не БПФ) только в тех частотных каналах, которые подлежат выделению. 

Если частоты каналов низкие, то реальное время не потребует высокой производительности от вычислителя. И тут будут совершенно ни к чему труДСП процессоры уровня "тигровых акул".

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

10 часов назад, jcxz сказал:

Годах этак в 2013-2014 делал я проект кардиографа, с анализом снимаемых данных. Было там множество и фильтров и БПФ. И отдельно - анализаторы этих данных,

Вот не выдержал. Вот зачем домысливать. Я никаких таких задач даже в помине не ставил. Речь идет о миниатюрном изделии которое выводит на экран сигналы отражающие эмоциональное состояние пациента. Не потянет ардуина, так много других МК. Расбери пока трогать не хочется. Вы наверно даже и не слышали о таком микрокомпютере, а там уже выросло целое поколение разных компьютеров помещающихся на долони. А Вам ардуна все жить не дает.

Зачем это все? Задачу поставили психологи а не врачи. Человек смотрит на экран и сам или под руководством тренера учится управлять своими эмоциями. Об обратной связи здесь должны знать. Это давно практикуют прогрессивные психологи при "психоанализе", не будем в суе вспоминать замшелых плешивых психоаналитиков, которые годами анализируют своих клиентов по Фройду и тянут бабку с сурйозным выражением на лице.

Геймеры тренируют свои мозги на нейроинтерфейсах для своих схваток на деньги. Это брейнфитнес называется. 

Инвалидов учат управлять своими протезами. Поклонники альцгеймера и паркинсона испражняют свои мозги до полного изнеможения в надежде не впасть в маразм

Китайцы и Запад уже нагородил кучу софта для тренировок на нейроинтерфейсах для коррекции поведения и терапии.

Что еще о нейроинтерфейсах..., а. Особо параноидальным не читать,  кондрашка она такая. Только тихо... набрали воздух... есть и военное применение. Еще тише... пиндосы на нейроинтерфейсах тренируют операторов боевых дронов, и знаете что.... Только это никому, ни слова. Тсс...

Короче очень много применения, и если ничего больше не читать кроме матанализа Смирнова или что вы там читаете- то и живите с миром.

Времена господа, меняются.

Выше переведенный текст не относится к прогрессивному человечеству и прошу не искать совпадений с реальными событиями.

Все, набрался злости, иду на тренировку а дальше и по пивку. Колотитесь 

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

1 hour ago, IhorU said:

Вы наверно даже и не слышали о таком микрокомпютере или что вы там читаете- то и живите с миром.

Времена господа, меняются.

Вы феерично смешны в своей первобытной наивности, пытаясь рассказать об ушастом запорожце людям, которые профессионально занимаются разработкой бентли...

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

16 часов назад, Xenia сказал:

Скажем, мочу по пробиркам разливают, лекарства больным дозируют :).

То есть эти устройства не должны быть critical safe и проходить соответствующие сертификации? Мне казалось, что сама элементная база, что написание кода в стиле Arduino не соответствуют данным требованиям.

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Присоединяйтесь к обсуждению

Вы можете написать сейчас и зарегистрироваться позже. Если у вас есть аккаунт, авторизуйтесь, чтобы опубликовать от имени своего аккаунта.

Гость
Ответить в этой теме...

×   Вставлено с форматированием.   Вставить как обычный текст

  Разрешено использовать не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отображать как обычную ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставлять изображения напрямую. Загружайте или вставляйте изображения по ссылке.

×
×
  • Создать...