$\gdef\e{\mathrm{e}}\gdef\d{\mathrm{d}}\gdef\i{\mathrm{i}}\gdef\N{\mathbb{N}}\gdef\Z{\mathbb{Z}}\gdef\Q{\mathbb{Q}}\gdef\R{\mathbb{R}}\gdef\C{\mathbb{C}}\gdef\F{\mathbb{F}}\gdef\E{\mathbb{E}}\gdef\P{\mathbb{P}}\gdef\M{\mathbb{M}}\gdef\O{\mathrm{O}}\gdef\b#1{\boldsymbol{#1}}\gdef\ker{\operatorname{Ker}}\gdef\im{\operatorname{Im}}\gdef\r{\operatorname{rank}}\gdef\id{\mathrm{id}}\gdef\span{\operatorname{span}}\gdef\spec{\operatorname{spec}}\gdef\mat#1{\begin{bmatrix}#1\end{bmatrix}}\gdef\dat#1{\begin{vmatrix}#1\end{vmatrix}}\gdef\eps{\varepsilon}\gdef\arcsinh{\operatorname{arcsinh}}\gdef\arccosh{\operatorname{arccosh}}\gdef\arctanh{\operatorname{arctanh}}\gdef\arccoth{\operatorname{arccoth}}\gdef\arcsech{\operatorname{arcsech}}\gdef\arccsch{\operatorname{arccsch}}\gdef\sgn{\operatorname{sgn}}\gdef\sech{\operatorname{sech}}\gdef\csch{\operatorname{csch}}\gdef\arccot{\operatorname{arccot}}\gdef\arcsec{\operatorname{arcsec}}\gdef\arccsc{\operatorname{arccsc}}\gdef\tr{\operatorname{tr}}\gdef\unit#1{\mathop{}!\mathrm{#1}}\gdef\re{\operatorname{Re}}\gdef\aut{\operatorname{Aut}}\gdef\diag{\operatorname{diag}}\gdef\D{\mathrm{D}}\gdef\p{\partial}\gdef\eq#1{\begin{align*}#1\end{align*}}\gdef\Pr{\mathsf{Pr}}\gdef\Ex{\mathsf{E}}\gdef\Var{\mathsf{Var}}\gdef\ip#1{\left\langle #1\right\rangle}\gdef\J{\mathrm{J}}\gdef\w{\wedge}$
封面 from wxf 的课件。
这次选择以问答的方式做笔记帮助复习。不过事实证明基本定义清楚之后刷题就够了。括号里面的编号表示出自 wxf ppt,如果想知道答案可以评论问。
多元函数极限与连续
-
距离,三个公理。对于线性空间,加平移不变性。
哪三个公理?(1-16)
-
范数,三个公理。线性空间下的距离和范数一一对应。
哪三个公理?(1-18)
$p$-范数,必须要求 $p\ge 1$。
证明 $p\ge 1$ 可以构成范数,$p<1$ 不行。
答案
https://www.zhihu.com/question/667693249/answer/3633863935
对于 $p<1$,$\|(0,1)+(1,0)\|>\|(0,1)\|+\|(1,0)\|$。 -
内积,三个公理。内积可以定义范数,范数不一定可以由内积定义。$p$-范数中只有 $2$-范数可以。范数可以由内积得到,当且仅当 $\|\b x+\b y\|^2+\|\b x-\b y\|^2=2(\|\b x\|^2+\|\b y\|^2)$。(1-25)
证明?
-
拓扑相关的定义,略。
当我们讨论连续性和极限时,对定义域的限定不同,具体是怎样的?(1-29)
答案
极限要求 $\b x_0$ 是聚点,连续性要求 $\b x_0$ 属于定义域。如果直接把连续性定义为极限等于该点的取值,可是可以,但是如果 $\b x_0$ 是个孤立点,就没法说了。在 wxf 讲的连续性中,孤立点处必连续。 -
任意范数在以下意义下与 $\|\cdot\|_{\infty}$ 等价:$\exists M_1,M_2$,使得 $\forall\b x$,$M_1\|\b x\|_{\infty}\le\|\b x\|\le M_2\|\b x\|_{\infty}$。(1-32)从而与有界、收敛、连续等相关的定理都可以只证明无穷范数下的情况。
证明?(1-35)
答案
$\|\b x\|\le C\|\b x\|_{\infty}$:写出 $\b x$ 的基的表示,用三角不等式拆开。
另一个方向:考虑反证,如果有一列 $\set{\b x_n}$,$\|\b x_n\|_{\infty}$ 比 $\|\b x_n\|$ 大得越来越离谱,那么可以找到一个范数的收敛子列,从而是柯西列,从而对应的 $\b x$ 子列也收敛(这个没有循环论证,可以先证柯西列收敛的),这个极限点处的两种范数就爆了。 -
一元微积分里的一些性质在多元情况也成立:柯西收敛定理、致密性定理、最值定理、闭集连续 = 一致连续。多元函数除了没有一元函数的单调性,别的都可以讨论。
证明?
-
证明:$S\subseteq V$ 是有界闭集,当且仅当 $S$ 中的任何点列的极限也在 $S$ 中。(1-37)
-
证明:如果 $S\subseteq\R^n$ 既是开集又是闭集,则它只能是全集或空集。
提示
找到一条从 $x\in S$ 到 $y\notin S$ 的路径 $\pi$,构造 $f:\pi\to\set{0,1}$。 -
证明多元情况下的闭区间套定理。
-
在证明多元连续函数的含参积分对参数的连续性时,必须搞出一致连续性才行。(1-41)
现在我们说 $(x,y_1,\cdots,y_n)\in[a,b]\times A$,$x$ 是被积变量,$\b y_0=(y_1,\cdots,y_n)$ 是参数。如何找到一个闭集?
答案
如果 $\b y_0$ 是内点,那选个完全包含在 $A$ 内的闭邻域就行。否则,如果 $A$ 是闭集,那随便选个闭邻域就行。否则,证明起来有点麻烦。我想到的一个思路是,对于一个 $\eps$,对于每个 $x\in [a,b]$,找到一个 $\delta(x)$ 为最大的($\sup$),不超过某个阈值的,使邻域内的函数值都 $\le\eps/2$ 的 $\delta$。可以通过测度的三角形不等式证明 $\delta(x)$ 连续,从而一致连续,从而 $\inf\delta(x)>0$ 就是选择的 $\delta$。 -
如何证明,$\R^2\setminus\set{(x,x^2)}$ 是开的?
提示
思路 1:通过闭集的充要条件证明 $\set{(x,x^2)}$ 是闭的;思路 2:对于上方的一个点 $(x_0,y_0)$,考虑 $f(x)=y_0-x^2$,它一定在 $x_0$ 附近 $>0$,就行了。思路 3:构造把抛物线拉平的连续函数 $\R^2\to\R^2$。 -
对于多元函数,无论按多少路径逼近是连续/极限相等,都不能说连续/极限存在。累次极限/大量路径极限存在 $\not\Rightarrow$ 连续;因定义域奇形怪状导致的:连续 $\not\Rightarrow$ 累次极限存在(如果存在一定等于)、在一个趋近于完整定义域的定义域子集内连续 $\not\Rightarrow$ 连续(2-11,$f(x,y)=x^y$)。
-
证明连续:用定义证;证明不连续:找到一条路径,极限不对。如何找?如果 $\to(0,0)$,就分析最低次项;如果不是,把 $x_0$、$y_0$ 减去后同理。注意还是要具体情况具体分析,多元函数中单项式的次数并不绝对说明增长速度,因为 $x$ 和 $y$ 之间可以有不同的关系(2-8,$xy/(x+y)$)。
-
无穷小量、渐进关系、同阶,这些概念都与范数的选择无关。wxf 讲的同阶不要求极限存在。(2-13)
-
证明:$x,y\to 0$ 时 $ax^2+2bxy+cy^2$ 与 $x^2+y^2$ 同阶,当且仅当 $b^2<ac$。(2-17)
提示
一个方向可以用三角换元。 -
双线性映射都是连续的。可以通过拆成基之和的方式来放缩,得到 $\|B(\b x,\b y)\|\le C\|(\b x,\b y)\|^2$,然后再利用双线性性。(2-28)
-
二元分式极限的求法:
- 极坐标。用于次数不是很大的情况。
- 换元 $y/x$。齐次的时候会比较好搞,不是齐次也不是不行。
- 不等式放缩。一般比较难,因为得考虑正负的问题。
- 渐进估计。用于出现 $\ln$、$\sin$、指数等的情况。
- 取 $y=f(x)$,把 $x$ 的成分减掉。这样可以证明极限不存在。
一定要注意,次数没有意义。哪怕分子的次数比分母大很多,极限仍有可能不存在——因为可以相减抵消,增高次数。
比较难的例子:$\lim_{(x,y)\to (0,0)}\frac{x^5+y^5}{x^3+y^3}$。这个用 1 可以搞但得分析三角函数的分式;用 2 的话会好处理。注意到 $\lim_{(x,y)\to (0,0)}\frac{x^p+y^p}{x^3+y^3}$ 在 $p=3,5,7,\cdots$ 时存在,$p=4,6,8,\cdots$ 时不存在。
-
定义矩阵的范数为 $\sqrt{\tr(A^\top A)}$ 或者 $\max\set{\lvert A_{i,j}\rvert}$ whatever,证明所有行列式 $>0$ 的实矩阵构成开集、道路连通集。
答案
行列式是 $\M_{n\times n}(\R)\to\R$ 的连续函数,这个证明和上面一样。于是就是开集。
道路连通集:思路 1(with psy):所有初等变换都可以“慢动作”做成连续的,所以可以把 $A$ 变成 $\diag(\pm 1)$ 的形式,然后可以用 $\rm SO(2)$ 里旋转的方式,将任意两个元素同时取负号。所有操作都是可逆的。思路 2(deepseek):把矩阵极分解啰,正定矩阵是凸集,行列式为 $1$ 的正交矩阵可以通过同上的旋转方法,严谨证明可以见 https://math.stackexchange.com/q/711492。 -
可逆矩阵虽然不是道路连通的,但是是稠密开的。考虑 $|A-\varepsilon I|$。
-
证明:设 $A$ 为正定矩阵,证明 $\b x^\top A\b x$ 与 $\b x^\top\b x$ 同阶。(2-32)
提示
考虑相减后利用正定的等价定义,相除单位化,或归为范数问题。 -
证明实对称矩阵可正交对角化。线性代数里的证明是利用代数基本定理搞出一个特征值,再证明这个特征值是实的(扩域的又一个应用)。
不通过这个来证明,能找到一个特征值。(3-13)
提示
对于所有 $\b x$ 满足 $\|\b x\|_2=1$,找到 $\max\set{\b x^\top A\b x}=\lambda$。利用反证法,证明 $A\b x^*-\lambda\b x^*=\b 0$。 -
写出压缩不动点的使用条件及证明。(3-18)
-
定义矩阵算子范数 $\|A\|_p=\max_{\|\b x\|_p=1}{\|A\b x\|_p}$,每种向量范数都可以导出一种算子范数,例如 $p=1$ 时,对应的算子范数求的是最大的【列的绝对值之和】。证明:https://zhuanlan.zhihu.com/p/30674132。算子范数满足 $\|AB\|\le\|A\|\|B\|$。
证明当 $\|B\|<1$ 时,$I-B$ 可逆;$\|(I-B)^{-1}-I\|\le\|B\|/(1-\|B\|)$;$\|(I-B)^{-1}-I-B\|\le\|B\|^2/(1-\|B\|)$。(3-21)
提示
用压缩不动点定理。先证唯一。
多元函数的微分与偏导
-
陈述公式 $\D\ip{,}(\b a,\b b)(\b x,\b y)=\ip{\b x,\b b}+\ip{\b a,\b y}$ 的含义。(4-8)
-
求 $\D\det(A)(B)$。(4-14)
答案
行列式是多重线性函数,展开。
答案为 $\tr({A^*}^\top B)$。 -
求 $\D(A^{-1})(B)$。(4-14)
-
回顾 $\b y^\prime=A(x)\b y$ 的解系构成的 $Y(x)$,它的行列式称为 Wronsky 行列式。其满足 $|Y(x)|^\prime=\tr A(x)|Y(x)|$。毕竟它说到底是个一元的函数,所以上个学期给出了证明,就是直接展开分析。
现在,通过 2 的答案,给出一个更高大上的推导。注意多元情况下,导数的除法极限定义会很不好搞,用渐进表示可能更容易。当然最终还是对一个变量 $x$ 求导。(4-11)
-
证明多元情况下,微分的链索法则 $\D(g\circ f)(\b x)=\D g(f(\b x))\circ\D f(\b x)$。(4-16)
注意这里要求 $g$ 和 $f$ 在对应点都可微,而不仅仅是偏导存在。
-
当 $f$ 是 $\to\R$ 时,$\D f$ 会写成 $\d f$。(4-3)
-
这里的链索法则,看似是“一元”的,实际上它是普适的。例如,我们考虑 $\d\ip{f(\b x),g(\b y)}$,它本质上是笛卡尔积上的映射(4-15) $$ (\b x,\b y)\xrightarrow{(f(\cdot_1),g(\cdot_2))}(f(\b x),g(\b y))\xrightarrow{\ip{\cdot_1,\cdot_2}}\ip{f(\b x),g(\b y)} $$ 所以 $$ \d\ip{f(\b x),g(\b y)}(\b u,\b v)=(\ip{f(\b x),\cdot_2}+\ip{\cdot_1,g(\b y)})\circ(\d f(\b x)(\cdot_1),\d g(\b y)(\cdot_2))(\b u,\b v)=\ip{f(\b x),\d g(\b y)(\b v)}+\ip{\d f(\b x)(\b u),g(\b y)} $$
-
如果 $\b x=\sum\xi_i\b e_i$,那么 $$ \d f(\b x_0)(\b x)=\sum\xi_i\d f(\b x_0)(\b e_i)=\sum\xi_i\frac{\p f}{\p\b e_i}(\b x_0)=\sum\frac{\p f}{\p x_i}(\b x_0)\d x_i(\b x) $$ 其中 $\d x_i$ 也是线性映射,即为取第 $i$ 维坐标的投影。因此 $$ \d f(\b x_0)=\sum\frac{\p f}{\p x_i}(\b x_0)\d x_k=\mat{\frac{\p f}{\p x_1}(\b x_0)&\cdots&\frac{\p f}{\p x_n}(\b x_0)}\mat{\d x_1\\ \vdots\\ \d x_n} $$ 这个东西和梯度没有任何关系,但是可以以矩阵形式处理复合运算——复合的微分,就是线性变换的复合。例如要求 $\D(g\circ f)(\b x_0)$,那么 $$ \begin{array}{ccccc} \D(g\circ f)(\b x_0)&=&\D g(f(\b x_0))&\circ&\D f(\b x_0)\\ \mat{\frac{\p(g\circ f)_1}{\p x_1}(\b x_0)&\cdots&\frac{\p(g\circ f)_1}{\p x_k}(\b x_0)\\ \vdots&\ddots&\vdots\\ \frac{\p(g\circ f)_n}{\p x_1}(\b x_0)&\cdots&\frac{\p(g\circ f)_n}{\p x_k}(\b x_0)}&=&\mat{\frac{\p g_1}{\p x_1}(f(\b x_0))&\cdots&\frac{\p g_1}{\p x_m}(f(\b x_0))\\ \vdots&\ddots&\vdots\\ \frac{\p g_n}{\p x_1}(f(\b x_0))&\cdots&\frac{\p g_n}{\p x_m}(f(\b x_0))}&\cdot&\mat{\frac{\p f_1}{\p x_1}(\b x_0)&\cdots&\frac{\p f_1}{\p x_k}(\b x_0)\\ \vdots&\ddots&\vdots\\ \frac{\p f_m}{\p x_1}(\b x_0)&\cdots&\frac{\p f_m}{\p x_k}(\b x_0)}\\ \J(g\circ f)(\b x_0)&=&\J(g)(f(\b x_0))&\cdot&\J(f)(\b x_0)\\ \displaystyle\frac{\p(g\circ f)_i}{\p x_j}(\b x_0)&=&\displaystyle\sum_k\frac{\p g_i}{\p t_k}(f(\b x_0))&\cdot&\displaystyle\frac{\p f_k}{\p x_j}(\b x_0) \end{array} $$
-
我们说线性泛函 $L$ 的梯度向量,是指向量 $\b b$ 使得 $\forall\b x$,$L(\b x)=\ip{\b b,\b x}$。现在对于一般的(不一定单位正交的)基底 $\b e_1,\cdots,\b e_n$,已知 $L(\b e_1),\cdots,L(\b e_n)$,不一定有 $\b b=L(\b e_1)\b e_1+\cdots+L(\b e_n)\b e_n$。
用 $\b e_i$ 和 $L(\b e_i)$ 表示出 $\b b$。(4-25、5-26)
上面那个式子一般而言并不是内积。仅在单位正交基时,内积等于行向量乘列向量。一般而言,设 $\b b=\sum c_i\b e_i$,根据内积的定义, $$ \ip{\b b,\b x}=\sum_{i,j}c_i\xi_j\ip{\b e_i,\b e_j}=\mat{c_1&\cdots&c_n}\mat{\ip{\b e_1,\b e_1}&\cdots&\ip{\b e_1,\b e_n}\\ \vdots&\ddots&\vdots\\ \ip{\b e_n,\b e_1}&\cdots&\ip{\b e_n,\b e_n}}\mat{\xi_1\\ \vdots\\ \xi_n} $$ 但另一方面,上面我们看到,根据微分的定义, $$ \d f(\b x_0)(\b x)=\mat{\frac{\p f}{\p x_1}(\b x_0)&\cdots&\frac{\p f}{\p x_n}(\b x_0)}\mat{\xi_1\\ \vdots\\ \xi_n} $$ 因此 $$ \b b=\mat{c_1&\cdots&c_n}^\top=\mat{\ip{\b e_1,\b e_1}&\cdots&\ip{\b e_1,\b e_n}\\ \vdots&\ddots&\vdots\\ \ip{\b e_n,\b e_1}&\cdots&\ip{\b e_n,\b e_n}}^{-1}\mat{\frac{\p f}{\p x_1}(\b x_0)\\ \vdots\\ \frac{\p f}{\p x_n}(\b x_0)} $$
-
利用梯度的定义证明 $\nabla(fg)=g\nabla f+f\nabla g$。(4-28)
-
$f$ 在某点的各方向导数都存在,哪怕各方向导数能统一写成内积的形式,也不意味着它在该点连续或可微。例如 $f(\rho,\theta)=[\rho\le\theta]$($\theta\in(0,2\pi]$),这个“平台”的分界不是“一致”的,在原点处无法取得一个全 $1$ 邻域。也可以定义路径导数,当可微时,路径 $\gamma$ 的导数 $\d f(\gamma(t))/\d t=\nabla f\cdot\gamma^\prime$(链索法则);一般情况下也行,如果要求方向导数(单位化),就取弧长参数。与连续性相同的是,可能出现各方向导数都一致,但沿某个曲线的方向导数爆炸,例如 $y/x^3(x=0\rightarrow 0)$。
求单位球面上的函数 $f$,沿经线与纬线分别的方向导数。(5-16)
-
区分沿向量的导数与方向导数,主要是是否单位化的区别。“一个函数沿坐标轴方向的方向导数就是该函数对这个坐标求的偏导数”这句话是错的,因为坐标的基向量可能不是单位的。用公式来表示,就是 $$ \frac{\p f}{\p\b{\hat e}_i}\ne\frac{\p f}{\p\b e_i}=\frac{\p f}{\p x_i} $$
-
Jacobi 矩阵行列记忆方法:线性泛函的微分是行向量,因此每行对应值域的一个维度。
-
利用复合函数的链索法则,解释多元函数的一阶微分不变性。(5-31)
-
对于一个函数 $f$,如果原来它的自变量为 $x,y$,然后换元后成了 $u,v$,那么求相关导数的思路是:将 $x,y$ 视作中间变量——也就是说,要求出 $x,y$ 关于 $u,v$ 的表达式,这样就可以求出 $f_u,f_v$ 等了。如果要将 $f_x,f_y$ 表示成 $u,v$ 相关的,那么要么用 $u,v$ 表示出 $x,y$,要么把第一个问题的东西反解出来(矩阵求逆)。
多元函数的高阶微分
-
$\mathscr{C}^k$ 类函数的定义是?(6-7)
-
证明矩阵求逆是 $\mathscr{C}^\infty$ 的。(6-8)
-
注意高阶偏导数“分母”上顺序与求导的顺序。(6-12)
-
简述以下命题的证明思路:
- $f$ 关于 $x$ 连续、$y$ 连续单调 $\Longrightarrow$ $f$ 连续。
- $f$ 在单点处 $x$ 偏导存在,邻域内 $y$ 偏导存在且连续 $\Longrightarrow$ $f$ 单点可导。(6-5)
- Clairaut 定理(6-33)
-
Clairaut 定理条件:等式两边的导数都是连续的,允许减弱到一侧连续。反例:https://math.stackexchange.com/a/4603178。
-
高阶微分: $$ \eq{ \d^kf(\b x_0)(\b u_1,\cdots,\b u_k)&=\left.\frac{\d}{\d t}\left(\d^{k-1}f(\b x_0+t\b u_k)(\b u_1,\cdots,\b u_{k-1})\right)\right|_{t=0}\\ &=\left.\frac{\d}{\d t}\sum_{i_1,\cdots,i_{k-1}}\frac{\p^{k-1}f(\b x_0+t\b u_k)}{\p x_{i_{k-1}}\cdots\p x_{i_1}}\b u_{1,i_1}\cdots\b u_{k-1,i_{k-1}}\right|_{t=0}\\ &=\sum_{i_1,\cdots,i_k}\frac{\p^kf(\b x_0)}{\p x_{i_k}\cdots\p x_{i_1}}\b u_{1,i_1}\cdots\b u_{k,i_k} } $$
-
推导微分方程 $u_{tt}=u_{xx}$ 的解。(6-17)
-
简述求拉普拉斯算子在极坐标下表述的思路。(6-25)
答案
我们要使用 $r$ 和 $\theta$ 表示 $u_x$ 和 $u_y$,这个东西如果直接用链式求导法则,会涉及到 $r$ 关于 $x$ 的偏导之类的,会比较麻烦。合适的方式是: $$ \frac{\p u}{\p(r,\theta)}=\frac{\p u}{\p(x,y)}\frac{\p(x,y)}{\p(r,\theta)} $$ 这样之后,矩阵求逆反解出 $\p u/\p(x,y)$ 即可。
然后会得到 $\p_x$ 和 $\p_y$ 被 $\p_r$ 和 $\p_\theta$ 表示的式子,再作用一次即可。对于 $\p_{yy}$,可以将坐标系关于 $x=y$ 对称一下,直接从 $\p_{xx}$ 的表达式推得。 -
偏导与积分交换的条件:求完偏导的函数中,求偏导的变量关于被积变量一致连续。(6-36)一个充分条件是,函数本身的偏导均连续。不可换序的反例:https://www.zhihu.com/question/661186049。
-
证明:多元函数 $f\in\mathscr{C}^r\Longrightarrow g=\int_tf(t)\d t\in\mathscr{C}^r$。(6-37)
-
不用定义证明变上下限积分的偏导。(6-39)
提示
考虑把积分写成一个多元(中间变量)函数。
注意说明链索法则的使用条件。 -
对于二重积分的换序,一个充分条件是被积函数是连续的。
证明这种情况下可以换序。(6-43)
提示
考虑利用上一个问题。注意偏导与积分可交换的条件。 -
积分换序条件:Fubini 定理。https://wuli.wiki/online/Fubin0.html。
-
泰勒展开推导。对于 $f\in\mathscr{C}^r$,我们考虑近似 $f(\b x_0+\b v)$。
根据一元泰勒展开, $$ \eq{ f(\b x_0+t\b v)&=\sum_{i=0}^{r-1}\frac{1}{i!}\left.\frac{\d^i}{\d t^i}f(\b x_0+t\b v)\right|_{t=0}t^i+\frac1{(r-1)!}\int_0^t\left.\frac{\d^r}{\d t^r}f(\b x_0+t\b v)\right|_{t=s}(t-s)^{r-1}\d s\\ f(\b x_0+\b v)&=\sum_{i=0}^{r-1}\frac{1}{i!}\left.\frac{\d^i}{\d t^i}f(\b x_0+t\b v)\right|_{t=0}+\frac{1}{(r-1)!}\int_0^1\left.\frac{\d^r}{\d t^r}f(\b x_0+t\b v)\right|_{t=s}(1-s)^{r-1}\d s\\ \text{(Integral)}&=\sum_{i=0}^{r-1}\frac{1}{i!}\sum_{j_1,\cdots,j_i}\frac{\p^i f}{\p x_{j_i}\cdots\p x_{j_1}}(\b x_0)\b v^{j_1}\cdots\b v^{j_i}+\frac{1}{(r-1)!}\int_0^1\sum_{j_1,\cdots,j_r}\frac{\p^r f}{\p x_{j_r}\cdots\p x_{j_1}}(\b x_0+s\b v)\b v^{j_1}\cdots\b v^{j_r}(1-s)^{r-1}\d s\\ &=\sum_{i=0}^{r-1}\frac{1}{i!}\sum_{j_1,\cdots,j_i}\frac{\p^i f}{\p x_{j_i}\cdots\p x_{j_1}}(\b x_0)\b v^{j_1}\cdots\b v^{j_i}+\frac{1}{(r-1)!}\sum_{j_1,\cdots,j_r}\frac{\p^r f}{\p x_{j_r}\cdots\p x_{j_1}}(\b x_0+\theta\b v)\b v^{j_1}\cdots\b v^{j_r}\int_0^1(1-s)^{r-1}\d s\\ \text{(Lagrange)}&=\sum_{i=0}^{r-1}\frac{1}{i!}\sum_{j_1,\cdots,j_i}\frac{\p^i f}{\p x_{j_i}\cdots\p x_{j_1}}(\b x_0)\b v^{j_1}\cdots\b v^{j_i}+\frac{1}{r!}\sum_{j_1,\cdots,j_r}\frac{\p^r f}{\p x_{j_r}\cdots\p x_{j_1}}(\b x_0+\theta\b v)\b v^{j_1}\cdots\b v^{j_r}\\ &=\sum_{i=0}^{r}\frac{1}{i!}\sum_{j_1,\cdots,j_i}\frac{\p^i f}{\p x_{j_i}\cdots\p x_{j_1}}(\b x_0)\b v^{j_1}\cdots\b v^{j_i}+\frac{1}{r!}\sum_{j_1,\cdots,j_r}\left(\frac{\p^r f}{\p x_{j_r}\cdots\p x_{j_1}}(\b x_0+\theta\b v)-\frac{\p^r f}{\p x_{j_r}\cdots\p x_{j_1}}(\b x_0)\right)\b v^{j_1}\cdots\b v^{j_r}\\ \text{(Peano)}&=\sum_{i=0}^{r-1}\frac{1}{i!}\sum_{j_1,\cdots,j_i}\frac{\p^i f}{\p x_{j_i}\cdots\p x_{j_1}}(\b x_0)\b v^{j_1}\cdots\b v^{j_i}+\omicron(\|\b v\|^r) } $$ 这里 Peano 余项是通过 $r$ 阶偏导的连续性得到(所以实际上 $r$ 阶可导的 Lagrange 余项,加上连续条件,就可以推出 $r+1$ 阶的 Peano 余项。这个上学期没提到)。
由于 $f\in\mathscr{C}^r$,故偏导的顺序无所谓。可以枚举每个自变量求了几次偏导,这样, $$ \frac{1}{i!}\sum_{j_1,\cdots,j_i}\frac{\p^if}{\p x_{j_1}\cdots\p x_{j_i}}(\b x_0)\b v^{j_1}\cdots\b v^{j_i}=\sum_{\alpha_1+\cdots+\alpha_k=i}\frac{\p^if}{\alpha_1!\cdots\alpha_k!\p^{\alpha_1}x_1\cdots\p^{\alpha_k}x_k}(\b x_0)(\b v^1)^{\alpha_1}\cdots(\b v^k)^{\alpha_k} $$ 方便计算。
证明:$f\in\mathscr{C}^r$ 在 Peano 余项形式的泰勒展开中的 $r$ 次多项式唯一。(7-10)从而有更方便的分多步展开方法。
提示
相减后通过一些简化的手段,变成单变量的情况。 -
默写二元函数极值点的分类讨论,并基于泰勒公式证明。(7-16)
-
说明为什么非退化驻点附近的等高线是椭圆或双曲线。
-
一个二元函数在某一点上,各个方向该点都是极小值点,该点不一定是极小值点。(7-23)
进一步地,甚至对于可微的函数,唯一的驻点是极值点,也不一定是最值点:https://math.stackexchange.com/a/3730294。
-
一些泰勒展开的套路:
- $f(x,y)=\left\{\begin{aligned}&\frac{1-\e^{x(x^2+y^2)}}{x^2+y^2},&(x,y)\ne(0,0)\\ &0,&(x,y)=(0,0)\end{aligned}\right.$
- $f(x,y)=\sin(xy)\e^{x+y}$
- $f(x,y)=\arctan\frac{1+x+y}{1-x+y}$
- $f(x,y)=x^y((x,y)\to(1,0))$
-
Hesse 矩阵在 $f$ 不二阶连续时,$(H_f)_{i,j}=f_{i,j}$,先对 $i$ 元求偏导。存疑
-
关于多元函数极值的推导(7-14~7-18):
-
从方向导数的角度描述驻点的必要条件,并写成代数语言。
-
为什么二阶部分不退化时,Hesse 矩阵的性质是极值点的充要条件?
-
写出极大、极小、鞍点的 Hesse 矩阵等价表述。
-
证明正定与顺序主子式 $>0$ 的等价性。
-
-
用数学语言说明:沿梯度反方向能走到极小值点。进一步说明为什么不太可能走到鞍点,或者说极小值是吸引不动点。进一步给出吸引的速度。(7-34、8-16、8-17)
-
共轭梯度法的思路:延续最速下降法的二次截断解极值点,然后只改方向——改为反伸缩时的正交方向。如果有多维,正交方向取负梯度被已有向量 Gram-Schmidt $H_f$-正交化后的结果。
-
说明牛顿迭代法是二阶收敛的。
-
再次强调凸是向下凸,凹是向上凸。
-
多元函数的凸性和任意线段上的凸性是等价的。
-
证明 $\mathscr{C}^2$ 函数在非退化极小值点的一个凸邻域内是严格凸函数。(9-9)
-
证明 $\mathscr{C}^2$ 函数是凸的 $\Longleftrightarrow$ $H$ 处处半正定;$H$ 处处正定 $\Longrightarrow$ 严格凸。(9-10)
-
证明凸函数的驻点是最小值点,以及严格凸 $\Longrightarrow$ 最小值点唯一。(9-12)
-
写出 Legendre 变换的定义,用 $f$、$\nabla f$ 及其反函数表示 $f^*$。(9-15~9-17)
-
Legendre 变换是对合的。
证明(deepseek)
Legendre 变换在多元情形下的对合性证明如下:设 $ f(\mathbf{x}) $ 为严格凸且二阶连续可导的 $ n $ 元函数,其 Legendre 变换定义为
$$ g(\mathbf{p}) = \sup_{\mathbf{x}} \left[ \mathbf{p} \cdot \mathbf{x} - f(\mathbf{x}) \right], $$ 其中 $\mathbf{p} = \nabla f(\mathbf{x})$ (梯度)。由严格凸性,Hessian 矩阵 $\nabla^2 f(\mathbf{x})$ 正定,故 $\mathbf{x}$ 与 $\mathbf{p}$ 一一对应,记 $\mathbf{x} = \mathbf{x}(\mathbf{p})$ ,则 $ g(\mathbf{p}) = \mathbf{p} \cdot \mathbf{x}(\mathbf{p}) - f(\mathbf{x}(\mathbf{p})) $ 。
对 $ g(\mathbf{p}) $ 再次应用 Legendre 变换,得到
$$ h(\mathbf{y}) = \sup_{\mathbf{p}} \left[ \mathbf{y} \cdot \mathbf{p} - g(\mathbf{p}) \right]. $$ 计算 $ g(\mathbf{p}) $ 的梯度:
$$ \nabla g(\mathbf{p}) = \mathbf{x}(\mathbf{p}) + \left[ \mathbf{p} \cdot \nabla_{\mathbf{p}} \mathbf{x}(\mathbf{p}) - \nabla f(\mathbf{x}(\mathbf{p})) \cdot \nabla_{\mathbf{p}} \mathbf{x}(\mathbf{p}) \right]. $$ 由于 $\mathbf{p} = \nabla f(\mathbf{x}(\mathbf{p}))$ ,第二项抵消,得 $\nabla g(\mathbf{p}) = \mathbf{x}(\mathbf{p})$ ,即 $\mathbf{y} = \nabla g(\mathbf{p}) = \mathbf{x}(\mathbf{p})$ ,故 $\mathbf{p} = \mathbf{p}(\mathbf{y})$ 为 $\mathbf{y} = \mathbf{x}(\mathbf{p})$ 的逆映射。
将 $ g(\mathbf{p}) = \mathbf{p} \cdot \mathbf{x}(\mathbf{p}) - f(\mathbf{x}(\mathbf{p})) $ 代入 $ h(\mathbf{y}) $ ,得
$$ h(\mathbf{y}) = \mathbf{y} \cdot \mathbf{p}(\mathbf{y}) - \left[ \mathbf{p}(\mathbf{y}) \cdot \mathbf{y} - f(\mathbf{y}) \right] = f(\mathbf{y}). $$ 因此,两次 Legendre 变换后恢复原函数,即 $ h(\mathbf{y}) = f(\mathbf{y}) $ ,证毕。 -
Young 不等式:$x^\alpha/\alpha+x^\beta/\beta\ge xy\;(x,y>0)$,$1/\alpha+1/\beta=1$,当且仅当 $x=y^{1/(\alpha-1)}$ 时取等。(9-19)
隐函数相关
-
隐函数定理:$\mathscr{C}^r\ni F:\R^m\times\R^n\to\R^n$ 写成 $F(\b x,\b y)$,$\b x\in\R^m$,$\b y\in\R^n$。已知 $F(\b x_0,\b y_0)=\b 0$ 且 $\J_{\b y}F(\b x_0,\b y_0)$ 可逆,那么存在 $\b x_0$ 的一个邻域 $U$ 及 $\b y_0$ 的一个邻域 $V$,存在 $\mathscr{C}^r\ni g:U\to V$,使得邻域内 $F(\b x,\b y)=0\Longleftrightarrow g(\b x)=\b y$。
证明:(10-9~10-20)
- 当 $n=1$ 时(记 $\b y$ 为 $y$,$F$ 为 $f$),存在 $U$、$V$,使得 $\b x\in U$ 时有唯一零点 $y(\b x)\in V$。
- 这样的 $y(\b x)$ 连续。
- 这样的 $y(\b x)$ 可微。
- $\J y(\b x)=\dfrac{\J_{\b x}f(\b x,y(\b x))}{f_y(\b x,y(\b x))}$,其中上下两个偏导都是对 $f$ 的变量求导而不是对整个复合函数求导。
- 当 $n>1$ 时结论也成立,且构造得到的 $G\in\mathscr{C}^r$。
提示
1. 利用 $\mathscr{C}^1$ 连续性,介值定理。只需 bound 在两个界之间。 2. 增加 $y$ 的邻域限制。 3. 先构造 4 中的这个结果,然后考虑直接证可微的定义。利用 $\mathscr{C}^1$,再进行一些移项,反解出 $y(\b x)-y_0$ 的界。 4. 链索法则。 5. 通过乘上 $\left[\J_{\b y}F(\b x_0,\b y_0)\right]^{-1}$,将 Jacobi 矩阵变成 $I_n$。这样一来可以直接确定 $y$ 的第 $n$ 个分量,依赖于 $\b x$ 以及前 $n-1$ 个分量。归纳。 -
注意隐函数定理的逆命题不一定成立。例如 $\sqrt[3]x$。
-
逆映射定理:$\mathscr{C}^r\ni H:\R^n\to\R^n$,若 $\b y_0=H(\b x_0)$ 且 $\J H(\b x_0)$ 可逆,则在 $\b x_0$ 的小邻域内,$\b x$ 可以写成关于 $\b y$ 的 $\mathscr{C}^r$ 函数。
证明隐函数定理和逆映射定理等价。(10-21~10-22)
-
整体微分同胚定理:$U$ 是 $\R^n$ 的开集。$\mathscr{C}^r\ni F:U\to\R^n$。$F$ 是同胚 $\Longleftrightarrow$ $F$ 单,且 $\D F$ 处处可逆。
这个定理的意思大概说的是,如果处处存在连续逆映射,那么整体就存在。
-
写出 $x^3+ax^2+bx-1=0$ 的解的近似表达式($a\to 0$ 且 $b\to 0$)。
-
写出 $n$ 维空间中 $m$ 维 $\mathscr{C}^r$ 曲面的定义。(11-6)$\set{(x,y,z)\mid x,y\ge 0\land z=x+y}$ 是二维曲面吗?$\set{(x,y,z)\mid |x|=|y|}$ 呢?
答案
不是,$x=0$ 处没法取邻域;不是,$x=y=0$ 处没法写成函数。 -
$\Sigma=\set{\b x\mid \b x_1^2+\cdots+\b x_n^2=1}\subset\R^n$ 是一个 $n-1$ 维球面。这个说法的证明是,每一处一定存在一个非零的 $\b x_i$,这个 $\b x_i$ 可以写成关于其他量的根式。
证明若 $F:\R^n\to\R$,$F^{-1}(0)$ 非空且不含 $F$ 的驻点,则它为一个 $n-1$ 维曲面。(11-10)推广到多个限制的情况。(11-11) -
写出并证明参数曲面相应的命题。(11-15)
-
注意以上所有命题都要在开集上——隐函数的 $F$、参数方程中的 $f_i$ 都要定义在开集上。
-
当 $\left(\sqrt{x^2+y^2}-R\right)^2+z^2=r^2$ 在 $r\ge R>1$ 时不再是个曲面。之前用隐函数和参数化对其的证明在哪里出了问题?
答案
隐函数:$F$ 在 $x=y=0$ 不可微;参数化:不是单射。 -
对于参数化曲面 $\b x(\b t)$,证明 $\span\J\b x(\b t_0)$ 恰好是所有过 $\b t_0$ 的曲线的切向量。思路是:每个向量都是切向量——构造 $\b t$ 各分量的线性组合(直线),然后代入 $\b x$;每个切向量都是向量——对于一个以另一个 $s$ 参数化的曲线,先把它写成一部分 $\b x$ 分量关于另一部分(一个)$\b x$ 分量的表达式,再根据后一部分 $\b x$ 和 $\b t$ 的可逆关系得出 $\b t$。现在就得到 $\b t$ 关于 $s$ 的表达式,再对 $\b x(\b t(s))$ 求导。
-
参数化曲面的切平面方程:$\b x=\b x(\b t_0)+\span\Set{\frac{\p\b x}{\p t_1}(\b t_0),\cdots,\frac{\p\b x}{\p t_n}(\b t_0)}$,即(可以理解为 $\b x$ 的一阶泰勒展开) $$ \b x-\b x_0=(\J\b x)\b t\quad\text{i.e.}\quad\left\{\begin{aligned}x-x_0&=\lambda\frac{\p x}{\p u}+\mu\frac{\p x}{\p v}+\cdots\\ y-y_0&=\lambda\frac{\p y}{\p u}+\mu\frac{\p y}{\p v}+\cdots\\ \vdots\end{aligned}\right. $$ 隐函数曲面的切平面方程:考虑任意一条曲线 $F(\b x(t))=0\Rightarrow\sum F_i\frac{\p\b x_i}{\p t}=0$,即 $\nabla F\cdot\Delta\b x=0$。对于多个限制(可以理解为 $F$ 的一阶泰勒展开等于零) $$ (\J F)(\b x-\b x_0)=0\quad\text{i.e.}\quad\left\{\begin{aligned}\frac{\p F_1}{\p x}(x-x_0)+\frac{\p F_1}{\p y}(y-y_0)+\cdots&=0\\ \frac{\p F_2}{\p x}(x-x_0)+\frac{\p F_2}{\p y}(y-y_0)+\cdots&=0\\ \vdots\end{aligned}\right. $$
-
条件极值($\min F(\b x)\text{ s.t. } \b g(\b x)=\b 0$)的思路:其实很简单。就是要求梯度向量在约束曲面上的投影为 $\b 0$,也就相当于梯度向量在约束曲面的切空间上投影为 $\b 0$,即 $\nabla F$ 可以表示为切空间的法向量,即 $\set{\nabla\b g_i}$ 的线性组合。即 $$ \frac{\p}{\p x_i}\left(F(\b x)-\sum_{j=1}^n\lambda_jg_j(\b x)\right)=0 $$ 对于定义域无限的情况,需要考虑无穷远处;否则需要考虑边界。一般处理都是说,当 $\|\b x\|>C$ 时 $F(\b x)$ 肯定过大。
-
条件极值的充分条件:$\mathcal{L}(\b x)=F(\b x)-\b\lambda\cdot\b g(\b x)$,在其驻点,如果固定 $\b\lambda$ 后关于 $\b x$ 的 $H_{\mathcal{L}}$ 在约束曲面的切空间上正定,那么 $\b x$ 就是极小值点。
证明(deepseek)
在多元条件极值问题中,若拉格朗日函数 $ L(x, \lambda) = F(x) - \lambda \cdot g(x) $ 在驻点处关于 $ x $ 的 Hessian 矩阵限制在约束曲面的切空间上正定,则 $ F(x) $ 在该点处取得局部极小值。以下是严格证明:
步骤 1:设定与隐函数定理的应用
设 $ x^* $ 是满足 $ g(x^*) = 0 $ 的可行点,且存在 $ \lambda^* $ 使得 $ \nabla_x L(x^*, \lambda^*) = 0 $。根据隐函数定理,若 $ g $ 关于后 $ n $ 个变量的 Jacobian 矩阵 $ J_v g(x^*) $ 非奇异,则存在局部参数化 $ v = G(u) $,其中 $ u \in \mathbb{R}^m $ 为前 $ m $ 个变量,$ v \in \mathbb{R}^n $ 为后 $ n $ 个变量,使得 $ g(u, G(u)) = 0 $。
步骤 2:参数化目标函数
将约束问题转化为无约束优化问题,定义参数化函数: $$ f(u) = F(u, G(u)). $$ 目标转化为证明 $ u^* $ 是 $ f(u) $ 的局部极小点,对应 $ x^* = (u^*, G(u^*)) $。
步骤 3:一阶条件
由隐函数定理,参数化后的梯度为: $$ \nabla f(u) = \nabla_u F + \nabla_v F \cdot DG(u), $$ 其中 $ DG(u) = -[J_v g]^{-1} J_u g $。在驻点 $ x^* $,由于 $ \nabla_x L = 0 $,即: $$ \nabla F(x^*) = \lambda^{*T} J_g(x^*), $$ 代入 $ J_g = [J_u g, J_v g] $,可得 $ \nabla_u F = \lambda^{*T} J_u g $ 和 $ \nabla_v F = \lambda^{*T} J_v g $。由此验证: $$ \nabla f(u^*) = \lambda^{*T} J_u g + \lambda^{*T} J_v g \cdot (-[J_v g]^{-1} J_u g) = 0, $$ 即 $ u^* $ 是 $ f(u) $ 的临界点。
步骤 4:二阶条件与 Hessian 分析
计算 $ f(u) $ 的二阶导数: $$ H_f(u^*) = \nabla_{uu}^2 F + \nabla_{uv}^2 F \cdot DG + DG^T \cdot \nabla_{vu}^2 F + DG^T \cdot \nabla_{vv}^2 F \cdot DG + \nabla_v F \cdot D^2 G(u^*). $$ 其中 $ D^2 G(u^*) $ 由隐函数定理的二阶展开确定: $$ D^2 G = -[J_v g]^{-1} \left( \nabla_{uu}^2 g + 2 \nabla_{uv}^2 g \cdot DG + DG^T \cdot \nabla_{vv}^2 g \cdot DG \right). $$ 代入 $ \nabla_v F = \lambda^{*T} J_v g $,可得: $$ \nabla_v F \cdot D^2 G = -\lambda^{*T} \left( \nabla_{uu}^2 g + 2 \nabla_{uv}^2 g \cdot DG + DG^T \cdot \nabla_{vv}^2 g \cdot DG \right). $$ 因此,$ H_f(u^*) $ 可整理为: $$ H_f(u^*) = \left[ \nabla_{xx}^2 F - \lambda^{*T} \nabla_{xx}^2 g \right] \bigg|_{\text{切空间}}. $$ 这恰好是拉格朗日函数 $ L $ 的 Hessian $ \nabla_{xx}^2 L(x^*, \lambda^*) $ 在切空间上的投影。
步骤 5:正定性结论
若 $ \nabla_{xx}^2 L(x^*, \lambda^*) $ 在切空间上正定,则 $ H_f(u^*) $ 正定,故 $ u^* $ 是 $ f(u) $ 的严格局部极小点,对应 $ x^* $ 为 $ F(x) $ 在约束 $ g(x) = 0 $ 下的严格局部极小点。
结论
当拉格朗日函数 $ L(x, \lambda) $ 在驻点处的 Hessian 矩阵限制于约束曲面的切空间上正定时,原函数 $ F(x) $ 在该点取得局部极小值。这由二阶充分条件保证,其核心在于参数化后的无约束问题中,Hessian 的正定性与拉格朗日 Hessian 的投影正定性一致。
(详细步骤解释)在约束优化问题中,关键是将带约束的优化转化为无约束优化,并分析参数化后函数的二阶导数。下面详细解释为什么参数化后的Hessian矩阵 $ H_f(u^*) $ 等于拉格朗日函数Hessian在切空间上的投影。
1. 参数化函数的二阶导数
设约束曲面 $ g(x) = 0 $ 在 $ x^* = (u^*, v^*) $ 附近局部可表示为 $ v = G(u) $,其中 $ u \in \mathbb{R}^m $,$ v \in \mathbb{R}^n $。参数化目标函数为: $$ f(u) = F(u, G(u)). $$ 其Hessian矩阵 $ H_f(u) $ 的完整表达式为: $$ H_f(u) = \underbrace{\nabla_{uu}^2 F}_{\text{直接项}} + \underbrace{\nabla_{uv}^2 F \cdot DG + DG^T \cdot \nabla_{vu}^2 F}_{\text{交叉项}} + \underbrace{DG^T \cdot \nabla_{vv}^2 F \cdot DG}_{\text{二次项}} + \underbrace{\nabla_v F \cdot D^2 G}_{\text{隐函数二阶项}}. $$ 其中 $ DG = \frac{\partial G}{\partial u} $,$ D^2 G $ 为 $ G $ 的二阶导数。
2. 隐函数定理的二阶导数公式
由于 $ g(u, G(u)) = 0 $,对两边求导可得: $$ \nabla_u g + \nabla_v g \cdot DG = 0 \implies DG = -[\nabla_v g]^{-1} \nabla_u g. $$ 进一步对 $ u $ 求导,得到 $ D^2 G $: $$ D^2 G = -[\nabla_v g]^{-1} \left( \nabla_{uu}^2 g + \nabla_{uv}^2 g \cdot DG + DG^T \cdot \nabla_{vu}^2 g + DG^T \cdot \nabla_{vv}^2 g \cdot DG \right). $$ 注意这里用到了链式法则和隐函数定理的二阶展开。
3. 代入拉格朗日条件
在驻点 $ x^* $,一阶条件为: $$ \nabla_x L = \nabla F - \lambda^* \nabla g = 0 \implies \nabla F = \lambda^* \nabla g. $$ 具体分块为: $$ \nabla_u F = \lambda^* \nabla_u g, \quad \nabla_v F = \lambda^* \nabla_v g. $$ 利用这一关系,我们可以将参数化Hessian中的 $ \nabla_v F \cdot D^2 G $ 项替换为 $ \lambda^* \nabla_v g \cdot D^2 G $。
4. 显式展开参数化Hessian
将 $ D^2 G $ 的表达式代入 $ H_f(u^*) $,并利用一阶条件: $$ \nabla_v F \cdot D^2 G = -\lambda^* \left( \nabla_{uu}^2 g + 2 \nabla_{uv}^2 g \cdot DG + DG^T \cdot \nabla_{vv}^2 g \cdot DG \right). $$ 因此,参数化Hessian可写为: $$ H_f(u^*) = \nabla_{uu}^2 F + 2 \nabla_{uv}^2 F \cdot DG + DG^T \cdot \nabla_{vv}^2 F \cdot DG - \lambda^* \left( \nabla_{uu}^2 g + 2 \nabla_{uv}^2 g \cdot DG + DG^T \cdot \nabla_{vv}^2 g \cdot DG \right). $$ 观察到: $$ H_f(u^*) = \left[ \nabla_{uu}^2 F + 2 \nabla_{uv}^2 F \cdot DG + DG^T \cdot \nabla_{vv}^2 F \cdot DG \right] - \lambda^* \left[ \nabla_{uu}^2 g + 2 \nabla_{uv}^2 g \cdot DG + DG^T \cdot \nabla_{vv}^2 g \cdot DG \right]. $$ 这正是: $$ H_f(u^*) = \left[ \nabla_{xx}^2 F - \lambda^* \nabla_{xx}^2 g \right] \bigg|_{(u, G(u))} \text{在切空间上的投影}. $$
5. 切空间与Hessian投影
- 切空间定义:约束曲面 $ g(x) = 0 $ 在 $ x^* $ 处的切空间由满足 $ \nabla g(x^*) \cdot h = 0 $ 的向量 $ h \in \mathbb{R}^{m+n} $ 构成。
- 参数化切向量:在参数化 $ x = (u, G(u)) $ 下,切向量可表示为 $ h = (du, DG \cdot du) $,其中 $ du \in \mathbb{R}^m $。
- Hessian的作用:对任意切向量 $ h $,有: $$ h^T \cdot \nabla_{xx}^2 L \cdot h = h^T \cdot (\nabla_{xx}^2 F - \lambda^* \nabla_{xx}^2 g) \cdot h. $$ 这正是参数化Hessian $ H_f(u^*) $ 的作用方式。
6. 关键结论
参数化后的Hessian $ H_f(u^*) $ 与拉格朗日Hessian $ \nabla_{xx}^2 L $ 在切空间上的投影完全一致。因此:
- 若 $ \nabla_{xx}^2 L $ 在切空间上正定,则 $ H_f(u^*) $ 正定。 - 由无约束优化的二阶充分条件,$ u^* $ 是 $ f(u) $ 的严格局部极小点,对应 $ x^* $ 为原问题的严格局部极小点。
总结图示 $$ \boxed{ \begin{aligned} &\text{参数化Hessian } H_f(u^*) \\\ &= \text{拉格朗日Hessian } \nabla_{xx}^2 L \text{ 在切空间上的投影} \\\ &= \left[ \nabla_{xx}^2 F - \lambda^* \nabla_{xx}^2 g \right] \bigg|_{\text{切向量 } h = (du, DG \cdot du)}. \end{aligned} } $$这里有另一个问题是,既然已知“$\mathcal{L}$ 在某个驻点处关于 $\b x$ 的 Hessian 局限在该点约束曲面的切空间上正定,则对应 $F(\b x)$ 是约束曲面上的极小值点”,那么是否有“$F(\b x)$ 在 $\cal L$ 驻点对应的 $\b x$ 处的 Hessian 局限在该点约束曲面的切空间上正定,则对应 $F(\b x)$ 是约束曲面上的极小值点”?
答案应该是否定的,具体反例不太会构造。
广义含参积分
-
几个换序定理的条件(继累次极限交换以及 Clairaut):
-
$f$ 连续,$g(\b x)=\int_a^bf(t,\b x)\d t$,则 $g$ 连续。
-
$f_{x_i}$ 连续,$g(\b x)=\int_a^bf(t,\b x)\d t$,则 $g_{x_i}(\b x)=\int_a^bf_{x_i}(t,\b x)\d t$ 且连续。
-
$f$ 连续,则 $\int_c^d\int_a^bf(t,x)\d t\d x=\int_a^b\int_c^df(t,x)\d x\d t$。
-
$f$ 连续,$g(\b x)=\int_a^\omega f(t,\b x)\d t$ 一致收敛,则 $g$ 连续。
用连续的定义,将 $\int_a^\omega f(t,\b x)\d t-\int_a^\omega f(t,\b x_0)\d t$ 拆成 $\int_a^b(f(t,\b x)-f(t,\b x_0))\d t$、$\int_b^\omega f(t,\b x_0)\d t$、$\int_b^\omega f(t,\b x)\d t$ 三部分放缩。
-
$f$ 连续,$\int_a^\omega f(t,x)\d t$ 一致收敛,则 $\int_c^d\int_a^\omega f(t,x)\d t\d x=\int_a^\omega\int_c^df(t,x)\d x\d t$。
$\int_c^d\int_a^\omega\gets\int_c^d\int_a^b=\int_a^b\int_c^d\to\int_a^\omega\int_c^d$。
-
$f_{x_i}$ 连续,$g(\b x)=\int_a^\omega f(t,\b x)\d t$(某单点)收敛,$\int_a^\omega f_{x_i}(t,\b x)\d t$ 一致收敛,则 $g_{x_i}(\b x)=\int_a^\omega f_{x_i}(t,\b x)\d t$ 且连续。
用定义拆导数,$f(t,\b x+s\b e_i)-f(t,\b x)$ 写成积分,用 5 换序。
-
-
叙述广义积分的一致 Cauchy 条件,并证明 Weierstrass 判别法。(14-20~14-22)
-
说明含参积分不一致收敛的几个方法:定义的否命题、一致 Cauchy 条件的否命题、一致收敛 $\Rightarrow$ 连续的逆否命题。
-
关于几个一致收敛判别法的对比(严格证明:https://zhuanlan.zhihu.com/p/700675944):
展开
1. Cauchy 准则(引理)
定理:$\int_a^{+\infty} f(x) \d x $ 收敛 $\Longleftrightarrow$ $\forall \epsilon > 0$,$\exists N > a $,$\forall b, c > N $,$\left| \int_b^c f(x) \d x \right| < \epsilon $。
一致版本:$\int_a^{+\infty} f(x, \theta) \d x $ 对 $\theta$ 一致收敛 $\Longleftrightarrow$ $\forall \epsilon > 0$,$\exists N > a $,$\forall b, c > N $,$\forall\theta$,$\left| \int_b^c f(x, \theta) \d x \right| < \epsilon $。
2. 比较判别法。
定理:若 $ |f(x)| \leq M(x) $ 且 $\int_a^{+\infty} M(x) \d x $ 收敛,则 $\int_a^{+\infty} f(x) \d x $ 绝对收敛。
证明:对任意 $ b > a $,有 $\int_a^b |f(x)| \d x \leq \int_a^b M(x) \d x \leq \int_a^{+\infty} M(x) \d x $,故 $\int_a^{+\infty} |f(x)| \d x $ 收敛。
一致版本(Weierstrass 判别法):若 $ |f(x, \theta)| \leq M(x,\theta) $ 且 $\int_a^{+\infty} M(x,\theta)\d x $ 对 $\theta$ 一致收敛,则 $\int_a^{+\infty} f(x, \theta)\d x $ 对 $\theta$ 一致收敛。
证明:由一致收敛的 Cauchy 准则,$\forall\epsilon > 0$,$\exists N > a $,$\forall b, c > N $,$\forall \theta$,$\int_b^c M(x,\theta) \d x < \epsilon $。从而 $\left| \int_b^c f(x, \theta) \d x \right| \leq \int_b^c |f(x, \theta)| \d x \leq \int_b^c M(x,\theta) \d x < \epsilon $,与 $\theta$ 无关。
3. Dirichlet 判别法
定理:设 $ f(x) $ 和 $ g(x) $ 满足 $ F(b) = \int_a^b f(x) \d x $ 有界,$g(x)$ 单调且 $\to 0$,则 $\int_a^{+\infty} f(x)g(x) \d x $ 收敛。
证明:分部积分得:$\int_a^b f(x)g(x) \d x = F(b)g(b) - F(a)g(a)-\int_a^b F(x)g^\prime(x) \d x $。当 $ b \to +\infty $,第一项 $ F(b)g(b) \to 0 $(因 $ g(b) \to 0 $),第二项因 $ |F(x)| \leq C $ 且 $ \int_a^{+\infty} g^\prime(x) \d x=-g(a)$,故绝对收敛。
一致版本:设 $ f(x, \theta) $ 和 $ g(x, \theta) $ 满足 $F(b,\theta)=\int_a^b f(x, \theta) \d x $ 对 $\theta$ 一致有界,$ g(x, \theta) $ 关于 $ x $ 单调且 $ g(x, \theta) \to 0 $ 对 $\theta$ 一致收敛,则 $\int_a^{+\infty} f(x, \theta)g(x, \theta) \d x $ 一致收敛。
证明:$F(b,\theta)g(b,\theta)$ 是一致有界乘一致收敛,所以一致收敛;$\int_a^bF(x,\theta)g^\prime(x,\theta)\d x$ 把 $F$ 提出来也是一致有界乘一致收敛。
4. Abel 判别法
定理:设 $ f(x) $ 和 $ g(x) $ 满足 $\int_a^{+\infty} f(x) \d x $ 收敛,$ g(x) $ 单调有界,则 $\int_a^{+\infty} f(x)g(x) \d x $ 收敛。
证明:$\int_a^{+\infty}f(x)g(x)\d x=g(+\infty)\int_a^{+\infty}f(x)\d x+\int_a^{+\infty}f(x)(g(x)-g(+\infty))\d x$,第一部分收敛,第二部分由 Dirichlet 判别法收敛。
一致版本:设 $ f(x, \theta) $ 和 $ g(x, \theta) $ 满足 $\int_a^{+\infty} f(x, \theta) \d x $ 对 $\theta$ 一致收敛,$ g(x, \theta) $ 关于 $ x $ 单调且对 $\theta$ 一致有界,则 $\int_a^{+\infty} f(x, \theta)g(x, \theta) \d x $ 对 $\theta$ 一致收敛。
证明:$g(+\infty,\theta)\int_a^{+\infty}f(x,\theta)\d x$ 是一致有界乘一致收敛,后面是一样的。
重积分
-
高维矩形 $I$ 上的 Lebesgue 准则:$f$ Riemann 可积 $\Longleftrightarrow$ $f$ 有界且 Darboux 可积 $\Longleftrightarrow$ $f$ 有界且间断点测度为零。
写出高维情况下,零测度的定义。(15-12)
-
证明:对 $I=[a,b]\times[\alpha,\beta]$ 上的连续函数 $f$,(15-13) $$ \iint_I f(x,y)\d x\d y=\int_a^b\int_\alpha^\beta f(x,y)\d y\d x=\int_\alpha^\beta\int_a^b f(x,y)\d x\d y $$
提示
将重积分写成 $$\lim\sum_{i,j}f(\xi_{i,j},\eta_{i,j})\Delta x_i\Delta y_j$$ 的形式,然后利用 $\sum$ 的交换律得到后两种积分。一般情况(Fubini 定理):https://zhuanlan.zhihu.com/p/695338328、https://zhuanlan.zhihu.com/p/668687961。这里值得注意的点是,广义积分的情况要求绝对可积。
-
证明:$\p A$ 是闭的。(15-23)
-
对于一般区域 $A$ 上的积分,是通过乘上示性函数,转化为矩形上的积分来定义的。如果有界 $f:A\to\R$ 的间断点为零测集,且 $\p A$ 为零测集,则 $\int_A f\d x$ 有定义。(15-24)
-
通过换元的方式,将以下积分变成 $[0,1]^3$上的积分:(15-34) $$ \iiint_{x,y,z\ge 0,\|(x,y,z)\|\le 1}xyz\d x\d y\d z $$
-
理一下重积分换元的逻辑:假设有一个映射(微分同胚)$f:U\to V$,那么一个小区域 $I$ 经过 $f$ 映射后,体积变为 $\lvert\det\J f\rvert$ 倍。这意味着,在 $V$ 中的区域的积分 $\int_V w(\b x)\d v$,如果要换到 $U$ 中,就需要逐点乘上这个伸缩倍数,也就是 $\int_U w(f(\b t))\lvert\det\J f(\b t)\rvert\d v$。换句话说,这个换元的逻辑关系是逆过来的,或者说可以记忆为,重积分换元的标准形式是第二类换元法,其中 $f^\prime(t)$ 变为了 $\lvert\det\J\rvert$。如果要用第一类换元法,就是把同胚逆过来,用第二类换元法。
给出以下重积分的换元,使得区域变为矩形:(16-21) $$ \iint_{(x^2-y^2,2xy)\in[a,b]\times[\alpha,\beta]}f(x,y)\d x\d y $$
-
求 $\iint_D(x^2+y^2+1)^{-2}\d x\d y$,其中 $D$ 为圆心在 $(1,1)$ 且经过原点的圆盘。(16-26)
-
计算 $\int_{[0,1]^m}\min\set{x_1,\cdots,x_n}\d x_1\cdots\d x_n$。
曲线曲面积分
-
证明曲线积分与参数化方法无关。即 $\gamma$ 的两种参数化方法 $\b x(s)$、$\b x(t)$,有(17-5) $$ \int_\gamma f(\b x)\d l=\int_a^bf(\b x(s))\|\b x^\prime(s)\|\d s=\int_\alpha^\beta f(\b x(t))\|\b x^\prime(t)\|\d t $$
-
根据叉积公式写出二维曲面的面积公式。(17-8)证明与参数选择无关。(17-11)
-
证明 $n$ 维空间中的 $k$ 维体的体积,可以由它 $k$ 条边的向量的度量矩阵开根算出。(17-12~17-15)
提示
考虑 Gram-Schmidt 正交化,归纳法。 -
求 $n$ 维超曲面 $x^{n+1}=f(x^1,\cdots,x^n)$ 的面积公式,利用 3。(17-19~17-20)
-
利用上面的方法,严格说明三维旋转体的表面积公式 $\int 2\pi y\sqrt{1+(y^\prime)^2}\d x$。(17-26~17-27)
-
说明曲面维数和背景空间维数相同时,积分系数即为换元系数。(16-16)
-
总结:第一类曲线/曲面积分: $$ \int_Df(\b x)\d v=\int_Df(\b x)\sqrt{\det\left(\ip{\frac{\p\b x}{\p t_i},\frac{\p\b x}{\p t_j}}\right)}\d t_1\cdots\d t_k $$
-
说明做功、环量、通量都可以写成以下形式:(18-7~18-9) $$ \int_\gamma\b F\cdot\d\b l=\int_a^b\b F(\b x(t))\cdot\b x^\prime(t)\d t=\int_a^b\sum_i\b F^i(\b x(t))\d\b x^i(t) $$
-
一阶微分形式的定义?一阶微分形式的积分的定义?(18-12)
答案
一阶微分形式 $\omega$ 是一个从向量场到函数的映射。在区域(严格来说是微分流形)$R$ 上,$\omega$ 将一个定义在 $R$ 上各点的向量函数(可以理解为每个点上的切向量)映射到定义在 $R$ 上各个点的实函数,这里对于单个点来说,向量到实数的映射是线性的。如果要取向量场 $\b v$ 经过 $\omega$ 后某一点 $\b x$ 处的结果,写作 $\omega(\b v)(\b x)$ 或 $\omega(\b x)\b v$。
$\int_\gamma\omega$ 定义为 $\int_a^b\omega(\b x(t))(\b x^\prime(t))\d t$,其中 $\b x(t)$ 是 $\gamma$ 的参数化表示。
值得注意的是,对于 $f:\R^n\to\R$,$\D f$ 也是个一阶微分形式,但 $\D f(\b x_0)$ 只是一个线性泛函。 -
说明 $\int_\gamma\omega$ 与(保方向)参数化表达选取无关。(18-13)
-
保守场的定义?(18-18)无旋场的定义?(18-24)
- 证明梯度场是保守场。(18-16)
- 证明保守场是梯度场。(18-19)
- 证明保守场是无旋场。(18-23)
- 说明什么时候无旋场不是保守场。(18-31)
-
求第二类曲线积分的方法是:首先先尝试凑全微分。剩下凑不了的部分,参数化之后暴力积。
注意凑出全微分后还要保证区域单连通。
求 $\int_\gamma(z\d x+x\d y+y\d z)$,其中 $\gamma$ 为 $\set{x^2+y^2+z^2=R^2}\cap\set{x+y+z=0}$ 在 $xy$ 平面上的投影,逆时针。(18-28)
-
证明 $\int_\gamma x\d y=-\int_\gamma y\d x=\frac12\int_\gamma(x\d y-y\d x)$。(18-29)
-
拉回的概念:对于 $\Phi:A(\subseteq\R^m)\to B(\subseteq\R^n)$,它的微分 $\D\Phi$(记为 $\Phi_*$)在 $A$ 的每一个点上定义了一个 $\R^m\to\R^n$ 的线性映射(或者也可以理解为将 $A$ 上的向量场映到 $B$ 上的向量场),就是切向量的线性映射。$\Phi$ 的拉回 $\Phi^*$ 定义了 $\Omega^1(B)\to\Omega^1(A)$ 的线性映射,这个线性映射将 $B$ 上的一阶微分形式 $\omega$ 映射成这样的一个 $\R^m$ 上的微分形式:拿到一个向量场,将它先依据 $\Phi_*$ 映到 $\R^m$ 上,再应用 $\omega$。直观地说,就是 $\omega$ 对应 $B$ 上积分时的“权重函数”,$\Phi^*(\omega)$ 就是 $A$ 上积分时对应的权重。它首先将 $A$ 上切向量微元对应到 $B$ 上,再应用原来的权重。用数学语言来说就是,对于一个 $A$ 上的向量场 $\b v$, $$ \Phi^*(\omega)(\b v)(\b x)=\omega(\Phi_*(\b v))(\b x) $$ (这里 $\Phi_*$ 用的是括号里的定义,就是向量场之间的映射,注意笔记和课件的其他地方一般不会用这种方式看待微分,所以微分后面的括号里一般写的是个位置而不是向量场)
拉回严格定义了第二类曲线积分的换元方式: $$ \int_{\Phi(\gamma)}\omega=\int_\gamma\Phi^*(\omega) $$ 证明: $$ \begin{align*} \int_\gamma\Phi^*(\omega)&=\int_a^b\Phi^*(\omega)(\b x^\prime)(\b x(t))\d t\\ &=\int_a^b\omega(\Phi_*(\b x^\prime))(\b x(t))\d t\\ &=\int_a^b\sum_i F^i(\Phi(\b x(t)))(\D\Phi)(\b x^\prime)(\b x(t))^i\d t\\ &=\int_a^b\sum_i F^i(\b y(t))\b y^\prime(t)^i\d t\\ &=\int_a^b\omega(\b y^\prime)(\b y(t))\d t\\ &=\int_{\Phi(\gamma)}\omega \end{align*} $$ 其中 $\b y(t)=\Phi(\b x(t))$,而 $\b y^\prime(t)=\D\Phi(\b x(t))\cdot\b x^\prime(t)$,如果把 $\D\Phi$ 理解成向量场之间的映射,就是 $\D\Phi(\b x^\prime)(\b x(t))$。
对于 18-39~18-45 的例子,$\Phi$ 是从 $(\xi,\eta,\zeta)$ 到 $(x,y,z)$ 的映射,为了计算 $\Phi^*(z\d x+x\d y+y\d z)$,把 $\xi,\eta,\zeta$ 代入展开即可。
-
证明一个 $\R^2\to\R^2$ 的线性函数($\R^2$ 上的 $\R^2$ 向量场)可以拆成一个无旋场和一个无源场之和。(19-5)
-
证明:对于一个矩形 $R$ 而言,$\int_{\p R}X\d x+Y\d y=\iint_R(Y_x-X_y)\d x\d y$。(19-9)
-
同理可以证明 $\int_{\p R}(-Y\d x+X\d y)=\iint_R(X_x+Y_y)\d x\d y$。写出由这两个式子给出的旋度和散度的另一种表达式。
-
二阶微分形式将两个向量场映到一个标量场。写出斜积的定义。(19-21)
-
超曲面定义为维数比背景空间维数小 $1$ 的曲面。对于隐函数 $F(\b x)=0$($0$ 是 $F$ 的正则值),它定义了一个可定向的超曲面。法向量场为 $\nabla F$ 单位化的结果。利用这个结果证明:
-
可微函数的图像是一个可定向曲面。
答案
$y=f(\b x)$ 可以写成 $y-f(\b x)=0$,它的梯度向量为 $\mat{-\nabla f\\\ 1}$。 -
正则参数曲面是一个可定向曲面。回顾正则指的是单射且 $\p\b x/\p\b t$ 列满秩。
答案
构造函数 $L(\b v)=\dat{\b v&\frac{\p\b x}{\p\b t}}$,我们声称 $\nabla L$ 就是法向量。这是因为,如果代入某个 $\frac{\p\b x}{\p t_i}=\b v$,$L(\b v)$ 就等于 $0$ 了,在笛卡尔坐标系下,这就说明 $\nabla L$ 与所有 $\frac{\p\b x}{\p t_i}$ 内积为 $0$,即垂直。$\nabla L$ 是线性函数,从而是连续的。
-
-
计算高维超曲面通量的方法:首先记 $\b A=\nabla L$,面积微元在数值上等于单位高度的体积,也就是 $$ L(\b{\hat A})=\dat{\frac{\b A_1}{|\b A|}&\frac{\p\b x_1}{\p t_1}&\cdots&\frac{\p\b x_1}{\p t_m}\\ \vdots&\vdots&\ddots&\vdots\\ \frac{\b A_{m+1}}{|\b A|}&\frac{\p\b x_{m+1}}{\p t_1}&\cdots&\frac{\p\b x_{m+1}}{\p t_m}}=\sum_{i=1}^{m+1}\frac{\b A_i^2}{|\b A|}=|\b A| $$
$$ \begin{align*} \int_\Omega\b V\cdot\b n\d S&=\int_T\b V\cdot\frac{\b A}{|\b A|}(|\b A|\d t_1\cdots\d t_m)\\ &=\int_T\b V\cdot\b A\d t_1\cdots\d t_m\\ &=\int_T\dat{\b V_1&\frac{\p\b x_1}{\p t_1}&\cdots&\frac{\p\b x_1}{\p t_m}\\ \vdots&\vdots&\ddots&\vdots\\ \b V_1&\frac{\p\b x_{m+1}}{\p t_1}&\cdots&\frac{\p\b x_{m+1}}{\p t_m}}\d t_1\cdots\d t_m \end{align*} $$
其中 $T$ 为参数的定义域。证明这个计算与参数化表达的选择无关。(20-17)
-
三维中二阶微分形式的定义:
对于一个二阶微分形式 $\omega=X\d y\w\d z+Y\d z\w\d x+Z\d x\w\d y$,对于定向曲面 $\Sigma$ 上 $\omega$ 的积分,这里通过参数化表达来定义。设参数化表达为 $(x,y,z)=(x(u,v),y(u,v),z(u,v))$,其中通过 $\nabla\det\left(\b v,\frac{\p(x,y,z)}{\p(u,v)}\right)$ 得到的法向量和 $\Sigma$ 的定向一致(否则需要加负号),那么 $$ \iint_\Sigma\omega=\iint_D\left(X\det\frac{\p(y,z)}{\p(u,v)}+Y\det\frac{\p(z,x)}{\p(u,v)}+Z\det\frac{\p(x,y)}{\p(u,v)}\right)\d u\d v $$ 如果觉得难理解的话,可以以换元的角度来看。例如 $$ \d y\w\d z(u,v)=\det\frac{\p(y,z)}{\p(u,v)} $$ 以及 $$ \d y\w\d z=(y_u\d u+y_v\d v)\w(z_u\d u+z_v\d v)=(y_uz_v-y_vz_u)\d u\w\d v $$ 而在一致定向的情况下,$\iint\d u\w\d v=\iint\d u\d v$,也就是说,我们将微分形式的积分转换成黎曼积分。
二维曲面通量积分与二阶微分形式积分的关系:注意到 $$ \iint_\Sigma\b V\cdot\b n\d S=\iint_D\dat{\b V_x&x_u&x_v\\ \b V_y&y_u&y_v\\ \b V_z&z_u&z_v}\d u\d v=\iint_\Sigma\omega $$
-
证明:任意微分形式 $\omega$,$\d(\d\omega)=0$。
答案
$$ \begin{align*} \d\omega&=\sum_{i_1,\cdots,i_m}\left(\sum_{i=1}^n\frac{\p F_{i_1,\cdots,i_m}}{\p x_i}\d x_i\right)\w\d x_{i_1}\w\cdots\w\d x_{i_m}\\\ \d(\d\omega)&=\sum_{i_1,\cdots,i_m}\left(\sum_{i=1}^n\sum_{j=1}^n\frac{\p F_{i_1,\cdots,i_m}}{\p x_j\p x_i}\d x_j\w\d x_i\right)\w\d x_{i_1}\w\cdots\w\d x_{i_m} \end{align*} $$ 其中 $x_i$ 和 $x_j$ 可互换,正负抵消;$i=j$ 的时候 $\d x_j\w\d x_i=0$。顺带一提,对于 $n-1$ 阶微分形式,它的外微分是 $\tr\frac{\p\b F}{\p\b x}x_1\w\cdots\w x_n$。
-
三维的 Gauss 公式证明大致思路:将体拆成一堆小块,每个小块侧面是悬崖,每个 $(x,y)$ 对应一个 $z$ 区间,不过俯视图不一定是矩形。然后分三个部分,分别会得到 $X_x$,$Y_y$ 和 $Z_z$。做的过程中会得到形似 $\iint(Z(x,y,g(x,y))-Z(x,y,f(x,y)))\d x\d y$ 的,硬把里面一层写成积分即可。
-
三维的 Stokes 公式证明思路:考虑把参数化曲面上的一个小正方形,换元变到参数自变量的 $[0,1]^2$,然后利用二维 Green 公式,再变回去。用微分形式的语言来说,就是 $$ \iint_{\p\Sigma_k}\omega=\iint_{\p[0,1]^2}\Phi^*(\omega)\xlongequal{\text{Green}}\iint_{[0,1]^2}\d\Phi^*(\omega)\xlongequal{?}\iint_{[0,1]^2}\Phi^*(\d\omega)=\iint_{\Sigma_k}\d\omega $$ 证明 $\d\Phi^*=\Phi^*\d$ 可以直接展开 $\Phi$ 和 $\d$ 的定义,问题不大。(21-25)
-
微分形式统一 Stokes 公式: $$ \begin{array}{} \Omega^0&{}\xrightarrow{\quad\d\quad}{}&\Omega^1&{}\xrightarrow{\quad\d\quad}{}&\Omega^2&{}\xrightarrow{\quad\d\quad}{}&\Omega^3\\ f&&X\d x+Y\d y+Z\d z&&X\d y\w\d z+Y\d z\w\d x+Z\d x\w\d y&&f\d x\w\d y\w\d z\\ f&{}\xrightarrow{\quad\nabla\quad}{}&\b V&{}\xrightarrow{\quad\nabla\times\quad}{}&\b V&{}\xrightarrow{\quad\nabla\cdot\quad}{}&f\\ &\displaystyle f|_{\p\Gamma}=\int_\Gamma\nabla f\cdot\d\b x&&\displaystyle\oint_{\p\Sigma}\b V\cdot\d\b l=\iint_\Sigma(\nabla\times\b V)\cdot\d\b s&&\displaystyle\oiint_{\p\Omega}\b V\cdot\d\b s=\iiint_\Omega\nabla\cdot\b V\d v \end{array} $$
级数
-
注意以下涉及到绝对值的地方都可以换成范数。
-
写出级数收敛的 Cauchy 准则。(22-18)推论:收敛的必要条件是 $|u_n|\to 0$,注意这个对积分不一定成立。
-
证明 $\sum n^{-x}$ 收敛当且仅当 $x>1$。(22-21)
-
写出积分判别法,并证明。(22-23)
-
写出比较判别法,并证明。(22-32)
-
幂级数的收敛是怎样刻画的?(22-33)
提示
${<{}}r$ 收敛,$=r$ 不知道,$>r$ 发散。考虑使用比较判别法证明。 -
$a_n\le b_n\le c_n$,若 $\sum a_n$ 和 $\sum c_n$ 都收敛,则 $\sum b_n$ 收敛。进一步,若 $\sum a_n$ 或 $\sum c_n$ 绝对收敛,则全都绝对收敛。证明?(22-36)
提示
考虑 $b_n-a_n$、$c_n-a_n$,以及比较判别法、收敛的线性性。 -
写出 D'Alembert 判别法,并证明。(22-37)
-
写出 Cauchy 判别法,并证明。(22-38)
第一种的优势是比值容易求,第二种的优势在于可以处理 $1/2,1/3,1/2^2,1/3^2,\cdots$ 这种有震荡的情况。
-
利用 Cauchy 判别法写出幂级数的收敛半径的表达式。(32-44)幂级数求导、积分后的幂级数收敛半径不变。
-
写出 Raabe 判别法,并证明。(22-48)Raabe 判别法的思想是:已知 $\sum n^{-x}$ 的收敛半径,注意到 $(n+1)^x/n^x=1+x/n+\omicron(1/n)$,所以我们可以关注 $n^{-1}$ 项的系数。
-
写出 Leibniz 判别法,并证明。(23-3)
-
分析 $\sum\frac{(-1)^{n-1}}{n^\alpha+(-1)^{n-1}}$ 在不同 $\alpha$ 下的收敛情况。(23-7)
-
写出 Dirichlet 与 Abel 判别法,并证明。(23-11~23-13)
-
分析 $\sum\frac{z^n}{n^\alpha}$ 的收敛情况。(23-14)
-
证明条件收敛的实数级数不满足交换律,换句话说可以通过重排项得到任何收敛值(包括 $\pm\infty$)或震荡。(23-17~23-21)
-
证明绝对收敛的级数满足交换律。(23-22~23-23)
-
如果一个级数加括号之后收敛,且它单项 $\to 0$,则它本身也不一定收敛。考虑 $1-1+\frac12+\frac12-\frac12-\frac12+\frac13+\frac13+\frac13-\frac13-\frac13-\frac13+\cdots$。只有每个括号内的符号相同才能保证原级数收敛。
-
注意条件收敛指的就是收敛不绝对收敛。
-
绝对收敛和一致收敛是无关的。(ylj 26-55)
-
对于函数项级数 $\sum u_n(x)\to S(x)$,证明如果在 $x\to a$ 时 $S(x)\to+\infty$,那它在 $(a,b]$ 上不一致收敛。
然而逆命题是不对的,给出一个反例,$S(x)$ 在 $[a,b]$ 上连续,然而在 $(a,b]$ 上不一致收敛。
insane counterexample: ylj 26-67
-
如何说明函数项级数在某个区间上不一致收敛?(24-3)
-
函数项级数 $\sum u_n(x)\overset I\rightrightarrows S(x)$,证明 $S(x)$ 保持有界性与连续性;若每项都可积,则 $S(x)$ 可积。(24-14~24-15)
-
证明:若导数的级数一致收敛且原级数单点收敛,则原函数一致收敛且 $S(x)$ 可导。(24-16~24-17)
-
如何说明每个 Taylor 级数都对应多个 $\mathscr{C}^\infty$ 函数?
提示
考虑函数 $f(x)=\begin{cases}\e^{-1/x^2},&x\ne 0\\\ 0,&x=0\end{cases}$ -
说明一个幂级数的封闭形式是某个函数的方法:构造微分方程,利用初值问题解的唯一性。
傅里叶级数
-
$\sum_n\frac{\sin nx}n$ 这类级数的求法:首先由于 $$ \sum_{n\le N}\sin nx=\frac{1}{\sin\frac x2}\sum_{n\le N}\sin nx\sin\frac x2=\frac{1}{2\sin\frac x2}\sum_{n\le N}\left(\cos\left(n-\frac12\right)x-\cos\left(n+\frac12\right)x\right)=\frac{\cos\frac x2-\cos(N+\frac 12)x}{2\sin\frac x2}=\frac{\sin\frac{(N+1)x}2\sin\frac{Nx}2}{\sin\frac x2} $$ 有界,故由 Dirichlet,它收敛,且在 $x$ 的取值范围不包含(接近)$2k\pi$ 时一致收敛。其次 $$ \sum_{n\le N}\frac{\sin nx}{n}=\sum_{n\le N}\int_\pi^x\cos nt\d t=\int_\pi^x\sum_{n\le N}\cos nt\d t=\int_\pi^x\frac{\sin(N+\frac12)t-\sin\frac t2}{2\sin\frac t2}\d t\to\frac{\pi-x}{2} $$ 其中 $\int\sin(N+\frac 12)t\d t$ 这个东西可以通过一次分部积分让分母出来一个 $2N+1$,所以 $\to 0$。
对于 $x=2k\pi$ 显然和为 $0$。
上面那个 $\to 0$ 有个一般的定理,就是 Riemann-Lebesgue 引理: $$ \lim_{n\to\infty}\int_a^bf(x)g(nx)\d x=\frac1T\left(\int_a^bf(x)\d x\right)\left(\int_0^Tg(x)\d x\right) $$
-
写出 $\mathscr{L}^p(I)$ 函数的定义。(26-3)
-
傅里叶级数的定义: $$ a_n=\frac1\pi\int_{-\pi}^{+\pi}f(x)\cos nx\d x,b_n=\frac1\pi\int_{-\pi}^{+\pi}f(x)\sin nx\d x $$
$$ f(x)\textcolor{red}\sim\frac{a_0}2+\sum_{n=1}^{\infty}(a_n\cos nx+b_n\sin nx) $$
-
定义内积:$\ip{f,g}=\int_{-\pi}^{+\pi}f(x)g(x)\d x$。
Bessel 不等式:根据正交性—勾股定理,已知 $\sum\frac{(\b v\cdot\b e_i)^2}{\b e_i\cdot\b e_i}\le\|\b v\|^2$。对于傅里叶级数的情况,即为(取极限): $$ \frac{a_0^2}{2}+\sum_n(a_n^2+b_n^2)\le\frac1\pi\int_{-\pi}^{+\pi}f(x)^2\d x $$ Parseval 等式即为上式取等。Parseval 等式与内积形式(两个不同函数)等价。
对于 $2\pi$ 周期函数 $f$,Parseval 等式成立等价于 $f$ 的傅里叶级数均方收敛: $$ \lim_{n\to\infty}\left\lVert S_n-f\right\rVert_2<\infty $$ 而均方收敛当且仅当 $f\in\mathscr{L}^2$。
-
Riemann 可积函数一定满足 Parseval 等式,从而均方收敛。证明思路:
-
首先证明方波 $f(x)=[x\bmod 2\pi\le L]$ 满足:$a_0=\frac{L}{\pi}$,$a_n=\frac{\sin nL}{n\pi}$,$b_n=\frac{1-\cos nL}{n\pi}$。于是 $$ \frac{a_0^2}{2}+\sum_n(a_n^2+b_n^2)=\frac{L^2}{2\pi^2}+\sum_n\frac{2-2\cos nL}{n^2\pi^2}\qquad\cdots(1) $$ 这里需要计算一个关键的式子: $$ \sum_n\frac{\cos nx}{n^2}=-\sum_n\left(\frac{1}{n^2}+\int_0^x\frac{\sin nt}n\d t\right)=\sum_n\frac1{n^2}-\int_0^x\frac{\pi-t}2\d t=\sum_n\frac1{n^2}-\frac{\pi x}{2}+\frac{x^2}4 $$ 神秘处理:两边同时关于 $x$ 积分: $$ 0=\sum_n\frac1{n^2}\int_0^{2\pi}\cos nx\d x=2\pi\sum_n\frac{1}{n^2}-\frac{\pi^3}3 $$ 故 $\zeta(2)=\pi^2/6$, $$ \sum_n\frac{\cos nx}{n^2}=\frac14x^2-\frac\pi2x+\frac{\pi^2}6 $$ 代回 $(1)$ 得 $a_0^2/2+\sum_n(a_n^2+b_n^2)=L/\pi$,这恰好等于 $\int_0^{2\pi}f(x)^2\d x$。
-
满足 Parseval 等式的,正交的两个函数之和也满足。这个通过展开 $f+g$ 和 $f-g$ 积分,然后代回 Bessel 不等式,可得到交叉乘积项 $=0$。
-
分段常值函数满足 Parseval 不等式。
-
使用分段常值函数逼近 $f$。通过达布积分的方式得到 $f$ 的一个下界分段函数 $g$。于是 $$ \|f-S_n(f)\|_2\le\|f-S_n(g)\|_2\le\|f-g\|_2+\|g-S_n(g)\|_2 $$ 第一个不等号是最小二乘性质。最后两个距离都可以任意小。
-
-
上述定理说明了 $\mathscr{L}^2$ 函数和 $\ell^2$ 函数通过傅里叶级数构成一个等距同构(这里 $\mathscr{L}^2$ 中两个函数相等当且仅当差的模为 $0$),求积分和求级数可以互相转换。
-
证明 $f$ 连续且一致收敛,则傅里叶级数收敛到 $f$。(27-22~27-24)
提示
用 $S_n$ 联系 $f$ 和极限 $S$。一方面用均方收敛,一方面用一致收敛。过程中可能要用 Cauchy 不等式在平方和绝对值积分之间转换。最后得到 $\int|f-S|=0$,再用连续性。 -
证明 $\sum_n n^p(|a_n|+|b_n|)$ 收敛则傅里叶级数是 $\mathscr{C}^p$ 的。(27-37)
提示
归纳。利用导数、极限与函数项级数的交换条件。 -
分段连续的情况下(不要求傅里叶级数收敛), $$ \int_\alpha^\beta f(x)\d x=\frac{a_0}2(\beta-\alpha)+\sum_n\left(a_n\frac{\sin n\alpha-\sin n\beta}{n}+b_n\frac{\cos n\beta-\cos n\alpha}{n}\right) $$
重要定理汇总
判定条件 | Cauchy | 比较 | Dirichlet | Abel |
---|---|---|---|---|
广义积分收敛 | 1 | 2 | 3 | 4 |
广义含参积分一致收敛 | 5 | 6 | 7 | 8 |
数项级数收敛 | 9 | 10 | 11 | 12 |
函数项级数一致收敛 | 13 | 14 | 15 | 16 |
- $\int_a^\omega f(x)\d x$ 收敛当且仅当 $\forall\eps>0$,$\exists b<\omega$,$\forall b_2\ge b_1>b$,$\left\lvert\int_{b_1}^{b_2}f(x)\d x\right\rvert<\eps$。
- 若 $\int_a^\omega g(x)\d x$ 收敛且 $\exists b<\omega$,$x\in(b,\omega)$ 时 $|f(x)|\le g(x)$,则 $\int_a^\omega f(x)\d x$ 收敛。
- 若 $\int_a^bf(x)\d x$(关于 $b$)有界且 $g(x)$ 单调 $\to 0$,则 $\int_a^\omega f(x)g(x)\d x$ 收敛。
- 若 $\int_a^\omega f(x)\d x$ 收敛且 $g(x)$ 单调有界,则 $\int_a^\omega f(x)g(x)\d x$ 收敛。
- $\int_a^\omega f(x,y) \d x $ 在 $I$ 上一致收敛当且仅当 $\forall \eps> 0$,$\exists b<\omega$,$\forall y\in I$,$\forall b_2\ge b_1 > b$,$\left| \int_{b_1}^{b_2} f(x,y) \d x \right| < \eps$。
- 若 $\int_a^\omega g(x,y)\d x$ 一致收敛且 $\exists b<\omega$,$x\in(b,\omega)$ 时 $|f(x,y)|\le g(x,y)$,则 $\int_a^\omega f(x,y)\d x$ 一致收敛。
- 若 $\int_a^bf(x,y)\d x$ 对 $b$ 和 $y$ 一致有界,$g(x,y)$ 关于 $x$ 单调且 $x\to\omega$ 时一致 $\to 0$,则 $\int_a^\omega f(x,y)g(x,y)\d x$ 一致收敛。
- 若 $\int_a^\omega f(x,y)\d x$ 一致收敛,$g(x,y)$ 关于 $x$ 单调且对 $x$ 和 $y$ 一致有界,则 $\int_a^\omega f(x,y)g(x,y)\d x$ 一致收敛。
- $\sum_na_n$ 收敛当且仅当 $\forall\eps>0$,$\exists N$,$\forall N_2\ge N_1>N$,$\left\lVert\sum_{n=N_1}^{N_2}a_n\right\rVert<\eps$。
- 若正项级数 $\sum_nb_n$ 收敛且 $\|a_n\|=\Omicron(b_n)$,则 $\sum_na_n$ 绝对收敛。
- 若 $\sum_na_n$ 部分和有界,$b_n$ 单调且 $\to 0$,则 $\sum_na_nb_n$ 收敛。
- 若 $\sum_na_n$ 收敛,$b_n$ 单调有界,则 $\sum_na_nb_n$ 收敛。
- $\sum_nf_n(x)$ 在 $I$ 上一致收敛当且仅当 $\forall\eps>0$,$\exists N$,$\forall x\in I$,$\forall N_2\ge N_1>N$,$\left\lvert\sum_{n=N_1}^{N_2}f_n(x)\right\rvert<\eps$。
- 若恒正的函数项级数 $\sum_ng_n(x)$ 一致收敛且 $\forall x\in I$,$|f_n(x)|=\Omicron(g_n(x))$,则 $\sum_nf_n(x)$ 一致收敛。
- 若 $\sum_nf_n(x)$ 部分和一致有界,$g_n(x)$ 关于 $n$ 单调且一致 $\to 0$,则 $\sum_nf_n(x)g_n(x)$ 一致收敛。
- 若 $\sum_nf_n(x)$ 一致收敛,$g_n(x)$ 关于 $n$ 单调且一致有界,则 $\sum_nf_n(x)g_n(x)$ 一致收敛。
正项级数. 若部分和有界则收敛。
积分判别法. 若 $a_n$ 能扩展为减函数 $f:[1,+\infty)\to\R$,则 $\sum_na_n=\sum_nf(n)$ 收敛 $\Longleftrightarrow$ $\int_1^{+\infty}f(x)\d x$ 收敛。
“夹逼”判别法. 若 $a_n\le b_n\le c_n$ 且 $\sum_na_n$ 与 $\sum_nc_n$ 都收敛,则 $\sum_nb_n$ 收敛。若 $\sum_na_n$ 或 $\sum_nc_n$ 绝对收敛,则另外两个也绝对收敛。
D'Alembert. 若 $\overline\lim_{n\to\infty}\left\lvert\frac{a_{n+1}}{a_n}\right\rvert<1$ 则 $\sum_na_n$ 绝对收敛,若 $\underline\lim_{n\to\infty}\left\lvert\frac{a_{n+1}}{a_n}\right\rvert>1$,则 $\sum_na_n$ 发散。
Cauchy. 若 $\overline\lim_{n\to\infty}\sqrt[n]{|a_n|}<1$ 则 $\sum_na_n$ 绝对收敛,若 $\overline\lim_{n\to\infty}\sqrt[n]{|a_n|}>1$,则 $\sum_na_n$ 发散。
幂级数. $\frac1R=\overline\lim_{n\to\infty}\sqrt[n]{|a_n|}$。
Raabe. 对于正项级数,若 $\underline\lim_{n\to\infty}n\left(\frac{a_n}{a_{n+1}}-1\right)>1$ 则 $\sum_na_n$ 收敛,若 $\overline\lim_{n\to\infty}n\left(\frac{a_n}{a_{n+1}}-1\right)<1$,则 $\sum_na_n$ 发散。
Leibniz. 对于单调不增的 $a_n\ge 0$,$\sum_n(-1)^na_n$ 收敛仅要求 $a_n\to 0$。
Hardy. 若 $\sum_nf_n(x)$ 部分和一致有界,$\sum_n|g_{n+1}(x)-g_n(x)|$ 一致收敛且 $g_n(x)$ 一致 $\to 0$,则 $\sum_nf_n(x)g_n(x)$ 一致收敛。
换序条件 | 极限 | 求导 | 积分 | 广义积分 | 级数 |
---|---|---|---|---|---|
极限 | 1 | ||||
求导 | 2 | 3 | |||
积分 | 4 | 5 | 6 | ||
广义积分 | 7 | 8 | 9 | 10 | |
级数 | 11 | 12 | 13 | 14 | 15 |
- 若 $\lim_{y\to y_0}f(x,y)$ 在 $U(x_0)$ 内一致收敛,$\lim_{x\to x_0}f(x,y)$ 在 $U(y_0)$ 内收敛,则 $\lim_{x\to x_0}\lim_{y\to y_0}f(x,y)$ 与 $\lim_{y\to y_0}\lim_{x\to x_0}f(x,y)$ 存在且相等。
- 若 $\frac{\p f}{\p y}(x,y_0)$ 在 $U(x_0)$ 内存在且求导的极限一致收敛,$\lim_{x\to x_0}f(x,y)$ 在 $U(y_0)$ 内收敛,则 $\frac{\p}{\p y}\left(\lim_{x\to x_0}f(x,y)\right)\big|_{y=y_0}$ 与 $\lim_{x\to x_0}\frac{\p f}{\p y}(x,y_0)$ 存在且相等。
- 若 $\frac{\p f}{\p x\p y}$ 在 $U(x_0,y_0)$ 内存在且连续,$\frac{\p f}{\p y\p x}(x_0,y_0)$ 存在,则 $\frac{\p f}{\p x\p y}(x_0,y_0)=\frac{\p f}{\p y\p x}(x_0,y_0)$。(加强:https://en.wikipedia.org/wiki/Symmetry_of_second_derivatives#Sufficiency_of_twice-differentiability)
- 若 $f$ 在 $[a,b]\times U(y_0)$ 内连续,则 $\lim_{y\to y_0}\int_a^bf(x,y)\d x=\int_a^bf(y_0)\d x$。
- 若 $\frac{\p f}{\p y}$ 在 $[a,b]\times U(y_0)$ 内连续,则 $\frac{\p}{\p y}\left(\int_a^bf(x,y)\d x\right)\big|_{y=y_0}=\int_a^b\frac{\p f}{\p y}(x,y_0)\d x$。
- 若 $f$ 在 $[a,b]\times[c,d]$ 内连续,则 $\int_a^b\int_c^df(x,y)\d y\d x=\int_c^d\int_a^bf(x,y)\d x\d y$。
- 若 $f$ 在 $[a,\omega)\times U(y_0)$ 内连续,且 $\int_a^\omega f(x,y)\d x$ 一致收敛,则 $\lim_{y\to y_0}\int_a^\omega f(x,y)\d x =\int_a^\omega f(x,y_0)\d x$。
- 若 $\frac{\p f}{\p y}$ 在 $[a,\omega)\times U(y_0)$ 内连续,且 $\int_a^\omega\frac{\p f}{\p y}(x,y)\d x$ 一致收敛,且 $\int_a^\omega f(x,y_0)\d x$ 收敛,则 $\int_a^\omega f(x,y)\d x$ 在闭区间上一致收敛,且 $\frac{\p}{\p y}\left(\int_a^\omega f(x,y)\d x\right)=\int_a^\omega\frac{\p f}{\p y}(x,y)\d x$。
- 若 $f$ 在 $[a,\omega)\times[c,d]$ 内连续,且 $\int_a^\omega f(x,y)\d x$ 一致收敛,则 $\int_a^\omega\int_c^df(x,y)\d y\d x=\int_c^d\int_a^\omega f(x,y)\d x\d y$。
- 若 $f$ 在 $[a,\omega)\times[c,\psi)$ 内连续,且 $\int_a^\omega f(x,y)\d x$ 一致收敛,$\int_c^\psi f(x,y)\d y$ 一致收敛,$\int_a^\omega\int_c^\psi |f(x,y)|\d y\d x$ 和 $\int_c^\psi\int_a^\omega |f(x,y)|\d x\d y$ 至少一个收敛,则 $\int_a^\omega\int_c^\psi f(x,y)\d y\d x$ 与 $\int_c^\psi\int_a^\omega f(x,y)\d x\d y$ 存在且相等。
- 若 $\sum_nf_n(x)$ 在 $U(x_0)$ 上一致收敛,且 $\lim_{x\to x_0}f_n(x)$ 存在,则 $\sum_n\lim_{x\to x_0}f_n(x)$ 与 $\lim_{x\to x_0}\sum_nf_n(x)$ 存在且相等。
- 若 $\sum_nf_n(x)$ 满足在 $U(x_0)$ 上 $f^\prime_n(x)$ 存在且连续,且 $\sum_nf^\prime_n(x)$ 一致收敛,且 $\sum_nf(x_0)$ 收敛,则 $\sum_nf_n(x)$ 在闭区间上一致收敛,且 $\left(\sum_nf_n(x)\right)^\prime=\sum_nf^\prime_n(x)$。
- 若 $\sum_nf_n(x)$ 在 $[a,b]$ 上一致收敛,且 $\int_a^bf_n(x)\d x$ 可积,则 $\sum_n\int_a^bf_n(x)\d x$ 与 $\int_a^b\sum_nf_n(x)\d x$ 存在且相等。
- 若在 $[a,\omega)$ 上 $\sum_nf_n(x)=f(x)$,且存在可积函数 $F(x)$ 使 $\forall N$,$\left\lvert\sum_{n=1}^Nf_n(x)\right\rvert\le F(x)$,则 $\sum_n\int_a^\omega f_n(x)\d x$ 与 $\int_a^\omega\sum_nf_n(x)\d x$ 存在且相等。(控制收敛定理)
- 若 $\sum_{m,n}a_{m,n}$ 能通过安排顺序变为绝对收敛的级数,则 $\sum_{m,n}a_{m,n}$ 绝对收敛,且等于 $\sum_m\sum_na_{m,n}=\sum_n\sum_ma_{m,n}$。
傅里叶级数的收敛:
- 一致收敛 $+$ 连续 $\implies$ 逐点收敛到 $f$。
- 参数级数绝对收敛 $\implies$ 一致收敛 $+$ 连续。如果 $f$ 也连续则逐点收敛到 $f$(注意 $f$ 可能有可去间断点)。
- 参数单调 $\to 0$ $\implies$ 内闭一致收敛 $+$ 开区间连续。
- 连续 $+$ 分段 $\mathscr{C}^1$ $\implies$ 一致逐点收敛到 $f$。
- 分段 $\mathscr{C}^1$ $+$ 只有第一类间断点 $\implies$ 逐点收敛到两侧极限的平均值。
- 有界 $+$ 分段单调连续 $\implies$ 逐点收敛到两侧极限的平均值。
- $\mathscr{L}^1$ $+$ 各点左右极限存在 $+$ $\int_{0^+}|(f(x-t)-f(x-)+f(x+t)-f(x+))/t|\d t<\infty$ $\implies$ 逐点收敛到两侧极限的平均值。
- (存疑)分段可微 $\implies$ 逐点收敛到两侧极限的平均值。
- $\mathscr{L}^{p>1}$ $\implies$ 几乎处处收敛。
- $\sum_n n^p(|a_n|+|b_n|)$ 收敛 $\implies$ $S(n)\in\mathscr{C}^p$。