Главная » 2011 » Сентябрь » 23 » Ресемплинг - моделирование
20:27 Ресемплинг - моделирование |
Чтобы так сказать расставить точки над и решил смоделировать два подхода к ресемплингу (линейная интерполяция и заполнение новых отсчетов нулями) и сравнить как они влияют на частотный спектр сигнала. Для минимизации времени вычислений взял исходный сигнал в частотном диапазоне - 0 - 1024 Гц при частоте семплирования 2048 Гц - это рабочий диапазон частот. Использовался вектор в 2048 отсчетов. Разрешающая способность частотного спектра при этом получается 1 Гц. Соответсвенно ресемплинг делал в 4096 Гц. Для того чтобы оценить как оба метода модифицируют частотный спектр использовал два сигнала -
Исходный сигнал синтезировал из предварительно созданного частотного спектра. Для варианта а) все значения частот, кроме заданных равны нулю, для варианта б) все амплитудные значения имеют одну величину, при этом значения фазы заполнялось случайно в диапазоне от -Пи до +Пи. Далее к полученному спектру применялось обратное преобразование Фурье, на выходе соответственно имеем зависимость амплитуды от времени. Далее обработка состояла из шагов:
Вычитание спектров означает что из амплитудного значения каждой частоты ресемпленного сигнала вычиталось амплитудное значение той же частоты исходного сигнала. Также нужно понимать что ресемпленный сигнал имеет спектр в два раза шире базового, тоесть потенциально может содержать частоты до 2048 Гц в моем случае. Соответственно вычитание производится только в той области ресемпленого спектра где есть общие частоты. область 1025-2048 Гц модифицированного сигнала показана без изменений. Результат показан на графиках ниже. Вариант A: Я увеличил частоту 950 Гц для сравнения: Здесь видно, что:
Вот. Теперь факты
|
|
Всего комментариев: 6 | |||||||
| |||||||