版本11.2--追求极致的极限

Asdmath2020

█ 本文译自算法R&D，内核开发工程师 Devendra Kapadia 于2017年11月9日的博客文章： Limits without Limits in Version 11.2.

0 D$ W- @. s9 i4 |: g
这是一个序列的10个项：
: |3 H4 _" o3 M1 B

这是它们的数值：
0 t& q* I) a& M3 W

这个数列的极值是什么？如果这个数列无限继续下去，我们会得到什么？
" j% n' D( j* \8 y

在 Mathematica 和 Wolfram 语言中，有一个直接求解这一问题的函数：

! ~) R) u. E0 c+ a+ r* X

7 f7 j7 U) S" y8 W3 F

极限是包括数论、几何学和计算复杂性理论在内的许多领域的中心概念。它们还是微积分的核心，不仅是因为它们被用来定义导数和积分的概念。

7 V" Z% }8 m% ~

计算极限是 Mathematica 和 Wolfram 语言一直以来就具有的功能。在版本11 .2 中，这一功能被大幅扩展。我们综合利用 Wolfram 语言的许多特色功能，并发明了一些全新的算法来实现这一点。为了确保我们所做的一切满足用户的需要，我们从Wolfram|Alpha 抽选了100多万个极限问题。

让我们来讨论一下哈代和拉马努金（Hardy 和 Ramanujan：http://blog.stephenwolfram.com/2016/04/who-was-ramanujan/）在1918年所做的一个极限问题。首先，定义如下序列 ：

8 O4 R& D* b; ^! L. q! S
3 M1 t5 K5 {: v6 s  C" T
这是该数列前10个值的表格：
# ]% a4 b- K* m  c9 |' P

下图表明，当 n 趋向于无穷时，数列收敛于 0。
) L, f9 _% y* j- I6 Q1 U6 i
( z7 N' p1 Y+ T, g0 N. d

在版本11 .2 中引入的  DiscreteLimit  函数证实了这个数列的极限是0。

# F. L3 X' s0 R3 O( J在现实生活中出现的许多序列 (例如，在信号通信中) 是周期性的，它们的值在规则的间隔时间内重复出现。此类间隔的长度称为序列的周期。例如，考虑一个使用 Mod 定义的序列。
/ w2 x, p% b0 V  z) h
' J7 k9 r! |7 j" f
* y+ T- q  j5 H/ a5 K+ @$ }序列图表明该序列是周期性的，周期为6。

( Y' ]7 c& R- J. D' U9 u5 x/ l/ o
* L% B$ K) w' I8 k: L" {/ v与我们的第一个例子相反，这个序列不收敛，它在0和5之间振荡。因此，在这种情况下， DiscreteLimit  返回 Indeterminate（不定）。
/ n5 M" a8 b. }) n0 {8 q

5 c+ k) Q2 ~% e0 K# {在这种情况下，版本11.2 的新增函数 DiscreteMinLimit 和 DiscreteMaxLimit 可以分别用于计算振荡的下限和上限。因此，我们有:

! t3 i9 }2 A& a, \( y
0 f: G0 x" [8 x7 o$ o% g
在数学文献中，DiscreteMinLimit 和 DiscreteMaxLimit 通常被称为 "lim inf" 和 "lim sup"。这些极限的传统底线和顶线符号是可用的，如下所示。

: k- G) z5 d, `% I7 C
' V) Y9 r, v/ X8 [" w
: A8 c- h$ K" k7 ?* [; X我们的下一个例子是由三角函数 Sin 和 Cos 得到的一个振荡序列，定义如下。
3 c  K6 @7 `3 I# m
  W2 d5 H" W9 s

尽管在实数范围内 Sin 和 Cos 作为函数是周期性的，这个整数序列的行为很奇特，远远不是一个周期序列，可以通过下图证实。

因此，该序列的极限不存在。
" w. _' |9 C# V2 J; i
+ o! _5 x/ A% R# `  K' E0 P0 M3 @
然而，事实证明，对于这种 "密集型非周期序列"，极值可以通过将它们作为实际函数来计算。不出所料，DiscreteMinLimit使用此方法返回该示例的答案0。


