高等微积分笔记 Lesson 29

高等微积分 Lesson 29

换元公式

/Theorem/

φ:JI\varphi:J\to IC1C^1的,设fC(I)f\in C(I),则a,bI\forall a,b\in IA,BJA,B\in J满足(端点对齐)φ(A)=a\varphi(A)=aφ(B)=b\varphi(B)=b,有

abf(x)dx=ABf(φ(t))φ(t)dt\int_a^bf(x)\text{d}x=\int_A^Bf(\varphi(t))\varphi'(t)\text{d}t

此版本的换元公式对于ff要求比较多,必须连续;而对φ\varphi的要求很少,只要求C1C^1(甚至允许φ\varphi不单). 之后在高维版本的换元公式中,对φ\varphi的要求会变高,要求φ\varphi是双射,且C1C^1;而对于ff的要求降低为要求可积.

我们这个一维版本的证明中用到了 Newton - Leibniz 公式,绕开了ff中间细致的内容,只需要端点对齐就可以换元.

下面来看一些应用.

/Claim/ (奇函数和偶函数的积分化简)

  1. 奇函数在对称区间上的积分消失.

    f(x)=f(x)f(-x)=-f(x)x\forall xI=I+II=I_+\cup I_-I=I+I_-=-I_+. 则有

    I+f+If=0\int_{I_+}f+\int_{I_-}f=0

  2. 偶函数在对称区间上的积分是单侧的两倍.

    f(x)=f(x)f(-x)=f(x)x\forall x,则

    If=I+f\int_{I_-}f=\int_{I_+}f

从几何上看这个命题是显然的,当然也可以使用换元公式给出一个证明:

/Proof/

I+=[a,b]I_+=[a,b]I=[b,a]I_-=[-b,-a]. 对奇函数有

If=baf(x)dx=x=yabf(y)dy=abf(y)dy=I+f\int_{I_-}f=\int_{-b}^{-a}f(x)\text{d}x\overset{x=-y}{=}\int_a^bf(-y)\text{d}y=-\int_a^bf(y)\text{d}y=-\int_{I_+}f

得证. 偶函数反之.

可以发现,在一维下的对称性是很容易看出来的,但是等到我们学习高维的积分时,会发现对称性难以看出来,这时候就需要利用换元公式来计算. 观察积分的对称性时,用换元公式更易实现,也更加机械化.

/Claim/ (周期函数积分)

周期函数在一个周期内的积分不依赖于起点的位置.

假设ff的周期为TT,也即f(x+T)=f(x)f(x+T)=f(x)x\forall x. 则

aa+Tf(x)dx=0Tf(x)dx\int_a^{a+T}f(x)\text{d}x=\int_0^Tf(x)\text{d}x

/Proof/

几乎是显然的. 有

LHS=aTf+Ta+Tf=x=T+yaTf+0af(T+y)dy=aTf+0af(y)dy=0Tf(x)dx=RHS\begin{aligned} \text{LHS}&=\int_a^Tf+\int_T^{a+T}f\overset{x=T+y}{=}\int_a^Tf+\int_0^af(T+y)\text{d}y\\ &=\int_a^Tf+\int_0^af(y)\text{d}y=\int_0^Tf(x)\text{d}x=\text{RHS} \end{aligned}

证毕.

一个更加奇怪的例子:

/Example/

ff是偶函数,证明:

aa(f(x)1+ex)dx=0af(x)dx\int_{-a}^a(\frac{f(x)}{1+e^x})\text{d}x=\int_0^af(x)\text{d}x

/Proof/

LHS=a0f(x)1+exdx+0af(x)1+exdx=x=ya0f(y)1+ey(dy)+0af(x)1+exdx=0aeyf(y)1+eydy+0af(x)1+exdx=0af(x)dx=RHS\begin{aligned} \text{LHS}&=\int_{-a}^0\frac{f(x)}{1+e^x}\text{d}x+\int_0^a\frac{f(x)}{1+e^x}\text{d}x\\\\ &\overset{x=-y}{=}\int_a^0\frac{f(-y)}{1+e^{-y}}(-\text{d}y)+\int_0^a\frac{f(x)}{1+e^x}\text{d}x\\\\ &=\int_0^a\frac{e^yf(y)}{1+e^y}\text{d}y+\int_0^a\frac{f(x)}{1+e^x}\text{d}x\\\\ &=\int_0^af(x)\text{d}x=\text{RHS} \end{aligned}

