Вопрос про децимацию Опции

[:-)]

Возникла вот какая задача. Требуется перейти от частоты дискретизации 4 МГц к частоте дискретизации 2.048МГц. Сейчас для выполнения децимации просто выкидываю лишние данные (беру только i*2048/4000 1<i<=4000). Но возник вопрос насколько это коректно?

Каким алгоритмом следует воспользоваться? (Два главных критерия: максимальная "корректоность" операции и максимальная скорость, т.к. требуется выполнять эту операцию параллельно с другими в реальном времени на ПК).

[анатолий]

Вы наверно хотели сказать что сигнал в самых младших битах и естественный ДД очень маленький (т.к. еще от 30(1 период) до 43(20 периодов) дб сверху вытягивает ПСП), поэтому подавление паразитов от некратной децимации не обязательно делать сильным (типа 80 дб)? В такой постановке я соглашусь...

Совершенно верно.
Сигнал сильно замешан в шумах - его уровень может быть в десятки раз ниже шумов, т.е. он может быть "закопан" в случайности ошибок округления.
По этому случаю в таких приемниках, напр, опорная синусоида гетеродина - всего 2-4 разряда,
а последующий ФНЧ - 6-8 разрядный.
Просто, вся информация потом собирается в корреляторе, для которого достаточно 1-4 входных правильных разряда, т.е. напр, один лишь знаковый разряд.
Но коррелятор должен быть очень синхронным, для чего частота дискретизации играет роль.

[AlikM]

Совершенно верно.
Сигнал сильно замешан в шумах - его уровень может быть в десятки раз ниже шумов, т.е. он может быть "закопан" в случайности ошибок округления.
По этому случаю в таких приемниках, напр, опорная синусоида гетеродина - всего 2-4 разряда,
а последующий ФНЧ - 6-8 разрядный.
Просто, вся информация потом собирается в корреляторе, для которого достаточно 1-4 входных правильных разряда, т.е. напр, один лишь знаковый разряд.
Но коррелятор должен быть очень синхронным, для чего частота дискретизации играет роль.

Я как-то занимался приемом и демодуляцией GPS сигналов, попытаюсь поделиться своим опытом. Сигнал GPS действительно ниже уровня теплового шума. Из документации по GPS знаем, что расчетный уровень сигнала на поверхности земли около -160 dBW. Для простоты будем считать что спектральная плотность мощности теплового шума -204 dBW/Hz. Для полосы 2МГц уровень шума -140 dBW, т.е. на 20 dB ниже уровня шумов. Это надо обязательно учитывать при согласовании приемного тракта и АЦП. При аналого-цифровом преобразовании уровень сигнала на входе АЦП должен быть выше шумов квантования. А "закапывая" сигнал под шум квантования Вы сами себе выкручиваете руки.
Тут еще разгорелся спор про интерполяцию полифазником или полиномом Логранжа. Логранж можно считать FIR фильром, а если посчитать его коэффициенты для фиксированных значений параметра(равные промежутки задержки), то получим тоже что-то вроде полифазника. А сказать какой является оптимальным можно рассматривая АЧХ FIR фильтров (полифазного и полинома Логранжа). Например, если речь идет об оптимальности с точки зрения максимального подавления за полосой и минимальной переходной полосы, то полифазник полученный по Ремезу будет оптимальным и полином Логранжа для получения тех же характеристик потребует больше вычислений. Можно предположить, что при других критериях оптимальности выигрыш будет у интерполятора Логранжа, но для меня исторически сложилось, что мне удобней полифазник. В аккумулятор задаю код требуемой на выходе частоты, остаток при переполнении аккумулятора показывает нужный номер фильтра. Поэтому, конечно же, ни о каких гигагерцах речь не идет.