我对西奥多森的功能有以下表达：
C(k)= (besselh(1,2,k))./(besselh(1,2,k)+1i*besselh(0,2,k));
其中k是降低的频率
我想找到一个带有极点和零点的西奥多森函数的有理分数近似。
我已经尝试过pade fnc，但没有得到想要的东西（如此work所示）：
在上面的链接中，他们进行了最小二乘拟合，获得了以下二阶近似值（2个极点-2个零）：
C(k) =0.5(ik+0.135)(ik+0.651)/(ik+0.0965)(ik+0.4555)
s = ik
我试图在不使用最小二乘方的情况下得出泰奥多尔（Theodorsen）C（k）函数的泰勒或帕德近似值，因为这种方法减小了整个频率范围的误差。所以，我想做一个泰勒级数，使其非常适合低频 然后尝试建立更高的频率（在波德图中）。