| |
 Интерполяция сигнала, вопрос... |
|
|
|
|
 Dec 1 2009, 10:05
|
Знающий
Группа: Свой
Сообщений: 688
Регистрация: 4-09-09
Пользователь №: 52 195
|
Хочу минимизировать эфект "растекания" спектра. Сигнал содержит 40 гармоник кратных основной частоте.
Известен интервал изменения основной гармоники (50+ 1Гц)
Частота дискретизации задана жёстко 16000Гц.
По расчётам кол-во точек сигнала равно N = 8192.
Затем по некоторому алгоритму получаем частоту основной гармоники (пусть частота равна fосн. = 50.55Гц).
Далее, по идее, для минимизации "растекания" нужно интерполировать сигнал с учётом полученной частоты основной гармоники.
Расчётные данные показали, что нужно пересчитать на частоту fd = (327.68*fосн) = 16564,224 Гц.
________________________________________________
Вот здесь самое интересное начинается!!!
Чтобы идеализировать интерполяцию просто изначально задаю частоту 16564,224 Гц на том же интервале времени, который равен интервалу времени выборки 8192 отчётов на частоте 16000Гц (т.е. (1/16000)*8192 = 0,512 с.).
Другими словами окно временное оставляю тем же, что и для 16000Гц.
Конечно же в спектре сигнала наблюдается сильное искажение...
Вот такой вопрос:
Я сделал вывод что при интерполяции каким-то образом меняется величина временного окна... иначе никак((
КАК ЭТО ДЕЛАЕТСЯ????
Встал на этом вопросе прирасчётах, очень нужна помощь...
Заранее спасибо всем ответившим!
|
|
|
|
|
|
|
 |
Ответов
|
|
Dec 1 2009, 15:13
|
Эксперт
Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183
|
Цитата(TigerSHARC @ Dec 1 2009, 17:51)  Вопрос в том, как построить адаптивную систему, чтобы менять длину временного окна (почитайте отрывок из диссертации выше))
Я вот думаю что с помощью интерполяции... Для прецизионности системы нужна интерполяция синком А что спектр надо оценивать не только для гармоник основной частоты, но и между ними? Если нет - то зачем БПФ? Достаточно ДПФ и быстрее. Существует для любого окна N Как говорил petrov - банк фильтров на кратных частотах.
|
|
|
|
|
|
|
|
Dec 1 2009, 15:25
|
Знающий
Группа: Свой
Сообщений: 688
Регистрация: 4-09-09
Пользователь №: 52 195
|
Цитата(fontp @ Dec 1 2009, 18:13) Для прецизионности системы нужна интерполяция синком
А что спектр надо оценивать не только для гармоник основной частоты, но и между ними?
Если нет - то зачем БПФ? Достаточно ДПФ и быстрее. Существует для любого окна N
Как говорил petrov - банк фильтров на кратных частотах. Интерполяция меняет размер окна??? Разве на 40 гармоник ДПФ будет быстрее??
|
|
|
|
|
|
|
|
Dec 1 2009, 15:29
|
Эксперт
Группа: Свой
Сообщений: 1 467
Регистрация: 25-06-04
Пользователь №: 183
|
Цитата(TigerSHARC @ Dec 1 2009, 18:25) Разве на 40 гармоник ДПФ будет быстрее?? 40 гармоник ДПФ - это примерно 80 * N действительных операций (40 для синуса, 40 для косинуса MAC) FFT - это 2*N*logN = 26*N комплексных операций. Вы же вычисляете не 40, а 4096 гармоник, из которых Вас интересует, возможно, только каждая 10-я. Зависит от постановки задачи. Конечно, если Вам нужен красивый спектроанализатор тогда другое дело... Если не быстрее, то примерно одинаково Цитата(TigerSHARC @ Dec 1 2009, 18:25) Интерполяция меняет размер окна??? Какая интерполяция?  Если Вы пересчитаете свои 8192 отсчета на другую дискретную сетку, то понятно, что размер окна в отсчетах изменится. А в физических единицах - нет. Ширина полосы фильтра, который можно получить по этим точкам в герцах не изменится
|
|
|
|
|
|
Сообщений в этой теме
 TigerSHARC Интерполяция сигнала Dec 1 2009, 10:05 thermit ЦитатаХочу минимизировать эфект "растекания... Dec 1 2009, 11:33  mdmitry Цитата(thermit @ Dec 1 2009, 14:33) Таког... Dec 1 2009, 13:54 TigerSHARC Цитата(mdmitry @ Dec 1 2009, 16:54) Знает... Dec 1 2009, 14:12 alex_os Цитата(TigerSHARC @ Dec 1 2009, 17:12) вс... Dec 1 2009, 14:44 petrov ДПФ это банк фильтров с прямоугольной ИХ, АЧХ кажд... Dec 1 2009, 11:53 alex_os Цитата(TigerSHARC @ Dec 1 2009, 13:05) ..... Dec 1 2009, 12:04 petrov Цитата(alex_os @ Dec 1 2009, 15:04) А зач... Dec 1 2009, 12:09 TigerSHARC Цитата(thermit @ Dec 1 2009, 14:33) Таког... Dec 1 2009, 13:52 thermit Цитатаmdmitry:
Знает, в книге А.Б. Сергиенко ... Dec 1 2009, 14:30 mdmitry Цитата(thermit @ Dec 1 2009, 17:30) Не чи... Dec 1 2009, 15:37 TigerSHARC как по научному-то??? у Сергиенко тоже "расте... Dec 1 2009, 14:36 alex_os Цитата(TigerSHARC @ Dec 1 2009, 17:51) во... Dec 1 2009, 15:30 thermit ЦитатаTigerSHARC:
как по научному-то??? у Сергиенк... Dec 1 2009, 15:37 TigerSHARC А как объяснить такое.
Читаю:
"... С учётом б... Dec 1 2009, 15:44 fontp Цитата(TigerSHARC @ Dec 1 2009, 18:44) А ... Dec 1 2009, 16:04 TigerSHARC Цитата(fontp @ Dec 1 2009, 19:04) Наприме... Dec 1 2009, 16:14 alex_os Цитата(TigerSHARC @ Dec 1 2009, 18:44) ь.... Dec 1 2009, 16:05 thermit Цитатаmdmitry:
Укажите какие.
Рабинер Голд Теория... Dec 1 2009, 15:48 mdmitry Цитата(thermit @ Dec 1 2009, 18:48) Рабин... Dec 1 2009, 22:07 thermit ЦитатаTigerSHARC:
преобразую в сигнал с частотой 1... Dec 1 2009, 16:15 Alex11 Да не нужен Вам никакой ресемплер, наложите Гауссо... Dec 1 2009, 16:21 TigerSHARC Цитата(Alex11 @ Dec 1 2009, 19:21) Да не ... Dec 1 2009, 18:00 SPACUM Привет SHARK. По Твоей тематике можешь посмотреть ... Dec 1 2009, 21:33 TigerSHARC Цитата(SPACUM @ Dec 2 2009, 00:33) Привет... Dec 2 2009, 11:04 SPACUM Цитата(TigerSHARC @ Dec 2 2009, 14:04) Мы... Dec 2 2009, 16:55 TigerSHARC Цитата(SPACUM @ Dec 2 2009, 19:55) Прошу ... Dec 2 2009, 19:28 fontp Цитата(TigerSHARC @ Dec 2 2009, 22:28) а ... Dec 3 2009, 07:47 Tanya Цитата(fontp @ Dec 3 2009, 10:47) Там как... Dec 3 2009, 08:14 SPACUM Цитата(Tanya @ Dec 3 2009, 11:14) Вот эта... Dec 4 2009, 06:56 TigerSHARC Цитата(SPACUM @ Dec 4 2009, 09:56) Большо... Dec 4 2009, 10:13 SPACUM Цитата(TigerSHARC @ Dec 4 2009, 13:13) Го... Dec 4 2009, 13:31 TigerSHARC Цитата(SPACUM @ Dec 4 2009, 16:31) Честно... Dec 4 2009, 15:01 fontp Цитата(TigerSHARC @ Dec 4 2009, 18:01) Та... Dec 4 2009, 15:05 hobgoblin Цитата(fontp @ Dec 4 2009, 18:05) Так он ... Dec 4 2009, 15:31 fontp Цитата(hobgoblin @ Dec 4 2009, 18:31) Фил... Dec 4 2009, 15:54 SPACUM Цитата(TigerSHARC @ Dec 4 2009, 18:01) ..... Dec 4 2009, 16:07 TigerSHARC Цитата(SPACUM @ Dec 4 2009, 19:07) Прошу ... Dec 4 2009, 17:25 fontp Вот та статья VU Jiekang из IEEE более соответству... Dec 3 2009, 19:36 thermit ЦитатаTigerSHARC:
Так вот главный вопрос всё ещё а... Dec 4 2009, 15:34 hobgoblin Сам по себе фильтр фильтр фарроу это просто элемен... Dec 4 2009, 15:52 hobgoblin Это не философия, а вопрос терминологии В после... Dec 4 2009, 16:13 fontp Цитата(hobgoblin @ Dec 4 2009, 19:13) Это... Dec 4 2009, 17:42 TigerSHARC Цитата(fontp @ Dec 4 2009, 20:42) Для пос... Jan 8 2010, 21:14 hobgoblin ЦитатаДля постоянной дробной задержки я предпочту ... Dec 4 2009, 20:12 alexr22b Можно применить сдвиг частоты так, как это делаетс... Dec 4 2009, 21:47 ivan219 Ну как вареант можно ещё и здвинуть сигнал использ... Jan 9 2010, 19:00 TigerSHARC по-моему тогда сильно возрастёт погрешность... Jan 10 2010, 14:48
2 чел. читают эту тему (гостей: 2, скрытых пользователей: 0)
Пользователей: 0
|
|
|
