Сдвиг частот для подавления акустической обратной связи

@4eaeNeoStudy · Регистрация: 14.05.2013

Author24 — интернет-сервис помощи студентам

Есть звук - делаю БФП, спектр сдвигаю на несколько Герц вправо, делаю обратное БФП - получается какой-то шум резкий. Подскажите пожалуйста, как правильно сдвинуть спектр. Звук 32-битный, 44100Гц, БФП делается с помощью fftw

raxper · 14.05.2013, 16:29

...а ежели не сдвигать, после IFFT звук оки? Вероятно где-то ошиблись, где, не видно.

@4eaeNeoStudy · 15.05.2013, 10:24 **[ТС]**

Сообщение от raxp

...а ежели не сдвигать, после IFFT звук оки? Вероятно где-то ошиблись, где, не видно.

если не сдвигать, то нормально. Вообще, звук обрабатывается порциями по 512 семплов, это может влиять?

raxper · 15.05.2013, 12:12

...вот потому и говорю, не видно как. Вы что же, FFT сразу над всеми сэмплами делаете или поочередно все-таки?

...под сэмплом я подразумеваю - выборку за один раз, набор отсчетов в заданном буфере.

@vital792 · 15.05.2013, 13:51

Сообщение от 4eaeNeoStudy

делаю БФП, спектр сдвигаю на несколько Герц вправо, делаю обратное БФП

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

raxper · 15.05.2013, 13:58

...могу предположить: сдвиг на пару герц используется для подавления обратной связи для акустических систем при возбуждении. Ранее на сельсинах делали.

@4eaeNeoStudy · 15.05.2013, 15:13 **[ТС]**

Сообщение от raxp

...могу предположить: сдвиг на пару герц используется для подавления обратной связи для акустических систем при возбуждении. Ранее на сельсинах делали.

да, сдвиг для подавления обратной связи. Я под семплом имел виду, то что вы называете отсчётом. С микрофона с частотой 44100, поступает звук буфером 512 отсчётов (32 битный звук), я делаю БФП - получаю спектр на 512 частот, последняя половина нули. Расстояние между двумя соседнями частотами получается (44100/512 ~ 86 Гц) Чтобы сдвинуть на примерно 5 Гц я добавляю между соседнями частотами 15 нулей, получем спектр на 512 * 16 = 8192 частот, растояние между которыми чуть больше 5 Гц. Затем весь спектр, кроме нулевой частоты, сдвигаю на один отсчет вправо и делаю обратное БПФ, получаю 8192 отсчётов, из которых беру первые 512. Вроде правильно рассуждаю - просто, вроде бы. Получается совсем не очень)

raxper · 15.05.2013, 16:02

...вот в конце ошибка: изначально в исходном сигнале 512 отсчетов взято для заданной fd, на выходе имеем буфер отсчетов того же сигнала, но нарезанного как бы чаще (но на той же длине), а вы вдруг выкидываете большую часть восстановленного сигнала. Для записи с тем же размером буфера сырых (но уже сдвинутых по частоте) отсчетов следует сделать прореживание, т.е. взять каждый 17-й из набора 8820 (у меня для нарезки по 5 герц получилось).

@4eaeNeoStudy · 15.05.2013, 16:29 **[ТС]**

что-то не очень понимаю. Когда в спектр добавляю между каждыми частотами нули - это эквивалетно тому, что я взял 512 отсчётов и повторил их 15 раз, итого получил 512*16 отсчётов с той же частотой 44100. Поэтому и при обратном нужно взять 512 первых - частота-та как раз не меняется - или я что-то не понимаю? Или вообще не так делаю?

raxper · 15.05.2013, 16:47

...после обратного БПФ вы работаете с огибающей сигнала, вот это то вы и забываете: когда вы добавляете виртуальные гармоники (нули) в исходный спектр - вы имитируете взятие куска огибающей с большим количеством отсчетов, но в конце это та же огибающая на том же временном участке, только нарезана чаще. Просто возьмите листочек бумаги и прорисуйте все пошагово. Ферштеен?

@4eaeNeoStudy · 16.05.2013, 10:01 **[ТС]**

Начинаю понимать, спасибо большое. Буду разбираться

