Детектирование "нелинейных" звуковых чирпов в аудио потоке., чем сдетектировать? Опции

[tmtlib]

Подскажите, чем определять синусоидальный сигнал с нелинейно растущей частотой? Есть шумы и эхо. Спектрограмма БПФ по времени даёт весьма размазанную картину, так как частота растёт быстро, амплитуда падает в шум, т.е. БПФ с большим окном попросту не подходит для таких писков. В спектрограмме конечно что-то видно, но алгоритмически это выделять уже сложно - и нужна память на хранение спектрограммы. На частоте дискретизации 22.5кГц в контроллер влазит примерно 2048 отсчёта (упаковываю по два 8-битных отстчёта в одну 16-битную ячейку памяти).

[tmtlib]

ну что же, стоит признать, что мой метод оказался бредом, хотя и не полнейшим. Та маленькая картинка сверху - это гистограмма, а не спектр. Спектр выглядит намного хуже - всё в шумах.

Теперь подумываю над предложением sup-sup засунуть эталон в FIR фильтр и натравить его на децимированный сигнал.

p.s. Всем спасибо, сделать открытие в области ЦОС мне пока не удалось

У меня ещё вопрос! Я хочу выбрать эталон, с которым потом будет проводиться корреляция. Но при этом в звуковом сигнале присутствует многократное эхо, из-за чего синусоиды этого звука перемешались. Будет ли в этом случае у меня работать корреляция? Я ещё могу представить как это будет работать на "аналитическом" чирпе, нечто вроде ДПФ для одной частоты. Но здесь свёртка с такой фиговиной. Я конечно попробую, потом отпишу о результатах...

[sup-sup]

.. и натравить его на децимированный сигнал.
.. многократное эхо

Я предложил децимацию после свертки сигнала с комплексным эталоном. То есть, получить огибающую, с шагом в N выборок. Шаг децимации зависит от максимальной частоты спектра огибающей.

[tmtlib]

Я предложил децимацию после свертки сигнала с комплексным эталоном. То есть, получить огибающую, с шагом в N выборок. Шаг децимации зависит от максимальной частоты спектра огибающей.

ok. А я думал для ускорения свёртки.....

[sup-sup]

ok. А я думал для ускорения свёртки.....

И еще, нужно не забыть про комплексное дополнение для ядра фильтра - можно сделать через преобразование Гильберта.

Сообщение отредактировал sup-sup - Aug 3 2011, 11:34