* [! R( b* p& z) C使用相同的方法，DiscreteMaxLimit返回一个含有根式对象的相当凌乱的结果。
+ N6 |7 [! K$ P( N* L6 `

1 Y9 W5 Q% P; q) \  K这个结果的数值近似于 0.8，正如人们可能从图中猜到的那样。


8 L, i. q" I/ M0 b" T当我们试图计算无限嵌套的根式值时，离散极限也会以自然的方式出现。例如，考虑以下嵌套根式的求值问题。

" X2 g# U. P; U+ x5 r! |

由于该序列满足非线性重复性，因此可以使用 RSolveValue 来生成根式展开式的连续项。例如，展开式的第三个项可以通过以下所示得到。

通过 RecurrenceTable 所生成的下图可以看出，无限嵌套的根式值似乎是 2。

使用版本 11.2，通过求解 RSolveValue 中的值 r(∞)，我们可以确认极值确实是 2，如下所示。
2 ?7 N1 f* }( f/ m
7 \4 H$ T2 q+ {6 ^) F

极限的研究是数学的一个分支，称为渐近分析。渐近分析提供了在特定值 (如0或无穷大) 附近获取问题近似解的方法。事实证明，在实际中，渐进逼近的效率通常恰恰会在相应的精确计算变得困难的情况下得到提高！这一现象的一个著名的例子是在研究整数分区时，随着数量的增加，分区会增长非常快。例如，数字6可以使用IntegerPartitions 以11种不同的方式进行分区，如下所示。

$ c# u' [6 U# h* q/ ?- V) F. K8 h

可以直接使用 PartitionsP 找到不同分区的数量，如下所示。

如前所述，分区的数量随着整数的大小而快速增长。例如，数字200有近4兆分区。

) R/ c0 ^9 C! G( U$ o
在 1918年，哈代和拉马努金提供了这个数字的渐近逼近，由以下公式给出。

+ c& h' i/ o" L; v, |, h

根据上式估计，数字200的答案相当接近4兆。

8 B3 Z% }7 D4 S
提供有一个更大的整数，我们几乎可以瞬间得到一个更好的分区数量近似，如下例所示。
. p4 l* j  O( A6 {

最后，我们可以证实，使用 DiscreteLimit，当 n 趋于无穷时，渐近估计接近于分区数，这正印证了上面讨论的哈迪-拉马努金公式。
1 Y8 l/ g+ w% u% [' U
. W, B( e  |5 P3 d

正式的表述是，当n 接近无穷时，精确和近似公式的分区数是渐近等价的。

, Y  \$ `( Q+ L, x

渐近概念在函数极限的研究中也起着重要的作用。例如，三角学中的小角度逼近断言 "在 x 取很小的值时，sin (x) 近似等于 x 。" 还可以表述成， "当 x 趋向于 0 时，sin (x) 与x渐近等价。" 这个结果可以使用计算函数极值的 Limit 正式表述如下。

0 G9 g+ F5 ]' k" Z) q2 R2 k. R
这个图形直观地证实，极限确实是1。
8 x6 F( }( f$ ~  H, b5 ?9 G

) {% l+ W: N$ ?2 O上述极限也可以使用洛必达法则，通过分别计算分子和分母的导数 cos(x) 和1 得到，如图所示。

9 Y* b$ B& s1 u: I4 V0 a- W" E

洛必达法则为计算许多极限提供了一个强有力的方法。然而，在得出答案之前，可能需要大量的步骤。例如，考虑以下极限。

( q) H" r9 n3 x

由于所有的中间计算都给出了不确定的结果，这个极限要求重复应用六次洛必达法则规则才能得到答案0。

8 m( |" v5 @" ~8 @& V  ?
  d8 u) \3 ^, ]$ j  r

我们看到在求函数极值时，洛必达法则作为一个实用的算法，仍存在一定的局限性，因为我们不可能确定算法应该在什么时候停止！而内置函数 Limit 使用级数展开和现代算法的组合，可以对涉及指数和对数的输入（即所谓的"exp-log"类）很好地工作。事实上，在版本11 .2 中，Limit 得到了大幅更新，可以以相当全面的方式，对各种难题迎刃而解，如下所示（最后两个示例仅适用于最新版本）。

& W, {; ^6 R: \* y2 [9 W5 ~

( n0 u7 Z* d- f, ?- \. e4 f

与序列的情况相同，周期函数和振荡函数的极限往往不存在。你可以使用 MaxLimit 和MinLimit，就像它们对应的离散情况一样，给函数的振荡设置严格的边界，就像下面这个经典的例子所示。

6 L+ Q7 [$ y+ A/ Y. j
. n# p: T6 b5 m0 \3 g如图所示，函数在-1和1之间接近0处快速振荡。这些界限由 MaxLimit 和 MinLimit确认，而 Limit 本身返回 Indeterminate。

在前面的示例中，由于函数在原点附近剧烈振荡，极限不存在。不连续函数还可以提供在某点处极限可能不存在的其他类型示例。我们这里考虑这样一个函数，它在原点和其他值上有跳跃不连续性。该函数由 SquareWave 和 FresnelS  定义，如下所示。

从这个图可以看到跳跃不连续性，这是由于函数定义中存在方波所导致。
! b, l- \* U6 _$ h
- h% s# c2 ?- i; J* F

我们看到，取决于我们接近原点的方向，函数在0点有极值。可以使用 Direction （方向）选项计算右侧 ("FromAbove")的极值。

; i, Z: ]" [4 |* Q. Q

同样，左侧的极值可以计算如下。

极值如果存在的话，是函数的 "双面" 极值，在本例中不存在。

2 M7 o1 l, A( K2 n* a& T

默认情况下，在版本11 .2 中 Limit 计算双面极值。这与早期版本不同，后者在默认情况下，从上面计算极值。因此，在不设定 Direction 选项时， 我们从Limit 得到一个不确定的结果 Indeterminate。
3 z; k" u% Y% }$ O

方向极值在多变量情况下更有意义，因为在更高的维度中有许多可能的方向接近给定点。例如，考虑定义如下的二元函数 f (x，y)。

; o$ J3 U4 ^& z" |9 q6 W

如果我们沿着 x 轴（由 y=0 给出）求解，则这个函数在原点的极值是 0，因为函数沿着这条线有常数值0。

同样，如果我们沿着 y 轴（即x=0），则函数在原点的极值为0。
) t2 k3 i: [" X3 P( G9 i
但是，如果我们沿着线 y=x 的方向接近原点，则极限为 1/2，如下所示。

" S5 C% B- ]' Z; W7 Q1 v6 Y% D更一般地，当我们沿着不同的线 y=m x 接近原点时，极值会发生变化。
. @& W0 A) \3 a1 i! v, ^' x
2 @5 a4 r0 n, N0 @  V% |极值的方向依赖性意味着不存在真正的多元极值。在版本11 .2 中，Limit 可以轻松处理多变量情况，并且快速返回该函数在原点上极值的预期答案 Indeterminate。
3 W2 P, _. h! X& {/ ^. c. z
曲面z=f(x，y)的图形确认函数在原点附近的行为。

答案可以通过对函数直接应用 Limit 来确定。

) i# T0 o7 u$ V5 x5 E) B+ }用户每天从 Wolfram|Alpha 所输入的查询语句是我们获取多变量极值示例的一个丰富来源。我们获得了大约10万个匿名 Wolfram|Alpha 查询，然后用版本11 .2 运算。这个示例是最新版本中 Limit 能轻松求解的浩如烟海的问题之一。
5 m* z  ^( x2 Y9 _- ?

& C. p' Z# Z7 w! _8 A6 R浏览 Wolfram|Alpha 示例本身就乐趣无穷，我们在此与你分享其中1000个精彩示例（在后台发送“极限”，便可获取这1000个范例的下载链接）。示例图像如下所示。

% |7 o9 F7 @, y, ?% p1 Q7 `

版本11 .2 能够计算整个集合的90 %，这是相当了不起的，因为多元极限是本版本推出的最新功能。

, f( J) t$ K4 b: \5 _! U对于更大的集合——来自 Wolfram|Alpha 的100万个单变量极限问题，版本11 .2 能够求解的比例更高（96%），而版本11 .1为94 %。这两个版本之间的微小百分比差异可以解释为，大多数Wolfram|Alpha单变量极限查询与大学微积分中的第一或第二个课程有关，在任一版本中都很容易用 Limit 计算。
3 P1 \* |% h: L0 y; t
" G0 J! @* k! L2 ?9 [- l% F自 Wolfram 语言首次发布版本 1 开始（1988年），Limit 一直是最可靠的函数之一。版本11 .2 中对该函数的改进，以及 DiscreteLimit 和其他新函数的推出，为我们求解极限世界的难题提供了便利。希望您能喜欢本文简介，并欢迎您对新功能提出任何评论或建议。
) e- s% z9 v( e+ z& `
# \; `7 g, @) C8 B9 h