跳转至

6.定积分的应用

总结

  • 旋转体的模型:圆柱体 \(\Delta V = \pi \Delta (r^2 \cdot h)\)
    • f(x) 绕 y = c 旋转 \(\implies\) \(\Delta V \approx \pi r^2 \Delta h \implies \int \pi (f_x-c)^2~d_x\)
    • f(x) 绕 x = c 旋转 \(\implies\) \(\Delta V \approx \pi \Delta r^2 h \approx 2\pi r \Delta r h \implies \int 2\pi (x-c)f_x~d_x\)
    • 对于 \(f(x)\) 的“对称函数” \(f^{-1}(y)\),偷偷地把 x 和 y 互换 (通常要用到这种方法往往是因为关于 f(x) 的被积函数难以求解)
    • 注:正确性不能保证;对于一道题,最多只有两种解法,不能更多了
  • 平面面积的模型:矩形 \(\Delta S = \Delta (wh)\)
    • f(x) 和 x轴 围成的面积:
    • f(x) 和 g(x) 围成的面积:
    • f(x) 和 直线g(x)=mx+b 以及 两个垂直于g(x)的直线 围成的面积:\(\frac 1{1+m^2}\int_a^b [f(x) - mx - b][1+mf'(x)]~d_x\) (详见 p574)

本章讨论定积分计算曲线间的面积,立体体积,压力做功 的应用

这些方法与求曲线下的面积相似:现将数据 Q 分割成若干份,然后\(\Delta x f(x_i^*)\) 近似每一部分,利用黎曼和求近似 Q 再求极限,将Q表示成积分的形式。最后利用 微积分基本定理 和 langrange中值定理 求积分的估计值

Note

  • 原理都一样,最关键的莫过于 \(\Delta x f(x_i)\),然后给 f 和 x 赋予某些含义

1.两条曲线间的面积

Warning

  • 面积 是一个正值,而非 "面积的代数和"

计算 f(x) 和 g(x) 之间的面积,其黎曼和为:\(A = \sum\limits_{i=1}^n \Delta x [f(x_i^*) - g(x_i^*)]\)

两条曲线间的面积

A 是曲线 y=f(x),y = g(x) 和 直线 x = a,x = b 围成的区域的面积,其中 f 和 g 在区间 [a, b] 上连续,且 f(x) \(\ge\) g(x),则:

\(A = \int_a^b [f(x) - g(x)] d_x\)

\(f(x) - g(x)\) 可正可负,则 \(A = \int_a^b |f(x) - g(x)| d_x\)

例子

  1. 假设 \(S=\{y=e^x,y=x,x=0,x=1\}\),那么 \(\|S\|=\int_0^1|e^x-x|d_x=\left[e^x-\frac12x^2\right]_0^1=e-\frac32\)
  2. 假设 \(S=\{y=x^2,y=2x-x^2\}\),求解方程 \(y=x^2,y=2x-x^2\)\(x=0,x=1\),于是 \(\|S\|=\int_0^1|x^2-(2x-x^2)|d_x=\int_0^1(2x-2x^2)d_x=\Big[x^2-\frac23x^3\Big]_0^1=\frac13\)
  3. 假设 \(S=\{y=x/\sqrt{x^2+1},y=x^4-x\}\),求解方程 \(x/\sqrt{x^2+1}=x^4-x\)\(x=0,x\approx1.18\),于是 \(\|S\|\approx\int_0^{1.18}|x/\sqrt{x^2+1}-(x^4-x)|d_x=\int_0^{1.18}|x/\sqrt{x^2+1}-(x^4-x)|d_x=\left[\sqrt{x^2+1}-\frac15x^5+\frac12x^2\right]_0^{1.18}\approx0.785\)
  4. 假设 \(S=\{y=\sin x,y=\cos x,x=0,x=\frac\pi2\}\),求解方程 \(\sin x=\cos x,0\le x\le\frac\pi2\)\(x=\frac\pi4\),于是 \(\|S\|=\int_0^{\frac\pi2}|\sin x-\cos x|d_x=\left(-\int_0^{\frac\pi4}+\int_{\frac\pi4}^{\frac\pi2}\right)(\sin x-\cos x)d_x=-\left[-\cos x-\sin x\right]_0^{\frac\pi4}+\left[-\cos x-\sin x\right]_{\frac\pi4}^{\frac\pi2}=2\sqrt2-2\)
  5. 假设 \(S=\{y=x-1,y^2=2x+6\}\),求解方程 \(y=x-1,y^2=2x+6\)\((x,y)=(-1,-2)\)\((x,y)=(5,4)\)\(\|S\|=\int_{-2}^4\left|(y+1)-\frac12(y^2-6)\right|d_y=\int_{-2}^4(-\frac12y^2+y+4)d_y=\left[-\frac16y^3+\frac12y^2+4y\right]_{-2}^4=18\)

Tip

  • 可以注意到,上述的面积并不严格等于 f(x) 和 g(x) 之间的封闭区域面积
  • 计算封闭区域面积:
    1. f(x),g(x) 均是 单值函数
      • 这种情况可以用两种方法求 (y 是 x 的函数 或 x 是 y 的函数)
      • 计算 \(h = f - g\)\(h = f^{-1} - g^{-1}\) 的零点,并判断 左右 或 上下 两边的正负 (这里的 逆函数 是形式上的,数学定义上这是不正确的说法,这里是把 x 和 y 互换了)
      • 计算 最左端 和 最右端 的零点之间区间的封闭面积
      • 假设区间为 [a, b],则面积为 \(\int_a^b |y_T - y_B| d_x\)\(\int_a^b |x_R - x_L| d_y\)
    2. f(x) 是 单值函数,而 g(x) 是 多值函数 (比如 函数方程含有 关于y的不小于2次的多项式)
      • 只有一种合适的求法:x 是 y 函数 (除非 y 是 x 的函数时,分类讨论)
      • 计算 \(h = f^{-1} - g^{-1}\) 的零点
      • 计算 最左端 和 最右端 的零点之间区间的封闭面积
      • 假设区间为 [a, b],则面积为 \(\int_a^b |x_R - x_L| d_y\)
    3. f(x),g(x) 均是 多值函数

总结

  1. 曲边梯形:\(S=\{(x,y)|~x\in[a,b],\min\{f(x),g(x)\}\le y\le\max\{f(x),g(x)\}\}\) 定义为曲边梯形,简记为 \(\{y=f(x),y=g(x),a\le x\le b\}\);区域 S 的面积记为 \(\|S\|\)
  2. 曲边梯形的面积:假设曲边矩形 \(S=\{x\in[a,b],y=f(x),y=g(x)\}\),那么其面积为 \(\|S\|=\int_a^b|f(x)-g(x)|d_x\)
  3. 函数封闭区域:\(y=f(x)\)\(y=g(x)\) 之间围成的区域称为函数封闭区域,记为 \(S=\{y=f(x),y=g(x)\}\)
  4. 函数封闭区域的面积:假设 \(S=\{y=f(x),y=g(x)\}\),若 \(\forall i=1..n,[a_i,b_i]\)\(y=f(x),y=g(x)\) 的封闭区间,那么 S 的面积为 \(\|S\|=\sum\limits_{i=1}^n\int_{a_i}^{b_i}|f(x)-g(x)|d_x\)
  5. 曲线封闭区域:一般曲线 \(f(x,y)=0\)\(g(x,y)=0\) 围成的区域称为曲线封闭区域,记为 \(S=\{f(x,y)=0,g(x,y)=0\}\)

一级结论

  1. 假设曲边矩形 S 为 \(\{x\in[a,b],y=f(x),y=g(x)\}\),若 \(f,g\)\([a,b]\) 上单射,那么 S 可记为 \(S=\{\min\{f^{-1}(a),f^{-1}(b)\}\le y\le\max\{f^{-1}(a),f^{-1}(b)\},x=f^{-1}(y),x=g^{-1}(y)\}\)

2.体积

将我们要测量体积的立体穿过一条轴线,我们通常取其为 x 轴,用一个垂直于x轴平面 S 截断 或 "扫描" 立体,得到 横截面 的面积 A,实际上这个面积是关于 x 的函数,记为 A(x)

我们用n-1个平面 \(P_{x_i}\) 将立体(或 要测量的区间 [a, b] 之间) 划分为 n 个 高 (或 宽) 为 \(\Delta x\) 的 区域 或 "薄片",在此基础上再讨论体积

Note

  • 计算面积时我们 "切条",计算体积时我们 "切片",计算??时我们 "??"

体积定义

令 S 是介于 \(x=a\)\(x=b\) 固体,如果经过点 x 垂直于 x 轴的平面 \(P_x\) 在 S 上的横截面的面积为 \(A(x)\),其中 A 为连续函数,那么体积为:

\(V = \lim\limits_{n\to +∞} \Delta x A(x_i^*) = \int_a^b A(x) d_x\)

例子

  1. 计算半径为 r 的球体体积:球体记为 \(V=\{-r\le x\le r,A(x)=\pi(r^2-x^2)\}\)(利用几何和代数计算 \(A(x)\)),于是 \(\|V\|=\int_{-r}^r\pi(r^2-x^2)d_x=\left[\pi(r^2x-\frac13x^3)\right]_{-r}^r=\frac43\pi r^3\)
  2. 计算 \(y=\sqrt x\) 从 0 到 1 关于 x 轴旋转得到的区域体积:该区域为旋转立体 \(V=\{a\le x\le b,A(x)=\pi(\sqrt x)^2\}\),于是 \(\|V\|=\int_0^1\pi|x|d_x=\frac12\pi x^2\Big|_0^1=\frac12\pi\)
  3. 计算 \(y=x^3,y=8,x=0\) 围成的区域绕 y 轴旋转得到的立体体积:该区域为旋转立体 \(V=\{0\le y\le 8,A(y)=\pi(\sqrt[3]y)^2\}\),于是 \(\|V\|=\int_0^8\pi y^{2/3}d_y=\left[\frac35\pi y^{5/3}\right]_0^8=\frac{96}5\pi\)(注意区分区域 \(\{0\le y\le 8,A(y)=\pi|(\sqrt[3]y)^2-2^2|\}\),其体积为 \(\|V\|=\int_0^8\pi(4-y^{2/3})d_y=\left[4\pi y-\frac35\pi y^{5/3}\right]_0^8=\frac{64}5\pi\)
  4. 计算函数封闭区域 \(\{y=x,y=x^2\}\) 绕 x 轴旋转得到的旋转立体 V 的体积:求解方程 \(x=x^2\)\(x=0,x=1\),那么 \(V=\{0\le x\le 1,A(x)=\pi|x^2-(x^2)^2|\}\),于是 \(\|V\|=\int_0^1\pi|x^2-x^4|d_x=\left[\frac13\pi x^3-\frac15\pi x^5\right]_0^1=\frac2{15}\pi\)
  5. 计算函数封闭区域 \(\{y=x,y=x^2\}\)\(y=2\) 旋转得到的一般旋转立体 V 的体积:\(V=\{x\in[0,1],A(x)=\pi|(x-x^2)(x+x^2-4)|\}\)\(\|V\|=\int_0^1\pi|(x-x^2)(x+x^2-4)|d_x=\pi\int_0^1(x^4-5x^2+4x)d_x=\pi\left[\frac15x^5-\frac53x^3+2x^2\right]_0^1=\frac8{15}\pi\)
  6. 计算函数封闭区域 \(S=\{y=x,y=x^2\}\)\(x=1\) 旋转得到的旋转立体 V 的体积:根据[6.1一级结论],S 可变换为 \(S=\{y\in[0,1],x=y,x=\sqrt y\}\),蕴含 \(V=\{y\in[0,1],A(y)=\pi|(y-\sqrt y)(y+\sqrt y-2)|\}\),于是 \(\|V\|=\int_0^1\pi|(y-\sqrt y)(y+\sqrt y+2)|d_y=\pi\int_0^1(y^2+y-2\sqrt y)d_y=\pi\left[\frac13y^3+\frac12y^2-\frac13y^{3/2}\right]_0^1=\frac12\pi\)
  7. 计算底面为半径为 r 的圆,垂直于底面的横截面为等边三角形的立体 V 的体积:在底面上建立 \(xoy\)\(\forall x_0\in[-r,r]\)\(x=x_0\) 与 V 相交的等边三角形边长为 \(d=2\sqrt{r^2-x^2}\),蕴含其面积为 \(A(x)=\frac{\sqrt3}4d^2=\frac{\sqrt3}4(2\sqrt{r^2-x^2})^2=\sqrt3(r^2-x^2)\),于是 \(\|V\|=\int_{-r}^r\sqrt3(r^2-x^2)d_x=2\int_0^r\sqrt3(r^2-x^2)d_x=2\sqrt3\left[r^2x-\frac13x^3\right]_0^r=\frac{4\sqrt3}3r^3\)(注意到 \(A(x)\) 的偶函数)
  8. 计算底边长为 L,高为 h 的棱锥 V 的体积:\(\forall x_0\in[0,h]\)\(x=x_0\) 与 V 相交的正方形边长的一半 d 满足 \(\frac h{L/2}=\frac x{d}\),蕴含 \(d=\frac{xL}{2h}\),于是 \(A(x)=(2d)^2=\frac{x^2L^2}{h^2}\),于是 \(\|V\|=\int_0^h\frac{L^2}{h^2}x^2d_x=\frac{L^2}{h^2}\left[\frac13x^3\right]_0^h=\frac13L^2h\)
  9. (锲的体积)假设足够高的圆柱体 W 半径为 r,平面 \(V_1\) 垂直于 W 的轴,平面 \(V_2\)\(\theta\) 角交于 \(V_1\cap W\) 的直径,计算 \(W,V_1,V_2\) 的封闭区域 V 的体积:
    1. 法1:\(\forall x_0\in[-r,r]\)\(x=x_0\) 与 V 截得一个底边为 \(d=\sqrt{r^2-x^2}\) 高为 \(d\tan\theta\) 的直角三角形,蕴含 \(A(x)=\frac12d^2\tan\theta=\frac12(r^2-x^2)\tan\theta\),于是 \(\|V\|=\int_{-r}^r\frac12(r^2-x^2)\tan\theta~d_x=\frac12\tan\theta\left[r^2x-\frac13x^3\right]_{-r}^r=\frac23r^3\tan\theta\)
    2. 法2:\(\forall x_0\in[0,r]\)\(x=x_0\) 与 V 截得一个底边为 \(d=2\sqrt{r^2-x^2}\) 高为 \(x\tan\theta\) 的矩形,蕴含 \(A(x)=2\sqrt{r^2-x^2}x\tan\theta\),于是 \(\int_0^r2\sqrt{r^2-x^2}x\tan\theta~d_x=2\tan\theta\int_{r^2}^0\sqrt{r^2-x^2}x/(-2x)~d(r^2-x^2)=-\tan\theta\int_{r^2}^0\sqrt t~d_t=\tan\theta\left[\frac23t^{3/2}\right]_0^{r^2}=\frac23r^3\tan\theta\)

Tip

  • “体积”这个概念是三维的,而“一维积分”只能处理 “一个自变量x和其函数 \(f(x)\) 的乘积” 的 “两维乘法问题”,此时就需要将 \(f(x)\) 复合一下:\(f=g\circ h\)
  • 在笛卡尔坐标系中
    1. 选定一个合适的x轴来 迭代 立体的每个横截面(坐标中心与立体的相对位置要合适),然后 y 轴迭代 当前横截面的每条线;于是y轴与线长度函数积分得到当前横截面面积,然后 z轴 与 横截面面积函数 积分得到立体体积 (\(\int S(x) d_x = \int \int L(x, y)~d_y~d_x\)) (并不是划分得越细越好)
    2. 选定一个合适的x轴来 迭代 立体的每个横截面,但是此时横截面的面积直接依赖于x (横截面的面积是x的函数);于是x轴与横截面面积函数积分得到立体体积 (\(\int S(x) d_x\))
    3. 但是有时候立体体积直接依赖于r(如球体体积)
  • 旋转问题:
    • \(f(x)\)\(g(x)\) 围成n个封闭区域(有n+1零点),再将其绕 x 轴旋转,求生成的立体体积:讨论某个封闭区域 \(x\in[a, b]\),那么横截面函数为 \(S(x) = \pi |g(x)^2 - f(x)^2|\),于是立体体积为 \(V = \int_a^b \pi |g(x)^2 - f(x)^2| d_x\)
    • y=(f-g)(x) 绕 y=c 轴旋转:立体体积为 \(V = \int_a^b \pi |(c - g(x))^2 - (c - f(x))^2| d_x\)
    • y=(f-g)(x) 绕 x=c 轴旋转:立体体积为 \(V = \int_a^b \pi |(c - g^{-1}(x))^2 - (c - f^{-1}(x))^2| d_x\)

练习

  1. 旋转体
    1. \(y=\sqrt x\) 从 x=0..1 旋转所得立体体积
    2. \(y=x^3, y=8, x=0\) 绕 y 轴旋转所得立体体积 (如何通过先x切割横截面得到体积?)
    3. \(y=x, y=x^2\) 分别求 绕 x 轴旋转,绕 y=2 旋转,绕 x=-1 旋转 时所得体积
  2. 非旋转体
    1. 横截面为三角形
    2. 棱锥体积:底边长为 L,高为 H
    3. 锲的体积:一个平面垂直于圆柱的轴 并 切割圆柱,另一个平面以 \(\theta\) 角切割圆柱 并且两平面交线与圆柱截面半径的半径相交

总结

  1. 曲面椎体:假设一个立体的上下表面是平面,侧面是曲面,若 x 轴垂直于上下表面,横截面积函数 \(A(x)\) 是连续的,那么该立体称为曲面立方,记为 \(V=\{(x,y,z)|~a\le x\le b,A(x)=\|S(x)\|\}\)
  2. 曲面椎体的体积:\(V=\{(x,y,z)|~a\le x\le b,A(x)=\|S(x)\|\}\) 的体积为 \(\|V\|=\int_a^bA(x)d_x\)
  3. 旋转立体(标准曲面椎体):假设一个曲边梯形为 \(S=\{x\in[a,b],y=f(x),y=g(x)\}\),立体 V 由 S 旋转得到
    1. 若 V 由 S 绕 x 轴旋转一周而成,那么 V 称为关于 x 轴的旋转立体,记为 \(V=\{x\in[a,b],A(x)=\pi|f^2(x)-g^2(x)|\}\)
    2. 若 V 由 S 绕 \(y=c\) 旋转一周而成,那么 V 称为关于 \(y=c\) 的旋转立体一般旋转立体,记为 \(V=\{x\in[a,b],A(x)=\pi|(f(x)-c)^2-(g(x)-c)^2|\}\)\(V=\{x\in[a,b],A(x)=\pi|(f(x)-g(x))(f(x)+g(x)-2c)|\}\)
    3. 若 V 由 S 绕 y 轴旋转一周而成,那么 V 称为关于 y 轴的旋转立体,记为 \(\{x\in[a,b],A(x)=2\pi|x|\cdot|f(x)-g(x)|\}\)
    4. 若 V 由 S 绕 \(x=c\) 旋转一周而成,那么 V 称为关于 \(x=c\) 的旋转立体,记为 \(\{x\in[a,b],A(x)=2\pi|x-c|\cdot|f(x)-g(x)|\}\)

注:旋转轴不应与区域 S 相交,否则体积会重复计算;(3.1)和(3.2)是第一类旋转,(3.3)和(3.4)是第二类旋转

一级结论

  1. 若关于 x 轴的曲边矩形 \(S=\{a\le x\le b,y=f(x),y=g(x)\}\) 中,\(f,g\)\([a,b]\) 上单射,那么 S 是可绕 \(x=c\) 旋转的

二级结论

  1. 假设 \(xoy\) 上的区域 S 与 \(l:ax+by+c=0\) 不相交,若将 S 绕 l 旋转得到曲面椎体 V,那么 \(\|V\|=\iint_S2\pi\frac{|ax+by+c|}{\sqrt{a^2+b^2}}d_\sigma\)(参考:曲线的旋转立体的体积

三级结论

  1. 斜边为 a 的等边三角形面积为 \(\frac{\sqrt3}4a^2\)
  2. 斜边为 a 的等边直角三角形面积为 \(\frac14a^2\)

练习

  1. 证明:内外半径分别为 \(r,R\) 的圆环 V 的体积为 \(\|V\|=2\pi^2r^2R\)
  2. 若立体 V 以椭圆 \(\frac{x^2}{a^2}+\frac{y^2}{b^2}=1\) 为底面,对于 \(\forall c\in[-a,b]\)\(x=c\) 与 V 交于一个等边直角三角形,证明 \(\|V\|=\frac43ab^2\)
  3. 证明:半径为 r 的球体中,高为 h 的球帽体积为 \(\pi h^2(r-\frac13h)\)
  4. 证明:高为 h,上下底半径分别为 \(R,r\) 的正圆锥的平截头体 V 的体积为 \(\frac13\pi h(R^2+Rr+r^2)\)
  5. 假设下列积分表示立体体积,反推这些立体 V 的表达式
    1. \(\pi\int_0^{\pi/2}\cos^2x~d_x\)
    2. \(\pi\int_2^5y~d_y\)
    3. \(\pi\int_0^1(y^4-y^8)d_y\)
    4. \(\pi\int_0^{\pi/2}[(1+\cos x)^2-1^2]d_x\)
  6. 证明:两个半径均为 r 的圆柱正交的公共部分 V 的体积为 \(\frac{16}3r^3\)
  7. 证明:两个半径均为 r 的球体相交,并且球心在另一个球体表面,两球的公共部分 V 的体积为 \(\frac54\sqrt3+4\pi\)

提示

(1)

法1:V 是由区域 \(\{y=R+\sqrt{r^2-x^2},y=R-\sqrt{r^2-x^2}\}\) 绕 x 轴旋转得到的,即 \(V=\{x\in[-r,r],A(x)=\pi|(R+\sqrt{r^2-x^2})^2-(R-\sqrt{r^2-x^2})^2|\}\)

其中 \(A(x)=\pi(2R)(2\sqrt{r^2-x^2})=4\pi R\sqrt{r^2-x^2}\)

于是 \(\|V\|=4\pi R\int_{-r}^r\sqrt{r^2-x^2}d_x=4\pi R\cdot\frac12\pi r^2=2\pi^2r^2R\)

法2(柱面法):V 是由区域 \(\{x\in[R-r,R+r],y=0,y=\pm\sqrt{r^2-(R-x)^2}\}\) 绕 y 轴旋转得到的,即 \(V=\{x\in[R-r,R+r],A(x)=2\pi x\cdot2\sqrt{r^2-(R-x)^2}\}\)

于是 \(\|V\|=4\pi\int_{R-r}^{R+r}x\sqrt{r^2-(R-x)^2}d_x=-4\pi\int_r^{-r}x\sqrt{r^2-(R-x)^2}d(R-x)\) \(=4\pi\int_{-r}^r(R-t)\sqrt{r^2-t^2}d_t=4\pi R\int_{-r}^r\sqrt{r^2-t^2}d_t=4\pi R\cdot\frac12\pi r^2=2\pi^2r^2R\)

(2) \(\forall x_0\in[-a,a]\)\(x=x_0\) 与 V 交于斜边长为 \(\left(2\frac ba\sqrt{a^2-x^2}\right)\) 的等边直角三角形,其面积为 \(A(x)=\frac14\left(2\frac ba\sqrt{a^2-x^2}\right)^2=\frac{b^2}{a^2}(a^2-x^2)\)

于是 \(V=\{x\in[-a,a],A(x)=\frac{b^2}{a^2}(a^2-x^2)\}\),那么 \(\|V\|=\frac{b^2}{a^2}\int_{-a}^a(a^2-x^2)d_x=2\frac{b^2}{a^2}\left[a^2x-\frac13x^3\right]_0^a=\frac43ab^2\)

(3)(绕 x 轴)

注意到球帽由区域 \(\{x\in[r-h,r],y\in[0,\sqrt{r^2-x^2}]\}\) 绕 x 轴旋转所得

于是 \(\|V\|=\int_{r-h}^r\pi(r^2-x^2)d_x=\pi\left[r^2x-\frac13x^3\right]_{r-h}^r\)

\(=\pi\left[r^2(r-(r-h))-\frac13(r^3-(r-h)^3)\right]=\pi[r^2h-\frac13h(r^2+r(r-h)+(r-h)^2)]=\pi h^2(r-\frac13h)\)

(4) 记 u 为平截头体的母线的延长线与圆锥轴交点 P 到 V 的上表面的距离,满足 \(\frac ur=\frac{u+h}R\),即 \(u=\frac{hr}{R-r}\)

V 为区域 \(\{x\in[u,u+h],y=0,y=\frac rux\}\) 绕 x 轴旋转所得的立体,

于是 \(\|V\|=\int_u^{u+h}\pi(\frac rux)^2d_x=\pi\frac{r^2}{u^2}\left[\frac13x^3\right]_u^{u+h}=\frac13\pi\frac{r^2}{u^2}\cdot h(3u^2+3uh+h^2)=\frac{\pi r^2h}3(3+3h/u+h^2/u^2)\)

\(=\frac{\pi r^2h}3(3+3\frac{R-r}r+\frac{(R-r)^2}{r^2})=\frac13\pi h(3Rr+R^2-2Rr+r^2)=\frac13\pi h(R^2+Rr+r^2)\)

(5) (以下旋转均是第一类旋转)

  1. \(V=\{x\in[0,\pi/2],y=0,y=\cos^2x\}\)\(\{y=0,y=\cos^2x\}\) 绕 x 轴旋转)
  2. \(V=\{y\in[2,5],x=0,x=\sqrt y\}\)\(\{x=0,x=y\}\) 绕 y 轴旋转)
  3. \(V=\{y\in[0,1],x=y^4,x=y^8\}\)\(\{x=y^4,x=y^8\}\) 绕 y 轴旋转)
  4. \(V=\{y\in[0,\pi/2],y=1+\cos x,y=1\}\)\(\{y=1+\cos x,y=1\}\) 绕 x 轴旋转)