证毕.

要想看出对称性,基本上都使用换元公式.

/Example/

rr(r2x2)n2dx=x=rcosθπ0(r2r2cos2θ)n2(rsinθ)dθ=0πrnsinnθrsinθdθ=rn+10πsinn+1θdθ=rn+10π2sinn+1θdθ+rn+1π2πsinn+1θdθ=θ=παn+10π2sinn+1θdθ+rn+1π2πsinn+1(πα)(dα)=2rn+1In+1\begin{aligned} \int_{-r}^r(r^2-x^2)^{\frac{n}{2}}\text{d}x&\overset{x=r\cos\theta}{=}\int_\pi^0(r^2-r^2\cos^2\theta)^{\frac{n}{2}}(-r\sin\theta)\text{d}\theta\\ &=\int_0^\pi r^n\sin^n\theta\cdot r\sin\theta\text{d}\theta\\ &=r^{n+1}\int_0^\pi\sin^{n+1}\theta\text{d}\theta\\ &=r^{n+1}\int_0^\frac{\pi}{2}\sin^{n+1}\theta\text{d}\theta+r^{n+1}\int_\frac{\pi}{2}^\pi\sin^{n+1}\theta\text{d}\theta\\ &\overset{\theta=\pi-\alpha}{=}^{n+1}\int_0^\frac{\pi}{2}\sin^{n+1}\theta\text{d}\theta+r^{n+1}\int_\frac{\pi}{2}^\pi\sin^{n+1}(\pi-\alpha)(-\text{d}\alpha)\\ &=2r^{n+1}I_{n+1} \end{aligned}

回忆之前的递推公式,Im=m1mIm2I_m=\frac{m-1}{m}I_{m-2}m2\forall m\geq2. 得到结果.

/Example/ (Gauss)

电荷分布均匀的球面产生的电场E\vec{E}. 考虑球面S2={(x,y,z)x2+y2+z2=1}S^2=\{(x,y,z)|x^2+y^2+z^2=1\},所以z\vec{z}轴上(0,0,a)(0,0,a)点处的电场是

ES2dS(u2+v2+(wa)2)rr=S2(u,v,aw)(u2+v2+(wa)2)3/2dS\vec{E}\propto\iint_{S^2}\frac{\text{d}S}{(u^2+v^2+(w-a)^2)}\frac{\vec{r}}{r}=\iint_{S^2}\frac{(-u,-v,a-w)}{(u^2+v^2+(w-a)^2)^{3/2}}\text{d}S

因为xxyy方向是奇函数,积分区间是对称区间,所以这两个方向的分量消失;考虑对zz方向的分量的一维积分,应该有

(E)z0π(cosxa)sinx(1+a22acosx)3/2dx=11(ua)du(1+a22au)3/2\begin{aligned} (\vec{E})_z&\propto\int_0^\pi\frac{(\cos x-a)\sin x}{(1+a^2-2a\cos x)^{3/2}}\text{d}x=\int_{-1}^1\frac{(u-a)\text{d}u}{(1+a^2-2au)^{3/2}} \end{aligned}

如果我假装我笨一点的话……哦根本不用假装.
——艾神

考虑把分子上的uu用完,写成

=112au+1+a22a+1a22a(1+a22au)3/2du=1112a(1+a22au)1/2du+111a22a(1+a22au)3/2du=12a(1+a22au)1/212(2a)11+1a22a(1+a22au)1/212(2a)11=\begin{aligned} &=\int_{-1}^1\frac{\frac{-2au+1+a^2}{-2a}+\frac{1-a^2}{2a}}{(1+a^2-2au)^{3/2}}\text{d}u\\\\ &=\int_{-1}^1\frac{1}{-2a}(1+a^2-2au)^{-1/2}\text{d}u+\int_{-1}^1\frac{1-a^2}{2a}(1+a^2-2au)^{-3/2}\text{d}u\\\\ &=-\frac{1}{2a}\frac{(1+a^2-2au)^{1/2}}{\frac{1}{2}(-2a)}|^1_{-1}+\frac{1-a^2}{2a}\frac{(1+a^2-2au)^{-1/2}}{-\frac{1}{2}(-2a)}|^1_{-1}\\\\ &=\cdots \end{aligned}

这里一定要小心计算,不能算错,因为在考试的时候没有时间检查.

你能做完就不错了!
——艾神