Новые блоги и статьи Все статьи Все блоги /
Kubernetes с Apache Flink для обработки данных в реальном времени Mr. Docker 17.05.2025 Kubernetes — это целая философия управления распределёнными приложениями. В отличие от "примитивных" решений вроде Docker Swarm, K8s (как его ласково называют в тусовке DevOps-инженеров) предлагает. . .	Использование декораторов в Python py-thonny 17.05.2025 Если вы когда-нибудь задумывались о том, как красиво расширить функциональность кода без лишней возни и дублирования, декораторы в Python — та самая волшебная палочка, которую вы искали. По сути, это. . .	Как некоторые отстают в навыках поиска информации в сети. Дублирование страницы открытого адреса в сети. Hrethgir 17.05.2025 . . . берём бразуер, для многих он является инструментом поиска информации в сети. Понятное дело, что чтобы переходить по сслыкам есть такая фишка - "открыть в новой вкладке", но, представим себе что вы. . .	Реализация многопоточных сетевых серверов на Python py-thonny 16.05.2025 Когда сталкиваешься с необходимостью писать высоконагруженные сетевые сервисы, выбор технологии имеет критическое значение. Python, со своей элегантностью и высоким уровнем абстракции, может. . .	C# и IoT: разработка Edge приложений с .NET и Azure IoT UnmanagedCoder 16.05.2025 Мир меняется прямо на наших глазах, и интернет вещей (IoT) — один из главных катализаторов этих перемен. Если всего десять лет назад концепция "умных" устройств вызывала скептические улыбки, то. . .
Гибридные квантово-классические вычисления: Примеры оптимизации EggHead 16.05.2025 Гибридные квантово-классические вычисления — это настоящий прорыв в подходах к решению сложнейших вычислительных задач. Представьте себе союз двух разных миров: классические компьютеры, с их. . .	Использование вебсокетов в приложениях Java с Netty Javaican 16.05.2025 HTTP, краеугольный камень интернета, изначально был спроектирован для передачи гипертекста с минимальной интерактивностью. Его главный недостаток в контексте современных приложений — это. . .	Реализация операторов Kubernetes Mr. Docker 16.05.2025 Концепция операторов Kubernetes зародилась в недрах компании CoreOS (позже купленной Red Hat), когда команда инженеров искала способ автоматизировать управление распределёнными базами данных в. . .	Отражение в C# и динамическое управление типами stackOverflow 16.05.2025 Reflection API в . NET — это набор классов и интерфейсов в пространстве имён System. Reflection, который позволяет исследовать и манипулировать типами, методами, свойствами и другими элементами. . .	Настройка гиперпараметров с помощью Grid Search и Random Search в Python AI_Generated 15.05.2025 В машинном обучении существует фундаментальное разделение между параметрами и гиперпараметрами моделей. Если параметры – это те величины, которые алгоритм "изучает" непосредственно из данных (веса. . .

@4eaeNeoStudy 0 / 0 / 0 Регистрация: 14.05.2013 Сообщений: 5

	Сдвиг частот для подавления акустической обратной связи 14.05.2013, 15:41. Показов 9979. Ответов 10 Метки нет (Все метки) Есть звук - делаю БФП, спектр сдвигаю на несколько Герц вправо, делаю обратное БФП - получается какой-то шум резкий. Подскажите пожалуйста, как правильно сдвинуть спектр. Звук 32-битный, 44100Гц, БФП делается с помощью fftw 0

raxper 10236 / 6614 / 498 Регистрация: 28.12.2010 Сообщений: 21,154 Записей в блоге: 1
	14.05.2013, 16:29
	...а ежели не сдвигать, после IFFT звук оки? Вероятно где-то ошиблись, где, не видно. 0

raxper 10236 / 6614 / 498 Регистрация: 28.12.2010 Сообщений: 21,154 Записей в блоге: 1
	15.05.2013, 12:12
	...вот потому и говорю, не видно как. Вы что же, FFT сразу над всеми сэмплами делаете или поочередно все-таки? ...под сэмплом я подразумеваю - выборку за один раз, набор отсчетов в заданном буфере. 0

raxper 10236 / 6614 / 498 Регистрация: 28.12.2010 Сообщений: 21,154 Записей в блоге: 1
	15.05.2013, 13:58
	...могу предположить: сдвиг на пару герц используется для подавления обратной связи для акустических систем при возбуждении. Ранее на сельсинах делали. 1

raxper 10236 / 6614 / 498 Регистрация: 28.12.2010 Сообщений: 21,154 Записей в блоге: 1
	15.05.2013, 16:02
	...вот в конце ошибка: изначально в исходном сигнале 512 отсчетов взято для заданной fd, на выходе имеем буфер отсчетов того же сигнала, но нарезанного как бы чаще (но на той же длине), а вы вдруг выкидываете большую часть восстановленного сигнала. Для записи с тем же размером буфера сырых (но уже сдвинутых по частоте) отсчетов следует сделать прореживание, т.е. взять каждый 17-й из набора 8820 (у меня для нарезки по 5 герц получилось). 0

@4eaeNeoStudy 0 / 0 / 0 Регистрация: 14.05.2013 Сообщений: 5
	15.05.2013, 16:29 [ТС]
	что-то не очень понимаю. Когда в спектр добавляю между каждыми частотами нули - это эквивалетно тому, что я взял 512 отсчётов и повторил их 15 раз, итого получил 512*16 отсчётов с той же частотой 44100. Поэтому и при обратном нужно взять 512 первых - частота-та как раз не меняется - или я что-то не понимаю? Или вообще не так делаю? 0

raxper 10236 / 6614 / 498 Регистрация: 28.12.2010 Сообщений: 21,154 Записей в блоге: 1
	15.05.2013, 16:47
	...после обратного БПФ вы работаете с огибающей сигнала, вот это то вы и забываете: когда вы добавляете виртуальные гармоники (нули) в исходный спектр - вы имитируете взятие куска огибающей с большим количеством отсчетов, но в конце это та же огибающая на том же временном участке, только нарезана чаще. Просто возьмите листочек бумаги и прорисуйте все пошагово. Ферштеен? 2

@4eaeNeoStudy 0 / 0 / 0 Регистрация: 14.05.2013 Сообщений: 5
	16.05.2013, 10:01 [ТС]
	Начинаю понимать, спасибо большое. Буду разбираться 0

Опции темы