[an error occurred while processing this directive]
Специалисты подскажите, пожалуйста, как строить спектр сигнала, если его частота попадает между “палок”, т.е. если частоты в спектре размещены через 10 Гц, а частота сигнала, например 1005 Гц.
(«Телесистемы»: Конференция «Цифровые сигнальные процессоры (DSP) и их применение»)
|
Отправлено
msn 29 июля 2003 г. 19:36
|
|
|
|
Если кто сталкивался с такой задачей, подскажите, как правильно построить спектр. Например, если есть двух частотный сигнал 1*sin(2*pi*900*t)+1*sin(2*pi*1005*t) то после БПФ N=4096, например с Fдискретизации = 40960, частотные составляющие будут идти через 10 Гц, тогда в результате получается, две “палки” на 900 и 1000-1010 Гц, но амплитуда второй почти в полтора раза меньше, чем первой (“размыта” возле частоты 1005), но в сигнале их амплитуды одинаковые. Существует ли способ сделать так чтобы соотношение амплитуд составляющих сигнала в частотной и временной области было одинаково, при этом, не увеличивая количество точек БПФ и не уменьшая частоту дискретизации.
Спасибо.
Составить ответ
|||
Конференция
|||
Архив
Ответы
- Не очень понятен вопрос? Цель тут измерение частоты/уровня/фазы или потом потребуется восстановление во временной области? (+) — SM (29.07.2003 20:39, 131 байт)
- А он хочет, чтобы меня плавно частоту сигнала, например от 995 Гц до 1005 Гц, амплитуда спектральной составляющей не менялась. — msn (29.07.2003 21:24, пустое)
- Тогда результаты FFT проинтерполируйте например сплайном кубическим, и постройте на ЖКИ график функций-интерполянтов. Получится то, что надо Вашему заказчику :) — SM (29.07.2003 23:35, пустое)
- Если не очень накладно, могли бы Вы более подробно объяснить, как это сделать, или посоветуйте, пожалуйста, литературу или ссылки где можно почитать. — msn (29.07.2003 23:50, пустое)
- Ответ (+) — SM (30.07.2003 00:35, 2463 байт, ссылка)
- Еще есть интересный и простой по вычислительным затратам вариант - сплайн Катмулла-Рома. Широко используется при интерполяции при обработки изображений. — SM (30.07.2003 00:38, пустое)
- А вот и ссылочка на кое-какую скромную информацию про этого Катмулла с ромом... — SM (30.07.2003 00:48, пустое, ссылка)
- Огромное спасибо!!! Я тут давича почитал, так для решения этой задачи хорошо подходит наложение на оцифрованный сигнал перед FFT сглаживающей оконной функции, но сплайновую интерполяцию обязательно попробую. — msn (30.07.2003 01:20, пустое)
- Э-э - мухи отдельно, а котлеты отдельно. Окно нужно для другого! — SM (30.07.2003 07:50, пустое)
- Может как раз окно с широким плоским главным лепестком и проще для данной задачи? — Отделение 13 (30.07.2003 11:01, пустое)
- В данной задаче (нарисовать график на ЖКИ) окно имеет преимущество в скорости вычислений (если само окно посчитано заранее), а погрешность в амплитуде скроет ЖКИ :) — Отделение 13 (31.07.2003 11:41, пустое)
- Ответ: Окно могло бы быть проще, если бы не было нужно разрешение. Вроде упоминались 900 и 1005 Гц одновременно. К тому же, окно уменьшит модуляцию спектра, но не уберет ее полностью. Вместо полутора раз можно получить где-то 1.1 раза. Например 92 дб Блэкмана-Харриса. — ВН (30.07.2003 21:15, пустое)
- Какое окно не применяй, а частоту оно не сдвинет. И как она попадала "меж палок", так и попадет. — SM (30.07.2003 20:57, пустое)
- Это точно, но, применяя окно Блэкмана-Харриса получается, что “провалы” амплитуды не превышают -0,88 дБ, например при 995 Гц – 8300 единиц, а при 997,5 Гц – 7550 единиц но на 1000 Гц. Вроде бы хорошо. А если применять интерполяцию то, например при частоте сигнала 997,5 Гц на 995 получается около 6500, на 1000 – 7900, а на 1005 – 6600, то есть после интерполяции так и не получится 8300 на 997,5 Гц. — msn (31.07.2003 00:05, пустое)
- Цель – удовлетворить заказчика -)) и построить спектр сигнала на ЖКИ. — msn (29.07.2003 21:23, пустое)
Перейти к списку ответов
|||
Конференция
|||
Архив
|||
Главная страница
|||
Содержание
|||
Без кадра
E-mail:
info@telesys.ru