主题
Search

施图姆函数


给定函数 f(x)=f_0(x),记 f_1=f^'(x) 并通过下式定义施图姆函数

 f_n(x)=-{f_(n-2)(x)-f_(n-1)(x)[(f_(n-2)(x))/(f_(n-1)(x))]},
(1)

其中 [P(x)/Q(x)] 是多项式商。然后构造以下施图姆函数链,

f_0=q_0f_1-f_2
(2)
f_1=q_1f_2-f_3
(3)
f_2=q_2f_3-f_4
(4)
|
(5)
f_(s-2)=q_(s-2)f_(s-1)-f_s,
(6)

被称为 施图姆链。当获得常数 -f_s(x) 时,链终止。

施图姆函数为求代数方程在给定区间内实根的数量提供了一种便捷的方法。具体来说,在两个点 x=ax=b 处评估的施图姆函数之间的符号变化数之差,给出了区间 (a,b) 内的实根数。这个强大的结果被称为 施图姆定理。然而,当数值应用该方法时,在计算多项式商时必须小心,以避免由于舍入误差而产生虚假结果。

SturmFunction

作为施图姆函数在寻找 多项式 的一个具体应用,考虑函数 f_0(x)=x^5-3x-1,如上图所示,它有根 -1.21465-0.3347340.0802951+/-1.32836i 和 1.38879(其中三个是实数)。导数f^'(x)=5x^4-3 给出,然后 施图姆链 由下式给出

f_0=x^5-3x-1
(7)
f_1=5x^4-3
(8)
f_2=1/5(12x+5)
(9)
f_3=(59083)/(20736).
(10)

下表显示了 f_i 的符号和为间隔 Deltax=2 分隔的点获得的符号变化数 Delta

xf_0f_1f_2f_3Delta
-2-11-113
0-1-1111
211110

这表明 3-1=2 个实根位于 (-2,0) 中,并且 1-0=1 个实根位于 (0,2) 中。将间隔减小到 Deltax=0.5 得到下表。

xf_0f_1f_2f_3Delta
-2.0-11-113
-1.5-11-113
-1.011-112
-0.51-1-112
0.0-1-1111
0.5-1-1111
1.0-11111
1.511110
2.011110

该表隔离了三个实根,并表明它们位于区间 (-1.5,-1.0)(-0.5,0.0)(1.0,1.5) 中。如果需要,可以进一步缩小根所在的区间。

施图姆函数满足以下条件

1. 在区间内的任何点,两个相邻的函数不会同时消失。

2. 在施图姆函数的零点,其两个相邻函数具有不同的符号。

3. 在 f_0(x) 的零点周围的足够小区间内,f_1(x) 处处大于零或处处小于零。


另请参阅

笛卡尔符号法则, 施图姆链, 施图姆定理

使用 Wolfram|Alpha 探索

参考文献

Acton, F. S. 数值方法(第二版印刷本)。 Washington, DC: Math. Assoc. Amer., p. 334, 1990.Dörrie, H. "施图姆关于根数的问题。" §24 in 初等数学的 100 个伟大问题:其历史和解决方案。 New York: Dover, pp. 112-116, 1965.Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. FORTRAN 数值食谱:科学计算的艺术,第二版。 Cambridge, England: Cambridge University Press, p. 469, 1992.Rusin, D. "已知数学。" http://www.math.niu.edu/~rusin/known-math/96/sturm.Sturm, C. "关于数值方程解法的回忆录。" Bull. des sciences de Férussac 11, 1929.

在 Wolfram|Alpha 上被引用

施图姆函数

请引用为

Weisstein, Eric W. "施图姆函数。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/SturmFunction.html

学科分类