学习目标与公理起点
从目标误差向后规划 epsilon 证明、追踪条件依赖,并完整证明递推数列收敛,而不是只解不动点方程。假设实数构成完备有序域:满足域的恒等式与逆元、相容的全序,且每个非空有上界子集都有上确界。有序域公理本身并不蕴含完备性。
完备性推出单调有界收敛定理。若 xn 非递减且有上界,令 L=sup{xn:n≥1}。对任意 ε>0,上确界的逼近性质给出指标 N,使 L−ε<xN≤L。对所有 n≥N,由单调性得 L−ε<xn≤L,因此 xn→L。非递增情形可对数列取负后得到。这里清楚显示了完备性用在何处。
1. 乘积误差规划需要先建立缓冲界
设 un→a 且 vn→b。为证明 unvn→ab,将误差拆成
∣unvn−ab∣≤∣un∣∣vn−b∣+∣b∣∣un−a∣.
第一个系数起初不是固定常数。先用容差一及收敛性选取 N0,使 n≥N0 时 ∣un∣≤∣a∣+1=:M。此时 M>0 与 n 无关。
给定任意 ε>0,选取 N1 使 ∣vn−b∣<ε/(2M),再选取 N2 使 ∣un−a∣<ε/[2(∣b∣+1)]。当 n≥max(N0,N1,N2) 时,上述误差小于 ε/2+ε/2=ε。使用 ∣b∣+1 可避免 b=0 时除以零。
依赖次序很重要:先选固定界 M,再选依赖它的误差阈值;最终阈值依赖 ε,但不依赖随后取的指标 n。若在选 N 时使用尚未受到控制的指标上的 un,就颠倒了量词次序。
2. 递推定义的五项独立任务
对递推式 xn+1=F(xn),检查:
- **良定义性:**所有迭代值都留在 F 的定义域内。
- **不变量:**某个区域条件或不等式经更新后仍成立。
- **极限存在:**例如使用单调有界收敛或 Cauchy 估计。
- **识别极限:**只有在极限存在后,才能用有效的连续性结论取极限。
- **误差:**若需要近似,应推导收敛速率,而非只断言收敛。
解 L=F(L) 只涉及第四项;即便如此,它也只是给出候选值,并未证明任何数列趋向该值。递推 xn+1=1−xn、x0=0 有唯一不动点 1/2,但迭代值在零和一之间交替。不动点方程本身不能证明收敛。
3. Newton 平方根迭代:从定义到极限
设 A>0、x0>0,定义
xn+1=21(xn+xnA).
令 α>0 为 A 的唯一正实平方根,其存在性可独立由完备性或介值定理得到。这里不假设该递推本身已经构造了平方根。
良定义性。 正输入产生正输出,故由归纳法得每个 xn>0,所有除法均有意义。
不变量界。 直接代数运算得到
xn+1−α=2xn(xn−α)2≥0.
因此从 x1 开始的每个迭代值都不小于 α,即使 x0<α 也如此。后一情形中,声称从任意正初值出发整个数列都递减是错误的。
第一步之后的单调性。 对 n≥1,
xn+1−xn=2xnA−xn2≤0,
因为 xn≥α。因此尾部数列非递增且以下界 α 为界。由完备性得到极限 L≥α>0。
识别极限。 非零点处倒数函数的连续性允许取极限:L=(L+A/L)/2,故 L2=A。由正性及正平方根唯一性,得 L=α。先证明收敛并确认极限分母非零,再解此方程才是有效的。
4. 推导误差速率与精确变换递推
对 n≥1 令 en=xn−α。不变量恒等式给出
en+1=2xnen2≤2αen2,en+1≤21en.
第二个不等式来自 en=xn−α<xn;当 en=0 时取等。它给出初等的几何误差估计;第一个不等式则表明误差相对于 α 足够小时具有二次改善。
还可作一个精确变换。对每个 n≥0,令 rn=(xn−α)/(xn+α);分母为正,且此比值的绝对值小于一。代数运算给出
xn+1+α=2xn(xn+α)2,rn+1=rn2.
因此 rn=r02n 且 xn=α(1+rn)/(1−rn)。这提供了另一条收敛途径,并解释精度迅速提高的原因。它不是计算未知 α 的实用公式,因为变换本身使用了 α;其作用是证明与误差分析。
经过验证的例子。 对 A=2,x0=1,精确有理数迭代得到 x1=3/2、x2=17/12 和 x3=577/408。由 (7/5)2<2<(3/2)2 得 0<e1<1/10 和 α>1,因此 e2<1/200、e3<1/80000。有理区间 [577/408−1/80000,577/408] 包含二的正平方根。误差界由不等式证明,而非由小数显示相符推断。
5. 带必要假设的子列判别
若有界实数列只有一个可能的子列极限 L,则它收敛到 L。假设不收敛,由收敛的否定可取 ε0>0,并取一子列使其每项到 L 的距离至少为 ε0。由有界性及 Bolzano–Weierstrass 定理,该子列有收敛子列。绝对值的连续性说明其极限仍至少相距 ε0,与可能的子列极限唯一矛盾。
有界性不可缺少。偶数项为零、奇数项为 n 的数列,唯一有限子列极限为零,但它不收敛,因为奇数项无界。有限聚点唯一并不能单独保证所有项最终靠近它。
分级证明规划练习
- 修正量词。 为什么不能选择依赖随后指标 n 的收敛阈值?用先行固定界和只依赖 ε 与给定数列的阈值,重写乘积极限证明计划。
- 不动点分析。 对任意实数 x0 分析 xn+1=1−xn,准确确定哪些初值使数列收敛,并解释为什么只有代数候选值不够。
- 第一步例外。 在 Newton 迭代中,若 0<x0<α,证明第一步跨到 α 的另一侧,并说明哪一单调性结论仍成立。
- 误差传播。 设 α≥m>0 且 0≤eN≤η,仅用 m,η 给出 eN+1、eN+2 的上界,解释绝对精度与相对精度的区别。
- 子列证明。 证明有界且仅有一个可能子列极限的判别法,并用明确的指标记法给出无界反例。
- Cauchy 分析。 证明 xn=n 相邻增量趋于零,但数列不是 Cauchy 数列。解释为什么 Cauchy 定义控制充分靠后的所有项对,而不仅是相邻项。
完整解答
- 极限定义要求同一个 N 适用于每个 n≥N;依赖指标的阈值改变了量词次序,可能没有验证任何共同尾部。先固定 N0 使 ∣un∣≤∣a∣+1=M。给定 ε,由收敛性选取使 ∣vn−b∣<ε/(2M) 和 ∣un−a∣<ε/[2(∣b∣+1)] 成立的阈值。将它们与 N0 取最大值,所得阈值与后续的 n 无关,并控制整个乘积误差。
- 偶数项等于 x0,奇数项等于 1−x0。两者相等时且仅此时收敛,即 x0=1/2,此时数列为常数;否则两个常值子列具有不同极限。不动点方程只识别可能的极限,不能决定任意迭代的行为。
- 恒等式 x1−α=(x0−α)2/(2x0) 在 x0=α 时严格为正,因此尽管 x0<α,仍有 x1>α。数列不能从指标零开始非递增;从指标一开始则非递增且以下界 α 为界,分别由不变量和差值计算得到。
- 二次误差恒等式给出 eN+1≤η2/(2m),进而得 eN+2≤η4/(8m3)。这些绝对误差界与平方根的单位相同。相对误差要除以 α;其大小至多为对应绝对误差界除以 m,但实际分母是 α。只有把数值误差界与指定舍入容差比较后,才能确定小数位数。
- 若不收敛,则存在无穷子列始终与极限相距至少 ε0。有界性提供收敛的子列之子列,其极限仍与原候选极限分离,矛盾。反例取 x2k=0、x2k−1=2k−1。任何收敛子列只能包含有限多个奇指标,因为对应值趋于无穷,因此其极限为零;但整个数列无界且不收敛。
- 有理化得到 n+1−n=1/(n+1+n)→0。然而对指标 n 与 4n,差为 4n−n=n,在任何指定尾部之后都可任意大。故以固定容差一即可否定 Cauchy 条件。相邻小变化可在多步中累积;Cauchy 条件排除了所有充分靠后的项对之间出现这类累积。
总结与术语
不变量在每次迭代后仍成立。不动点在收敛得到证明前只是候选极限。缓冲界在选取 epsilon 阈值前把可变系数控制为固定常数。子列极限是某个子列的极限;有界性通过完备性提供这样的极限。严谨近似论证需分别建立迭代值有效、不变量控制、极限存在、极限识别及误差界。