参见§1.12连分式的相关性质,包括以下定义:
C类n个是n个第个近似的或收敛的到C类.
每个收敛、渐近或形式级数
可以转换为连分数C类类型为(3.10.1),并且具有n个第个收敛的C类n个=A类n个/B类n个到C类是等于n个中级数的第个部分和(3.10.3)也就是说,
例如,如果u个n个消失,然后我们可以定义
然而,具有相同极限的其他连续分数可能在复平面的面积大于(3.10.4)和(3.10.5). 例如,通过转换麦克劳林展开阿卡坦z(z)(4.24.3),我们获得具有相同收敛区域的连分式(|z(z)|≤1,z(z)≠±我),而连分数(4.25.4)收敛为所有人z(z)∈ℂ除了在树枝上我到我∞和−我到−我∞.
形式的连分数
称为Stieltjes分数(S公司-分数)我们这么说对应到形式幂级数
如果它的扩展n个第个收敛的C类n个以提升的力量z(z)同意(3.10.7)直至并包括术语z(z)n个−1,n个=1,2,三,….
对于几个特殊功能S公司-分数是明确的,但在在任何情况下,系数一n个始终可以从幂级数计算系数,通过商微分算法; 看见表3.10.1.
此表中的前两列由定义
其中c(c)n个(≠0)出现在(3.10.7). 我们继续的手段菱形规则
然后是系数一n个的S公司-分数(3.10.6)是由提供
商微分算法通常不稳定,可能需要高精度算法或精确算法。更稳定的版本算法在中进行了讨论斯托克斯(1980)适用于贝塞尔函数和Whittaker函数(第章10和13),请参阅Gargantini和Henrici(1967).
称为雅可比分数(J-分数)我们说是的相关用形式幂级数(f)(z(z))在里面(3.10.7)如果它的扩展n个第个收敛的C类n个在里面提升的力量z(z),同意(3.10.7)到和将该术语包含在z(z)2n个−1,n个=1,2,三,….对于同一功能(f)(z(z)),收敛的C类n个雅可比分数的(2011年3月3日)等于收敛的C类2n个Stieltjes分数(3.10.6).
有关基本函数,请参见§§4.9和4.35.
有关特殊功能,请参见§5.10(伽马函数),§7.9(错误函数),§8.9(不完整的伽马函数),§8.17(v)(不完整的β函数),§8.19(七)(广义指数积分),§§10.10和10.33(贝塞尔函数的商),§13.6(合流超几何函数的商),§13.19(Whittaker函数的商),以及§15.7(超几何函数的商)。
有关更多信息和示例,请参见洛伦岑和瓦德兰(1992第292-330、560-599页)和库伊特等。(2008).
这个A类n个和B类n个第页,共页(3.10.2)可以通过以下方式计算三项递推关系(1.12.5). 然而,这可能是不稳定;评估时也可能发生溢出和下溢A类n个和B类n个(需要不时重新缩放)。
要计算C类n个第页,共页(3.10.2)我们执行迭代除法
然后u个0=C类n个.为了达到规定的精确度先验的知识的价值是需要的n个,或n个通过试验确定错误。一般来说,该算法比前向算法更稳定;看见琼斯和王位(1974).
连续分数
可以写在表格中
哪里
这个n个第次部分和t吨0+t吨1+⋯+t吨n个−1等于n个第个收敛的第页,共页(3.10.13),n个=1,2,三,…与前一项相比本小节中的算法不会出现缩放问题,也不会先验的需要信息。
在高茨基(1979年)正向序列算法用于评估不完整伽马函数的连续分数(参见§8.9).
该正向算法在计算收敛点C类n个=A类n个/B类n个,与正向系列相关递归算法。同样,不会出现缩放问题,也不会先验的需要信息。
让
(∇是后向差分算子.)然后针对n个≥2,
复发持续到(∇C类n个)/C类n个在规定范围内相对精度。
Steed算法的替代方案是伦茨算法 伦茨(1976)和改进的Lentz算法 汤普森和巴奈特(1986).
有关上述算法的更多信息,包括复数平面和加速收敛的方法,参见布兰奇(1964)和洛伦岑和瓦德兰(1992,第3章)。对于用连分式计算特殊函数库伊特等。(2008),高奇(1967, §1),吉尔等。(2007年a,第6章)、和Wimp公司(1984,第4章,§5)。请参阅也包括§§6.18(i),7.22(i),8.25(iv),10.74(v),14.32,28.34(ii),29.20(i),30.16(i),33.23(v).