我们来看看有没有算错……这看起来几乎是对的,我们就满足于此. 我漏抄了一个1a21-a^2?那这没办法,我抄漏了.
——艾神

反常积分

有时候我们想要在一个无穷的区间中做一个积分,那么这时就需要先积分再取极限. 或者有时在某点处被积函数无界,需要做截断,之后取极限.

这里我们就讲到了几种积分:

Riemann integral=boundedbounded functioninfinite integral=unboundedunbounded functionimproper integral=boundedunbounded functiom\begin{aligned} &\text{Riemann integral}=\int_{\text{bounded}}\text{bounded function}\\\\ &\text{infinite integral}=\int_{\text{unbounded}}\text{unbounded function}\\\\ &\text{improper integral}=\int_{\text{bounded}}\text{unbounded functiom} \end{aligned}

其中后面两种被称为反常积分,或者广义积分.

无穷积分

先截断、再取极限:

/Definition/

  1. a+f(x)dx=limb+abf(x)dx\int_a^{+\infty}f(x)\text{d}x=\lim_{b\to+\infty}\int_a^bf(x)\text{d}x

    称无穷积分收敛 / 发散     \iff 上述极限存在 / 不存在.

  2. bf(x)dx=limaabf(x)dx\int_{-\infty}^bf(x)\text{d}x=\lim_{a\to-\infty}\int_a^bf(x)\text{d}x

    同理.

  3. R\R上的无穷积分:

    Rf=f(x)dx=defined as0f+0f\int_\R f=\int_{-\infty}^\infty f(x)\text{d}x\overset{\text{defined as}}{=}\int_{-\infty}^0f+\int_0^\infty f

    称积分收敛     \iff 两个积分均收敛;积分发散     \iff 至少有一个积分发散.

如果每一次判断都需要用定义,未免太麻烦,所以我们写下一些命题方便判断.

/Claim/

R\R上无穷积分收敛,则可用对称截断计算:即

f=limR+RRf\int_{-\infty}^\infty f=\lim_{R\to+\infty}\int_{-R}^Rf

证明:求和即可

注:从limR+RRf\underset{R\to+\infty}{\lim}\int_{-R}^Rf无法得到Rf\int_\R f收敛,比如对任意奇函数都有极限值为00,但是f(x)=xf(x)=x的无穷积分在单侧不收敛,所以整体也不收敛.

/Example/

0exdx=limb+0bexdx=limbex0b=1\int_0^\infty e^{-x}\text{d}x=\lim_{b\to+\infty}\int_0^be^{-x}\text{d}x=\lim_{b\to\infty}-e^{-x}|^b_0=1

/Example/

