学习目标与准备
本章假定读者已掌握代数、函数、对数及微积分基础。我们从导数的意义出发,研究边际量、连续最优化和数值方程,并区分局部检验与全局结论、连续最优值与可行整数决策,以及数值估计与经过验证的误差界。
1.1 导数是变化率,不是有限差量
在内部点 a a a ,若以下有限极限存在,导数定义为
f ′ ( a ) = lim h → 0 f ( a + h ) − f ( a ) h , f'(a)=\lim_{h\to0}\frac{f(a+h)-f(a)}h, f ′ ( a ) = h → 0 lim h f ( a + h ) − f ( a ) ,
可微性等价于局部展开
f ( a + h ) = f ( a ) + f ′ ( a ) h + r ( h ) , r ( h ) / h → 0. f(a+h)=f(a)+f'(a)h+r(h),\qquad r(h)/h\to0. f ( a + h ) = f ( a ) + f ′ ( a ) h + r ( h ) , r ( h ) / h → 0.
f ′ ( a ) h f'(a)h f ′ ( a ) h 近似有限变化,但通常不等于该变化。若 C ( q ) C(q) C ( q ) 是成本,则 C ′ ( q ) C'(q) C ′ ( q ) 是边际成本 ,单位为每单位数量的成本。平均成本 C ( q ) / q C(q)/q C ( q ) / q 是另一个函数,仅在 q > 0 q>0 q > 0 时定义。
例 1:边际成本与下一单位的成本。 设 C ( q ) = 100 + 5 q + 0.2 q 2 C(q)=100+5q+0.2q^2 C ( q ) = 100 + 5 q + 0.2 q 2 ,则 C ′ ( q ) = 5 + 0.4 q C'(q)=5+0.4q C ′ ( q ) = 5 + 0.4 q 。在 q = 10 q=10 q = 10 时,边际成本为 9 9 9 ,而产量从 10 10 10 增至 11 11 11 的精确成本为
C ( 11 ) − C ( 10 ) = 9.2. C(11)-C(10)=9.2. C ( 11 ) − C ( 10 ) = 9.2.
多出的 0.2 0.2 0.2 是单位变化的二次余项。对较小变化 h h h ,精确差量为 ( 5 + 0.4 q ) h + 0.2 h 2 (5+0.4q)h+0.2h^2 ( 5 + 0.4 q ) h + 0.2 h 2 。
使用常见求导规则:在适当定义域上 ( x n ) ′ = n x n − 1 (x^n)'=nx^{n-1} ( x n ) ′ = n x n − 1 ;( e x ) ′ = e x (e^x)'=e^x ( e x ) ′ = e x ;当 x > 0 x>0 x > 0 时 ( ln x ) ′ = 1 / x (\ln x)'=1/x ( ln x ) ′ = 1/ x ;以及和、积、商与链式法则。求导必须保留原定义域:对数表达式求导并不会使非正输入变为允许的输入。
1.2 单调性与均值定理
为何实数不只有算术规则
有理数支援熟悉的有序算术,但集合 { q ∈ Q : q > 0 , q 2 < 2 } \{q\in\mathbb Q:q>0,\ q^2<2\} { q ∈ Q : q > 0 , q 2 < 2 } 没有有理数最小上界。实数系统加入完备性:非空有上界集合都有最小上界。这不只是另一条不等式运算规则,而是分析中许多存在性结论的来源。
定义
上界与上确界 对非空 A ⊂ R A\subset\mathbb R A ⊂ R ,若每个 a ∈ A a\in A a ∈ A 均满足 a ≤ u a\le u a ≤ u ,则 u u u 是上界。若 s s s 是上界,并对每个上界 u u u 都有 s ≤ u s\le u s ≤ u ,则 s = sup A s=\sup A s = sup A 。最大值是既属于 A A A 又是上界的元素。最大值若存在便是上确界;上确界则不一定属于 A A A 。下界、下确界和最小值以反向不等式定义。
定理
完备性公理 每个非空且有上界的实数子集都有实数上确界。等价地,每个非空且有下界的子集都有实数下确界,因为 inf A = − sup ( − A ) \inf A=-\sup(-A) inf A = − sup ( − A ) 。
例如 sup ( 0 , 1 ) = 1 \sup(0,1)=1 sup ( 0 , 1 ) = 1 ,但 1 1 1 不是该集合的最大值;相反,sup [ 0 , 1 ] = 1 \sup[0,1]=1 sup [ 0 , 1 ] = 1 且是最大值。依此定义,空集及没有上界的集合没有实数上确界。
定理
上确界近似判据 非空集合 A A A 的上界 s s s 是 sup A \sup A sup A ,当且仅当对每个 ε > 0 \varepsilon>0 ε > 0 ,存在 a ∈ A a\in A a ∈ A 使 s − ε < a ≤ s s-\varepsilon<a\le s s − ε < a ≤ s 。
证明
证明两个方向 若 s = sup A s=\sup A s = sup A ,但没有元素大于 s − ε s-\varepsilon s − ε ,则 s − ε s-\varepsilon s − ε 是更小的上界,矛盾。反过来,假设 s s s 是上界且具近似性质。若另一上界 u < s u<s u < s ,取 ε = s − u > 0 \varepsilon=s-u>0 ε = s − u > 0 ,近似性质给出 a > u a>u a > u ,违反 u u u 是上界。因此 s s s 是最小上界。
例题
上确界证明的两部分 令 A = { n / ( n + 1 ) : n ∈ N , n ≥ 1 } A=\{n/(n+1):n\in\mathbb N, n\ge1\} A = { n / ( n + 1 ) : n ∈ N , n ≥ 1 } 。每个元素小于 1 1 1 ,所以 1 1 1 是上界。给定 ε > 0 \varepsilon>0 ε > 0 ,由下述阿基米德性质选 n > 1 / ε n>1/\varepsilon n > 1/ ε 。则 1 − n / ( n + 1 ) = 1 / ( n + 1 ) < ε 1-n/(n+1)=1/(n+1)<\varepsilon 1 − n / ( n + 1 ) = 1/ ( n + 1 ) < ε ,故 sup A = 1 \sup A=1 sup A = 1 。没有元素等于 1 1 1 ,所以 A A A 没有最大值。只证上界并不足以证明它是最小上界:2 2 2 也是上界。
阿基米德性质与稠密性
定理
阿基米德性质 对每个实数 x x x ,存在正整数 n > x n>x n > x 。因此对每个 ε > 0 \varepsilon>0 ε > 0 ,存在 n ∈ N n\in\mathbb N n ∈ N 使 1 / n < ε 1/n<\varepsilon 1/ n < ε 。
证明
完备性迫使整数无上界 若正整数有上界,完备性给 s = sup N s=\sup\mathbb N s = sup N 。因 s − 1 s-1 s − 1 不是上界,存在 n ∈ N n\in\mathbb N n ∈ N 满足 n > s − 1 n>s-1 n > s − 1 。但 n + 1 > s n+1>s n + 1 > s 且 n + 1 ∈ N n+1\in\mathbb N n + 1 ∈ N ,违反 s s s 是上界。对 1 / ε 1/\varepsilon 1/ ε 应用该性质即得第二个结论。
任意实数 x < y x<y x < y 之间都有有理数。选 n n n 使 n ( y − x ) > 1 n(y-x)>1 n ( y − x ) > 1 ;由阿基米德性质及整数部分性质,可选整数 m m m 使 n x < m < n y nx<m<ny n x < m < n y ,因此 x < m / n < y x<m/n<y x < m / n < y 。整数部分性质可由整数在正负两方向均无界及自然数良序性得到。x x x 和 y y y 之间也有无理数:选有理数 q q q 使 x − 2 < q < y − 2 x-\sqrt2<q<y-\sqrt2 x − 2 < q < y − 2 ,则 q + 2 q+\sqrt2 q + 2 是位于区间内的无理数。稠密性不表示有理数与无理数集合相同。
例题
以反射求下确界 对 B = { 2 + 1 / n : n ≥ 1 } B=\{2+1/n:n\ge1\} B = { 2 + 1/ n : n ≥ 1 } ,所有元素均大于 2 2 2 。给定 ε > 0 \varepsilon>0 ε > 0 ,选 n > 1 / ε n>1/\varepsilon n > 1/ ε 可得 2 < 2 + 1 / n < 2 + ε 2<2+1/n<2+\varepsilon 2 < 2 + 1/ n < 2 + ε 。所以 2 2 2 是下确界,但不是最小值。最大值为 3 3 3 ,由 n = 1 n=1 n = 1 取到。
收敛中的量词
定义
实数数列的收敛 数列 ( x n ) (x_n) ( x n ) 收敛至 L ∈ R L\in\mathbb R L ∈ R ,指对每个 ε > 0 \varepsilon>0 ε > 0 ,存在 N ∈ N N\in\mathbb N N ∈ N ,使每个 n ≥ N n\ge N n ≥ N 均有 ∣ x n − L ∣ < ε |x_n-L|<\varepsilon ∣ x n − L ∣ < ε 。记作 x n → L x_n\to L x n → L 。N N N 可以依赖 ε \varepsilon ε ,但必须对其后所有下标同时有效。
「收敛至 L L L 」的否定式为:存在 ε 0 > 0 \varepsilon_0>0 ε 0 > 0 ,使对每个 N N N 都有某个 n ≥ N n\ge N n ≥ N 满足 ∣ x n − L ∣ ≥ ε 0 |x_n-L|\ge\varepsilon_0 ∣ x n − L ∣ ≥ ε 0 。少数初始项不影响收敛;在任意后期仍反复出现的失败则会影响。
例题
构造 epsilon–N 证明 对 x n = ( 3 n + 1 ) / ( n + 2 ) x_n=(3n+1)/(n+2) x n = ( 3 n + 1 ) / ( n + 2 ) ,先猜 L = 3 L=3 L = 3 。因为
∣ 3 n + 1 n + 2 − 3 ∣ = 5 n + 2 < 5 n , \left|\frac{3n+1}{n+2}-3\right|=\frac5{n+2}<\frac5n, n + 2 3 n + 1 − 3 = n + 2 5 < n 5 , 选 N > 5 / ε N>5/\varepsilon N > 5/ ε ,则每个 n ≥ N n\ge N n ≥ N 都满足 ∣ x n − 3 ∣ < ε |x_n-3|<\varepsilon ∣ x n − 3∣ < ε 。猜测只是动机;不等式及带量词的选取才建立证明。
定理
唯一性、有界性及极限定律 收敛实数数列的极限唯一,且数列有界。若 x n → x x_n\to x x n → x 、y n → y y_n\to y y n → y ,则 x n + y n → x + y x_n+y_n\to x+y x n + y n → x + y 、x n y n → x y x_ny_n\to xy x n y n → x y ;若 y ≠ 0 y\ne0 y = 0 且商最终有定义,则 x n / y n → x / y x_n/y_n\to x/y x n / y n → x / y 。若最终有 x n ≤ z n ≤ y n x_n\le z_n\le y_n x n ≤ z n ≤ y n ,而两侧均收敛至 L L L ,则 z n → L z_n\to L z n → L 。
证明
唯一性与有界性 若 x n → L x_n\to L x n → L 且 x n → M x_n\to M x n → M ,对任意 ε > 0 \varepsilon>0 ε > 0 ,充分后期的项同时满足 ∣ x n − L ∣ < ε / 2 |x_n-L|<\varepsilon/2 ∣ x n − L ∣ < ε /2 及 ∣ x n − M ∣ < ε / 2 |x_n-M|<\varepsilon/2 ∣ x n − M ∣ < ε /2 。三角不等式给 ∣ L − M ∣ < ε |L-M|<\varepsilon ∣ L − M ∣ < ε ;此式对所有正 ε \varepsilon ε 成立,故 L = M L=M L = M 。要证有界,选 N N N 使 n ≥ N n\ge N n ≥ N 时 ∣ x n − L ∣ < 1 |x_n-L|<1 ∣ x n − L ∣ < 1 。有限个初始项与尾部界 ∣ x n ∣ ≤ ∣ L ∣ + 1 |x_n|\le|L|+1 ∣ x n ∣ ≤ ∣ L ∣ + 1 一起给出整个数列的有限界。
有界性的逆命题不成立:( − 1 ) n (-1)^n ( − 1 ) n 有界但不收敛。若收敛至 L L L ,其偶数子列收敛至 1 1 1 、奇数子列收敛至 − 1 -1 − 1 ,违反极限唯一性。
证明
乘积定律的机制 收敛数列 ( y n ) (y_n) ( y n ) 有界,设 ∣ y n ∣ ≤ M |y_n|\le M ∣ y n ∣ ≤ M 。由
∣ x n y n − x y ∣ ≤ ∣ y n ∣ ∣ x n − x ∣ + ∣ x ∣ ∣ y n − y ∣ , |x_ny_n-xy|\le|y_n||x_n-x|+|x||y_n-y|, ∣ x n y n − x y ∣ ≤ ∣ y n ∣∣ x n − x ∣ + ∣ x ∣∣ y n − y ∣ , 两项均趋零。对商,y n → y ≠ 0 y_n\to y\ne0 y n → y = 0 蕴含最终有 ∣ y n ∣ ≥ ∣ y ∣ / 2 |y_n|\ge|y|/2 ∣ y n ∣ ≥ ∣ y ∣/2 ;分母因而远离零,得 1 / y n → 1 / y 1/y_n\to1/y 1/ y n → 1/ y 。
子列
定义
子列 ( x n ) (x_n) ( x n ) 的子列是 ( x n k ) (x_{n_k}) ( x n k ) ,其中 n 1 < n 2 < ⋯ n_1<n_2<\cdots n 1 < n 2 < ⋯ 为正整数。若 x n → L x_n\to L x n → L ,每个子列都收敛至 L L L ,因为 n k ≥ k n_k\ge k n k ≥ k ,最终会超过任意指定下标。
证明
套叠区间证明 把所有项放入闭有界区间 I 1 I_1 I 1 ,并把它平分。至少一半含无穷多个数列项;选该半区间为 I 2 I_2 I 2 。继续进行,得到长度趋零的套叠闭区间 I k I_k I k ,每个都含无穷多项。选严格递增的下标 n k n_k n k 使 x n k ∈ I k x_{n_k}\in I_k x n k ∈ I k ;因每个区间含无穷多项,故可做到。完备性给出所有区间的一个共同点 L L L :左端点集合的上确界位于每个区间内,因为每个左端点都不超过每个右端点。由 x n k , L ∈ I k x_{n_k},L\in I_k x n k , L ∈ I k ,二者距离不超过区间长度,故趋零。
结论是关于子列,而非整个数列。例如 ( − 1 ) n (-1)^n ( − 1 ) n 有分别收敛至 1 1 1 及 − 1 -1 − 1 的常数子列,但自身不收敛。
实直线闭集:序列判据与开补集
在通常的经典集合论背景下讨论实直线,并允许可数次见证选择。集合 U U U 为开集,是指每个 y y y 属于 U U U 都有某个 r > 0 r>0 r > 0 ,使 ( y − r , y + r ) (y-r,y+r) ( y − r , y + r ) 包含于 U U U 。使用阿基米德性质的推论:1 / n 1/n 1/ n 趋于零。对任意实直线子集 K K K ,下列条件等价:
K K K 中任意在实直线上收敛的序列,其极限仍属于 K K K 。
K K K 的补集为开集。
由 2 推出 1。 设 x n x_n x n 属于 K K K 且 x n x_n x n 趋于 x x x 。若 x x x 不属于 K K K ,补集的开放性给出 r > 0 r>0 r > 0 ,使 ( x − r , x + r ) (x-r,x+r) ( x − r , x + r ) 与 K K K 不相交。收敛性保证充分靠后时 ∣ x n − x ∣ < r |x_n-x|<r ∣ x n − x ∣ < r ,与 x n x_n x n 属于 K K K 矛盾,故 x x x 属于 K K K 。
由 1 推出 2。 固定不属于 K K K 的 y y y 。若以 y y y 为中心的每个正半径区间都遇到 K K K ,则对每个正整数 n n n ,都有 x n x_n x n 属于 K K K 且 ∣ x n − y ∣ < 1 / n |x_n-y|<1/n ∣ x n − y ∣ < 1/ n 。将这些见证选成一个序列。给定 ε > 0 \varepsilon>0 ε > 0 ,取 N N N 使 1 / N < ε 1/N<\varepsilon 1/ N < ε ;当 n ≥ N n\ge N n ≥ N ,有 ∣ x n − y ∣ < 1 / n ≤ 1 / N < ε |x_n-y|<1/n\le1/N<\varepsilon ∣ x n − y ∣ < 1/ n ≤ 1/ N < ε ,故 x n x_n x n 趋于 y y y 。条件 1 会迫使 y y y 属于 K K K ,矛盾。因此 y y y 周围存在避开 K K K 的区间。对补集每点均如此,故补集开放。选择见证序列这一步使用已声明的经典背景;这里不声称在不带选择的基础体系中证明同一结论。若 K K K 为空,其中没有序列,其补集为整个开实直线;若 K K K 为整个实直线,其空补集的开放性真空成立。两种边界情形也满足等价关系。
学习目标与定义
集合 K ⊆ R K\subseteq\mathbb R K ⊆ R 称为紧致的,如果任意一族开集只要其并包含 K K K ,就存在有限个成员仍覆盖 K K K 。集合若包含其点构成的每个收敛序列的极限,则称为闭集;等价地,其补集是开集。序列紧致性指 K K K 中每个序列都有收敛到 K K K 中某点的子列。
空集是紧致、闭且有界的,也满足序列紧致性;这些性质均以空真方式成立。在需要选取点或无穷序列的证明中,要单独处理空集。我们使用先前数列内容中的实数完备性及其闭区间套和 Bolzano–Weierstrass 推论。
定理:实数直线上的 Heine–Borel 定理
对 K ⊆ R K\subseteq\mathbb R K ⊆ R ,下列条件等价:开覆盖意义下的紧致性;闭且有界;序列紧致性。
紧致推出闭且有界。 区间 ( − n , n ) (-n,n) ( − n , n ) (n ≥ 1 n\ge1 n ≥ 1 )覆盖 R \mathbb R R ,因而覆盖 K K K 。有限子覆盖给出一个界。为证明闭性,取 y ∉ K y\notin K y ∈ / K 。对每个 x ∈ K x\in K x ∈ K ,区间 B ( x , ∣ x − y ∣ / 2 ) B(x,|x-y|/2) B ( x , ∣ x − y ∣/2 ) 不含 y y y 。取其中覆盖 K K K 的有限子族,令 δ \delta δ 为这些区间中心到 y y y 的正距离的一半所组成的有限集合的最小值。所选以 x x x 为中心的区间内每点到 y y y 的距离严格大于 ∣ x − y ∣ / 2 ≥ δ |x-y|/2\ge\delta ∣ x − y ∣/2 ≥ δ 。因此 B ( y , δ ) B(y,\delta) B ( y , δ ) 与 K K K 不相交,因此对每个集合外的点,都找到不与集合相交的开邻域,补集是开集。若 K K K 为空,结论立即成立。
闭且有界推出紧致。 将非空的 K K K 放入闭有界区间 I 0 I_0 I 0 。假设 K K K 的某个开覆盖不存在有限子覆盖。将 I 0 I_0 I 0 二等分;至少一个闭半区间与 K K K 的交集不能被有限覆盖,否则两组有限覆盖合起来便覆盖 K K K 。选取这样的半区间并重复,得到长度趋于零的闭区间套 I n I_n I n ,且 K ∩ I n K\cap I_n K ∩ I n 不能被有限覆盖。每个交集非空,选取 x n ∈ K ∩ I n x_n\in K\cap I_n x n ∈ K ∩ I n 。闭区间套性质给出唯一公共点 p p p ;因两者均在 I n I_n I n 中,故 x n → p x_n\to p x n → p 。由闭性得 p ∈ K p\in K p ∈ K 。覆盖中某个成员 O O O 含 p p p ;由其开放性,O O O 包含以 p p p 为中心的某个开球。当 n n n 足够大时,整个区间 I n I_n I n 均在该球内,于是单个开集 O O O 就覆盖 K ∩ I n K\cap I_n K ∩ I n ,矛盾。
闭且有界推出序列紧致。 K K K 中每个序列有界,故 Bolzano–Weierstrass 定理给出收敛子列。由闭性,其极限仍在 K K K 中。
序列紧致推出闭且有界。 若无界,可选取 x n ∈ K x_n\in K x n ∈ K ,使 ∣ x n ∣ > n |x_n|>n ∣ x n ∣ > n ;它没有收敛到有限实数的子列,矛盾。若 p p p 是 K K K 中点的极限,则可选取 x n ∈ K x_n\in K x n ∈ K 使 ∣ x n − p ∣ < 1 / n |x_n-p|<1/n ∣ x n − p ∣ < 1/ n 。由序列紧致性,某子列收敛到 K K K 中一点;它也收敛到 p p p ,由实数极限唯一性得 p ∈ K p\in K p ∈ K 。因此 K K K 是闭集,所有蕴涵均已证明。
实数线上的紧致子集
集合 K ⊂ R K\subset\mathbb R K ⊂ R 若每个开覆盖都有有限子覆盖,便称紧致。在 R \mathbb R R 中,Heine–Borel 定理指出,这等价于闭且有界,也等价于 K K K 内每个数列都有收敛至 K K K 中某点的子列。闭性使子列极限留在集合内;有界性则容许抽取子列。以下证明使用序列形式。
定理
极值定理 若 K ⊂ R K\subset\mathbb R K ⊂ R 非空且紧致,f : K → R f:K\to\mathbb R f : K → R 连续,则 f f f 有界并在 K K K 上取得最大值及最小值。
证明
抽取逼近最大值的子列 若 f f f 无上界,选 x n ∈ K x_n\in K x n ∈ K 使 f ( x n ) > n f(x_n)>n f ( x n ) > n 。存在子列 x n k → x ∈ K x_{n_k}\to x\in K x n k → x ∈ K ;连续性给 f ( x n k ) → f ( x ) f(x_{n_k})\to f(x) f ( x n k ) → f ( x ) ,与其趋向无穷矛盾。对 − f -f − f 作同样论证可得下界。令 s = sup f ( K ) s=\sup f(K) s = sup f ( K ) ,选 x n ∈ K x_n\in K x n ∈ K 使 s − 1 / n < f ( x n ) ≤ s s-1/n<f(x_n)\le s s − 1/ n < f ( x n ) ≤ s 。某收敛子列趋向 x ∈ K x\in K x ∈ K ,由连续性有 f ( x ) = s f(x)=s f ( x ) = s 。对 − f -f − f 同理得最小值。
例题
哪个假设失效? 连续函数 f ( x ) = x f(x)=x f ( x ) = x 在 ( 0 , 1 ) (0,1) ( 0 , 1 ) 上有界,但没有最大值或最小值,因为 0 0 0 、1 1 1 不在定义域。此集合有界但不闭,故不紧致。在 [ 0 , 1 ] [0,1] [ 0 , 1 ] 上两端点都存在,因此取得极值。紧致性是充分的全局假设,不是每个个别函数取得极值的必要条件。
费马必要条件。 在可微的内点局部最大值或最小值处,f ′ ( a ) = 0 f'(a)=0 f ′ ( a ) = 0 。
证明。 在局部最大值处,充分小的正增量给出非正差商,负增量给出非负差商。双侧导数的极限要同时满足两种符号,只能为零。对函数取负处理最小值。□ \square □
内点条件很重要:f ( x ) = x f(x)=x f ( x ) = x 在 [ 0 , 1 ] [0,1] [ 0 , 1 ] 上的极值出现在端点,尽管导数为 1。逆命题不成立:x 3 x^3 x 3 在零处导数为零,但那里不是极值。
Rolle 定理。 若 a < b a<b a < b ,f f f 在 [ a , b ] [a,b] [ a , b ] 上连续、在 ( a , b ) (a,b) ( a , b ) 上可微,且 f ( a ) = f ( b ) f(a)=f(b) f ( a ) = f ( b ) ,则存在 c ∈ ( a , b ) c\in(a,b) c ∈ ( a , b ) 使 f ′ ( c ) = 0 f'(c)=0 f ′ ( c ) = 0 。
证明。 极值定理保证最大值和最小值存在。若函数为常数,任意内点均符合要求。否则某值不同于相等的端点值,因此最大值或最小值之一出现在内部。可在该处应用费马条件。□ \square □
中值定理。 在相同的连续与可微条件下,不要求端点值相等,则存在 c ∈ ( a , b ) c\in(a,b) c ∈ ( a , b ) 满足 f ( b ) − f ( a ) = f ′ ( c ) ( b − a ) f(b)-f(a)=f'(c)(b-a) f ( b ) − f ( a ) = f ′ ( c ) ( b − a ) 。
从 f f f 减去割线 f ( a ) + [ f ( b ) − f ( a ) ] ( x − a ) / ( b − a ) f(a)+[f(b)-f(a)](x-a)/(b-a) f ( a ) + [ f ( b ) − f ( a )] ( x − a ) / ( b − a ) ,再应用 Rolle 定理即可。这一证明不假定 f ′ f' f ′ 连续。推论包括:导数在区间上为零则函数为常数;导数非负则函数不减;区间上 ∣ f ′ ∣ ≤ M |f'|\le M ∣ f ′ ∣ ≤ M 则其中各点满足 ∣ f ( y ) − f ( x ) ∣ ≤ M ∣ y − x ∣ |f(y)-f(x)|\le M|y-x| ∣ f ( y ) − f ( x ) ∣ ≤ M ∣ y − x ∣ 。不能由一点的导数值推出整个区间上的导数界。
由此,f ′ > 0 f'>0 f ′ > 0 在区间内处处成立蕴含严格递增,f ′ < 0 f'<0 f ′ < 0 蕴含严格递减。在线段上有界 ∣ f ′ ∣ ≤ M |f'|\le M ∣ f ′ ∣ ≤ M 可得 ∣ f ( b ) − f ( a ) ∣ ≤ M ∣ b − a ∣ |f(b)-f(a)|\le M|b-a| ∣ f ( b ) − f ( a ) ∣ ≤ M ∣ b − a ∣ 。这些是从局部导数信息得到的区间结论;不能跨过定义域缺失点应用定理。
若 f ′ = 0 f'=0 f ′ = 0 在一个区间上处处成立,则 f f f 在该区间上为常数。
这是必要条件而非充分条件:x 3 x^3 x 3 在零处导数为零,却没有局部极值。不可微点和端点也须检查,它们不属于 Fermat 条件覆盖的内部可微情形。
1.3 求出并验证最优值
当候选点只有有限多个时:
找出内部 f ′ = 0 f'=0 f ′ = 0 或导数不存在的点;
加入两个端点;
计算并比较所有候选值。
一阶导数从正变负给出局部极大值,从负变正给出局部极小值。若函数在驻点附近有连续二阶导数,则 f ′ ′ ( a ) > 0 f''(a)>0 f ′′ ( a ) > 0 给出严格局部极小值,f ′ ′ ( a ) < 0 f''(a)<0 f ′′ ( a ) < 0 给出严格局部极大值;f ′ ′ ( a ) = 0 f''(a)=0 f ′′ ( a ) = 0 无法判定。例如 x 4 x^4 x 4 、− x 4 -x^4 − x 4 和 x 3 x^3 x 3 在零点一、二阶导数都为零,但极值行为不同。
例 2:约束利润。 设价格 p ( q ) = 60 − q p(q)=60-q p ( q ) = 60 − q ,成本 C ( q ) = 100 + 10 q + q 2 C(q)=100+10q+q^2 C ( q ) = 100 + 10 q + q 2 ,其中 0 ≤ q ≤ 20 0\le q\le20 0 ≤ q ≤ 20 。利润为
Π ( q ) = q p ( q ) − C ( q ) = 50 q − 2 q 2 − 100. \Pi(q)=qp(q)-C(q)=50q-2q^2-100. Π ( q ) = qp ( q ) − C ( q ) = 50 q − 2 q 2 − 100.
导数为 50 − 4 q 50-4q 50 − 4 q ,在 12.5 12.5 12.5 以下为正、以上为负。驻点可行,且 Π ( 12.5 ) = 212.5 \Pi(12.5)=212.5 Π ( 12.5 ) = 212.5 、Π ( 0 ) = − 100 \Pi(0)=-100 Π ( 0 ) = − 100 、Π ( 20 ) = 100 \Pi(20)=100 Π ( 20 ) = 100 。因此该区间上的全局最大值为 212.5 212.5 212.5 ,在 q = 12.5 q=12.5 q = 12.5 处取得。若产能限制为 q ≤ 10 q\le10 q ≤ 10 ,利润在可行区间上递增,最大值改为 Π ( 10 ) = 200 \Pi(10)=200 Π ( 10 ) = 200 。
若数量必须是整数,连续最优点未必可行。这里配方得 Π ( q ) = 212.5 − 2 ( q − 12.5 ) 2 \Pi(q)=212.5-2(q-12.5)^2 Π ( q ) = 212.5 − 2 ( q − 12.5 ) 2 ,所以 q = 12 q=12 q = 12 与 q = 13 q=13 q = 13 都使整数利润最大,最大值为 212 212 212 。此处取邻近整数有精确二次形状作依据,不能推广为对所有连续最优点都可直接四舍五入。
1.4 弹性比较相对变化
对正价格 p p p 处的正可微需求 q = q ( p ) q=q(p) q = q ( p ) ,有符号的相对响应为 p q ′ ( p ) / q ( p ) p q'(p)/q(p) p q ′ ( p ) / q ( p ) 。需求递减时,通常报告非负价格弹性
E ( p ) = − p q ′ ( p ) q ( p ) . E(p)=-\frac{p q'(p)}{q(p)}. E ( p ) = − q ( p ) p q ′ ( p ) .
对小幅价格变化,Δ q / q ≈ − E Δ p / p \Delta q/q\approx-E\,\Delta p/p Δ q / q ≈ − E Δ p / p ;这是一阶近似,不是有限百分比变化的精确公式。
价格对应的收入为 R ( p ) = p q ( p ) R(p)=pq(p) R ( p ) = pq ( p ) ,故
R ′ ( p ) = q ( p ) + p q ′ ( p ) = q ( p ) ( 1 − E ( p ) ) . R'(p)=q(p)+pq'(p)=q(p)(1-E(p)). R ′ ( p ) = q ( p ) + p q ′ ( p ) = q ( p ) ( 1 − E ( p )) .
在上述正性假设下,E < 1 E<1 E < 1 时收入随价格局部递增,E > 1 E>1 E > 1 时局部递减。
例 3:幂函数需求。 设 q ( p ) = 100 p − 2 q(p)=100p^{-2} q ( p ) = 100 p − 2 ,p > 0 p>0 p > 0 ,则 q ′ = − 200 p − 3 q'=-200p^{-3} q ′ = − 200 p − 3 ,故各价格处 E = 2 E=2 E = 2 。收入为 100 / p 100/p 100/ p ,随价格递减,与弹性判别一致。
1.5 二分法:带根区间的估计
设实连续函数 f f f 在 [ a , b ] [a,b] [ a , b ] 上满足 f ( a ) f ( b ) < 0 f(a)f(b)<0 f ( a ) f ( b ) < 0 。介值定理保证至少有一个根。二分法计算中点,保留仍有异号端点的半区间并重复。经过 n n n 次二分,区间宽度为 ( b − a ) / 2 n (b-a)/2^n ( b − a ) / 2 n ;其中点到该区间内某个根的距离不超过 ( b − a ) / 2 n + 1 (b-a)/2^{n+1} ( b − a ) / 2 n + 1 。唯一性还需要严格单调等额外假设。
取 m = ( a + b ) / 2 m=(a+b)/2 m = ( a + b ) /2 ;若 f ( m ) = 0 f(m)=0 f ( m ) = 0 ,就停止。
例 4:开始区间计算。 对 f ( x ) = x 3 + x − 1 f(x)=x^3+x-1 f ( x ) = x 3 + x − 1 ,有 f ( 0 ) = − 1 f(0)=-1 f ( 0 ) = − 1 、f ( 1 ) = 1 f(1)=1 f ( 1 ) = 1 。导数 3 x 2 + 1 > 0 3x^2+1>0 3 x 2 + 1 > 0 保证根唯一。在 0.5 0.5 0.5 处 f = − 0.375 f=-0.375 f = − 0.375 ,第一个区间为 [ 0.5 , 1 ] [0.5,1] [ 0.5 , 1 ] ;在 0.75 0.75 0.75 处 f = 0.171875 f=0.171875 f = 0.171875 ,得到 [ 0.5 , 0.75 ] [0.5,0.75] [ 0.5 , 0.75 ] 。中点 0.625 0.625 0.625 与根相距至多 0.125 0.125 0.125 。继续二分可缩小保证误差,而不要求根有方便的精确表达式。
符号计算必须可靠:在近似运算中,接近零而符号不确定的值不能认证为正或负;在声称严格包围前应提高精度或采用区间界。
1.6 Newton 法:快速提议而非保证
在满足 f ′ ( x n ) ≠ 0 f'(x_n)\ne0 f ′ ( x n ) = 0 的迭代点 x n x_n x n ,切线为
y = f ( x n ) + f ′ ( x n ) ( x − x n ) . y=f(x_n)+f'(x_n)(x-x_n). y = f ( x n ) + f ′ ( x n ) ( x − x n ) .
其零点给出 Newton 迭代
x n + 1 = x n − f ( x n ) f ′ ( x n ) . x_{n+1}=x_n-\frac{f(x_n)}{f'(x_n)}. x n + 1 = x n − f ′ ( x n ) f ( x n ) .
它在简单根附近往往收敛很快,但初值不佳、导数为零或迭代点超出定义域都可能导致失败。步长小本身不能证明接近根;若要获得可靠结论,应保留包围区间或建立误差界。
例 5:平方根与残差界。 对 f ( x ) = x 2 − 2 f(x)=x^2-2 f ( x ) = x 2 − 2 ,取 x 0 = 3 / 2 x_0=3/2 x 0 = 3/2 ,则
x 1 = 1 2 ( 3 2 + 2 3 / 2 ) = 17 12 . x_1=\frac12\left(\frac32+\frac{2}{3/2}\right)=\frac{17}{12}. x 1 = 2 1 ( 2 3 + 3/2 2 ) = 12 17 .
第二步得到 x 2 = 577 / 408 x_2=577/408 x 2 = 577/408 。根 α = 2 \alpha=\sqrt2 α = 2 和迭代值都在 [ 1 , 2 ] [1,2] [ 1 , 2 ] 内,而该区间上 f ′ ( x ) = 2 x ≥ 2 f'(x)=2x\ge2 f ′ ( x ) = 2 x ≥ 2 。由均值定理,
∣ x n − α ∣ ≤ ∣ f ( x n ) ∣ 2 . |x_n-\alpha|\le\frac{|f(x_n)|}{2}. ∣ x n − α ∣ ≤ 2 ∣ f ( x n ) ∣ .
经验证的导数下界将残差转化为误差保证。一般而言,单独的 ∣ f ( x n ) ∣ |f(x_n)| ∣ f ( x n ) ∣ 不是误差估计:把方程整体乘以很小的常数可改变残差而不移动根。
1.7 分级练习
基础
若 C ( q ) = 80 + 4 q + q 2 / 2 C(q)=80+4q+q^2/2 C ( q ) = 80 + 4 q + q 2 /2 ,求 q > 0 q>0 q > 0 时的边际成本和平均成本,并求从 q q q 变到 q + h q+h q + h 的精确变化。
对 f ( x ) = ln ( 1 + x ) + x e − x f(x)=\ln(1+x)+x e^{-x} f ( x ) = ln ( 1 + x ) + x e − x 求导并写出实数定义域。
求 f ( x ) = x 3 − 3 x f(x)=x^3-3x f ( x ) = x 3 − 3 x 的驻点,并用导数符号分类。
核心
在 [ 0 , 10 ] [0,10] [ 0 , 10 ] 上最大化 Π ( q ) = 24 q − q 2 − 20 \Pi(q)=24q-q^2-20 Π ( q ) = 24 q − q 2 − 20 ,说明为何无约束驻点不能解决该约束问题。
对 0 < p < 30 0<p<30 0 < p < 30 上的 q ( p ) = 60 − 2 p q(p)=60-2p q ( p ) = 60 − 2 p ,计算 E ( p ) E(p) E ( p ) 并求收入最大处。
举出一个可微函数,使 f ′ ( 0 ) = 0 f'(0)=0 f ′ ( 0 ) = 0 但零点不是局部极值点,并说明理由。
初始根区间长度为 1 1 1 ,二分多少次可保证中点误差至多为 10 − 3 10^{-3} 1 0 − 3 ?
综合
从 x 0 = 2 x_0=2 x 0 = 2 出发,对 f ( x ) = x 2 − 3 f(x)=x^2-3 f ( x ) = x 2 − 3 进行两次 Newton 迭代,并利用 [ 1 , 2 ] [1,2] [ 1 , 2 ] 上的导数界给第二次迭代值一个精确误差界。
对 f ( x ) = x 3 − 2 x + 2 f(x)=x^3-2x+2 f ( x ) = x 3 − 2 x + 2 ,Newton 法从 x 0 = 0 x_0=0 x 0 = 0 开始。计算 x 1 x_1 x 1 和 x 2 x_2 x 2 ,并说明这对无条件收敛说法的意义。
设 f f f 在 [ a , b ] [a,b] [ a , b ] 上连续、内部可微,在区间内有根 α \alpha α ,且内部处处 f ′ ( x ) ≥ m > 0 f'(x)\ge m>0 f ′ ( x ) ≥ m > 0 。证明根唯一,并对每个 x ∈ [ a , b ] x\in[a,b] x ∈ [ a , b ] 证明 ∣ x − α ∣ ≤ ∣ f ( x ) ∣ / m |x-\alpha|\le|f(x)|/m ∣ x − α ∣ ≤ ∣ f ( x ) ∣/ m 。
1.8 完整解答
C ′ = 4 + q C'=4+q C ′ = 4 + q ,C / q = 80 / q + 4 + q / 2 C/q=80/q+4+q/2 C / q = 80/ q + 4 + q /2 。展开得 C ( q + h ) − C ( q ) = ( 4 + q ) h + h 2 / 2 C(q+h)-C(q)=(4+q)h+h^2/2 C ( q + h ) − C ( q ) = ( 4 + q ) h + h 2 /2 ;有限差量包含边际预测中没有的二阶项。
定义域为 x > − 1 x>-1 x > − 1 ,导数为 1 / ( 1 + x ) + e − x − x e − x = 1 / ( 1 + x ) + ( 1 − x ) e − x 1/(1+x)+e^{-x}-xe^{-x}=1/(1+x)+(1-x)e^{-x} 1/ ( 1 + x ) + e − x − x e − x = 1/ ( 1 + x ) + ( 1 − x ) e − x 。
f ′ = 3 ( x − 1 ) ( x + 1 ) f'=3(x-1)(x+1) f ′ = 3 ( x − 1 ) ( x + 1 ) :在 x < − 1 x<-1 x < − 1 为正,在 − 1 < x < 1 -1<x<1 − 1 < x < 1 为负,在 x > 1 x>1 x > 1 为正。因此 − 1 -1 − 1 是局部极大点,1 1 1 是局部极小点。
在 [ 0 , 10 ] [0,10] [ 0 , 10 ] 上 Π ′ = 24 − 2 q > 0 \Pi'=24-2q>0 Π ′ = 24 − 2 q > 0 ,故最大值在 q = 10 q=10 q = 10 ,值为 240 − 100 − 20 = 120 240-100-20=120 240 − 100 − 20 = 120 。无约束驻点 12 12 12 在可行区间外。
E = 2 p / ( 60 − 2 p ) = p / ( 30 − p ) E=2p/(60-2p)=p/(30-p) E = 2 p / ( 60 − 2 p ) = p / ( 30 − p ) 。收入 R = 60 p − 2 p 2 = 450 − 2 ( p − 15 ) 2 R=60p-2p^2=450-2(p-15)^2 R = 60 p − 2 p 2 = 450 − 2 ( p − 15 ) 2 ,因此在 p = 15 p=15 p = 15 最大,此时 E = 1 E=1 E = 1 ;该价格在所给定义域内。
取 f ( x ) = x 3 f(x)=x^3 f ( x ) = x 3 ,其导数 3 x 2 3x^2 3 x 2 给出 f ′ ( 0 ) = 0 f'(0)=0 f ′ ( 0 ) = 0 ,但任意接近零的负输入对应负值、正输入对应正值,故零点既非局部极大也非局部极小。
二分 n n n 次后中点误差至多为 1 / 2 n + 1 1/2^{n+1} 1/ 2 n + 1 。由于 2 9 = 512 < 1000 2^9=512<1000 2 9 = 512 < 1000 、2 10 = 1024 ≥ 1000 2^{10}=1024\ge1000 2 10 = 1024 ≥ 1000 ,取 n = 9 n=9 n = 9 ;八次仅给出较大的界 1 / 512 1/512 1/512 。
迭代为 x n + 1 = ( x n + 3 / x n ) / 2 x_{n+1}=(x_n+3/x_n)/2 x n + 1 = ( x n + 3/ x n ) /2 ,故 x 1 = 7 / 4 x_1=7/4 x 1 = 7/4 、x 2 = 97 / 56 x_2=97/56 x 2 = 97/56 。根与 x 2 x_2 x 2 都在 [ 1 , 2 ] [1,2] [ 1 , 2 ] 内,且 f ′ ≥ 2 f'\ge2 f ′ ≥ 2 。残差为 ( 97 / 56 ) 2 − 3 = ( 9409 − 9408 ) / 3136 = 1 / 3136 (97/56)^2-3=(9409-9408)/3136=1/3136 ( 97/56 ) 2 − 3 = ( 9409 − 9408 ) /3136 = 1/3136 ,所以误差至多为 1 / 6272 1/6272 1/6272 。
在零处 f = 2 f=2 f = 2 、f ′ = − 2 f'=-2 f ′ = − 2 ,故 x 1 = 1 x_1=1 x 1 = 1 ;在一处 f = 1 f=1 f = 1 、f ′ = 1 f'=1 f ′ = 1 ,故 x 2 = 0 x_2=0 x 2 = 0 。迭代在 0 0 0 与 1 1 1 间循环而不收敛,说明每步导数非零仍不能保证收敛。
均值定理给出 f ( v ) − f ( u ) ≥ m ( v − u ) > 0 f(v)-f(u)\ge m(v-u)>0 f ( v ) − f ( u ) ≥ m ( v − u ) > 0 (u < v u<v u < v ),故 f f f 严格递增,至多有一个根。若 x ≠ α x\ne\alpha x = α ,在二者之间应用同一定理,得到某个内部 c c c 使 f ( x ) − f ( α ) = f ′ ( c ) ( x − α ) f(x)-f(\alpha)=f'(c)(x-\alpha) f ( x ) − f ( α ) = f ′ ( c ) ( x − α ) 。取绝对值并用 f ( α ) = 0 f(\alpha)=0 f ( α ) = 0 ,得 ∣ f ( x ) ∣ ≥ m ∣ x − α ∣ |f(x)|\ge m|x-\alpha| ∣ f ( x ) ∣ ≥ m ∣ x − α ∣ ;x = α x=\alpha x = α 时结论显然。
总结
边际量是导数,不是精确的下一单位变化。求最优值须检查定义域和端点,驻点本身并不足够。弹性比较相对变化率,须说明符号约定。二分法保留根区间并给出透明误差界;Newton 法提供有效迭代,但须有保障或分析才能断言收敛与精度。