欧拉积分和曲线积分是数学分析中计算技巧性非常强、概念也很抽象的部分。
本节主要讲了两类特殊的函数(用来算复杂积分的工具)和一种新的积分(在弯曲的线上做积分)。
之前学的积分,比如 ∫x2dx,都是初等函数。但很多积分算不出来具体的初等函数,于是数学家定义了两个"模具",只要积分长得像这个模具,就能直接算出答案。
工具 A:伽玛函数 Γ(s) (可以理解为“广义阶乘”)
- 模具样子:积分区间是 (0,+∞),里面是“幂函数 × 指数函数”。
Γ(s)=∫0+∞xs−1e−xdx(s>0)
- 核心性质(必须背下来):
- 阶乘性质:Γ(n+1)=n! (当 n 是整数)。例如 Γ(3)=2!=2。
- 递推公式:Γ(s+1)=sΓ(s)。
- 特殊值:Γ(21)=π (这个非常重要,做题常考)。
工具 B:贝塔函数 B(p,q)
- 模具样子:积分区间是 (0,1),里面是“x 的幂 × (1−x) 的幂”。
B(p,q)=∫01xp−1(1−x)q−1dx
- 三角形式(常用于正弦余弦的积分):
B(p,q)=2∫02πsin2p−1θcos2q−1θdθ
- A与B的关系(万能公式):
B(p,q)=Γ(p+q)Γ(p)Γ(q)
- 余元公式(用来处理奇怪的分数):
Γ(p)Γ(1−p)=sin(pπ)π(0<p<1)
- 直观理解:计算一段弯曲铁丝的质量——已知铁丝每一处的密度 ρ(x,y),沿着铁丝把所有微小的质量加起来。
- 计算方法(口诀:“一换二代三积分”):
- 换:把曲线 L 写成参数方程 {x=x(t)y=y(t)。
- 代:把被积函数里的 x,y 都换成 t。
- 算弧长微分:ds=(x′)2+(y′)2dt。
- 最后变成一个关于 t 的普通定积分:∫αβf(x(t),y(t))(x′)2+(y′)2dt。
解答原则:能化成 Γ 或 B 函数的,尽量往那上面凑。
(1) ∫01x−x2dx
- 分析:看到积分区间 (0,1) 和根号下的二次项,先提取公因式凑 B 函数定义。
- 步骤:
- 变形被积函数:x−x2=x(1−x)=x1/2(1−x)1/2。
- 这就完全符合 B(p,q) 的定义:∫01xp−1(1−x)q−1dx。
- 这里 p−1=1/2⇒p=3/2,q−1=1/2⇒q=3/2。
- 原式 =B(23,23)。
- 利用关系式:=Γ(3/2+3/2)Γ(3/2)Γ(3/2)=Γ(3)(21Γ(21))2。
- 代入 Γ(1/2)=π,Γ(3)=2!=2。
- 原式 =2(21π)2=241π=8π。
- 答案:8π
(2) ∫0π3−cosxdx
- 分析:这道题比较特殊,它是一个椭圆积分,通常无法用简单的 Γ 函数直接算出数值,除非这里有特殊技巧。但标准解法是化简为椭圆积分形式。如果是在考查欧拉积分的章节,可能存在特定参数技巧,但我们按标准步骤化简。
- 步骤:
- 利用半角公式:cosx=1−2sin22x。
- 分母变成 3−(1−2sin22x)=2+2sin22x=21+sin22x。
- 令 t=2x,则 dx=2dt。当 x=0,t=0;当 x=π,t=2π。
- 原式 =∫02π21+sin2t2dt=2∫02π1+sin2tdt。
- 注:这是一个第一类椭圆积分,通常记作 K(k) 形式,但无法算出像 π 这样的简单数。
- 修正思路:如果这是必做题,可能题目有误(例如分母是 2(1−cosx) 就能算)。但如果必须按原题算,只能保留在椭圆积分形式或使用级数展开。鉴于这是普通练习题,这道题可能是让你练习化简,或者题目抄写时有偏差(比如可能是 ∫0π/2...)。
- 我们这里给出一个最简化的积分类形式:2F(2π,i) (涉及到虚数模,通常本科阶段不要求算出具体数)。
- 建议:如果老师没有讲过椭圆积分,请检查题目是否抄错(例如是否为 ∫0π1−cosxdx? 那样能算)。
- 假设题目无误的回答:该积分为第一类完全椭圆积分形式,值为 2K(1/2)≈2.622。
(3) ∫01n1−xndx(n>1)
- 分析:把 n1−xn 处理成标准 B 函数形式。
- 步骤:
- 令 t=xn,那么 x=t1/n,dx=n1tn1−1dt。
- 区间:x∈[0,1]⇒t∈[0,1]。
- 被积函数 (1−xn)1/n1=(1−t)−1/n。
- 代入积分:∫01(1−t)−1/n⋅n1tn1−1dt=n1∫01tn1−1(1−t)−1/ndt。
- 这是 B 函数:p=n1,q=1−n1。
- 原式 =n1B(n1,1−n1)。
- 利用余元公式 Γ(p)Γ(1−p)=sinpππ:
B(p,1−p)=Γ(p)Γ(1−p)=sin(pπ)π。
- 结果 =n1⋅sin(nπ)π。
- 答案:nsinnππ
(4) ∫0+∞1+xnxm−1dx(n>m>0)
- 分析:这是 B 函数的另一种形式:∫0∞1+yya−1dy=sinaππ(其实就是余元公式的变形)。
- 步骤:
- 令 u=xn,则 x=u1/n,dx=n1un1−1du。
- 区间仍是 (0,+∞)。
- 代入:∫0+∞1+u(u1/n)m−1⋅n1un1−1du。
- 整理指数:nm−1+n1−1=nm−1。
- 原式 =n1∫0+∞1+uunm−1du。
- 利用公式 ∫0∞1+uua−1du=sinaππ,这里 a=nm。
- 因为 0<m<n,所以 0<a<1,公式成立。
- 答案:nsinnmππ
(5) ∫0+∞(1+x)24xdx
- 分析:这符合 B 函数的第三种定义:B(p,q)=∫0+∞(1+x)p+qxp−1dx。
- 步骤:
- 被积函数为 (1+x)2x1/4。
- 对比 xp−1,得 p−1=1/4⇒p=5/4。
- 对比分母 (1+x)p+q,得 p+q=2⇒5/4+q=2⇒q=3/4。
- 原式 =B(45,43)。
- 利用 B(p,q)=Γ(p+q)Γ(p)Γ(q):
=Γ(2)Γ(5/4)Γ(3/4)。
- 分母 Γ(2)=1!=1。分子用递推公式 Γ(5/4)=41Γ(1/4)。
- 原式 =41Γ(41)Γ(43)。
- 利用余元公式 Γ(z)Γ(1−z),这里 z=1/4。
Γ(1/4)Γ(3/4)=sin(π/4)π=2/2π=2π。
- 最终结果 =41⋅2π。
- 答案:42π
(6) ∫02πsin7xcos21xdx
- 分析:直接套用三角形式的 B 函数公式。
- 步骤:
- 公式:∫0π/2sinαxcosβxdx=21B(2α+1,2β+1)。
- 这里 α=7,β=1/2。
- p=27+1=4,q=21/2+1=43。
- 原式 =21B(4,43)=21Γ(4.75)Γ(4)Γ(3/4)。
- 利用递推公式 Γ(n+s)=(n+s−1)…sΓ(s) 化简分母:
Γ(19/4)=415⋅411⋅47⋅43Γ(3/4)。
- 分子 Γ(4)=3!=6。
- 代入:21⋅4415⋅11⋅7⋅3Γ(3/4)6⋅Γ(3/4)=34653⋅256=1155256。
- 答案:1155256
(7) ∫0+∞xme−xndx(m,n>0)
- 分析:Γ 函数的标准凑法。目的是把 e−xn 变成 e−t。
- 步骤:
- 令 t=xn⇒x=t1/n,dx=n1tn1−1dt。
- 代入积分:∫0+∞(t1/n)me−t⋅n1tn1−1dt。
- 整理:n1∫0+∞tnm⋅tn1−1e−tdt=n1∫0+∞tnm+1−1e−tdt。
- 根据 Γ 定义 ∫0∞ts−1e−tdt=Γ(s),这里 s=nm+1。
- 答案:n1Γ(nm+1)
(8) ∫01xp−1(1−xn)q−1dx(p,q,n>0)
- 分析:这是 B 函数的变体,需要把 (1−xn) 变成 (1−t)。
- 步骤:
- 令 t=xn⇒x=t1/n,dx=n1tn1−1dt。
- 代入:∫01(t1/n)p−1(1−t)q−1⋅n1tn1−1dt。
- 整理 t 的指数:np−1+n1−1=np−1。
- 原式 =n1∫01tnp−1(1−t)q−1dt。
- 直接对照 B 函数定义。
- 答案:n1B(np,q)
- 证明等式:
- 利用 T1(7) 的结论,令 m=0(即 x0=1)。
- 直接代入 T1(7) 结果:n1Γ(n0+1)=n1Γ(n1)。得证。
- 求极限 limn→∞n1Γ(n1):
- 令 t=n1,当 n→∞ 时,t→0+。
- 原式变为 limt→0+tΓ(t)。
- 利用性质 Γ(t+1)=tΓ(t)。
- 所以 limt→0+tΓ(t)=limt→0+Γ(t+1)=Γ(1)=0!=1。
这是含参变量积分求导的问题。
- 证明:
- Γ(s)=∫0+∞xs−1e−xdx。
- 对 s 求导,根据“积分符号下求导法则”,我们可以把导数 dsd 移到积分号里面。
- 被积函数 f(x,s)=xs−1e−x=e(s−1)lnxe−x。
- 对 s 求偏导:∂s∂(xs−1e−x)=xs−1e−x⋅lnx (指数函数求导法则 (as)′=aslna)。
- 所以 Γ′(s)=∫0+∞xs−1e−xlnxdx。
- 进一步:继续求导,n 阶导数就是乘 n 次 lnx。
- Γ(n)(s)=∫0+∞xs−1e−x(lnx)ndx。
- 直观思路:Γ(n)=(n−1)!,当 n 很大时,阶乘趋向无穷大。
- 严格证明:
- 取 x∈[1,+∞) 时,积分部分 ∫1+∞xs−1e−xdx。
- 当 s 很大时,选取一个子区间 [2,3]。在该区间上 xs−1≥2s−1。
- Γ(s)>∫232s−1e−xdx=2s−1(e−2−e−3)。
- 当 s→+∞ 时,2s−1→+∞,常数 (e−2−e−3)>0。
- 所以整体趋向 +∞。
这是一个经典的高难度积分。
- 步骤:
- 利用对称性,设 I=∫01lnΓ(x)dx。
- 令 t=1−x,则 I=∫01lnΓ(1−x)dx。
- 两式相加:2I=∫01[lnΓ(x)+lnΓ(1−x)]dx=∫01ln[Γ(x)Γ(1−x)]dx。
- 利用余元公式 Γ(x)Γ(1−x)=sinπxπ。
- 2I=∫01lnsinπxπdx=∫01(lnπ−lnsinπx)dx。
- 第一部分 ∫01lnπdx=lnπ。
- 第二部分:令 u=πx,∫01lnsinπxdx=π1∫0πlnsinudu。
- 已知经典积分 ∫0πlnsinudu=−πln2。
- 所以 2I=lnπ−(−ln2)=lnπ+ln2=ln(2π)。
- I=21ln(2π)。
- 答案:21ln(2π)
- 题目:∭Ω1+xα+yβ+zγdxdydz
- 方法:广义球面坐标变换(Dirichlet 积分技巧)。
- 令 u=xα,v=yβ,w=zγ。
- x=u1/α,dx=α1u1/α−1du。同理处理 y,z。
- 区域 Ω 变为 u,v,w>0。
- I=αβγ1∭1+u+v+wuα1−1vβ1−1wγ1−1dudvdw。
- 收敛条件:
- 分母增长阶数为 1,分子总体阶数为 α1+β1+γ1−3。
- 对于无穷远处收敛,需要总阶数小于 −1(类似三维空间 rk1 需要 k>3 收敛,这里分母是一次,所以需要分子非常小?不,其实可以转化为 Gamma 函数)。
- 该积分为 Liouville 积分形式。利用公式:
∭⋯=Γ(1/α+1/β+1/γ)Γ(1/α)Γ(1/β)Γ(1/γ)∫0∞1+tt1/α+1/β+1/γ−1dt。
- 最后的单积分要收敛,需要指数 s=α1+β1+γ1<1。
- 计算值:
- 当 α1+β1+γ1<1 时,
- I=αβγ1Γ(S)Γ(1/α)Γ(1/β)Γ(1/γ)⋅Γ(S)Γ(1−S) (其中 S=1/α+1/β+1/γ)。
- 化简得 I=αβγ1Γ(1/α)Γ(1/β)Γ(1/γ)Γ(1−(α1+β1+γ1))。
- 步骤:
- tanαx=sinαxcos−αx。
- 利用 T1(6) 的 B 函数三角形式公式。
- 原式 =21B(2α+1,2−α+1)。
- 令 p=21+α,则 1−p=1−21+α=21−α。恰好对应。
- 利用余元公式:21⋅sin(21+απ)π。
- sin(2π+2απ)=cos(2απ)。
- 得证:2cos2αππ。
例子 1.1: ∫01xa(1−xb)cdx
- 令 t=xb⇒x=t1/b。
- dx=b1tb1−1dt。
- 代入:∫01ta/b(1−t)cb1t1/b−1dt=b1∫01tba+1−1(1−t)cdt。
- 结果:b1B(ba+1,c+1)。
例子 1.2: ∫01x81−x3dx
- 对应上题,a=8,b=3,c=1/2。
- 结果:31B(38+1,21+1)=31B(3,23)。
- 计算:31Γ(9/2)Γ(3)Γ(3/2)=3127252321π2⋅21π=31161051=31516。
例子 1.3: 三叶玫瑰线面积
- 面积公式(极坐标):S=21∫ρ2dθ。
- 这里 x,y 给出的是参数方程,我们将其转为极坐标。注意到 x=ρcosθ,y=ρsinθ,对比题目 ρ=acos3θ。
- 因为有三片叶子,我们可以算一片(θ∈[−π/6,π/6])然后乘 3。或者直接利用周期性。
- θ 从 0 到 2π 会画出完整的花。
- S=21∫0π(acos3θ)2dθ (注意区间重复性,通常三叶玫瑰线 0→π 即可画完,再往后是重复。如果是 2π,可能算了两次,题目给定 2π 则需确认是否覆盖)。
- 一般计算:S=21a2∫0πcos23θdθ。
- 利用 Walli's 公式或倍角公式,∫0πcos2kxdx=π/2。
- 结果 S=21a2⋅2π=4πa2。
(1) 三角形 O(0,0),A(1,0),B(0,1)
- 积分 ∫L(x+y)ds 分三段算。
- OA(y=0,x:0→1): ds=dx。∫01xdx=1/2。
- AB(y=1−x,x:1→0): 直线 y=−x+1,ds=12+(−1)2dx=2∣dx∣。注意第一类曲线积分下限小上限大。或者直接参数化 x=t,y=1−t,t:1→0 是错的,应该用长度参数。简单说,把线段拉直积分。
在 AB 上,函数值 x+y=x+(1−x)=1。线段长是 2。积分 = 1×2=2。
- BO(x=0,y:1→0): ds=dy。函数值 x+y=y。∫01ydy=1/2。
- 总和:1/2+2+1/2=1+2。
(2) 单位圆 x2+y2=1,求 ∫L∣y∣ds
- 参数化:x=cost,y=sint,t∈[0,2π]。
- ds=1dt。
- 原式 =∫02π∣sint∣dt。
- 利用对称性,=4∫0π/2sintdt=4[−cost]0π/2=4(0−(−1))=4。
(3) 星形线 x2/3+y2/3=a2/3,求 ∫L∣x∣1/3ds
- 参数化:x=acos3t,y=asin3t,t∈[0,2π]。
- 算 ds:
x′=−3acos2tsinty′=3asin2tcostds=x′2+y′2dt=3a∣sintcost∣dt。
- 利用对称性,只算第一象限 (t:0→π/2) 乘 4。
- 函数 ∣x∣1/3=(acos3t)1/3=a1/3cost。
- 积分 =4∫0π/2a1/3cost⋅3asintcostdt=12a4/3∫0π/2cos2tsintdt。
- 令 u=cost,=12a4/3∫01u2du=12a4/3⋅31=4a4/3。
(4) 双纽线 (x2+y2)2=x2−y2,求 ∫L∣x∣ds
- 极坐标方程:r2=cos2θ。对称性,只算第一象限半个叶子 (θ:0→π/4) 乘 4。
- ds 公式(极坐标):r2+(r′)2dθ。
2rr′=−2sin2θ⇒r′=−rsin2θ。
ds=r2+r2sin22θdθ=r2r4+sin22θdθ=r2cos22θ+sin22θdθ=r1dθ。
或者更简单的:双纽线弧长微分 ds=cos2θ1dθ。
- ∣x∣=rcosθ=cos2θcosθ。
- 积分 =4∫0π/4(cos2θcosθ)⋅cos2θ1dθ=4∫0π/4cosθdθ。
- =4[sinθ]0π/4=4⋅22=22。
(5) 螺旋线 x=acost,y=asint,z=bt
- x′=−asint,y′=acost,z′=b。
- ds=a2+b2dt (常数)。
- 被积函数 x2+y2+z2=a2+(bt)2。
- 积分 ∫02π(a2+b2t2)a2+b2dt。
- =a2+b2[a2t+3b2t3]02π=a2+b2(2πa2+38π3b2)。
(6) 曲线 x=t,y=322t3/2,z=21t2
- x′=1,y′=2t1/2,z′=t。
- ds=1+2t+t2dt=(1+t)2dt=(1+t)dt (因为 t>0)。
- 被积函数 xyz=t⋅322t3/2⋅21t2=32t9/2。
- 积分 ∫0132t9/2(1+t)dt=32∫01(t9/2+t11/2)dt。
- =32(112+132)=32214313+11=32214324=143162。
(7) 球面 x2+y2+z2=a2 与平面 x+y+z=0 交线
- 这是一个过原点的平面截球面,得到的是大圆,半径为 a。
- 被积函数 (xy+yz+zx)。
- 利用代数恒等式:(x+y+z)2=x2+y2+z2+2(xy+yz+zx)。
- 在曲线上,满足 x+y+z=0 和 x2+y2+z2=a2。
- 代入:0=a2+2(xy+yz+zx)。
- 所以被积函数 xy+yz+zx=−2a2 (常数!)。
- 积分 =∫L(−2a2)ds=−2a2⋅(大圆周长)。
- 周长 =2πa。
- 结果 =−2a2⋅2πa=−πa3。
- 参数化:x=acost,y=bsint。
- ds=a2sin2t+b2cos2tdt。
- 质量 M=∫02π∣bsint∣a2sin2t+b2cos2tdt。
- 利用对称性,上下左右对称,算第一象限 (0→π/2) 乘 4。
- M=4∫0π/2bsinta2sin2t+b2cos2tdt。
- 技巧:把根号里的 sin2t 换成 1−cos2t。
=4b∫0π/2sinta2−(a2−b2)cos2tdt。
- 令 u=cost,du=−sintdt。
=4b∫01a2−(a2−b2)u2du。
- 设 c2=a2−b2 (若 a>b)。积分形式 ∫a2−c2u2du。
- 这是一个标准积分:∫A2−x2dx=2xA2−x2+2A2arcsinAx。
- 稍微复杂,但这道题如果 a=b (圆) 就很简单。如果 a=b,结果含 arcsin 或对数。
- 计算结果(假设 a>b):
M=2b2+a2−b22ab2arcsinaa2−b2。
(如果老师没要求这么细,列出积分式子通常也能拿大部分分)。