(6) \(\forall x\in[-r,r]\)

V 的俯视图形是两个正交圆柱的侧面图形的相交的封闭图形,即一个矩形,其边长计算为 \(2\sqrt{r^2-x^2}\)

蕴含 \(A(x)=(2\sqrt{r^2-x^2})^2=4(r^2-x^2)\)

于是 \(\|V\|=\int_{-r}^r4(r^2-x^2)d_x=8\int_0^r(r^2-x^2)d_x=\frac{16}3r^3\)

(7) V 由 \(S_1=\{x\in[0,r/2],y=0,y=\sqrt{r^2-(x+r/2)^2}\}\)\(S_2=\{x\in[-r/2,0],y=0,y=\sqrt{r^2-(x-r/2)^2}\}\) 绕 x 轴旋转所得,

而其中 \(S_1\)\(S_2\) 对称,于是 \(\|V\|=2\int_0^{r/2}\pi[r^2-(x+r/2)^2]d_x=\frac5{12}\pi r^3\)

3.柱面法求体积

如果 f(x) 绕 y 轴旋转,利用平面切割立体时难以计算横截面的面积,此时可用考虑 柱体薄壳法 来解决

柱体薄壳法 通过将 \(x_{i-1}\)\(x_{i}\) 之间的区域绕y轴旋转 得到 "薄壳",此时高度设为 \(h=f(x_i^*)\),计算此薄壳体积为:

\(V(x_i) = \pi x_{i}^2 f(x_i^*) - \pi x_{i-1}^2 f(x_i^*) = \pi (x_{i}^2 - x_{i-1}^2)f(x_i^*) = \pi (x_i-x_{i-1})(x_i+x_{i-1})f(x_i^*)\)

我们构造一个样本点 \(x_i^* = \frac {x_i-x_{i-1}}2 = \overline x_i\);而根据定义 \(\Delta x = x_i - x_{i-1}\),得到

\(V(x_i) = 2\pi \cdot \Delta x \cdot \overline x_i f(\overline x_i)\),黎曼和的极限:\(\lim\limits_{n\to +∞}\sum\limits_{i=1}^n 2\pi \cdot \Delta x \cdot \overline x_i f(\overline x_i) = \int_a^b 2\pi \cdot x f(x)~d_x\)

Tip

  • 实际上样本点取 左侧\(x_{i-1}\) 或 右侧\(x_i\),结果都与 中值 一样
  • 与上一节 通过平面切割立体 的方法不同,本节使用的是 柱体薄壳(一个立体) 来切割立体

柱体薄壳法

由曲线 y=f(x) 从 a 到 b 下的区域绕 y 轴旋转所得立体体积为 \(V = \int_a^b 2\pi \cdot x f(x)~d_x\) & \(0\le a < b\)

