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

[:-)]

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

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

[Singer]

В случае разумных коэффициентов изменения частоты дискретизации (типа 125/256) использую всегда полифазный FIR фильтр. Просто и доступно, легко ложится на архитектуру ПЛИС (DSP48 блоки), а самое главное позволяет легко добиться нужной АЧХ и подавления заворотов спектра .Например, можно сразу получить АЧХ корня из приподнятого косинуса с нужным Roll-off, что требуется в цифровой свзяи...В этом случае всякие чисто математические Лагранжи неуместны, поскольку искажают характеристику, (действительно что называется "не подходят для сигнальной обработки"). Единственная проблема - обычно получается табличка на несколько тысяч коэффициентов, т.е. нужно памяти заметно побольше чем на простой дециматор скажем. Но по операциям MAC это все равно эквивалентно обычному децимирующему фильтру, т.к. операции со вставленными нулями пропускаются, плюс отбрасываемые сэмплы не вычисляются. При желаний конечно можно и нули на гигагерцах перемножать, но к практике это отношения не имеет...
З.Ы. - Итого мое мнение - максимальная корректность операции - полноценнный FIR фильтр, все остальное - от лукавого.
З.Ы.З.Ы. Кстати, кто вам сказал, что за пределами полосы 2.046 в GPS ничего нет? Там же уши от функции sin(x)/x -идут с периодом 2.046 (см спектр прямоуголного сигнала). В GPS передатчиках узких фильтров не ставят - там еще P код надо передать, полоса там мегагерц 20 минимум. Если не применять хороший фильтр, вырезающий основные 2 мгц, то все эти уши при децимации благополучно завернутся в основную полосу.