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

[TigerSHARC]

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

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

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

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

[des00]

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

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

а работу можно посмотреть ? как то все странно это. единственное место потери точности у фарроу это умножение на mu/D. Но даже по крайней мере до третьего порядка это можно сделать вообще без потери точности.