关于项目
数值方法地区

§3.10连续分数

目录
  1. §3.10(i)介绍
  2. §3.10(ii)与幂级数的关系
  3. §3.10(iii)连分式的数值计算

§3.10(i)介绍

参见§1.12连分式的相关性质,包括以下定义:

3.10.1 C类=b条0+1b条1+2b条2+,
n个0,
3.10.2 C类n个=b条0+1b条1+2b条2+n个b条n个=A类n个B类n个.

C类n个n个第个近似的收敛的C类.

§3.10(ii)与幂级数的关系

每个收敛、渐近或形式级数

3.10.3 u个0+u个1+u个2+

可以转换为连分数C类类型为(3.10.1),并且具有n个第个收敛的C类n个=A类n个/B类n个C类等于n个中级数的第个部分和(3.10.3)也就是说,

3.10.4 A类n个B类n个=u个0+u个1++u个n个,
n个=0,1,.

例如,如果u个n个消失,然后我们可以定义

3.10.5 b条0 =u个0,
b条1 =1,
1 =u个1,
b条n个 =1+u个n个u个n个1,
n个 =u个n个u个n个1,
n个2.

然而,具有相同极限的其他连续分数可能在复平面的面积大于(3.10.4)和(3.10.5). 例如,通过转换麦克劳林展开阿卡坦z(z)(4.24.3),我们获得具有相同收敛区域的连分式(|z(z)|1,z(z)±),而连分数(4.25.4)收敛为所有人z(z)除了在树枝上.

Stieltjes分数

形式的连分数

3.10.6 C类=011z(z)12z(z)1

称为Stieltjes分数(S公司-分数)我们这么说对应到形式幂级数

3.10.7 (f)(z(z))=c(c)0+c(c)1z(z)+c(c)2z(z)2+

如果它的扩展n个第个收敛的C类n个以提升的力量z(z)同意(3.10.7)直至并包括术语z(z)n个1,n个=1,2,,.

商差算法

对于几个特殊功能S公司-分数是明确的,但在在任何情况下,系数n个始终可以从幂级数计算系数,通过商微分算法; 看见3.10.1.

表3.10.1:商差方案。
\开始{picture}(8.0,9.0)\put(0.0,7.0){$e_{0}^{1}$}\投入(0.0,5.0){$e_{0}^{2}$}\put(0.0,3.0){$e_{0}^{3}$}\put(0.0,1.0){$e_{0}^{4}$}\put(1.0,8.0){$q{1}^{0}$}\put(1.0,6.0){$q{1}^{1}$}\put(1.0,4.0){$q{1}^{2}$}\put(1.0,2.0){$q{1}^{3}$}\put(1.0,0.0){$\ddots$}\put(2.0,7.0){$e_{1}^{0}$}\put(2.0,5.0){$e_{1}^{1}$}\put(2.0,3.0){$e_{1}^{2}$}\put(2.0,1.0){$e_{1}^{3}$}\put(3.0,6.0){$q{2}^{0}$}\put(3.0,4.0){$q{2}^{1}$}\put(3.0,2.0){$q{2}^{2}$}\put(3.0,0.0){$\ddots$}\put(4.0,5.0){$e_{2}^{0}$}\put(4.0,3.0){$e_{2}^{1}$}\put(4.0,1.0){$e_{2}^{2}$}\put(5.0,4.0){$q{3}^{0}$}\put(5.0,2.0){$q{3}^{1}$}\put(5.0,0.0){$\ddots$}\put(6.0,3.0){$e_{3}^{0}$}\put(6.0,1.0){$e_{3}^{1}$}\put(7.0,2.0){$\ddots$}\put(7.0,0.0){$\ddots$}\结束{图片}

此表中的前两列由定义

3.10.8 e(电子)0n个 =0,
n个=1,2,,
q个1n个 =c(c)n个+1/c(c)n个,
n个=0,1,,

其中c(c)n个(0)出现在(3.10.7). 我们继续的手段菱形规则

3.10.9 e(电子)j个k个 =e(电子)j个1k个+1+q个j个k个+1q个j个k个,
j个1,k个0,
q个j个+1k个 =q个j个k个+1e(电子)j个k个+1/e(电子)j个k个,
j个1,k个0.

然后是系数n个S公司-分数(3.10.6)是由提供

3.10.10 0 =c(c)0,
1 =q个10,
2 =e(电子)10,
=q个20,
4 =e(电子)20,
.

商微分算法通常不稳定,可能需要高精度算法或精确算法。更稳定的版本算法在中进行了讨论斯托克斯(1980)适用于贝塞尔函数和Whittaker函数(第章1013),请参阅Gargantini和Henrici(1967).

雅可比分数

形式的连分数

3.10.11 C类=β01α0z(z)β1z(z)21α1z(z)β2z(z)21α2z(z)

称为雅可比分数(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).

以下示例S公司-和J-分数

有关基本函数,请参见§§4.94.35.

有关特殊功能,请参见§5.10(伽马函数),§7.9(错误函数),§8.9(不完整的伽马函数),§8.17(v)(不完整的β函数),§8.19(七)(广义指数积分),§§10.1010.33(贝塞尔函数的商),§13.6(合流超几何函数的商),§13.19(Whittaker函数的商),以及§15.7(超几何函数的商)。

有关更多信息和示例,请参见洛伦岑和瓦德兰(1992第292-330、560-599页)库伊特等。(2008).

§3.10(iii)连分式的数值计算

前向递归算法

这个A类n个B类n个第页,共页(3.10.2)可以通过以下方式计算三项递推关系(1.12.5). 然而,这可能是不稳定;评估时也可能发生溢出和下溢A类n个B类n个(需要不时重新缩放)。

向后递归算法

要计算C类n个第页,共页(3.10.2)我们执行迭代除法

2012年3月3日 u个n个 =b条n个,
u个k个 =b条k个+k个+1u个k个+1,
k个=n个1,n个2,,0.

然后u个0=C类n个.为了达到规定的精确度先验的知识的价值是需要的n个,或n个通过试验确定错误。一般来说,该算法比前向算法更稳定;看见琼斯和王位(1974).

前向序列递归算法

连续分数

2013年3月3日 C类=011121

可以写在表格中

3.10.14 C类=k个=0t吨k个,

哪里

2015年3月3日 t吨0 =0,
t吨k个 =ρk个t吨k个1,
ρ0 =0,
ρk个 =k个(1+ρk个1)1k个(1+ρk个1),
k个=1,2,,.

这个n个第次部分和t吨0+t吨1++t吨n个1等于n个第个收敛的第页,共页(3.10.13),n个=1,2,,与前一项相比本小节中的算法不会出现缩放问题,也不会先验的需要信息。

高茨基(1979年)正向序列算法用于评估不完整伽马函数的连续分数(参见§8.9).

Steed算法

该正向算法在计算收敛点C类n个=A类n个/B类n个,与正向系列相关递归算法。同样,不会出现缩放问题,也不会先验的需要信息。

3.10.16 C类0 =b条0,
D类1 =1/b条1,
C类1 =1D类1,
C类1 =C类0+C类1.

(后向差分算子.)然后针对n个2,

3.10.17 D类n个 =1D类n个1n个+b条n个,
C类n个 =(b条n个D类n个1)C类n个1,
C类n个 =C类n个1+C类n个.

复发持续到(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).