Math Analysis Huashi Xia
只挑了几点
函数项级数
不能乱交换
第一种是我们直接让导函数发散:
f n ( x ) = sin n x n \begin{gathered}
f_n(x)=\dfrac{\sin nx}{\sqrt n}
\end{gathered} f n ( x ) = n sin n x
第二种是构造一个不可导点,用极限把很小的弯曲区域变没:
f n = ∣ x ∣ 1 + 1 n \begin{gathered}
f_n=|x|^{1+\dfrac{1}{n} }
\end{gathered} f n = ∣ x ∣ 1 + n 1
首先可以考虑构造不可积,比如我们选定极限是迪利克雷函数d ( x ) d(x) d ( x ) .于是随便取一个比如
f n ( x ) = { 1 , x = p q , q < n 0 , 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} f n ( x ) = ⎩ ⎨ ⎧ 1 , x = q p , q < n 0 , otherwise , x ∈ [ 0 , 1 ]
然后它的间断点显然有限所以可积.
那么你自然会想如果补个可积的条件呢?
一个比较容易想到的是一个向右移动的面积为1 1 1 的矩形:
f n ( 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} f n ( x ) = { 1 , x ∈ [ n , n + 1 ] 0 , otherwise
现在我们在有限区域内考虑呢?
考虑一个恒定面积,把多余的宽度累到高度上,仍然使得函数列能逐点收敛到0 0 0 :
f n ( x ) = { 2 n , x ∈ [ 2 − n , 2 − n + 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} f n ( x ) = { 2 n , x ∈ [ 2 − n , 2 − n + 1 ] 0 , otherwise , x ∈ ( 0 , 1 ]
你分析发现此时不对是因为0 0 0 处发散了,lim n → ∞ lim x → 0 f n ( x ) = + ∞ \lim_{n \to \infty} \lim_{x\to 0} f_n(x)=+\infty lim n → ∞ lim x → 0 f n ( x ) = + ∞ .
于是我们再加条件,如果极限函数可积且函数列一致有界:
函数列f n → f f_n \to f f n → f ,f n , f ∈ R [ a , b ] , ∀ x ∈ [ a , b ] , n ∈ Z , ∣ f n ( x ) ∣ ≤ M f_n,f\in R[a,b],\forall x\in [a,b],n\in Z,|f_n(x)|\le M f n , f ∈ R [ a , b ] , ∀ x ∈ [ a , b ] , n ∈ Z , ∣ f n ( x ) ∣ ≤ M
则 lim n → ∞ ∫ a b f n ( x ) d x = ∫ a b f ( x ) d x \lim_{n \to \infty} \int_a^b f_n(x) dx = \int_a^b f(x) dx lim n → ∞ ∫ a b f n ( x ) d x = ∫ a b f ( x ) d x
你发现此时在勒贝格积分是显然的 控制收敛定理(甚至可以去掉黎曼可积的条件)
在黎曼积分下呢?不会做了,发现这个叫 阿泽拉有界收敛定理(Arzelà's Bounded Convergence Theorem)并找到一个初等证明:
首先简化到 f n : [ 0 , 1 ] → [ 0 , 1 ] , f n → 0 , f n ∈ R [ 0 , 1 ] ⟹ ∫ 0 1 f n d x → 0 f_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 f n : [ 0 , 1 ] → [ 0 , 1 ] , f n → 0 , f n ∈ R [ 0 , 1 ] ⟹ ∫ 0 1 f n d x → 0 .
反证,则存在ϵ \epsilon ϵ ,使得f n f_n f n 的一个子列(仍记为f n f_n f n )满足∫ 0 1 f n d x > 4 ϵ \int_0^1 f_ndx>4\epsilon ∫ 0 1 f n d x > 4 ϵ .
根据黎曼积分定义,对每个f n f_n f n ,可以取一个分割使得∑ f n ( x ) Δ x > 2 ϵ \sum f_n(x) \Delta x>2 \epsilon ∑ f n ( x ) Δ x > 2 ϵ .且我们可以让每个区间的f n ( x ) f_n(x) f n ( x ) 取最小值,则一定存在一些区间满足f n ( x ) > ϵ f_n(x)>\epsilon f n ( x ) > ϵ ,且总长度大于ϵ \epsilon ϵ 记这些区间的并为U n U_n U n .我们可以把端点扔了弄成开区间使得U n U_n U n 是开集.U n U_n U n 中的点都f n ( x ) > ϵ f_n(x)>\epsilon f n ( x ) > ϵ .我们其实想找一个点被包含在无限个U n U_n U n 中.
此时经典令V n = ⋃ i > n U i V_n=\bigcup_{i>n} U_i V n = ⋃ i > n U i .我们的目标是说明⋂ V n ≠ ∅ \bigcap V_n\ne \varnothing ⋂ V n = ∅ .(这一步好像是处理被无限个包含的经典方法,仿照上下极限)
显然V n ⊃ V n + 1 V_n\supset V_{n+1} V n ⊃ V n + 1 ,此时一个结论说,R R R 上的任意开集V V V 可以被唯一分解成至多可数个不交开区间的并,方法是考虑对任意两个点a , b a,b a , b ,定义等价关系 a ∼ b ⟺ [ a , b ] ⊂ V a\sim b \iff [a,b]\subset V a ∼ b ⟺ [ a , b ] ⊂ V ,则每个等价类是一个开区间,就确定了分解.
分解之后,我们定义每个开区间是一个节点,如果这个节点被上层的某个开区间包含则连边,可以得到一棵树.
此时你可以让树每层只有有限个节点(如果是可数个区间扔掉一些)并保持每层总长度大于ϵ 2 \dfrac \epsilon 2 2 ϵ ,你可以扔掉没有后代的节点,因为由下层的长度保证不会扔空.
此时因为每层都有长度大于ϵ 2 \dfrac \epsilon 2 2 ϵ ,所以一定存在一条无限延伸的路.现在:
如果存在一条路最终长度不趋近到0 0 0 ,那么显然是有交的(比如考虑区间中点序列的极限点一定是交里面的)
否则路最终长度趋近0 0 0 .注意开区间不能直接用闭区间套.
考虑如果一个点在某一层至少有三个后代(不一定是儿子),那么只要我下一步走的是中间那个点,我就可以把这个开区间收缩成闭区间(因为去掉的端点被左右侧儿子挡着)
而如果这个点在后继所有点长度都维持在两个以内,又已知这两条路长度都趋近于0 0 0 ,所以它对总长度的贡献趋近于0 0 0 .
那么形式化的,我们定义坏点是从它出发的所有路的总长度趋近于0 0 0 ,这样的点及它的所有后代都可以删了,而不影响我们能找到路的情况.
则现在树上的点都是不坏的,拿出任意一个点,要么有最终也有长度的通路,要么一定有三个后代可以收缩成闭区间.
于是你可以构造一个闭区间套解决问题.
这个先用开区间性质构造树,再用树的结构构造闭区间套感觉很厉害啊!
Dini定理
若闭区间上的函数列 f n → f f_n\to f f n → f ,f n ∈ C , f ∈ C , f n < f n + 1 f_n\in C,f\in C,f_n<f_{n+1} f n ∈ C , f ∈ C , f n < f n + 1 ,则f n ⇉ f f_n\rightrightarrows f f n ⇉ f
因为f n f_n f n 单增,所以f − f n f-f_n f − f n 一定是单调减的.
所以很直觉的,对每个点x 0 x_0 x 0 ,存在N N N 使得 ∀ n > N , ∣ f n ( x 0 ) − f ( x 0 ) ∣ < δ 0 \forall n>N,|f_n(x_0)-f(x_0)|<\delta_0 ∀ n > N , ∣ f n ( x 0 ) − f ( x 0 ) ∣ < δ 0 ,因为f , f n f,f_n f , f n 连续,这个点一定有一个邻域U U U 有∀ x ∈ U , ∣ f N ( x ) − f N ( x 0 ) ∣ < δ 0 , ∣ f − f ( x 0 ) < δ 0 ∣ \forall x\in U,|f_N(x)-f_N(x_0)|<\delta_0,|f-f(x_0)<\delta_0| ∀ x ∈ U , ∣ f N ( x ) − f N ( x 0 ) ∣ < δ 0 , ∣ f − f ( x 0 ) < δ 0 ∣ 于是这个邻域内∀ n > N , ∣ f n ( x ) − f ( x ) ∣ ≤ ∣ f N ( x ) − f ( x ) ∣ ≤ ∣ f N ( x ) − f N ( x 0 ) ∣ + ∣ f N ( x 0 ) − f ( x 0 ) ∣ + ∣ f ( x ) − f ( x 0 ) ∣ = 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 ∀ n > N , ∣ f n ( x ) − f ( x ) ∣ ≤ ∣ f N ( x ) − f ( x ) ∣ ≤ ∣ f N ( x ) − f N ( x 0 ) ∣ + ∣ f N ( x 0 ) − f ( x 0 ) ∣ + ∣ f ( x ) − f ( x 0 ) ∣ = 3 δ 0 .
这就很显然了,叠一个有限覆盖就结束了.对开区间我们就只能内闭一致收敛.
一致连续推连续和积分的两个性质证起来都很显然吧.注意积分要证明可积,用不连续点是零测集可以秒.
交换极限和求导
{ f n ( x ) ∈ C 1 [ a , b ] f n ′ ( x ) ⇉ g ( x ) ∃ c ∈ [ a , b ] , lim n → ∞ f n ( c ) ∃ ⟹ { f n ( x ) ⇉ f ( x ) ∈ C 1 [ 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} ⎩ ⎨ ⎧ f n ( x ) ∈ C 1 [ a , b ] f n ′ ( x ) ⇉ g ( x ) ∃ c ∈ [ a , b ] , lim n → ∞ f n ( c ) ∃ ⟹ { f n ( x ) ⇉ f ( x ) ∈ C 1 [ a , b ] f ′ ( x ) = g ( x )
看书的时候一开始觉得这不就把积分写成导数抄一遍吗.但是课上讲的怎么这么复杂.
第一反应是因为导数连续所以可以f n ( x ) = f n ( c ) + ∫ c x f n ′ ( t ) d t f_n(x)=f_n(c)+\int_c^x f'_n(t)dt f n ( x ) = f n ( c ) + ∫ c x f n ′ ( t ) d t ,然后对任意x x x ,由刚才的定理可以得到f n ( x ) → f ( x ) f_n(x)\to f(x) f n ( x ) → f ( x ) ,所以这里得到的是f n → f f_n\to f f n → f .但没一致.
你再思考一下,一致收敛的积分性质中,是对ϵ \epsilon ϵ ,存在N N N 使得∀ n > N , ∫ a b ∣ f n ( x ) − f ( x ) ∣ d x < ϵ ( b − a ) \forall n>N,\int_a^b |f_n(x)-f(x)|dx<\epsilon(b-a) ∀ n > N , ∫ a b ∣ f n ( x ) − f ( x ) ∣ d x < ϵ ( b − a ) .所以写出来你发现这玩意真的是一致的.
因为点态收敛的结论就已经是f ( x ) = f ( c ) + ∫ c x g ( t ) f(x)=f(c)+\int_c^x g(t) f ( x ) = f ( c ) + ∫ c x g ( t ) 了求个导就做完了.
所以和把积分抄一遍的区别就在于中间那个证一致收敛啊.
Weierstrass 第一逼近定理
∀ f ∈ C [ 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 ∀ f ∈ C [ 0 , 1 ] , ϵ > 0 , ∃ polynomial P ( x ) , s . t . ∀ x , ∣ P ( x ) − f ( x ) ∣ < ϵ
定义
B n ( f , x ) = ∑ i = 0 n f ( i n ) x i ( 1 − x ) n − i ( n i ) \begin{gathered}
B_n(f,x)=\sum_{i=0}^n f(\dfrac i n)x^i(1-x)^{n-i} \binom ni
\end{gathered} B n ( f , x ) = i = 0 ∑ n f ( n i ) x i ( 1 − x ) n − i ( i n )
多元函数微分
可微的条件
我们知道可微是说你可以被平面近似.
方向导数均存在且满足条件d f d v = ∑ i f x i cos ( θ i ) \dfrac{df}{dv}=\sum_i f_{x_i} \cos(\theta_i) d v df = ∑ i f x i cos ( θ i ) 不保证函数可微
首先如果你可微的话,这个式子是一定满足的,考虑cos θ i \cos \theta_i cos θ i 实际上就是分量v i v_i v i :
d f d v = lim k → 0 f ( x + k v ) − f ( x ) k = A k v + o ( k v ) k = A v = ∑ i A i , . v = ∑ i ( d f i d x ) ⋅ v j \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} d v df = k → 0 lim k f ( x + k v ) − f ( x ) = k A k v + o ( k v ) = A v = i ∑ A i , . v = i ∑ ( d x d f i ) ⋅ v j
但是这个满足不一定可微,因为方向导数说的是你看他某一条直线上的一元函数,但直线上对不能让平面上对,实际上他甚至不一定连续啊(经典 x 2 y x 4 + y 2 \dfrac{x^2y}{x^4+y^2} x 4 + y 2 x 2 y ),如果我们让他连续,举例:
f ( x , y ) = { 0 , x = y = 0 x 2 y x 4 + y 2 ⋅ x , 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} f ( x , y ) = ⎩ ⎨ ⎧ 0 , x = y = 0 x 4 + y 2 x 2 y ⋅ x , otherwise
这么写是因为,x 2 y x 4 + y 2 \dfrac{x^2y}{x^4+y^2} x 4 + y 2 x 2 y 是一个常见的说明直线上连续不等于连续的例子,特点是y = x 2 y=x^2 y = x 2 上是 1 2 \dfrac12 2 1 ,而其他地方是0 0 0 ,乘完了之后发现:
d f d v ⃗ = f ( k v ⃗ ) − f ( 0 ) k = 1 k = 0 \begin{gathered}
\dfrac{df}{d\vec v}= \dfrac{f(k\vec v)-f(0)}{k}=\dfrac{1}{k}=0
\end{gathered} d v df = k f ( k v ) − f ( 0 ) = k 1 = 0
但f ( x , x 2 ) = 1 2 x ≠ o ( ∣ ( x , x 2 ) ∣ ) f(x,x^2)=\dfrac{1}{2}x\ne o(|(x,x^2)|) f ( x , x 2 ) = 2 1 x = o ( ∣ ( x , x 2 ) ∣ ) ,所以不对.
假设f ( x , y ) f(x,y) f ( x , y ) 的偏导f x ( x , y ) , f y ( x , y ) ∈ C ( x 0 , y 0 ) f_x(x,y),f_y(x,y)\in C(x_0,y_0) f x ( x , y ) , f y ( x , y ) ∈ C ( x 0 , y 0 ) ,则:
f ( x 0 + Δ x , y 0 + Δ y ) − f ( x , y ) = f ( x 0 + Δ x , y 0 + Δ y ) − f ( x 0 + Δ x , y 0 ) + f ( x 0 + Δ x , y 0 ) − f ( x 0 , y 0 ) = f y ( x 0 + Δ x , ξ 1 ) Δ y + o ( Δ y ) + f x ( ξ 2 , y 0 ) Δ x + o ( Δ x ) = f x ( ξ 2 , y 0 ) Δ x + f y ( x 0 + Δ x , ξ 1 ) Δ y + o ( ( Δ x ) 2 + ( Δ y ) 2 ) → f x ( x 0 , y 0 ) Δ x + f y ( x 0 , y 0 ) Δ 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} f ( x 0 + Δ x , y 0 + Δ y ) − f ( x , y ) = f ( x 0 + Δ x , y 0 + Δ y ) − f ( x 0 + Δ x , y 0 ) + f ( x 0 + Δ x , y 0 ) − f ( x 0 , y 0 ) = f y ( x 0 + Δ x , ξ 1 ) Δ y + o ( Δ y ) + f x ( ξ 2 , y 0 ) Δ x + o ( Δ x ) = f x ( ξ 2 , y 0 ) Δ x + f y ( x 0 + Δ x , ξ 1 ) Δ y + o ( ( Δ x ) 2 + ( Δ y ) 2 ) → f x ( x 0 , y 0 ) Δ x + f y ( x 0 , y 0 ) Δ y
这里还会讲一个中值定理,但一开始看他感觉很平凡,大概说的是如果在包含A , B A,B A , B 这条直线的区域上可微就有f ( A ) − f ( B ) = ∇ f ( ξ ) ⋅ A B ⃗ f(A)-f(B)=\nabla f(\xi)\cdot \vec{AB} f ( A ) − f ( B ) = ∇ f ( ξ ) ⋅ A B .你考虑在这条线上的一元拉格朗日中值,然后套上方向导数那个式子.所以实际上也很平凡.
一阶偏导在( a , b ) (a,b) ( a , b ) 邻域内连续且在( 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) 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) φ ( t ) = f ( t , b + y ) − f ( t , b ) ,则
W ( x , y ) = ϕ ′ ( a + ξ 1 ) x = x f x ( a + ξ 1 , b + y ) − x f x ( a + ξ 1 , b ) = f x y ( a , b ) x y + o ( x x 2 + y 2 ) \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 ) = ϕ ′ ( a + ξ 1 ) x = x f x ( a + ξ 1 , b + y ) − x f x ( a + ξ 1 , b ) = f x y ( a , b ) x y + o ( x x 2 + y 2 )
同理还可以得到
W ( x , y ) = f y x ( a , b ) x y + o ( y x 2 + y 2 ) \begin{gathered}
W(x,y)=f_{yx}(a,b)xy+o(y\sqrt{x^2+y^2})
\end{gathered} W ( x , y ) = f y x ( a , b ) x y + o ( y x 2 + y 2 )
那么极限
lim x → 0 W ( x , x ) x 2 = f x y ( a , b ) = f y x ( a , b ) \begin{gathered}
\lim_{x \to 0} \dfrac{W(x,x)}{x^2} =f_{xy}(a,b)=f_{yx}(a,b)
\end{gathered} x → 0 lim x 2 W ( x , x ) = f x y ( a , b ) = f y x ( a , b )
得证.
混合偏导中任意一个连续,且混合偏导都存在,则全都相等
多元的情况,可以每次交换两元.所以等价于二元情况.
不妨设是f x y f_{xy} f x y 连续.
f y x ( x 0 , y 0 ) = lim Δ x → 0 lim Δ y → 0 f ( x 0 + Δ x , y 0 + Δ y ) − f ( x 0 + Δ x , y 0 ) − f ( x 0 , y 0 + Δ y ) + f ( x 0 , y 0 ) Δ x Δ y = use lagrange for φ ( t ) = f ( x , y 0 + Δ y ) − f ( x , y 0 ) lim Δ x → 0 lim Δ y → 0 f x ( x 0 + θ 1 Δ x , y 0 + Δ y ) − f x ( x 0 + θ 1 Δ x , y 0 ) Δ y = lim Δ x → 0 lim Δ y → 0 f x y ( x 0 + θ 1 Δ x , y 0 + θ 2 Δ y ) = f x y ( x 0 , y 0 ) \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 y x ( x 0 , y 0 ) = Δ x → 0 lim Δ y → 0 lim Δ x Δ y f ( x 0 + Δ x , y 0 + Δ y ) − f ( x 0 + Δ x , y 0 ) − f ( x 0 , y 0 + Δ y ) + f ( x 0 , y 0 ) use lagrange for φ ( t ) = f ( x , y 0 + Δ y ) − f ( x , y 0 ) Δ x → 0 lim Δ y → 0 lim Δ y f x ( x 0 + θ 1 Δ x , y 0 + Δ y ) − f x ( x 0 + θ 1 Δ x , y 0 ) = Δ x → 0 lim Δ y → 0 lim f x y ( x 0 + θ 1 Δ x , y 0 + θ 2 Δ y ) = f x y ( x 0 , y 0 )
隐函数定理
{ f ( x ⃗ 0 , u ⃗ 0 ) = 0 ⃗ 在边长为 d 的闭超立方体D中,f连续且有连续偏导 det J f , u ⃗ ( ( x ⃗ 0 , u ⃗ 0 ) ) ≠ 0 ⟹ { f ( x ⃗ , u ⃗ ) = 0 ⃗ defines the function u ⃗ ( x ⃗ ) u ⃗ ( x ⃗ ) ∈ C ( D ) J u ⃗ , x ⃗ ( p ) = J f , u ⃗ ( p ) − 1 × J f , 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} ⎩ ⎨ ⎧ f ( x 0 , u 0 ) = 0 在边长为 d 的闭超立方体 D 中 ,f 连续且有连续偏导 det J f , u (( x 0 , u 0 )) = 0 ⟹ ⎩ ⎨ ⎧ f ( x , u ) = 0 defines the function u ( x ) u ( x ) ∈ C ( D ) J u , x ( p ) = J f , u ( p ) − 1 × J f , x ( p ) u has continuous partial derivative
其中J a , b ( p ) J_{a,b}(p) J a , b ( p ) 表示p p p 处函数a a a 对b b b 的 Jaccob Matrix.
要求u u u 和f f f 的维数相同.
这里有两个思路:
一是基于单调性和降维归纳的做法:
我们先考虑最简单的情况:u , x u,x u , x 是一维,f : R 2 → R f:R^2\to R f : R 2 → R .
此时由于x 0 x^0 x 0 处d f d u ≠ 0 \dfrac {df}{du}\ne 0 d u df = 0 ,不妨设d f d u > 0 \dfrac {df}{du}>0 d u df > 0 ,则由连续性在x x x 的小邻域内都有d f d u > 0 \dfrac {df}{du}>0 d u df > 0 ,只看这个小邻域,f f f 关于u u u 单调.可以找到δ \delta δ 使得f ( x 0 , u 0 + δ ) > ϵ , f ( x 0 , u 0 − δ ) < − ϵ f(x^0,u^0+\delta)>\epsilon,f(x^0,u^0-\delta)<-\epsilon f ( x 0 , u 0 + δ ) > ϵ , f ( x 0 , u 0 − δ ) < − ϵ .又由连续性,又能找到δ 1 \delta_1 δ 1 使得∀ x ∈ [ x 0 − δ 1 , x 0 + δ 1 ] , f ( x , u 0 + δ ) > ϵ 2 > 0 , f ( x , u 0 − δ ) < − ϵ 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 ∈ [ x 0 − δ 1 , x 0 + δ 1 ] , f ( x , u 0 + δ ) > 2 ϵ > 0 , f ( x , u 0 − δ ) < − 2 ϵ < 0 .从而∀ x ∈ [ x 0 − δ 1 , x 0 + δ 1 ] \forall x\in[x^0-\delta_1,x^0+\delta_1] ∀ x ∈ [ x 0 − δ 1 , x 0 + δ 1 ] ,在u u u 上使用介值定理都存在唯一一个u 1 u_1 u 1 使得f ( x , u 1 ) = 0 f(x,u_1)=0 f ( x , u 1 ) = 0 .
上面走了一个 ( x 0 , u 0 ) → ( x 0 , u 0 ± δ ) → ( [ x 0 − δ 1 , x 0 + δ 1 ] , u 0 ± δ ) (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 0 , u 0 ) → ( x 0 , u 0 ± δ ) → ([ x 0 − δ 1 , x 0 + δ 1 ] , u 0 ± δ ) 的过程,完成了对函数和导数局部性质的应用.
对于连续性只需要再跑一遍上面这个流程:对刚才那个范围内的任意一点( x , u ) (x,u) ( x , u ) ,对任意足够小的ϵ \epsilon ϵ ,考虑它也处在单调范围内,所以可以找( x , u ± ϵ ) (x,u\pm \epsilon) ( x , u ± ϵ ) 的函数值异号,然后分别用连续性拓展到x x x 上的一段,用介值定理说明存在x 1 ∈ N ( x , δ ) , ∣ f ( x 1 ) − f ( x ) ∣ < ϵ x^1\in N(x,\delta),|f(x^1)-f(x)|<\epsilon x 1 ∈ N ( x , δ ) , ∣ f ( x 1 ) − f ( x ) ∣ < ϵ .
对于导函数,考虑对这个范围内的一点( x 1 , u 1 ) (x^1,u^1) ( x 1 , u 1 ) ,和它邻域内趋近( x 1 , u 1 ) (x^1,u^1) ( x 1 , u 1 ) 的一点( x , u ) (x,u) ( x , u ) ,则由拉格朗日中值:
f ( x 1 , u 1 ) = f ( x , u ) = 0 ⟹ ∃ ξ s . t . d f d x ( x 1 − x ) + d f d u ( u 1 − u ) = 0 ⟹ u 1 − u x 1 − x = − d f d x ( ξ ) d f d u ( ξ ) ⟹ u ′ ( x ) = − d f d x d f d u ( 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} f ( x 1 , u 1 ) = f ( x , u ) = 0 ⟹ ∃ ξ s . t . d x df ( x 1 − x ) + d u df ( u 1 − u ) = 0 ⟹ x 1 − x u 1 − u = − d u df ( ξ ) d x df ( ξ ) ⟹ u ′ ( x ) = − d u df d x df ( x , u )
从而可导且连续
然后开始拓展:首先你注意到x ⃗ \vec x x 维度的拓展是容易的:证存在性和连续性的部分只是从原来的( [ x 0 − δ 1 , x 0 + δ 1 ] , u ) ([x^0-\delta_1,x^0+\delta_1],u) ([ x 0 − δ 1 , x 0 + δ 1 ] , u ) 这个区域改写成( N ( x ⃗ 0 ) , u ⃗ ) (N(\vec x^0),\vec u) ( N ( x 0 ) , u ) ,其实是一样的.而导函数的部分现在是对每维固定住其他的分别走这个流程,也是完全一样的.
要拓展u ⃗ \vec u u 的维度是比较复杂的:u ⃗ \vec u u 的维度和f f f 的维度(也就是方程的个数)是对应的,现在你需要处理多个方程.考虑雅可比矩阵行列式不为0 0 0 的意思是,这个地方线性化替代是成立的.所以直接仿照线性化的过程:应用一维的情况,把u ⃗ \vec u u 的第一分量u 1 u_1 u 1 看成因变量,其他的看成自变量,用一维的隐函数定理,可以得到u 1 u_1 u 1 关于其他所有变量的隐函数,然后用f ( u 1 ( … ) , … ) f(u_1(\ldots),\ldots) f ( u 1 ( … ) , … ) 代替原来的f f f 就能把维度减1 1 1 .
于是,因为雅可比矩阵不全为0 0 0 一定可以找到一个变量u i u_i u i 满足f 1 u i ≠ 0 \dfrac{f_1}{u_i}\ne 0 u i f 1 = 0 ,对f 1 , u i f_1,u_i f 1 , u i 使用结论,对新函数f ′ = f ( u i ( u ⃗ , x ⃗ ) , u ⃗ , x ⃗ ) f'=f(u_i(\vec u,\vec x),\vec u,\vec x) f ′ = f ( u i ( u , x ) , u , x ) 求导得到:
d f j ′ d u k = d f j d u i d u i d u k + d f j d u k \begin{gathered}
\dfrac{df'_j}{du_k}=\dfrac{df_j}{du_i}\dfrac{du_i}{du_k}+\dfrac{df_j}{du_k}
\end{gathered} d u k d f j ′ = d u i d f j d u k d u i + d u k d f j
所以矩阵的变化是删掉第一行第i i i 列,然后每列加上原来第i i i 列的某系数倍,一定仍然满秩.且验证其他条件也均是符合的,可以归纳下去.可以拿到所有的存在性,连续性,连续导数性质.
而要得到求导公式,考虑链式法则:
J f , u ⃗ J u ⃗ , x ⃗ = J f , x ⃗ \begin{gathered}
J_{f,\vec u} J_{\vec u,\vec x}=J_{f,\vec x}
\end{gathered} J f , u J u , x = J f , x
然后移过去就完了.
至高Gemini表示另一个做法是使用Banach不动点定理:即压缩映射一定有不动点,内容在250915-math-analysis里面.
所以可以构造压缩映射ϕ ( x ⃗ , u ⃗ ) = u ⃗ − J f , u ⃗ ( u ⃗ ) − 1 f ( x ⃗ , u ⃗ ) \phi(\vec x,\vec u)=\vec u-J_{f,\vec u}(\vec u)^{-1}f(\vec x,\vec u) ϕ ( x , u ) = u − J f , u ( u ) − 1 f ( x , u ) ,则有:
J ϕ , u ⃗ ( u ⃗ 1 ) = I − J f , u ⃗ ( u ⃗ 0 ) − 1 J f , u ⃗ ( u ⃗ 1 ) \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} J ϕ , u ( u 1 ) = I − J f , u ( u 0 ) − 1 J f , u ( u 1 )
则u ⃗ 1 = u ⃗ 0 \vec u^1=\vec u^0 u 1 = u 0 时是全0 0 0 ,则由偏导的连续性,存在一个小邻域使得矩阵的模(sup ∣ A u ⃗ ∣ ∣ u ⃗ ∣ \sup \dfrac{|A\vec u|}{|\vec u|} sup ∣ u ∣ ∣ A u ∣ )不超过1 2 < 1 \dfrac{1}{2}<1 2 1 < 1 ,于是在这个邻域内是压缩映射.则在这个范围内任意x ⃗ \vec x x 都会由不动点也就是方程的解.
对于连续性,注意不到
∣ u ⃗ ( x ⃗ ) − u ⃗ ( x ⃗ 1 ) ∣ = ∣ ϕ ( x ⃗ , u ⃗ ( x ⃗ ) ) − ϕ ( x ⃗ 1 , u ⃗ ( x ⃗ 1 ) ) ∣ ≤ ∣ ϕ ( x ⃗ , u ⃗ ( x ⃗ ) ) − ϕ ( x ⃗ , u ⃗ ( x ⃗ 1 ) ) ∣ + ∣ ϕ ( x ⃗ , u ⃗ ( x ⃗ 1 ) ) − ϕ ( x ⃗ 1 , u ⃗ ( x ⃗ 1 ) ) ∣ ≤ 1 2 ∣ u ⃗ ( x ⃗ ) − u ⃗ ( x ⃗ 1 ) ∣ + ∣ ϕ ( x ⃗ , u ⃗ ( x ⃗ 1 ) ) − ϕ ( x ⃗ 1 , u ⃗ ( x ⃗ 1 ) ) ∣ \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} ∣ u ( x ) − u ( x 1 ) ∣ = ∣ ϕ ( x , u ( x )) − ϕ ( x 1 , u ( x 1 )) ∣ ≤ ∣ ϕ ( x , u ( x )) − ϕ ( x , u ( x 1 )) ∣ + ∣ ϕ ( x , u ( x 1 )) − ϕ ( x 1 , u ( x 1 )) ∣ ≤ 2 1 ∣ u ( x ) − u ( x 1 ) ∣ + ∣ ϕ ( x , u ( x 1 )) − ϕ ( x 1 , u ( x 1 )) ∣
最后一步是压缩映射的性质.最后移项取极限即可得到连续性.
对于可微性,可以对每个分量用中值定理,用思路一做.
总之很麻烦啊!
重积分
面积定义
对于一个平面点集T T T ,对它任意进行划分格子,然后定义格子的集合A A A 为被完全包含的格子的集合,B B B 为与点集有交的格子的集合,C C C 为B − A B-A B − A .对每个集合内部的格子的面积累加得到S ( A ) , S ( B ) , S ( C ) S(A),S(B),S(C) S ( A ) , S ( B ) , S ( C ) ,则T T T 为可求面积的当且仅当lim S ( C ) = 0 \lim S(C)=0 lim S ( C ) = 0 ,定义其面积为S ( A ) = S ( B ) S(A)=S(B) S ( A ) = S ( B ) .
我们应该给他起个黎曼测度之类的名字的,因为这个可求面积确实只是在这种情况下的定义.
区别在于 黎曼测度只允许你用有限个格子逼近,而勒贝格测度允许你用可数个格子逼近
Fubini
闭超立方区域D = × i [ a i , b i ] D=\times_i [a_i,b_i] D = × i [ a i , b i ] 上的f f f 是重积分可积的,且某个次序的累次积分存在,则它们相等.
定义∫ n = ( ∫ … ∫ ) \int^n=(\int \ldots \int) ∫ n = ( ∫ … ∫ ) .
不妨设关于x 1 x_1 x 1 是可积的.
∫ n f ( d x 1 ∧ … ∧ d x n ) = ∫ n − 1 ∫ f d x 1 ( d x 2 ∧ … ∧ d x n ) \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} ∫ n f ( d x 1 ∧ … ∧ d x n ) = ∫ n − 1 ∫ f d x 1 ( d x 2 ∧ … ∧ d x n )
寻找一个分割组合的列( T n 1 , T n 2 ) (T^1_n,T^2_n) ( T n 1 , T n 2 ) 分别是作用再第一维的和剩下其他维空间的分割,满足 lim n → ∞ ∥ T n i ∥ = 0 \lim_{n \to \infty} \|T^i_n\|=0 lim n → ∞ ∥ T n i ∥ = 0 .
则只需证
lim n → ∞ ∑ T n 1 , T n 2 f ( x ) ∏ i Δ x i = lim n → ∞ ∑ T n 2 ∑ T n 1 f ( x ) ∏ i Δ x i = lim n → ∞ ∑ T n 2 ∫ a 1 b 1 f ( x ) d x 1 ∏ i ≠ 1 Δ x i \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} n → ∞ lim T n 1 , T n 2 ∑ f ( x ) i ∏ Δ x i = n → ∞ lim T n 2 ∑ T n 1 ∑ f ( x ) i ∏ Δ x i = n → ∞ lim T n 2 ∑ ∫ a 1 b 1 f ( x ) d x 1 i = 1 ∏ Δ x i
因为我们有
lim n → ∞ ∑ T n 2 ∫ a 1 b 1 f ( x ) d x 1 ∏ i ≠ 1 Δ x i = lim n → ∞ ∑ T n 2 ∑ i ∫ T n 1 , i T n 1 , i + 1 f ( x ) ∏ i ≠ 1 Δ x i ∈ [ lim n → ∞ ∑ T n 2 ∑ i M Δ x 1 ∏ i ≠ 1 Δ x i , lim n → ∞ ∑ T n 2 ∑ i m Δ x 1 ∏ i ≠ 1 Δ x i ] \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} n → ∞ lim T n 2 ∑ ∫ a 1 b 1 f ( x ) d x 1 i = 1 ∏ Δ x i = n → ∞ lim T n 2 ∑ i ∑ ∫ T n 1 , i T n 1 , i + 1 f ( x ) i = 1 ∏ Δ x i ∈ [ n → ∞ lim T n 2 ∑ i ∑ M Δ x 1 i = 1 ∏ Δ x i , n → ∞ lim T n 2 ∑ i ∑ m Δ x 1 i = 1 ∏ Δ x i ]
其中m m m 和M M M 分别是区域内的最大值和最小值.
然后你因为后面这两个都是收敛的直接就证明完了.
我真是糖丸了.
Peano曲线
peano curve是一类曲线而不是一个.
书上的种类看起来比wikipeida上更容易证明:构造说的是,对于一个等边三角形,连接它的中心和两个顶点构造一条折线,这样联通了它的两个顶点.
现在把一个边长为1 1 1 的定边三角形,取所有边的中点,并过它们做平行于三条边的所有的直线,可以把三角形分成四个小的,定义这个操作为一次分割.
定义曲线c n c_n c n 为进行了n n n 次分割后得到4 n 4^n 4 n 个三角形,在每个三角形内部进行一次连接,使得最后连出一个经过全部4 n 4^n 4 n 个小三角形的曲线.f n f_n f n 为[ 0 , 1 ] [0,1] [ 0 , 1 ] 到这条曲线上长度比例位置点的映射.
现在要证明这个构造是连续且都映射上的.
观察这个递归结构你会发现对任意m ≤ n m\le n m ≤ n 和t ∈ [ 0 , 1 ] t\in [0,1] t ∈ [ 0 , 1 ] ,可以找到边长为1 2 m \dfrac 1{2^m} 2 m 1 的三角形同时包含f m ( t ) , f n ( t ) f_m(t),f_n(t) f m ( t ) , f n ( t ) ,所以∣ f m ( t ) − f n ( t ) ∣ ≤ 1 2 m |f_m(t)-f_n(t)|\le \dfrac 1{2^m} ∣ f m ( t ) − f n ( t ) ∣ ≤ 2 m 1 ,一致收敛,连续.
要证明满,只需注意对三角形内任意一点,到c n c_n c n 这条线的最短距离一定小于边长1 2 n \dfrac 1 {2^n} 2 n 1 ,又有∣ f n − f ∣ < ϵ |f_n-f|<\epsilon ∣ f n − f ∣ < ϵ ,所以可以取n n n 拿到f f f 像上任意近的点.也就是说三角形内的点都是f f f 的像集的聚点.又因为连续映射把紧集映射成紧集,所以f ( [ 0 , 1 ] ) f([0,1]) f ([ 0 , 1 ]) 是闭集,包含它的所有聚点.
变量代换
U U U 是R n R^n R n 上的开集,D D D 是其中边界分段光滑的区域,则
{ f ( x ) ∈ C ( D ) g ( x ) is a bijection on D g ( x ) ∈ C 1 ( D ) ⟹ ∬ D f ( 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 ) ∈ C ( D ) g ( x ) is a bijection on D g ( x ) ∈ C 1 ( D ) ⟹ ∬ D f ( g ( x )) d Σ = ∬ g ( D ) f ( x ) d Σ
本源映射
若f ( x ) f(x) f ( x ) 与x x x 最多只有一维分量不同,则f ( x ) f(x) f ( x ) 是本源映射.
若g g g 是本源映射,则一个小超立方体区域D 1 = × i = 1 n [ a i , b i ] D_1=\times_{i=1}^n [a_i,b_i] D 1 = × i = 1 n [ a i , b i ] 有
m g ( D 1 ) = det J g , x ( x 0 ) m ( D 1 ) , x 0 ∈ D 1 \begin{gathered}
m g(D_1)=\det J_{g,x}(x_0) m(D_1),x_0\in D_1
\end{gathered} m g ( D 1 ) = det J g , x ( x 0 ) m ( D 1 ) , x 0 ∈ D 1
因为g g g 是双射,我们可以说它det J \det J det J 不为0 0 0 .我们让它在这个小区域内J J J 不变号.
考虑g g g 的J J J 是什么,假设它只改变第k k k 维,那么J J J 删掉第k k k 行/列一定是I n − 1 I_{n-1} I n − 1 ,而第k k k 行是d g n d x i \dfrac{dg_n}{dx_i} d x i d g n ,第k k k 列其他地方都是0 0 0 ,而行列式就直接是d g n d x n \dfrac{dg_n}{dx_n} d x n d g n .所以它也不为0 0 0 ,函数关于第k k k 维单调.
不妨设g g g 在这一维是单调增,设D 1 D_1 D 1 去掉第k k k 维的n − 1 n-1 n − 1 维区域是 D 1 ′ D_1' D 1 ′ .考虑:
m g ( D 1 ) = ∫ a n b n … ∫ a 1 b 1 ( 1 d x 1 … d x n ) = ∫ … ∫ ( ∫ g ( x 1 , … , x k − 1 , a k , x k + 1 , … , x n ) x 1 , … , x k − 1 , b k , x k + 1 , … , x n 1 d x k ) ( d x 1 … d x k − 1 d x k + 1 … d x n ) = ( ∏ i ≠ k ( b i − a i ) ) ( g ( x 1 , … , x k − 1 , b k , x k + 1 , … , x n ) − g ( x 1 , … , x k − 1 , a k , x k + 1 , … , x n ) ) = ∏ i ( b i − a i ) d g d x k ( x 1 , … , x k − 1 , ξ , x k + 1 , … , x n ) , ξ ∈ [ a k , b k ] = d g d x k ( x 1 , … , x k − 1 , ξ , x k + 1 , … , x n ) m g ( D 1 ) \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} m g ( D 1 ) = ∫ a n b n … ∫ a 1 b 1 ( 1 d x 1 … d x n ) = ∫ … ∫ ( ∫ g ( x 1 , … , x k − 1 , a k , x k + 1 , … , x n ) x 1 , … , x k − 1 , b k , x k + 1 , … , x n 1 d x k ) ( d x 1 … d x k − 1 d x k + 1 … d x n ) = ( i = k ∏ ( b i − a i )) ( g ( x 1 , … , x k − 1 , b k , x k + 1 , … , x n ) − g ( x 1 , … , x k − 1 , a k , x k + 1 , … , x n )) = i ∏ ( b i − a i ) d x k d g ( x 1 , … , x k − 1 , ξ , x k + 1 , … , x n ) , ξ ∈ [ a k , b k ] = d x k d g ( x 1 , … , x k − 1 , ξ , x k + 1 , … , x n ) m g ( D 1 )
我们证明的是小立方体,所以自然要考虑重积分的定义把区域切成小立方体和一些边界的不规则图形,记这些小区域的集合为{ D n } \{D_n\} { D n } .
假设现在所有小区域中,那些是小立方体的标号集合是A A A (对应了内部),剩下的小区域的标号集合是B B B .
用定义拆开:
( ∫ … ∫ ) D f ( g ( x ⃗ ) ) d Σ = lim ∥ T ∥ → 0 ∑ i ∈ A f ( ξ i ) m g ( D i ) + ∑ i ∈ B f ( ξ i ) m g ( D i ) , ξ i ∈ g ( D i ) \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} ( ∫ … ∫ ) D f ( g ( x )) d Σ = ∥ T ∥ → 0 lim i ∈ A ∑ f ( ξ i ) m g ( D i ) + i ∈ B ∑ f ( ξ i ) m g ( D i ) , ξ i ∈ g ( D i )
对第一部分可以直接刚刚证明的结论,因为我们可以保证结论中x 0 x_0 x 0 在原来的格子里.但第二部分可能会有x 0 x_0 x 0 落在D D D 区域外面,但边界很小,所以可以特殊处理:
= lim ∥ T ∥ → 0 ∑ A f ( ξ i ) m g ( D i ) + ∑ C f ( ξ i ) m g ( D i ) + ( ∑ C f ( ξ i ) m D i − ∑ C f ( ξ i ) m g ( D i ) ) = lim ∥ T ∥ → 0 ∑ A f ( ξ i ) m ( D i ) ∣ J g , x ( ξ i ′ ) ∣ + ∑ C f ( ξ i ) m ( D i ) ∣ J g , x ( ξ i ′ ) ∣ + ( ∑ C f ( ξ i ) m D i − ∑ C f ( ξ i ) m g ( D i ) ) \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} = ∥ T ∥ → 0 lim A ∑ f ( ξ i ) m g ( D i ) + C ∑ f ( ξ i ) m g ( D i ) + ( C ∑ f ( ξ i ) m D i − C ∑ f ( ξ i ) m g ( D i )) = ∥ T ∥ → 0 lim A ∑ f ( ξ i ) m ( D i ) ∣ J g , x ( ξ i ′ ) ∣ + C ∑ f ( ξ i ) m ( D i ) ∣ J g , x ( ξ i ′ ) ∣ + ( C ∑ f ( ξ i ) m D i − C ∑ f ( ξ i ) m g ( D i ))
注意ξ i \xi_i ξ i 是g ( D i ) g(D_i) g ( D i ) 中的点,而我们可以选取ξ i = g ( ξ i ′ ) \xi_i=g(\xi_i') ξ i = g ( ξ i ′ ) ,则就会变成
= ( ∫ … ∫ ) g ( D ) f ⋅ J g , x d Σ + lim ∥ T ∥ → 0 ( ∑ C f ( ξ i ) m D i − ∑ C f ( ξ i ) m g ( D i ) ) \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} = ( ∫ … ∫ ) g ( D ) f ⋅ J g , x d Σ + ∥ T ∥ → 0 lim ( C ∑ f ( ξ i ) m D i − C ∑ f ( ξ i ) m g ( D i ))
注意到因为f f f 有界,g g g 的雅可比矩阵因为连续所以有界,设有共同的界M M M ,所以后面这个括号的东西有:
∣ ∑ C f ( ξ i ) m D i − ∑ C f ( ξ i ) m g ( D i ) ∣ ≤ M ∑ C m ( D i ) + M 2 ∑ C m ( D i ) \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} ∣ C ∑ f ( ξ i ) m D i − C ∑ f ( ξ i ) m g ( D i ) ∣ ≤ M C ∑ m ( D i ) + M 2 C ∑ m ( D i )
而显然∑ C m ( D i ) \sum_C m(D_i) ∑ C m ( D i ) 是趋近于0 0 0 的,所以这部分最后极限是0 0 0 .
于是证毕.
还可以容易的用Fubini说明对分量置换(仅改变x i x_i x i 分量顺序的变换)是成立的.
g ∈ C 1 ( D ) g\in C^1(D) g ∈ C 1 ( D ) 在x x x 的邻域上可以被分解为有限个本源映射和分量置换映射的复
考虑因为雅可比矩阵可逆,可以进行置换使得d g n d x n ≠ 0 \dfrac{dg_n}{dx_n}\ne 0 d x n d g n = 0 ,于是在小邻域内构造
h ( x 1 , … , x n ) = ( x 1 , … , x n − 1 , g n ( x 1 , … , x n ) ) \begin{gathered}
h(x_1,\ldots,x_n)=(x_1,\ldots,x_{n-1},g_n(x_1,\ldots,x_n))
\end{gathered} h ( x 1 , … , x n ) = ( x 1 , … , x n − 1 , g n ( x 1 , … , x n ))
则由反函数定理,这个家伙的雅可比矩阵行列式是d g n d x n ≠ 0 \dfrac{dg_n}{dx_n}\ne 0 d x n d g n = 0 ,存在反函数h − 1 h^{-1} h − 1 .而且h − 1 h^{-1} h − 1 也是只改变x n x_n x n 的本源映射).
显然g ∘ h − 1 ∘ h = g g\circ h^{-1}\circ h=g g ∘ h − 1 ∘ h = g ,但h h h 过后x n x_n x n 的分量就已经是g ( x n ) g(x_n) g ( x n ) 的值了,所以g ∘ h − 1 g\circ h^{-1} g ∘ h − 1 是不改变x n x_n x n 的变换.
考虑这意味着你可以对g g g 再进行置换,考虑如果g i g_i g i 是不改变其中i i i 个分量的变换,那么把一个会改变的扔到x n x_n x n ,得到h i h_i h i ,然后因为x n ∘ h i − 1 ∘ h i x_n\circ h_i^{-1}\circ h_i x n ∘ h i − 1 ∘ h i ,一定是不改变i + 1 i+1 i + 1 个元素的,这样你不断重复就能拆完了.
现在的问题是对本源映射的分解也都是小区域,所以考虑进行切分:对U U U 中的每个点可以找一个开集邻域,再找到一个D D D 的有限覆盖.设有限覆盖中开集的最小半径是δ \delta δ .
那么再去划分网格,当网格中区域的直径最大值小于δ 2 \dfrac\delta 2 2 δ 时,一定每个格都完全被一个邻域完全包含,可以用这个邻域对应的分解套用刚才的证明.最后就证明完了.
可积与绝对可积
对定义在无限区域U U U 上的函数f f f ,若可以用简单且面积为0 0 0 的轨迹l l l 在U U U 内画出有限区域D l D_l D l ,定义其直径为d ( D l ) d(D_l) d ( D l ) ,则积分的值定义为
lim d ( D l ) → + ∞ ∫ D l f ( x ) d x \begin{gathered}
\lim_{d(D_l) \to +\infty} \int_{D_l} f(x)dx
\end{gathered} d ( D l ) → + ∞ lim ∫ D l f ( x ) d x
注意这要求它对任意划分l l l 都成立.
在R n ( n ≥ 2 ) R^n(n\ge 2) R n ( n ≥ 2 ) 空间以上,如果一个无限区域上的函数可积则它绝对可积.
含参积分
含参正常积分
对形如
F ( x ) = ∫ c ( x ) d ( x ) f ( x , y ) d y , x ∈ [ a , b ] \begin{gathered}
F(x)=\int_{c(x)}^{d(x)} f(x,y)dy,x\in [a,b]
\end{gathered} F ( x ) = ∫ c ( x ) d ( x ) f ( x , y ) d y , x ∈ [ a , b ]
的函数,若对任意x x x 都可积有定义,则F F F 是关于x x x 的函数,是含参量的正常积分.
若f f f 连续,c ( x ) , d ( x ) c(x),d(x) c ( x ) , d ( x ) 连续,则F F F 连续.
先考虑c ( x ) , d ( x ) c(x),d(x) c ( x ) , d ( x ) 为常数的情况.
则对任意x ∈ [ a , b ] x\in [a,b] x ∈ [ a , b ] ,对任意ϵ > 0 \epsilon>0 ϵ > 0 ,考虑
∣ F ( x ) − F ( x + δ ) ∣ = ∫ c d ( f ( x , y ) − f ( x + δ , y ) ) d y \begin{gathered}
|F(x)-F(x+\delta)|=\int_c^d (f(x,y)-f(x+\delta,y))dy
\end{gathered} ∣ F ( x ) − F ( x + δ ) ∣ = ∫ c d ( f ( x , y ) − f ( x + δ , y )) d y
因为f f f 连续,而闭区间上连续推一致连续,所以存在δ \delta δ 使得f ( x , y ) − f ( x + δ , y ) < ϵ 0 f(x,y)-f(x+\delta,y)<\epsilon_0 f ( x , y ) − f ( x + δ , y ) < ϵ 0 .于是就成了
∣ F ( x ) − F ( x + δ ) ∣ ≤ ϵ 0 ( d − c ) \begin{gathered}
|F(x)-F(x+\delta)|\le \epsilon_0(d-c)
\end{gathered} ∣ F ( x ) − F ( x + δ ) ∣ ≤ ϵ 0 ( d − c )
那么取ϵ 0 = ϵ d − c \epsilon_0=\dfrac{\epsilon}{d-c} ϵ 0 = d − c ϵ 即可.
现在考虑c ( x ) , d ( x ) c(x),d(x) c ( x ) , d ( x ) 不为常数,对F ( x ) = ∫ c ( x ) d ( x ) f ( x , y ) d y F(x)=\int_{c(x)}^{d(x)}f(x,y)dy F ( x ) = ∫ c ( x ) d ( x ) f ( x , y ) d y ,换元y = c ( x ) + ( d ( x ) − c ( x ) ) t y=c(x)+(d(x)-c(x))t y = c ( x ) + ( d ( x ) − c ( x )) t :
∫ c ( x ) d ( x ) f ( x , y ) d y = ∫ 0 1 f ( x , c ( x ) + ( d ( x ) − c ( x ) ) t ) ( d ( x ) − c ( x ) ) d t \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} ∫ c ( x ) d ( x ) f ( x , y ) d y = ∫ 0 1 f ( x , c ( x ) + ( d ( x ) − c ( x )) t ) ( d ( x ) − c ( x )) d t
此时对这个x , t x,t x , t 的含参积分用刚才的结论即知连续.
若闭区域上f ( x , y ) , d f d x , c ( x ) , d ( x ) f(x,y),\dfrac{df}{dx},c(x),d(x) f ( x , y ) , d x df , c ( x ) , d ( x ) 连续,则F ( x ) = ∫ c ( x ) d ( x ) f ( x , y ) d y F(x)=\int_{c(x)}^{d(x)}f(x,y)dy F ( x ) = ∫ c ( x ) d ( x ) f ( x , y ) d y 可微,且
d F d x = ∫ c ( x ) d ( x ) d f d x d y + f ( x , d ( x ) ) ⋅ d d d x − f ( x , c ( x ) ) ⋅ d c d x \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} d x d F = ∫ c ( x ) d ( x ) d x df d y + f ( x , d ( x )) ⋅ d x dd − f ( x , c ( x )) ⋅ d x d c
仍然先考虑c ( x ) , d ( x ) c(x),d(x) c ( x ) , d ( x ) 为常数的情况.
lim δ → 0 F ( x + δ ) − F ( x ) δ = lim δ → 0 1 δ ∫ c d ( f ( x + δ , y ) − f ( x , y ) ) d y = lim δ → 0 1 δ ∫ c d δ d f d x ( ξ ( y ) , y ) d y \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} δ → 0 lim δ F ( x + δ ) − F ( x ) = δ → 0 lim δ 1 ∫ c d ( f ( x + δ , y ) − f ( x , y )) d y = δ → 0 lim δ 1 ∫ c d δ d x df ( ξ ( y ) , y ) d y
于是
lim δ → 0 F ( x + δ ) − F ( x ) δ − ∫ c d d f d x d y = lim δ → 0 1 δ ∫ c d δ d f d x ( ξ ( y ) , y ) d y − ∫ c d d f d x ( x , y ) d y = lim δ → 0 ∫ c d ( d f d x ( ξ ( y ) , y ) − d f d x ( x , y ) ) d y \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} δ → 0 lim δ F ( x + δ ) − F ( x ) − ∫ c d d x df d y = δ → 0 lim δ 1 ∫ c d δ d x df ( ξ ( y ) , y ) d y − ∫ c d d x df ( x , y ) d y = δ → 0 lim ∫ c d ( d x df ( ξ ( y ) , y ) − d x df ( x , y )) d y
因为d f d x \dfrac{df}{dx} d x df 连续所以闭区域上一致连续,所以当δ \delta δ 足够小的时候,一定有∣ d f d x ( ξ ) − d f d x ( y ) ∣ < ϵ |\dfrac{df}{dx}(\xi)-\dfrac{df}{dx}(y)|<\epsilon ∣ d x df ( ξ ) − d x df ( y ) ∣ < ϵ ,于是可以证明这个极限为0 0 0 .
对非常数的情况,设H ( x , c , d ) = ∫ c d f ( x , y ) d y , 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)) H ( x , c , d ) = ∫ c d f ( x , y ) d y , F = H ( x , c ( x ) , d ( x )) ,用上面的可以得到H H H 对x x x 的导数,所以对H H H 全微分就做完了.
gemini补充:应该先使用H H H 偏导连续说明它全微分存在
闭区域上f f f 连续的话F F F 连续所以F F F 可积,此时会考虑能否交换积分号:
若f ( x , y ) ∈ C ( [ a , b ] × [ c , d ] ) f(x,y)\in C([a,b]\times [c,d]) f ( x , y ) ∈ C ([ a , b ] × [ c , d ]) ,则:
∫ a b ∫ c d f ( x , y ) d x d y = ∫ c d ∫ a b f ( x , y ) d y d x \begin{gathered}
\int_a^b \int_c^d f(x,y)dxdy=\int_c^d \int_a^b f(x,y)dydx
\end{gathered} ∫ a b ∫ c d f ( x , y ) d x d y = ∫ c d ∫ a b f ( x , y ) d y d x
这也是显然的,连续说明这些累次积分都存在+Fubini.
含参反常积分
两类欧拉积分