例子

  1. 计算区域 \(\{y=0,y=2x^2-x^3\}\) 绕 y 轴旋转所得立体 V 的体积:\(V=\{x\in[0,2],A(x)=2\pi x(2x^2-x^3)\}\),于是 \(\|V\|=\int_0^22\pi x(2x^2-x^3)d_x=2\pi\left[\frac12x^4-\frac15x^5\right]_0^2=\frac{16}5\pi\)
  2. 计算区域 \(\{y=x,y=x^2\}\) 绕 y 轴旋转所得立体 V 的体积:\(V=\{x\in[0,2],A(x)=2\pi x(x-x^2)\}\),于是 \(\|V\|=\int_0^12\pi x(x-x^2)d_x=2\pi\int_0^1(x^2-x^3)d_x=2\pi\left[\frac13x^3-\frac14x^4\right]_0^1=\frac16\pi\)
  3. 计算区域 \(\{x\in[0,1],y=0,y=\sqrt x\}\) 绕 x 轴旋转所得立体 V 的体积:
    1. 法1:\(V=\{x\in[0,1],A(x)=\pi(\sqrt x)^2\}\),于是 \(\|V\|=\int_0^1\pi xd_x=\frac12\pi\)
    2. 法2:\(V=\{y\in[0,1],A(y)=2\pi|y|\cdot|y^2-0|\}\),于是 \(\|V\|=\int_0^12\pi y\cdot y^2d_y=\frac12\pi\)
  4. 计算区域 \(\{y=0,y=x-x^2\}\)\(x=2\) 旋转所得立体的体积:\(V=\{x\in[0,1],A(x)=2\pi|x-2|\cdot|(x-x^2)-0|\}\),于是 \(\|V\|=\int_0^12\pi|x-2|\cdot(x-x^2)d_x=2\pi\int_0^1(x^3-3x^2+2x)d_x=2\pi\left[\frac14x^4-x^3+x^2\right]_0^1=\frac12\pi\)

