Как приблизить функцию Теодорсена, используя ряды Паде или Тейлора?
У меня есть следующее выражение для функции Теодорсена:
C(k)= (besselh(1,2,k))./(besselh(1,2,k)+1i*besselh(0,2,k));
где k это уменьшенная частота
Я хотел бы найти рациональное дробное приближение к функции Теодорсена с полюсами и нулями. я пробовал padeFNC, но я не получаю то, что я ищу (как найдено в этой работе):
В приведенной выше ссылке они применили метод наименьших квадратов, получив следующее приближение 2-го порядка (2 полюса - 2 нуля):
C(k) =0.5(ik+0.135)(ik+0.651)/(ik+0.0965)(ik+0.4555)
где s = ik
Я пытаюсь придумать приближение Тейлора или Паде к функции C(k) Теодорсена без использования метода наименьших квадратов, поскольку этот подход уменьшает ошибку во всем частотном диапазоне. Итак, я хотел бы сделать, как ряд Тейлора, чтобы он действительно хорошо подходил к более низким частотам, а затем попытаться построить на более высоких частотах (на графике Боде).