1xpdx=limb1bxpdx={limbxp+1p+11b=limbbp+11p+1={divergent,p+1>0convergent,p+1<1,p1limblnx1b=limblnb=divergent,p=1\begin{aligned} \int_1^\infty x^p\text{d}x&=\lim_{b\to\infty}\int_1^bx^p\text{d}x\\\\ &=\left\{\begin{array}{lr} \underset{b\to\infty}{\lim}\frac{x^{p+1}}{p+1}|^b_1=\underset{b\to\infty}{\lim}\frac{b^{p+1}-1}{p+1}=\left\{\begin{array}{lr} \text{divergent}\,,\quad p+1>0\\\\ \text{convergent}\,,\quad p+1<1 \end{array}\right.\,,p\ne-1\\\\ \underset{b\to\infty}{\lim}\ln x|^b_1=\underset{b\to\infty}{\lim}\ln b=\text{divergent}\,,\quad p=-1 \end{array}\right. \end{aligned}

瑕积分

/Definition/

  1. 右端点bb是瑕点(ffbb点无界)

    abf=limε0+abεf=limcbacf\int_a^bf=\lim_{\varepsilon\to0+}\int_a^{b-\varepsilon}f=\lim_{c\to b-}\int_a^cf

  2. 左端点aa是瑕点(ffaa点无界)

    abf=limε0+a+εbf\int_a^bf=\lim_{\varepsilon\to0+}\int_{a+\varepsilon}^bf

  3. 内部一点cc为瑕点.

    abf=acf+cbf\int_a^bf=\int_a^cf+\int_c^bf

一个标准的例子是下面的:

/Example/

ab1(bx)qdx={divergent,q1convergent,q<1\int_a^b\frac{1}{(b-x)^q}\text{d}x=\left\{\begin{array}{lr} \text{divergent}\,,\quad q\geq1\\\\ \text{convergent}\,,\quad q<1 \end{array}\right.

判断无穷积分的收敛性

区分非负函数的无穷积分对应的方法和一般函数对应的方法.

非负函数无穷积分收敛性

也就是f0f\ge0,在[a,+)[a,+\infty)上. 则

a+f is convergent    limbF(b)\int_a^{+\infty}f\text{ is convergent}\iff\lim_{b\to\infty}F(b)

同时根据一致收敛原理,等价于{F(b)ba}\{F(b)|b\geq a\}有上界.

/Claim/

ff非负,则

af is convergent    {F(b)=abfba} has sup\int_a^\infty f\text{ is convergent}\iff\{F(b)=\int_a^bf|b\ge a\}\text{ has sup}

根据上面的命题,我们可以提出一个极为重要的定理,

/Theorem/ (比较定理)

0fg0\leq f\leq g,则

  1. ag\int_a^\infty g收敛,则af\int_a^\infty f收敛;
  2. (上面的逆否命题)若af\int_a^\infty f发散,则ag\int_a^\infty g发散.

/Proof/

只证明第 1 条即可.

ag\int_a^\infty g收敛 \Longrightarrow {G(b)=abg}\{G(b)=\int_a^bg\}有上界KK,所以F(b)G(b)KF(b)\leq G(b)\leq Kb\forall b,因此F(b)F(b)有上界,所以af\int_a^\infty f收敛.

构建不等式难于计算极限,所以我们考虑使用算极限的方式给出比较定理的另外一个形式:

f,g0f,g\geq0,且limxf(x)g(x)=LR0{\underset{x\to\infty}{\lim}\frac{f(x)}{g(x)}=L\in\R_{\geq0}\cup\{符号+}+\infty\}.

  1. ag\int_a^\infty g收敛,且LR0L\in\R_{\geq0},则af\int_a^\infty f收敛;
  2. ag\int_a^\infty g发散且L0L\neq0,则af\int_a^\infty f发散;
  3. LR+L\in\R_+,则两个积分收敛性一致.

/Proof/

只证明第 1 条:设limfg=L<+\lim\frac{f}{g}=L<+\infty,则M\exist M使得xM\forall x\geq Mf(x)g(x)<L+1\frac{f(x)}{g(x)}<L+1.

此时f(x)g(x)(L+1)f(x)\leq g(x)(L+1),由比较定理知道,因为M(L+1)g\int_M^\infty(L+1) g收敛,所以Mf\int_M^\infty f收敛.

/Example/

分析下面积分的收敛性:

I=2+dxxklnxI=\int_2^{+\infty}\frac{\text{d}x}{x^k\ln x}

  1. k>1k>1,则

    1xklnx1xkln2\frac{1}{x^k\ln x}\leq\frac{1}{x^k\ln2}

    回忆1/xk1/x^k的收敛性,由比较定理知道积分收敛.

  2. k<1k<1,则

    limx+1xklnx1xkxε=limx+xεlnx=+\lim_{x\to+\infty}\frac{\frac{1}{x^k\ln x}}{\frac{1}{x^k x^\varepsilon}}=\lim_{x\to+\infty}\frac{x^\varepsilon}{\ln x}=+\infty

    0<ε1k0<\varepsilon\leq1-kII发散.

  3. k=1k=1,就是ln(lnx)\ln(\ln x),发散.

一般函数的无穷积分

af\int_a^\infty f收敛     \iff limbF(b)\underset{b\to\infty}{\lim}F(b)存在,可以类比数列的 Cauchy 准则,得到前面的表述等价于ε>0\forall\varepsilon>0K>0\exist K>0使得b2>b1K\forall b_2>b_1\geq KF(b2)F(b1)<ε|F(b_2)-F(b_1)|<\varepsilon.

/Theorem/ (Cauchy 准则)

af\int_a^\infty f收敛     \iff ε>0\forall\varepsilon>0K>0\exist K>0使得b2>b1K\forall b_2>b_1\geq KF(b2)F(b1)<ε|F(b_2)-F(b_1)|<\varepsilon.

推论:若af\int_a^\infty|f|收敛,则af\int_a^\infty f收敛.

/Proof/ (推论的证明)

已知af\int_a^\infty|f|收敛,由 Cauchy 准则知ε>0\forall\varepsilon>0K>0\exist K>0使得b2>b1K\forall b_2>b_1\geq Kb1b2f<ε\int_{b_1}^{b_2}|f|<\varepsilon,所以

b1b2fb1b2f<ε|\int_{b_1}^{b_2}f|\leq\int_{b_1}^{b_2}|f|<\varepsilon

再用 Cauchy 准则知道af\int_a^\infty f收敛.

由上面的推论可以给出一个新的概念:

/Definition/ (绝对收敛)

af\int_a^\infty f绝对收敛,若af\int_a^\infty|f|收敛;条件收敛,若af\int_a^\infty f收敛且af\int_a^\infty|f|发散.

绝对收敛显然强于收敛.

无穷积分收敛有可能是绝对收敛,也有可能是条件收敛,一个条件收敛的例子:

/Example/

考虑

1sinxxdx\int_1^\infty\frac{\sin x}{x}\text{d}x

首先证明收敛:对于复杂的广义积分,只能分部计算.

1bsinxxdx=1b(cosx)xdx=cosxx1b1b(cosx)(1)x2dx\begin{aligned} \int_1^b\frac{\sin x}{x}\text{d}x&=\int_1^b\frac{(-\cos x)'}{x}\text{d}x=-\frac{\cos x}{x}|^b_1-\int_1^b(-\cos x)(-1)x^{-2}\text{d}x \end{aligned}

所以

1sinxxdx=cos11cosxx2dx\int_1^\infty\frac{\sin x}{x}\text{d}x=\cos 1-\int_1^\infty\frac{\cos x}{x^2}\text{d}x

因为1cosxx2dx\int_1^\infty|\frac{\cos x}{x^2}|\text{d}x收敛,所以去掉绝对值也收敛,因此要证明的积分也收敛.

其次证明绝对值积分发散:反证法,设收敛,而

sinxxsin2xx\frac{|\sin x|}{x}\geq\frac{\sin^2x}{x}

因此 RHS 收敛,但是

sin2xxdx=1cos2x2xdx\int\frac{\sin^2x}{x}\text{d}x=\int\frac{1-\cos2x}{2x}\text{d}x

拆分成两项,后一项可证明收敛,所以112xdx\int_1^\infty\frac{1}{2x}\text{d}x收敛,这不成立,所以假设不成立,绝对值积分发散.

瑕积分收敛性判断

艾神:建议大家看去年的课程回放.

\uparrow显然是讲不完了

计算广义积分

/Theorem/ (分部积分)

0fg=fg00fg\int_0^\infty f'g=fg|^\infty_0-\int_0^\infty fg'

/Proof/

LHS=limb(abfg)=limb(fgababfg)=RHS\text{LHS}=\lim_{b\to\infty}(\int_a^bf'g)=\lim_{b\to\infty}(fg|^b_a-\int_a^bfg')=\text{RHS}

得证.

证毕.

/Example/ (Gauss 积分)

实际上我们还不会计算,因为要用到二维的积分(唉这个实在是刻在 DNA 里面了,我觉着也没必要讲),将其作为已知条件:

ex2dx=π\int_{-\infty}^\infty e^{-x^2}\text{d}x=\sqrt{\pi}

在此基础上,求:

eax2bxcdx,(a>0)\int_{-\infty}^\infty e^{-ax^2-bx-c}\text{d}x\,,\quad(a>0)

上式得到

=ea(x+b2a)2+b24ac4adx=eb24ac4ae(a(x+b2a))2d(a(x+b2a))a=eb24ac4aey2dy/a=eb24ac4aπa\begin{aligned} &=\int_{-\infty}^\infty e^{-a(x+\frac{b}{2a})^2+\frac{b^2-4ac}{4a}}\text{d}x\\ &=e^{\frac{b^2-4ac}{4a}}\int_{-\infty}^\infty e^{-(\sqrt{a}(x+\frac{b}{2a}))^2}\frac{\text{d}(\sqrt{a}(x+\frac{b}{2a}))}{\sqrt{a}}\\ &=e^{\frac{b^2-4ac}{4a}}\int_{-\infty}^\infty e^{-y^2}\text{d}y/\sqrt{a}\\ &=e^{\frac{b^2-4ac}{4a}}\sqrt{\frac{\pi}{a}} \end{aligned}


高等微积分笔记 Lesson 29
https://physnya.top/2024/12/27/integral29/
作者
菲兹克斯喵
发布于
2024年12月27日
许可协议