Warning

  • 对于绕 x=c 的情况,就有必要重新建模了,假设 \(c\ge b\)
    • \(\displaystyle V(x_i) = \pi((c-x_{i-1})^2 - (c-x_i)^2)f(x_i^*) = \pi[2c(x_i - x_{i-1}) + x_{i-1}^2-x_i^2]f(x_i^*) = \pi (x_i-x_{i-1})[2c-(x_i + x_{i-1})]f(x_i^*)\)
    • \(\displaystyle = 2\pi\Delta x\cdot (c-\frac {x_i + x_{i-1}}2) f(x_i^*) = 2\pi\Delta x\cdot (c-\overline x_i) f(\overline x_i)\)
    • 因此 \(\displaystyle \lim\limits_{n\to +∞}\sum\limits_{i=1}^n 2\pi\Delta x\cdot (c-\overline x_i) f(\overline x_i) = \int_a^b 2\pi(c-x)f(x)d_x\)
  • 对于 y=f(x) 绕 x=c 旋转,体积为 \(\int_a^b 2\pi(c-x)f(x)d_x\)

Info

  • 柱面法的本质:\(V_i = \lim\limits_{\Delta h\to 0} \pi(r+\Delta r)^2(h+\Delta h)-\pi r^2h = \lim\limits_{\Delta h\to 0}\pi[(2r\Delta r+\Delta r^2)h + (r^2+2r\Delta r+\Delta r^2)\Delta h]\)
  • \(= 2\pi r\Delta r h+\pi \Delta r^2h = 2\pi \Delta rh(r+\frac12\Delta r)\approx 2\pi \Delta rhr\) (取 \(r\approx \overline r=r+\frac12\Delta r\)

总结

  1. 参见[6.2总结]

练习

  1. 使用薄壳法(柱面法)计算 \(y=e^{-x^2},y=0\),x=0,x=1$ 所围成的区域绕 y 轴旋转所得立体的体积
  2. 计算区域 \(\{y=x^2,y=0,x=1,x=2\}\)\(x=4\) 所得立体体积
  3. 用两种方法表示区域 \(S=\{y\in[0,\pi],x=0,x=\sqrt{\sin y}\}\)\(y=4\) 所得立体体积
  4. \(\int_0^32\pi x^5d_x\) 是薄壳法表示的立体体积,逆推该立体的表达式
  5. 计算底半径为 r,高为 h 的正圆锥体积
  6. 用半径为 r 的圆筒削去半径为 R 的球体(\(R>r\))的中心,计算球体剩下部分的体积

提示

(1) \(\int_0^12\pi x\cdot e^{-x^2}d_x=\int_0^{-1}2\pi xe^{-x^2}/(-2x)d(-x^2)=\int_{-1}^0\pi e^td_t=\pi(1-e^{-1})\)

(2) \(\int_1^22\pi|x-4|\cdot x^2d_x=2\pi\int_1^2(4x^2-x^3)d_x=2\pi\left[\frac43x^3-\frac14x^4\right]_1^2=\frac{67}6\pi\)

(3)

  1. \(S=\{y\in[0,\pi],x=0,x=\sqrt{\sin y}\}\) 时,\(\|V\|=\int_0^\pi2\pi|y-4|\cdot\sqrt{\sin y}~d_y\)
  2. \(\|V\|=\int_0^\pi2\pi|y-4|\cdot\sqrt{\sin y}~d_y\) 时,\(\|V\|=\int_0^1\pi|(\sin^{-1}x^2-4)^2-(0-4)^2|d_x\)

(4) \(V=\{x\in[0,3],A(x)=2\pi x\cdot x^4\}\)

(5) \(\int_0^r2\pi x\cdot|h-\frac hr x|d_x=2\pi\frac hr\int_0^r(rx-x^2)d_x=2\pi\frac hr\left[\frac12rx^2-\frac13x^3\right]_0^r=\frac13\pi r^2h\)

(6) 球体剩下的部分为轴中心对称立体,于是可以方便地使用旋转法计算体积(区域按 x 轴旋转或按 y 轴旋转)

\(V=\{x\in[-\sqrt{R^2-r^2},\sqrt{R^2-r^2}],A(x)=\pi|\sqrt{(R^2-x^2})^2-r^2|\}\)

于是 \(\|V\|=\int_{-\sqrt{R^2-r^2}}^{\sqrt{R^2-r^2}}\pi|\sqrt{(R^2-x^2})^2-r^2|d_x=2\pi\int_0^{\sqrt{R^2-r^2}}(R^2-r^2-x^2)d_x=2\pi\left[(R^2-r^2)x-\frac13x^3\right]_0^{\sqrt{R^2-r^2}}=\frac43\pi(R^2-r^2)^{3/2}\)

4.功

5.函数的均值

函数 f 在 [a, b] 之间的平均值为:\(\overline f = \frac 1{b-a}\int_a^b f(x) d_x\)

积分中值定理

如果 f 在 [a, b] 之间连续,那么 [a, b] 内存在 c 使得: \(\int_a^b f(x) d_x = f(c) (b-a)\)

Tip

  1. \(f(x)=\cos^4x\sin x\)\([0,\pi]\) 上的均值为 \(\frac{\int_0^\pi\cos^4x\sin x~d_x}{\pi-0}=\frac{-\frac15\cos^5x\Big|_0^\pi}{\pi}=\frac2{5\pi}\)