А в чем собственно говоря проблема? Берете и строите Паде-аппроксимацию (в любой матоболочке типа Maple) для функции
в точке
, где
- вещественное положительное число (выбирается в зависимости от того, в какой частотной области вам необходимо обеспечить хорошую аппроксимацию). При этом обращаете внимание на то, чтобы полюсы были устойчивыми. Экспресс-тест показывает, что для данной функции только Паде-аппроксимации с одинаковыми степенями числителя и знаменателя обладают устойчивыми полюсами. Далее добавляете фильтр Баттерворта достаточно высокого порядка, чтобы обрубить высокие частоты в нужном месте (можно поиграть частотой среза, чтоб в итоге получить в точности то, что вам надо).