адаптивный ресеплинг, нужна помощь Опции

[TigerSHARC]

Требуется построить адаптивную систему изменения частоты дискретизации сигнала.

Сначала в сторону фильтра Фарроу смотрел, там как я понял в основе лежит апроксимация полиномами Лагранжа.
Но тут вдруг недавно наткнулся на научную работу, где прочитал следующие строки:
"Реализация формулы Лагранжа в вычислительном контексте неудобна, так как произведение n сомножителей приводит с существенным погрешностям при округлениях в разрядной сетке вычислительного устройства. Поэтому на практике применяют интерполяционные формулы, в основе которых лежит идея конечных разностей."

Далее приводится вывод о том, что наиболее подходит интерполяция Стирлинга 4-го порядка...

Где можно найти пример реализации такого подхода?
И какие ещё существуют методы изменения частоты дискретизации?

[EvgenyNik]

В вашей задаче (если речь про Фурье всех гармоник) как раз Фарроу 3-го порядка и достаточен и вот почему:
1. Частота дискретизации, по-любому, из соображений точности должна накрывать сигнал так, чтобы на интервале между соседними 3-4 выборками сигнал мало отличался от 3-го порядка, т.к. иначе Вы с достаточной степенью точности этот сигнал не отфильтруете из-за вносимых искажений от гармоник более высоких порядков.
2. Т.к. исходный сигнал, всё-таки, содержит небольшой шум АЦП, к включению которого в огибающую интерполятор будет стремиться в меру своей степени, то наращивание порядка интерполяции приведёт к генерации несуществующих компонент.