2026-03-02

Math Analysis Huashi Xia

Math Analysis Huashi Xia

只挑了几点

函数项级数

不能乱交换

构造:求极限不保证连续性

fn=xnf_n=x^n

求导和求极限不交换

第一种是我们直接让导函数发散:

fn(x)=sinnxn\begin{gathered} f_n(x)=\dfrac{\sin nx}{\sqrt n} \end{gathered}

第二种是构造一个不可导点,用极限把很小的弯曲区域变没:

fn=x1+1n\begin{gathered} f_n=|x|^{1+\dfrac{1}{n} } \end{gathered}

积分也不行

首先可以考虑构造不可积,比如我们选定极限是迪利克雷函数d(x)d(x).于是随便取一个比如

fn(x)={1,x=pq,q<n0,otherwise,x[0,1]\begin{gathered} f_n(x)=\begin{cases} 1,x=\dfrac pq,q<n \\ 0,\text{otherwise} \end{cases},x\in [0,1] \end{gathered}

然后它的间断点显然有限所以可积.

那么你自然会想如果补个可积的条件呢?

一个比较容易想到的是一个向右移动的面积为11的矩形:

fn(x)={1,x[n,n+1]0,otherwise\begin{gathered} f_n(x)=\begin{cases} 1,x\in [n,n+1] \\ 0,\text{otherwise} \end{cases} \end{gathered}

现在我们在有限区域内考虑呢?

考虑一个恒定面积,把多余的宽度累到高度上,仍然使得函数列能逐点收敛到00:

fn(x)={2n,x[2n,2n+1]0,otherwise,x(0,1]\begin{gathered} f_n(x)=\begin{cases} 2^n,x\in [2^{-n},2^{-n+1}] \\ 0,\text{otherwise} \end{cases},x\in (0,1] \end{gathered}

你分析发现此时不对是因为00处发散了,limnlimx0fn(x)=+\lim_{n \to \infty} \lim_{x\to 0} f_n(x)=+\infty.

于是我们再加条件,如果极限函数可积且函数列一致有界:

函数列fnff_n \to f,fn,fR[a,b],x[a,b],nZ,fn(x)Mf_n,f\in R[a,b],\forall x\in [a,b],n\in Z,|f_n(x)|\le M

limnabfn(x)dx=abf(x)dx\lim_{n \to \infty} \int_a^b f_n(x) dx = \int_a^b f(x) dx

你发现此时在勒贝格积分是显然的 控制收敛定理(甚至可以去掉黎曼可积的条件)

在黎曼积分下呢?不会做了,发现这个叫 阿泽拉有界收敛定理(Arzelà's Bounded Convergence Theorem)并找到一个初等证明:

首先简化到 fn:[0,1][0,1],fn0,fnR[0,1]    01fndx0f_n:[0,1]\to [0,1],f_n\to 0,f_n\in R[0,1] \implies \int_0^1 f_n dx\to 0.

反证,则存在ϵ\epsilon,使得fnf_n的一个子列(仍记为fnf_n)满足01fndx>4ϵ\int_0^1 f_ndx>4\epsilon.

根据黎曼积分定义,对每个fnf_n,可以取一个分割使得fn(x)Δx>2ϵ\sum f_n(x) \Delta x>2 \epsilon.且我们可以让每个区间的fn(x)f_n(x)取最小值,则一定存在一些区间满足fn(x)>ϵf_n(x)>\epsilon,且总长度大于ϵ\epsilon记这些区间的并为UnU_n.我们可以把端点扔了弄成开区间使得UnU_n是开集.UnU_n中的点都fn(x)>ϵf_n(x)>\epsilon.我们其实想找一个点被包含在无限个UnU_n中.

此时经典令Vn=i>nUiV_n=\bigcup_{i>n} U_i.我们的目标是说明Vn\bigcap V_n\ne \varnothing.(这一步好像是处理被无限个包含的经典方法,仿照上下极限)

显然VnVn+1V_n\supset V_{n+1},此时一个结论说,RR上的任意开集VV可以被唯一分解成至多可数个不交开区间的并,方法是考虑对任意两个点a,ba,b,定义等价关系 ab    [a,b]Va\sim b \iff [a,b]\subset V,则每个等价类是一个开区间,就确定了分解.

分解之后,我们定义每个开区间是一个节点,如果这个节点被上层的某个开区间包含则连边,可以得到一棵树.

此时你可以让树每层只有有限个节点(如果是可数个区间扔掉一些)并保持每层总长度大于ϵ2\dfrac \epsilon 2,你可以扔掉没有后代的节点,因为由下层的长度保证不会扔空.

此时因为每层都有长度大于ϵ2\dfrac \epsilon 2,所以一定存在一条无限延伸的路.现在:

如果存在一条路最终长度不趋近到00,那么显然是有交的(比如考虑区间中点序列的极限点一定是交里面的)

否则路最终长度趋近00.注意开区间不能直接用闭区间套.

考虑如果一个点在某一层至少有三个后代(不一定是儿子),那么只要我下一步走的是中间那个点,我就可以把这个开区间收缩成闭区间(因为去掉的端点被左右侧儿子挡着)

而如果这个点在后继所有点长度都维持在两个以内,又已知这两条路长度都趋近于00,所以它对总长度的贡献趋近于00.

那么形式化的,我们定义坏点是从它出发的所有路的总长度趋近于00,这样的点及它的所有后代都可以删了,而不影响我们能找到路的情况.

则现在树上的点都是不坏的,拿出任意一个点,要么有最终也有长度的通路,要么一定有三个后代可以收缩成闭区间.

于是你可以构造一个闭区间套解决问题.

这个先用开区间性质构造树,再用树的结构构造闭区间套感觉很厉害啊!

Dini定理

若闭区间上的函数列 fnff_n\to f,fnC,fC,fn<fn+1f_n\in C,f\in C,f_n<f_{n+1},则fnff_n\rightrightarrows f

因为fnf_n单增,所以ffnf-f_n一定是单调减的.

所以很直觉的,对每个点x0x_0,存在NN使得 n>N,fn(x0)f(x0)<δ0\forall n>N,|f_n(x_0)-f(x_0)|<\delta_0,因为f,fnf,f_n连续,这个点一定有一个邻域UUxU,fN(x)fN(x0)<δ0,ff(x0)<δ0\forall x\in U,|f_N(x)-f_N(x_0)|<\delta_0,|f-f(x_0)<\delta_0|于是这个邻域内n>N,fn(x)f(x)fN(x)f(x)fN(x)fN(x0)+fN(x0)f(x0)+f(x)f(x0)=3δ0\forall n>N,|f_n(x)-f(x)| \le |f_N(x)-f(x)|\le |f_N(x)-f_N(x_0)|+|f_N(x_0)-f(x_0)|+|f(x)-f(x_0)|=3\delta_0.

这就很显然了,叠一个有限覆盖就结束了.对开区间我们就只能内闭一致收敛.

一致连续推连续和积分的两个性质证起来都很显然吧.注意积分要证明可积,用不连续点是零测集可以秒.

交换极限和求导

{fn(x)C1[a,b]fn(x)g(x)c[a,b],limnfn(c)    {fn(x)f(x)C1[a,b]f(x)=g(x)\begin{gathered} \begin{cases} f_n(x)\in C^1[a,b] \\ f_n'(x)\rightrightarrows g(x) \\ \exists c\in [a,b],\lim_{n \to \infty} f_n(c)\exists \end{cases} \\ \implies \begin{cases} f_n(x)\rightrightarrows f(x)\in C^1[a,b] \\ f'(x)=g(x) \end{cases} \end{gathered}

看书的时候一开始觉得这不就把积分写成导数抄一遍吗.但是课上讲的怎么这么复杂.

第一反应是因为导数连续所以可以fn(x)=fn(c)+cxfn(t)dtf_n(x)=f_n(c)+\int_c^x f'_n(t)dt,然后对任意xx,由刚才的定理可以得到fn(x)f(x)f_n(x)\to f(x),所以这里得到的是fnff_n\to f.但没一致.

你再思考一下,一致收敛的积分性质中,是对ϵ\epsilon,存在NN使得n>N,abfn(x)f(x)dx<ϵ(ba)\forall n>N,\int_a^b |f_n(x)-f(x)|dx<\epsilon(b-a).所以写出来你发现这玩意真的是一致的.

因为点态收敛的结论就已经是f(x)=f(c)+cxg(t)f(x)=f(c)+\int_c^x g(t)了求个导就做完了.

所以和把积分抄一遍的区别就在于中间那个证一致收敛啊.

Weierstrass 第一逼近定理

fC[0,1],ϵ>0, polynomial P(x), s.t. x,P(x)f(x)<ϵ\forall f\in C[0,1],\epsilon>0,\exists \text{ polynomial } P(x), \ s.t.\ \forall x,|P(x)-f(x)|< \epsilon

定义

Bn(f,x)=i=0nf(in)xi(1x)ni(ni)\begin{gathered} B_n(f,x)=\sum_{i=0}^n f(\dfrac i n)x^i(1-x)^{n-i} \binom ni \end{gathered}

多元函数微分

可微的条件

我们知道可微是说你可以被平面近似.

方向导数均存在且满足条件dfdv=ifxicos(θi)\dfrac{df}{dv}=\sum_i f_{x_i} \cos(\theta_i)不保证函数可微

首先如果你可微的话,这个式子是一定满足的,考虑cosθi\cos \theta_i实际上就是分量viv_i:

dfdv=limk0f(x+kv)f(x)k=Akv+o(kv)k=Av=iAi,.v=i(dfidx)vj\begin{gathered} \dfrac{df}{dv} =\lim_{k \to 0} \dfrac{f(x+kv)-f(x)}{k} =\dfrac{Akv+o(kv)}{k}=Av=\sum_i A_{i,.}v=\sum_i (\dfrac{df_i}{dx}) \cdot v_j \end{gathered}

但是这个满足不一定可微,因为方向导数说的是你看他某一条直线上的一元函数,但直线上对不能让平面上对,实际上他甚至不一定连续啊(经典 x2yx4+y2\dfrac{x^2y}{x^4+y^2}),如果我们让他连续,举例:

f(x,y)={0,x=y=0x2yx4+y2x,otherwise\begin{gathered} f(x,y)=\begin{cases} 0, x=y=0 \\ \dfrac{x^2y}{x^4+y^2} \cdot x,\text{otherwise} \end{cases} \end{gathered}

这么写是因为,x2yx4+y2\dfrac{x^2y}{x^4+y^2}是一个常见的说明直线上连续不等于连续的例子,特点是y=x2y=x^2上是 12\dfrac12,而其他地方是00,乘完了之后发现:

dfdv=f(kv)f(0)k=1k=0\begin{gathered} \dfrac{df}{d\vec v}= \dfrac{f(k\vec v)-f(0)}{k}=\dfrac{1}{k}=0 \end{gathered}

f(x,x2)=12xo((x,x2))f(x,x^2)=\dfrac{1}{2}x\ne o(|(x,x^2)|),所以不对.

偏导数连续则函数可微

假设f(x,y)f(x,y)的偏导fx(x,y),fy(x,y)C(x0,y0)f_x(x,y),f_y(x,y)\in C(x_0,y_0),则:

f(x0+Δx,y0+Δy)f(x,y)=f(x0+Δx,y0+Δy)f(x0+Δx,y0)+f(x0+Δx,y0)f(x0,y0)=fy(x0+Δx,ξ1)Δy+o(Δy)+fx(ξ2,y0)Δx+o(Δx)=fx(ξ2,y0)Δx+fy(x0+Δx,ξ1)Δy+o((Δx)2+(Δy)2)fx(x0,y0)Δx+fy(x0,y0)Δy\begin{gathered} f(x_0+\Delta x,y_0+\Delta y)-f(x,y) \\ =f(x_0+\Delta x,y_0+\Delta y)-f(x_0+\Delta x,y_0)+f(x_0+\Delta x,y_0)-f(x_0,y_0) \\ =f_y(x_0+\Delta x,\xi_1)\Delta y+o(\Delta y)+f_x(\xi_2,y_0)\Delta x+o(\Delta x) \\ =f_x(\xi_2,y_0)\Delta x+f_y(x_0+\Delta x,\xi_1)\Delta y+o({\sqrt{(\Delta x)^2+(\Delta y)^2}}) \\ \to f_x(x_0,y_0)\Delta x+f_y(x_0,y_0)\Delta y \end{gathered}

这里还会讲一个中值定理,但一开始看他感觉很平凡,大概说的是如果在包含A,BA,B这条直线的区域上可微就有f(A)f(B)=f(ξ)ABf(A)-f(B)=\nabla f(\xi)\cdot \vec{AB}.你考虑在这条线上的一元拉格朗日中值,然后套上方向导数那个式子.所以实际上也很平凡.

一阶偏导在(a,b)(a,b)邻域内连续且在(a,b)(a,b)处可微则(a,b)(a,b)处混合偏导相等.

W(x,y)=f(a+x,b+y)f(a,b+y)f(a+x,b)+f(a,b)W(x,y)=f(a+x,b+y)-f(a,b+y)-f(a+x,b)+f(a,b).

然后为了一致,对φ(t)=f(t,b+y)f(t,b)\varphi(t)=f(t,b+y)-f(t,b),则

W(x,y)=ϕ(a+ξ1)x=xfx(a+ξ1,b+y)xfx(a+ξ1,b)=fxy(a,b)xy+o(xx2+y2)\begin{gathered} W(x,y)=\phi'(a+\xi_1)x=xf_x(a+\xi_1,b+y)-xf_x(a+\xi_1,b) \\ =f_{xy}(a,b)xy+o(x\sqrt{x^2+y^2}) \end{gathered}

同理还可以得到

W(x,y)=fyx(a,b)xy+o(yx2+y2)\begin{gathered} W(x,y)=f_{yx}(a,b)xy+o(y\sqrt{x^2+y^2}) \end{gathered}

那么极限

limx0W(x,x)x2=fxy(a,b)=fyx(a,b)\begin{gathered} \lim_{x \to 0} \dfrac{W(x,x)}{x^2} =f_{xy}(a,b)=f_{yx}(a,b) \end{gathered}

得证.

混合偏导中任意一个连续,且混合偏导都存在,则全都相等

多元的情况,可以每次交换两元.所以等价于二元情况.

不妨设是fxyf_{xy}连续.

fyx(x0,y0)=limΔx0limΔy0f(x0+Δx,y0+Δy)f(x0+Δx,y0)f(x0,y0+Δy)+f(x0,y0)ΔxΔy=use lagrange for φ(t)=f(x,y0+Δy)f(x,y0)limΔx0limΔy0fx(x0+θ1Δx,y0+Δy)fx(x0+θ1Δx,y0)Δy=limΔx0limΔy0fxy(x0+θ1Δx,y0+θ2Δy)=fxy(x0,y0)\begin{gathered} f_{yx}(x_0,y_0)=\lim_{\Delta x \to 0}\lim_{\Delta y \to 0} \dfrac{f(x_0+\Delta x,y_0+\Delta y)-f(x_0+\Delta x,y_0)-f(x_0,y_0+\Delta y)+f(x_0,y_0)}{\Delta x\Delta y} \\ \xlongequal{ \text{use lagrange for } \varphi(t)=f(x,y_0+\Delta y)-f(x,y_0) } \\ \lim_{\Delta x \to 0} \lim_{\Delta y \to 0} \dfrac{f_x(x_0+\theta_1\Delta x,y_0+\Delta y)-f_x(x_0+\theta_1\Delta x,y_0)}{\Delta y} \\ =\lim_{\Delta x \to 0} \lim_{\Delta y \to 0} f_{xy}(x_0+\theta_1\Delta x,y_0+\theta_2\Delta y) \\ =f_{xy}(x_0,y_0) \end{gathered}

隐函数定理

{f(x0,u0)=0在边长为d的闭超立方体D中,f连续且有连续偏导detJf,u((x0,u0))0    {f(x,u)=0 defines the function u(x)u(x)C(D)Ju,x(p)=Jf,u(p)1×Jf,x(p)u has continuous partial derivative\begin{gathered} \begin{cases} f(\vec x^0,\vec u^0)=\vec 0 \\ \text{在边长为$d$的闭超立方体D中,f连续且有连续偏导} \\ \det J_{f,\vec u}((\vec x^0,\vec u^0) ) \ne 0 \end{cases} \\ \implies \begin{cases} f(\vec x,\vec u)=\vec 0 \text{ defines the function } \vec u(\vec x) \\ \vec u(\vec x) \in C(D) \\ J_{\vec u,\vec x}(p)=J_{f,\vec u}(p)^{-1}\times J_{f,\vec x}(p) \\ \vec u \text{ has continuous partial derivative} \end{cases} \end{gathered}

其中Ja,b(p)J_{a,b}(p)表示pp处函数aabb的 Jaccob Matrix. 要求uuff的维数相同.

这里有两个思路:

一是基于单调性和降维归纳的做法:

我们先考虑最简单的情况:u,xu,x是一维,f:R2Rf:R^2\to R.

此时由于x0x^0dfdu0\dfrac {df}{du}\ne 0,不妨设dfdu>0\dfrac {df}{du}>0,则由连续性在xx的小邻域内都有dfdu>0\dfrac {df}{du}>0,只看这个小邻域,ff关于uu单调.可以找到δ\delta使得f(x0,u0+δ)>ϵ,f(x0,u0δ)<ϵf(x^0,u^0+\delta)>\epsilon,f(x^0,u^0-\delta)<-\epsilon.又由连续性,又能找到δ1\delta_1使得x[x0δ1,x0+δ1],f(x,u0+δ)>ϵ2>0,f(x,u0δ)<ϵ2<0\forall x\in [x^0-\delta_1,x^0+\delta_1],f(x,u^0+\delta)>\dfrac \epsilon 2>0,f(x,u^0-\delta)<-\dfrac \epsilon2<0.从而x[x0δ1,x0+δ1]\forall x\in[x^0-\delta_1,x^0+\delta_1],在uu上使用介值定理都存在唯一一个u1u_1使得f(x,u1)=0f(x,u_1)=0.

上面走了一个 (x0,u0)(x0,u0±δ)([x0δ1,x0+δ1],u0±δ)(x^0,u^0)\to (x^0,u^0\pm \delta) \to ([x^0-\delta_1,x^0+\delta_1],u^0\pm \delta)的过程,完成了对函数和导数局部性质的应用.

对于连续性只需要再跑一遍上面这个流程:对刚才那个范围内的任意一点(x,u)(x,u),对任意足够小的ϵ\epsilon,考虑它也处在单调范围内,所以可以找(x,u±ϵ)(x,u\pm \epsilon)的函数值异号,然后分别用连续性拓展到xx上的一段,用介值定理说明存在x1N(x,δ),f(x1)f(x)<ϵx^1\in N(x,\delta),|f(x^1)-f(x)|<\epsilon.

对于导函数,考虑对这个范围内的一点(x1,u1)(x^1,u^1),和它邻域内趋近(x1,u1)(x^1,u^1)的一点(x,u)(x,u),则由拉格朗日中值:

f(x1,u1)=f(x,u)=0    ξ s.t. dfdx(x1x)+dfdu(u1u)=0    u1ux1x=dfdx(ξ)dfdu(ξ)    u(x)=dfdxdfdu(x,u)\begin{gathered} f(x^1,u^1)=f(x,u)=0 \\ \implies \exists \xi \ s.t.\ \dfrac{df}{dx}(x^1-x)+\dfrac{df}{du}(u^1-u)=0 \\ \implies \dfrac{u^1-u}{x^1-x}=-\dfrac{\dfrac{df}{dx}(\xi) }{\dfrac{df}{du}(\xi) } \\ \implies u'(x)=-\dfrac{\dfrac{df}{dx} }{\dfrac{df}{du} } (x,u) \end{gathered}

从而可导且连续

然后开始拓展:首先你注意到x\vec x维度的拓展是容易的:证存在性和连续性的部分只是从原来的([x0δ1,x0+δ1],u)([x^0-\delta_1,x^0+\delta_1],u)这个区域改写成(N(x0),u)(N(\vec x^0),\vec u),其实是一样的.而导函数的部分现在是对每维固定住其他的分别走这个流程,也是完全一样的.

要拓展u\vec u的维度是比较复杂的:u\vec u的维度和ff的维度(也就是方程的个数)是对应的,现在你需要处理多个方程.考虑雅可比矩阵行列式不为00的意思是,这个地方线性化替代是成立的.所以直接仿照线性化的过程:应用一维的情况,把u\vec u的第一分量u1u_1看成因变量,其他的看成自变量,用一维的隐函数定理,可以得到u1u_1关于其他所有变量的隐函数,然后用f(u1(),)f(u_1(\ldots),\ldots)代替原来的ff就能把维度减11.

于是,因为雅可比矩阵不全为00一定可以找到一个变量uiu_i满足f1ui0\dfrac{f_1}{u_i}\ne 0,对f1,uif_1,u_i使用结论,对新函数f=f(ui(u,x),u,x)f'=f(u_i(\vec u,\vec x),\vec u,\vec x)求导得到:

dfjduk=dfjduiduiduk+dfjduk\begin{gathered} \dfrac{df'_j}{du_k}=\dfrac{df_j}{du_i}\dfrac{du_i}{du_k}+\dfrac{df_j}{du_k} \end{gathered}

所以矩阵的变化是删掉第一行第ii列,然后每列加上原来第ii列的某系数倍,一定仍然满秩.且验证其他条件也均是符合的,可以归纳下去.可以拿到所有的存在性,连续性,连续导数性质.

而要得到求导公式,考虑链式法则:

Jf,uJu,x=Jf,x\begin{gathered} J_{f,\vec u} J_{\vec u,\vec x}=J_{f,\vec x} \end{gathered}

然后移过去就完了.

至高Gemini表示另一个做法是使用Banach不动点定理:即压缩映射一定有不动点,内容在250915-math-analysis里面.

所以可以构造压缩映射ϕ(x,u)=uJf,u(u)1f(x,u)\phi(\vec x,\vec u)=\vec u-J_{f,\vec u}(\vec u)^{-1}f(\vec x,\vec u),则有:

Jϕ,u(u1)=IJf,u(u0)1Jf,u(u1)\begin{gathered} J_{\phi,\vec u}(\vec u^1)=I-J_{f,\vec u}(\vec u^0)^{-1}J_{f,\vec u}(\vec u^1) \end{gathered}

u1=u0\vec u^1=\vec u^0时是全00,则由偏导的连续性,存在一个小邻域使得矩阵的模(supAuu\sup \dfrac{|A\vec u|}{|\vec u|})不超过12<1\dfrac{1}{2}<1,于是在这个邻域内是压缩映射.则在这个范围内任意x\vec x都会由不动点也就是方程的解.

对于连续性,注意不到

u(x)u(x1)=ϕ(x,u(x))ϕ(x1,u(x1))ϕ(x,u(x))ϕ(x,u(x1))+ϕ(x,u(x1))ϕ(x1,u(x1))12u(x)u(x1)+ϕ(x,u(x1))ϕ(x1,u(x1))\begin{gathered} |\vec u(\vec x)-\vec u(\vec x^1)|=|\phi(\vec x,\vec u(\vec x))-\phi(\vec x^1,\vec u(\vec x^1))| \\ \le |\phi(\vec x,\vec u(\vec x))-\phi(\vec x,\vec u(\vec x^1))|+|\phi(\vec x,\vec u(\vec x^1))-\phi(\vec x^1,\vec u(\vec x^1))| \\ \le \dfrac{1}{2}|\vec u(\vec x)-\vec u(\vec x^1)|+|\phi(\vec x,\vec u(\vec x^1))-\phi(\vec x^1,\vec u(\vec x^1))| \end{gathered}

最后一步是压缩映射的性质.最后移项取极限即可得到连续性.

对于可微性,可以对每个分量用中值定理,用思路一做.

总之很麻烦啊!

重积分

面积定义

对于一个平面点集TT,对它任意进行划分格子,然后定义格子的集合AA为被完全包含的格子的集合,BB为与点集有交的格子的集合,CCBAB-A.对每个集合内部的格子的面积累加得到S(A),S(B),S(C)S(A),S(B),S(C),则TT为可求面积的当且仅当limS(C)=0\lim S(C)=0,定义其面积为S(A)=S(B)S(A)=S(B).

我们应该给他起个黎曼测度之类的名字的,因为这个可求面积确实只是在这种情况下的定义.

区别在于 黎曼测度只允许你用有限个格子逼近,而勒贝格测度允许你用可数个格子逼近

Fubini

闭超立方区域D=×i[ai,bi]D=\times_i [a_i,b_i]上的ff是重积分可积的,且某个次序的累次积分存在,则它们相等.

定义n=()\int^n=(\int \ldots \int).

不妨设关于x1x_1是可积的.

nf(dx1dxn)=n1fdx1(dx2dxn)\begin{gathered} \int^n f(dx_1\wedge\ldots\wedge dx_n)=\int^{n-1} \int fdx_1 (dx_2\wedge\ldots\wedge dx_n) \end{gathered}

寻找一个分割组合的列(Tn1,Tn2)(T^1_n,T^2_n)分别是作用再第一维的和剩下其他维空间的分割,满足 limnTni=0\lim_{n \to \infty} \|T^i_n\|=0.

则只需证

limnTn1,Tn2f(x)iΔxi=limnTn2Tn1f(x)iΔxi=limnTn2a1b1f(x)dx1i1Δxi\begin{gathered} \lim_{n \to \infty} \sum_{T^1_n,T^2_n} f(x)\prod_i \Delta x_i \\ =\lim_{n \to \infty} \sum_{T^2_n}\sum_{T^1_n}f(x)\prod_i\Delta x_i \\ =\lim_{n \to \infty} \sum_{T^2_n}\int_{a_1}^{b_1}f(x)dx_1\prod_{i\ne 1}\Delta x_i \end{gathered}

因为我们有

limnTn2a1b1f(x)dx1i1Δxi=limnTn2iTn1,iTn1,i+1f(x)i1Δxi[limnTn2iMΔx1i1Δxi,limnTn2imΔx1i1Δxi]\begin{gathered} \lim_{n \to \infty} \sum_{T^2_n}\int_{a_1}^{b_1}f(x)dx_1\prod_{i\ne 1}\Delta x_i \\ =\lim_{n \to \infty} \sum_{T^2_n} \sum_{i} \int_{T_n^{1,i}}^{T_n^{1,i+1}} f(x)\prod_{i\ne 1}\Delta x_i \\ \in [ \lim_{n \to \infty} \sum_{T^2_n} \sum_{i} M\Delta x_1\prod_{i\ne 1}\Delta x_i ,\\ \lim_{n \to \infty} \sum_{T^2_n} \sum_{i} m\Delta x_1\prod_{i\ne 1}\Delta x_i ] \end{gathered}

其中mmMM分别是区域内的最大值和最小值.

然后你因为后面这两个都是收敛的直接就证明完了.

我真是糖丸了.

Peano曲线

构造一个线到一个面的区域的连续满映射.

peano curve是一类曲线而不是一个.

书上的种类看起来比wikipeida上更容易证明:构造说的是,对于一个等边三角形,连接它的中心和两个顶点构造一条折线,这样联通了它的两个顶点.

现在把一个边长为11的定边三角形,取所有边的中点,并过它们做平行于三条边的所有的直线,可以把三角形分成四个小的,定义这个操作为一次分割.

定义曲线cnc_n为进行了nn次分割后得到4n4^n个三角形,在每个三角形内部进行一次连接,使得最后连出一个经过全部4n4^n个小三角形的曲线.fnf_n[0,1][0,1]到这条曲线上长度比例位置点的映射.

现在要证明这个构造是连续且都映射上的.

观察这个递归结构你会发现对任意mnm\le nt[0,1]t\in [0,1],可以找到边长为12m\dfrac 1{2^m}的三角形同时包含fm(t),fn(t)f_m(t),f_n(t),所以fm(t)fn(t)12m|f_m(t)-f_n(t)|\le \dfrac 1{2^m},一致收敛,连续.

要证明满,只需注意对三角形内任意一点,到cnc_n这条线的最短距离一定小于边长12n\dfrac 1 {2^n},又有fnf<ϵ|f_n-f|<\epsilon,所以可以取nn拿到ff像上任意近的点.也就是说三角形内的点都是ff的像集的聚点.又因为连续映射把紧集映射成紧集,所以f([0,1])f([0,1])是闭集,包含它的所有聚点.

变量代换

UURnR^n上的开集,DD是其中边界分段光滑的区域,则

{f(x)C(D)g(x) is a bijection on Dg(x)C1(D)    Df(g(x))dΣ=g(D)f(x)dΣ\begin{gathered} \begin{cases} f(x)\in C(D) \\ g(x) \text{ is a bijection on } D \\ g(x)\in C^1(D) \end{cases} \\ \implies \iint_{D} f(g(x))d\Sigma=\iint_{g(D)}f(x)d\Sigma \end{gathered}

本源映射

f(x)f(x)xx最多只有一维分量不同,则f(x)f(x)是本源映射.

gg是本源映射,则一个小超立方体区域D1=×i=1n[ai,bi]D_1=\times_{i=1}^n [a_i,b_i]

mg(D1)=detJg,x(x0)m(D1),x0D1\begin{gathered} m g(D_1)=\det J_{g,x}(x_0) m(D_1),x_0\in D_1 \end{gathered}

因为gg是双射,我们可以说它detJ\det J不为00.我们让它在这个小区域内JJ不变号.

考虑ggJJ是什么,假设它只改变第kk维,那么JJ删掉第kk行/列一定是In1I_{n-1},而第kk行是dgndxi\dfrac{dg_n}{dx_i},第kk列其他地方都是00,而行列式就直接是dgndxn\dfrac{dg_n}{dx_n}.所以它也不为00,函数关于第kk维单调.

不妨设gg在这一维是单调增,设D1D_1去掉第kk维的n1n-1维区域是 D1D_1'.考虑:

mg(D1)=anbna1b1(1dx1dxn)=(g(x1,,xk1,ak,xk+1,,xn)x1,,xk1,bk,xk+1,,xn1dxk)(dx1dxk1dxk+1dxn)=(ik(biai))(g(x1,,xk1,bk,xk+1,,xn)g(x1,,xk1,ak,xk+1,,xn))=i(biai)dgdxk(x1,,xk1,ξ,xk+1,,xn),ξ[ak,bk]=dgdxk(x1,,xk1,ξ,xk+1,,xn)mg(D1)\begin{gathered} mg(D_1)=\int_{a_n}^{b_n}\ldots\int_{a_1}^{b_1} (1dx_1\ldots dx_n) \\ =\int\ldots \int(\int_{g(x_1,\ldots,x_{k-1},a_k,x_{k+1},\ldots,x_n)}^{x_1,\ldots,x_{k-1},b_k,x_{k+1},\ldots,x_n}1dx_k)(dx_1\ldots dx_{k-1}dx_{k+1}\ldots dx_n) \\ =(\prod_{i\ne k}(b_i-a_i))(g(x_1,\ldots ,x_{k-1},b_k,x_{k+1},\ldots,x_n)-g(x_1,\ldots,x_{k-1},a_k,x_{k+1},\ldots,x_n)) \\ =\prod_i (b_i-a_i) \dfrac {dg}{dx_k} (x_1,\ldots ,x_{k-1},\xi,x_{k+1},\ldots,x_n),\xi \in [a_k,b_k] \\ =\dfrac {dg}{dx_k} (x_1,\ldots ,x_{k-1},\xi,x_{k+1},\ldots,x_n) mg(D_1) \\ \end{gathered}

gg是本源映射,则原命题成立

我们证明的是小立方体,所以自然要考虑重积分的定义把区域切成小立方体和一些边界的不规则图形,记这些小区域的集合为{Dn}\{D_n\}.

假设现在所有小区域中,那些是小立方体的标号集合是AA(对应了内部),剩下的小区域的标号集合是BB.

用定义拆开:

()Df(g(x))dΣ=limT0iAf(ξi)mg(Di)+iBf(ξi)mg(Di),ξig(Di)\begin{gathered} \left(\int\ldots\int\right)_{D}f(g(\vec x))d\Sigma \\ =\lim_{\|T\| \to 0} \sum_{i\in A} f(\xi_i)mg(D_i) + \sum_{i\in B} f(\xi_i)mg(D_i),\xi_i\in g(D_i) \\ \end{gathered}

对第一部分可以直接刚刚证明的结论,因为我们可以保证结论中x0x_0在原来的格子里.但第二部分可能会有x0x_0落在DD区域外面,但边界很小,所以可以特殊处理:

=limT0Af(ξi)mg(Di)+Cf(ξi)mg(Di)+(Cf(ξi)mDiCf(ξi)mg(Di))=limT0Af(ξi)m(Di)Jg,x(ξi)+Cf(ξi)m(Di)Jg,x(ξi)+(Cf(ξi)mDiCf(ξi)mg(Di))\begin{gathered} =\lim_{\|T\| \to 0} \sum_A f(\xi_i) mg(D_i) + \sum_C f(\xi_i) mg(D_i) \\ + (\sum_{C} f(\xi_i) mD_i-\sum_{C}f(\xi_i)mg(D_i)) \\ =\lim_{\|T\| \to 0} \sum_A f(\xi_i) m(D_i)|J_{g,x}(\xi_i')| + \sum_C f(\xi_i) m(D_i)|J_{g,x}(\xi_i')| \\ + (\sum_{C} f(\xi_i) mD_i-\sum_{C}f(\xi_i)mg(D_i)) \\ \end{gathered}

注意ξi\xi_ig(Di)g(D_i)中的点,而我们可以选取ξi=g(ξi)\xi_i=g(\xi_i'),则就会变成

=()g(D)fJg,xdΣ+limT0(Cf(ξi)mDiCf(ξi)mg(Di))\begin{gathered} =\left(\int\ldots\int\right)_{g(D)}f\cdot J_{g,x} d\Sigma+ \lim_{\|T\| \to 0} (\sum_{C} f(\xi_i) mD_i-\sum_{C}f(\xi_i)mg(D_i)) \end{gathered}

注意到因为ff有界,gg的雅可比矩阵因为连续所以有界,设有共同的界MM,所以后面这个括号的东西有:

Cf(ξi)mDiCf(ξi)mg(Di)MCm(Di)+M2Cm(Di)\begin{gathered} |\sum_{C} f(\xi_i) mD_i-\sum_{C}f(\xi_i)mg(D_i)| \\ \le M\sum_C m(D_i)+M^2\sum_C m(D_i) \end{gathered}

而显然Cm(Di)\sum_C m(D_i)是趋近于00的,所以这部分最后极限是00.

于是证毕.

还可以容易的用Fubini说明对分量置换(仅改变xix_i分量顺序的变换)是成立的.

gC1(D)g\in C^1(D)xx的邻域上可以被分解为有限个本源映射和分量置换映射的复

考虑因为雅可比矩阵可逆,可以进行置换使得dgndxn0\dfrac{dg_n}{dx_n}\ne 0,于是在小邻域内构造

h(x1,,xn)=(x1,,xn1,gn(x1,,xn))\begin{gathered} h(x_1,\ldots,x_n)=(x_1,\ldots,x_{n-1},g_n(x_1,\ldots,x_n)) \end{gathered}

则由反函数定理,这个家伙的雅可比矩阵行列式是dgndxn0\dfrac{dg_n}{dx_n}\ne 0,存在反函数h1h^{-1}.而且h1h^{-1}也是只改变xnx_n的本源映射).

显然gh1h=gg\circ h^{-1}\circ h=g,但hh过后xnx_n的分量就已经是g(xn)g(x_n)的值了,所以gh1g\circ h^{-1}是不改变xnx_n的变换.

考虑这意味着你可以对gg再进行置换,考虑如果gig_i是不改变其中ii个分量的变换,那么把一个会改变的扔到xnx_n,得到hih_i,然后因为xnhi1hix_n\circ h_i^{-1}\circ h_i,一定是不改变i+1i+1个元素的,这样你不断重复就能拆完了.

现在的问题是对本源映射的分解也都是小区域,所以考虑进行切分:对UU中的每个点可以找一个开集邻域,再找到一个DD的有限覆盖.设有限覆盖中开集的最小半径是δ\delta.

那么再去划分网格,当网格中区域的直径最大值小于δ2\dfrac\delta 2时,一定每个格都完全被一个邻域完全包含,可以用这个邻域对应的分解套用刚才的证明.最后就证明完了.

可积与绝对可积

对定义在无限区域UU上的函数ff,若可以用简单且面积为00的轨迹llUU内画出有限区域DlD_l,定义其直径为d(Dl)d(D_l),则积分的值定义为

limd(Dl)+Dlf(x)dx\begin{gathered} \lim_{d(D_l) \to +\infty} \int_{D_l} f(x)dx \end{gathered}

注意这要求它对任意划分ll都成立.

Rn(n2)R^n(n\ge 2)空间以上,如果一个无限区域上的函数可积则它绝对可积.

含参积分

含参正常积分

对形如

F(x)=c(x)d(x)f(x,y)dy,x[a,b]\begin{gathered} F(x)=\int_{c(x)}^{d(x)} f(x,y)dy,x\in [a,b] \end{gathered}

的函数,若对任意xx都可积有定义,则FF是关于xx的函数,是含参量的正常积分.

ff连续,c(x),d(x)c(x),d(x)连续,则FF连续.

先考虑c(x),d(x)c(x),d(x)为常数的情况.

则对任意x[a,b]x\in [a,b],对任意ϵ>0\epsilon>0,考虑

F(x)F(x+δ)=cd(f(x,y)f(x+δ,y))dy\begin{gathered} |F(x)-F(x+\delta)|=\int_c^d (f(x,y)-f(x+\delta,y))dy \end{gathered}

因为ff连续,而闭区间上连续推一致连续,所以存在δ\delta使得f(x,y)f(x+δ,y)<ϵ0f(x,y)-f(x+\delta,y)<\epsilon_0.于是就成了

F(x)F(x+δ)ϵ0(dc)\begin{gathered} |F(x)-F(x+\delta)|\le \epsilon_0(d-c) \end{gathered}

那么取ϵ0=ϵdc\epsilon_0=\dfrac{\epsilon}{d-c}即可.

现在考虑c(x),d(x)c(x),d(x)不为常数,对F(x)=c(x)d(x)f(x,y)dyF(x)=\int_{c(x)}^{d(x)}f(x,y)dy,换元y=c(x)+(d(x)c(x))ty=c(x)+(d(x)-c(x))t:

c(x)d(x)f(x,y)dy=01f(x,c(x)+(d(x)c(x))t)(d(x)c(x))dt\begin{gathered} \int_{c(x)}^{d(x)}f(x,y)dy \\ =\int_0^1 f(x,c(x)+(d(x)-c(x))t)(d(x)-c(x))dt \\ \end{gathered}

此时对这个x,tx,t的含参积分用刚才的结论即知连续.

若闭区域上f(x,y),dfdx,c(x),d(x)f(x,y),\dfrac{df}{dx},c(x),d(x)连续,则F(x)=c(x)d(x)f(x,y)dyF(x)=\int_{c(x)}^{d(x)}f(x,y)dy可微,且

dFdx=c(x)d(x)dfdxdy+f(x,d(x))dddxf(x,c(x))dcdx\begin{gathered} \dfrac{dF}{dx}=\int_{c(x)}^{d(x)}\dfrac{df}{dx}dy+f(x,d(x))\cdot \dfrac{dd}{dx}-f(x,c(x))\cdot \dfrac{dc}{dx} \end{gathered}

仍然先考虑c(x),d(x)c(x),d(x)为常数的情况.

limδ0F(x+δ)F(x)δ=limδ01δcd(f(x+δ,y)f(x,y))dy=limδ01δcdδdfdx(ξ(y),y)dy\begin{gathered} \lim_{\delta \to 0} \dfrac{F(x+\delta)-F(x)}{\delta} \\ =\lim_{\delta \to 0} \dfrac{1}{\delta}\int_c^d(f(x+\delta,y)-f(x,y))dy \\ =\lim_{\delta \to 0} \dfrac1\delta \int_c^d \delta \dfrac{df}{dx}(\xi(y),y)dy \\ \end{gathered}

于是

limδ0F(x+δ)F(x)δcddfdxdy=limδ01δcdδdfdx(ξ(y),y)dycddfdx(x,y)dy=limδ0cd(dfdx(ξ(y),y)dfdx(x,y))dy\begin{gathered} \lim_{\delta \to 0} \dfrac{F(x+\delta)-F(x)}{\delta} - \int_c^d \dfrac{df}{dx}dy \\ =\lim_{\delta \to 0} \dfrac1\delta \int_c^d \delta \dfrac{df}{dx}(\xi(y),y)dy-\int_c^d \dfrac{df}{dx}(x,y)dy \\ \\ =\lim_{\delta \to 0} \int_c^d (\dfrac{df}{dx}(\xi(y),y)-\dfrac{df}{dx}(x,y))dy \end{gathered}

因为dfdx\dfrac{df}{dx}连续所以闭区域上一致连续,所以当δ\delta足够小的时候,一定有dfdx(ξ)dfdx(y)<ϵ|\dfrac{df}{dx}(\xi)-\dfrac{df}{dx}(y)|<\epsilon,于是可以证明这个极限为00.

对非常数的情况,设H(x,c,d)=cdf(x,y)dy,F=H(x,c(x),d(x))H(x,c,d)=\int_c^d f(x,y)dy,F=H(x,c(x),d(x)),用上面的可以得到HHxx的导数,所以对HH全微分就做完了.

gemini补充:应该先使用HH偏导连续说明它全微分存在

闭区域上ff连续的话FF连续所以FF可积,此时会考虑能否交换积分号:

f(x,y)C([a,b]×[c,d])f(x,y)\in C([a,b]\times [c,d]),则:

abcdf(x,y)dxdy=cdabf(x,y)dydx\begin{gathered} \int_a^b \int_c^d f(x,y)dxdy=\int_c^d \int_a^b f(x,y)dydx \end{gathered}

这也是显然的,连续说明这些累次积分都存在+Fubini.

含参反常积分

两类欧拉积分