| |
 Интерполяция сигнала, вопрос... |
|
|
|
|
 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, 11:33
|
Знающий
Группа: Участник
Сообщений: 781
Регистрация: 3-08-09
Пользователь №: 51 730
|
Цитата Хочу минимизировать эфект "растекания" спектра Такого эффекта на сегодняшний день наука не знает... Цитата Известен интервал изменения основной гармоники (50+ 1Гц)
Частота дискретизации задана жёстко 16000Гц.
По расчётам кол-во точек сигнала равно N = 8192.
Затем по некоторому алгоритму получаем частоту основной гармоники (пусть частота равна fосн. = 50.55Гц).
Далее, по идее, для минимизации "растекания" нужно интерполировать сигнал с учётом полученной частоты основной гармоники.
Расчётные данные показали, что нужно пересчитать на частоту fd = (327.68*fосн) = 16564,224 Гц. угу. Вы хотите сделать передискретизацию так, чтоб гармоника была кратна бину дпф. Цитата Чтобы идеализировать интерполяцию просто изначально задаю частоту 16564,224 Гц на том же интервале времени, который равен интервалу времени выборки 8192 отчётов на частоте 16000Гц (т.е. (1/16000)*8192 = 0,512 с.).
Другими словами окно временное оставляю тем же, что и для 16000Гц.
Конечно же в спектре сигнала наблюдается сильное искажение...
Вот такой вопрос:
Я сделал вывод что при интерполяции каким-то образом меняется величина временного окна... иначе никак(( Гы. Вы повышаете частоту дискретизации. На интервале 0.512 c получится больше отсчетов. В Вашем случае их станет 8480,882688. Никакой метафизики тут нету... Цитата КАК ЭТО ДЕЛАЕТСЯ???? Что именно? Как вычислить дпф на 8к точек по 8480,882688-и точкам? Дык вариантов много. Можно вычислить несколько дпф на 8к точек со смещением по выборке из 8480, а потом их усреднить. Есть еще более сложные способы. Искажения спектра скорее всего связаны с горбатой интерполяцией.
|
|
|
|
|
|
|
|
Dec 1 2009, 14:12
|
Знающий
Группа: Свой
Сообщений: 688
Регистрация: 4-09-09
Пользователь №: 52 195
|
Цитата(mdmitry @ Dec 1 2009, 16:54)  Знает, в книге А.Б. Сергиенко "Цифровая обработка сигналов" посмотрите.
По теме: может, использовать общепринятые окна или метод MUSIC? всем кто смотрит тему: "пожалуйста посмотрите файл с вырезкой из диссертации в моём посте выше" Стандартное окно (Хемминга и др...) не применимо из-за погрешности порядка процентов. Здесь строится прецизионная система. Нужно использовать именно интерполяцию (предположительно фильтром Фарроу 3-го порядка)
|
|
|
|
|
|
|
|
Dec 1 2009, 14:44
|
Знающий
Группа: Свой
Сообщений: 521
Регистрация: 12-05-06
Пользователь №: 17 030
|
Цитата(TigerSHARC @ Dec 1 2009, 17:12) всем кто смотрит тему: "пожалуйста посмотрите файл с вырезкой из диссертации в моём посте выше"
Стандартное окно (Хемминга и др...) не применимо из-за погрешности порядка процентов.
Здесь строится прецизионная система.
Нужно использовать именно интерполяцию (предположительно фильтром Фарроу 3-го порядка) Где окно погрешности вносит ? В измерении амплитуды гармоник? Насчет интерполятора, Вы смотрели какая будет АЧХ у полиномиального интерполятора 3го порядка? Наверное теже несколько процентов получится завал АЧХ на верхних частотах...
--------------------
ну не художники мы...
|
|
|
|
|
|
Сообщений в этой теме
 TigerSHARC Интерполяция сигнала Dec 1 2009, 10:05 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 TigerSHARC вопрос стоит не в том, что применять весовое окно ... Dec 1 2009, 14:51 fontp Цитата(TigerSHARC @ Dec 1 2009, 17:51) Во... Dec 1 2009, 15:13 TigerSHARC Цитата(fontp @ Dec 1 2009, 18:13) Для пре... Dec 1 2009, 15:25 fontp Цитата(TigerSHARC @ Dec 1 2009, 18:25) Ра... Dec 1 2009, 15:29 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
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|
